发明名称 一种基于压缩感知的FPM方法
摘要 一种基于压缩感知的FPM方法,包括以下步骤:1)用FPM平台采集不同光照下的低分辨率的图像r<sub>i</sub>(x,y);2)基于不同角度光照下所得到的图像其实是正射下的图像在频域进行平移得到的特点,对采集到的图像ri(x,y)建立约束,基于该约束,根据压缩感知构造求解最优问题的优化问题;3)通过迭代的方法求解所述优化问题得到稀疏系数α,将α与过完备字典相乘得到最终结果。本方法利用压缩感知技术的优势,将原有的FPM算法进行了数学抽象,将频域迭代的方法抽象成求解最优解的问题,从一个新的角度来解决图像超分辨率重建问题,提高了算法的重建效果。
申请公布号 CN104200449B 申请公布日期 2016.05.25
申请号 CN201410420927.0 申请日期 2014.08.25
申请人 清华大学深圳研究生院 发明人 张永兵;蒋伟鑫;戴琼海
分类号 G06T5/50(2006.01)I 主分类号 G06T5/50(2006.01)I
代理机构 深圳市汇力通专利商标代理有限公司 44257 代理人 李保明;张慧芳
主权项 一种基于压缩感知的FPM方法,其特征在于包括以下步骤:1)用FPM平台采集不同光照下的低分辨率的图像r<sub>i</sub>(x,y),包括:控制LED阵列中LED逐点发光,成像装置选择当前LED对应的曝光时间采集图像r<sub>i</sub>(x,y);2)基于不同角度光照下所得到的图像其实是正射下的图像在频域进行平移得到的特点,对采集到的图像ri(x,y)建立约束<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>r</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>=</mo><mo>|</mo><mo>|</mo><msup><mi>F</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mo>{</mo><msub><mi>C</mi><mi>i</mi></msub><mo>*</mo><mi>F</mi><mo>{</mo><mi>H</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>)</mo></mrow><mo>}</mo><mo>*</mo><msubsup><mi>C</mi><mi>i</mi><mi>T</mi></msubsup><mo>}</mo><mo>|</mo><msub><mo>|</mo><mn>2</mn></msub></mrow>]]></math><img file="FDA0000931276510000011.GIF" wi="863" he="78" /></maths>其中,H(x,y)为最终结果,C<sub>i</sub>为提取矩阵,<img file="FDA0000931276510000012.GIF" wi="56" he="65" />为提取矩阵C<sub>i</sub>的转置;基于该约束,根据压缩感知构造求解最优问题的优化问题:<maths num="0002" id="cmaths0002"><math><![CDATA[<mfenced open = "" close = ""><mtable><mtr><mtd><mrow><mi>m</mi><mi>i</mi><mi>n</mi><mo>|</mo><mo>|</mo><mi>&alpha;</mi><mo>|</mo><msub><mo>|</mo><mn>1</mn></msub></mrow></mtd><mtd><mrow><mi>s</mi><mo>.</mo><mi>t</mi><mo>.</mo><mfenced open = "{" close = ""><mtable><mtr><mtd><mrow><msub><mi>r</mi><mn>1</mn></msub><mrow><mo>(</mo><mrow><mi>x</mi><mo>,</mo><mi>y</mi></mrow><mo>)</mo></mrow><mo>=</mo><mo>|</mo><mo>|</mo><msup><mi>F</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mrow><mo>{</mo><mrow><msub><mi>C</mi><mn>1</mn></msub><mo>*</mo><mi>F</mi><mrow><mo>{</mo><mrow><mi>&phi;</mi><mi>&alpha;</mi></mrow><mo>}</mo></mrow><mo>*</mo><msubsup><mi>C</mi><mn>1</mn><mi>T</mi></msubsup></mrow><mo>}</mo></mrow><mo>|</mo><msub><mo>|</mo><mn>2</mn></msub></mrow></mtd></mtr><mtr><mtd><mrow><msub><mi>r</mi><mn>2</mn></msub><mrow><mo>(</mo><mrow><mi>x</mi><mo>,</mo><mi>y</mi></mrow><mo>)</mo></mrow><mo>=</mo><mo>|</mo><mo>|</mo><msup><mi>F</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mrow><mo>{</mo><mrow><msub><mi>C</mi><mn>2</mn></msub><mo>*</mo><mi>F</mi><mrow><mo>{</mo><mrow><mi>&phi;</mi><mi>&alpha;</mi></mrow><mo>}</mo></mrow><mo>*</mo><msubsup><mi>C</mi><mn>2</mn><mi>T</mi></msubsup></mrow><mo>}</mo></mrow><mo>|</mo><msub><mo>|</mo><mn>2</mn></msub></mrow></mtd></mtr><mtr><mtd><mtable><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr></mtable></mtd></mtr><mtr><mtd><mrow><msub><mi>r</mi><mi>L</mi></msub><mrow><mo>(</mo><mrow><mi>x</mi><mo>,</mo><mi>y</mi></mrow><mo>)</mo></mrow><mo>=</mo><mo>|</mo><mo>|</mo><msup><mi>F</mi><mrow><mo>-</mo><mn>1</mn></mrow></msup><mrow><mo>{</mo><mrow><msub><mi>C</mi><mi>L</mi></msub><mo>*</mo><mi>F</mi><mrow><mo>{</mo><mrow><mi>&phi;</mi><mi>&alpha;</mi></mrow><mo>}</mo></mrow><mo>*</mo><msubsup><mi>C</mi><mi>L</mi><mi>T</mi></msubsup></mrow><mo>}</mo></mrow><mo>|</mo><msub><mo>|</mo><mn>2</mn></msub></mrow></mtd></mtr></mtable></mfenced></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000931276510000013.GIF" wi="1211" he="359" /></maths>其中,L等于LED阵列中LED的数量,φ为过完备字典,α是稀疏系数;3)通过迭代的方法求解所述优化问题得到稀疏系数α,将α与过完备字典相乘得到最终结果;其中,采用以下方法确定LED阵列中各个LED对应的成像装置曝光时间:首先确定中心LED的位置:固定一个曝光时间,LED阵列上的灯逐点发光,用成像装置采集每一张图像,计算图像的亮度分布,估计中心LED的位置;然后确定不同区域的曝光时间:选定几个不同的曝光时间,LED阵列上的灯逐点发光,用成像装置采集每一张图像,筛选每个LED对应的图像,确定每个LED对应需要的曝光时间;主要的原则是:a)曝光时间尽可能选择短的,但是一定保证基本保留了图像信息;b)与中心LED距离相同的LED所需要的曝光时间应该一样;所述中心LED为样本正下方的LED。
地址 518055 广东省深圳市南山区西丽大学城清华校区