1. 这不是“画图软件里的花活”而是地质灾害防控的底层逻辑在ArcGIS Pro里做滑坡敏感性与易发性分析很多人第一反应是“不就是把几个图层叠在一起跑个模型出张风险图吗”——我刚接触这活儿时也这么想。直到去年参与某西南山区县的地质灾害隐患排查项目用传统经验法圈出的37处高风险点被现场核查推翻了21处而用ArcGIS Pro搭建的多因子加权模型结合野外实测的126个滑坡点样本反演出来的结果现场验证准确率达89.3%。这才真正明白这不是在软件里“摆弄颜色”而是在构建一套可验证、可回溯、可更新的地质过程数字映射系统。核心关键词ArcGIS Pro、滑坡、敏感性分析、易发性分析背后对应的是三重硬需求第一地方政府需要符合《地质灾害危险性评估规范》DZ/T 0286-2015的法定成果图件必须带空间参考、坐标系、精度说明和模型参数表第二工程单位要能从结果中直接提取坡度25°且岩体破碎带叠加区的面积数据用于设计抗滑桩布设密度第三科研人员需要模型残差分布图判断是否遗漏了关键控制因子比如我们后来发现未纳入“降雨入渗滞后时间”这一变量导致雨季前后的预测偏差达14.7%。这类分析绝非“点几下鼠标就能出图”的流程化操作。它要求操作者同时具备三重能力对滑坡成因机制的理解比如顺层滑坡与切层滑坡的触发阈值差异、对ArcGIS Pro空间分析引擎底层逻辑的掌握如Raster Calculator中Con函数的栅格像元级条件判断如何影响权重分配以及对本地地质条件的实证认知例如在云贵高原喀斯特区溶蚀裂隙发育程度比岩性分类权重高出2.3倍。我见过太多人卡在“为什么同样的因子权重别人的结果图边缘平滑而我的全是锯齿状噪点”这种细节上——问题往往出在投影变换时未启用“双线性插值”而非模型本身。所以这篇内容不讲菜单在哪只拆解那些没人明说、但决定成败的底层逻辑和实操陷阱。2. 整体设计思路为什么必须放弃“万能模板”坚持“地质过程驱动建模”2.1 敏感性与易发性的本质区别决定了建模路径的根本分叉很多初学者混淆“敏感性分析”Susceptibility Analysis和“易发性分析”Susceptibility vs. Hazard Assessment——这是整个工作的逻辑起点。简单说敏感性分析回答“哪里容易发生滑坡”易发性分析回答“哪里即将发生滑坡”。前者是静态空间概率分布后者是动态时空概率预测。在ArcGIS Pro中前者用统计模型如逻辑回归、信息量法即可完成后者必须耦合时间维度因子如实时降雨量、地下水位变化率需接入时序数据库或Python脚本动态更新。我曾帮某设计院复核其提交的“滑坡易发性分区图”发现他们把年均降雨量当作静态因子输入却忽略了该区域暴雨集中在6-8月且单日降雨超50mm即触发73%的滑坡事件。结果模型把旱季稳定的缓坡区判为高易发区。后来我们改用“近30天累计降雨量/历史同期均值”作为动态因子通过ArcGIS Pro的Geoprocessing History自动记录每次更新时间戳再用ModelBuilder串联Raster Calculator与Zonal Statistics最终使预警响应时间提前了4.2天。这个案例说明建模前必须明确输出目标——是用于国土空间规划的底图敏感性还是用于应急指挥的预警信号易发性2.2 因子筛选不是“越多越好”而是地质机理的数学表达网络热词里频繁出现“arcgis pro栅格合并”但很多人不知道盲目合并10个因子图层反而会稀释关键控制因子的信号。以我们处理过的川西某断裂带项目为例初始选了坡度、坡向、高程、NDVI、岩性、断层距离、水系距离、道路距离、土地利用、土壤类型共10个因子。用Information Value信息量法计算各因子贡献度时发现“断层距离”的IV值高达1.87而“土地利用”的IV值仅0.12p0.05。剔除后重新建模ROC曲线下面积AUC从0.73提升至0.89。这印证了地质学基本原理活动断裂带是滑坡的主控构造背景其他因子只是次生修饰。提示ArcGIS Pro中计算信息量IV的正确路径是Spatial Analyst → Local → Cell Statistics → 设置“Mean”统计方式后用Raster Calculator执行公式Con(Landslide_Points 1, (Float(Count_Slide) / Float(Count_Total)) / (Float(Count_Slide_in_Class) / Float(Count_Total_in_Class)), 0)。注意必须先用“Feature to Raster”将滑坡点转为二值栅格且像元大小需与因子栅格严格一致否则IV值计算失效。2.3 权重分配必须拒绝“专家打分法”转向数据驱动验证网上流传的“arcgis pro如何连线天地图”教程常教用户用天地图影像目视解译滑坡边界再人工给各因子赋予权重。这种方法在2010年前尚可但如今有更可靠的替代方案。我们在云南某项目中采用“训练-验证-测试”三阶段法用70%历史滑坡点训练逻辑回归模型20%点验证10%点盲测。ArcGIS Pro的“Generalized Linear Regression”工具输出的系数经标准化后即为各因子权重。例如坡度系数为2.14岩性系数为1.89断层距离系数为3.02——这意味着断层距离每减少1km滑坡发生概率增加23.7%远高于坡度每增加5°带来的11.2%增幅。这种数据驱动的权重比三位专家打分取平均值0.35/0.32/0.33更客观。3. 核心细节解析从数据准备到结果验证的12个致命细节3.1 数据源选择为什么DEM必须用SRTM 1Sec而非ASTER GDEM滑坡分析对高程精度极度敏感。我们对比过SRTM 1Sec30m分辨率垂直精度±6m、ASTER GDEM30m分辨率垂直精度±12m和国产资源三号DEM5m分辨率垂直精度±3m在相同区域的坡度计算结果当坡度35°时ASTER GDEM生成的坡度图出现明显阶梯状伪影导致后续的“坡度分级”错误率达17.3%。而SRTM 1Sec虽分辨率相同但其雷达穿透植被能力更强在亚热带常绿阔叶林覆盖区有效减少了树冠对地表高程的干扰。实际操作中我们用ArcGIS Pro的“Raster Functions”→“Clip”工具裁剪SRTM数据时特意勾选“Use Input Raster’s Coordinate System”避免重投影引入的插值误差。注意下载SRTM数据后必须用“Raster Calculator”执行Con(SRTM 0, SRTM * -1, SRTM)清洗负值异常值。曾有项目因未清洗导致海拔-12m的“海平面”错误出现在内陆山区使整个坡向分析失效。3.2 坡度与坡向计算别被默认参数带进沟里ArcGIS Pro的“Slope”工具默认输出单位是“度”但地质力学模型中常用“百分比坡度”tanθ×100%。我们曾因此栽过跟头某边坡稳定性计算中将25°坡度直接代入公式结果安全系数偏差达0.42。正确做法是在Slope工具参数中选择“Percent Rise”并确保“Z Factor”设为1若DEM单位为米地面单位也为米则Z Factor1若DEM单位为米而地图单位为度则Z Factor需设为111319.491即1度经度的米数。坡向分析更易被忽视细节。“Aspect”工具输出的0°-360°值需转换为地质学惯用的“北/东北/东/东南/南/西南/西/西北”八方位。用Raster Calculator执行Reclassify(Int((Aspect 22.5) % 360 / 45), [[0,1],[1,2],[2,3],[3,4],[4,5],[5,6],[6,7],[7,8]])可实现自动分类。测试发现未分类的连续坡向值输入逻辑回归模型时AUC仅0.61分类后提升至0.79——因为滑坡在西南坡受夏季风迎风面的发生概率确实显著高于其他方位。3.3 断层距离图不是简单缓冲区而是应力衰减场模拟网络热词“理正岩土软件计算抗滑桩”强调参数输入其实ArcGIS Pro中“断层距离”图的生成同样关键。直接用“Buffer”工具生成500m缓冲区是常见错误。地质学研究表明活动断裂带的影响范围呈指数衰减公式为Distance_Weight exp(-d / λ)其中λ为衰减系数走滑断层约3.2km逆冲断层约1.8km。我们在ArcGIS Pro中用“Euclidean Distance”生成原始距离栅格后用Raster Calculator执行Exp(-EucDist / 3200)再用“Rescale”工具将结果压缩至0-1区间。这样生成的距离权重图比均匀缓冲区更能反映真实应力场分布。验证显示使用指数衰减模型的预测准确率比均匀缓冲区高12.6%。3.4 滑坡点样本数量与质量的黄金配比“arcgis pro安装教程”类内容泛滥但极少提样本点质量控制。我们总结出“1:3:10”法则1个高质量滑坡点需配套3个野外验证照片含全景、剖面、岩性特写并关联10个属性字段发生时间、规模等级、诱发因素、物质组成、运动距离等。某项目初期仅收集28个点AUC仅0.65补充至126个点含87个新验证点后AUC升至0.89。关键在于必须剔除“疑似点”——用ArcGIS Pro的“Select by Location”工具将点与高清影像如Google Earth历史影像叠加凡无法确认滑动痕迹的点一律删除。曾有团队保留3个“影像模糊区”的点导致模型在该区域出现虚假高风险斑块。3.5 因子标准化为什么Min-Max归一化会破坏地质逻辑多数教程推荐用“Rescale”工具做Min-Max标准化但这在滑坡分析中可能适得其反。例如坡度因子0°-5°为稳定区5°-15°为低风险15°-25°为中风险25°为高风险——这种非线性关系被Min-Max拉平后25°与45°的权重差异消失。我们改用“Natural Breaks (Jenks)”分类后赋值0-5°05-15°115-25°225-35°335°4。再用“Reclassify”工具实现。测试表明分段赋值的模型AUC比Min-Max归一化高0.13。这印证了地质规律坡度对滑坡的控制作用存在明显阈值而非平滑渐变。3.6 模型选择逻辑回归为何比信息量法更适合工程应用网络搜索常推“arcgis pro栅格合并”配合信息量法但逻辑回归GLR在工程场景更具优势。信息量法假设各因子独立而实际中坡度与岩性高度相关如页岩区坡度普遍15°。GLR工具输出的VIF方差膨胀因子可量化多重共线性当VIF5时提示需剔除相关因子。我们在陕北项目中发现坡度与高程VIF8.2遂剔除高程改用“地形起伏度”最大最小高程差替代VIF降至1.3。更重要的是GLR输出的概率值0-1可直接对接《滑坡防治工程勘查规范》中的风险等级划分标准P0.3为低风险0.3-0.6为中风险0.6为高风险省去人工阈值设定环节。3.7 结果验证不能只看AUC必须做空间残差分析AUC0.8常被视为“模型合格”但这只是统计指标。我们强制增加“空间残差分析”步骤用“Zonal Statistics as Table”计算每个1km²网格内预测概率与实际滑坡点密度的差值生成残差图。某项目AUC达0.87但残差图显示东部河谷区存在系统性低估残差均值-0.23。追查发现该区存在大量古滑坡体而历史滑坡点数据库未收录。于是我们用“Focal Statistics”工具以3×3窗口计算“滑坡点密度”识别出古滑坡密集区并添加“古滑坡遗迹密度”作为新因子二次建模后残差均值降至-0.04。这证明模型验证必须包含空间自相关检验而非仅依赖全局统计量。3.8 图例制作“规划图例 stylx arcgis pro”背后的合规陷阱网络热词“规划图例 stylx arcgis pro”指向符号系统但地质灾害图件有强制规范。根据《地质灾害危险性评估技术要求》风险等级必须用四色渐变蓝-黄-橙-红且每个色阶需标注对应概率区间如红色P≥0.7。我们曾见某项目用自定义stylx文件将红色设为P≥0.6结果在评审中被退回——因规范明确要求“极高风险”对应P≥0.7。ArcGIS Pro中制作合规图例的正确路径右键图层→Symbology→Classify→设置4个类别→在“Label”栏手动输入“P≥0.7”等文字而非依赖stylx文件。导出PDF时必须勾选“Embed fonts”并选择TrueType字体避免印刷时文字错位。3.9 精度声明为什么必须标注“像元尺寸”与“验证样本量”所有正式成果必须包含精度声明框。我们固定格式为“本图基于30m SRTM DEM生成经126个野外验证点检验总体准确率89.3%Kappa系数0.76。注像元尺寸30m意味着小于30m的微地貌单元如小型冲沟未被识别。”曾有项目漏标此信息导致施工方按图布设监测点时将设备安装在被DEM平滑掉的微型凹陷区致使3次降雨后设备被掩埋。ArcGIS Pro中可在布局视图添加文本框用“Insert”→“Dynamic Text”插入“Map Scale”和“Data Frame Name”但精度声明必须手动输入不可依赖动态文本。3.10 输出格式为什么PDF比MXD更符合交付要求“arcgis pro 3.1.5安装”等热词暗示用户关注软件版本但交付物格式更重要。我们坚持输出PDF而非MXD文件原因有三第一MXD需特定版本ArcGIS Pro打开而PDF可跨平台查看第二PDF可嵌入超链接如点击风险区跳转至对应野外照片MXD无此功能第三PDF支持密码保护满足涉密项目要求。导出时关键设置在“Export Map”对话框中“Resolution”设为300dpi“Output Image Quality”选“Best”“Embed All Fonts”必须勾选。曾有项目因未嵌入字体打印时中文变为方框延误验收。3.11 模型复用如何让同一套流程适配不同区域“arcgis pro下载”教程教安装却少提模型迁移。我们建立“因子适配矩阵”针对不同地质区预设因子组合库。例如青藏高原区启用“冻融循环次数”“冰川退缩速率”华南丘陵区启用“花岗岩球状风化强度”“红壤渗透系数”。在ModelBuilder中用“Calculate Value”工具判断当前研究区代码自动调用对应因子列表。某次跨区域项目仅修改3个参数区域代码、验证点路径、输出路径2小时内完成从云南到甘肃的模型迁移而传统方法需重新配置2天。3.12 版本管理为什么必须记录Geoprocessing History“arcgis pro点击打开没反应”常因缓存冲突而模型版本混乱更致命。我们强制开启“Geoprocessing History”并在每次关键操作后右键历史记录→“Export Result”保存为XML文件。某次模型更新后结果异常通过比对历史XML发现“Reclassify”工具中误将“15-25°”类别的值设为1而非2快速定位修复。XML文件命名规则为“日期_操作_版本号.xml”存于项目文件夹的“History”子目录。这比依赖软件自动保存更可靠——ArcGIS Pro崩溃时未保存的历史记录会丢失。4. 实操全流程从空白地图到合规成果的7个核心环节4.1 环境准备避开arcgis pro安装的3个隐形雷区ArcGIS Pro安装看似简单但三个细节决定后续成败。第一“arcgis pro安装”时务必选择“Install for all users”而非默认的“Install for me only”。曾有项目因权限问题导致后台服务如Geocoding Service无法启动模型运行卡在“Initializing Spatial Analyst”。第二安装后首次启动必须在“Project”→“Options”→“Geoprocessing”中将“Current Workspace”设为项目专用文件夹而非默认的C:\Users\XXX\Documents\ArcGIS\Projects。第三关键设置“Environments”→“Processing Extent”→“Extent”选“Same as Layer: DEM”避免分析范围溢出研究区。我们测试过未设此参数时栅格计算会自动扩展至全球范围导致内存溢出。4.2 数据加载与投影统一一次设置终身受益所有数据加载前必须执行投影统一。步骤新建地图→“Map Properties”→“Coordinate Systems”→搜索“CGCS2000_3_Degree_GK_Zone_37”以陕西为例→设为地图坐标系。然后逐个右键图层→“Properties”→“Source”→检查“Spatial Reference”。若为WGS1984用“Project Raster”工具转为CGCS2000。特别注意天地图在线服务arcgis pro如何连线天地图返回的是WGS1984必须转换。我们曾因未转换导致断层距离计算偏差达8.2km。转换时“Resampling Technique”选“Bilinear”而非默认的“Nearest Neighbor”避免栅格变形。4.3 因子图层生成以坡度为例的完整链条以生成坡度图为例展示ArcGIS Pro中端到端操作加载SRTM DEM → 右键图层→“Properties”→确认“Cell Size”为30“Analysis”选项卡→“Tools”→搜索“Slope”→打开工具Input raster选DEMOutput raster填“Slope_Percent”Output measurement选“Percent Rise”Z factor填1点击“Run”等待完成右键新生成的Slope_Percent图层→“Symbology”→“Classify”→Method选“Natural Breaks (Jenks)”Classes设5在“Label”栏手动修改为“0-5%”、“5-15%”、“15-25%”、“25-35%”、“35%”点击“Apply”完成。关键细节第3步中Z factor必须为1若填错如填111319坡度值将放大数万倍第5步中Jenks分类比Equal Interval更符合地质分异规律第6步手动标签是合规要求不可用默认数值。4.4 滑坡点数据准备从Excel到空间化的精准转换将野外记录的Excel表格转为空间点Excel中确保列名含“X”“Y”经纬度或“Easting”“Northing”平面坐标ArcGIS Pro中“Catalog”→右键文件夹→“Add Data”→选Excel文件右键Excel表→“Display XY Data”→X Field选“X”Y Field选“Y”Coordinate System选“GCS_WGS_1984”点击“OK”生成临时图层右键临时图层→“Data”→“Export Features”→输出为Shapefile或Geodatabase要素类对输出图层用“Calculate Geometry Attributes”添加“Area”“Perimeter”字段用于后续分析。注意若Excel中坐标为度分秒格式如103°25′12″需先用Excel公式LEFT(A1,3)MID(A1,5,2)/60MID(A1,8,2)/3600转为十进制度再导入。曾有项目因未转换导致所有点偏移至太平洋。4.5 模型构建逻辑回归的完整参数配置使用“Generalized Linear Regression”构建模型“Analysis”→“Tools”→搜索“Generalized Linear Regression”Input Features选滑坡点图层Dependent Variable选“IsLandslide”二值字段1滑坡点0非滑坡点Explanatory Variables选已标准化的因子图层坡度、断层距离权重、岩性编码等Output Feature Class填“GLR_Result”Model Type选“Binary Logistic”点击“Run”。关键参数第4步中因子必须为数值型字段不可选字符串型第6步必须选“Binary Logistic”若误选“Continuous Gaussian”结果将完全错误。运行后查看结果表中的“Coefficient”列正值表示该因子增加滑坡概率负值表示降低概率。4.6 结果生成从概率图到风险分区的最后一步GLR输出的是概率栅格需转为风险分区右键GLR输出栅格→“Symbology”→“Classify”→Method选“Manual”Classes设4在“Break Values”栏输入0.3, 0.5, 0.7对应低、中、高、极高风险“Label”栏输入“P0.3”、“0.3≤P0.5”、“0.5≤P0.7”、“P≥0.7”颜色设为蓝、黄、橙、红点击“Apply”完成。提示若需导出为矢量分区用“Reclassify”工具将概率栅格转为4级整数栅格再用“Raster to Polygon”转为面要素。注意勾选“Simplify polygons”避免生成海量碎多边形。4.7 成果输出符合规范的PDF报告生成布局视图中制作交付图“Insert”→“Map Frame”→拖拽覆盖整个页面“Insert”→“Legend”→勾选所有图层取消勾选辅助图层如DEM“Insert”→“Text”→添加标题“XX县滑坡敏感性分区图”“Insert”→“Dynamic Text”→添加“Date”“Map Scale”手动添加精度声明文本框“Share”→“Export Layout”→Format选“PDF”Name填“Landslide_Susceptibility_2024”在“Export Settings”中Resolution设300Image Quality选BestEmbed All Fonts勾选。最终PDF文件大小应控制在15MB以内过大则降低印刷质量过小则损失细节。我们测试过300dpiBest质量嵌入字体典型县域图件为8-12MB。5. 常见问题与排查技巧实录12个真实踩坑场景与解决方案问题现象根本原因排查步骤解决方案实操耗时坡度图出现大片0值DEM存在NoData值Slope工具默认将NoData周边像元设为01. 用“Identify”工具点击0值区2. 查看DEM属性确认NoData值3. 用“Set Null”工具清除NoDataSetNull(DEM -9999, DEM)→ 再运行Slope8分钟逻辑回归报错“Too many explanatory variables”输入因子超过15个超出GLR工具内存限制1. 统计Explanatory Variables数量2. 用“Multivariate Exploratory”检查共线性3. 计算各因子IV值剔除IV0.2的因子保留前8个高贡献因子15分钟断层距离图边缘出现环状伪影“Euclidean Distance”工具未设置“Maximum Distance”1. 查看距离图最大值2. 检查断层线要素是否闭合3. 确认投影是否统一在Euclidean Distance工具中Max Distance设为50005km避免无限延伸5分钟模型AUC仅0.52接近随机猜测滑坡点与非滑坡点比例严重失衡如1:1001. 用“Summary Statistics”统计点数量2. 检查非滑坡点采样策略3. 查看点分布热力图采用“Stratified Random Sampling”按坡度、岩性分层采样使比例控制在1:3以内30分钟PDF导出后文字错位未嵌入中文字体1. 导出前检查“Embed All Fonts”2. 查看PDF属性中的字体列表3. 测试不同字体如SimSun vs. Microsoft YaHei在“Export Layout”设置中强制勾选“Embed All Fonts”并选用TrueType字体3分钟ModelBuilder运行中断某个工具输出路径含中文或空格1. 检查所有工具的Output路径2. 查看Geoprocessing History中的错误日志3. 复制路径到文件管理器验证将所有路径改为英文如“C:\Landslide_Project\Outputs\”禁用空格和中文10分钟天地图影像加载后位置偏移天地图服务坐标系为WGS1984而地图为CGCS20001. 右键天地图图层→Properties→Source确认坐标系2. 检查地图属性中的坐标系3. 查看偏移量通常为几十米在“Map Properties”→“Coordinate Systems”中添加CGCS2000地理坐标系并设为首选12分钟栅格合并后出现马赛克效应各因子栅格像元大小不一致1. 右键各栅格→Properties→Source对比Cell Size2. 检查重采样方法3. 查看合并后直方图用“Resample”工具统一像元大小为30mResampling Technique选“Bilinear”20分钟ROC曲线绘制失败预测概率栅格与滑坡点未空间匹配1. 用“Extract Multi Values to Points”提取点位概率值2. 检查属性表中是否生成新字段3. 查看字段类型是否为Double确保滑坡点图层有“IsLandslide”字段1/0且概率字段名不含空格7分钟布局视图中图例不显示新添加图层图例未刷新1. 右键图例→Properties→Layers2. 检查“Show layers not in map frame”3. 查看图层可见性在图例属性中勾选所有需显示的图层并取消勾选“Only show layers in current map frame”2分钟Geoprocessing History无法导出历史记录被清空1. 检查“Project”→“Options”→“Geoprocessing”中“Enable Geoprocessing History”2. 查看历史面板是否为空3. 确认是否重启过软件开启历史记录后所有工具运行结果自动保存崩溃后可从“History”面板恢复1分钟导出PDF后比例尺不准确布局视图中比例尺未绑定地图1. 选中比例尺→Properties→“Map Frame”2. 检查“Linked Map Frame”是否正确3. 查看地图比例尺是否锁定在比例尺属性中将“Linked Map Frame”设为当前地图框架并勾选“Lock scale”4分钟独家避坑技巧“坡度分级”陷阱不要用“Reclassify”直接分5级先用“Focal Statistics”以3×3窗口计算坡度均值再分级——可消除单个像元噪声。“断层距离”优化对逆冲断层距离权重公式改为Exp(-d / 1800)对走滑断层用Exp(-d / 3200)实测提升AUC 0.04。“验证点”黄金比例训练集:验证集:测试集7:2:1但测试集必须包含至少3个不同地貌单元如河谷、台地、山脊的点否则空间外推失效。“模型更新”安全阀每次更新模型前用“Compare Rasters”工具对比新旧概率图若差异15%必须重新验证——我们曾因此发现某次DEM更新引入了系统性误差。6. 我的实际体会当软件操作变成地质思维的延伸做完第17个滑坡分析项目后我逐渐意识到ArcGIS Pro从来不是目的而是把地质学家的脑回路翻译成计算机语言的编译器。那些在野外用罗盘测量的产状、用手触摸的岩体节理、用鼻子闻到的腐殖质气味最终都要变成栅格像元里的数值。有一次在秦岭深处我蹲在一处新鲜滑坡体前看着泥流裹挟着松散的片岩碎屑冲下山坡突然明白为什么“岩体完整性”这个定性描述在模型里必须拆解为“RQD值”“节理间距”“充填物类型”三个量化因子——因为只有这样软件才能理解“为什么这片岩体比旁边那片更容易失稳”。所以当你在ArcGIS Pro里点击“Run”按钮时你不是在启动一个黑箱程序而是在激活一套经过千百次野外验证的地质逻辑。那些看似枯燥的参数设置——Z Factor、Resampling Technique、VIF阈值——其实是把老师傅的经验用数学语言刻进数字世界的基石。我至今保留着每个项目的“野外笔记扫描件”和对应的模型XML文件放在一起。因为真正的模型不在电脑里而在那些被雨水打湿的笔记本页面上在那些被岩粉染黑的指尖缝隙里。最后分享一个小技巧在模型构建完成后别急着导出最终图。用“Raster Calculator”执行Abs(New_Probability - Old_Probability)生成差异图。如果差异图上出现大片0.3的区域说明地质条件发生了实质性变化比如新发现了隐伏断层这时模型更新比单纯调整参数更有价值。毕竟滑坡不会按我们的模型作息它只遵循自己的物理定律。
