结论
投影对栅格的统计结果影响较大。比较不同区域统计结果时需要统一投影。
地理坐标系和投影坐标系下,统计结果相差较大。如下表的第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 |
|