亚洲在线久爱草,狠狠天天香蕉网,天天搞日日干久草,伊人亚洲日本欧美

為了賬號安全,請及時綁定郵箱和手機立即綁定
已解決430363個問題,去搜搜看,總會有你想問的

在 Python 底圖上鏡像的 GeoTIFF 柵格

在 Python 底圖上鏡像的 GeoTIFF 柵格

達令說 2023-07-27 15:42:29
我正在嘗試使用底圖在地圖上繪制 .tif 柵格。使用 QGIS,我看到了它應該是的柵格層:?QGIS 圖像然而,當使用 python 底圖繪制時,顏色關閉,并且投影以某種方式旋轉 180 度并進行鏡像,在圖表的左側投影一條隨機藍線:有關 .tif 文件的一些信息(通過 rasterio 和 Earthpy 包獲?。?lt;osgeo.gdal.Dataset; proxy of <Swig Object of type 'GDALDatasetShadow *' at 0x7fcf9bdbef90> >{'driver': 'GTiff', 'dtype': 'float32', 'nodata': -3.4028234663852886e+38, 'width': 2760, 'height': 1350, 'count': 1, 'crs': CRS.from_epsg(4326), 'transform': Affine(0.01, 0.0, 3.7,? ? ? ?0.0, -0.01, 71.2)}EPSG:4326BoundingBox(left=3.7, bottom=57.7, right=31.3, top=71.2)+proj=longlat +datum=WGS84 +no_defs原始柵格數據在此處下載(森林恢復潛力)并使用 QGIS 裁剪到正確的緯度和經度。我究竟做錯了什么?import gdalfrom numpy import linspacefrom numpy import meshgridimport matplotlib.pyplot as pltfrom mpl_toolkits.basemap import Basemaplowlong = 3.7 #lower left corner of longitudelowlat = 57.7 #lower left corner of latitudeupplong = 31.3 #upper right corner of longitudeupplat = 71.2 #upper right corner of latitudepathToRaster = r'~/Data/Shapefile/NorwayPotential.tif'raster = gdal.Open(pathToRaster,1)print(raster)geo = raster.GetGeoTransform()geo = raster.ReadAsArray()mp = Basemap(projection='merc',? ? ? ? ? ? ?llcrnrlon=lowlong,? ? ? ? ? ? ?llcrnrlat=lowlat,? ? ? ? ? ? ?urcrnrlon=upplong,? ? ? ? ? ? ?urcrnrlat=upplat,? ? ? ? ? ? ?resolution='i')mp.drawcoastlines()mp.drawcountries()x = linspace(0,mp.urcrnrx,geo.shape[1])y = linspace(0,mp.urcrnry,geo.shape[0])xx,yy = meshgrid(x,y)mp.pcolormesh(xx,yy,geo)plt.show()
查看完整描述

1 回答

?
寶慕林4294392

TA貢獻2021條經驗 獲得超8個贊

代碼行之后:

geo = raster.ReadAsArray()

您可以通過以下方式翻轉數據數組

geo =  geo[::-1,:]

并應該得到正確的結果。


查看完整回答
反對 回復 2023-07-27
  • 1 回答
  • 0 關注
  • 125 瀏覽
慕課專欄
更多

添加回答

舉報

0/150
提交
取消
微信客服

購課補貼
聯系客服咨詢優惠詳情

幫助反饋 APP下載

慕課網APP
您的移動學習伙伴

公眾號

掃描二維碼
關注慕課網微信公眾號