mayvai使用tif绘制地形图

Question

问题遇到的现象和发生背景

使用在国家地理云数据平台下载的DEM高程数据绘制地形图时，x与y轴似乎合在了一起

问题相关代码，请勿粘贴截图

import numpy as np
from mayavi import mlab
from osgeo import gdal
gdal.UseExceptions()

def read(filename):
    ds = gdal.Open(filename)
    elev = ds.ReadAsArray()

    # True x, y coordinates
    x0, dx, dxdy, y0, dydx, dy = ds.GetGeoTransform()
    i, j = np.mgrid[:elev.shape[0], :elev.shape[1]]
    x = x0 + dx * j + dxdy * i
    y = y0 + dy * i + dydx * j

    return ds.ReadAsArray(), x, y

z, x, y = read('ASTGTMV003_N02W064_dem.tif')


fig = mlab.figure()
mlab.mesh(x, y, z, colormap='gist_earth')
mlab.axes()
mlab.show()

运行结果及报错内容

我的解答思路和尝试过的方法

使用原GitHub中的tif文件则可绘出地形图 https://github.com/joferkington/scipy2015-3d_printing

我想要达到的结果

问题出在哪

Answer 1

代码是用于读取tif图形文件立体呈现地形，在用数据文件data\alaska\clipped_elev.tif正常显示，使用其他数据文件显示不正常，说明代码没有问题，是数据文件本身问题，检查数据格式是否符合要求。