从地理区域和光栅文件获取多边形的工作流程是什么？

Question

鉴于与欧洲相对应的TIFF文件，并转换为形状轮廓：

gdal_contour -a elev Europe_input.tif Europe_output.shp -fl -1000 -500 -200 -50 0 50 100 200 500 1000 2000 4000 6000

给定对应法国的定义区域：

West : 005° 48' W
East : 010° E
North : 051° 30' N
South : 041° N

如何仅使用gdal获取几何与我定义的区域相交的要素？

Answer 1

Clip the raster到gdal_translate感兴趣的区域，例如：

gdal_translate -projwin ulx uly lrx lry Europe_input.tif Europe_subregion.tif

其中ulx uly lrx lry是感兴趣区域的上/下x / y范围，单位为Europe_input.tif（也许是米，也许是弧秒，也许是度数，我无法分辨）。剪切光栅后，使用问题中的命令处理Europe_subregion.tif gdal_contour。