发明名称 基于全变分范数图像分块梯度计算的图像重构方法
摘要 本发明涉及一种基于全变分范数图像分块梯度计算的图像重构方法,包括图像分块处理、初始化、测量值获取、图像块重构、合并图像块等步骤;在对分块图像进行重构时,除了其最上面一横排图像块和最左面的一竖排图像块,其余的各小块充分地利用已完成重构的上方和左方的图像块的边缘信息,因此,本发明可以有效地去除块效应,减少重建图像所消耗的时间,提高图像的重构质量。
申请公布号 CN106097263A 申请公布日期 2016.11.09
申请号 CN201610392351.0 申请日期 2016.06.03
申请人 江苏大学 发明人 宋雪桦;张鑫晟;曹曦文;孔尧;孙海威;叶润武;高亚红;孙旭;化瑞
分类号 G06T5/00(2006.01)I 主分类号 G06T5/00(2006.01)I
代理机构 南京苏高专利商标事务所(普通合伙) 32204 代理人 罗敏
主权项 一种基于全变分范数图像分块梯度计算的图像重构方法,其特征在于包括如下步骤:1)图像分块:将大小为M×N的图像u分成大小为n×n的图像块,其中,n为M、N的公约数,满足s=M/n,t=N/n,即待处理图像被分成H=s×t个图像块,图像块从左到右,从上到下依次记为K<sub>1</sub>,K<sub>2</sub>…K<sub>H</sub>;2)初始化图像块变量c=1;3)利用测量矩阵对图像块K<sub>c</sub>进行随机测量,得到图像块K<sub>c</sub>随机像素点的像素值矩阵y;4)计算图像块K<sub>c</sub>水平方向和垂直方向的梯度值,其中,若图像块u<sub>c</sub>位于第1行第1列,水平方向和垂直方向的梯度值计算方法如下:<maths num="0001"><math><![CDATA[<mrow><mtable><mtr><mtd><mrow><msub><mi>D</mi><mi>h</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mi>K</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>-</mo><msub><mi>K</mi><mrow><mi>i</mi><mo>-</mo><mn>1</mn><mo>,</mo><mi>j</mi></mrow></msub><mo>,</mo></mrow></mtd><mtd><mrow><mn>1</mn><mo>&lt;</mo><mi>i</mi><mo>&le;</mo><mi>n</mi></mrow></mtd></mtr><mtr><mtd><mrow><mn>0</mn><mo>,</mo></mrow></mtd><mtd><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow></mtd></mtr></mtable></mfenced></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>D</mi><mi>v</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mi>K</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>-</mo><msub><mi>K</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>,</mo></mrow></mtd><mtd><mrow><mn>1</mn><mo>&lt;</mo><mi>i</mi><mo>&le;</mo><mi>n</mi></mrow></mtd></mtr><mtr><mtd><mrow><mn>0</mn><mo>,</mo></mrow></mtd><mtd><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow></mtd></mtr></mtable></mfenced></mrow></mtd></mtr></mtable><mo>;</mo></mrow>]]></math><img file="FDA0001008871290000011.GIF" wi="694" he="302" /></maths>其中,D<sub>h</sub>(i,j)表示像素点(i,j)的水平梯度值,D<sub>v</sub>(i,j)表示像素点(i,j)的垂直梯度值,K<sub>i,j</sub>表示像素点(i,j)的像素值,i和j满足i∈[1,n],j∈[1,n],转步骤8),否则转步骤5);5)若图像块K<sub>c</sub>位于第1行,水平方向和垂直方向的梯度值的计算方法如下:<maths num="0002"><math><![CDATA[<mrow><mtable><mtr><mtd><mrow><msub><mi>D</mi><mi>h</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mi>K</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>-</mo><msub><mi>K</mi><mrow><mi>i</mi><mo>-</mo><mn>1</mn><mo>,</mo><mi>j</mi></mrow></msub><mo>,</mo></mrow></mtd><mtd><mrow><mn>1</mn><mo>&lt;</mo><mi>i</mi><mo>&le;</mo><mi>n</mi></mrow></mtd></mtr><mtr><mtd><mrow><mn>0</mn><mo>,</mo></mrow></mtd><mtd><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow></mtd></mtr></mtable></mfenced></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>D</mi><mi>v</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mi>K</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>-</mo><msub><mi>K</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>,</mo></mrow></mtd><mtd><mrow><mn>1</mn><mo>&lt;</mo><mi>i</mi><mo>&le;</mo><mi>n</mi></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>K</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>-</mo><msub><mi>U</mi><mrow><mi>i</mi><mo>,</mo><mi>n</mi></mrow></msub><mo>,</mo></mrow></mtd><mtd><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow></mtd></mtr></mtable></mfenced></mrow></mtd></mtr></mtable><mo>;</mo></mrow>]]></math><img file="FDA0001008871290000012.GIF" wi="709" he="311" /></maths>其中,U<sub>i,n</sub>用来表示待重构块K<sub>c</sub>的上方图像块的边缘信息,转步骤8),否则转步骤6);6)若图像块K<sub>c</sub>位于第1列,水平方向和垂直方向的梯度值的计算方法如下:<img file="FDA0001008871290000013.GIF" wi="694" he="295" />转步骤8),否则转步骤7);7)水平方向和垂直方向的梯度值的计算方法如下:<maths num="0003"><math><![CDATA[<mrow><mtable><mtr><mtd><mrow><msub><mi>D</mi><mi>h</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mi>K</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>-</mo><msub><mi>K</mi><mrow><mi>i</mi><mo>-</mo><mn>1</mn><mo>,</mo><mi>j</mi></mrow></msub><mo>,</mo></mrow></mtd><mtd><mrow><mn>1</mn><mo>&lt;</mo><mi>i</mi><mo>&le;</mo><mi>n</mi></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>K</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>-</mo><msub><mi>L</mi><mrow><mi>n</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>,</mo></mrow></mtd><mtd><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow></mtd></mtr></mtable></mfenced></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>D</mi><mi>v</mi></msub><mrow><mo>(</mo><mi>i</mi><mo>,</mo><mi>j</mi><mo>)</mo></mrow><mo>=</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mi>K</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>-</mo><msub><mi>K</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi><mo>-</mo><mn>1</mn></mrow></msub><mo>,</mo></mrow></mtd><mtd><mrow><mn>1</mn><mo>&lt;</mo><mi>i</mi><mo>&le;</mo><mi>n</mi></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>K</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>-</mo><msub><mi>U</mi><mrow><mi>i</mi><mo>,</mo><mi>n</mi></mrow></msub><mo>,</mo></mrow></mtd><mtd><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow></mtd></mtr></mtable></mfenced></mrow></mtd></mtr></mtable><mo>;</mo></mrow>]]></math><img file="FDA0001008871290000014.GIF" wi="694" he="301" /></maths>其中,L<sub>n,j</sub>用来表示待重构块K<sub>c</sub>的左方图像块的边缘信息,转步骤8);8)对图像块K<sub>c</sub>进行压缩传感重构;9)修改图像块变量:c=c+1,转步骤10);10)如果c&lt;=H,转步骤3),否则转步骤11);11)合并图像块,按照拆分顺序将重构后的图像块合并,整个重构过程结束。
地址 212013 江苏省镇江市学府路301号