发明名称 一种结构自适应和结构保持的图像局部扭曲方法
摘要 本发明公开了一种结构自适应和结构保持的图像局部扭曲方法。该方法在单尺度下,首先分析图像的局部边缘方向和一致性,据此设计方向平滑核函数,使用平滑核函数对初始偏移场的边缘分量做方向平滑,最后使用平滑后的偏移场扭曲图像。为同时保持大尺度和小尺度结构,构建了多尺度的图像金字塔和偏移场金字塔,自上而下平滑偏移场,在每层偏移场上使用迭代平滑提高平滑效果,将上一层平滑后的偏移场作为下一层的初始偏移场。本发明通过分析和利用图像中的感知显著结构信息,达到不改变此类重要结构而扭曲图像的效果,提高了扭曲后图像的主观感知质量,并且不影响其对水印系统等的攻击效果。
申请公布号 CN103793880B 申请公布日期 2016.09.21
申请号 CN201410075833.4 申请日期 2014.03.04
申请人 山东科技大学 发明人 颜斌;杨红梅;崔鑫;郭银景;王卓鹏;郝建军;王凤瑛;张同军;张仁彦
分类号 G06T3/00(2006.01)I 主分类号 G06T3/00(2006.01)I
代理机构 济南舜源专利事务所有限公司 37205 代理人 王连君
主权项 一种结构自适应和结构保持的图像局部扭曲方法,其特征在于,包括以下处理步骤:a分析图像局部结构构建图像金字塔并分析各层局部结构;b平滑偏移场构建偏移场金字塔并自顶向下对偏移场做方向平滑;c扭曲图像使用平滑后的偏移场扭曲图像;所述步骤a具体包括以下处理步骤:a1构造图像的高斯金字塔;a2对每层图像分析其局部边缘方向并计算出一个一致性度量;所述步骤a2具体包括以下处理步骤:a21首先计算图像的梯度<img file="FDA0000992962960000011.GIF" wi="562" he="63" />并在一个小邻域内计算结构张量J:<maths num="0001"><math><![CDATA[<mrow><mi>J</mi><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><mrow><munder><mo>&Sigma;</mo><mrow><mo>(</mo><mi>k</mi><mo>,</mo><mi>l</mi><mo>)</mo><mo>&Element;</mo><mi>N</mi><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow></munder><msub><mi>I</mi><mi>x</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>,</mo><mi>l</mi><mo>)</mo></mrow><mo>&CenterDot;</mo><msub><mi>I</mi><mi>x</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>,</mo><mi>l</mi><mo>)</mo></mrow></mrow></mtd><mtd><mrow><munder><mo>&Sigma;</mo><mrow><mo>(</mo><mi>k</mi><mo>,</mo><mi>l</mi><mo>)</mo><mo>&Element;</mo><mi>N</mi><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow></munder><msub><mi>I</mi><mi>x</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>,</mo><mi>l</mi><mo>)</mo></mrow><mo>&CenterDot;</mo><msub><mi>I</mi><mi>y</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>,</mo><mi>l</mi><mo>)</mo></mrow></mrow></mtd></mtr><mtr><mtd><mrow><munder><mo>&Sigma;</mo><mrow><mo>(</mo><mi>k</mi><mo>,</mo><mi>l</mi><mo>)</mo><mo>&Element;</mo><mi>N</mi><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow></munder><msub><mi>I</mi><mi>x</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>,</mo><mi>l</mi><mo>)</mo></mrow><mo>&CenterDot;</mo><msub><mi>I</mi><mi>y</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>,</mo><mi>l</mi><mo>)</mo></mrow></mrow></mtd><mtd><mrow><munder><mo>&Sigma;</mo><mrow><mo>(</mo><mi>k</mi><mo>,</mo><mi>l</mi><mo>)</mo><mo>&Element;</mo><mi>N</mi><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow></munder><msub><mi>I</mi><mi>y</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>,</mo><mi>l</mi><mo>)</mo></mrow><mo>&CenterDot;</mo><msub><mi>I</mi><mi>y</mi></msub><mrow><mo>(</mo><mi>k</mi><mo>,</mo><mi>l</mi><mo>)</mo></mrow></mrow></mtd></mtr></mtable></mfenced><mo>,</mo></mrow>]]></math><img file="FDA0000992962960000012.GIF" wi="1046" he="223" /></maths>其中,(i,j)表示像素位置,I(i,j)表示灰度值,I<sub>x</sub>(i,j)和I<sub>y</sub>(i,j)分别是沿着水平和垂直方向的方向导数,N(i,j)表示以像素(i,j)为中心的一个小邻域;对结构张量J做特征值分解,即:<maths num="0002"><math><![CDATA[<mrow><mi>J</mi><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><mi>e</mi></mtd><mtd><mi>f</mi></mtd></mtr></mtable></mfenced><mfenced open = "[" close = "]"><mtable><mtr><mtd><msub><mi>&lambda;</mi><mi>e</mi></msub></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><msub><mi>&lambda;</mi><mi>f</mi></msub></mtd></mtr></mtable></mfenced><mfenced open = "[" close = "]"><mtable><mtr><mtd><msup><mi>e</mi><mi>T</mi></msup></mtd></mtr><mtr><mtd><msup><mi>f</mi><mi>T</mi></msup></mtd></mtr></mtable></mfenced><mo>,</mo><msub><mi>&lambda;</mi><mi>e</mi></msub><mo>&gt;</mo><msub><mi>&lambda;</mi><mi>f</mi></msub><mo>,</mo></mrow>]]></math><img file="FDA0000992962960000013.GIF" wi="710" he="150" /></maths>得到边缘方向向量e;其中,单位向量e称为边缘方向,而单位向量f称为流向,这两个方向正交,λ<sub>e</sub>表示对应e的特征值,而λ<sub>f</sub>表示对应f的特征值;a22如下计算局部结构的一致性度量κ:<maths num="0003"><math><![CDATA[<mrow><mi>&kappa;</mi><mo>=</mo><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mfrac><mi>&tau;</mi><mrow><msub><mi>&lambda;</mi><mi>e</mi></msub><mo>-</mo><msub><mi>&lambda;</mi><mi>f</mi></msub></mrow></mfrac><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0000992962960000014.GIF" wi="422" he="159" /></maths>其中,τ是阈值参数;最终获得局部结构,即每个像素位置的边缘方向e(i,j)和一致性度量κ(i,j)。
地址 266590 山东省青岛市经济技术开发区前湾港路579号
您可能感兴趣的专利