一种基于卷积神经网络的海冰类型遥感分类方法技术

本发明专利技术公开了一种基于卷积神经网络的海冰类型遥感分类方法。属于遥感地学应用技术领域。本方法采用卫星雷达高度计数据对海冰进行分类（分为多年冰、一年冰和开阔水域三种类型），其步骤包括从高度计数据中获取测量点的经纬度坐标及其雷达回波波形。下载海冰类型数据，并提取相应测量点位处的海冰类型信息。将上述提取的数据信息进行空间位置匹配，通转换为带经纬度坐标的矢量点数据。所获取的雷达回波波形和对应的海冰类型作为训练数据，对卷积神经网络进行训练，对待分类的雷达回波波形进行海冰类型的识别。分类后将海冰类型标记到矢量数据上；标记后的矢量数据在通过投影变换、栅格转换和空间重采样等处理，获得整个研究区的海冰类型。

【技术实现步骤摘要】
一种基于卷积神经网络的海冰类型遥感分类方法
本专利技术涉及一种基于卷积神经网络的海冰类型遥感分类方法，属于遥感应用
技术背景海冰类型的变化是气候变化的重要指示器，其变化会影响海冰厚度的反演、船只航行路线的划定以及其他人类重要的极地活动。早期海冰类型数据只能通过人工实地调查得到，之后利用光学和SAR遥感影像可以实现海冰的分类，但是由于遥感影像空间覆盖范围小，光学影像受光线、云的影响较大，SAR影像获取成本比较高等困难的限制，难以实现大尺度的海冰分类。随着卫星雷达高度计的发展，可利用高度计数据对海冰进行分类。虽然高度计常常被用来反演海冰厚度，但自1990年以来，多位学者相继发现不同类型的海冰其雷达回波波形具有一定的差异性，并且这种差异可以用来对海冰进行分类。同时由于高度计的空间覆盖范围更大，时间分辨率也较高，因此可采用高度计数据实现的海冰分类。然而要实现更为准确的海冰分类，需要更加先进的雷达高度计。2010年4月，欧空局发射的CryoSat-2卫星，是目前为止最为先进的卫星雷达高度计，观测数据范围可达南北纬88°N的区域，且拥有369天的全周期和30天的小周期，一个月即可实现对南极海冰的全部观测。卫星上携带的合成孔径干涉雷达高度计(SIRAL)的垂直测量精度达到1～3cm，并采用延迟多普勒雷达高度计(DDA)技术将卫星地面足迹减小到沿轨约为0.3km，跨轨约为1.5km，分辨率优于传统的高度计。此外，SIRAL对地表点进行多视处理以减少雷达斑点引起的噪声，其测量海平面高度数据精度约为传统雷达高度计的2倍。美国国家海冰中心(NIC)每周提供一次南极海域的冰况图。该海冰类型产品是使用光学、近红外、SAR卫星数据和船舶走航数据得到的。采用其shapefile格式的数据产品，产品中主要包含三种海冰类型数据：一年冰、多年冰和开阔水域。因此，我们以美国国家海冰中心的冰况图作为参考对照数据，针对全南极进行一年冰、多年冰和开阔水域的分类。
技术实现思路
本专利技术要解决的技术问题是：提供一种基于CryoSat-2卫星测高数据的海冰分类方法。考虑到高度计数据较高的空间分辨率和时间分辨率，该数据大大提高了海冰分类的时效性，同时提出了一种快速、大尺度的海冰分类方法。为了解决上述技术问题，本专利技术提出的技术方案是：一种基于卷积神经网络的海冰类型遥感分类方法，包括以下步骤：第一步：准备训练数据和待分类数据,具体包括如下几个方面的内容：a.下载某一时期内和待分类的CryoSat-2卫星SAR模式L1b级数据，将所述某一时期内的数据作为训练数据；读取L1b级数据DBL格式的原始文件，从中提取各测量点的经纬度和雷达回波波形；b.下载与训练数据相同时期内的美国国家海冰中心的冰况图数据使用多年冰、一年冰和开阔水域数据，获得相应的海冰类型信息；第二步、获取雷达回波波形图像和提取波形特征：根据训练数据和待分类数据中所有测量点的雷达回波波形生成雷达回波波形图像像，并提取波形特征值用于后续的海冰分类；所述波形特征值包含脉冲宽度、前缘宽度和栈标准差；第三步、分别将训练数据和待分类数据中提取的波形特征值与对应的经纬度坐标进行空间匹配；将匹配好的数据转换成带经纬度坐标的矢量点数据，同时记录每个雷达回波波形图像所代表的测量点处的经纬度；第四步、分类前处理：对于训练数据和待分类数据，剔除前缘宽度大于14的雷达回波波形图像，再剔除脉冲宽度小于0.3且栈标准差大于4的雷达回波波形图像；第五步、生成训练样本：将训练数据中带有经纬度坐标和波形特征值的矢量点与同时期内的冰况图数据进行空间匹配，进而得到各矢量点的海冰类型信息，并根据矢量点的经纬度信息得到各自雷达回波波形图像所代表的海冰类型，将该雷达回波波形图像集作为训练样本。第六步、训练分类器并利用分类器进行分类：并利用训练样本对卷积神经网络进行训练，使用训练完毕后的卷积神经网络对待分类的雷达回波波形图像进行分类，分类完成后根据经纬度信息对相应的矢量点数据进行类型标记；第七步、生成海冰类型图：将矢量点数据投影至栅格图内，每个栅格内所有矢量点的海冰类型的众数作为该栅格的海冰类型，并根据栅格的海冰类型对栅格赋值。本专利技术利用了CryoSat-2数据和NIC海冰类型数据，数据获取相对简单，操作便捷。CryoSat-2可以提供精确的波形信息，NIC海冰类型数据可以提供时空分辨率极高的海冰类型数据。利用卷积神经网络将两者有效地结合，从而可以用于大尺度的海冰分类，为后续海冰厚度的估算和气候变化的研究提供了较好的理论基础。本专利技术相关数据提取和矢量转换过程均通过MATLAB编程和ArcGIS平台实现，减少了人工参与，提高了分类效率。综上，本专利技术方法的执行步骤简单易行，分类效果较好。目前大尺度海冰类型资料很少，基于传统的实地调查和光学、SAR影像的分类方法获取的海冰类型范围十分有限。本专利技术使用高精度、大空间覆盖范围的卫星雷达高度计CryoSat-2数据，利用一种经典的深度学习方法-卷积神经网络探究雷达回波波形与海冰类型之间的关系，实现大空间尺度上的海冰分类。对精确量化南极海冰物质平衡变化、研究南极海冰类型对全球气候的影响有深远的科学意义。附图说明下面结合附图对本专利技术作进一步的说明。图1是基于卷积神经网络的海冰分类方法流程图。图2为2017年9月NIC冰况图。图3为提取的矢量点的经纬度坐标和3种波形特征数据。图4与NIC冰况图匹配后包含海冰类型的矢量点的有关信息。图5根据经纬度信息得到雷达回波波形图像的海冰类型。图6为2017年9月海冰的分类结果。具体实施方式下面根据附图详细阐述本专利技术，使本专利技术的技术路线和操作步骤更加清晰。本专利技术实例采用的数据是CryoSat-2卫星L1b级SAR模式基线C数据和NICshapfile格式数据。示例数据获取时间为2017年9月1日——2017年9月30日，NIC海冰类型数据获取时间为2017年9月14日。图1为基于卷积神经网络的海冰遥感分类方法流程图，具体步骤如下:第一步：准备训练数据和待分类数据，并对数据进行信息读取。具体包括如下几个方面的内容：a.下载某一时期内的CryoSat-2卫星SAR模式L1b级数据，获得示例数据(2017年9月1日——2017年9月30日)，按照训练样本和测试样本3∶1的比例随机生成训练数据和待分类数据，从两个数据中提取每个测量点的经纬度信息和雷达回波波形。b.读取同时期NIC冰况图(海冰类型)shapefile格式文件(2017年9月)，获得研究区海冰类型分布图。2017年9月NIC冰况图如图2所示。第二步、获取雷达回波波形图像和提取波形特征值。根据训练数据和待分类数据中所有测量点的雷达回波波形生成雷达回波波形图像，并提取波形特征值用于后续的海冰分类。波形特征值主要包含脉冲宽度、前缘宽度、栈标准差，其中栈标准差可以直接从Cryosat-2L1b数据中获取，其他特征计算式如下：脉冲宽度PP是最大能量与累计能量和的比值：式中，Pi是在第i个距离门处的回波能量值，Pmax是回波波形中最大回波能量值；前缘宽度是能量值等于第一个1％Pmax的至能量值等于第一个最大能量值99％Pmax之间的距离门数量。第三步、分别将训练数据和待分类数据中提取的波形特征值与对应的经纬度坐标进行本文档来自技高网...

