ArcGIS填补修复DEM数据缺失值

问题:下载的DEM数据存在斑点状NoData空值;

     

方法(一):栅格计算器

 在ArcMap10.8中,ArcToolbox > Spatial Analyst工具 > 地图代数 > 栅格计算器;

输入公式:

Con(IsNull("inputraster"), FocalStatistics("inputraster", NbrRectangle(5,5, "CELL"), "MEAN"), "inputraster")

"inputraster"根据实际输入处理栅格图像;

方法(二):焦点统计

在ArcMap10.8中,ArcToolbox > Spatial Analyst工具 > 邻域分析 > 焦点统计;

可根据实际需求调整邻域分析区域(矩形/圆形/环形等)、大小(像元高度和宽度)、统计类型(MEAN或MAXIMUM等);

两种方法的区别:栅格计算器处理时保留非NoData值像元的原值,而焦点统计时根据邻域范围和大小对图像像元进行重新统计计算。

参考文章:

https://blog.youkuaiyun.com/qq_39546299/article/details/120980737

DEM数据下载来源:海南省12.5mDEM免费下载 (qq.com)[微信公众号:地质掘墓人]

### DEM 数据处理流程 DEM(Digital Elevation Model,数字高程模型)的数据处理通常涉及多个环节,这些环节涵盖了从原始数据获取到最终应用的全过程。以下是基于GIS和遥感技术的典型DEM数据处理流程: #### 1. 原始数据采集 DEM 的构建依赖于高质量的地形数据源。常见的数据来源包括卫星影像、航空摄影测量、激光雷达扫描等[^1]。 例如,SRTM(Shuttle Radar Topography Mission)、ASTER GDEM 和 LiDAR 是常用的全球或局部区域高精度DEM数据。 #### 2. 数据预处理 在实际使用前,需要对原始数据进行必要的校正和清理工作: - **辐射校正**:消除传感器响应差异的影响。 - **几何校正**:调整图像投影误差,确保地理坐标的准确性。 - **去云处理**:对于光学影像,去除因云层遮挡造成的无效像素。 - **重采样**:统一不同分辨率数据的空间尺度[^2]。 #### 3. 地形提取与建模 通过专业的GIS软件(如ArcGIS 或 QGIS),可以从经过预处理后的数据中提取地形特征并生成DEM文件。具体方法如下: - 利用立体像对计算高度差; - 使用LiDAR点云生成TIN(Triangulated Irregular Network)三角网格结构后再转换成规则格网形式的DEM; - 应用InSAR干涉测量技术获得地面起伏信息。 #### 4. 后期编辑与质量控制 为了提高DEM的质量,需执行一系列后期编辑操作: - 平滑滤波减少噪声干扰; - 插值填补缺失值; - 边界裁剪限定研究范围; - 叠加其他辅助数据验证其合理性。 #### 5. 空间分析与应用 完成上述步骤后即可开展各种空间分析任务,比如坡度坡向计算、流域划分、视线分析等。此外还可以结合降雨分布预测潜在泥石流发生概率或者规划最佳交通线路避开陡峭地段等问题解决方案提供依据。 ```python import arcpy from arcpy.sa import * # 设置工作环境 arcpy.env.workspace = r"C:\path\to\data" # 输入DEM数据 dem = Raster("input_dem.tif") # 计算坡度 slope = Slope(dem, "DEGREE", 1) # 输出结果 slope.save(r"C:\path\to\output\slope_output.tif") ``` ---
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值