做个小记录自用,具体操作步骤是没有记录完整的,谨慎参考。
1 看看就行 pro插值
看看过程就行。有详细过程,类似tutorial
https://learn.arcgis.com/zh-cn/projects/interpolate-temperatures-using-the-geostatistical-wizard/
2 一些注意
2.1 注意用投影坐标系
wgs84转投影坐标系。一般进行arcgis的空间分析操作都建议用投影坐标系
。常用为UTM-WGS1984下的。
https://blog.csdn.net/qq_40323256/article/details/120751710
2.2 修复几何
做插值时有时会报错,例如空数据之类,附带一个报错问题类似 000358: Invalid expression(注:仅示意报错问题样式,非报错问题)。
点进去提示需要检查几何或者修复几何
.
建议在插值数据正式插值前,用修复几何处理下,因为数据量的关系,有时会有个别数据经纬度等存在问题但没有发现。
3 请使用-探索性插值
提供多种插值方法,最后选出精度较好的一种
Exploratory Interpolation
https://pro.arcgis.com/zh-cn/pro-app/latest/tool-reference/geostatistical-analyst/exploratory-interpolation.htm
3.1 确保字段为数值型
可在数据中查看。如果不是数值型,可通过新建一个字段。设置为数值型,比如长整型,再通过字段计算器将新字段的值=原字段值。
3.2 数据服从正态分布
右击图层,创建图表,生产直方图。
如果不符合正态部分,需要自己调整变化。平均值和中值接近,偏度较小即可
但这里视数据情况定,
3.3 趋势分析-arcmap
注意这个只有arcmap里面才有,
详见https://zhuanlan.zhihu.com/p/636377714
arcpro目前没有引入这个功能的,可以见esri社区的这个讨论帖:
https://community.esri.com/t5/arcgis-pro-ideas/add-trend-analysis-tool-in-arcgis-pro/idi-p/1124000
3.4 插值
3.5 如果数据无法覆盖所需范围
一个小问题就是,如果你的点无法覆盖你想要的范围,比如下图,那么你可能需要自己将前面探索的最好的一种,自己在工具箱里手工插值
举个例子,比如我这里提示我最佳方法是:经验贝叶斯克里金法。
就自己在工具箱里,找到该方法,差值。
在环境中提供了范围选项。
这是上一步探索性差值没有的
另,实际上链接https://pro.arcgis.com/zh-cn/pro-app/latest/tool-reference/geostatistical-analyst/exploratory-interpolation.htm的参数,指的是下图中用地统计向导的界面参数。
但是,这个界面也没有空间范围可以选!!
插值出来还是,覆盖不了所需范围的。
3.6 如果要导出插值结果书属性表
经常做分析不是在矢量或者栅格要素上操作,而是属性表数据。
但存在一个问题,就是有时候数据量过大,比如上百万,常规导出的方法,比如复制属性表(会复制不完),text to table(2025年了还限制65533行) 等,是行不通的。
而直接读取shp文件对应的dbf文件,由于arcgis的编码问题,会读取失败,比如我就遇到了编码为Windows-1254这种,能正常显示,就是无法python读取的编码。
经试验,推荐以下方法,可以 快速、无痛导出含有大规模数据的属性表,如下图