DEM制作和表面分析

DEM制作和表面分析地理信息系统核心技术：数字高程模型构建与地形表面分析Contents课程目录数字高程模型与空间分析核心知识体系概览01DEM基本概念与理论02TIN模型与DEM制作方法03空间表面分析技术04空间内插方法与算法CHAPTER01DEM基本概念与理论数字高程模型的定义、分类与应用场景DIGITALTERRAIN什么是DEM数字高程模型（DEM）是用有序数值阵列表示地面高程的实体模型，作为DTM的核心分支，专注于地形高程信息的数字化表达，是现代地理信息系统中不可或缺的基础数据。DEM三维地形可视化渲染01DEM全称为DigitalElevationModel，是用一组有序数值阵列形式表示地面高程的实体地面模型02DEM是数字地形模型（DTM）的一个分支，DTM描述高程、坡度、坡向、坡度变化率等多种地貌因子03核心功能是将连续地形表面离散化为规则或不规则的高程点集合，便于计算机处理和分析04作为4D产品（DOM、DEM、DRG、DLG）之一，是目前最主要的基础地理信息产品之一DATAMODELDEM的两种数据模型DEM采用Grid（栅格）和TIN（不规则三角网）两种数据模型，前者结构简单高效适合宏观分析，后者精确灵活适合精细建模，两者互补构成完整的地形表达体系。Grid（栅格）模型—规则矩形网格结构Grid（栅格）模型将地形划分为规则的矩形网格，每个网格单元存储一个高程值，数据结构简单统一计算效率高，存储需求相对较小，适合区域性、小比例尺的宏观地形分析小比例尺在平坦区域可能造成数据冗余，在复杂地形区域可能丢失细节信息TIN（不规则三角网）模型将空间上不规则分布的点连接成三角形网络，能更精确地表达复杂地形特征数据量大、计算复杂、耗内存多，常用于精细的、大比例尺的地形建模大比例尺在地形变化剧烈区域能自适应加密节点，在平坦区域可减少节点数量TIN（不规则三角网）模型—自适应三角网络结构SpatialDataModelTIN模型详解TIN（不规则三角网）通过将空间不规则分布的点连接成三角形网络来表达连续表面，作为特殊的矢量数据格式，它在地形精细建模中具有不可替代的优势。01TIN将具有(X,Y,Z)坐标值且在空间分布上不规则的点连接成三角形，相邻三角形形成网络表示地形特征02尽管TIN属于矢量数据格式，但由于不能表达单个要素，因此是一种特殊的矢量格式03TIN用于表示连续的表面数据（如地形高程），而不能表示离散的数据对象04TIN能够自适应地形复杂度，在陡峭区域加密节点、在平坦区域稀疏节点，实现数据的最优分布05TIN模型在大比例尺精细地形建模、工程测量、三维可视化等领域有广泛应用TIN三角网地形建模渲染效果SPECIFICATIONDEM坐标与尺寸规定DEM制作遵循严格的坐标排列规则和格网尺寸标准，通过高程中误差等精度指标确保数据质量，为不同比例尺应用提供规范化的技术保障。01DEM格网坐标按由西向东、由北向南排列，左上角第一个格网点对应高斯平面坐标系的起始坐标021:500至1:2000比例尺DEM格网尺寸不大于千分之一比例尺分母，1:5000至1:10万不大于万分之五比例尺分母03DEM成果按精度分为三级，格网高程精确到小数后两位，用格网点的高程中误差表示04高程中误差的2倍为采样点数据的最大误差，内插点高程中误差是格网点高程中误差的1.2倍不同比例尺DEM格网尺寸标准比例尺格网尺寸要求典型应用1:500~1:2000≤1/1000比例尺分母工程设计、城市规划1:5000~1:10万≤5/10000比例尺分母区域规划、资源调查1:25万~1:100万根据精度等级确定宏观分析、科研教学不同比例尺DEM对应不同的格网尺寸要求，比例尺越大，格网尺寸越小，精度越高APPLICATIONSDEM的主要应用领域DEM作为基础地理信息产品，在工程测量、城市规划、资源管理、防灾减灾等领域发挥着关键作用，是现代空间分析不可或缺的核心数据源。工程测量与规划查询单点高程、计算土石方量、地表面积计算、绘制断面图等工程测量基础应用与DOM结合制作三维景观图，支持城市规划、建筑设计、道路选线等可视化需求Engineering&Planning资源环境管理水资源分析、流域划分、洪水模拟，支持水利工程建设与防洪减灾决策土壤侵蚀评估、植被分布研究、生态环境监测，为资源可持续利用提供数据支撑ResourceManagement防灾减灾滑坡风险评估、泥石流模拟、地质灾害预警，保障人民生命财产安全地形稳定性分析、边坡工程评估，支持矿山开采与基础设施建设安全管理DisasterPreventionCHAPTER02TIN模型与DEM制作方法从数据采集到DEM生成的完整技术流程TECHNIQUESDEM主要制作方法DEM制作方法涵盖航空摄影、LiDAR、地形图矢量化、野外测量等多种技术路线，不同方法在精度、效率、成本上各有优劣，需根据实际需求选择。航空摄影测量专业设备01航空摄影法：通过航空影像进行立体测图，量测等高线并内插得到DEM，是传统主流方法02雷达/LiDAR：快速获取大范围高精度点云数据，通过滤波和内插生成DEM03地形图矢量化：将纸质地形图数字化，经扫描、配准、矢量化提取高程信息04全野外数据采集：使用全站仪、GPS实地测量高程点，精度最高但效率较低05DSM滤波法：对数字表面模型滤波去除植被、建筑等非地面要素，得到DEMProductionWorkflowDEM生产基本流程DEM生产遵循从资料准备到成果提交的标准化流程，涵盖定向建模、特征采集、TIN构建、数据编辑、接边镶嵌等关键环节，确保数据质量和规范合规。DEM生产流程各阶段要点阶段主要任务输出成果01资料准备收集控制点、影像、地形图基础数据集02定向建模建立影像内外方位元素定向参数03特征采集提取特征点、等高线、断裂线矢量特征数据04TIN构建构建TIN并内插生成DEMTIN和DEM数据05数据编辑修正错误、填补空洞编辑后DEM06接边镶嵌图幅拼接、边缘处理完整区域DEM07裁切输出按需裁切、格式转换最终成果数据08质量检查精度检验、完整性检查质检报告与成果DEM生产流程包含8个主要阶段，从资料准备到质量检查，每个阶段都有明确的任务和输出成果ArcGIS·3DAnalyst创建TIN操作步骤TIN创建需要正确准备等值线和离散点数据，通过ArcGIS的3DAnalyst工具实现，Hardline、Softline、MassPoints等要素类型对TIN质量有重要影响。数据准备01准备等值线数据，确保包含高程字段，用于构建地形的基本骨架02准备离散点数据（MassPoints），必须包含高程值字段，用于补充地形细节03准备Hardline（硬线）数据，用于表示表面不连续性，如悬崖、堤坝、河岸等等值线·离散点操作步骤01加载等值线数据，在属性中正确指定高程字段，确保高程信息被正确读取02打开3DAnalyst扩展模块，选择CreateTIN工具，设置输入要素和参数03指定输出TIN的存储路径和名称，执行工具生成TIN模型3DAnalyst关键概念01Hardline：提供Z值的硬线，在TIN中生成节点，表示表面不连续性02Softline：不提供Z值的软线，表示线性要素但不改变表面坡度，如区域边界03MassPoints：离散高程点，用于补充地形细节，生成的TIN没有线要素Hard·Soft·MassTINOPERATIONSTIN裁剪与高级编辑TIN裁剪和高级编辑功能（HardErase、HardReplace等）是精细化地形建模的关键技术，能够有效处理局部数据更新、错误修正和区域裁剪等复杂需求。01HardErase：删除TIN中指定区域的数据，用于去除不需要的区域或修正错误数据02HardReplace：用新的数据替换TIN中的指定区域，字段值会被用于更新高程信息03连续性处理：两块等高线创建TIN时可能出现不连续问题，需要通过裁剪或编辑确保表面连续性04裁剪边界：裁剪操作需要指定裁剪边界，可以是多边形要素或矩形范围，确保输出TIN符合研究区域05拓扑验证：编辑后的TIN需要重新检查拓扑关系和表面连续性，确保没有悬挂三角形或空洞GIS软件中TIN地形编辑工作场景ArcGIS·SpatialAnalysisTIN转DEM操作TIN转DEM是地形建模的关键步骤，通过ConvertTINtoRaster工具实现，CellSize参数的合理设置直接影响DEM的精度、数据量和后续分析效果。使用ConvertTINtoRaster工具将TIN转换为栅格DEM，便于后续空间分析Elevation字段按原值转换，确保高程信息完整保留CellSize默认值为外接矩形最短边÷250，可按精度需求调整单位取决于坐标系：地理坐标为度，投影坐标为米CellSize越大粒度越大、行列数越少，需权衡精度与数据量CellSize可为浮点数或小于1，适用高精度小范围地形建模CellSize对DEM数据的影响CellSize精度数据量适用场景较小（如1m）高精度大数据量精细工程测量中等（如10m）中等精度中等数据量区域规划分析较大（如30m）较低精度小数据量宏观地形分析CellSize是DEM精度的关键参数，需在精度需求和数据量之间找到平衡RemoteSensing·DEM遥感卫星影像制作DEM遥感卫星影像是现代DEM制作的重要数据源，具有覆盖范围广、获取周期短、成本相对较低的优势，特别适合大范围地形数据的快速获取和更新。常用卫星数据源Landsat系列：影像质量较高，空间分辨率30m，获取周期16天，适合中长期地形监测Sentinel-1：合成孔径雷达卫星，全天候获取能力，空间分辨率5-20m，获取周期6天Sentinel-2：多光谱光学卫星，分辨率10-60m，获取周期5天，适合植被覆盖区域30m/6天影像预处理辐射定标：将卫星影像的辐射量转换为反射率，消除传感器响应差异大气校正：消除大气散射、吸收对影像的影响，恢复地表真实反射特性几何校正：纠正影像的几何变形，确保与实际地形坐标精确对应3步骤高程信息提取像元直接提取：分析像元光谱信息直接提取高程，适用于地形平坦区域地形模型间接提取：利用ASTERGDEM等参考模型配准提取，适用于复杂地形立体像对提取：利用同一地区不同角度的影像进行立体匹配，生成高精度DEM3方法REMOTESENSINGSentinel-2制作DEM实例以Sentinel-2卫星影像为例，从数据获取到DEM生成需经历辐射定标、大气校正、几何校正、高程提取、格式转换等关键步骤，每步都需要严格的质量控制。01获取影像数据从欧空局Copernicus数据中心免费下载Sentinel-2卫星影像数据Copernicus02辐射定标与大气校正使用Sen2Cor等专用工具对影像进行辐射定标和大气校正，消除大气影响Sen2Cor03几何校正利用ASTERGDEM作为参考，对Sentinel-2影像进行几何校正，确保坐标精度ASTERGDEM04高程信息提取采用立体匹配或与参考DEM配准的方法，提取像元的高程信息ELEVATION05DEM数据生成将高程信息整理编辑，生成GeoTIFF或ENVI格式的DEM数据GeoTIFF质量控制检查DEM的精度、完整性、连续性，确保与实际地形相符合QA/QCCHAPTER03空间表面分析技术基于DEM的地形因子提取与空间分析方法TERRAINANALYSIS坡度分析坡度是地形分析最基本的因子，表示地表倾斜程度，通过DEM计算获得，在农业、工程建设、城市规划等领域有广泛应用。坡度分级与应用场景不同坡度等级对应不同的土地利用方式和开发限制，坡度分析是土地规划的重要依据坡度等级坡度范围适宜用途限制因素平地0–5°城市建设、农田耕作排水需人工设计缓坡5–15°果园、梯田、低密度建设需考虑水土保持斜坡15–25°林地、牧业、景观开发建设成本高、水土流失风险陡坡25°以上生态保护、水源涵养禁止开发建设农业应用坡度影响土壤侵蚀程度、灌溉方式选择、机械化作业可行性工程建设坡度决定道路选线、建筑选址、边坡稳定性评估、土石方量计算城市规划坡度影响用地分类、开发强度、排水系统设计、绿地布局TerrainAnalysis坡向分析坡向表示坡面朝向，以0-360度方位角表示，对光照、温度、植被分布有重要影响，是生态研究、农业规划和建筑设计中不可或缺的地形分析因子。方位角系统用0–360°表示坡面朝向：0°正北、90°正东、180°正南、270°正西0°–360°ArcGIS计算Aspect工具从DEM计算坡向，基于每个像元与相邻像元的高程变化方向Aspect坡向分析地形渲染图阳坡特征南坡与西南坡光照充足、温度高、蒸发量大，适合喜阳植物生长南坡·西南坡阴坡特征北坡与东北坡光照少、温度低、湿度大，适合耐阴植物生长北坡·东北坡应用领域影响农业作物选择与灌溉，以及建筑采光设计与能耗评估农业·建筑TerrainAnalysis山体阴影分析山体阴影（Hillshade）通过模拟太阳光照射地形的阴影效果，将二维DEM转化为具有三维立体感的可视化图像，是地形展示和地图制图中最常用的晕渲方法。模拟太阳光照射地形产生的阴影效果，让二维地图呈现三维立体感ArcGISHillshade工具需设置方位角（Azimuth）和高度角（Altitude）两个关键参数方位角表示光源水平方向，通常设为315°（西北方向），产生最佳视觉效果高度角表示光源垂直角度，通常设为45°，平衡阴影深度和细节可见性不包含具体高程数值，但能直观展示地形起伏、山脊、山谷等地貌特征常用于地图底图、与坡度坡向叠加显示、三维场景背景，增强地形可视化效果地形晕渲效果示意·Hillshade渲染TERRAINANALYSIS剖面分析剖面分析从DEM中提取沿特定路径的地形起伏信息，生成高程变化曲线，是道路设计、管线铺设、地质勘探等工程设计中不可或缺的地形分析工具。剖面分析在地形设计中的典型应用从DEM中提取沿特定路径的地形起伏信息，生成高程变化曲线ArcGIS工具：3DAnalyst/Surface生成剖面图，支持直线或折线路径道路设计：计算填挖方量、设计纵坡、评估视距、优化路线管线铺设：确定管道坡度、计算泵站位置、评估施工难度地质勘探：分析地层结构、评估地质条件、规划钻探位置水利工程：河道纵坡分析、水库库容计算、堤防设计SPATIALANALYSIS可视性分析可视性分析确定观察点的可见区域或两点间的通视性，包括视域分析和通视分析两种类型，在通信选址、军事侦察、景观规划等领域有重要应用价值。视域分析Viewshed确定从观察点可以看到的区域，输出栅格图层表示每个像元是否可见考虑地球曲率和大气折射影响，可设置观察点高度和目标高度参数通信基站选址、瞭望塔布局、风电场视觉影响评估、观景点规划栅格输出通视分析LOS确定两点之间是否相互可见，输出线段标注可见与不可见部分考虑地形起伏对视线遮挡，可设置观察点和目标点的高度偏移军事侦察、射击阵地选择、输电线路设计、景观视线廊道分析线段标注关键参数设置观察点高度—观察点相对于地面的高度，如塔高、身高目标高度—目标对象相对于地面的高度，如天线、建筑高度分析半径—视域分析的最大距离，超出范围不进行分析OFFSETHYDROLOGY·GIS流域分析流域分析是基于DEM的水文分析核心方法，通过填洼、流向、流量、流域划分等步骤提取河网和流域边界，是水资源管理、洪水预报、生态评估的重要基础。流域分析主要步骤与工具步骤ArcGIS工具输入数据输出结果1.填洼Fill原始DEM无洼陷DEM2.流向FlowDirection无洼陷DEM流向栅格3.流量FlowAccumulation流向栅格流量栅格4.河网RasterCalculator流量栅格+阈值河网栅格5.流域Watershed流向栅格+出水点流域多边形NOTE五步流程形成完整的水文分析链路，广泛应用于水资源管理、洪水预报、水质监测、生态评估与水利工程设计。DEMAnalysis等高线生成从DEM生成等高线是地形图制作的基本操作，通过设置等高距和起始高程，将栅格DEM转化为矢量等高线，既可用于地图制图，也可用于DEM质量检验。CONTOURTOOLArcGIS的Contour工具从DEM生成等高线，需要设置等高距（ContourInterval）和起始等高线SCALE等高距选择取决于地形复杂度和地图比例尺：大比例尺用小等高距，小比例尺用大等高距VECTOR生成的等高线是矢量线要素，每个线要素包含高程属性，可用于地形图制图和三维建模INDEXLINE计曲线（IndexContour）：每隔4-5条基本等高线加粗一条，便于读图和计算高程SMOOTH等高线平滑：使用SmoothLine工具对生成的等高线进行平滑处理，提高地图美观度QACHECKDEM质量检验：通过检查等高线的平滑度、连续性、闭合性，发现DEM中的错误和异常Chapter04空间内插方法与算法DEM制作中核心内插算法的原理与应用SPATIALINTERPOLATION空间内插基本概念空间内插是根据已知点数据估算未知点数据的方法，基于地理学第一定律，分为确定性方法和随机性方法两大类，是DEM制作从离散点到连续表面的核心技术。空间内插曲面拟合三维渲染01根据已知点数据估算未知点，实现从离散点到连续表面的转换02基于地理学第一定律：空间上相近的事物比相距远的更相似03确定性方法：反距离权重法（IDW）、样条函数法（Spline），直接计算未知点04随机性方法：克里金法（Kriging），考虑空间自相关并提供不确定性度量05DEM制作：通过等高线和高程点建立TIN，再基于内插算法生成DEM06方法选择需考虑数据密度、地形复杂度、精度要求与计算效率SPATIALINTERPOLATION线性与双线性内插线性内插和双线性多项式内插是基础的空间内插方法，原理简单、计算效率高，适合大规模数据处理，但在地形复杂区域精度有限，常用于DEM重采样等快速处理场景。线性内插假设未知点的值是已知点的线性组合，在一维情况下是两点间的直线连接二维情况下基于三角形或矩形网格，在网格内部进行线性插值原理简单、计算速度快，但表面不够平滑，可能出现折痕一维/二维双线性多项式内插在矩形网格四个角点之间进行两次线性插值，先沿X方向，再沿Y方向数学表达式Z=a+bX+cY+dXY，包含四个待定系数比线性内插更平滑，计算仍较简单，但地形复杂区域精度有限Z=a+bX+cY+dXY应用场景DEM重采样：将DEM从一个分辨率转换到另一个分辨率影像配准：在影像几何校正中估算新位置像元的值快速预览：在需要快速处理大规模数据时作为初步估算方法DEM·配准·预览InterpolationMethods双三次卷积内插双三次卷积内插使用16个邻近点进行估算，基于三次样条函数保证表面一阶和二阶导数连续，生成平滑连续的表面，适合高精度DEM制作，但计算量较大。01使用4×4窗口共16个邻近点估算未知值，比双线性内插考虑更多邻近信息02数学基础为三次样条函数，保证表面一阶与二阶导数连续，表面非常平滑03生成平滑连续表面，适合光滑变化地形，视觉效果优于双线性内插04计算量约为双线性内插4倍，数据稀疏区可能出现过冲现象05过冲时估算值超出已知点范围，需后处理裁剪；广泛用于高精度DEM与遥感处理不同内插方法对比内插方法使用点数平滑度计算速度最近邻1个不连续最快双线性4个一阶连续快双三次卷积16个二阶连续中等样条函数全局非常平滑较慢不同内插方法在使用点数、平滑度和计算速度上有显著差异，需根据应用需求选择InterpolationMethods移动拟合法内插移动拟合法在每个待插值点周围选择局部窗口拟合多项式曲面，能够适应地形的局部变化，适合地形复杂度变化大的区域。01在待插值点周围选取局部窗口，用窗口内已知点拟合多项式曲面，估算待插值点的值。02拟合窗口随待插值点移动而移动，能够适应地形的局部变化特征。03通常使用二次多项式：Z=a+bX+cY+dX²+eXY+fY²，含6个待定系数。04至少需要6个已知点进行最小二乘求解，实际应用中通常使用更多点以提高精度。05平坦区域可用大窗口简单拟合，复杂区域可用小窗口精细拟合，灵活适应地形。06在DEM制作中广泛应用，特别适合地形复杂度变化较大的区域。移动拟合法地形曲面拟合效果示意SpatialInterpolation反距离权重法（IDW）反距离权重法（IDW）用已知点值的加权平均估算未知点，权重与距离的幂成反比，原理简单、计算快速，但可能产生牛眼效应，适合数据分布均匀、地形平缓的区域。01基本原理：距离待插值点越近的已知点影响越大，距离越远影响越小02加权平均：权重是距离幂函数的倒数：W=1/dp03幂参数：控制距离衰减速度，幂值越大近距离点影响越大，表面越不平滑04优点：原理简单、计算快速、结果直观、易于理解和实现05缺点：可能产生牛眼效应，围绕数据点形成圆形等高线，不考虑地形趋势06适用场景：数据分布均匀、地形变化平缓区域，或快速预览的初步估算IDW反距离权重插值结果示意SPATIALINTERPOLATION克里金法（Kriging）克里金法是一种考虑空间自相关性的高级随机性内插方法，通过变异函数量化空间结构，可提供估算的不确定性度量。01考虑空间自相关性，不仅用距离还用数据之间的空间关系进行估算02通过变异函数量化空间自相关性，描述数据差异随距离变化的规律03可提供估算不确定性（方差），评估可靠性；可考虑趋势和方向性04计算复杂，需建立变异函数模型，对参数设置敏感，需要专业知识05变异函数模型：球面、指数、高斯模型等，描述空间自相关性衰减方式06应用于矿产资源评估、环境污染物分布、土壤属性估算等高精度场景模型类型特征适用场景球面模型存在明确的相关范围大多数自然现象指数模型相关性逐渐衰减连续渐变的属性高斯模型原点处平滑过渡非常平滑的现象线性模型线性衰减无上限趋势明显的现象不同变异函数模型描述不同的空间自相关特征，需根据数据特性选择合适的模型INTERPOLATIONMETHODS样条函数法（Spline）样条函数法用分段多项式拟合已知点，保证表面高阶导数连续，生成非常平滑的表面，能精确通过所有已知点，适合