使用Mapnik comp-op使用矢量图层屏蔽栅格图层

时间:2015-01-04 13:50:25

标签: python gis gdal mapnik cartography

我正在尝试使用单独的陆地和海洋数据集制作Mapnik地图。

Land是Landsat图像中的栅格数据集。 Ocean是使用NOAA水深测量的gdaldem创建的栅格数据集。

Landsat数据集包括海洋图像,但它很难看。我想敲掉它,所以我们看到下面的NOAA测深层,创建一个类似于蓝色大理石地图的外观。

我想我可以使用包含基于世界边界矢量数据集的PolygonSymbolizer的第三层作为掩码,方法是将符号化程序的comp-op设置为dst_in并将其放在图层堆栈的最顶层,如下所示:

  • VECTOR MASK(COMP-OP DST_IN)
  • LAND RASTER
  • SEA RASTER

我已尝试过每个alpha comp模式,无法让它工作。我错过了什么?或者也许有一种更简单的方法?

这是我的代码的要点...... https://gist.github.com/john-sanders/800b973858bc13e4417c

1 个答案:

答案 0 :(得分:0)

退后一步,似乎我正在使用错误的工具来完成工作。或者至少不是最简单的工具。

gdalwarp有一个-cutline参数,它接受一个shapefile来删除不需要的栅格数据。

https://gis.stackexchange.com/questions/45053/gdalwarp-cutline-along-with-shapefile

此外,似乎在gis.stackexchange.com上更好地询问了这个问题

相关问题