打开一景遥感影像,看到的是漂亮的色块,但真正决定它能不能量面积、判岩性的,是从辐射值到反射率、再到地质信息的这几步转换有没有做对。这篇文章把链条拆成五段——成像、定标、校正、提取、验证,给出四种分辨率的选型判据、预处理的固定顺序,以及哪些地方遥感确实到不了。适用于区域尺度的岩性识别、构造解译与地表变化监测;不适用于厘米级形变测量和植被穿透探测,这两类要另配雷达干涉或地面手段。
一、从光子到数字:卫星如何把地表"读"成数据
一颗在七百多公里高空运行的卫星,一天能回传上百景遥感卫星数据。这些数据记录了地表每一平方米的反射强度,但它既不是颜色,也不是岩石名称。从原始数值到一张能勾出岩性界线的图件,中间隔着定标、校正、提取、验证四道工序,省掉任何一道,图上的色块就只是色块。
太阳辐射穿过大气时,水汽和气溶胶先散射、吸收掉一部分,剩下的照到地面。不同矿物对各波段光的吸收和反射并不一样:铁氧化物在 0.85 至 0.9 微米附近有明显吸收,含羟基的黏土矿物在 2.2 微米附近有一道特征谷(据公开资料,常见矿物光谱库)。这些"指纹"才是遥感能区分岩性的物理依据,与画面好不好看没有关系。
传感器把每个波段接收到的辐射量量化成整数,8 位对应 0 至 255,12 位对应 0 至 4095。这一步得到的是 DN 值,一个相对亮度。同一片蚀变带,夏天正午成像和冬天下午成像,DN 值能差出一大截,因为太阳高度角变了,穿过大气的路径长度也变了。
所以第一件事要认清:图上的亮不等于反射率高,图上的红也不等于岩石发红。多数时候那只是三个波段被分别赋成红绿蓝的结果,属于制图习惯,不是物理结论。

二、四个分辨率:选数据之前先把这四个词分开
提到分辨率,多数人只想到"能看清多大的东西"。遥感里的四个分辨率经常被混着用,选数据时混一个,后面全要返工。
表1 四种分辨率的含义、常见量级与代价
| 分辨率类型 | 含义 | 常见量级(据公开资料) | 地质上的直接作用 | 代价 |
|---|---|---|---|---|
| 空间分辨率 | 单个像元对应的地面尺寸 | 陆地卫星 8 号多光谱 30 米;哨兵 2 号 10 至 60 米;高分二号 0.8 米 | 决定最小可识别地质体的尺度 | 空间分辨率越高,幅宽越窄,单景价格越高 |
| 光谱分辨率 | 波段数量与各波段带宽 | 多光谱 4 至 13 个波段;高光谱 100 个以上 | 决定能否区分矿物种类而非只分大类 | 波段越多,数据量与处理复杂度同步上升 |
| 时间分辨率 | 同一地点重复观测的间隔 | 陆地卫星 8 号约 16 天;哨兵 2 号约 5 天 | 决定变化监测能力与去云补空的余地 | 高频重访通常以牺牲空间分辨率为代价 |
| 辐射分辨率 | 量化的灰度级数 | 8 位 256 级;12 位 4096 级 | 决定弱差异地物能否被分开 | 高位深需要配套定标产品与更大存储 |
判据很干脆:目标地质体的边长至少要覆盖 3 至 5 个像元,边界才能勾得稳;低于这个尺度,看到的是混合像元,不是一个地质体【行业测算值】。算面积同理——一个 30 米像元是 900 平方米,边界上判错一个像元,就是 900 平方米的误差,这个量级在矿区圈定里不能忽略。

