CN102231190B - 冲洪积扇信息的自动提取方法 - Google Patents
冲洪积扇信息的自动提取方法 Download PDFInfo
- Publication number
- CN102231190B CN102231190B CN 201110191094 CN201110191094A CN102231190B CN 102231190 B CN102231190 B CN 102231190B CN 201110191094 CN201110191094 CN 201110191094 CN 201110191094 A CN201110191094 A CN 201110191094A CN 102231190 B CN102231190 B CN 102231190B
- Authority
- CN
- China
- Prior art keywords
- image
- remote sensing
- ground object
- index
- 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.)
- Active
Links
- 238000000605 extraction Methods 0.000 title claims abstract description 34
- 238000000034 method Methods 0.000 claims abstract description 27
- 230000002194 synthesizing effect Effects 0.000 claims abstract description 4
- 230000009466 transformation Effects 0.000 claims abstract description 4
- 238000002310 reflectometry Methods 0.000 claims description 29
- 230000011218 segmentation Effects 0.000 claims description 8
- 238000001914 filtration Methods 0.000 claims description 6
- 239000002131 composite material Substances 0.000 claims description 3
- 238000012937 correction Methods 0.000 claims description 3
- 238000005260 corrosion Methods 0.000 claims description 3
- 230000007797 corrosion Effects 0.000 claims description 3
- 230000000877 morphologic effect Effects 0.000 claims description 3
- 238000000638 solvent extraction Methods 0.000 abstract 1
- 230000000007 visual effect Effects 0.000 description 4
- 238000011835 investigation Methods 0.000 description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 3
- 239000000284 extract Substances 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 239000007787 solid Substances 0.000 description 2
- 238000005267 amalgamation Methods 0.000 description 1
- 210000004556 brain Anatomy 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000019771 cognition Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 239000011499 joint compound Substances 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 230000007480 spreading Effects 0.000 description 1
- 239000004575 stone Substances 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
Images
Landscapes
- Investigating Or Analysing Materials By Optical Means (AREA)
- Image Processing (AREA)
Abstract
本发明公开了一种冲洪积扇信息的自动提取方法,包括如下步骤:(1)计算遥感影像中地物的遥感提取指数、植被指数和湿度指数,并分别生成遥感提取指数灰度图像、植被指数灰度图像和湿度指数灰度图像;(2)将所述遥感提取指数灰度图像、所述植被指数灰度图像和所述湿度指数灰度图像合成彩色合成图像;(3)对所述彩色合成图像进行主成分变换并提取第一主成分生成第一主成分灰度图像;(4)对所述第一主成分灰度图像进行阈值分割,得到冲洪积扇的二值图像,从而提取出冲洪积扇信息。本发明能准确将冲洪积扇信息提取出来,且提取的精度高、提取的速度快。
Description
技术领域
本发明涉及地物信息遥感提取领域,尤其涉及一种冲洪积扇信息的自动提取方法。
背景技术
冲洪积扇包括冲积扇和洪积扇,它们是由常年性河流、暂时性山地水流携带的沙、石、泥等固体物质在出山后,以山麓谷口为顶点向开阔地展布而形成的扇状堆积地貌,是冲积平原的组成部分。冲洪积扇在水资源、环境调查和评价中具有重要的应用和意义,同时冲洪积扇由于成因、组成颗粒的分选程度不同,因此在工程地质调查中非常重要。冲洪积扇是一个空间地理实体,其发育范围多在几十公里甚至更大,单靠地面调查研究其演变过程比较困难,利用遥感影像从宏观角度观察和分析其空间演变具有很好的效果。
目前,对冲洪积扇的遥感提取方法主要是以人工目视解译为主,辅以人机交互分类的方法。人工目视解译是一种人工提取信息的方法,使用眼睛目视观察,同时借助一些光学仪器火灾计算机显示屏幕上,凭借丰富的解译经验,扎实的专业知识和手头的相关资料,通过人脑的分析、推理和判断,提取有用的信息。
然而,一方面,由于采用人工目视解译的方法解译遥感图像是一个复杂的认知过程,往往需要经历几次反复解译才能得到正确的结果,一景数据多种灾害信息的解译往往需要几天甚至更长的时间,因而,该方法的效率比较低。另一方面,采用人工目视解译的方法解译遥感图像时,只有在遥感图像上存在着颜色或者明度的差异时,并且这种差异能为解译者视觉所区分,才有可能将地物目标与背景区别开,且这种方法的解译精度因人而异,会成解译成果的拼接,因而,该方法的译精度比较低。
因此,有必要提供一种冲洪积扇信息的自动提取方法克服上述缺陷。
发明内容
本发明的目的是提供一种冲洪积扇信息的自动提取方法,能准确将冲洪积扇信息提取出来,且提取的精度高、提取的速度快。
为了实现上述目的,本发明提供了一种冲洪积扇信息的自动提取方法,包括如下步骤:(1)获取通过卫星检测的地物的红光波长、蓝光波长、近红外波长和绿光波长,并根据公式计算遥感影像中地物的遥感提取指数并生成遥感提取指数灰度图像,其中,IAFI为遥感提取指数,λRED为遥感影像中地物的红光波长,λNIR为遥感影像中地物的近红外波长,λBLUE为遥感影像中地物的蓝光波长,λCREEN为遥感影像中地物的绿光波长;(2)获取通过卫星检测的地物的近红外波段反射率和红光波段反射率,并根据公式计算遥感影像中地物的归一化植被指数并生成植被指数灰度图像,其中,INDVI为归一化植被指数,ρNIR为遥感影像中地物的近红外波段反射率,ρRED为遥感影像中地物的红光波段反射率;(3)获取通过卫星检测的地物的中红外波段反射率和绿光波段反射率,并根据公式计算遥感影像中地物的湿度指数并生成湿度指数灰度图像,其中,INDMI为遥感影像中地物的湿度指数,ρGREEN为遥感影像中地物的绿光波段反射率,ρMIR为遥感影像中地物的中红外波段反射率;(4)将所述遥感提取指数灰度图像、所述植被指数灰度图像和所述湿度指数灰度图像合成彩色合成图像;(5)对所述彩色合成图像进行主成分变换并提取第一主成分生成第一主成分灰度图像;(6)对所述第一主成分灰度图像进行灰度值拉伸;(7)获取拉伸后的所述第一主成分灰度图像的灰度直方图,通过迭代法计算所述灰度直方图的最大灰度值和最小灰度值,并选择两个最大灰度值之间的最小灰度值作为分割阈值;(8)根据所述分割阈值,将所述第一主成分灰度图像进行二值化,得到冲洪积扇的二值图像,从而提取出冲洪积扇信息。
较佳地,所述步骤(1)之前还包括:对通过卫星获得的遥感影像的五个波段进行大气校正预处理。
较佳地,所述卫星为TM卫星,遥感影像的所述五个波段为近红外波段、红光波段、绿光波段、蓝光波段和中红外波段,遥感影像为冲积扇的遥感影像或者洪积扇的遥感影像。
较佳地,所述步骤(6)具体为:将所述第一主成分灰度图像中的灰度值范围拉伸到0-255。
较佳地,所述步骤(8)之后还包括:对所述冲洪积扇的二值图像进行形态学膨胀滤波、腐蚀滤波以及去噪处理。
与现有技术相比,本发明的方法通过计算遥感图像中的遥感提取指数、植被指数和湿度指数,从而增强了遥感图像中的冲洪积扇的灰度值,使冲洪积扇的灰度值与其他地物的灰度值的差异变大,消除了地形阴影的影响,因而,可通过阈值自动选取算法精确地选取分割阈值,从而将冲洪积扇信息准确地提取出来。另外,由于冲洪积扇信息的提取都是基于多峰直方图阈值自动选取算法实现自动提取的,因而,提取的精度高且提取的速度快。
通过以下的描述并结合附图,本发明将变得更加清晰,这些附图用于解释本发明的实施例。
附图说明
图1为本发明冲洪积扇信息的自动提取方法的主流程图。
图2为经过本发明冲洪积扇信息的自动提取方法的步骤S1-步骤S6处理后得到的冲积扇的第一主成分灰度图像。
图3为采用本发明冲洪积扇信息的自动提取方法的得到的冲积扇的二值图像。
图4为采用本发明冲洪积扇信息的自动提取方法的得到的洪积扇的二值图像。
具体实施方式
现在参考附图描述本发明的实施例,附图中类似的元件标号代表类似的元件。
如图1,本实施例冲洪积扇信息的自动提取方法包括如下步骤:
步骤S1,通过TM卫星拍摄目标地的冲积扇的遥感影像或者洪积扇的遥感影像,并对所述冲积扇的遥感影像或者洪积扇的遥感影像的近红外波段、红光波段、绿光波段、蓝光波段和中红外波段进行大气校正预处理;
其中,IAFI为遥感提取指数,λRED为遥感影像中地物的红光波长,λNIR为遥感影像中地物的近红外波长,λBLUE为遥感影像中地物的蓝光波长,λCREEN为遥感影像中地物的绿光波长;
其中,INDVI为归一化植被指数,ρNIR为遥感影像中地物的近红外波段反射率,ρRED为遥感影像中地物的红光波段反射率;
其中,INDMI为遥感影像中地物的湿度指数,ρGREEN为遥感影像中地物的绿光波段反射率,ρMIR为遥感影像中地物的中红外波段反射率;
步骤S5,将所述遥感提取指数灰度图像、所述植被指数灰度图像和所述湿度指数灰度图像合成彩色合成图像;
步骤S6,对所述彩色合成图像进行主成分变换并提取第一主成分生成第一主成分灰度图像;
步骤S7,将所述第一主成分灰度图像中的灰度值范围拉伸到0-255;
步骤S8,提取拉伸后的所述第一主成分灰度图像的灰度直方图,通过迭代法计算所述灰度直方图的最大灰度值和最小灰度值,并选择两个最大灰度值之间的最小灰度值作为分割阈值;
步骤S9,根据所述分割阈值,将所述第一主成分灰度图像进行二值化,得到冲洪积扇的二值图像,从而提取出冲洪积扇信息;
步骤S10,对所述冲洪积扇的二值图像进行形态学膨胀滤波、腐蚀滤波以及去噪处理。
如图2-3所示,冲积扇的遥感影像经过步骤S1-步骤S6的处理后,冲积扇与植被、水体、耕地及人工建筑等地物的光谱差异进一步增大。因而,利用阈值自动选取算法可简单、准确的提取冲积扇的图斑。提取的结果中会存在极少的道路,在后续计算中加入地形数据,可极为方便的剔除,不影响提取的精度。
如图4所示,洪积扇的遥感影像经过本发明的冲洪积扇信息的自动提取方法处理后,通过最终获得的洪积扇的二值图像,洪积扇的信息被精确地提取出。
由上可知,本发明的方法通过计算遥感图像中的遥感提取指数、植被指数和湿度指数,从而增强了遥感图像中的冲洪积扇的灰度值,使冲洪积扇的灰度值与其他地物的灰度值的差异变大,消除了地形阴影的影响,因而,可通过阈值自动选取算法精确地选取分割阈值,从而将冲洪积扇信息准确地提取出来。另外,由于冲洪积扇信息的提取都是基于多峰直方图阈值自动选取算法实现自动提取的,因而,提取的精度高且提取的速度快。
以上结合最佳实施例对本发明进行了描述,但本发明并不局限于以上揭示的实施例,而应当涵盖各种根据本发明的本质进行的修改、等效组合。
Claims (5)
1.一种冲洪积扇信息的自动提取方法,包括如下步骤:
其中,IAFI为遥感提取指数,λRED为遥感影像中地物的红光波长,λNIR为遥感影像中地物的近红外波长,λBLUE为遥感影像中地物的蓝光波长,λCREEN为遥感影像中地物的绿光波长;
其中,INDVI为归一化植被指数,ρNIR为遥感影像中地物的近红外波段反射率,ρRED为遥感影像中地物的红光波段反射率;
其中,INDMI为遥感影像中地物的湿度指数,ρGREEN为遥感影像中地物的绿光波段反射率,ρMIR为遥感影像中地物的中红外波段反射率;
(4)将所述遥感提取指数灰度图像、所述植被指数灰度图像和所述湿度指数灰度图像合成彩色合成图像;
(5)对所述彩色合成图像进行主成分变换并提取第一主成分生成第一主成分灰度图像;
(6)对所述第一主成分灰度图像进行灰度值拉伸;
(7)提取拉伸后的所述第一主成分灰度图像的灰度直方图,通过迭代法计算所述灰度直方图的最大灰度值和最小灰度值,并选择两个最大灰度值之间的最小灰度值作为分割阈值;
(8)根据所述分割阈值,将所述第一主成分灰度图像进行二值化,得到冲洪积扇的二值图像,从而提取出冲洪积扇信息。
2.如权利要求1所述的冲洪积扇信息的自动提取方法,其特征在于,所述步骤(1)之前还包括:
对通过卫星获得的遥感影像的五个波段进行大气校正预处理。
3.如权利要求2所述的冲洪积扇信息的自动提取方法,其特征在于,所述卫星为TM卫星,遥感影像的所述五个波段为近红外波段、红光波段、绿光波段、蓝光波段和中红外波段,遥感影像为冲积扇的遥感影像或者洪积扇的遥感影像。
4.如权利要求1所述的冲洪积扇信息的自动提取方法,其特征在于,所述步骤(6)具体为:
将所述第一主成分灰度图像中的灰度值范围拉伸到0-255。
5.如权利要求1所述的冲洪积扇信息的自动提取方法,其特征在于,所述步骤(8)之后还包括:
对所述冲洪积扇的二值图像进行形态学膨胀滤波、腐蚀滤波以及去噪处理。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110191094 CN102231190B (zh) | 2011-07-08 | 2011-07-08 | 冲洪积扇信息的自动提取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110191094 CN102231190B (zh) | 2011-07-08 | 2011-07-08 | 冲洪积扇信息的自动提取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102231190A CN102231190A (zh) | 2011-11-02 |
CN102231190B true CN102231190B (zh) | 2012-10-31 |
Family
ID=44843753
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201110191094 Active CN102231190B (zh) | 2011-07-08 | 2011-07-08 | 冲洪积扇信息的自动提取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102231190B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105427305B (zh) * | 2015-11-19 | 2018-12-21 | 国家***北海预报中心 | 一种绿潮信息提取方法 |
CN105629325B (zh) * | 2015-12-23 | 2018-03-16 | 中国石油大学(北京) | 前陆盆地冲积扇精细刻画与预测方法 |
CN108154068A (zh) * | 2016-12-02 | 2018-06-12 | 航天星图科技(北京)有限公司 | 一种遥感影像中泥石流信息的提取方法 |
CN107256452A (zh) * | 2017-05-23 | 2017-10-17 | 国家地理空间信息中心 | 基于高分卫星遥感数据的城镇化水平评价*** |
CN108875615B (zh) * | 2018-06-07 | 2021-04-30 | 中国石油天然气股份有限公司 | 沉积区域遥感识别方法、装置、电子设备及存储介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1710380A (zh) * | 2005-07-01 | 2005-12-21 | 武汉大学 | 基于影像分割技术的立体影像行树三维信息提取方法 |
CN1924610A (zh) * | 2005-09-01 | 2007-03-07 | 中国林业科学研究院资源信息研究所 | 利用陆地卫星数据反演森林郁闭度和蓄积量的方法 |
CN102032903A (zh) * | 2010-12-01 | 2011-04-27 | 南京大学 | 基于Landsat数据源的珊瑚岛礁遥感信息自动提取方法 |
CN102054274A (zh) * | 2010-12-01 | 2011-05-11 | 南京大学 | 一种海岸带水体遥感信息全自动提取方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006266858A (ja) * | 2005-03-24 | 2006-10-05 | Railway Technical Res Inst | 硬質堆積岩地域の災害発生危険度評価方法 |
JP5481140B2 (ja) * | 2009-09-28 | 2014-04-23 | 株式会社日立ソリューションズ | スペクトル解析装置 |
-
2011
- 2011-07-08 CN CN 201110191094 patent/CN102231190B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1710380A (zh) * | 2005-07-01 | 2005-12-21 | 武汉大学 | 基于影像分割技术的立体影像行树三维信息提取方法 |
CN1924610A (zh) * | 2005-09-01 | 2007-03-07 | 中国林业科学研究院资源信息研究所 | 利用陆地卫星数据反演森林郁闭度和蓄积量的方法 |
CN102032903A (zh) * | 2010-12-01 | 2011-04-27 | 南京大学 | 基于Landsat数据源的珊瑚岛礁遥感信息自动提取方法 |
CN102054274A (zh) * | 2010-12-01 | 2011-05-11 | 南京大学 | 一种海岸带水体遥感信息全自动提取方法 |
Non-Patent Citations (3)
Title |
---|
JP特开2006-266858A 2006.10.05 |
JP特开2011-69757A 2011.04.07 |
杨树文等.一种利用TM图像自动提取洪积扇的方法.《国土资源遥感》.2011,(第2期),第65-69页. * |
Also Published As
Publication number | Publication date |
---|---|
CN102231190A (zh) | 2011-11-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
KR101728137B1 (ko) | 위성 영상과 gis를 사용한 토지피복 항목별 영상분류 방법 | |
CN111932567B (zh) | 基于卫星影像的冰湖轮廓自动提取方法 | |
CN108051371B (zh) | 一种面向生态环境参量遥感反演的阴影提取方法 | |
Yang et al. | River delineation from remotely sensed imagery using a multi-scale classification approach | |
CN102231190B (zh) | 冲洪积扇信息的自动提取方法 | |
CN113033670A (zh) | 一种基于Sentinel-2A/B数据的水稻种植面积提取方法 | |
CN105139396B (zh) | 一种全自动遥感影像云雾检测方法 | |
Aliero et al. | The usefulness of unmanned airborne vehicle (UAV) imagery for automated palm oil tree counting | |
CN112907587B (zh) | 一种基于GEE的Otsu和边缘检测算法的高山林线提取方法 | |
CN113077486A (zh) | 一种山区植被覆盖率监测方法及*** | |
Siregar et al. | Mapping shallow waters habitats using OBIA by applying several approaches of depth invariant index in North Kepulauan Seribu | |
Cruz et al. | Precise real-time detection of nonforested areas with UAVs | |
Kadhim et al. | Shadow detection from very high resoluton satellite image using grabcut segmentation and ratio-band algorithms | |
Zhao et al. | Image dehazing based on haze degree classification | |
CN107507151B (zh) | 一种多光谱遥感图像真实色彩还原方法及*** | |
Nanare et al. | Remote sensing satellite image analysis for deforestation in yavatmal district, maharashtra, india | |
CN113469104A (zh) | 基于深度学习的雷达遥感影像地表水体变化检测方法和设备 | |
Raja | Tree crown delineation from high resolution satellite images | |
CN109544558B (zh) | 一种城市复杂环境***影与水体分离方法 | |
Karydas et al. | Fine scale mapping of agricultural landscape features to be used in environmental risk assessment in an olive cultivation area | |
Reddy et al. | A novel method for water and water canal extraction from landsat-8 oli imagery | |
Bhaskaran et al. | Mapping shadows in very high-resolution satellite data using HSV and edge detection techniques | |
CN113128523B (zh) | 一种基于时间序列遥感影像自动提取珊瑚礁的方法 | |
Bouhennache et al. | Using landsat images for urban change detection, a case study in Algiers town | |
Hamid et al. | Low-illuminated SPOT-5 image improvement for density-based vegetation identification using three-layer colour manipulation approach |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |