发明名称 一种光谱和空间信息结合的高光谱数据非均一校正方法
摘要 一种光谱和空间信息结合的高光谱数据非均一校正方法,它有六大步骤:一、高光谱数据的读入;二、由边缘检测算子计算每一个单波段图像中非均一特征的边缘;三、由步骤二提取的非均一特征边缘实现图像分割,判定响应非均一探元的位置;四、由步骤三响应非均一探元位置确定参考探元位置,并计算相应的光谱维比率;五、由步骤四探元光谱维比率计算响应非均一探元的校正系数;六、由校正系数实现响应非均一像元值的更新,实现非均一性校正。本发明能稳定的分离图像中的非均一特征,并结合光谱和空间信息获得更可靠的校正结果,它在高光谱图像预处理领域里具有较好的实用价值和广阔的应用前景。
申请公布号 CN102682430B 申请公布日期 2014.07.09
申请号 CN201210125152.5 申请日期 2012.04.25
申请人 北京航空航天大学 发明人 赵慧洁;江澄;贾国瑞
分类号 G06T5/00(2006.01)I;G06T7/00(2006.01)I 主分类号 G06T5/00(2006.01)I
代理机构 北京慧泉知识产权代理有限公司 11232 代理人 王顺荣;唐爱华
主权项 一种光谱和空间信息结合的高光谱数据非均一校正方法,其特征在于:该方法具体步骤如下:步骤一:高光谱数据的读入;步骤二:由边缘检测算子计算每一个单波段图像中非均一特征的边缘;步骤三:由步骤二提取的非均一特征边缘实现图像分割,判定响应非均一探元的位置;步骤四:由步骤三响应非均一探元位置确定参考探元位置,并计算相应的光谱维比率;步骤五:由步骤四探元光谱维比率计算响应非均一探元的校正系数;步骤六:由校正系数实现响应非均一像元值的更新,实现非均一性校正;其中,步骤三中所述的实现图像分割,判定响应非均一探元的位置,判定准则如下:<img file="FDA0000486487560000011.GIF" wi="837" he="148" />式中,Loc为1表示此探元响应非均一,为0表示响应均一;R是非均一特征边缘上标记为边缘点的数目,P是最大连续边缘点的数目;M是图像的行数;Res是数据的空间分辨率;其中,步骤四中所述参考探元是指与非均一探元最近邻的均一探元,光谱维比率由不同波段的辐亮度计算得到:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>r</mi><mrow><mo>(</mo><msub><mi>x</mi><mn>0</mn></msub><mo>,</mo><msub><mi>y</mi><mi>j</mi></msub><mo>,</mo><msub><mi>k</mi><mn>0</mn></msub><mo>,</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mi>L</mi><mrow><mo>(</mo><msub><mi>x</mi><mn>0</mn></msub><mo>,</mo><msub><mi>y</mi><mi>j</mi></msub><mo>,</mo><msub><mi>k</mi><mn>0</mn></msub><mo>+</mo><mi>n</mi><mo>)</mo></mrow></mrow><mrow><mi>L</mi><mrow><mo>(</mo><msub><mi>x</mi><mn>0</mn></msub><mo>,</mo><msub><mi>y</mi><mi>j</mi></msub><mo>,</mo><msub><mi>k</mi><mn>0</mn></msub><mo>)</mo></mrow></mrow></mfrac></mrow>]]></math><img file="FDA0000486487560000021.GIF" wi="777" he="184" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mover><mi>r</mi><mo>&OverBar;</mo></mover><mrow><mo>(</mo><msub><mover><mi>x</mi><mo>&OverBar;</mo></mover><mn>0</mn></msub><mo>,</mo><msub><mi>y</mi><mi>j</mi></msub><mo>,</mo><msub><mi>k</mi><mn>0</mn></msub><mo>,</mo><mi>n</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mi>L</mi><mrow><mo>(</mo><msub><mover><mi>x</mi><mo>&OverBar;</mo></mover><mn>0</mn></msub><mo>,</mo><msub><mi>y</mi><mi>j</mi></msub><mo>,</mo><msub><mi>k</mi><mn>0</mn></msub><mo>+</mo><mi>n</mi><mo>)</mo></mrow></mrow><mrow><mi>L</mi><mrow><mo>(</mo><msub><mover><mi>x</mi><mo>&OverBar;</mo></mover><mn>0</mn></msub><mo>,</mo><msub><mi>y</mi><mi>j</mi></msub><mo>,</mo><msub><mi>k</mi><mn>0</mn></msub><mo>)</mo></mrow></mrow></mfrac></mrow>]]></math><img file="FDA0000486487560000022.GIF" wi="785" he="184" /></maths>式中,r(x<sub>0</sub>,y<sub>j</sub>,k<sub>0</sub>,n)表示探元(x<sub>0</sub>,k<sub>0</sub>)的光谱维比率,<img file="FDA0000486487560000023.GIF" wi="290" he="68" />表示参考探元<img file="FDA0000486487560000024.GIF" wi="154" he="71" />的光谱维比率;L(x<sub>0</sub>,y<sub>j</sub>,k<sub>0</sub>+n)、L(x<sub>0</sub>,y<sub>j</sub>,k<sub>0</sub>)分别表示探元(x<sub>0</sub>,k<sub>0</sub>)在第k<sub>0</sub>+n波段、第k<sub>0</sub>波段的辐亮度值,<img file="FDA0000486487560000025.GIF" wi="556" he="77" />分别表示探元<img file="FDA0000486487560000026.GIF" wi="150" he="71" />在第k<sub>0</sub>+n波段、第k<sub>0</sub>波段的辐亮度值,其中y<sub>j</sub>是空间维的行序号,j=1,2,……,M,并且M是图像的行数;其中,步骤五中求取响应非均一探元的校正系数计算公式如下:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>a</mi><mrow><mo>(</mo><msub><mi>x</mi><mn>0</mn></msub><mo>,</mo><msub><mi>k</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mi>Std</mi><mrow><mo>(</mo><msub><mi>x</mi><mn>0</mn></msub><mo>,</mo><msub><mi>k</mi><mn>0</mn></msub><mo>)</mo></mrow></mrow><mrow><mi>Std</mi><mrow><mo>(</mo><msub><mover><mi>x</mi><mo>&OverBar;</mo></mover><mn>0</mn></msub><mo>,</mo><msub><mi>k</mi><mn>0</mn></msub><mo>)</mo></mrow></mrow></mfrac></mrow>]]></math><img file="FDA0000486487560000027.GIF" wi="456" he="142" /></maths><maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>b</mi><mrow><mo>(</mo><msub><mi>x</mi><mn>0</mn></msub><mo>,</mo><msub><mi>k</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>=</mo><mi>Mean</mi><mrow><mo>(</mo><msub><mi>x</mi><mn>0</mn></msub><mo>,</mo><msub><mi>k</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>-</mo><mfrac><mrow><mi>Std</mi><mrow><mo>(</mo><msub><mi>x</mi><mn>0</mn></msub><mo>,</mo><msub><mi>k</mi><mn>0</mn></msub><mo>)</mo></mrow></mrow><mrow><mi>Std</mi><mrow><mo>(</mo><msub><mover><mi>x</mi><mo>&OverBar;</mo></mover><mn>0</mn></msub><mo>,</mo><msub><mi>k</mi><mn>0</mn></msub><mo>)</mo></mrow></mrow></mfrac><mi>Mean</mi><mrow><mo>(</mo><msub><mover><mi>x</mi><mo>&OverBar;</mo></mover><mn>0</mn></msub><mo>,</mo><msub><mi>k</mi><mn>0</mn></msub><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000486487560000028.GIF" wi="1014" he="142" /></maths>式中,a(x<sub>0</sub>,k<sub>0</sub>)、b(x<sub>0</sub>,k<sub>0</sub>)是非均一探元(x<sub>0</sub>,k<sub>0</sub>)的校正系数,分别表征增益和偏置;Std(x<sub>0</sub>,k<sub>0</sub>)、<img file="FDA0000486487560000029.GIF" wi="226" he="71" />分别表示非均一探元光谱维比率、参考探元光谱维比率的标准差;Mean(x<sub>0</sub>,k<sub>0</sub>)、<img file="FDA00004864875600000210.GIF" wi="270" he="71" />分别表示非均一探元光谱维比率、参考探元光谱维比率的均值。
地址 100191 北京市海淀区学院路37号