三、一份影像的预处理链条
拿到一景原始遥感卫星数据,工序顺序基本固定,颠倒任何一步都要从头再来。
第一步是辐射定标,用元数据里的定标系数把 DN 值换算成辐亮度或表观反射率。这一步不能跳,跳过去后面的所有比较都失去基准。第二步是大气校正,把表观反射率还原成地表反射率,去掉气溶胶和水汽的影响;跨时相做对比、做定量反演必须做,只做单时相制图可以省,但镶嵌时接缝会明显偏色。第三步是几何校正与正射纠正,用有理多项式系数配合数字高程模型处理,平地区可以放宽,山区必须做,否则山坡上的线性体会整体偏移,位置误差随坡度放大。
后面三步属于精加工。影像融合把全色波段的空间细节注入多光谱,但算法会改动光谱值,做蚀变提取最好回到融合前的原始多波段。镶嵌与匀色处理多景拼接,先做直方图匹配再拼,否则接缝线一眼可辨。最后按工作区裁剪,水体、云、阴影单独做掩膜,别让它们混进统计。
一句判据记住:只要涉及"比较"——比时间、比区域、比波段,就必须走到地表反射率这一层;只涉及"看",表观反射率够用。产品级别怎么选、哪些级别自带定标参数,可以对着[多光谱遥感数据产品级别说明]({{内链1}})逐条核对。
四、地质场景里的三种典型用法
岩性识别靠两样东西叠加:光谱特征和纹理。厚层灰岩纹理均一,花岗岩常见节理网状纹理,单靠颜色分不开的单元,用纹理和地形阴影补。判据是同一套光谱曲线在两类岩石上出现可分差异,且野外能对上,才算识别成立。
蚀变异常提取常用波段比值和主成分变换。铁染类异常多用短波红外与可见光的比值组合,羟基类异常用 2.2 微米附近的波段比值,具体波段编号随传感器不同要查对应产品的波段说明,不能照搬。配合主成分变换可以压掉植被干扰,但阈值定在哪个分位数,必须在工作区里用已知矿点反推,不能套用别的区域的经验值。
构造解译看线性体、环形构造和节理密集带。在数字高程模型的山影图上叠多方向光源,比只看影像更容易勾出隐伏断裂的走向;把解译出的线性体与区域地质图对照,方向一致才算可信,不一致要回头查是位置偏移还是判读错误。同一套方法在更大范围怎么铺开,可参照[区域地质调查中的遥感解译流程]({{内链2}})。
五、边界、风险与落地建议
遥感到不了的地方
云是最直接的限制。南方多云地区一年内有效观测天数有限,单景成图往往要靠多时相合成补空【行业测算值】。植被和厚层风化壳是第二层限制,地表被覆盖之后,下面的岩性信号传不上来,红土覆盖区的识别精度会明显下降。混合像元的问题在边界带和小规模地质体上最突出,一个像元里混了几种地物,反演出来的是平均值而不是任何一类。另外,不同传感器、不同版本的产品,反射率不能直接混用,定标系数一改,历史数据就要重处理。最后一条最容易被忽略:没有剖面、没有采样点,精度只能定性地说"看着像",给不出可核查的准确率。
一套能跑起来的最小流程
先定目标尺度,再按 3 至 5 个像元的判据反推需要什么级别的遥感卫星数据,不要先选数据再想用途。做完定标和大气校正,只提取一个要素,别一次上五六个。出图前用野外点或已有地质资料抽验,抽验比例按图斑数取一成到两成。最后出图时必须写明误差范围和数据时相,图件没有时相标注,等于不可复核。
判据是样本量:验证样本少于 30 个,不要给百分比精度,只写"定性一致"。
要点梳理
- DN 值不是反射率,跨时相比较必须做大气校正。
- 目标体至少覆盖 3 至 5 个像元,否则看到的是混合像元。
- 融合会改动光谱,蚀变提取用融合前的原始多波段。
- 山区必须正射纠正,否则线性体位置随坡度整体偏移。
- 蚀变异常是概率不是矿,没有地面验证就没有精度。
- 先跑通一个要素,再考虑叠加第二个。
六、结语
一颗卫星做的事,本质是把电磁波换成数字,剩下的判断全在地面这一端。遥感能提供的是空间上的连续性和尺度的整体感,给不出成分的确定值。把它当成一张可以量、可以比、可以反复核对的底图,比当成一张找矿图更贴近它的能力边界。真要落到图件上,先问清三个问题:目标体多大、要比较什么、有没有验证点。这三问答完,数据选型和流程基本就定死了。