1

我想制作一个光栅图像的子集并将其放在 800x600 的尺寸中。我正在查看 Rasterio 食谱，但它似乎不允许我输入诸如 800x600 之类的尺寸。这是我一直在看的：https ://mapbox.s3.amazonaws.com/playground/perrygeo/rasterio-docs/cookbook.html

另外，我看到了这个并认为它可能会起作用：https ://rasterio.readthedocs.io/en/latest/topics/windowed-rw.html

我使用了阅读器代码片段：

import rasterio
with rasterio.open('MyRasterImage.tif') as src:
    w = src.read(1, window=Window(0, 0, 800, 600))

print(w.shape)

但是，当我运行它时，它给了我错误消息：

w = src.read(1, window = Window(0, 0, 800, 600))

NameError: name 'Window' is not defined

不知道是什么导致了这个错误。我在想 Windows 是 rasterio 中的一个内置函数，我可以简单地调用它并调整图像大小以制作一个子集。

我还希望能够在屏幕上显示新的 800x600 图像（使用 Spyder），但不确定这是如何完成的。

任何和所有的帮助都将不胜感激，并将投票赞成。

谢谢

4

1 回答 1

1 
import numpy
import rasterio
from matplotlib import pyplot
from rasterio.windows import Window

width = 800
height = 600

with rasterio.open('MyRasterImage.tif') as src:
    w = src.read(1, window=Window(0, 0, width, height))
    profile = src.profile
    profile['width'] = width
    profile['height'] = height
    # Create output
    result = numpy.full((width, height), dtype=profile['dtype'], fill_value=profile['nodata'])

#writting
with rasterio.open('/tmp/sampled_image.tif', 'w', **profile) as dataset:
    dataset.write_band(1, result)

#plotting
with rasterio.open('/tmp/sampled_image.tif') as src:
    pyplot.imshow(src.read(1), cmap='pink')
    pyplot.show()
于 2020-07-17T22:36:05.087 回答