[发明专利]一种钠交代岩型铀矿关键蚀变矿物组合获取方法

权利要求书

1.一种钠交代岩型铀矿关键蚀变矿物组合获取方法，其特征在于：该方法包括以下步骤：

步骤1.对高分遥感数据进行预处理；

所述的步骤1中对Worldview-3卫星高分遥感数据进行预处理；

步骤2.对步骤1中预处理后的高分遥感数据大气校正，获取地物反射率图像；所述的步骤2中对步骤1中预处理后的Worldview-3卫星高分遥感数据进行大气校正，获取地物反射率图像；

步骤3.建立赤铁矿信息、绿泥石信息和碳酸盐信息的提取模型；所述的步骤3具体包括以下步骤：步骤3.1，依据步骤2中大气校正后的Worldview-3卫星高分遥感数据对应的赤铁矿反射光谱特征波段，建立提取赤铁矿信息的模型；所述的步骤3.1中的Worldview-3卫星遥感高分数据对应的赤铁矿反射光谱特征波段为可见-近红外第3、7两个波段位置的吸收峰和第6波段的强反射峰，可见-近红外第3、6、7波段的反射率值分别记为B3、B6和B7；高分遥感数据提取赤铁矿信息的模型为(B6 GT B7)AND(B6 GT B3)AND((B6+2*B3)LT 3*B4)*(2*B6-B3-B7)，B4为可见-近红外第4波段的反射率值，GT表示大于，LT表示小于；

步骤3.2，依据步骤2中大气校正后的Worldview-3卫星高分遥感数据应的绿泥石反射光谱特征波段，建立提取绿泥石信息的模型；所述的步骤3.2中的高分遥感数据对应的绿泥石反射光谱特征波段为短波红外第6、7、8三个波段，6、7、8三个波段的反射率值分别记为Bˊ6、Bˊ7和Bˊ8；高分遥感数据提取绿泥石信息的模型为(Bˊ6GT Bˊ7)AND(Bˊ6GT Bˊ8)AND((9*Bˊ6+14*Bˊ8)GT 23*Bˊ7)；

步骤3.3，依据步骤2中大气校正后的Worldview-3卫星高分遥感数据应的碳酸盐反射光谱特征波段，建立提取碳酸盐信息的模型；所述的步骤3.3中的高分遥感数据对应的碳酸盐信息反射光谱特征波段为短波红外第6、7、8三个波段，其反射率值分别记为Bˊ6、Bˊ7和Bˊ8；高分遥感数据提取碳酸盐信息的模型为(Bˊ6GT Bˊ7)AND(Bˊ7GT Bˊ8)AND((( 9*Bˊ6+14*Bˊ8)/23*Bˊ7)GT 0.9)AND((( 9*Bˊ6+14*Bˊ8)/23*Bˊ7)LT 1.1)；

步骤4.使用步骤3得到的赤铁矿信息、绿泥石信息和碳酸盐信息提取模型，对步骤2中获取的地物反射率图像进行波段运算，分别获取赤铁矿信息、绿泥石信息和碳酸盐分布信息，从而完成钠交代岩型铀矿关键蚀变矿物组合的获取。

2.根据权利要求1所述的一种钠交代岩型铀矿关键蚀变矿物组合获取方法，其特征在于：所述的步骤1中的高分遥感数据的预处理包括辐射校正、几何校正、图像裁剪和镶嵌。

3.根据权利要求2所述的一种钠交代岩型铀矿关键蚀变矿物组合获取方法，其特征在于：所述的步骤2中的高分遥感数据的大气校正，采用快速大气校正方法或者择FLAASH大气校正方法获得。

4.根据权利要求3所述的一种钠交代岩型铀矿关键蚀变矿物组合获取方法，其特征在于：所述的步骤3.1中的赤铁矿光谱特征为B6＞B3且B6＞B7；赤铁矿的反射率光谱曲线公式为(B6+2B3)＜3B4。

5.根据权利要求4所述的一种钠交代岩型铀矿关键蚀变矿物组合获取方法，其特征在于：所述的步骤3.2中的绿泥石光谱特征为Bˊ6＞Bˊ7且Bˊ6＞Bˊ8；绿泥石反射率光谱曲线公式为

6.根据权利要求5所述的一种钠交代岩型铀矿关键蚀变矿物组合获取方法，其特征在于：所述的步骤3.3中的碳酸盐信息光谱特征为Bˊ6＞Bˊ7且Bˊ7＞Bˊ8；碳酸盐的反射率光谱曲线公式为

7.根据权利要求6所述的一种钠交代岩型铀矿关键蚀变矿物组合获取方法，其特征在于：所述的步骤4的具体步骤如下：

步骤4.1，使用高分遥感数据可见-近红外第3、4、6波段，通过步骤3.1中得到的赤铁矿提取模型(B6+2B3)＜3B4，对步骤2中获取的地物反射率图像进行波段运算，获取赤铁矿分布信息；

步骤4.2，使用高分遥感数据短波红外第6、7、8波段，通过步骤3.2得到的绿泥石提取模型(Bˊ6GT Bˊ7)AND(Bˊ6GT Bˊ8)AND((9*Bˊ6+14*Bˊ8)GT 23*Bˊ7)，对步骤2中获取的地物反射率图像进行波段运算，获取绿泥石分布信息；

步骤4.3，使用高分遥感数据短波红外第6、7、8波段，通过步骤3.3中得到的碳酸盐提取模型公式(Bˊ6GT Bˊ7)AND(Bˊ7GT Bˊ8)AND((( 9*Bˊ6+14*Bˊ8)/23*Bˊ7)GT 0.9)AND((( 9*Bˊ6+14*Bˊ8)/23*Bˊ7)LT1.1)，对步骤2中获取的地物反射率图像进行波段运算，获取碳酸盐分布信息。