1
rasterioとgeopandasを使用してラスターをクロップする
歴史的な航空写真を切り抜いています。これらの写真は、端に大きな黒い領域があります(値0)。ただし、値が0の有効なデータもあります。私が使用しているワークフローは: rasterioでラスターを読み込む rasterio.features.shapes()を使用してラスターをポリゴン化します 値が0でサイズが5000平方メートルを超えるポリゴンを特定する 元の画像をポリゴンでマスクし、反転マスクを実行します これが単一の画像をマスクするための私の現在のコードです: import rasterio from rasterio import features from rasterio import mask import json import geopandas as gpd results = [] final_results = [] with rasterio.open(r"C:\1927_oahu\tif\_Line1_6to8_0.tif") as src: src_meta = src.meta src_affine = src_meta.get("transform") band = src.read(1) for geometry, raster_value in features.shapes(band, transform=src_affine): if raster_value == …