需求与目的

在批量将全球GeoTif裁剪得到某一区域之后(具体操作可参考我的另一篇博客使用shp文件批量裁剪GeoTif遥感图像),在ArcGIS打开图像发现,图像中有部分像元是空洞即值为Nodata(如下图所示)
2019-11-14 GeoTif中Nodata像元填充
可以看出,Nodata空洞一方面会影响ArcGIS中图像显示的美观效果,而且在使用ArcGIS进行其它计算分析的时候也会造成影响,故而需要将Nodata填充为某一个值。

运行环境

Windows 7、10(其它系统版本未尝试);
ArcGIS10.5自带Python 2.7.12,因为要调用ArcGIS中的ArcPy模块,所以只能使用ArcGIS自带的Python 2.7.12版本;
ArcGIS需要在扩展模块中勾选 Spatial Analyst模块。如下图所示。
2019-11-14 GeoTif中Nodata像元填充

Python代码实现

代码如下:
2019-11-14 GeoTif中Nodata像元填充
最核心的代码:data = Con(IsNull(data), FocalStatistics(data, NbrRectangle(3, 3, “CELL”), “MEAN”, “DATA”), data)

结果

2019-11-14 GeoTif中Nodata像元填充

用ArcGIS打开看一下:
2019-11-14 GeoTif中Nodata像元填充
OK,原图像中Nodata的像元已被其周围3*3像元值的平均值所填充!
最后,如果本文有帮助到你,请为我点个赞打个赏,十分感谢!

相关文章:

  • 2022-02-04
  • 2021-05-19
  • 2022-02-12
  • 2022-12-23
  • 2021-07-31
  • 2021-10-21
  • 2021-04-29
  • 2022-02-06
猜你喜欢
  • 2021-07-30
  • 2021-08-10
  • 2021-04-20
  • 2022-12-23
  • 2021-12-31
  • 2022-01-09
  • 2022-12-23
相关资源
相似解决方案