CN109725311A - 一种基于雷达卫星提取地质构造的方法 - Google Patents
一种基于雷达卫星提取地质构造的方法 Download PDFInfo
- Publication number
- CN109725311A CN109725311A CN201811588325.0A CN201811588325A CN109725311A CN 109725311 A CN109725311 A CN 109725311A CN 201811588325 A CN201811588325 A CN 201811588325A CN 109725311 A CN109725311 A CN 109725311A
- Authority
- CN
- China
- Prior art keywords
- atural object
- image
- object image
- radar satellite
- extraction
- 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
Landscapes
- Geophysics And Detection Of Objects (AREA)
- Image Processing (AREA)
Abstract
本发明属于地质构造技术领域,具体公开一种基于雷达卫星提取地质构造的方法,该方法的具体步骤如下:步骤1、对地物影像的原始数据进行辐射定标,并对定标后的影像进行多视处理;步骤2、对地物影像进行滤波和增强处理;步骤3、对上述步骤2中逆差距后的地物影像进行评价,并进行地质构造提取。本发明的方法能够得出地质构造提取最终影像,实现快速、准确地地质构造提取;选择传统光学影像构造识别方法与本识别方法进行目视解译比较,经过本方法增强后的影像,岩体内部纹理特征清晰,线性构造特征明显,小岩体得到增强,与背景区分明显,利于与目视分析。
Description
技术领域
本发明属于地质构造技术领域,具体涉及一种基于雷达卫星提取地质构造的方法。
背景技术
现有技术中的地质构造流程是,首先获取研究区可见光数据,接着根据地质背景情况,绘制地质构造图,最后现场踏勘进一步提高构造解译精度。传统方法的不足体现在:首先,易受不利天气、山体阴影、水汽含量高等因素影响,导致大量数据无法使用;其次,不同工作人员由于个人主观经验的不同,导致构造解译结果具有很大的主观性,缺乏科学性;最后,由于人力物力的限制,踏勘无法遍布所有研究区,限制了提取的精度。
因此,迫切需要研究一种新的方法克服传统技术的弊端。
发明内容
针对现有技术的缺陷,本发明提供一种基于雷达卫星提取地质构造的方法,该方法能够得出地质构造提取最终影像,实现快速、准确地地质构造提取;选择传统光学影像构造识别方法与本识别方法进行目视解译比较,经过本方法增强后的影像,岩体内部纹理特征清晰,线性构造特征明显,小岩体得到增强,与背景区分明显,利于与目视分析。
实现本发明目的的技术方案:一种基于雷达卫星提取地质构造的方法,该方法的具体步骤如下:
步骤1、对地物影像的原始数据进行辐射定标,并对定标后的影像进行多视处理;
步骤2、对地物影像进行滤波和增强处理;
步骤3、对上述步骤2中逆差距后的地物影像进行评价,并进行地质构造提取。
所述的步骤1中对地物影像的原始数据进行辐射定标的公式如下:
σ0 i,j=10lg(DNi,j 2/Aσ 2)。
所述的步骤1中地物影像进行多视处理中,通过输入方位向视数和距离向视数,就可以得到多视处理后的地物影像。
所述的步骤1中地距分辨率=距离向分辨率/sin(中心入射角);方位向视数=地距分辨率/方位向分辨率。
所述的步骤1中地物影像头文件中记录了方位向分辨率、距离向分辨率及中心入射角。
所述的步骤2中对地物影像进行Lee滤波。
所述的步骤2中滤波窗口大小5×5。
所述的步骤2中采用灰度共生矩阵对地物影影像进行增强处理。
所述的步骤2中强参数选择3×3、5×5、7×7。
所述的步骤3中逆采用差距方法对地物影影像进行地质构造提取。
本发明的有益技术效果在于:本发明的方法通过对所获取的雷达卫星数据进行一系列预处理,定标后雷达数据,通过Lee滤波,去除噪声与雷达影像的“椒盐现象”,对雷达影像进行逆差距处理,根据提取结果综合目视解译,进行精度对比;根据精度评价结果,得出地质构造提取最终影像,实现快速、准确地地质构造提取。雷达卫星技术主动遥感、不受天气影响、不受山体阴影遮挡等优势,能够提高构造识别提取。本发明的方法基于雷达卫星的电磁波C波段,具有可穿透性特点,通过增强图像灰度值的空间相关特性,突出具有同样亮度或接近亮度的像素之间的位置分布特征。本发明的方法使用了Lee滤波对影像进行处理后,再对影像使用灰度共生矩阵实现图像增强,经过这两个处理之后地质构造体的线性构造得到了增强,与背景信息得到了分离,提高了影像的可用性,增强了解译效果。本发明的方法选择雷达数据,建立了地质构造提取的新技术,这一方法能够在光学影像受天气影响的情况下,增强线性地质构造,实现地质构造提取工作。
附图说明
图1为本发明所提供的一种基于雷达卫星提取地质构造的方法的流程图。
具体实施方式
下面结合附图和实施例对本发明提供的一种基于雷达卫星提取地质构造的方法作进一步详细说明。
如图1所示,本发明提供的一种基于雷达卫星提取地质构造的方法,该方法的具体步骤如下:
步骤1、为获取地物影像的后向散射系数,对地物影像的原始数据进行辐射定标,并对定标后的影像进行多视处理。
步骤1.1、对地物影像的原始数据进行辐射定标,辐射定标公式如下:
σ0 i,j=10lg(DNi,j 2/Aσ 2)
其中,σ0为地物影像的后向散射系数,i和j表示第i行,第j列;DN是雷达图像的灰度值,Aσ是定标参数,从雷达卫星数据的头文件中读取A,结合雷达图像灰度DN值,通过辐射定标公式,可以求得地物的后向散射系数σ0。
步骤1.2、对地物影像进行多视处理
地物影像头文件中记录了方位向分辨率、距离向分辨率及中心入射角。
地距分辨率=距离向分辨率/sin(中心入射角)
方位向视数=地距分辨率/方位向分辨率
多视处理中,通过输入方位向视数和距离向视数(默认1),就可以得到多视处理后的地物影像。
多视处理抑止了斑点噪声,同时提高了辐射分辨率,增强了图像信息,使图像更具可读性。
步骤2、对地物影像进行滤波和增强处理
步骤2.1、对地物影像进行Lee滤波
在雷达卫星获取的地物影像上,选择Lee滤波作为滤波方法,选择滤波窗口大小5×5,对地物影像进行Lee滤波
步骤2.2、对地物影像进行增强处理
在滤波后的雷达卫星获取的地物影像上,选择灰度共生矩阵对影像进行增强处理。通过灰度共生矩阵可使影像中地质体信息得到增强,增强参数选择3×3、5×5、7×7。
步骤3、对上述步骤2中逆差距后的地物影像进行评价,并进行地质构造提取
通过对比逆差距窗口3×3、5×5、7×7的影像,结合目视解译对影像进行评价,选择相对最优的结果进行地质构造提取。
选择使用逆差距方法对地物影影像进行地质构造提取。
上面结合附图和实施例对本发明作了详细说明,但是本发明并不限于上述实施例,在本领域普通技术人员所具备的知识范围内,还可以在不脱离本发明宗旨的前提下作出各种变化。本发明中未作详细描述的内容均可以采用现有技术。
Claims (10)
1.一种基于雷达卫星提取地质构造的方法,其特征在于,该方法的具体步骤如下:
步骤1、对地物影像的原始数据进行辐射定标,并对定标后的影像进行多视处理;
步骤2、对地物影像进行滤波和增强处理;
步骤3、对上述步骤2中逆差距后的地物影像进行评价,并进行地质构造提取。
2.根据权利要求1所述的一种基于雷达卫星提取地质构造的方法,其特征在于:所述的步骤1中对地物影像的原始数据进行辐射定标的公式如下:
σ0 i,j=10lg(DNi,j 2/Aσ 2)。
3.根据权利要求2所述的一种基于雷达卫星提取地质构造的方法,其特征在于:所述的步骤1中地物影像进行多视处理中,通过输入方位向视数和距离向视数,就可以得到多视处理后的地物影像。
4.根据权利要求3所述的一种基于雷达卫星提取地质构造的方法,其特在于:所述的步骤1中地距分辨率=距离向分辨率/sin(中心入射角);方位向视数=地距分辨率/方位向分辨率。
5.根据权利要求4所述的一种基于雷达卫星提取地质构造的方法,其特征在于:所述的步骤1中地物影像头文件中记录了方位向分辨率、距离向分辨率及中心入射角。
6.根据权利要求5所述的一种基于雷达卫星提取地质构造的方法,其特征在于:所述的步骤2中对地物影像进行Lee滤波。
7.根据权利要求6所述的一种基于雷达卫星提取地质构造的方法,其特征在于:所述的步骤2中滤波窗口大小5×5。
8.根据权利要求7所述的一种基于雷达卫星提取地质构造的方法,其特征在于:所述的步骤2中采用灰度共生矩阵对地物影影像进行增强处理。
9.根据权利要求8所述的一种基于雷达卫星提取地质构造的方法,其特征在于:所述的步骤2中强参数选择3×3、5×5、7×7。
10.根据权利要求9所述的一种基于雷达卫星提取地质构造的方法,其特征在于:所述的步骤3中逆采用差距方法对地物影影像进行地质构造提取。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811588325.0A CN109725311A (zh) | 2018-12-25 | 2018-12-25 | 一种基于雷达卫星提取地质构造的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811588325.0A CN109725311A (zh) | 2018-12-25 | 2018-12-25 | 一种基于雷达卫星提取地质构造的方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN109725311A true CN109725311A (zh) | 2019-05-07 |
Family
ID=66296345
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811588325.0A Pending CN109725311A (zh) | 2018-12-25 | 2018-12-25 | 一种基于雷达卫星提取地质构造的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109725311A (zh) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104517284A (zh) * | 2014-09-29 | 2015-04-15 | 西安电子科技大学 | 基于深度置信网的极化sar图像分割 |
CN104700110A (zh) * | 2015-04-03 | 2015-06-10 | 电子科技大学 | 一种基于全极化sar图像的植被覆盖信息提取方法 |
CN106447079A (zh) * | 2016-08-31 | 2017-02-22 | 贵州师范大学 | 一种基于Radarsat‑2的喀斯特山区烟草产量预测方法 |
CN106680794A (zh) * | 2015-11-10 | 2017-05-17 | 核工业北京地质研究院 | 基于空间建模技术的雷达数据快速去噪方法 |
CN107273813A (zh) * | 2017-05-23 | 2017-10-20 | 国家地理空间信息中心 | 基于高分卫星遥感数据的地理空间要素提取系统 |
US20180128925A1 (en) * | 2008-09-09 | 2018-05-10 | Geooptics, Inc. | Cellular Interferometer for Continuous Earth Remote Observation (CICERO) |
-
2018
- 2018-12-25 CN CN201811588325.0A patent/CN109725311A/zh active Pending
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20180128925A1 (en) * | 2008-09-09 | 2018-05-10 | Geooptics, Inc. | Cellular Interferometer for Continuous Earth Remote Observation (CICERO) |
CN104517284A (zh) * | 2014-09-29 | 2015-04-15 | 西安电子科技大学 | 基于深度置信网的极化sar图像分割 |
CN104700110A (zh) * | 2015-04-03 | 2015-06-10 | 电子科技大学 | 一种基于全极化sar图像的植被覆盖信息提取方法 |
CN106680794A (zh) * | 2015-11-10 | 2017-05-17 | 核工业北京地质研究院 | 基于空间建模技术的雷达数据快速去噪方法 |
CN106447079A (zh) * | 2016-08-31 | 2017-02-22 | 贵州师范大学 | 一种基于Radarsat‑2的喀斯特山区烟草产量预测方法 |
CN107273813A (zh) * | 2017-05-23 | 2017-10-20 | 国家地理空间信息中心 | 基于高分卫星遥感数据的地理空间要素提取系统 |
Non-Patent Citations (3)
Title |
---|
XIANGHUI YUAN ; TAO LIU: "《Texture invariant estimation of equivalent number of looks based on log-cumulants in polarimetric radar imagery》", 《 JOURNAL OF SYSTEMS ENGINEERING AND ELECTRONICS》 * |
武鼎: "《基于Lee滤波的SAR图像噪声滤除方法》", 《全国铀矿大基地建设学术研讨会论文集(下)》 * |
陈勇国: "《基于雷达遥感的冰川泥石流灾害判识与信息提取》", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109919875B (zh) | 一种高时频遥感图像特征辅助的居民地提取与分类方法 | |
CN111242224B (zh) | 一种基于无人机提取分类样本点的多源遥感数据分类方法 | |
Baccini et al. | A first map of tropical Africa’s above-ground biomass derived from satellite imagery | |
CN105989322B (zh) | 一种基于高分辨率遥感影像的多指标融合滑坡检测方法 | |
CN107145891B (zh) | 一种基于遥感影像的水体提取方法及系统 | |
CN113850139B (zh) | 一种基于多源遥感的森林年际物候监测方法 | |
CN108051371B (zh) | 一种面向生态环境参量遥感反演的阴影提取方法 | |
Chen et al. | Digital mapping of the soil thickness of loess deposits over a calcareous bedrock in central France | |
Goodarzi Mehr et al. | Using the mixture-tuned matched filtering method for lithological mapping with Landsat TM5 images | |
CN108319923A (zh) | 一种云阴影识别方法及系统 | |
CN105825217B (zh) | 基于活动轮廓模型的高光谱影像感兴趣区域自动提取方法 | |
Oommen et al. | Documenting earthquake-induced liquefaction using satellite remote sensing image transformations | |
Zhan et al. | Quantitative analysis of shadow effects in high-resolution images of urban areas | |
Gadi et al. | Understanding soil surface water content using light reflection theory: a novel color analysis technique considering variability in light intensity | |
CN109738369B (zh) | 一种采用航空高光谱遥感玉器的考古探测方法 | |
CN111428582A (zh) | 一种利用互联网街景照片计算城市天空开阔度的方法 | |
Alperovich et al. | A new combined wavelet methodology: Implementation to GPR and ERT data obtained in the Montagnole experiment | |
CN106569209A (zh) | 一种利用全极化雷达数据提取土壤含水信息的方法 | |
Piga et al. | Integrated geophysical and aerial sensing methods for archaeology: a case history in the Punic Site of Villamar (Sardinia, Italy) | |
CN109725311A (zh) | 一种基于雷达卫星提取地质构造的方法 | |
Alexakis et al. | Remote sensing applications in archaeological research | |
CN114118231B (zh) | 多源遥感浅覆盖区地层岩性分类特征空间构建方法、岩性分类方法及系统 | |
CN106680794A (zh) | 基于空间建模技术的雷达数据快速去噪方法 | |
Seid et al. | Identification of lithology and structures in Serdo, Afar, Ethiopia using remote sensing and Gis techniques | |
CN108230419B (zh) | 一种遥感图像光谱层析剥离方法及系统 |
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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190507 |