本发明专利技术涉及基于神经网络模型的极地海冰融池提取方法，属于遥感信息全自动提取方法领域。考虑MODIS表面反射率数据包含了不同类型地表光谱组合的特点，根据极地地区冰雪、水体、融池三者光谱响应曲线的差异，结合混合像元分解原理建立各波段反射率与三者所占比例之间的多元线性方程组，并增加对方程组解的约束条件。利用神经网络方法优化并快速解算此线性方程组，得到各像元内冰雪、水体、融池所占比例。融池区域是相对于除水体外的海冰表面而言的，因此进一步提取出相对的融池覆盖比例，实现对极地地区海冰融池范围的全自动提取。
【技术实现步骤摘要】
本专利技术涉及一种，属于遥感信息全自动提取方法领域。
技术介绍
海冰融化过程中，融水在海冰表面低凹处汇集形成融池。由于融池的反照率介于海冰和海水之间，提取融池覆盖率对于准确计算太阳辐射在极区海洋的吸收、了解极区气-冰-海耦合系统的热力过程有重要意义。融池的重要性突出表现在它的物理特性上。融池就是冰上水池，所以融池具有海水的一切性质，尤其在于其反照率，反照率是决定海冰吸收太阳辐射能量多少的一个标准，反照率高，则吸收的太阳辐射能量就少，融化就慢，反之融化就快。而海水的表面反照率为10%——15%，海冰的反照率为90%，差别很大。如果某一块海冰上有一个小融池，该融池有可能是海水的动力学或者是热力学作用形成的，那么在该块海冰上由于反照率的大幅减少，会吸收更多的太阳辐射能量，尤其在夏天又很少出现冻结情况，就会加速该块海冰的融化。夏季极地遥感影像中往往含有大量的融池区域，相关的数据也表明，北极夏季海冰的反照率有一定程度的降低，印证了融池使海冰融化加剧的现象。获取融池的相关数据，对更好地分析海冰对全球气候变暖的响应有重大的科学意义。覆盖率是融池模拟中的一项重要参数，也是空间遥感在融池研究中的主要研究方向，然而目前有关的遥感算法仍以局地性的个例实验为主，没有适用于更大范围的融池遥感提取算法。在北极不同地区已经进行过众多的海冰表面观测和船基走航观测，以研究融池的反射率、光谱特征以及融池的大小和分布等。虽然有较高的观测准确率和精度，但是局限于航线，观测结果的时间和空间尺度都极其有限，仅能反映航线涉及区域的融池情况。需要一种在较长时间序列和空间范围内提取极地海冰融池的遥感方法。
技术实现思路
本专利技术要解决的技术问题是:克服现有技术上述缺陷，提供一种基于神经网络模型的海冰融池提取方法，在极地地区遥感图像中确定冰雪、水体和融池的覆盖比例，通过快速优化地解算基于混合像元分解模型的多元线性方程组，减少计算代价并且保证精度，完成极地海冰融池范围的提取。为了解决上述技术问题，本专利技术提出的技术方案是:一种，其步骤包括:第一步、将MODIS表面反射率产品M0D09A1由正弦投影转换为极地方位投影；第二步、结合M0D09A1产品中的陆地掩膜和云掩膜数据，提取目标区域的图像，并选取λ” λ2、λ 3波段的表面反射率值；第三步、根据冰雪、水体和融池的光谱响应曲线，得到三者分别在λ P λ2、λ3波段的标准光谱反射率值:在A1波段:冰雪的光谱反射率F1(A1) = 0.95，水体的光谱反射率L(A1)=0.08，融池的光谱反射率rj λ J = 0.16 ；在λ2波段:冰雪的光谱反射率Γι(λ2) = 0.87，水体的光谱反射率rwU2)=0.08，融池的光谱反射率Γη(λ2) = 0.07 ；在λ3波段:冰雪的光谱反射率Γι(λ3) = 0.95，水体的光谱反射率rwU3)=0.08，融池的光谱反射率4 ( λ 3) =0.22；第四步、对图像中的每个像元，根据其中冰雪、水体和融池所占比例Sp Sw、Sni与三者光谱反射率的关系，建立多元线性方程组；Si^ri ( λ D +Sw*rw( λ !) +Sm*rm( λ = R(X1)S1=Kr1 ( λ 2) +Sw*rw ( λ 2) +Sn*rn ( λ 2) = R ( λ 2)S1^r1 ( λ 3) +Sw*rw( λ 3) +Sn*rn( λ 3) = R(A3)SJSJSni= IScKr(X1)nR(A2)nR(A3)分别是像兀在λ2、λ 3波段的反射率值；第五步、从图像中选取5% -10%数量的像元，以作为人工神经网络的训练样本，利用选取的像元直接求解带约束条件的方程组，获得的解为每个像元内冰雪、水体和融池所占比例，作为人工神经网络的训练样本集；第六步、构建神经网络模型，输入层为MODIS反射率数据的λ 1、λ2、λ3波段；输出层为冰雪、水体和融池所占比例的图像；第七步、运行神经网络模型，快速优化求解方程组，获得结果为冰雪、水体和融池分别在像元内所占的比例Sp sw、sn；第八步、通过公式Sni* (1-Sw)，计算得到目标区图像内的融池覆盖相对比例，根据融池覆盖相对比例绘制获得目标区融池范围分布图。本专利技术方法使用的数据源为公开下载的中分辨率成像光谱仪(MODIS)的表面反射率8天合成数据(M0D09A1)。M0D09A1包括7个波段，波长459_2155nm，空间分辨率500m。数据范围在O至I之间，表示了该波段上的表面反射率值。在极地海冰区域，图像中每个像元的表面反射率由像元内的冰雪、融池和水体的反射率值共同组成。本专利技术基于神经网络模型的海冰融池提取方法，还具有如下进一步的特征:1、第一步遥感影像预处理前，先在遥感图像处理软件中确定MODIS表面反射率产品M0D09A1影像的投影信息和地理位置信息，将原影像的正弦曲线投影转换为极地方位投影。2、第二步中A1波段的波长范围为620-670nm，λ 2波段的波长范围为841_876nm，λ 3波段的波长范围为459-479nm03、第五步中对于少量像元使用拟牛顿法(BFGS算法，SPBroyden-Fletcher-Goldfarb-Shanno算法)求解方程组，得到各个像元中冰雪、水体和融池分别所占的比例。4、第六步中神经网络模型为多层前馈感知器，由一个输入层、两层隐藏层和一个输出层组成，第一层隐藏层包含9个节点，第二层隐藏层包含27个节点。5、第七步中运行神经网络模型所使用的参数设置为:训练停止条件的RMS误差值为0.1 ;训练的迭代次数为10000次。6、第八步计算融池覆盖相对比例时，将融池在像元内所占比例换算成相对非水体面积而言的值，绘制融池范围分布图时，按照融池覆盖比例值所在区间划分为若干等级，并赋以不同等级的灰度形成图像。本专利技术的有益效果是:获取极地海冰融池的时空分布信息是一项意义重大的研究与技术工作。本专利技术实现了，利用不同类型地表光谱特征的差异，通过神经网络模型快速求解像元分解的多元线性方程组，确定图像中冰雪、水体和融池所占比例，提取融池范围分布。与目视解译和传统监督分类方法相比，本专利技术的自动化程度、分类速度和精度都大为提高。具体有益效果如下:第一，本专利技术成功提取了极地海冰融池范围分布信息，可进一步应用于极地海冰变化监测与趋势预测、不同消融状态下的冰雪反照率变化研究、海冰对全球气候变化的响应机制等研究。第二，本专利技术使用的数据获取便捷，通过申请下载公开的MODIS数据，即可进行极地海冰融池研究工作，无须专门购买数据。第三，本专利技术利用冰雪、水体和融池在遥感图像上光谱特征的差异，基于混合像元分解理论建立各自所占比例与反射率的定量关系，全自动完成混合像元模型下的遥感图像分类，精确提取融池覆盖比例信息。全过程无需人工介入，在确定神经网络模型及其参数后即可直接进行分类。第四，本专利技术构建了神经网络模型以完成混合像元分解计算，将原始MODIS波段数据输入神经网络，自动求解带约束条件的多元线性方程组，即可输出冰雪、水体和融池所占比例的图像。过程优化，运算处理快速，计算结果的精度符合要求。该方法还有一重意义:由于优化并加速了运算过程，可以完成大量遥感数据的分类处理，实现大范围、长时间序列的海冰融池提取。【附图说明】本文档来自技高网...
【技术保护点】
一种基于神经网络模型的极地海冰融池提取方法，其步骤包括：第一步、将MODIS表面反射率产品MOD09A1由正弦投影转换为极地方位投影；第二步、结合MOD09A1产品中的陆地掩膜和云掩膜数据，提取目标区域的图像，并选取λ1、λ2、λ3波段的表面反射率数据；第三步、根据冰雪、水体和融池的光谱响应曲线，得到三者分别在λ1、λ2、λ3波段的标准光谱反射率值：在λ1波段：冰雪的光谱反射率ri(λ1)＝0.95，水体的光谱反射率rw(λ1)＝0.08，融池的光谱反射率rm(λ1)＝0.16；在λ2波段：冰雪的光谱反射率ri(λ2)＝0.87，水体的光谱反射率rw(λ2)＝0.08，融池的光谱反射率rm(λ2)＝0.07；在λ3波段：冰雪的光谱反射率ri(λ3)＝0.95，水体的光谱反射率rw(λ3)＝0.08，融池的光谱反射率rm(λ3)＝0.22；第四步、对图像中的每个像元，根据其中冰雪、水体和融池所占比例Si、Sw、S,与三者光谱反射率的关系，建立多元线性方程组；Si*ri(λ1)+Sw*rw(λ1)+Sm*rm(λ1)＝R(λ1)Si*ri(λ2)+Sw*rw(λ2)+Sm*rm(λ2)＝R(λ2)Si*ri(λ3)+Sw*rw(λ3)+Sm*rm(λ3)＝R(λ3)Si+Sw+Sm＝1式中，R(λ1)、R(λ2)、R(λ3)分别是像元在λ1、λ2、λ3波段的反射率值；第五步、从图像中选取5％‑10％数量的像元，以作为人工神经网络的训练样本，利用选取的像元直接求解带约束条件的方程组，获得的解为每个像元内冰雪、水体和融池所占比例，作为人工神经网络的训练样本集；第六步、构建神经网络模型，输入层为MODIS反射率数据的λ1、λ2、λ3波段；输出层为冰雪、水体和融池所占比例的图像；第七步、运行神经网络模型,快速优化求解方程组，获得结果为冰雪、水体和融池分别在像元内所占的比例Si、Sw、Sm；第八步、通过公式Sm*(1‑Sw)，计算得到目标区图像内的融池覆盖相对比例，根据融池覆盖相对比例绘制获得目标区融池范围分布图。...
【技术特征摘要】
【专利技术属性】
技术研发人员：柯长青，金鑫，邵珠德，钱昊，
申请(专利权)人：南京大学，
类型：发明
国别省市：江苏;32
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。