【技术保护点】
1.一种基于卷积神经网络的海冰类型遥感分类方法，包括以下步骤：第一步：准备训练数据和待分类数据,具体包括如下几个方面的内容：a.下载某一时期内和待分类的CryoSat‑2卫星SAR模式L1b级数据，将所述某一时期内的数据作为训练数据；读取L1b级数据DBL格式的原始文件，从中提取各测量点的经纬度和雷达回波波形；b.下载与训练数据相同时期内的美国国家海冰中心的冰况图数据使用多年冰、一年冰和开阔水域数据，获得相应的海冰类型信息；第二步、获取雷达回波波形图像和提取波形特征值：根据训练数据和待分类数据中所有测量点的雷达回波波形生成雷达回波波形图像，并提取波形特征值用于后续的海冰分类；所述波形特征值包含脉冲宽度、前缘宽度和栈标准差；第三步、分别将训练数据和待分类数据中提取的波形特征值与对应的经纬度坐标进行空间匹配；将匹配好的数据转换成带经纬度坐标的矢量点数据，同时记录每个雷达回波波形图像所代表的测量点处的经纬度；第四步、分类前处理：对于训练数据和待分类数据，剔除前缘宽度大于14的雷达回波波形图像，再剔除脉冲宽度小于0.3且栈标准差大于4的雷达回波波形图像；第五步、生成训练样本：将训练数据中带有经纬度坐标和波形特征值的矢量点与同时期内的冰况图数据进行空间匹配，进而得到各矢量点的海冰类型信息，并根据矢量点的经纬度信息得到各自雷达回波波形图像所代表的海冰类型，将该雷达回波波形图像集作为训练样本。第六步、训练分类器并利用分类器进行分类：并利用训练样本对卷积神经网络进行训练，使用训练完毕后的卷积神经网络对待分类的雷达回波波形图像进行分类，分类完成后根据经纬度信息对相应的矢量点数据进行类型标记；第七步、生成海冰类型图：将矢量点数据投影至栅格图内，每个栅格内所有矢量点的海冰类型的众数作为该栅格的海冰类型，并根据栅格的海冰类型对栅格赋值。...

