我正在开发一个程序,该程序将使用GDAL以不同的分辨率处理各种用户选择的全球栅格数据集,以此作为将预计的国家级污染排放分配给网格正方形的基础。
我一直尝试使用的方法如下:
所有重新采样的数据集都将比原始的全局代理数据集(单元大小更大)更粗糙,我想使用总和聚合方法。
我已经能够使用GDAP warp来执行步骤1,但是还不能使步骤2正常工作。我知道进行聚合的正确方法是使用“ -r average”重采样方法,然后将结果乘以原始网格与重采样网格中的平方数之比。但是我不清楚如何创建都具有完全相同的行/列和xll / yll角的重采样网格。
有人可以建议使用gdal_warp,gdal_translate或类似方法执行此操作的最佳做法吗?
注意:我正在从另一个应用程序内调用这些实用程序-因此,确实需要一种无需用户交互即可自动(即从命令行)调用的解决方案。生成的网格将显示在用Delphi编写的Windows应用程序的屏幕上,该应用程序使用MapWinGIS来显示地图。
提前感谢您的任何建议!