擅长:python、mysql、java
<p>我同意这个解决方案,但我想补充一点,如果你和我一样,我通常有一个<code>rasterio datasetreader</code>对象(使用<code>rasterio.open</code>读取georeffed光栅数据的结果),而不仅仅是一个原始numpy数组。因此,对于RasterioV1.1.8,我必须执行额外的步骤,从datasetreader对象提取numpy数组。例如,对于单个波段:</p>
<pre><code>dem = rasterio.open("GIS/anaPlotDEM.tif")
fig, ax = plt.subplots(figsize=(10,10))
image_hidden = ax.imshow(dem.read()[0])
fig.colorbar(image_hidden, ax=ax)
rasterio.plot.show(dem, ax=ax)
</code></pre>
<p>(我想加上这句话作为评论,但没有声誉点)</p>