CN104796609B - 基于最优哈达玛编码的大视场高分辨率显微成像方法 - Google Patents
基于最优哈达玛编码的大视场高分辨率显微成像方法 Download PDFInfo
- Publication number
- CN104796609B CN104796609B CN201510186030.0A CN201510186030A CN104796609B CN 104796609 B CN104796609 B CN 104796609B CN 201510186030 A CN201510186030 A CN 201510186030A CN 104796609 B CN104796609 B CN 104796609B
- Authority
- CN
- China
- Prior art keywords
- mrow
- resolution
- msub
- mtd
- optimal
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Landscapes
- Transforming Light Signals Into Electric Signals (AREA)
- Image Processing (AREA)
Abstract
本发明公开了一种基于最优哈达玛编码的大视场高分辨率显微成像方法,首先利用计算机生成所需的最优哈达玛编码的矩阵,然后采用LED阵列作为显微镜光源,显示最优哈达玛编码图案,并用相机拍摄一系列待测样品在不同照明角度下的低分辨率显微图像,再根据最优哈达玛矩阵求解出待测样品受到单个角度照明时的低分辨率显微图像,最后利用频域合成孔径技术将多幅大视场低分辨率图像合成一幅大视场高分辨率图像。本发明既能降低相机所需的曝光时间,提高系统的图像采集速度,又能抑制噪声,提高系统采集的图像质量。
Description
技术领域
本发明属于光学成像技术,特别是一种基于最优哈达玛编码的大视场高分辨率显微成像方法。
背景技术
在显微成像领域,更高的分辨率一直是追求的目标,但是在提高分辨率的同时存在一个关键性问题,那就是并没有随分辨率一起提高的显微镜的空间带宽积,换言之即传统显微镜存在分辨率与视场大小难以同时兼顾的矛盾。因为,传统显微镜使用低倍物镜进行成像时视场大但是分辨率低,而使用高倍物镜进行成像时分辨率提高了但是相应的视场就会缩得很小。目前,为了突破分辨率与视场大小难以同时兼顾的矛盾,常见的方法是采用常规显微镜系统配合高精度机械扫描和后期空域图像拼接方法将多个小视场高分辨率图像拼接融合生成一幅大视场高分辨率图像([1]孙西钊,王震,李攀,李建勋,王文合.适用于结核杆菌抗酸染色图像拼接的装置:中国,2013205777012[P].2013-09-17.)。但是由于引入了机械移动装置,所以系统成像时的稳定性和成像速度又成为一对难以调和的矛盾,提高扫描速度必将影响成像稳定性。所以,想要突破分辨率与视场大小难以同时兼顾的矛盾又不引入了机械移动装置,必须采用近年来提出的计算成像的方法,比如基于合成孔径的成像方法。
基于合成孔径成像原理的扫描成像方法最早是由Hoppe为了研究晶体结构所提出的,并通过研究晶体和非晶体的扫描透射电子衍射显微成像,验证了此方法的有效性。Rodenburg和Faulkner等结合相位恢复算法将此方法多次改进,目前这种成像方法已在可见光域、X射线、电子显微镜等不同波段得到了实验证实,并发展出若干种技术以提高成像质量以及分辨率,该技术显示了在大幅面成像和高分辨成像方面的巨大潜力。传统的合成孔径成像技术是通过移动一个全透的小孔(或待测样品本身)使入射平面波照射到待测样品的不同部位,即由小孔控制照明光束尺寸、几何形状及位置,并利用由此得到的一系列衍射强度图样重构出待测样品的振幅与位相信息([2]王雅丽,史祎诗,李拓,等.可见光域叠层成像中照明光束的关键参量研究[J].物理学报,2013,Vol.62,No.6.064206-1-9)。合成孔径成像术的关键在于:每次照射待测样品的一个“子孔径”也就是待测样品的某一部分时,都要和至少一个其他的“子孔径”发生交叠。这样就可建立一种重构算法,在分别重构每“子孔径”的复振幅时也要同时满足其他“子孔径”衍射分布的约束,使得最后的待测样品的整体复振幅信息是所有“子孔径”的共同解,从而由各个“子孔径”拼接合成一幅大视场高分辨率的待测样品的图像。合成孔径成像术可以说是一种稳健而简约的显微成像技术,但目前最突出的问题在于其成像效率与质量之间的矛盾,因为使用合成孔径成像技术拍摄一幅完整的超分辨率的图像往往需要花费十几或数十分钟,因此如何提高成像效率就成为了合成孔径成像技术必须克服的一个技术难题。
发明内容
本发明的目的在于提供一种基于哈达玛编码的大视场高分辨率显微成像方法,既能降低相机所需的曝光时间,提高系统的图像采集速度,又能抑制噪声,提高系统采集的图像质量,以解决显微系统分辨率与视场大小难以同时兼顾的矛盾。
实现本发明目的的技术解决方案为:一种基于最优哈达玛编码的大视场高分辨率显微成像方法,其特征在于步骤如下:
步骤一:利用计算机编程生成初始哈达玛编码矩阵;
步骤二:对初始哈达玛编码矩阵进行优化,得到最优哈达玛编码矩阵;
步骤三:采用LED阵列作为显微镜光源,显示最优哈达玛编码矩阵,并用相机拍摄待测样品在不同照明角度下的一系列低分辨率图像,再根据最优哈达玛矩阵通过解方程组求解出待测样品受到单个角度照明时的低分辨率图像;
步骤四:利用频域合成孔径技术将多幅大视场低分辨率图像合成一幅大视场高分辨率图像。
本发明与现有技术相比,其显著优点:(1)本发明利用LED阵列产生不同角度的入射光照射待测样品,再用相机拍摄多幅低分辨率图像,最后采用频域合成孔径技术将多幅大视场低分辨率图像合成一幅大视场高分辨率图像,无需任何机械扫描装置,系统成像稳定性高。(2)本方法利用哈达玛矩阵来编码照明光,产生多个角度的入射光同时照射待测样品,这样既能降低相机所需的曝光时间,提高系统的图像采集速度,又能抑制噪声,提高系统采集的图像质量,所以本方法可以非常稳定并且精确的重建出大视场高分辨率图像。
下面结合附图对本发明作进一步详细描述。
附图说明
图1为基于最优哈达玛编码的大视场高分辨率显微成像方法步骤流程示意图。
图2为基于合成孔径超分辨率重构算法来重构待测样品的大视场高分辨率的光强和相位分布步骤流程示意图。
图3为本发明生成的N幅最优化哈达玛编码LED图案中的一幅。
图4为测试对象为1951USAF分辨率板时相机拍摄的N幅图像中的一幅。
图5为测试对象为1951USAF分辨率板时解出的单照明角度下的N幅图像中的一幅。
图6为测试对象为1951USAF分辨率板时重构的大视场高分辨率显微图像;图6(a)是4倍物镜拍摄到的原始图像;图6(b1)是重构后的原始图像图6(a)中选框区域的高分辨率显微图像的放大图像;图6(b2)是图6(b1)中选框区域的放大图像;图6(b3)是图6(b2)中选框区域的放大图像;图6(c1)、图6(c2)、图6(c3)分别是重构前的原始低分辨率图像在图6(b1)、图6(b2)、图6(b3)图像对应位置处的插值放大图。
图7为测试对象为淋巴细胞切片时重构的大视场高分辨率显微图像;图7(a)是4倍物镜拍摄到的原始图像;图7(b1)是重构后的原始图像图7(a)中选框区域的高分辨率显微图像的放大图像;图7(b2)是图7(b1)中选框区域的放大图像;图7(b3)是图7(b2)中选框区域的放大图像;图7(c1)、图7(c2)、图7(c3)分别是重构前的原始低分辨率图像在图7(b1)、图7(b2)、图7(b3)图像对应位置处的插值放大图。
具体实施方式
结合图1,本发明首先利用计算机编程生成初始哈达玛编码矩阵。之后对初始哈达玛编码矩阵进行优化,得到最优哈达玛编码矩阵。然后采用LED阵列作为显微镜光源,显示最优哈达玛编码矩阵,并用相机拍摄待测样品在不同照明角度下的一系列低分辨率图像。再根据最优哈达玛矩阵通过解方程组求解出待测样品受到单个角度照明时的低分辨率图像。最后利用频域合成孔径技术将多幅大视场低分辨率图像合成一幅大视场高分辨率图像。具体实现步骤如下:
步骤一:利用计算机Matlab软件(或其他类似编程软件)编程生成初始哈达玛编码矩阵。生成公式为:
R=Rand(1,N)
Rsort=Sort(R)
W0(m,n)=H(1,n+m-1-kN),n=1,...,N,m=1,...,N
其中,R为包含有N个随机数的行向量,Rand(·)为产生随机数的函数,N表示LED的总个数,Rsort为R经过排序后由大到小排列的行向量,Sort(·)为由大到小排序的排序函数。H表示以向量Rsort中第C个元素的值为阈值,对R进行0-1阈值分割后得到的行向量,C表示每次同时点亮的LED的个数,H(1,n)表示行向量H中的第n个元素,取值范围为1-N。W0表示初始哈达玛编码矩阵,W0(m,n)为矩阵W0中第m行n列的元素,k为一个中间变量,取值为0或1。
步骤二:对初始哈达玛编码矩阵进行优化,得到最优哈达玛编码矩阵。优化的公式为:
Woptimal=Wn,Wn-Wn-1<ε
其中,MSE表示含噪声图像相比原始图像的均方误差,trace表示求矩阵的迹,Wi-1表示第i-1次优化后得到的哈达玛编码矩阵,i为优化迭代次数,i=1,2,3,...,表示矩阵Wi-1的转置矩阵。Γ为均方误差MSE相对哈达玛编码矩阵W的导数,表示求导,(·)-1表示求矩阵的逆矩阵。表示对原始矩阵Wi-1经过优化后得到的非整数矩阵,ones(N,1)表示包含N个1的列向量,ones(1,N)表示包含N个1的行向量。Wi表示第i次优化后得到的哈达玛编码矩阵,Round(·)表示进行四舍五入运算,Woptimal为最优哈达玛编码矩阵,矩阵中所有元素都是0或1,Wn为第n次优化后得到的哈达玛编码矩阵,ε表示优化过程中判断停止迭代的阈值。
步骤三:采用LED阵列作为显微镜光源,显示最优哈达玛编码矩阵,并用相机拍摄待测样品在不同照明角度下的一系列低分辨率图像。再根据最优哈达玛矩阵通过解方程组求解出待测样品受到单个角度照明时的低分辨率图像。求解公式为:
其中,x和y分别为图像中的x轴和y轴坐标,表示单个角度照明下的所有N个低分辨率图像中坐标为(x,y)的N个像素点灰度值所组成的列向量,为N行N列的最优哈达玛编码矩阵的逆矩阵,为相机拍摄到的多角度照明下的所有N个低分辨率图像中坐标为(x,y)的N个像素点灰度值所组成的列向量。
步骤四:利用频域合成孔径技术将多幅大视场低分辨率图像合成一幅大视场高分辨率图像。如图2所示,具体实现流程为:
第一步,从不同照明角度下的衍射图像Ip,q中选取正入射照明视角下的衍射光强图像I0,0(其中下标p,q表示LED阵列中第p行q列的LED所对应的照明角度,p,q=-7,-6,...,6,7,下标0,0表示LED阵列中第0行0列的LED所对应的照明角度),对它进行双线性插值计算,得到一幅高分辨率的光强图像(其中上标h表示高分辨率),开根号求出其振幅分布然后随意构造一幅高分辨率的相位分布比如假设利用振幅分布与相位分布即可得到待测样品的高分辨率的复振幅分布对其进行傅立叶变换则可得到待测样品的高分辨率的频谱分布Fh。
第二步,使用一个圆形孔径低通滤波器,在待测样品的高分辨率的频谱分布Fh中截取一个低分辨率的“子孔径”频谱(其中上标l表示低分辨率),该低分辨率频谱的中心位于空间频率(kp,kq)(其中空间频率kp,kq分别表示LED阵列中第p行q列的LED照明时入射光对应的x,y方向的空间频率),圆形孔径低通滤波器的半径由物镜的数值孔径决定。对“子孔径”频谱做逆傅立叶变换求出空间频率(kp,kq)下待测样品的低分辨率复振幅
第三步,用相机实际拍摄得到的待测样品的低分辨率振幅去替换计算得到的低分辨率复振幅中的振幅部分,从而得到新约束条件下的低分辨率复振幅(其中上标c表示振幅经过替换后的结果)。然后再对其做傅立叶变换得到新约束条件下的低分辨率频谱最后用低分辨率频谱去替换第二步中圆形孔径低通滤波器位置处的高分辨率的频谱分布Fh中的“子孔径”频谱分量。
第四步,选取另一个照明角度,重复第二和第三步,直到所有N个照明角度都进行了运算,或者说频域中N个“子孔径”都完成了扫描拼接。
第五步,把第二步到第四步迭代重复计算2到3次,以保证最后合成孔径后的超分辨率光强分布和相位分布距离真实值误差足够小。迭代完成后,将最终高分辨率的频谱分布Fh变换到空域中即得到了待测样品大视场高分辨率的光强分布Ih和相位分布
为了测试基于最优哈达玛编码的大视场高分辨率显微系统重建图像的分辨率,本发明分别选取1951USAF分辨能力测试板和淋巴细胞切片进行了合成孔径大视场高分辨率成像。实验中,使用的LED阵列包含15行15列共225个LED,并利用其产生225个不同角度的照明光,LED之间间距为3mm,发出的红光波长为632.8nm。另外,我们选取每次同时点亮的LED个数为33个。生成的最优哈达玛编码LED图案如图3所示,相机拍摄到的多角度照明下的低分辨率显微图像如图4所示,从中解出的单角度照明时的低分辨率显微图像如图5所示。如图6所示,图6(a)表示的是使用数值孔径0.1、放大倍率为4倍的物镜拍摄得到的原始图像。选取图6(a)方框中的图像直接进行放大,得到图6(c1)。从中再选择更小的区域直接进行放大,依次得到图6(c2)和图6(c3)。针对图6(a)中同样的方框中的图像进行合成孔径大视场高分辨率成像重构,可以重构出图6(b1)。从中再选择更小的区域直接进行放大,依次得到图6(b2)和图6(b3)。比较图6(b3)和图6(c3),明显可以看出前者图像的分辨率远高于后者图像的分辨率。由此可以看出本方法合成孔径大视场高分辨率成像的优越性。
1951 USAF分辨能力测试板的物理参数表
具体的结果看图6(b3),分辨能力测试板上第9组第3个元素依然能够清晰分辨,那么根据1951USAF分辨能力测试板的物理参数(见上表),可以估算出本方法合成孔径大视场高分辨率成像的最高分辨率高于1000um/(645*2)=0.775um。
此外,本发明对淋巴细胞切片进行了合成孔径大视场高分辨率成像,其中一组实验结果图如图7所示。图7(a)是4倍物镜拍摄到的原始图像;图7(b1)是重构后的原始图像图7(a)中选框区域的高分辨率显微图像的放大图像;图7(b2)是图7(b1)中选框区域的放大图像;图7(b3)是图7(b2)中选框区域的放大图像;图7(c1)、图7(c2)、图7(c3)分别是重构前的原始低分辨率图像在图7(b1)、图7(b2)、图7(b3)图像对应位置处的插值放大图。对比其中图7(b3)和图7(c3)可以明显看出本发明能够在低数值孔径大视场情况下实现大视场高分辨率成像。
Claims (4)
1.一种基于最优哈达玛编码的大视场高分辨率显微成像方法,其特征在于步骤如下:
步骤一:利用计算机编程生成初始哈达玛编码矩阵;
步骤二:对初始哈达玛编码矩阵进行优化,得到最优哈达玛编码矩阵;
步骤三:采用LED阵列作为显微镜光源,显示最优哈达玛编码矩阵,并用相机拍摄待测样品在不同照明角度下的一系列低分辨率图像,再根据最优哈达玛矩阵通过解方程组求解出待测样品受到单个角度照明时的低分辨率图像;
步骤四:利用频域合成孔径技术将多幅大视场低分辨率图像合成一幅大视场高分辨率图像;
在步骤四中,频域合成孔径技术实现的流程为:
第一步,从不同照明角度下的衍射图像Ip,q中选取正入射照明视角下的衍射光强图像I0,0,对衍射光强图像I0,0进行双线性插值计算,得到一幅高分辨率的光强图像开根号求出其振幅分布然后随意构造一幅高分辨率的相位分布利用振幅分布与相位分布即可得到待测样品的高分辨率的复振幅分布对其进行傅立叶变换则可得到待测样品的高分辨率的频谱分布Fh;
第二步,使用一个圆形孔径低通滤波器,在待测样品的高分辨率的频谱分布Fh中截取一个低分辨率的“子孔径”频谱该低分辨率频谱的中心位于空间频率(kp,kq),圆形孔径低通滤波器的半径由物镜的数值孔径决定,对“子孔径”频谱做逆傅立叶变换求出空间频率(kp,kq)下待测样品的低分辨率复振幅
第三步,用相机实际拍摄得到的待测样品的低分辨率振幅去替换计算得到的低分辨率复振幅中的振幅部分,从而得到新约束条件下的低分辨率复振幅然后再对其做傅立叶变换得到新约束条件下的低分辨率频谱最后用低分辨率频谱去替换第二步中圆形孔径低通滤波器位置处的高分辨率的频谱分布Fh中的“子孔径”频谱分量;
第四步,选取另一个照明角度,重复第二和第三步,直到所有N个照明角度都进行了运算,或者说频域中N个“子孔径”都完成了扫描拼接;
第五步,把第二步到第四步迭代重复计算2-3次,以保证最后合成孔径后的超分辨率光强分布和相位分布距离真实值误差足够小;迭代完成后,将最终高分辨率的频谱分布Fh变换到空域中即得到了待测样品大视场高分辨率的光强分布Ih和相位分布
2.根据权利要求1所述的基于最优哈达玛编码的大视场高分辨率显微成像方法,其特征在于步骤一中,生成初始哈达玛编码矩阵的生成公式为:
R=Rand(1,N)
Rsort=Sort(R)
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>H</mi>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>,</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mn>1</mn>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>R</mi>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>,</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo>&GreaterEqual;</mo>
<msub>
<mi>R</mi>
<mrow>
<mi>s</mi>
<mi>o</mi>
<mi>r</mi>
<mi>t</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>C</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>0</mn>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>R</mi>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>,</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo><</mo>
<msub>
<mi>R</mi>
<mrow>
<mi>s</mi>
<mi>o</mi>
<mi>r</mi>
<mi>t</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>C</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<mi>N</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
W0(m,n)=H(1,n+m-1-kN),n=1,...,N,m=1,...,N
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mn>0</mn>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>n</mi>
<mo>+</mo>
<mi>m</mi>
<mo>-</mo>
<mn>1</mn>
<mo>&le;</mo>
<mi>N</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>1</mn>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>n</mi>
<mo>+</mo>
<mi>m</mi>
<mo>-</mo>
<mn>1</mn>
<mo>></mo>
<mi>N</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>n</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<mi>N</mi>
<mo>,</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<mi>N</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
其中,R为包含有N个随机数的行向量,Rand(·)为产生随机数的函数,N表示LED的总个数,Rsort为R经过排序后由大到小排列的行向量,Sort(·)为由大到小排序的排序函数,H表示以向量Rsort中第C个元素的值为阈值,对R进行0-1阈值分割后得到的行向量,C表示每次同时点亮的LED的个数,H(1,n)表示行向量H中的第n个元素,取值范围为1-N,W0表示初始哈达玛编码矩阵,W0(m,n)为矩阵W0中第m行n列的元素,k为一个中间变量,取值为0或1。
3.根据权利要求1所述的基于最优哈达玛编码的大视场高分辨率显微成像方法,其特征在于步骤二中,最优哈达玛编码矩阵的优化公式为:
<mrow>
<msub>
<mi>MSE</mi>
<mrow>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mi>N</mi>
</mfrac>
<mi>t</mi>
<mi>r</mi>
<mi>a</mi>
<mi>c</mi>
<mi>e</mi>
<mo>&lsqb;</mo>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>W</mi>
<mrow>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mi>t</mi>
</msubsup>
<msub>
<mi>W</mi>
<mrow>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>&rsqb;</mo>
<mo>,</mo>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mn>3</mn>
<mo>,</mo>
<mn>...</mn>
</mrow>
<mrow>
<msub>
<mi>&Gamma;</mi>
<mrow>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mi>d</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>MSE</mi>
<mrow>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>d</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>W</mi>
<mrow>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<mn>2</mn>
<mi>N</mi>
</mfrac>
<msup>
<mrow>
<mo>(</mo>
<msubsup>
<mi>W</mi>
<mrow>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mi>t</mi>
</msubsup>
<msub>
<mi>W</mi>
<mrow>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<msubsup>
<mi>W</mi>
<mrow>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mi>t</mi>
</msubsup>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
</mrow>
<mrow>
<msubsup>
<mi>W</mi>
<mrow>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>u</mi>
<mi>b</mi>
</mrow>
</msubsup>
<mo>=</mo>
<msub>
<mi>W</mi>
<mrow>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>+</mo>
<mfrac>
<mn>1</mn>
<mn>5</mn>
</mfrac>
<mo>{</mo>
<mfrac>
<mn>1</mn>
<mi>N</mi>
</mfrac>
<mo>&lsqb;</mo>
<mi>C</mi>
<mo>-</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>W</mi>
<mrow>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>&Gamma;</mi>
<mrow>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>*</mo>
<mi>o</mi>
<mi>n</mi>
<mi>e</mi>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mi>N</mi>
<mo>,</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mo>*</mo>
<mi>o</mi>
<mi>n</mi>
<mi>e</mi>
<mi>s</mi>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>,</mo>
<mi>N</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>&Gamma;</mi>
<mrow>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mrow>
<msub>
<mi>W</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<mi>R</mi>
<mi>o</mi>
<mi>u</mi>
<mi>n</mi>
<mi>d</mi>
<mrow>
<mo>(</mo>
<msubsup>
<mi>W</mi>
<mrow>
<mi>i</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mrow>
<mi>u</mi>
<mi>b</mi>
</mrow>
</msubsup>
<mo>)</mo>
</mrow>
</mrow>
Woptimal=Wn,Wn-Wn-1<ε
其中,MSE表示含噪声图像相比原始图像的均方误差,trace表示求矩阵的迹,Wi-1表示第i-1次优化后得到的哈达玛编码矩阵,i为优化迭代次数,i=1,2,3,...,表示矩阵Wi-1的转置矩阵,Γ为均方误差MSE相对哈达玛编码矩阵W的导数,表示求导,(·)-1表示求矩阵的逆矩阵,表示对原始矩阵Wi-1经过优化后得到的非整数矩阵,C表示每次同时点亮的LED的个数,ones(N,1)表示包含N个1的列向量,ones(1,N)表示包含N个1的行向量,Wi表示第i次优化后得到的哈达玛编码矩阵,Round(·)表示进行四舍五入运算,Woptimal为最优哈达玛编码矩阵,矩阵中所有元素都是0或1,Wn为第n次优化后得到的哈达玛编码矩阵,ε表示优化过程中判断停止迭代的阈值。
4.根据权利要求1所述的基于最优哈达玛编码的大视场高分辨率显微成像方法,其特征在于步骤三中,根据哈达玛编码矩阵求解单个角度照明下的低分辨率图像的公式为:
<mrow>
<mover>
<mi>I</mi>
<mo>&RightArrow;</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msubsup>
<mi>W</mi>
<mrow>
<mi>o</mi>
<mi>p</mi>
<mi>t</mi>
<mi>i</mi>
<mi>m</mi>
<mi>a</mi>
<mi>l</mi>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msubsup>
<mover>
<mi>M</mi>
<mo>&RightArrow;</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
</mrow>
其中,x和y分别为图像中的x轴和y轴坐标,表示单个角度照明下的所有N个低分辨率图像中坐标为(x,y)的N个像素点灰度值所组成的列向量,为N行N列的最优哈达玛编码矩阵的逆矩阵,为相机拍摄到的多角度照明下的所有N个低分辨率图像中坐标为(x,y)的N个像素点灰度值所组成的列向量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510186030.0A CN104796609B (zh) | 2015-04-17 | 2015-04-17 | 基于最优哈达玛编码的大视场高分辨率显微成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510186030.0A CN104796609B (zh) | 2015-04-17 | 2015-04-17 | 基于最优哈达玛编码的大视场高分辨率显微成像方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104796609A CN104796609A (zh) | 2015-07-22 |
CN104796609B true CN104796609B (zh) | 2018-01-05 |
Family
ID=53561099
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510186030.0A Active CN104796609B (zh) | 2015-04-17 | 2015-04-17 | 基于最优哈达玛编码的大视场高分辨率显微成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104796609B (zh) |
Families Citing this family (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105225202B (zh) * | 2015-09-23 | 2018-06-29 | 清华大学深圳研究生院 | 一种基于fpm算法的图像超分辨率重建的方法 |
JP2018537709A (ja) * | 2015-11-11 | 2018-12-20 | スコピオ ラボズ リミテッド | 異なる照明条件下で画像を生成するためのコンピュータ顕微鏡及び方法 |
CN108604285B (zh) * | 2015-12-24 | 2021-09-21 | 得利捷Ip 科技有限公司 | 编码信息读取器 |
CN105911692B (zh) * | 2016-06-24 | 2018-06-12 | 南京理工大学 | 一种傅立叶叠层显微成像系统的最优系统参数选择方法 |
CN106204434B (zh) * | 2016-06-24 | 2019-04-16 | 南京理工大学 | 一种面向大视场高分辨率显微成像的图像迭代重构方法 |
CN106092319A (zh) * | 2016-08-26 | 2016-11-09 | 中国科学院寒区旱区环境与工程研究所 | 光谱直读型阿达玛变换光谱成像系统设计 |
CN108616686B (zh) * | 2016-12-12 | 2020-09-11 | 南京理工大学 | 一种全哈达玛矩阵编码成像方法 |
CN106707486B (zh) * | 2017-01-24 | 2019-07-26 | 清华大学 | 基于fpm的宽视场显微成像方法及系统 |
JP2020098984A (ja) * | 2018-12-17 | 2020-06-25 | キヤノン株式会社 | 画像符号化装置及び画像復号装置及びそれらの制御方法及びプログラム |
CN111031264B (zh) * | 2019-11-29 | 2021-10-08 | 南京理工大学 | 一种基于透射式红外孔径编码成像系统及其超分辨方法 |
CN111462026B (zh) * | 2020-03-02 | 2023-04-07 | 清华大学 | 基于编码曝光的合成图像相位恢复方法及装置 |
CN112580808B (zh) * | 2020-12-30 | 2024-04-26 | 杭州电子科技大学 | 一种基于改进型led矩阵无透镜成像方法 |
CN113709325B (zh) * | 2021-07-20 | 2023-09-15 | 辽宁大学 | 基于Hadamard频域变换矩阵阈值滤波的单像素成像方法 |
CN113917677B (zh) * | 2021-09-09 | 2023-05-05 | 北京纳析光电科技有限公司 | 一种三维超分辨率光片显微成像方法和显微镜 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1308230A (zh) * | 2001-03-08 | 2001-08-15 | 武汉大学 | 高分辨阿达玛变换显微图像分析仪 |
CN101915920A (zh) * | 2010-07-02 | 2010-12-15 | 北京航空航天大学 | 一种地球同步轨道合成孔径雷达卫星的高分辨率成像方法 |
CN102879782A (zh) * | 2012-09-25 | 2013-01-16 | 北京理工大学 | 基于分数阶傅里叶变换的压缩感知sar成像方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8358739B2 (en) * | 2010-09-03 | 2013-01-22 | The University Of North Carolina At Chapel Hill | Systems and methods for temporal multiplexing X-ray imaging |
US8941061B2 (en) * | 2012-03-28 | 2015-01-27 | Uchicago Argonne, Llc | Compressive passive millimeter wave imager |
-
2015
- 2015-04-17 CN CN201510186030.0A patent/CN104796609B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1308230A (zh) * | 2001-03-08 | 2001-08-15 | 武汉大学 | 高分辨阿达玛变换显微图像分析仪 |
CN101915920A (zh) * | 2010-07-02 | 2010-12-15 | 北京航空航天大学 | 一种地球同步轨道合成孔径雷达卫星的高分辨率成像方法 |
CN102879782A (zh) * | 2012-09-25 | 2013-01-16 | 北京理工大学 | 基于分数阶傅里叶变换的压缩感知sar成像方法 |
Non-Patent Citations (1)
Title |
---|
Coded multi-angular illumination for Fourier ptychography based on Hadamard codes;Jiasong Sun et al.;《International Conference on Optical and Photonic Engineering(icOPEN 2015)》;20150109;第9524卷;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN104796609A (zh) | 2015-07-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104796609B (zh) | 基于最优哈达玛编码的大视场高分辨率显微成像方法 | |
CN107065159B (zh) | 一种基于大照明数值孔径的大视场高分辨率显微成像装置及迭代重构方法 | |
CN107966801A (zh) | 一种基于环形照明的高速傅立叶叠层成像装置及重构方法 | |
US11368608B2 (en) | Compressed sensing based object imaging system and imaging method therefor | |
Luo et al. | Pixel super-resolution for lens-free holographic microscopy using deep learning neural networks | |
CN108761752A (zh) | 傅立叶重叠关联成像系统、设备和方法 | |
CN105911692B (zh) | 一种傅立叶叠层显微成像系统的最优系统参数选择方法 | |
CN111343376B (zh) | 一种基于透射式双缝孔径编码成像系统及其超分辨方法 | |
CN106204434B (zh) | 一种面向大视场高分辨率显微成像的图像迭代重构方法 | |
CN104459971A (zh) | 一种基于集成光子芯片的结构光照明显微成像系统 | |
CN109035168A (zh) | 一种基于翻折掩模构造的计算成像方法及装置 | |
CN110378981A (zh) | 基于神经网络的傅里叶叠层显微镜光瞳恢复方法 | |
CN109188881B (zh) | 一种大视场高分辨率太赫兹波数字全息成像方法和系统 | |
CN110942423A (zh) | 一种基于傅里叶叠层成像的远场超分辨率重建方法 | |
CN110455834A (zh) | 基于光强传输方程的x射线单次曝光成像装置及方法 | |
CN109211790A (zh) | 一种基于傅里叶功率谱探测的单像素相位成像方法 | |
Li et al. | Fast confocal microscopy imaging based on deep learning | |
CN108537862A (zh) | 一种自适应降噪的傅里叶衍射扫描显微镜成像方法 | |
CN110470667A (zh) | 基于压缩感知和随机多焦点扫描的三维体成像方法及装置 | |
CN111179371A (zh) | 一种扫描电子显微成像方法 | |
WO2020007761A1 (en) | Programmable multiple-point illuminator, confocal filter, confocal microscope and method to operate said confocal microscope | |
Sun et al. | Double‐flow convolutional neural network for rapid large field of view Fourier ptychographic reconstruction | |
CN105259668A (zh) | 基于黑色支撑的叠层成像技术 | |
Ma et al. | Scanning-position error-correction algorithm in dual-wavelength ptychographic microscopy | |
CN116841028A (zh) | 一种基于镜面干涉增强的多焦点宽场超分辨显微成像方法及装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
EXSB | Decision made by sipo to initiate substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |