我正在使用arcpy进行脚本编写,我有一个由作物按县组成的NumPy矩阵。该矩阵用浮点数填充。
我还有一个单独的栅格图像显示作物以及一个显示县的形状文件。
我想将我的NumPy矩阵内容保存到新的栅格图像中,使用矩阵中给定的正确值,因为栅格像素是特定作物在特定县的情况。
我的问题是:我如何改进将矩阵转换为栅格图像的方式?
目前我做类似以下的事情:
我还有一个单独的栅格图像显示作物以及一个显示县的形状文件。
我想将我的NumPy矩阵内容保存到新的栅格图像中,使用矩阵中给定的正确值,因为栅格像素是特定作物在特定县的情况。
我的问题是:我如何改进将矩阵转换为栅格图像的方式?
目前我做类似以下的事情:
# initialize my output
outraster = 0 * cropraster
for crop in crops
for county in counties
outraster = outraster + \
conditional statements checking if cropraster pixel is in
current county and crop, and if so return matrix value
outraster.save(filename)
条件语句是嵌套的空间分析“Con”调用,其语法我将为您省略。
这种方法可以工作,但非常丑陋,我觉得应该有更好的方法来实现我的目标。当我尝试太多县或作物时,上述方法会崩溃。它在保存步骤时失败。
我想知道是否有更简单的方法来实现我想要的功能。有什么建议吗?