Получить вертикальное (z) разрешение растрового файла матрицы высот (используя gdal & python)? - PullRequest
Новая
       12

Получить вертикальное (z) разрешение растрового файла матрицы высот (используя gdal & python)?

0 голосов
/ 30 октября 2018

Я хотел бы найти разрешение по вертикали (z) растра DEM (цифровая модель рельефа).

Используя gdal, я легко могу найти разрешение x и y и ширину пикселя: 

dem = gdal.Open('DEMraster50.tif')    #DEM is a rectangular tif 
ncol = dem.RasterXSize                #number of columns (aka number of cells along x axis)
nrow = dem.RasterYSize                #number of rows (aka number of cells along y axis)
ulx, pixelwidthx, xskew, uly, yskew, pixelheighty = dem.GetGeoTransform() 
#get resolution and coordinate info (for some reason the order of skew and pixel size is flipped for y axis?!)

Я также могу найти минимальное и максимальное значения z (высота в метрах): 

srcband = dem.GetRasterBand(1)        #get the first band (with the elevation data)
zmin,zmax,zmean,zstdv = srcband.GetStatistics(True, True)  #get stats for this band (automatically ignores nan values)

Но то, что я хочу, это эквивалент nrows и ширины пикселя в направлении z. Это возможно? Если так, то как? Спасибо!

EDIT: Извините, у меня дерьмовая камера, но это может помочь. Я ищу Nlay и Z толщины. rastergrid

...