从 Shapefile 中提取栅格值

Question

我有一个光栅文件，其中包含我为映射要素绘制的多边形或线条。我现在要做的是沿着这些 polygons/lines 从栅格数据中提取值，并绘制沿像素的高程图。如 How to extract an arbitrary line of values from a numpy array?。只是这次是多边形不是直线

Answer 1

您必须决定采样间隔。您可以按所需间隔沿 line/polygon 边添加点，然后提取这些点处的栅格值（使用 gdal/numpy）。

您需要注意光栅分辨率和采样间隔之间的关系，以避免 "skipping" 像素产生伪影或在同一像素中采集两个样本，并且您可能需要应用某种 filtering/interpolation 在生成的配置文件上。