取平均值重采样栅格

时间:2019-04-17 18:56:05

标签: python raster resampling

我有一个总体分辨率为100m的人口栅格,我想将其重新采样为更精细的分辨率75m。 我已经在栅格中拥有一个人口数据,其中每个像素值就是人数,即每100m * 100m以内的人数。当我需要将其重新采样到75m时,新的像素范围(min_r_x,max_r_x,min_r_y,max_r_y)可能会跨越100m栅格的几个像素(此处我使用红色,蓝色,绿色和黄色四个部分作为示例)。我的想法是,首先获取这个所需的75m栅格中每个像素的四个顶点的位置,然后获取顶点所在的填充值,然后计算面积百分比以生成重采样结果值。
您可以参考下面的图片作为简单说明。

X-MY-CSRF

我尝试将1行2列4 m分辨率的栅格重新采样为2 m分辨率的新栅格。

Original_R = 100 
 R = 75 
 Area_red = (max_x - min_r_x)*(max_r_y-min_y)
 Area_blue = (R – (max_x –min_r_x))*(max_r_y-min_y)
 Area_green = (max_x - min_r_x)*(R – (max_r_y-min_y))
 Area_yellow = (R – (max_x –min_r_x))*(R – (max_r_y-min_y))
 Value_resample = (10 *area_red     +20*area_blue+40*area_green+50*area_yellow)/(original_R* original_R)

0 个答案:

没有答案