CN104007443A - 一种利用光谱特征提取极地海冰融池的方法 - Google Patents

一种利用光谱特征提取极地海冰融池的方法 Download PDF

Info

Publication number
CN104007443A
CN104007443A CN201410169032.4A CN201410169032A CN104007443A CN 104007443 A CN104007443 A CN 104007443A CN 201410169032 A CN201410169032 A CN 201410169032A CN 104007443 A CN104007443 A CN 104007443A
Authority
CN
China
Prior art keywords
melt pool
ice
snow
water body
polar region
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
Application number
CN201410169032.4A
Other languages
English (en)
Inventor
柯长青
邓娟
林辉
马东辉
李创
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing University
Original Assignee
Nanjing University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Nanjing University filed Critical Nanjing University
Priority to CN201410169032.4A priority Critical patent/CN104007443A/zh
Publication of CN104007443A publication Critical patent/CN104007443A/zh
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S17/00Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
    • G01S17/88Lidar systems specially adapted for specific applications
    • G01S17/89Lidar systems specially adapted for specific applications for mapping or imaging
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/48Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S17/00

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Electromagnetism (AREA)
  • Image Processing (AREA)

Abstract

本发明涉及一种利用光谱特征提取极地海冰融池的方法,属于遥感地学应用领域。以MODIS反射率产品MOD09A1为数据源,先进行投影转换、重采样、赋予坐标信息等预处理。再根据水体、冰雪与融池波谱曲线的差异,建立反射率与三者覆盖率之间的多元线性方程组。采用拟牛顿法(Quasi-NewtonMethod)解算此超定线性方程组。最后输出的结果为水体、冰雪和融池在各个像元上所占的比例。冰雪的比例是相对于水体以外的部分,而融池分布于冰雪表面。所以,总的比例1减去水体的比例,再乘以融池的计算比例,就得到融池覆盖的区域,进而提取出海冰融池。

Description

