带有GDAL和PostGIS的Python替代了ArcGIS

几次前，我处理了阿拉斯加大小的每日地面温度卫星图像。

我必须处理大约一千张图像。我下载了它们，进行镶嵌，投影，修剪到shapefile边界，从修剪的栅格中提取均值，然后在时间-温度图上绘制均值。时间-温度图是最终输出。我在ArcGIS中完成了所有这些工作，即使我创建了有助于自动完成一部分工作的ArcPy脚本，也仍然很艰辛。

现在，我想构建一个执行相同操作但自动100％完成的程序。我只想将拐角的坐标输入到程序中，以便定义感兴趣的区域并获得生成的最终输出，即时间-温度图。

您认为使用Python与GDAL库和PostGIS作为栅格数据存储库是否可行？

是的-简短的答案。我认为您将需要NumPy和SciPy作为Python解决方案的一部分。看一下scipy.ndimage模块，它涉及到计算均值（这些区域性均值吗？），因为这比仅使用NumPy进行计算要快得多。另外，通过使用Python 多处理模块，您将获得显着的速度提升。但是，请减少子进程的数量，因为栅格处理可能会占用大量处理器，并且最终可能会降低计算机的速度而不是加快速度。我建议比您的计算机上的内核少1个进程。