嗨我需要为每个纬度和经度值的Alt值!然而,我从Geotiff提取这些值有问题,不知道这个数据是公开在Excel或CSV文件格式?不幸的是我还没有找到任何我自己!

我的数据:https://wetransfer.com/downloads/8e5ba7aa4cd8c406068a60543b353b0420220622084939/d3c744对于那些知道Python:

从mpl_toolkits.basemap进口基础图进口matplotlib。pyplot作为plt从进口gdal osgeo numpy linspace进口从numpy进口meshgrid熊猫作为pd导入基础图=基础图(投影= tmerc, lat_0 = 0,lon_0 = 3,llcrnrlon = 1。819757266426611,llcrnrlat = 41。583851612359275,urcrnrlon = 1。841589961763497,urcrnrlat = 41。598674173123)ds=gdal.Open(的路径。tif”)数据=ds.ReadAsArray ()数据=数据(::1:)x= linspace (0,Basemap.urcrnrx,data.shape [1])y= linspace (0,Basemap.urcrnry,data.shape [0])xx, yy=meshgrid(x, y) xx,yy =基础图(xx, yy,逆= True)print (xx)分= pd。数据Frame({ 'Latitude':yy.reshape(1),“经度”:xx.reshape(1),“高度”:数据.reshape (1)})# dfl.to_csv(“any-file.csv”)

这是我写的代码到目前为止,也许会有人在那里找到一个缺陷

嗨我需要为每个纬度和经度值的Alt值!然而,我从Geotiff提取这些值有问题,不知道这个数据是公开在Excel或CSV文件格式?不幸的是我还没有找到任何我自己!

我的数据:https://wetransfer.com/downloads/8e5ba7aa4cd8c406068a60543b353b0420220622084939/d3c744对于那些知道Python:

进口numpy作为np进口matplotlib。pyplot plt从netCDF4进口数据集熊猫作为pd导入数据=数据集(“C: /用户/奥利弗weis /桌面/本科/程序/ Daten / Daten / ETOPO1_Bed_g_gdal.grd”,“r”)打印(数据。variables.keys ())lon_range=数据。变量[' x_range '] [:]lat_range=数据。变量[' y_range '] [:]topo_range=数据。变量[' z_range '][:]间距=数据。变量(“间距”][:]维度=数据。变量(“维度”][:]z=数据。变量[' z '] [:]lon_num=维度[0]lat_num=维度[1]朗=np。linspace (lon_range [0],lon_range [1],维度[0])纬度=np。linspace (lat_range [0],lat_range [1],维度[1])威尼斯平底渔船=np.reshape(z,(lat_num,lon_num))分= pd.DataFrame ({“纬度”:纬度.reshape (1),“经度”:朗.reshape (1),“高度”:威尼斯平底渔船.reshape (1)})打印(公司分)

来源:ETOPO1区域选择(在python中)

这是我写的代码到目前为止,也许会有人在那里找到一个缺陷

嗨我需要为每个纬度和经度值的Alt值!然而,我从Geotiff提取这些值有问题,不知道这个数据是公开在Excel或CSV文件格式?不幸的是我还没有找到任何我自己!

嗨我需要为每个纬度和经度值的Alt值!然而,我从Geotiff提取这些值有问题,不知道这个数据是公开在Excel或CSV文件格式?不幸的是我还没有找到任何我自己!

我的数据:https://wetransfer.com/downloads/8e5ba7aa4cd8c406068a60543b353b0420220622084939/d3c744对于那些知道Python:

从mpl_toolkits。技术进口技术导入matplotlib。pyplot plt从osgeo import gdal from numpy import linspace from numpy import meshgrid import pandas as pd Basemap = Basemap(projection='tmerc', lat_0=0, lon_0=3, llcrnrlon=1.819757266426611, llcrnrlat=41.583851612359275, urcrnrlon=1.841589961763497, urcrnrlat=41.598674173123) ds = gdal.Open('Path.tif') data = ds.ReadAsArray() data = data[::-1,:] x = linspace(0, Basemap.urcrnrx, data.shape[1]) y = linspace(0, Basemap.urcrnry, data.shape[0]) xx, yy = meshgrid(x, y) xx, yy = Basemap(xx,yy,inverse=True) print(xx) dfl = pd.DataFrame({ 'Latitude': yy.reshape(-1), 'Longitude': xx.reshape(-1), 'Altitude': data.reshape(-1) }) #dfl.to_csv('any-file.csv')

这是我写的代码到目前为止,也许会有人在那里找到一个缺陷