区域统计3类产品中不同植被类型的土壤碳密度

流程如下：
第一步：
R语言估算完成后，先把hkhcarbon数据中NA值的区域作为掩膜：

1. 通过isnull函数，将na值设置为0，其余值设置为1；Con(IsNull('hkhcarbon250.tif'),0,1)
2. 通过setnull函数，将0值再设置成na，其余值不变仍为1。SetNull('hkhcarbonNa0.tif' == 0,1)

第二步：
在R中将MOD12土地利用分类数据向hkhcarbon重采样，得到边界范围像元完全匹配的土地利用标准数据。
第三步：
需要明确SoilGrid数据和HWSD数据的差异。HWSD数据是将江河湖泊设置为了NA，而SoilGrid数据是将Na值设为0。所以在统计之前一定将HWSD数据的N值转化为0.
HWSD数据操作流程：
1.裁剪HKH正常边界更大的区域
2.重投影为HKHAlbers坐标系，并将分辨率设置为250m
3.通过is将Na值转化为0，其余值设为1. Con(IsNull('HWSD大范围albersCublic去异常.tif(4去除负数异常值)'),0,'HWSD大范围albersCublic去异常.tif(4去除负数异常值)')或者是 Con(IsNull('HWSD大范围albersCublic去异常.tif(4去除负数异常值实际到这一步就可以了)'),0,'HWSD大范围albersCublic去异常.tif(4去除负数异常值实际到这一步就可以了)')
4.将hkhcarbon的NA掩膜乘以HWSD数据中Na值转化为0的数据; 'hkhcarbonNaNaother1.tif' * 'HWSD大范围AlbersCubic去异常Na0.tif(5为匹配hkhcarbon将其Na值转化0)'
5.进行区域统计。空间分析工具—区域—区域