主权项 |
一种基于神经网络模型的极地海冰融池提取方法,其步骤包括:第一步、将MODIS表面反射率产品MOD09A1由正弦投影转换为极地方位投影;第二步、结合MOD09A1产品中的陆地掩膜和云掩膜数据,提取目标区域的图像,并选取λ<sub>1</sub>、λ<sub>2</sub>、λ<sub>3</sub>波段的表面反射率数据;第三步、根据冰雪、水体和融池的光谱响应曲线,得到三者分别在λ<sub>1</sub>、λ<sub>2</sub>、λ<sub>3</sub>波段的标准光谱反射率值:在λ<sub>1</sub>波段:冰雪的光谱反射率r<sub>i</sub>(λ<sub>1</sub>)=0.95,水体的光谱反射率r<sub>w</sub>(λ<sub>1</sub>)=0.08,融池的光谱反射率r<sub>m</sub>(λ<sub>1</sub>)=0.16;在λ<sub>2</sub>波段:冰雪的光谱反射率r<sub>i</sub>(λ<sub>2</sub>)=0.87,水体的光谱反射率r<sub>w</sub>(λ<sub>2</sub>)=0.08,融池的光谱反射率r<sub>m</sub>(λ<sub>2</sub>)=0.07;在λ<sub>3</sub>波段:冰雪的光谱反射率r<sub>i</sub>(λ<sub>3</sub>)=0.95,水体的光谱反射率r<sub>w</sub>(λ<sub>3</sub>)=0.08,融池的光谱反射率r<sub>m</sub>(λ<sub>3</sub>)=0.22;第四步、对图像中的每个像元,根据其中冰雪、水体和融池所占比例S<sub>i</sub>、S<sub>w</sub>、S<sub>,</sub>与三者光谱反射率的关系,建立多元线性方程组;S<sub>i*</sub>r<sub>i</sub>(λ<sub>1</sub>)+S<sub>w</sub>*r<sub>w</sub>(λ<sub>1</sub>)+S<sub>m</sub>*r<sub>m</sub>(λ<sub>1</sub>)=R(λ<sub>1</sub>)S<sub>i</sub>*r<sub>i</sub>(λ<sub>2</sub>)+S<sub>w*</sub>r<sub>w</sub>(λ<sub>2</sub>)+S<sub>m</sub>*r<sub>m</sub>(λ<sub>2</sub>)=R(λ<sub>2</sub>)S<sub>i</sub>*r<sub>i</sub>(λ<sub>3</sub>)+S<sub>w</sub>*r<sub>w</sub>(λ<sub>3</sub>)+S<sub>m</sub>*r<sub>m</sub>(λ<sub>3</sub>)=R(λ<sub>3</sub>)S<sub>i</sub>+S<sub>w</sub>+S<sub>m</sub>=1式中,R(λ<sub>1</sub>)、R(λ<sub>2</sub>)、R(λ<sub>3</sub>)分别是像元在λ<sub>1</sub>、λ<sub>2</sub>、λ<sub>3</sub>波段的反射率值;第五步、从图像中选取5%‑10%数量的像元,以作为人工神经网络的训练样本,利用选取的像元直接求解带约束条件的方程组,获得的解为每个像元内冰雪、水体和融池所占比例,作为人工神经网络的训练样本集;第六步、构建神经网络模型,输入层为MODIS反射率数据的λ<sub>1</sub>、λ<sub>2</sub>、λ<sub>3</sub>波段;输出层为冰雪、水体和融池所占比例的图像;第七步、运行神经网络模型,快速优化求解方程组,获得结果为冰雪、水体和融池分别在像元内所占的比例S<sub>i</sub>、S<sub>w</sub>、S<sub>m</sub>;第八步、通过公式S<sub>m</sub>*(1‑S<sub>w</sub>),计算得到目标区图像内的融池覆盖相对比例,根据融池覆盖相对比例绘制获得目标区融池范围分布图。 |