我在UTM32中有一个GeoTiff,在UTM32中也有一个矩形坐标。 (这种预测可能并非总是如此,但预测总是相同的)
我只需要使用矩形裁剪图像。
矩形由下式给出:(xmin,xmax,ymin,ymax)
699934.584491, 700160.946739, 6168703.00544, 6169364.0093
我知道如何从点制作多边形,如何从多边形制作shapefile,我知道如何使用点创建蒙版的numpy数组。但是,我不知道如何使用多边形,shapefile或蒙版来实际裁剪图像。
我已经看过以下描述: https://pcjericks.github.io/py-gdalogr-cookbook/raster_layers.html#clip-a-geotiff-with-shapefile
然而,我并不是真的理解它,它似乎过于复杂。 (就像我不知道除了混乱之外应该做什么直方图拉伸)
答案 0 :(得分:2)
尝试使用bbox =(xmin,ymin,xmax,ymax)
from osgeo import gdal
bbox = (xmin,ymin,xmax,ymax)
gdal.Translate('output_crop_raster.tif', 'input_raster.tif', projWin = bbox)