使用GDAL进行汇总/重采样?

时间:2018-12-19 18:22:31

标签: gis gdal

我正在开发一个程序,该程序将使用GDAL以不同的分辨率处理各种用户选择的全球栅格数据集,以此作为将预计的国家级污染排放分配给网格正方形的基础。

我一直尝试使用的方法如下:

  1. 首先,使用给定国家/地区的多边形形状的形状文件裁剪全局栅格数据集。
  2. 接下来,对数据集重新采样,以使它们都具有相同的分辨率(像元大小),相同的XLL / YLL角和相同数量的行/列。

所有重新采样的数据集都将比原始的全局代理数据集(单元大小更大)更粗糙,我想使用总和聚合方法。

我已经能够使用GDAP warp来执行步骤1,但是还不能使步骤2正常工作。我知道进行聚合的正确方法是使用“ -r average”重采样方法,然后将结果乘以原始网格与重采样网格中的平方数之比。但是我不清楚如何创建都具有完全相同的行/列和xll / yll角的重采样网格。

有人可以建议使用gdal_warp,gdal_translate或类似方法执行此操作的最佳做​​法吗?

注意:我正在从另一个应用程序内调用这些实用程序-因此,确实需要一种无需用户交互即可自动(即从命令行)调用的解决方案。生成的网格将显示在用Delphi编写的Windows应用程序的屏幕上,该应用程序使用MapWinGIS来显示地图。

提前感谢您的任何建议!

0 个答案:

没有答案