Маска растрово в питоне с использованием rasterio ⇐ Python

[Anonymous]

Я хочу обрезать один растер, основываясь на степени другого (меньшего) растра. Сначала я определяю координаты углов меньшего растра, используя < /p>

Код: Выделить всё

import rasterio as rio
import gdal
from shapely.geometry import Polygon

src = gdal.Open(smaller_file.tif)
ulx, xres, xskew, uly, yskew, yres  = src.GetGeoTransform()
lrx = ulx + (src.RasterXSize * xres)
lry = uly + (src.RasterYSize * yres)
geometry = [[ulx,lry], [ulx,uly], [lrx,uly], [lrx,lry]]

Это дает мне следующую выходную геометрию = [-174740.0, 592900.0], [-174740.0, 2112760,0], [900180.0, 2112760.0], [900180.0, 592900.0]

Подробнее здесь: https://stackoverflow.com/questions/717 ... g-rasterio