【技术特征摘要】
1.一种基于卷积神经网络的海冰类型遥感分类方法，包括以下步骤：第一步：准备训练数据和待分类数据,具体包括如下几个方面的内容：a.下载某一时期内和待分类的CryoSat-2卫星SAR模式L1b级数据，将所述某一时期内的数据作为训练数据；读取L1b级数据DBL格式的原始文件，从中提取各测量点的经纬度和雷达回波波形；b.下载与训练数据相同时期内的美国国家海冰中心的冰况图数据使用多年冰、一年冰和开阔水域数据，获得相应的海冰类型信息；第二步、获取雷达回波波形图像和提取波形特征值：根据训练数据和待分类数据中所有测量点的雷达回波波形生成雷达回波波形图像，并提取波形特征值用于后续的海冰分类；所述波形特征值包含脉冲宽度、前缘宽度和栈标准差；第三步、分别将训练数据和待分类数据中提取的波形特征值与对应的经纬度坐标进行空间匹配；将匹配好的数据转换成带经纬度坐标的矢量点数据，同时记录每个雷达回波波形图像所代表的测量点处的经纬度；第四步、分类前处理：对于训练数据和待分类数据，剔除前缘宽度大于14的雷达回波波形图像，再剔除脉冲宽度小于0.3且栈标准差大于4的雷达回波波形图像；第五步、生成训练样本：将训练数据中带有经纬度坐标和波形特征值的矢量点与同时期内的冰况图数据进行空间匹配，进而得到各矢量点的海冰类型信息，并根据矢量点的经纬度信息得到各自雷达回波波形图像所代表的海冰类型，将该雷达回波波形图像集作为训练样本。第六步、训练分类器并利用分类器进行分类：并利用训练样本对卷积神经网络进行训练，使用训练完毕后的卷积神经网络对待分类的雷达回波波形图像进行分类，分类完成后根据经纬度信息对相应的矢量点数...

【专利技术属性】
技术研发人员：柯长青，沈校熠，李萌萌，蔡宇，李海丽，
申请(专利权)人：南京大学，
类型：发明
国别省市：江苏,32