一种利用光谱特征提取极地海冰融池的方法
技术领域
本发明涉及一种提取极地海冰融池的方法,属于遥感地学应用技术领域。 
背景技术
极地海冰的融池会降低海冰表面反射率,并且会加快极地海冰的融化。夏季北极海冰表面有50%-60%的融池。一年冰覆盖的平坦区域,融池甚至达到了90%。当吸收的太阳短波辐射增加,海冰表面气温大于冰点时,海冰开始融化,在表面形成融池。融池使海冰反射率从大约0.8降到0.5。融池的出现极大地影响到海冰、大气与海水之间的辐射平衡,同时对海冰反射回馈的强度影响也很大。因此,极地海冰融池的空间分布研究有重要意义。 
在北极不同地区开展的众多海冰表面观测和船基走航观测,其目的是研究融池的反射率、光谱特征以及融池的大小和分布等。由俄罗斯绘制的海冰图集就是基于航空和海冰表面观测的北极部分融池分布图。 
不同表面类型的总体光谱和反射率可以从光谱曲线库中查到。水体的值在0.06左右,融池为0.29,而冰雪为0.87。 
融冰期与地表吸收的太阳辐射有关。在北冰洋南部,如白令海和哈德逊湾,4月中旬可以观测到海冰融化过程。在北极地区中部,6月开始有海冰融化,6月中旬,海冰表面已经有很多融池。 
融池是随着雪的融化开始的,冰或雪融化后的水聚集在表面凹陷处或者其它地表形变处。相比于多年冰不规则的表面形态,表面平整的一年冰更易于形成大的融池。多年冰上形成的融池面积较小,但个数多,而一年冰上形成的融池面积较大,个数少。由于水的对流,其热量传递要超过冰的热量传递。由于融池的反射率低,也使得更多的热量会被融池周围的冰吸收。这两方面的因素会导致融池周围冰的融化速度是裸冰融化速度的2-3倍。这样,融池会逐渐变深,直到融透整个冰层。但是随着融池深度的增加,融池的直径逐步减小。融池中的水温是高于冰点的,随着热水的下沉,热量传递会将底部的冰融化。在 每年的8月末到9月初,随着气温的降低,融池中的水会重新冻结,融池也逐渐减少。到北极降雪时,融池基本被雪覆盖,这时海冰表面反射率会增加。 
发明内容
本发明的目的是为提取极地海冰融池的覆盖率而提出的,一种过程简单、运算速度快的融池提取方法。 
本发明利用光谱特征提取极地海冰融池的方法,其步骤包括: 
第一步、将MODIS遥感影像的MOD09A1反射率数据由正弦曲线投影转换为极地方位投影; 
第二步、从完成投影转换的影像中提取目标区域的图像数据,并提取每个像元在λ3、λ1、λ2波段的反射率值; 
第三步、根据水体、冰雪和融池的光谱曲线,得到水体、冰雪和融池分别在λ3、λ1、λ2波段的光谱反射率; 
第四步、针对目标区域内的每个像元,根据光谱反射率与水体、冰雪和融池三者覆盖率的关系,建立多元线性方程组,并求解获得水体、融池、冰雪分别在像元内所占的百分比Sw、Sm、Si; 
第五步、根据公式计算得到目标区域内所有像元的融池覆盖比例; 
第六步、以像元的融池覆盖比例为要素对目标区像元进行渲染制图,获得目标区融池比例分布图。 
本发明方法使用的数据源为中分辨率成像光谱仪数据(MODIS),MODIS数据包括36个波段,波长为0.4-14.4μm,空间分辨率为250-1000m。根据光谱特征,水体、冰雪和融池在MODIS数据中可以分离。将MODIS的象元值归一化到0至1之间,计算得到融池的覆盖度。 
本发明利用光谱特征提取极地海冰融池的方法,还具有如下进一步的特征: 
1、所述第一步投影转换前,先查看MODIS反射率产品MOD09A1影像的地理 定位等信息,在遥感图像处理软件中对头文件信息进行适应性修改,将原影像的正弦曲线投影转换为极地方位投影。 
2、所述λ3波段的波长范围为459-479nm,λ1波段的波长范围为620-670nm,λ2波段的波长范围为841-876nm。 
3、所述第四步中,所述多元线性方程组为: 
Swrw3)+Smrm3)+Siri3)=R(λ3
Swrw1)+Smrm1)+Siri1)=R(λ1
Swrw2)+Smrm2)+Siri2)=R(λ2
Sw+Sm+Si=1 
式中,R(λ3)、R(λ1)、R(λ2)分别是像元在λ3、λ1、λ2波段的反射率值。 
4、所述第四步中,使用拟牛顿法求解方程组,得到各个像元中水体、融池和冰雪分别所占的比例。 
5、第六步制作融池比例分布图时,将融池覆盖比例划分为若干区间,根据像元的融池覆盖比例所在区间进行不同灰度等级的渲染。 
6、在λ3波段:水体的光谱反射率rw3)=0.08,冰雪的光谱反射率ri3)=0.95,融池的光谱反射率rm3)=0.22;在λ1波段:水体的光谱反射率rw1)=0.08,冰雪的光谱反射率ri1)=0.95,融池的光谱反射率rm1)=0.16;在λ2波段:水体的光谱反射率rw2)=0.08,冰雪的光谱反射率ri1)=0.87,融池的光谱反射率rm1)=0.07。 
本发明的有益效果是: 
融池是影响夏季极地海冰反射率的一个重要因素,但融池的时空分布情况尚不完全清楚。所以获取融池的时空分布信息是一项有意义的研究工作。本发明利用水体、冰雪和融池光谱特征的差异,计算得出融池的覆盖范围。 
第一,本发明方法操作简单,通过下载免费的MODIS数据,就可以实现对极地海冰融池的研究。 
第二,通过解算超定方程组得到融池的范围。该方法有两重意义,一是通过MODIS数据提取融池,可与现场观测数据进行对比,相互验证,进而获取大范围可靠的融池覆盖信息;二是提取的大范围融池覆盖信息,可以为极地气候变化研究提供基础数据。 
附图说明
下面结合附图对本发明作进一步的说明: 
图1本发明融池提取流程图。 
图2在完整的一景影像中选取大小为1000×1000个象元的实验区。 
图3提取1000×1000个象元的图像数据。 
图4水体、海冰和融池的光谱曲线。 
图5融池覆盖率图。 
具体实施方式
本实施例所利用的MODIS反射率产品MOD09A1影像从网站下载获得,其经纬度范围为:0°W-180°W,70°N-90°N,时间为2012年7月3日(儒略日第185天),下载网址为:http://ladsweb.nascom.nasa.gov/data/search.html。 
如图1所示,为本实施例提取极地海冰融池的方法流程图,具体包括如下步骤: 
第一步、将MODIS遥感影像的MOD09A1反射率数据由正弦曲线投影转换为极地方位投影。 
MOD09A1反射率数据是HDF格式,正弦曲线投影,但是很多遥感图像处理软件在打开HDF文件时,会有投影丢失现象。因此本例中利用遥感图像处理软件 查看HDF文件中的地理定位等信息(包含投影信息以及投影参数),在遥感图像处理软件中对头文件信息进行适应性修改,将原影像的正弦曲线投影转换为极地方位投影。 
第二步、从完成投影转换的影像中提取目标区域的图像数据,并提取每个像元在λ3、λ1、λ2波段的反射率值。 
如图2所示,选取1000×1000个象元的图像数据作为目标区域(图中矩形框内区域)提取海冰融池,该目标区域作为感兴趣区从整景影像中截取下来(如图3所示),再分别提出每个像元在λ3、λ1、λ2波段的的反射率值。 
λ3波段的波长范围为459-479nm,λ1波段的波长范围为620-670nm,λ2波段的波长范围为841-876nm。 
第三步、在遥感图像处理软件的光谱库中可以查看海冰、水体和融池的光谱曲线(如图4所示),L1为水体的光谱曲线,L2为冰雪的光谱曲线,L3为融池光谱曲线。根据水体、冰雪和融池的光谱曲线,得到水体、冰雪和融池分别在λ3、λ1、λ2波段的光谱反射率: 
在λ3波段:水体的光谱反射率rw3)=0.08,冰雪的光谱反射率ri3)=0.95,融池的光谱反射率rm3)=0.22; 
在λ1波段:水体的光谱反射率rw1)=0.08,冰雪的光谱反射率ri1)=0.95,融池的光谱反射率rm1)=0.16; 
在λ2波段:水体的光谱反射率rw2)=0.08,冰雪的光谱反射率ri1)=0.87,融池的光谱反射率rm1)=0.07。 
第四步、针对目标区域内的每个像元,根据光谱反射率与水体、冰雪和融池三者覆盖率的关系,建立多元线性方程组,并求解获得水体、融池、冰雪分别在像元内所占的百分比Sw、Sm、Si: 
Swrw3)+Smrm3)+Siri3)=R(λ3
Swrw1)+Smrm1)+Siri1)=R(λ1
Swrw2)+Smrm2)+Siri2)=R(λ2
Sw+Sm+Si=1 
式中,R(λ3)、R(λ1)、R(λ2)分别是像元在λ3、λ1、λ2波段的反射率值,该反射率值在第二步中已经获得。 
上述方程组为超定方程组,本例中使用拟牛顿法求解方,得到各个像元中水体、融池和冰雪分别所占的比例。 
第五步、融池的比例是相当于冰雪表面来说的,所以1减去水体的比例得到冰雪表面的比例,再乘以第五步中的融池比例,就得到了融池的覆盖范围。因此,本例中根据公式计算得到目标区域内所有像元的融池覆盖比例; 
第六步、以像元的融池覆盖比例为要素对目标区像元进行渲染制图,获得目标区融池比例分布图。 
本步骤制作融池比例分布图时,将融池覆盖比例划分为5个区间:0-0.2,0.2-0.4,0.4-0.6,0.6-0.8,0.8-1.0,根据像元的融池覆盖比例所在区间进行不同灰度等级的渲染,最终获得的目标区融池比例分布图(图5)。 
除上述实施例外,本发明还可以有其他实施方式。凡采用等同替换或等效变换形成的技术方案,均落在本发明要求的保护范围。 

Claims (7)

1.一种利用光谱特征提取极地海冰融池的方法,其步骤包括: 
第一步、将MODIS遥感影像的MOD09A1反射率数据由正弦曲线投影转换为极地方位投影; 
第二步、从完成投影转换的影像中提取目标区域的图像数据,并提取每个像元在λ3、λ1、λ2波段的反射率值; 
第三步、根据水体、冰雪和融池的光谱曲线,得到水体、冰雪和融池分别在λ3、λ1、λ2波段的光谱反射率; 
第四步、针对目标区域内的每个像元,根据光谱反射率与水体、冰雪和融池三者覆盖率的关系,建立多元线性方程组,并求解获得水体、融池、冰雪分别在像元内所占的百分比Sw、Sm、Si; 
第五步、根据公式计算得到目标区域内所有像元的融池覆盖比例; 
第六步、以像元的融池覆盖比例为要素对目标区像元进行渲染制图,获得目标区融池比例分布图。 
2.根据权利要求1所述的利用光谱特征提取极地海冰融池的方法,其特征在于:所述第一步投影转换前,先查看MODIS反射率产品MOD09A1影像的地理定位等信息,在遥感图像处理软件中对头文件信息进行适应性修改,将原影像的正弦曲线投影转换为极地方位投影。 
3.根据权利要求1所述的利用光谱特征提取极地海冰融池的方法,其特征在于:所述λ3波段的波长范围为459-479nm,λ1波段的波长范围为620-670nm,λ2波段的波长范围为841-876nm。 
4.根据权利要求1所述的利用光谱特征提取极地海冰融池的方法,其特征在于:第四步中,所述多元线性方程组为: 
Swrw3)+Smrm3)+Siri3)=R(λ3
Swrw1)+Smrm1)+Siri1)=R(λ1
Swrw2)+Smrm2)+Siri2)=R(λ2
Sw+Sm+Si=1 
式中,R(λ3)、R(λ1)、R(λ2)分别是像元在λ3、λ1、λ2波段的反射率值。 
5.根据权利要求4所述的利用光谱特征提取极地海冰融池的方法,其特征在于:所述第四步中,使用拟牛顿法求解方程组,得到各个像元中水体、融池和冰雪分别所占的比例。 
6.根据权利要求1所述的利用光谱特征提取极地海冰融池的方法,其特征在于:第六步制作融池比例分布图时,将融池覆盖比例划分为若干区间,根据像元的融池覆盖比例所在区间进行不同灰度等级的渲染。 
7.根据权利要求1所述的利用光谱特征提取极地海冰融池的方法,其特征在于:在λ3波段:水体的光谱反射率rw3)=0.08,冰雪的光谱反射率ri3)=0.95,融池的光谱反射率rm3)=0.22;在λ1波段:水体的光谱反射率rw1)=0.08,冰雪的光谱反射率ri1)=0.95,融池的光谱反射率rm1)=0.16;在λ2波段:水体的光谱反射率rw2)=0.08,冰雪的光谱反射率ri1)=0.87,融池的光谱反射率rm1)=0.07。 
CN201410169032.4A 2014-04-24 2014-04-24 一种利用光谱特征提取极地海冰融池的方法 Pending CN104007443A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410169032.4A CN104007443A (zh) 2014-04-24 2014-04-24 一种利用光谱特征提取极地海冰融池的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410169032.4A CN104007443A (zh) 2014-04-24 2014-04-24 一种利用光谱特征提取极地海冰融池的方法

