发明名称 一种单极化SAR假彩色影像生成方法
摘要 本发明公开了一种单极化SAR假彩色影像生成方法,本发明针对单极化SAR通常用灰度显示,无法像光学或全极化SAR一样彩色显示,从而限制了解译效果而用户认可度较低的问题,以地物类型丰富区域的某一单极化,HH、HV或VV的SAR影像的原始幅度、幅度均值、标准差等特征为输入,并以同区域的全极化SAR假彩色影像作为训练的先验知识,采用最小二乘拟合的方法,求取了从单极化影像特征到红绿蓝三个通道的三组多项式参数,使用该三组参数,能将同一传感器该单极化模式下的影像生成类似全极化的假彩色影,效提升单极化SAR影像的可视化效果。
申请公布号 CN106355625A 申请公布日期 2017.01.25
申请号 CN201610779249.6 申请日期 2016.08.31
申请人 甘宗平 发明人 邓少平;梁志华;孙盛;吴泽洪;甘宗平;刘学林;冷海芹
分类号 G06T11/00(2006.01)I 主分类号 G06T11/00(2006.01)I
代理机构 中山市科创专利代理有限公司 44211 代理人 凌信景;胡犇
主权项 一种单极化SAR假彩色影像生成方法,其特征在于包括以下步骤:步骤1:对全极化影像进行Pauli分解<img file="FDA0001104319110000011.GIF" wi="379" he="222" />步骤2:计算Pauli矢量各分量的幅度:A<sub>HH+VV</sub>=|s<sub>HH</sub>+s<sub>VV</sub>|、A<sub>HH‑VV</sub>=|s<sub>HH</sub>‑s<sub>VV</sub>|和A<sub>HV</sub>=|s<sub>HV</sub>|,其中|·|表示复数的模,即SAR后向散射的幅度,按照R=A<sub>HH‑VV</sub>,G=A<sub>HV</sub>,B=A<sub>HH+VV</sub>组合即可生成全极化假彩色影像;步骤3:计算全极化SAR影像中,与待处理极化方式X相同的极化通道的幅度影像:A<sub>X</sub>=|s<sub>X</sub>|,X=HH,HV或VV;步骤4:在局部邻域加权统计A<sub>X</sub>的均值M<sub>X</sub>和标准差V<sub>X</sub>作为新的特征,其中7×7邻域窗口的权矩阵W如下:<tables num="0001" wi="165"><table><tgroup cols="7"><colspec colname="c001" colwidth="14%" /><colspec colname="c002" colwidth="14%" /><colspec colname="c003" colwidth="14%" /><colspec colname="c004" colwidth="14%" /><colspec colname="c005" colwidth="14%" /><colspec colname="c006" colwidth="14%" /><colspec colname="c007" colwidth="16%" /><tbody><row><entry morerows="1">0.5</entry><entry morerows="1">0.5</entry><entry morerows="1">1.0</entry><entry morerows="1">1.5</entry><entry morerows="1">1.0</entry><entry morerows="1">0.5</entry><entry morerows="1">0.5</entry></row><row><entry morerows="1">0.5</entry><entry morerows="1">1.0</entry><entry morerows="1">1.5</entry><entry morerows="1">2.0</entry><entry morerows="1">1.5</entry><entry morerows="1">1.0</entry><entry morerows="1">0.5</entry></row><row><entry morerows="1">1.0</entry><entry morerows="1">1.5</entry><entry morerows="1">2.0</entry><entry morerows="1">2.5</entry><entry morerows="1">2.0</entry><entry morerows="1">1.5</entry><entry morerows="1">1.0</entry></row><row><entry morerows="1">1.5</entry><entry morerows="1">2.0</entry><entry morerows="1">2.5</entry><entry morerows="1">3.0</entry><entry morerows="1">2.5</entry><entry morerows="1">2.0</entry><entry morerows="1">1.5</entry></row><row><entry morerows="1">1.0</entry><entry morerows="1">1.5</entry><entry morerows="1">2.0</entry><entry morerows="1">2.5</entry><entry morerows="1">2.0</entry><entry morerows="1">1.5</entry><entry morerows="1">1.0</entry></row><row><entry morerows="1">0.5</entry><entry morerows="1">1.0</entry><entry morerows="1">1.5</entry><entry morerows="1">2.0</entry><entry morerows="1">1.5</entry><entry morerows="1">1.0</entry><entry morerows="1">0.5</entry></row><row><entry morerows="1">0.5</entry><entry morerows="1">0.5</entry><entry morerows="1">1.0</entry><entry morerows="1">1.5</entry><entry morerows="1">1.0</entry><entry morerows="1">0.5</entry><entry morerows="1">0.5</entry></row></tbody></tgroup></table></tables>均值和标准差的计算分别如下:<maths num="0001"><math><![CDATA[<mrow><msub><mi>M</mi><mi>X</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mrow><munderover><mi>&Sigma;</mi><mrow><msub><mi>i</mi><mn>0</mn></msub><mo>=</mo><mn>1</mn></mrow><mn>7</mn></munderover><munderover><mi>&Sigma;</mi><mrow><msub><mi>j</mi><mn>0</mn></msub><mo>=</mo><mn>1</mn></mrow><mn>7</mn></munderover><mi>W</mi><mrow><mo>(</mo><msub><mi>i</mi><mn>0</mn></msub><mo>,</mo><msub><mi>j</mi><mn>0</mn></msub><mo>)</mo></mrow></mrow></mfrac><munderover><mi>&Sigma;</mi><mrow><msub><mi>i</mi><mn>0</mn></msub><mo>=</mo><mn>1</mn></mrow><mn>7</mn></munderover><munderover><mi>&Sigma;</mi><mrow><msub><mi>j</mi><mn>0</mn></msub><mo>=</mo><mn>1</mn></mrow><mn>7</mn></munderover><mrow><mo>(</mo><mi>W</mi><mo>(</mo><mrow><msub><mi>i</mi><mn>0</mn></msub><mo>,</mo><msub><mi>j</mi><mn>0</mn></msub></mrow><mo>)</mo><msub><mi>A</mi><mi>X</mi></msub><mo>(</mo><mrow><mi>i</mi><mo>+</mo><msub><mi>i</mi><mn>0</mn></msub><mo>-</mo><mn>4</mn><mo>,</mo><mi>j</mi><mo>+</mo><msub><mi>j</mi><mn>0</mn></msub><mo>-</mo><mn>4</mn></mrow><mo>)</mo><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0001104319110000012.GIF" wi="1358" he="215" /></maths><img file="FDA0001104319110000021.GIF" wi="1600" he="312" />步骤5:将A<sub>HH‑VV</sub>、A<sub>HV</sub>、A<sub>HH+VV</sub>线性拉伸到0‑63的整数,并使得拉伸值为0和63像素比例为2%,分别记拉伸后的影像为N<sub>q,R</sub>、N<sub>q,G</sub>和N<sub>q,B</sub>;步骤6:利用全极化假彩色影像最小二乘法求解单极化到假彩色各通道的参数;步骤7:重复步骤6t次,并计算步骤6得到转换参数的数学平均值,记为<img file="FDA0001104319110000022.GIF" wi="291" he="64" />步骤8:提取单极化SAR影像上的特征,并变化到假彩色影像的RGB三个通道;步骤9:基于主成分分析PCA技术,并利用原始单极化影像增强假彩色影像的细节信息;步骤10:对<img file="FDA0001104319110000023.GIF" wi="51" he="73" />拉伸到0‑255的整数,并使得RGB每个通道值为0和255的值,且所占的比例为2%,即完成了假彩色影像的生成。
地址 528437 广东省中山市东区东裕路顺景花园19栋304房