CN103487804B - 基于冰雷达技术的南极冰下湖半自动识别方法 - Google Patents
基于冰雷达技术的南极冰下湖半自动识别方法 Download PDFInfo
- Publication number
- CN103487804B CN103487804B CN201310481963.3A CN201310481963A CN103487804B CN 103487804 B CN103487804 B CN 103487804B CN 201310481963 A CN201310481963 A CN 201310481963A CN 103487804 B CN103487804 B CN 103487804B
- Authority
- CN
- China
- Prior art keywords
- interface
- ice
- lake
- radar
- semi
- 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.)
- Expired - Fee Related
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
- G01S7/411—Identification of targets based on measurements of radar reflectivity
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明涉及基于冰雷达技术的南极冰下湖半自动识别方法,考虑冰下湖与冰的界面在冰雷达断面图像上的特征具有回波信号强度大、界面平坦的特点,利用表面粗糙度与雷达反射功率相结合的自动提取方法,实现了冰下湖面的计算机初步识别,大大提高了识别效率和识别精度;再结合人工修改,对冰下湖的位置进行修正,既不会增加太大工作量,又可以弥补自动提取的误差,使得最终的识别结果更加准确,速度和可靠性也更高。
Description
技术领域
本发明涉及一种基于冰雷达技术的南极冰下湖半自动识别方法,属于地球物理与遥感地学应用领域。
背景技术
冰下湖是人类的未知领域,也是地球上最后的重要科学前沿领域之一。其不但具有研究全球和区域气候变化机制的大量连续的古气候沉积记录,而且其水系循环对于冰动力过程也有重要影响,是未来海平面变化不确定性的重要影响因素之一,因此对其进行探测和识别研究具有重要的科学意义。
目前冰下湖的识别主要有两种方法。一种是利用卫星测高发现因冰下湖排水引起的冰盖表面起伏变化进行识别。另一种是利用冰雷达测定的冰下底部回波进行判断。相对于卫星测高方法,冰雷达发射的电磁波能够穿透几千米厚的冰层,获得冰底反射的回波信号,再根据信号特征识别冰下湖。这种观测方法更直观,而且能排除其他多种因素的干扰,所得结果准确性更高。相对于冰岩界面,冰下湖与冰的界面在冰雷达断面图像上具有回波信号强度大、界面平坦等特点。过去依据这些特征,多采用人工专业判读的方法从断面图像上识别冰下湖,具有工作量大且耗时长等缺点。
早期冰雷达回波信号记录在胶片上,是模拟信号,通过专业人员对胶片进行判读,根据上述特征找出可能有冰下湖的胶片,然后冲洗在相纸上。资深专家再进一步判读是否存在冰下湖。如果识别为冰下湖,再结合飞行导航的位置记录,即可推断冰下湖的位置。这种方式工作量大且耗时长。随着计算机和信号处理技术的发展,冰雷达数据已实现数字存储,冰下湖的半自动化识别已成为必然发展趋势。
发明内容
本发明要解决的技术问题是:克服现有技术缺点,提出一种基于冰雷达技术的南极冰下湖半自动识别方法,能够自动快速地从冰雷达数据中提取冰水和冰岩界面,通过人工删除界面中的孤立点,连接断点,完成冰下湖的识别。
为了解决上述技术问题,本发明提出的技术方案是:一种基于冰雷达技术的南极冰下湖半自动识别方法,包括以下步骤:
第一步、将读取的冰雷达数据进行预处理,获取可视化的冰雷达图像;
第二步、在冰雷达图像上,自下而上地寻找各道首次出现雷达反射功率突然增大的区域,该区域的界面包含有冰岩界面或冰水界面,结合GPS数据获得界面的坐标和高程;
第三步、根据所述区域界面的高程信息获取界面的粗糙度;
第四步、提取粗糙度小于0.1且雷达反射功率大于设定阈值的界面,获得冰下湖初步识别结果;
第五步、结合图像,人工删除界面中的孤立点,连接冰下湖的断点,完成冰下湖的识别。
本发明基于冰雷达技术的南极冰下湖半自动识别方法,还具有如下改进:
1、所述第一步中,冰雷达数据的预处理主要包括带通滤波降噪、手动或自动剔除坏道。
2、第二步中,自下而上首次出现雷达反射功率突然增大区域的界面提取方法如下:
步骤a、在垂直方向上获得反射功率突然增大区域的最大点,反射功率最大点构成准界面;
步骤b、对准界面进行人工修正,以确保得到的整个界面连续为准则,当准界面位置附近存在多个明显的反射层,且其中至少有一个反射层与两侧准界面连续时,以最上方与两侧准界面连续的反射层为准,最终完成界面半自动提取。
3、第三步中,对所述区域界面高程进行傅里叶变换计算粗糙度。
4、第四步中,设定的阈值范围为[0.75*a,0.85*a],a为整个区域界面中反射功率的最大值。
本发明采用半自动提取法获取界面,先获得反射功率突然增大区域内的最大点,然后根据规则对界面进行修正,保证了提取界面方法的统一性,使得提取界面的结果更客观,更准确。相比传统的人工数字化方法简单易行,效率和精度都比较高。
本发明利用表面粗糙度与雷达反射功率相结合的自动提取方法,实现了冰下湖面的计算机初步识别,大大提高了识别效率和识别精度;再结合人工修改,对冰下湖的位置进行修正,既不会增加太大工作量,又可以弥补自动提取的误差,使得最终的结果更加准确。
由于第二步中已经获得了区域界面,但是不能确定是冰水界面还是冰岩界面,考虑到冰水界面的雷达反射功率大于冰岩界面的雷达反射功率,经反复实验表明,当设定阈值位于[0.75*a,0.85*a](a为整个区域界面中反射功率的最大值)时,能够剔除冰岩界面,同时加上界面粗糙度作为另一个筛选条件,可以确保冰下湖面的提取精度。
本发明采用对冰岩或者冰水界面进行傅里叶变换计算粗糙度,这是国际上计算粗糙度的常用方法,成熟度高且运算简单,无需大量的编程工作。大量实践表明该方法获得的粗糙度能真实反映界面的起伏状况,且在因采样间隔大获得的点数稀少时仍能满足要求。
综上,本发明方法的执行步骤简单易行,由于界面提取和湖面识别都采用半自动化实现,因此冰下湖提取精度和效率都较高,使得本发明方法适合于从海量冰雷达数据中对冰下湖进行识别。
附图说明
下面结合附图对本发明作进一步的说明。
图1是本发明基于冰雷达技术的南极冰下湖半自动识别方法流程图。
图2是一幅包含冰下湖界面的冰雷达图像。
图3是提取的冰岩或者冰水界面的位置。
图4是得到的界面反射功率值。
图5是计算得到的界面粗糙度。
图6是粗糙度小于0.1的位置在界面上的表现。
图7是反射功率大于阈值的位置在界面上的表现。
图8是粗糙度小且反射功率大的位置在界面上的表现。
图9是经人工修整后最终识别出的冰下湖。
具体实施方式
下面根据附图详细阐述本发明,使本发明的目的和效果变得更加明显。
包含冰下湖的冰雷达数据目前在国际上公开较少,因此本发明的实例采用的数据主要为R.Jacobel等于2010年发布于美国国家冰雪数据中心(USA:National Snow and Ice Data Center)上名称为《Radar Studies of InternalStratigraphy and Bed Topography along the US ITASE-II Traverse》的冰雷达二级产品(编号为nsidc-0475)。该数据已经完成了本发明所述步骤1和2的操作。但为了更完整且更清晰地描述本发明的实施方式,下文仍将对这两步做详细说明。
图1所示,为本发明实施例基于冰雷达技术的南极冰下湖半自动识别方法的流程图,具体包括以下步骤:
第一步、将读取的冰雷达数据进行预处理,获取可视化的冰雷达图像;预处理主要包括带通滤波降噪、手动或自动剔除坏道。
本步骤为冰雷达数据读取和预处理的常规操作。在读取数据之后计算反射信号功率值,并进行去噪处理。由于数据中存在雷达信号异常引起的坏道和停车引起的停车道等异常情况,需要对其进行自动和手动删除。自动删除过程中计算每一道所有反射功率的标准偏差,如果小于某一阈值(由具体的雷达数据确定,小于该阈值说明没有明显的反射层,也即坏道),则删除。停车道之间连续且跨度较宽,采用手动删除的方式。最终获取的冰雷达图像如附图2所示。
第二步、在冰雷达图像上,自下而上地寻找各道首次出现雷达反射功率突然增大的区域,包含有冰岩界面或冰水界面,结合GPS数据获得界面的坐标和高程。
本步骤中,自下而上首次出现雷达反射功率突然增大区域的界面提取方法如下:
步骤a、在垂直方向上获得反射功率突然增大区域的最大点,反射功率最大点构成准界面;
步骤b、对所述准界面进行人工修正,以确保得到的整个界面连续为准则,当准界面位置附近存在多个明显的反射层,且其中至少有一个反射层与两侧准界面连续时,以最上方与两侧准界面连续的反射层为准,最终完成界面半自动提取。
上述人工修正之后,获取该连续界面距离冰盖表面的深度值(界面上覆冰层厚度),同时保存该界面位置的反射功率值。
本步骤中,结合GPS数据获得该区域界面的坐标和高程是这样实现的:将与冰雷达获取时间相对应的GPS数据的经纬度投影到极地坐标系中,所得投影点按冰雷达采样时间间隔(每道的时间)进行插值,同时对GPS测得的表面高程进行插值。将冰雷达每道与GPS数据对应起来,并对高程与冰厚求差得界面高程值。图3是获取的界面深度随距离的变化图,图4是该界面的反射功率值。
第三步、对所述区域界面的高程进行傅里叶变换计算粗糙度。
本步骤中,首先需要判断采样点的间距。如果间距不等,需重新采样以尽量保持间距均匀(设该均匀间距为d),对间隔较大的数据点还需要插值。根据Taylor等于2004年发表在Computers&Geosciences上第30卷第8期《Regional-scale bed roughness beneath ice masses:measurement andanalysis》一文中提出的有效插值不超过5个连续点的建议,在间隔小于等于均匀间距d的6倍的相邻两点之间进行插值。然后采用32个点为滑动窗口,计算窗口内的粗糙度。粗糙度的表示采用李鑫等于2010年发表在Journal ofGlaciology上第56卷第199期《Characterization of subglacial landscapesby a two-parameter roughness index》一文中提出的计算界面起伏程度系数ξ。界面粗糙度的具体计算方法属于现有技术范畴,本实施例不再详细赘述。
经本步骤,获得的界面粗糙度如图5所示。
第四步、从区域界面中提取粗糙度小于0.1且反射功率大于设定阈值的界面,获得冰下湖初步识别结果。
因为冰下湖表面极为平坦,粗糙度几乎可认为是0。所以可将界面上的粗糙度进行对比,识别粗糙度最小值所在位置。本实例根据对比选择小于0.1的值为粗糙度最小的情况,将符合该情况的所有位置标于图中,如图6所示。图中加粗的点、线表示界面中粗糙度小于0.1的部分。
相对于基岩表面的反射功率,冰下湖表面的反射功率更大。因此,可以通过反射功率值进一步识别冰下湖。本实施例中,设定阈值取80%的界面最大反射功率。第二步中提取的界面的最大反射功率为37500dB,设定阈值=30000dB。将反射功率值大于设定阈值(30000dB)的情况作为冰下湖的界面。本实例根据对比,将这些位置标于图中,如图7所示。图中加粗的点、线代表界面中反射功率大于设定阈值的部分。
本步骤中,将粗糙度与反射功率相结合(求交集)的提取结果,如图8所示,加粗的点、线代表冰下湖初步识别结果。
第五步、结合图像,人工删除界面(图8)中的孤立点,连接冰下湖的断点,完成冰下湖的识别。
本步骤中,将最终的结果标于图上,由于自动识别存在误判的情况,会产生一些错误的孤立点(孤立或者空间分布距离短),同时还会造成冰下湖界面的不连续。需要人工删除此类孤立点,并将空间分布较长且相距间隔较小的点连接,得到最终识别的结果,如图9所示。图9中,加粗的线代表冰下湖,细线代表基岩表面。
除上述实施例外,本发明还可以有其他实施方式。凡采用等同替换或等效变换形成的技术方案,均落在本发明要求的保护范围。
Claims (4)
1. 基于冰雷达技术的南极冰下湖半自动识别方法,包括以下步骤:
第一步、将读取的冰雷达数据进行预处理,获取可视化的冰雷达图像;
第二步、在冰雷达图像上,自下而上地寻找各道首次出现雷达反射功率突然增大的区域,该区域的界面包含有冰岩界面或冰水界面,结合GPS数据获得界面的坐标和高程;
第三步、根据所述区域界面的高程信息获取界面的粗糙度;
第四步、提取粗糙度小于0.1且雷达反射功率大于设定阈值的界面,获得冰下湖初步识别结果;设定阈值的范围为[0.75*a,0.85*a],a为整个区域界面中反射功率的最大值;
第五步、结合图像,人工删除界面中的孤立点,连接冰下湖的断点,完成冰下湖的识别。
2. 根据权利要求1所述的基于冰雷达技术的南极冰下湖半自动识别方法,其特征在于:第一步中,冰雷达数据的预处理主要包括带通滤波降噪、手动或自动剔除坏道。
3. 根据权利要求1所述的基于冰雷达技术的南极冰下湖半自动识别方法,其特征在于:第二步中,自下而上首次出现雷达反射功率突然增大区域的界面提取方法如下:
步骤a、在垂直方向上获得反射功率突然增大区域的最大点,反射功率最大点构成准界面;
步骤b、对准界面进行人工修正,以确保得到的整个界面连续为准则,当准界面位置附近存在多个明显的反射层,且其中至少有一个反射层与两侧准界面连续时,以最上方与两侧准界面连续的反射层为准,最终完成界面半自动提取。
4. 根据权利要求1所述的基于冰雷达技术的南极冰下湖半自动识别方法,其特征在于: 第三步中,对所述区域界面高程进行傅里叶变换计算粗糙度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310481963.3A CN103487804B (zh) | 2013-10-15 | 2013-10-15 | 基于冰雷达技术的南极冰下湖半自动识别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310481963.3A CN103487804B (zh) | 2013-10-15 | 2013-10-15 | 基于冰雷达技术的南极冰下湖半自动识别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103487804A CN103487804A (zh) | 2014-01-01 |
CN103487804B true CN103487804B (zh) | 2015-10-28 |
Family
ID=49828173
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310481963.3A Expired - Fee Related CN103487804B (zh) | 2013-10-15 | 2013-10-15 | 基于冰雷达技术的南极冰下湖半自动识别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103487804B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110441768A (zh) * | 2019-09-09 | 2019-11-12 | 北京工业大学 | 一种快速提取南极航空冰雷达数据中冰面与基岩界面的方法 |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103809179B (zh) * | 2014-02-27 | 2016-03-02 | 中国科学院电子学研究所 | 机载探冰雷达成像方法 |
CN109061647B (zh) * | 2018-07-04 | 2020-11-10 | 烟台雷奥电子科技有限公司 | 一种基于雷达的间断海冰区域检测方法、装置及设备 |
CN109283526B (zh) * | 2018-10-26 | 2020-04-24 | 浙江大学 | 一种基于全极化雷达的冰盖内部组构和冰流场分布特征识别方法 |
CN113671490B (zh) * | 2021-08-12 | 2023-06-02 | 同济大学 | 基于基底信号形态特征的南极冰下水探测方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101963509A (zh) * | 2009-05-08 | 2011-02-02 | 通用汽车环球科技运作公司 | 通过地形变化分析来检测畅通路径的方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2577532A4 (en) * | 2010-05-28 | 2016-04-06 | Geostellar Inc | SYSTEM AND METHOD FOR BASIC GEOMATIC MODELING OF VARIOUS RESOURCES ON LARGE LANDSCAPES |
-
2013
- 2013-10-15 CN CN201310481963.3A patent/CN103487804B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101963509A (zh) * | 2009-05-08 | 2011-02-02 | 通用汽车环球科技运作公司 | 通过地形变化分析来检测畅通路径的方法 |
Non-Patent Citations (1)
Title |
---|
机载SAR图像的海冰检测方法的研究;张珣;《中国优秀硕士学位论文全文数据库》;20121015;正文第5,6,9,18-21,46页 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110441768A (zh) * | 2019-09-09 | 2019-11-12 | 北京工业大学 | 一种快速提取南极航空冰雷达数据中冰面与基岩界面的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103487804A (zh) | 2014-01-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103487804B (zh) | 基于冰雷达技术的南极冰下湖半自动识别方法 | |
Jay et al. | Walrus areas of use in the Chukchi Sea during sparse sea ice cover | |
Baltsavias et al. | Digital surface modelling by airborne laser scanning and digital photogrammetry for glacier monitoring | |
Liu et al. | Wavelet‐transform based edge detection approach to derivation of snowmelt onset, end and duration from satellite passive microwave measurements | |
CN109804232A (zh) | 一种基于红外热像图分析的沥青路面裂缝发育程度检测方法 | |
CN102542551B (zh) | 极地冰盖边缘浮冰自动变化检测技术 | |
CN110568416B (zh) | 一种基于遥感影像的雷达有效检测区域提取方法 | |
CN108197583A (zh) | 基于图割优化和影像结构特征的建筑物变化检测方法 | |
CN104050474A (zh) | 一种基于LiDAR数据的海岛岸线自动提取方法 | |
Middleton et al. | Resolution and accuracy of an airborne scanning laser system for beach surveys | |
CN108537116B (zh) | 一种基于多尺度特征的海岸线二级类型提取方法及系统 | |
Drews et al. | Characterizing the glaciological conditions at Halvfarryggen ice dome, Dronning Maud Land, Antarctica | |
CN116645789A (zh) | 一种基于无人机测绘的土地自然灾害预警管控系统 | |
CN111639865A (zh) | 一种高铁沿线气象灾害发生风险分析方法 | |
CN102567726B (zh) | 极地冰盖边缘区域浮冰自动提取技术 | |
CN114280572A (zh) | 去除信号干扰杂波的单雷达回波质控方法、系统及终端 | |
Hendricks et al. | Sea-ice thickness variability in Storfjorden, Svalbard | |
Bollen et al. | Dynamic mass loss from Greenland's marine-terminating peripheral glaciers (1985–2018) | |
Zhang et al. | An automatic algorithm to extract nearshore bathymetric photons using pre-pruning quadtree isolation for ICESat-2 data | |
CN111881566B (zh) | 基于实景仿真的滑坡位移检测方法及装置 | |
Gallop et al. | Video observations of rip currents on an embayed beach | |
Malsam | Characterization of rockfall activity and identification of weather-rockfall relationships using high temporal resolution remote sensing methods | |
Li et al. | Snow stratigraphy observations from Operation IceBridge surveys in Alaska using S and C band airborne ultra-wideband FMCW (frequency-modulated continuous wave) radar | |
CN104297753A (zh) | 一种基于自适应缩减算子的导航雷达图像反演海面风向方法 | |
Tarmizi et al. | Shoreline data extraction from QuickBird satellite image using semi-automatic technique |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20151028 Termination date: 20161015 |
|
CF01 | Termination of patent right due to non-payment of annual fee |