栅格统计与投影的关系

结论

投影对栅格的统计结果影响较大。比较不同区域统计结果时需要统一投影。

地理坐标系和投影坐标系下，统计结果相差较大。如下表的第6条数据和第7、8条数据，像元数目约相差2000万个像元，占总像元数10%。DN总和约相差200万，占总DN总和的10%。

同一投影坐标系，不同投影中心（如4、5，6、7），统计结果有差别，但相差不大。洲尺度下（6、7），像元数目相差几千个，占总像元数0.001%。DN总和相差10万，占总DN总和的1%。

不同统计工具进行栅格的区域统计，统计结果略有差别，但差别不大。应该是在边界像元的处理上有不同。

统计结果

统计指定区域（矢量图形）内影像的影像值总和，分别尝试不同投影下的结果（重点关注Albers equal area conic投影在不同投影中心的差异）：

栅格分辨率15弧秒，投影后为500m。

较小的行政区域（area = 7.6km2）

小区域统计

序号	统计工具	坐标系统	像元数	像素值总和
1	python+gdal.Wrap cropToCutline=Fasle	地理坐标系 WGS84	43	5011.2
2	python+gdal.Wrap cropToCutline=True	地理坐标系 WGS84	36	3882.1
3	ArcTools  Zonal Statistics as Table	地理坐标系 WGS84	39	4652.24
4	ArcTools  Zonal Statistics as Table	投影坐标系：Asia_North_Albers_Equal_Area_Conic 地理坐标系：WGS84	32	3505.6
5	ArcTools  Zonal Statistics as Table	投影坐标系：Africa_Albers_Equal_Area_Conic 地理坐标系：WGS84	32	3532.1

非洲区域

大区域统计

序号	统计工具	坐标系统	像元数	像素值总和
6	ArcTools  Zonal Statistics as Table	地理坐标系 WGS84	147670192	11498671
7	ArcTools  Zonal Statistics as Table	投影坐标系：Asia_North_Albers_Equal_Area_Conic 地理坐标系：WGS84	119823937	8673634.4
8	ArcTools  Zonal Statistics as Table	投影坐标系：Africa_Albers_Equal_Area_Conic 地理坐标系：WGS84	119826271	8684811.6
9	python+gdal.Wrap cropToCutline=Fasle	地理坐标系 WGS84	147670167	11491753.7