流程如下：

第一步：

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.进行区域统计。空间分析工具—区域—区域统计作为表。选土地利用标准数据的value作为统计基础，对hkhcarbonNamaskmulHWSD4进行区域统计。

第四步：

SoilGrid数据只需要经过：

1.定义坐标系统，坐标系转化

2.将hkhcarbon的NA掩膜乘以上述数据; "hkhcarbonNaNaother1.tif" * "SoilGridGEEFromMollweideToAlbers.tif（2投影转化）"

3.将上述结果与土地利用标准数据连用获取不同土地利用的区域统计

第五步：

HKHcarbon250数据处理：

1.按照土地利用标准数据，制作水冰雪的掩膜。Con(("MosaicHKHAlbersFinalResampleInt.tif" == 15) | ("MosaicHKHAlbersFinalResampleInt.tif" == 17),0,1)

2.将上述数据乘HKHcarbon250数据。"hkhcarbon250.tif" * "MASKwatersnowice.tif"

3.再利用土地利用标准数据统计HKHCARBON。