Publications (1)

Publication Number Publication Date
CN104007443A true CN104007443A (zh) 2014-08-27

Family

ID=51368187

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410169032.4A Pending CN104007443A (zh) 2014-04-24 2014-04-24 一种利用光谱特征提取极地海冰融池的方法

Country Status (1)

Country Link
CN (1) CN104007443A (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105116464A (zh) * 2015-08-12 2015-12-02 南京大学 基于神经网络模型的极地海冰融池提取方法
CN106290183A (zh) * 2016-07-19 2017-01-04 国网天津市电力公司 基于光谱反射比的雪状态参数测量方法
CN110866363A (zh) * 2019-11-12 2020-03-06 南京信息工程大学 一种利用人工神经网络反演北极融池分布的方法
CN116597319A (zh) * 2023-04-19 2023-08-15 北京师范大学 一种海冰日数据综合处理方法、装置及计算设备

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040164235A1 (en) * 2003-01-21 2004-08-26 Miller Steven D. System and method for significant dust detection and enhancement of dust images over land and ocean
CN102636779A (zh) * 2012-05-07 2012-08-15 武汉大学 一种基于重采样回归分析的亚像素积雪覆盖率提取方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040164235A1 (en) * 2003-01-21 2004-08-26 Miller Steven D. System and method for significant dust detection and enhancement of dust images over land and ocean
CN102636779A (zh) * 2012-05-07 2012-08-15 武汉大学 一种基于重采样回归分析的亚像素积雪覆盖率提取方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
柯长青 等: "北极海冰的光谱特征分析", 《光谱学与光谱分析》 *
裴欢: "基于MODIS数据的北疆积雪信息提取及其应用研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105116464A (zh) * 2015-08-12 2015-12-02 南京大学 基于神经网络模型的极地海冰融池提取方法
CN106290183A (zh) * 2016-07-19 2017-01-04 国网天津市电力公司 基于光谱反射比的雪状态参数测量方法
CN106290183B (zh) * 2016-07-19 2018-12-04 国网天津市电力公司 基于光谱反射比的雪状态参数测量方法
CN110866363A (zh) * 2019-11-12 2020-03-06 南京信息工程大学 一种利用人工神经网络反演北极融池分布的方法
CN110866363B (zh) * 2019-11-12 2023-04-07 南京信息工程大学 一种利用人工神经网络反演北极融池分布的方法
CN116597319A (zh) * 2023-04-19 2023-08-15 北京师范大学 一种海冰日数据综合处理方法、装置及计算设备
CN116597319B (zh) * 2023-04-19 2024-02-02 北京师范大学 一种海冰日数据综合处理方法、装置及计算设备

Similar Documents

Publication Publication Date Title
Tachikawa et al. Characteristics of ASTER GDEM version 2
Santos et al. Differences in coastal and oceanic SST warming rates along the Canary upwelling ecosystem from 1982 to 2010
Ji et al. Estimating aboveground biomass in interior Alaska with Landsat data and field measurements
CN104198052B (zh) 基于海洋二号卫星扫描微波辐射计的海冰密集度获取方法
CN104535979A (zh) 一种陆地云光学厚度的遥感反演方法及***
CN104007443A (zh) 一种利用光谱特征提取极地海冰融池的方法
CN104615848A (zh) 一种集森林火险预报和地面火点监测于一体的综合应用方法
CN103927455A (zh) 基于高分一号卫星的陆地气溶胶光学性质反演方法
Dewidar Changes in the shoreline position caused by natural processes for coastline of Marsa Alam–Hamata, Red Sea, Egypt
CN105116464A (zh) 基于神经网络模型的极地海冰融池提取方法
Manyilizu et al. Interannual variability of sea surface temperature and circulation in the tropical western Indian Ocean
Liu et al. Characterization of recent glacier decline in the Cordillera Real by LANDSAT, ALOS, and ASTER data
CN103983360B (zh) 一种基于hj-1b irs卫星数据反演lst的方法
Hereher Environmental monitoring and change assessment of Toshka lakes in southern Egypt using remote sensing
Guo et al. A comparative study of coastline changes at Tampa Bay and Xiangshan Harbor during the last 30 years
Ma et al. Detection of thermal pollution from power plants on China’s eastern coast using remote sensing data
Rao et al. Near surface atmospheric response to the total solar eclipse at Dibrugarh on 22 July 2009
Kondrik et al. Tendencies in coccolithophorid blooms in some marine environments of the Northern Hemisphere according to the data of satellite observations in 1998–2013
Jun et al. Fog detection over China’s adjacent sea area by using the MTSAT geostationary satellite data
Le Hung et al. Evaluation of soil erosion risk using remote sensing and GIS data (a case study: Lang Chanh district, Thanh Hoa province, Vietnam)
Wen et al. Research on the features of chlorophyll-a derived from rapideye and EOS/MODIS data in chaohu lake
Khusrul Amin et al. Analysis of the wetland degradation around the vicinity of Dhaka city in Bangladesh
Yang et al. FY-3A satellite microwave data assimilation experiments in tropical cyclone forecast
Wiehle et al. Geophysical and morphological processes in the Arctic coastal regions and marginal ice zone by using remote sensing data
Pimentel et al. Snow evolution in a semi-arid mountainous area combining snow modelling and Landsat spectral mixture analysis

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
WD01 Invention patent application deemed withdrawn after publication
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20140827