CN114397251A - 一种植被覆盖区铀矿化蚀变信息遥感识别方法 - Google Patents
一种植被覆盖区铀矿化蚀变信息遥感识别方法 Download PDFInfo
- Publication number
- CN114397251A CN114397251A CN202111660772.4A CN202111660772A CN114397251A CN 114397251 A CN114397251 A CN 114397251A CN 202111660772 A CN202111660772 A CN 202111660772A CN 114397251 A CN114397251 A CN 114397251A
- Authority
- CN
- China
- Prior art keywords
- remote sensing
- vegetation
- vegetation coverage
- coverage area
- sensing image
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 230000004075 alteration Effects 0.000 title claims abstract description 56
- 229910052770 Uranium Inorganic materials 0.000 title claims abstract description 50
- JFALSRSLKYAFGM-UHFFFAOYSA-N uranium(0) Chemical compound [U] JFALSRSLKYAFGM-UHFFFAOYSA-N 0.000 title claims abstract description 50
- 230000033558 biomineral tissue development Effects 0.000 title claims abstract description 40
- 238000000034 method Methods 0.000 title claims abstract description 40
- 230000002159 abnormal effect Effects 0.000 claims abstract description 48
- 230000000873 masking effect Effects 0.000 claims abstract description 6
- 238000012847 principal component analysis method Methods 0.000 claims abstract description 5
- 238000007781 pre-processing Methods 0.000 claims abstract description 4
- 229910052500 inorganic mineral Inorganic materials 0.000 claims description 24
- 239000011707 mineral Substances 0.000 claims description 24
- 238000011160 research Methods 0.000 claims description 17
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 16
- 238000004364 calculation method Methods 0.000 claims description 6
- 238000012952 Resampling Methods 0.000 claims description 3
- 238000000605 extraction Methods 0.000 abstract description 9
- 230000000452 restraining effect Effects 0.000 abstract 1
- 241000132092 Aster Species 0.000 description 36
- 230000003595 spectral effect Effects 0.000 description 13
- 125000002887 hydroxy group Chemical group [H]O* 0.000 description 11
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 description 10
- 238000010521 absorption reaction Methods 0.000 description 6
- 230000014509 gene expression Effects 0.000 description 6
- 239000011435 rock Substances 0.000 description 6
- 238000001228 spectrum Methods 0.000 description 5
- 230000001629 suppression Effects 0.000 description 5
- 230000005856 abnormality Effects 0.000 description 4
- 239000000443 aerosol Substances 0.000 description 4
- 230000001186 cumulative effect Effects 0.000 description 4
- 230000005764 inhibitory process Effects 0.000 description 4
- 229910052742 iron Inorganic materials 0.000 description 4
- 241000195493 Cryptophyta Species 0.000 description 3
- 238000005065 mining Methods 0.000 description 3
- 238000002310 reflectometry Methods 0.000 description 3
- 241000196324 Embryophyta Species 0.000 description 2
- 229910052782 aluminium Inorganic materials 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 239000002734 clay mineral Substances 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 230000002401 inhibitory effect Effects 0.000 description 2
- 229910052749 magnesium Inorganic materials 0.000 description 2
- 239000011777 magnesium Substances 0.000 description 2
- 230000001089 mineralizing effect Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000000985 reflectance spectrum Methods 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 239000002689 soil Substances 0.000 description 2
- 238000010186 staining Methods 0.000 description 2
- 238000009966 trimming Methods 0.000 description 2
- OYPRJOBELJOOCE-UHFFFAOYSA-N Calcium Chemical compound [Ca] OYPRJOBELJOOCE-UHFFFAOYSA-N 0.000 description 1
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 1
- 206010013883 Dwarfism Diseases 0.000 description 1
- 102000006463 Talin Human genes 0.000 description 1
- 108010083809 Talin Proteins 0.000 description 1
- 238000000862 absorption spectrum Methods 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000005452 bending Methods 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 239000011575 calcium Substances 0.000 description 1
- 229910052791 calcium Inorganic materials 0.000 description 1
- 150000001768 cations Chemical class 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 239000013078 crystal Substances 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- GUJOJGAPFQRJSV-UHFFFAOYSA-N dialuminum;dioxosilane;oxygen(2-);hydrate Chemical compound O.[O-2].[O-2].[O-2].[Al+3].[Al+3].O=[Si]=O.O=[Si]=O.O=[Si]=O.O=[Si]=O GUJOJGAPFQRJSV-UHFFFAOYSA-N 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 238000004043 dyeing Methods 0.000 description 1
- 230000005670 electromagnetic radiation Effects 0.000 description 1
- 230000005274 electronic transitions Effects 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 229910052751 metal Inorganic materials 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 229910052901 montmorillonite Inorganic materials 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 238000000513 principal component analysis Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
- 229910001727 uranium mineral Inorganic materials 0.000 description 1
- 239000002023 wood Substances 0.000 description 1
- 239000010151 yanghe Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N2021/1793—Remote sensing
- G01N2021/1797—Remote sensing in landscape, e.g. crops
Landscapes
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Image Processing (AREA)
Abstract
本发明属于遥感地质技术领域,具体公开一种植被覆盖区火山岩型铀矿化信息遥感识别方法,如下:步骤1、在植被覆盖区进行遥感数据源选择;步骤2、对所选遥感影像预处理,得到植被覆盖区遥感影像图;步骤3、基于像元二分模型,对遥感影像进行植被覆盖度估算,确定高植被覆盖区和中低植被覆盖区;步骤4、对中低植被和高植被覆盖区,去除植被以外干扰之后,抑制中低植被覆盖区植被,掩膜高植被覆盖区植被,得到去除干扰因素的植被覆盖区影像图;步骤5、针对去除干扰因素的植被覆盖区影像图,采用主成分分析法提取铀矿化蚀变异常信息,完成植被覆盖区铀矿化蚀变信息遥感识别。本发明的方法解决了减轻植被干扰对遥感蚀变信息提取影响的技术问题。
Description
技术领域
本发明属于遥感地质技术领域,具体涉及一种植被覆盖区火山岩型铀矿化信息遥感识别方法。
背景技术
铀矿是具有放射性的国家重要的战略性矿产。火山岩型铀矿是我国乃至全球铀矿资源的主要类型,按照“多探少采”的总体要求,硬岩型铀矿勘查立足“增加储备”,重点围绕已有老铀矿田深部及外围探索。目前,遥感矿化蚀变信息提取的理论方法众多,在裸露地区的矿产资源勘查中遥感技术具有独特的优势。然而,大量植被覆盖区铀矿地质工作程度极低-空白区,其成矿地质构造复杂,成矿条件优越,矿产资源丰富,运用传统遥感蚀变信息提取难度较大。
植物对电磁波辐射的吸收、反射和散射作用构成了植物的特征光谱。研究发现,可见光波段的反射低谷分别位于0.45μm、0.67μm波长处的蓝光、红光;近红外波段形成的“反射高原”位于0.7~1.3μm波长处,波长1.3~2.4μm之间存在2个明显的吸收谷,分别位于1.45μm和1.94μm处。岩石矿物的吸收光谱特征主要由于电子跃迁或振动过程所产生的。矿物在可见光-近红外光谱带,电子跃迁形成的特征谱带主要是受到铁元素的影响,其中占主导地位的是Fe2+、Fe3+,在振动过程产生的吸收谱带特征主要受水、羟基基团和碳酸根离子的影响。羟基类是对以粘土矿物为特征矿物的蚀变类型的统称,粘土矿物是热液矿床形成过程中形成的蚀变矿物,主要由羟基、水和金属阳离子组成,光谱特征与羟基和水基团的晶体结构有关。该类蚀变的特征矿物较多,与铀矿化相关的有钙蒙脱石和碎片状水云母。该类蚀变矿物的光谱曲线中,强吸收谱带在1.4μm、1.9μm和2.2~2.4μm三个波长区,其中1.4μm的谱带是矿物组成中羟基团引起的,1.9μm的谱带则是矿物水分子的弯曲振动贡献,而2.2~2.4μm之间的谱带也是羟基团产生的。因此,植被光谱曲线与含有羟基基团蚀变矿物光谱曲线在1.4μm、1.9μm附近都有吸收谱带,在2.0~2.4μm波段处植被反射率较低,在1.6~2.4μm之间植被光谱曲线呈现出逐渐下降的特征,而在该区间上粘土矿物矿物的光谱曲线也呈现出了下降的特征,两者在总体特征上具有较大的相似性,所以传感器接收到的蚀变矿物反射光谱受到了植被反射光谱的很大影响。因此,在植被覆盖区提取铀矿化蚀变信息前对该地区的植被干扰进行抑制和去除是很有必要的。
国内外研究者利用遥感数据在低植被覆盖区或基岩裸露区蚀变信息提取已经做了大量的研究,并且取得了良好的效果。然而在植被覆盖区进行蚀变信息提取的研究和应用当前仍是一个难点。
发明内容
本发明的目的在于提供一种植被覆盖区火山岩型铀矿化信息遥感识别方法,该方法基于像元二分模型进行植被覆盖度估算,确定高植被覆盖区和中低植被覆盖区,釆用植被抑制算法抑制中低植被覆盖区,掩膜高植被覆盖区及其他干扰因素,在去除干扰因素的基础上,采用主成分分析方法进行植被覆盖区遥感蚀变信息提取,解决了减轻植被干扰对遥感蚀变信息提取影响的技术问题。
实现本发明目的的技术方案如下:一种植被覆盖区铀矿化蚀变信息遥感识别方法,该方法的具体步骤如下:
步骤1、在植被覆盖区进行遥感数据源选择;
步骤2、针对上述步骤1中得到的遥感影像进行预处理,得到遥感影像图;
步骤3、基于像元二分模型,对上述步骤2中遥感影像进行植被覆盖度估算,确定高植被覆盖区和中低植被覆盖区;
步骤4、针对上述步骤3中确定的中低植被覆盖区和高植被覆盖区,去除植被以外的干扰之后,抑制中低植被覆盖区植被,掩摸高植被覆盖区植被,得到去除干扰因素的植被覆盖区影像图;
步骤5、针对步骤4中得到的去除干扰因素的植被覆盖区影像图,采用主成分分析法提取铀矿化蚀变异常信息,从而完成植被覆盖区铀矿化蚀变信息的遥感识别。
所述步骤2的具体步骤如下:
步骤2.1、针对上述步骤1中的遥感影像进行重采样;
步骤2.2、针对上述步骤2.1中重采样后的遥感影像进行辐射定标,得到辐射定标后的光谱遥感影像;
步骤2.3、针对上述步骤2.2中辐射定标后的遥感影像去边框;
步骤2.4、针对上述步骤2.3中去边框后的遥感影像进行大气校正。
所述步骤3的具体步骤如下:
步骤3.1、归一化上述步骤2.4中大气校正后遥感影像的植被指数;
步骤3.2、去除上述步骤3.1中遥感影像归一化植被指数NDVI的异常值;
步骤3.3.1、对上述步骤3.2中去除NDVI异常值后的遥感影像归一化植被指数NDVI进行最值统计;
步骤3.3.2、将上述步骤3.3.1中最值统计后的遥感影像归一化植被指数NDVI生成NDVIsoil和NDVIveg参数文件;
步骤3.3.3、将上述步骤3.3.2中得到遥感影像归一化植被指数NDVI的参数文件,进行植被覆盖度计算;
步骤3.3.4、去除上述步骤3.3.3中得到植被覆盖度的异常值;
步骤3.3.5、针对上述步骤3.3.4中得到去除异常值的植被覆盖度进行等级划分,确定中低植被覆盖区和高植被覆盖区。
所述步骤4的具体步骤如下:
步骤4.1、去除水系、阴影、云雪和第四系覆盖区干扰因素;
步骤4.2、抑制中低植被覆盖区植被;
步骤4.3、去除高植被覆盖区植被;
步骤4.4、在上述步骤4.1、步骤4.2基础上,得到中低植被覆盖区去除干扰因素后的影像图;
步骤4.5、建立上述步骤4.1、步骤4.3高植被覆盖区的合膜,得到高植被覆盖区掩膜后的影像图。
所述的步骤5的具体步骤如下:
步骤5.1、针对步骤4.4中得到的去除干扰因素的植被覆盖区影像图,进行蚀变矿物异常主分量提取,得到植被覆盖区的异常主分量图像灰度值;
步骤5.2、针对步骤5.1中得到的植被覆盖区的异常主分量图像灰度值,得到研究区的蚀变信息,从而完成植被覆盖区铀矿化蚀变信息的遥感识别。
所述的步骤2.2中对遥感影像数据进行辐射定标采用ENVI5.3软件选择符合FLAASH大气校正要求的遥感影像数据参数,参数包括定标类型、存储顺序和缩放系数,输出得到辐射定标后的遥感影像。
所述的步骤2.3中采用掩摸方法去除遥感影像的边框。
所述的步骤2.3中的掩摸方法去除遥感影像的边框的全边框掩摸计算公式如下:b3>0=b3new,b6>0=b6new,b3new+b6new=b1,b5>0=b5new,b7>0=b7new,b5new+b7new=b2,b1+b2=bmask。
所述的步骤3.1中的遥感影像的归一化植被指数NDVI=(NIR-R)/(NIR+R)。
所述的步骤3.2中的去除异常值的公式为:-1>NDVI<1。
所述的步骤3.3.3中的植被覆盖度
VFC=(NDVI-NDVIsoil)/(NDVIveg-NDVIsoli)。
所述的步骤3.3.4中植被覆盖度去除异常值d1的公式为:0.0>d1<1.0。
本发明的有益技术效果在于:(1)本发明的方法引入像元二分模型进行植被覆盖度遥感估算,划分了植被覆盖度等级,准确定位于遥感影像图上,初步确定了研究区内植被覆盖情况。
(2)研究区内干扰蚀变信息提取的主要干扰因素包括植被、水系、阴影、云雪和第四系等。根据本发明的方法中像元二分模型设计的植被覆盖度估算,确定了研究区的中低植被覆盖区和高植被覆盖区。本发明在中低植被覆盖区釆用植被抑制算法最大程度抑制植被光谱信息,在高植被覆盖区和水系、阴影、云雪与第四系采用分别建立掩膜、相乘取得合膜的方法进行了去干扰操作。
(3)本发明的方法在去除干扰因素的基础上,采用“植被抑制算法抑制中低植被覆盖区+掩膜高植被覆盖区及其他干扰因素+主成分分析+异常主分量+无损线性拉伸+异常切割”的方法,提取植被覆盖区铀矿化蚀变异常信息,能够进行铀矿产资源寻找和预测。
(4)采用本发明的方法青龙铀矿田地区蚀变信息提取过程中植被是做为重要干扰因素被去除。
附图说明
图1为本发明所提供的青龙铀矿田及外围地区遥感异常分布图。
具体实施方式
下面结合附图和实施例对本发明作进一步详细说明。
下面以青龙铀矿田地区ASTER数据为例,描述对本发明所提的一种植被覆盖区铀矿化蚀变信息遥感识别方法,该方法的具体步骤如下:
步骤1、在植被覆盖区进行遥感数据源选择
在中等空间分辨率的多光谱遥感影像中Aster数据相对于ETM+数据增加了7个波段,ASTER数据在短波红外波段有6个波段,而ETM+数据在短波红外波段只有1个波段(如下表1所示),具有地面分辨率高、通道范围窄及波段多等优点,更利于矿物蚀变信息的提取。本实施例中选择覆盖研究区L1A级的ASTER数据,共1景,使用波段为可见光、近红外和短波红外波段。该遥感影像数据质量较好,云覆盖较少,图像清晰。
表1 Aster与ETM+遥感数据各波段对照
步骤2、针对上述步骤1中得到的遥感影像数据进行预处理,得到Aster遥感影像;所述的步骤2的具体步骤如下:
步骤2.1、针对上述步骤1中得到的Aster遥感影像数据进行重采样
由于Aster遥感影像数据中的可见光、近红外和短波红外分辨率不一致,将30米空间分辨率的Aster遥感影像数据中的SWIIR波段重采样至与VNIR波段相同的15m空间分辨率。然后与VNIR波段按波段顺序合并成新的文件,形成9个波段的遥感影像数据。
步骤2.2、针对上述步骤2.1中重采样后的的Aster遥感影像数据进行辐射定标,得到辐射定标后的Aster遥感影像
利用ENVI5.3软件中Radiometric Calibration工具单击Apply FLAASH Settings按钮,自动选择符合FLAASH大气校正要求的Aster遥感影像数据参数,包括定标类型(Calibration Type)、存储顺序(Type)和缩放系数(Scale Factor),输出得到辐射定标后的Aster遥感影像。
步骤2.3、针对上述步骤2.2中得到的辐射定标后的Aster遥感影像去边框
由于Aster传感器在可见光-近红外波长区间与短波红外区间成像时间略有差异,导致每个波段位置有所不同,主要体现在边界的一些信息不重合,因此,需要去Aster遥感影像的边框,剔除不重合部分,保留重合区域。Aster遥感影像去边框的具体步骤如下:
分别对ASTER数据3波段、6波段和5波段、第7波段分别进行波段“与”运算,构建全边框掩模,应用掩摸方法,使数据覆盖范围完全一致。其全边框掩摸计算公式如下:
b3>0=b3new (1)
b6>0=b6new (2)
b3new+b6new=b1 (3)
b5>0=b5new (4)
b7>0=b7new (5)
b5new+b7new=b2 (6)
b1+b2=bmask (7)
式中,b3、b5、b6、b7为ASTER数据的3波段、5波段、6波段、7波段;b3 new、b5 new、b6new、b7 new为ASTER数据新的3波段、5波段、6波段、7波段;b1、b2为东西向切边框数据、南北向切边框数据;bmask为全边框掩膜数据。
步骤2.4、针对上述步骤2.3中去边框后的Aster遥感影像进行大气校正
利用ENVI5.3软件中Flaash AtMospheric correction工具进行参数设置,参数包括:传感器基本信息、大气模型(Atmospheric Model)、气溶胶模型(Aerosol Model)、气溶胶反演方法(Aerosol Retrieval)、初始能见度(Initial Visibility),其余参数默认即可,得到大气校正后的Aster遥感影像。
步骤3、基于像元二分模型,对上述步骤2.4中大气校正后的Aster遥感影像进行植被覆盖度估算,确定高植被覆盖区和中低植被覆盖区;所述的步骤3的具体步骤如下:
步骤3.1、归一化上述步骤2.4中大气校正后的Aster遥感影像的植被指数
采用归一化植被指数(NDVI)做为计算研究区遥感影像的植被覆盖度指示因子,归一化植被指数NDVI值等于步骤2得到的Aster数据近红外波段与可见光红波段数值之差与波段数值之和之比。Aster遥感影像的归一化植被指数NDVI的公式为:
NDVI=(NIR-R)/(NIR+R) (8)
式中,NDVI为归一化植被指数;NIR为Aster遥感影像的近红外波段反射率;R为Aster遥感影像的红光波段反射率。
步骤3.2、去除上述步骤3.1中得到的Aster遥感影像归一化植被指数NDVI的异常值
由于大气校正后的Aster遥感影像结果有部分像元为负值,主要集中在阴影地区,这部分区域计算的NDVI在[-1,1]之外,对这部分像元进行处理,即NDVI值大于1的变为1,小于-1的变成-1,将步骤3.1得到的Aster遥感影像的归一化植被指数NDVI结果在[-1,1],去除异常值Band math表达式为:
-1>NDVI<1 (9)
式中,NDVI为归一化植被指数。
步骤3.3、基于像元二分模型,针对上述步骤3.2中得到的去除NDVI异常值后的Aster遥感影像归一化植被指数进行植被覆盖度估算
步骤3.3.1、对上述步骤3.2中去除NDVI异常值后的Aster遥感影像归一化植被指数NDVI进行最值统计
由于遥感成像受多种因素(如气象条件、观测角度、观测时间、太阳高度等)影响,因此在遥感图像中,难免会存在异常值。本发明中研究区地物类型分为林地、农业用地、城市用地、水体、其他。为了获取这5种地物类型对应的NDVI最大值和最小值,基于统计学原理,通过取一定置信度的值来确定步骤3.2得到的NDVI最大值的NDVImax和最小值NDVImin,这里取2%的置信度,即累积概率在2%和98%处的值。
取林地像元累积概率(ACC PCT)2%处的DN值作为NDVI最小值统计,累积概率98%处的DN值作为NDVI最大值。若累积概率无2%或98%,可取累积概率2%处前后DN值均值作为NDVI最小值。林地NDVImin=(0.365+0.380)/2;林地NDVImax同理。
同样方法得到其他地物覆盖类型的NDVI阈值,如下表2所示。
表2研究区地物类型NDVI阈值
水体没有植被(水藻不属于植被),认为这部分区域的植被覆盖度为0。
步骤3.3.2、将上述步骤3.3.1中最值统计后的Aster遥感影像归一化植被指数NDVI生成NDVIsoil和NDVIveg参数文件
根据步骤3.3.1得到的研究区土地覆盖类型的NDVI阈值,使用bandmath工具分别生成NDVIsoil和NDVIveg参数文件,其表达式如下:
NDVIsoil=c1*NDVImin1+c2*NDVImin2+c3*NDVImin3+b4*NDVImin4+b5*NDVImin5(10)
NDVIveg=c1*NDVImax1+c2*NDVImax2+c3*NDVImax3+c4*NDVImin4+c5*NDVImin5(11)
式中,c1为林地文件,c2为农业用地文件,c3为城市用地文件,c4为水体文件,c5为其他文件;NDVIsoil为完全是裸土或无植被覆盖区域像元的NDVI值,、NDVImin1为林地的无植被覆盖区域像元的NDVI最小值,、NDVImin2为农业用地的无植被覆盖区域像元的NDVI最小值,、NDVImin3为城市用地的无植被覆盖区域像元的NDVI最小值,通常认为水体没有植被(水藻不属于植被),认为这部分区域的植被覆盖度为0,即NDVImin4为0,NDVImin5为其他的无植被覆盖区域像元的NDVI最小值;NDVIveg为完全被植被所覆盖的像元的NDVI值,NDVImax1为林地的完全被植被所覆盖的像元的NDVI最大值,NDVImax2为农业用地的完全被植被所覆盖的像元的NDVI最大值,NDVImax3为城市用地的完全被植被所覆盖的像元的NDVI最大值,通常认为水体没有植被(水藻不属于植被),认为这部分区域的植被覆盖度为0,即NDVImax4为0,NDVImax5为其他的完全被植被所覆盖的像元的NDVI最大值。
步骤3.3.3、将上述步骤3.3.2中得到Aster遥感影像归一化植被指数NDVI的参数文件,进行植被覆盖度计算
NDVImin为影像土地覆盖类型对应的NDVI最小值,NDVImax为影像土地覆盖类型对应的NDVI最大值。利用步骤3.3.2得到的NDVIsoil和NDVIveg参数文件带入计算公式(12),使用Band Math表达式(13)进行植被覆盖度VFC计算。其表达式为
VFC=(NDVI-NDVIsoil)/(NDVIveg-NDVIsoli) (12)
式中:NDVI为“去除异常值的NDVI;NDVIsoil为完全是裸土或无植被覆盖区域像元的NDVI值;NDVIveg为完全被植被所覆盖的像元的NDVI值
步骤3.3.4、去除上述步骤3.3.3中得到植被覆盖度的异常值
分析步骤3.3.3得到的植被覆盖度估算结果,可能会存在一些异常值,即值在[0,1]之外。这些异常值是由于在NDVI置信度(2%—98%)之外的那部分像元产生的(也包括NDVI异常像元)。这些像元数量不多,大约占2%左右。还有背景及水体区域的植被覆盖度的值为-NaN,即无效值。因此利用Band Math表达式(13)去除这部分异常值,可以将小于0的值变成0,大于1的值变成1。
其表达式为:
0.0>d1<1.0 (13)
式中,d1为植被覆盖度。
步骤3.3.5、针对上述步骤3.3.4中得到去除异常值的植被覆盖度进行等级划分,确定中低植被覆盖区和高植被覆盖区
将步骤3.3.4得到的去除异常值的植被覆盖度(VFC)计算,利用ENVI5.3软件中的New Raster Color Slices工具进行分割,根据植被指数分级标准,将植被覆盖度划分为低植被覆盖、较低植被覆盖、较高植被覆盖、高植被覆盖共计4个等级(中低植被覆盖区和高植被覆盖区)(如下表3所示)。
表3植被指数分级标准
等级 | 等级名称 | 范围 |
Ⅰ | 低植被覆盖 | 0~0.10 |
Ⅱ | 较低植被覆盖 | 0.10~0.25 |
Ⅲ | 较高植被覆盖 | 0.25~0.40 |
Ⅳ | 高植被覆盖 | 0.40~0.55 |
步骤4、针对上述步骤3.3.5中划分研究区为中低植被覆盖区和高植被覆盖区,在去除植被以外的干扰之后,中低植被覆盖区采用植被抑制算法,高植被覆盖区进行掩摸去除植被,得到去除干扰因素的植被覆盖区的影像图;所述的步骤4的具体步骤如下:
步骤4.1、去除植被以外的干扰因素
针对水系、阴影、云雪和第四系覆盖物等干扰因素采用掩膜技术去除。
步骤4.2、去除中低植被覆盖区的干扰地物因素
根据步骤3.3.5得到的中低植被覆盖区,使用ENVI 5.3光谱分析软件中的植被抑制算法(Vegetation Suppression)去除中低植被覆盖区的干扰地物因素,突出地面光谱信息和增强裸露地面岩石光谱信息,最大程度抑制了植被光谱信息,减小对蚀变信息提取结果的影响。
表4中特征值NDVIQ表示中低植被覆盖区植被抑制前的NDVI特征,NDVIH表示植被抑制后的特征,其中NDVI的最大值从0.302326降到了0.30074,平均值也,0.106381降到了0.097708,说明了被抑制方法在本研宄区对抑制植被光谱信息,减小植被因素对烛变信息提取结果的影响是适用的和有效的。
表4中低植被覆盖区植被抑制前后统计特征
特征值 | 最小值 | 最大值 | 平均值 | 标准差 |
NDVI<sub>Q</sub> | -0.647059 | 0.302326 | 0.106381 | 0.109513 |
NDVI<sub>H</sub> | -0.519627 | 0.300474 | 0.097708 | 0.106364 |
步骤4.3、高植被覆盖区干扰地物因素去除
对上述步骤3.3.5中确定的高植被覆盖区进行掩膜处理,去除高植被覆盖区的干扰因素。
步骤4.4、中低植被覆盖区去干扰后影像图
在上述步骤4.1、步骤4.2基础上,得到中低植被覆盖区去除干扰因素后的影像图。
步骤4.5、建立高植被覆盖区的合膜,得到高植被覆盖区掩膜后的影像图
将步骤4.1、步骤4.3得到的高植被覆盖区掩膜、水系、阴影、云雪和第四系覆盖物掩膜一一相乘,得到高植被覆盖区的掩膜后的影像图。
步骤5、针对步骤4.4、步骤4.5中得到的去除干扰因素的植被覆盖区影像图,采用主成分分析法提取铀矿化蚀变异常信息,从而完成植被覆盖区铀矿化蚀变信息的遥感识别;所述的步骤5的具体步骤如下:
步骤5.1、针对步骤4.4、步骤4.5中得到的去除干扰因素的植被覆盖区影像图,进行蚀变矿物异常主分量提取,得到植被覆盖区的异常主分量图像灰度值
针对步骤4.4、步骤4.5中得到的去除干扰因素的植被覆盖区影像图,基于蚀变矿物的波谱特征选择主成分变换波段,分析变换后的特征向量载荷因子的大小和符号,确定每个波段对矿化蚀变矿物的光谱响应贡献,判别适合提取蚀变矿物异常的主分量图像。
通过对研究区蚀变分带模式中各蚀变带主要矿物波谱特征的分析,选取了蚀变矿物在ASTER近红外至短波红外具有光谱特征的波段,采用主成分分析法,分析变换后的特征向量载荷因子的大小和符号,确定每个波段对矿化蚀变矿物的光谱响应贡献,判别适合提取蚀变矿物异常的主分量图像。其中铝羟基蚀变提取选择波段1、3、4、6(表5),镁羟基蚀变提取选择波段1、3、4、8(表6);铁染提取选择波段1、2、3、4(表7)。
表5铝羟基蚀变PC1-PC4的特征向量表
特征向量 | Band1 | Band3 | Band4 | Band6 |
PC1 | 0.221349 | 0.472197 | 0.664537 | 0.535187 |
PC2 | 0.655237 | 0.569512 | -0.443477 | -0.222820 |
PC3 | 0.618439 | -0.630142 | -0.123135 | 0.453091 |
PC4 | -0.373099 | 0.235810 | -0.588690 | 0.677226 |
表6镁羟基蚀变PC1-PC4的特征向量表
特征向量 | Band1 | Band3 | Band4 | Band8 |
PC1 | 0.227621 | 0.484267 | 0.681135 | 0.499728 |
PC2 | 0.666167 | 0.512507 | -0.537614 | -0.067309 |
PC3 | 0.317209 | -0.564722 | -0.235985 | 0.724416 |
PC4 | -0.635444 | 0.428848 | -0.437421 | 0.470067 |
表7铁染蚀变PC1-PC4的特征向量表
特征向量 | Band1 | Band2 | Band3 | Band4 |
PC1 | 0.246988 | 0.349838 | 0.524998 | 0.735518 |
PC2 | 0.522297 | 0.568846 | 0.213519 | -0.598356 |
PC3 | 0.441137 | 0.231754 | -0.806775 | 0.317495 |
PC4 | 0.686732 | -0.707331 | 0.167037 | -0.013401 |
步骤5.2、针对步骤5.1中得到的植被覆盖区的异常主分量图像灰度值,得到研究区的蚀变信息;从而完成植被覆盖区铀矿化蚀变信息的遥感识别。
先将步骤5.1中得到的异常主分量图像灰度值范围无损线性拉伸至0到255之间,基于分量灰度图像元值的似正态分布的直方图,用σ表征正态分布曲线尺度进行数据切割,异常的分级利用标准离差的倍数Kσ。把主分量分析结果均值(X)理解为代表区域背景,利用计算切割时使用的公式(15)划分异常强度等级,然后进行中值滤波对噪声和孤立点驱除,得到研究区蚀变信息结果。研究区蚀变信息的计算切割公式为:
L=X+kσ (15)
式中,L代表密度分割时使用的低值;X为所选择分量的灰度图的均值;k值一般取在1.5—3左右;σ所选择分量的灰度图的标准差。
步骤5.3、根据步骤5.2中得到的研究区蚀变信息结果图中遥感异常分布的面积及强度,获取铀矿化蚀变信息分布特征
根据步骤5.2中研究区蚀变信息结果图中遥感异常分布的面积及强度,发现这些异常呈NE向展布,遥感蚀变异常部位往往受线性断裂、环状构造控制,发育于线性断裂、环状构造内。铀矿化(异常)点与遥感蚀变异常部位相互重叠(如图1所示)。铀矿化与蚀变关系密切。
根据遥感异常分布图与铀成矿部位进行分析,研究区羟基和铁染两类异常分布范围与目前已发现的铀矿床、铀矿点有对应关系。
羟基异常主要分布在木头凳、燕河营等区这些异常呈北东向展布,部分受火山岩系、构造控制,青龙铀矿化聚集区的干沟一带的羟基异常分布,受火山机构控制,与该区铀矿化可能有一定的关系。
铁染异常的分布与羟基异常分布区基本一致,分布范围稍大,主要分布在中下侏罗统中基性、中酸性火山岩系及基底老地层中。
上面结合附图和实施例对本发明作了详细说明,上述实施方式仅为本发明的最优实施例,但是本发明并不限于上述实施例,在本领域普通技术人员所具备的知识范围内,还可以在不脱离本发明宗旨的前提下做出各种变化。
Claims (12)
1.一种植被覆盖区铀矿化蚀变信息遥感识别方法,其特征在于,该方法的具体步骤如下:
步骤1、在植被覆盖区进行遥感数据源选择;
步骤2、针对上述步骤1中得到的遥感影像进行预处理,得到遥感影像图;
步骤3、基于像元二分模型,对上述步骤2中遥感影像进行植被覆盖度估算,确定高植被覆盖区和中低植被覆盖区;
步骤4、针对上述步骤3中确定的中低植被覆盖区和高植被覆盖区,在去除植被以外的干扰之后,抑制中低植被覆盖区植被,掩摸高植被区植被,得到去除干扰因素的植被覆盖区的影像图;
步骤5、针对步骤4中得到的去除干扰因素的植被覆盖区影像图,采用主成分分析法提取铀矿化蚀变异常信息,从而完成植被覆盖区铀矿化蚀变信息的遥感识别。
2.根据权利要求1所述的一种植被覆盖区铀矿化蚀变信息遥感识别方法,其特征在于:所述的步骤2的具体步骤如下:
步骤2.1、针对上述步骤1中的遥感影像进行重采样;
步骤2.2、针对上述步骤2.1中重采样后的遥感影像进行辐射定标,得到辐射定标后的遥感影像;
步骤2.3、针对上述步骤2.2中辐射定标后的遥感影像去边框;
步骤2.4、针对上述步骤2.3中去边框后的遥感影像进行大气校正。
3.根据权利要求2所述的一种植被覆盖区铀矿化蚀变信息遥感识别方法,其特征在于:所述的步骤3的具体步骤如下:
步骤3.1、归一化上述步骤2.4中大气校正后遥感影像的植被指数;
步骤3.2、去除上述步骤3.1中遥感影像归一化植被指数NDVI的异常值;
步骤3.3.1、对上述步骤3.2中去除NDVI异常值后的遥感影像归一化植被指数NDVI进行最值统计;
步骤3.3.2、将上述步骤3.3.1中最值统计后的遥感影像归一化植被指数NDVI生成NDVIsoil和NDVIveg参数文件;
步骤3.3.3、将上述步骤3.3.2中得到遥感影像归一化植被指数NDVI的参数文件,进行植被覆盖度计算;
步骤3.3.4、去除上述步骤3.3.3中得到植被覆盖度的异常值;
步骤3.3.5、针对上述步骤3.3.4中得到去除异常值的植被覆盖度进行等级划分,确定中低植被覆盖区和高植被覆盖区。
4.根据权利要求3所述的一种植被覆盖区铀矿化蚀变信息遥感识别方法,其特征在于:所述的步骤4的具体步骤如下:
步骤4.1、去除水系、阴影、云雪和第四系覆盖物干扰因素;
步骤4.2、抑制中低植被覆盖区植被;
步骤4.3、去除高植被覆盖区植被;
步骤4.4、在上述步骤4.1、步骤4.2基础上,得到中低植被覆盖区去除干扰因素后的影像图;
步骤4.5、建立上述步骤4.1、步骤4.3高植被覆盖区的的合膜,得到高植被覆盖区掩膜后的影像图。
5.根据权利要求4所述的一种植被覆盖区铀矿化蚀变信息遥感识别方法,其特征在于:所述的步骤5的具体步骤如下:
步骤5.1、针对步骤4.4、步骤4.5中得到的去除干扰因素的植被覆盖区影像图,进行蚀变矿物异常主分量提取,得到植被覆盖区的异常主分量图像灰度值;
步骤5.2、针对步骤5.1中得到的植被覆盖区的异常主分量图像灰度值,得到研究区的蚀变信息,从而完成植被覆盖区铀矿化蚀变信息的遥感识别。
6.根据权利要求5所述的一种植被覆盖区铀矿化蚀变信息遥感识别方法,其特征在于:所述的步骤2.2中对遥感影像数据进行辐射定标采用ENVI5.3软件选择符合FLAASH大气校正要求的遥感影像数据参数,参数包括定标类型、存储顺序和缩放系数,输出得到辐射定标后的遥感影像。
7.根据权利要求6所述的一种植被覆盖区铀矿化蚀变信息遥感识别方法,其特征在于:所述的步骤2.3中采用掩摸方法去除遥感影像的边框。
8.根据权利要求7所述的一种植被覆盖区铀矿化蚀变信息遥感识别方法,其特征在于:所述的步骤2.3中的掩摸方法去除遥感影像的边框的全边框掩摸计算公式如下:b3>0=b3new,b6>0=b6new,b3new+b6new=b1,b5>0=b5new,b7>0=b7new,b5new+b7new=b2,b1+b2=bmask。
9.根据权利要求8所述的一种植被覆盖区铀矿化蚀变信息遥感识别方法,其特征在于:所述的步骤3.1中的遥感影像的归一化植被指数NDVI=(NIR-R)/(NIR+R)。
10.根据权利要求9所述的一种植被覆盖区铀矿化蚀变信息遥感识别方法,其特征在于:所述的步骤3.2中的去除异常值的公式为:-1>NDVI<1。
11.根据权利要求10所述的一种植被覆盖区铀矿化蚀变信息遥感识别方法,其特征在于:所述的步骤3.3.3中的植被覆盖度VFC=(NDVI-NDVIsoil)/(NDVIveg-NDVIsoli)。
12.根据权利要求11所述的一种植被覆盖区铀矿化蚀变信息遥感识别方法,其特征在于:所述的步骤3.3.4中植被覆盖度去除异常值d1的公式为:0.0>d1<1.0。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111660772.4A CN114397251A (zh) | 2021-12-31 | 2021-12-31 | 一种植被覆盖区铀矿化蚀变信息遥感识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111660772.4A CN114397251A (zh) | 2021-12-31 | 2021-12-31 | 一种植被覆盖区铀矿化蚀变信息遥感识别方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN114397251A true CN114397251A (zh) | 2022-04-26 |
Family
ID=81229010
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111660772.4A Pending CN114397251A (zh) | 2021-12-31 | 2021-12-31 | 一种植被覆盖区铀矿化蚀变信息遥感识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114397251A (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103383348A (zh) * | 2013-05-28 | 2013-11-06 | 吉林大学 | 植被覆盖区高光谱遥感蚀变矿物提取方法 |
CN108717044A (zh) * | 2018-05-24 | 2018-10-30 | 青海师范大学 | 一种去除植被覆盖影响的表层土壤含水量卫星遥感估算方法 |
CN108985130A (zh) * | 2017-05-31 | 2018-12-11 | 核工业北京地质研究院 | 一种基于缓冲区干扰信息去除的蚀变信息提取方法 |
CN109165580A (zh) * | 2018-08-09 | 2019-01-08 | 东华理工大学 | 高密度植被覆盖区条件下的遥感影像蚀变异常信息提取方法 |
-
2021
- 2021-12-31 CN CN202111660772.4A patent/CN114397251A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103383348A (zh) * | 2013-05-28 | 2013-11-06 | 吉林大学 | 植被覆盖区高光谱遥感蚀变矿物提取方法 |
CN108985130A (zh) * | 2017-05-31 | 2018-12-11 | 核工业北京地质研究院 | 一种基于缓冲区干扰信息去除的蚀变信息提取方法 |
CN108717044A (zh) * | 2018-05-24 | 2018-10-30 | 青海师范大学 | 一种去除植被覆盖影响的表层土壤含水量卫星遥感估算方法 |
CN109165580A (zh) * | 2018-08-09 | 2019-01-08 | 东华理工大学 | 高密度植被覆盖区条件下的遥感影像蚀变异常信息提取方法 |
Non-Patent Citations (3)
Title |
---|
杨会峰 等: "《青海省共和盆地地下水资源及其环境问题研究》", 30 September 2021, 中国地质大学出版社, pages: 185 * |
蒋巍;: "福州地区植被覆盖遥感估算及变化分析", 内江师范学院学报, no. 08, 15 August 2008 (2008-08-15), pages 81 - 85 * |
郭昱彤: "基于遥感影像的济南市植被覆盖度提取研究", 《绿色科技》, 29 February 2016 (2016-02-29), pages 11 - 12 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Masoud et al. | Mapping soil salinity using spectral mixture analysis of landsat 8 OLI images to identify factors influencing salinization in an arid region | |
Helmer et al. | Cloud-free satellite image mosaics with regression trees and histogram matching | |
Hill | Image segmentation for humid tropical forest classification in Landsat TM data | |
Koike et al. | Lineament analysis of satellite images using a Segment Tracing Algorithm (STA) | |
Ourhzif et al. | Lithological mapping using landsat 8 oli and aster multispectral data in imini-ounilla district south high atlas of marrakech | |
Crippen et al. | Unveiling the lithology of vegetated terrains in remotely sensed imagery | |
CN106971156A (zh) | 一种基于面向对象分类的稀土开采区遥感信息提取方法 | |
CN108986040B (zh) | 一种基于遥感多光谱影像的ndvi阴影影响去除方法 | |
CN102622738A (zh) | 一种Landsat TM/ETM+图像中山体阴影区的光谱信息恢复方法 | |
AU2021100848A4 (en) | A Regional Extraction Method of Ecological Restoration Project in the Grassland Based on the High-resolution Remote Sensing Images | |
CN114596495A (zh) | 一种基于Sentinel-2A遥感影像溜砂坡判识与自动提取方法 | |
CN111476723B (zh) | 一种Landsat-7扫描线纠正器失效的遥感图像丢失像素恢复方法 | |
Lei et al. | A novel algorithm of individual tree crowns segmentation considering three-dimensional canopy attributes using UAV oblique photos | |
CN111310771A (zh) | 遥感影像的道路图像提取方法、装置、设备及存储介质 | |
CN117152634B (zh) | 基于色度指数的多源卫星影像浮水植物识别方法和系统 | |
CN114397251A (zh) | 一种植被覆盖区铀矿化蚀变信息遥感识别方法 | |
CN109871774B (zh) | 一种基于局域相近像元的混合像元分解方法 | |
De Jong et al. | Imaging spectrometry for surveying and modelling land degradation | |
CN112734819B (zh) | 一种适用于高分辨率遥感卫星dsm的地表滤波方法 | |
Pande et al. | High-resolution and hyperspectral data fusion for classification | |
CN111160182B (zh) | 一种识别古河道地表钙结岩的方法 | |
CN112033914A (zh) | 一种基于遥感图像的彩钢瓦厂房提取方法 | |
Smara et al. | Multisource ERS-1 and optical data for vegetal cover assessment and monitoring in a semi-arid region of Algeria | |
Hu et al. | Rapid dehazing algorithm based on large-scale median filtering for high-resolution visible near-infrared remote sensing images | |
Zylshal et al. | An object based image analysis approach to semi-automated karst morphology extraction |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination |