具体实施方式

下面结合实施例对本发明作进一步详细说明。

本发明一种适用于前寒武纪基底达马拉岩系地层的多源遥感识别方法，包括如下步骤：

步骤一、获取多源卫星遥感数据和DEM高程数据。

获取覆盖某地区地表出露前寒武纪基底达马拉岩系主要地层的多源卫星遥感数据和DEM高程数据；多源卫星遥感数据包括ASTER卫星遥感数据、QuickBird卫星遥感数据。ASTER卫星遥感数据最高空间分辨率为15米，具有可见-短波-热红外15个谱段；QuickBird卫星遥感数据最高空间分辨率为0.61米，具有可见-近红外5个谱段；上述两种遥感数据的采集时间为正午时分，且天空无云、信噪比高；DEM高程数据指空间分辨率为30米的数字高程模型数据。

覆盖某地区地表出露前寒武纪基底达马拉岩系主要地层包括可汗组、罗辛组、楚斯组、卡里毕比组。

步骤二、对步骤一中获取的多源卫星遥感数据进行预处理，获取预处理后的ASTER和QuickBird遥感数据。

对步骤一获取的多源卫星遥感数据，包括ASTER卫星遥感数据和QuickBird卫星遥感数据进行预处理，预处理包括辐射校正、几何校正、数据镶嵌和图像增强，获取预处理后的ASTER和QuickBird遥感数据。采用辐射回归分析法完成辐射校正，采用多项式纠正法完成几何校正，采用基于地理坐标的镶嵌法完成数据镶嵌，采用直方图匹配法完成图像增强等预处理。

步骤三、对步骤二中预处理后的多源卫星遥感数据的ASTER遥感数据进行假彩色合成和三维影像生成。

对步骤二预处理后的多源卫星遥感数据的ASTER遥感数据的第八波段、第六波段和第一波段进行红、蓝、绿彩色变换，并通过对比度拉伸形成假彩色合成影像；以假彩色合成影像为平面数据，步骤一中获取的DEM高程数据为Z轴数据，基于三维地质建模方法—趋势面模型法生成ASTER假彩色三维影像。

步骤四、对步骤二预处理后的多源卫星遥感数据的QuickBird遥感数据进行真彩色融合处理。

对步骤二中预处理后的多源卫星遥感数据的QuickBird遥感数据的第三、第二、第一3个波段进行红、绿、蓝彩色变换，并通过对比度拉伸形成真彩色影像，基于主成分变换法对合成的QuickBird真彩色数据与QuickBird遥感数据全色波段进行融合，获取空间分辨率为0.61米的QuickBird真彩色融合影像。

步骤五、主要地层岩石样品采集和光谱测量。

在前寒武纪基底达马拉岩系工作区采集主要地层岩石样品，样品采集尺寸为3cm(高)×6cm(宽)×9cm(长)，并用GPS记录样品采集地点经纬度格式的点位坐标；光谱测量采用可见-短波地面光谱测量仪，可见-短波地面光谱测量仪使用美国ASD公司生产的FieldSpec PRO FR光谱仪,在室内对岩石样品进行反射光谱测量，获取主要地层岩石样品的光谱曲线。

前寒武纪基底达马拉岩系工作区主要地层包括可汗组、罗辛组、楚斯组、卡里毕比组。

步骤六、主要地层岩石样品光谱预处理和诊断标志识别，获取前寒武纪基底达马拉岩系主要地层岩石的光谱诊断吸收峰对应的波长位置。

采用基线校正、填充零和包络线去除，对步骤五中获取的前寒武纪基底达马拉岩系主要地层岩石样品的光谱曲线进行预处理；采用连续统去除法光谱吸收特征分析方法对预处理后的光谱曲线进行诊断特征识别，获取前寒武纪基底达马拉岩系主要地层岩石的光谱诊断吸收峰对应的波长位置。

步骤七、构建可汗组地层遥感图谱识别标志。

若在步骤三中获取的ASTER假彩色三维影像上呈现浅蓝紫色夹浅黄绿色调，三维地貌为中低山，且地形切割强烈；在步骤四中获取的QuickBird真彩色融合影像上呈现灰黄色夹蓝紫色色调，爪状、弧状影纹；且地层岩石光谱诊断吸收峰波长为1098nm、2340nm，则该岩层为可汗组地层(图1)。

步骤八、构建罗辛组地层遥感图谱识别标志。

若在步骤三中获取的ASTER假彩色三维影像上呈现青蓝色深色调，三维地貌为高陡山地、耐风化，且地形切割强烈；在步骤四中获取的QuickBird真彩色融合影像上呈现浅灰色夹白色色调，辫状影纹；且地层岩石光谱诊断吸收峰波长为970nm、2342nm，则该岩层为罗辛组地层(图2)。

步骤九、构建楚斯组地层遥感图谱识别标志。

若在步骤三中获取的ASTER假彩色三维影像上呈现黄褐色夹深褐色调，三维地貌为中低山陡峭地貌、耐风化、地形切割强烈；在步骤四中获取的QuickBird真彩色融合影像上呈现浅灰色色调，局部呈黄绿色色调，梳状影纹；且地层岩石光谱无明显诊断吸收特征，则该岩层为楚斯组地层(图3)。

步骤十、构建卡里毕比组地层遥感图谱识别标志。

若在步骤三中获取的ASTER假彩色三维影像上呈现青蓝色深色调，三维地貌为低缓山坡、耐风化；在步骤四中获取的QuickBird真彩色融合影像上呈现较厚不连续浅白色色调，线状影纹；且地层岩石光谱诊断吸收峰波长为643nm、2150nm、2341nm，则该岩层为卡里毕比组地层(图4)。

步骤十一、识别前寒武纪基底达马拉岩系主要地层。

将步骤五获取的前寒武纪基底达马拉岩系主要地层岩石样品采集的GPS点位坐标，投影至步骤三获取的ASTER假彩色三维影像和步骤四获取的QuickBird真彩色融合影像上，结合步骤七、步骤八、步骤九、步骤十分别构建的主要地层遥感图谱识别标志(如下表1所示)，基于人工交互解译方法在ASTER假彩色三维影像和QuickBird真彩色融合影像上，识别并用红色线条勾勒出样品采集点位对应的主要地层分布范围，以该范围作为识别样本，采用监督分类法—最小距离法识别出全区前寒武纪基底达马拉岩系主要地层的分布范围。

表1前寒武纪基底达马拉岩系主要地层遥感图谱识别标志表

步骤十一中的人工解译指采用红色线条在ASTER假彩色三维影像和QuickBird真彩色融合影像上分别勾勒出前寒武纪基底达马拉岩系主要地层分布范围；训练样本指在ASTER假彩色三维影像和QuickBird真彩色融合影像上已确认的前寒武纪基底达马拉岩系主要地层的分布范围作为识别样本，基于该样本自动识别主要地层在工作区的全部分布范围；监督分类法指采用最小距离法进行分类。

上面对本发明的实施例作了详细说明，上述实施方式仅为本发明的最优实施例，但是本发明并不限于上述实施例，在本领域普通技术人员所具备的知识范围内，还可以在不脱离本发明宗旨的前提下做出各种变化。