1. 课程导入与实战目标
构造线性体是地质构造在遥感影像与地形数据上最直观的表达形式之一。断裂、节理密集带、褶皱轴迹、岩性界线以及韧性剪切带，往往在地表留下线性延伸的负地形、色调异常或纹理突变。在DEM（数字高程模型）上，这些构造表现为坡向突变、坡度脊线、地形剖面拐点与局部高程错断。把它们从连续的高程场中"剥离"出来，就是构造线性体提取的核心任务。它既是区域地质填图的辅助手段，也是活动构造调查、水文地质通道分析、工程选址与矿产预测的重要基础数据。
与目视解译相比，基于DEM的自动或半自动提取具有可重复、可量化、可批量处理的优势。但必须清醒认识到：算法提取的是"地形线性特征"，而非"地质构造"本身。一条提取出的线性体可能是断层，也可能是岩脉、冲沟、人工道路或地层界线。因此本课程强调"提取—筛选—地质验证"的完整链条，而不是把算法输出直接当作构造图。
为什么用DEM而不是只用影像
光学影像依赖光照与色调，阴影方向一变，线性体的可辨识度就大幅波动；植被覆盖区更是难以判读。DEM记录的是几何信息，与光照无关，能够稳定地反映地形起伏。尤其在植被茂密、影像纹理弱的地区，基于坡向、曲率、剖面形态的线性体增强往往比影像更可靠。实际工作中，二者应当互补：DEM负责几何骨架，影像负责岩性与色调佐证。
DEM数据源的选择
数据源直接决定可提取线性体的尺度下限。常用的几类数据源对比如下（参数为公开典型值，具体以实际数据说明为准）：
| 数据源 | 典型分辨率 | 高程精度（示例） | 适用尺度 | 主要限制 |
|---|---|---|---|---|
| SRTM 1弧秒 | 约30 m | 垂直约±16 m | 区域构造格架 | 植被/建筑偏差，无植被穿透 |
| ASTER GDEM | 约30 m | 垂直约±20 m | 区域构造格架 | 噪声较大，需去噪 |
| ALOS PALSAR DEM | 约12.5 m | 垂直约±5 m | 区带—局部构造 | 雷达侧视阴影与叠掩 |
| 机载LiDAR DEM | 0.5～2 m | 垂直约±0.1～0.3 m | 局部精细构造、节理带 | 成本高，覆盖有限 |
选择原则可以概括为"尺度匹配"：研究目标若是数十至数百千米的区域断裂格架，30 m级数据足够；若要识别数百米级的次级断裂或节理密集带，应优先选用12.5 m或更高分辨率数据。分辨率并非越高越好——过高分辨率会把微地貌噪声、树冠、田埂一并放大，反而干扰构造识别。此外还需关注数据是否经过"去植被"处理（DSM与DTM的差别），DSM中的树冠会形成大量伪线性体。
完整工作流概览
本课程贯穿的工作流可归纳为六个环节，下图给出了各环节的输入、处理与输出关系，以及它们之间的反馈路径：
需要特别说明的是图中那条虚线反馈：线性体提取很少一次成功。当统计结果显示方向杂乱、密度异常或与已知构造完全不符时，应回到地形增强环节调整滤波方向与窗口，甚至更换数据源，而不是在提取环节反复调阈值"凑结果"。
实战任务说明
本课程的实战任务设定为：以一幅覆盖典型断裂发育区的DEM为输入，完成从预处理到构造线性体图输出的全流程，并提交三项成果——线性体矢量图层、走向玫瑰图与密度等值图、以及一份简要的地质合理性说明。示例区可选择一条已知区域性断裂通过的山区，长度约40～60 km，以便用公开地质图做交叉验证。
为了让任务可执行、可检查，建议在开始前明确以下验收要点：
- 数据一致性：全流程使用同一投影坐标系（建议UTM或高斯-克吕格），像元大小与分辨率匹配，避免重采样引入虚假线性。
- 参数可追溯：记录滤波方向、窗口大小、长度阈值、最小间距等关键参数，保证结果可复现。
- 方向合理性：提取结果的主走向应与区域构造应力场或已知断裂走向大体一致，允许存在次级方向，但需能解释。
- 地貌响应：每条保留的线性体应能在山体阴影或剖面图上找到对应的地形证据，而非纯算法伪影。