基于自适应正则参数全变分算法的干涉高光谱图像分解方法与流程

文档序号:11833032阅读:来源:国知局

技术特征:

1.一种基于自适应正则参数全变分算法的干涉高光谱图像分解方法,其特征在于:包括以下步骤:

⑴读取一帧待分解的干涉高光谱图像X,设置正则参数取值的最小值和最大值;

⑵设置XB表示背景层,设置XI表示干涉条纹层,二者与X的关系是:XI=X-XB

⑶采用遗传算法来寻找最优正则化参数Opt_β;

⑷利用寻求最优背景层XB,然后根据XI=X-XB得到干涉条纹层;

⑸完成图像分解。

2.根据权利要求1所述的一种基于自适应正则参数全变分算法的干涉高光谱图像分解方法,其特征在于:步骤⑶所述Opt_β=GA(β,FitFun(β)),GA为遗传算法,FitFun为遗传算法的适应度函数。

3.根据权利要求2所述的一种基于自适应正则参数全变分算法的干涉高光谱图像分解方法,其特征在于:所述GA的适应度函数FitFun的计算包括以下步骤:

⑴设置迭代次数ini_num,且使XB=X;

⑵进行第一次迭代,计算

<mrow> <msub> <mi>X</mi> <mi>B</mi> </msub> <mo>=</mo> <msub> <mi>X</mi> <mi>B</mi> </msub> <mo>-</mo> <mi>&eta;</mi> <mfrac> <mrow> <mo>&part;</mo> <mo>{</mo> <munder> <mo>&Integral;</mo> <mrow> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>&Element;</mo> <mi>X</mi> </mrow> </munder> <msqrt> <msup> <mrow> <mo>(</mo> <mfrac> <mrow> <mo>&part;</mo> <mrow> <mo>(</mo> <mi>X</mi> <mo>(</mo> <mrow> <mi>x</mi> <mo>,</mo> <mi>y</mi> </mrow> <mo>)</mo> <mo>-</mo> <msub> <mi>X</mi> <mi>B</mi> </msub> <mo>(</mo> <mrow> <mi>x</mi> <mo>,</mo> <mi>y</mi> </mrow> <mo>)</mo> <mo>)</mo> </mrow> </mrow> <mrow> <mo>&part;</mo> <mi>y</mi> </mrow> </mfrac> <mo>)</mo> </mrow> <mn>2</mn> </msup> </msqrt> <mo>+</mo> <mi>&beta;</mi> <munder> <mo>&Integral;</mo> <mi>X</mi> </munder> <msqrt> <msup> <mrow> <mo>(</mo> <mfrac> <mrow> <mo>&part;</mo> <msub> <mi>X</mi> <mi>B</mi> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> </mrow> <mrow> <mo>&part;</mo> <mi>x</mi> </mrow> </mfrac> <mo>)</mo> </mrow> <mn>2</mn> </msup> </msqrt> <mo>}</mo> </mrow> <mrow> <mo>&part;</mo> <msub> <mi>X</mi> <mi>B</mi> </msub> </mrow> </mfrac> </mrow>

⑶迭代次数加一,将上一步骤计算出的结果带入公式继续进行迭代计算;

⑷判断迭代次数是否达到最大值,如果达到则进入下一步骤,如果未达到则返回上一步骤继续迭代计算;

⑸计算FitFun=abs(TVx(XB))+abs(TVy(X-XB))}。

4.根据权利要求1或2或3所述的一种基于自适应正则参数全变分算法的干涉高光谱图像分解方法,其特征在于:对中的XB求导,再使用最速下降法进行迭代计算,求出最优解XB,所述最速下降法公式是:

<mrow> <msub> <mi>X</mi> <mi>B</mi> </msub> <mo>=</mo> <msub> <mi>X</mi> <mi>B</mi> </msub> <mo>-</mo> <mi>&eta;</mi> <mfrac> <mrow> <mo>&part;</mo> <mrow> <mo>{</mo> <mrow> <mrow> <msub> <mo>&Integral;</mo> <mrow> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>&Element;</mo> <mi>X</mi> </mrow> </msub> <msqrt> <msup> <mrow> <mo>(</mo> <mfrac> <mrow> <mo>&part;</mo> <mrow> <mo>(</mo> <mi>X</mi> <mo>(</mo> <mrow> <mi>x</mi> <mo>,</mo> <mi>y</mi> </mrow> <mo>)</mo> <mo>-</mo> <msub> <mi>X</mi> <mi>B</mi> </msub> <mo>(</mo> <mrow> <mi>x</mi> <mo>,</mo> <mi>y</mi> </mrow> <mo>)</mo> <mo>)</mo> </mrow> </mrow> <mrow> <mo>&part;</mo> <mi>y</mi> </mrow> </mfrac> <mo>)</mo> </mrow> <mn>2</mn> </msup> </msqrt> </mrow> <mo>+</mo> <mi>O</mi> <mi>p</mi> <mi>t</mi> <mo>_</mo> <mi>&beta;</mi> <mrow> <msub> <mo>&Integral;</mo> <mi>x</mi> </msub> <msqrt> <msup> <mrow> <mo>(</mo> <mfrac> <mrow> <mo>&part;</mo> <msub> <mi>X</mi> <mi>B</mi> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> </mrow> <mrow> <mo>&part;</mo> <mi>x</mi> </mrow> </mfrac> <mo>)</mo> </mrow> <mn>2</mn> </msup> </msqrt> </mrow> </mrow> <mo>}</mo> </mrow> </mrow> <mrow> <mo>&part;</mo> <msub> <mi>X</mi> <mi>B</mi> </msub> </mrow> </mfrac> </mrow>

在最速下降法公式中:Opt_β为遗传算法中计算出的最优正则参数,η为步长,第一次迭代计算时,使XB的取值为X,第二次迭代时将第一次迭代的结果代入,然后下一次的迭代将上一次迭代的结果代入进行计算,循环后得到最优背景层XB,然后根据XI=X-XB得到干涉条纹层。

当前第2页1 2 3 
网友询问留言 已有0条留言
  • 还没有人留言评论。精彩留言会获得点赞!
1