terra::project中的加权求和:权重具体指什么?
terra::project() method="sum" 的权重机制说明
当使用terra包的project()函数并设置method="sum"时,加权和所采用的权重是输入栅格像元与输出栅格像元的重叠面积占输入像元总面积的比例。
具体逻辑如下:
- 投影坐标系转换时,输入栅格的单个像元可能被拆分,部分区域落入输出栅格的多个像元中;同理,输出的一个像元也可能包含多个输入像元的部分区域。
- 对于每个输出像元,会遍历所有与其有重叠的非NA输入像元:计算该输入像元在当前输出像元内的重叠面积,除以输入像元的总面积得到权重比例,再用这个比例乘以输入像元的原始值,得到该输入像元对当前输出像元的贡献值。
- 最终输出像元的值就是所有符合条件的输入像元贡献值的总和。
举个简单例子:
- 输入像元A面积为100,值为5,其中30%的区域落在输出像元X内,那么A对X的贡献为
5 × 0.3 = 1.5 - 输入像元B面积为80,值为10,其中50%的区域落在输出像元X内,那么B对X的贡献为
10 × 0.5 = 5 - 输出像元X的最终值就是
1.5 + 5 = 6.5
内容的提问来源于stack exchange,提问作者Jean-Luc Dupouey
相关产品推荐
相关产品推荐

