d3.js - 从地理区域和栅格文件中获取多边形的工作流程是什么？

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

使用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使用问题中的命令进行处理。