本发明属于数字全息图重建技术领域,涉及一种针对于大视野成像的大尺寸全息图的滑动窗分块重建方法。
背景技术:
数字全息是近年来新兴的一种精度可达微米级别的成像方法,由ccd等光学元器件直接记录获取样本图像,具有全视场、非接触、无损伤、实时性、定量化的优点。因为其无需对样本颜色即可清晰成像,故特别适合于活体生物样品的定量三维重建和快速跟踪,得以在生物医学应用领域尤其是细胞培养观测中发展迅速。
全息图成像视野直接对应于光学元器件ccd尺寸,所以其得天独厚具有大视野成像的优点。在能够对微米级别物体准确成像的前提下,全息图视野下所记录的物体信息远丰富于常规光学显微镜,在对大样本微观物体群的观察、记录方面颇具优势。
但是大尺寸全息图的重建无法回避的一个问题在于必须对图像进行分块后再逐一完成各图像块的卷积重建。经实验验证过大的分块或直接对全息图整图重建总会在重建时引入过多周围无关的光场信息,进而形成难以去除的横纵条带和背景斑块。此外考虑到卷积重建法需要大量正、逆傅里叶变换,而快速傅里叶变换对于长度为2的整数次幂的信号运算效率最高,故综合考量后认为单块尺寸选取256像素或512像素重建效果最佳,过大或过小均会造成不良效果。
伴随图像分割,一些物体势必处于分割的边缘位置,衍射环被人为割裂分布在不同图像子块上,不完整的衍射环重建后导致物体信息被损伤,难以反映物体的真实形貌,可能导致实验人员观察后误判、漏判,甚至造成实验失败。除此之外,受物体分布的影响,不同块之间重建后整体亮度差异往往十分明显,分块痕迹显著,不同图像子块内物体的可比对性受到影响,对于实验人员的观察和分析工作形成巨大干扰与挑战。因此,亟需一种针对于大视野成像的大尺寸全息图的滑动窗分块重建方法,用以提高重建图像质量和物体信息保真度。
技术实现要素:
有鉴于此,本发明的目的在于提供一种基于自适应滑动窗的分块处理策略,通过互相交叠的分块,确保全息图观察区域内位于分块边界上的物体总会在窗体移动后落入下一图像块内部,保证了其至少在一次重建过程中物体信息可以得到完整保留,且随滑动步长的细化,重建后物体信息得以完整保留的次数增多,达到提高重建图像质量和物体信息保真度的目的。
为达到上述目的,本发明提供如下技术方案:
一种用于提高大视野全息成像质量的自适应滑动窗重建方法,首先对原始全息图进行预处理和裁剪;然后使用常规分块方法对整图进行单次无交叠分块并完成各块重建,记录各块灰度均值;再计算各块灰度均值集合的离出率并据此设置窗体移动步长;最后按照设置步长滑动窗体加以分块重建,并将各块交叠对应像素点灰度赋为叠加平均值,使变化明显的背景干扰项得以抵消,使相对固定的物体信息得以保留;
该方法具体包括以下步骤:
s1:读入待重建全息图,进行预处理;
s2:设置分块尺寸,对图像进行裁剪;
s3:初步分块、重建,记录各块平均亮度,根据数据离出率自适应设置窗体滑动步长;
s4:按照s1中设置分块尺寸与s3中计算所得自适应步长移动窗体,对图像有效区域进行窗体遍历,并对各窗内图像块进行重建;
s5:将s4中处理所得的两个矩阵相同坐标处数据对应相除,得到最终经过自适应滑动窗分块处理后叠加平均的重建图像rec。
进一步,所述步骤s1具体包括:
将其从rgb彩色空间映射为灰度图像,图像rgb空间向灰度转换公式如下式所示,其中i、j为行列坐标,r、g、b分别为三原色中红、绿、蓝通道信息:
gray(i,j)=0.229×r(i,j)+0.587×g(i,j)+0.114×b(i,j)。
进一步,所述步骤s2具体包括:
s21:设图像块分块边长为l,整幅图像行列值为m和n,单位为像素(pixel);对应整除取其行、列所能排布图像块的个数row和col:
row=m|l
col=n|l
s22:取整幅图像中间row×l行与col×l列区域为有效重建区域,完成图像裁剪。
进一步,所述步骤s3具体包括:
s31:在有效区域内进行图像块的无交叠初步排布,即传统分块方法,对各图像块逐一重建,并记录各重建图像块平均亮度值,记录在row×col的矩阵中;
s32:对矩阵中所记录的平均亮度数据进行利群(outlier)分析,采用目前国际通用标准:以低于箱形图(boxplot)下箱体的1.5倍个四分位间距(inter-quartilerange,iqr),或是高于箱形图上箱体的1.5倍个四分位间距作为离群值界限;
s33:统计离群值个数numout,并计算各块平均灰度值的离群率ratioout:
s34:根据利群率设置自适应的窗体滑动步长step,其中[·]表示四舍五入取整:
该滑动步长的设置理由在于,图像各块平均灰度值离出率越高说明图像块重建后块间差异越大,需要更为细致的滑动窗进行处理。
进一步,所述步骤s4具体包括:
某一全息图像块,当其进行p次行方向滑动、q次列方向滑动后,以全息整图有效区域最左上角像素点为参考点,记为(1,1),则该图像块的行:l;列:l,对应覆盖在全息图有效区域中的范围为,行:(p-1)×l+1至p×l;列:(q-1)×l+1至q×l;
设置两个与全息图有效区域等大的空矩阵matrix1和matrix2,分别用于叠加各全息图块重建图recp,q对应位置处像素灰度值和叠加次数;
其中图像:x=(p-1)×l+i,y=(q-1)×l+j。
进一步,所述步骤s5中,所述重建图像rec为:rec(x,y)=marix1(x,y)/matrix2(x,y)。
进一步,对于所述步骤s1和s4中的全息图重建,均采用卷积重建法,具体流程如下:
卷积重建法是基于线性系统理论和瑞利-索默菲积分公式而构建的,能够准确反映光的空间衍射分布,由瑞利-索默菲衍射积分公式,全息图中衍射信息再现后对应光波:
其中,
r(x,y)=exp[-ima·k(xsinα+ysinβ)]
其中,k=2π/λ,α、β分别为平面参考光波与空间yoz平面和xoz平面的夹角;
而根据线性系统理论,再现像复振幅表示为如下卷积形式:
其中,g为自由空间脉冲响应;
再由卷积理论,使用傅里叶变换计算得到再现像,重建图:
u=ft-1[ft(h·c)·ft(g)]。
本发明的有益效果在于:本发明使大尺寸全息图重建结果的保真度得到了增强,所涵盖的物体信息得以完善,整图视觉感受更加趋于真实和自然,大大方便了实验人员的观察记录。
附图说明
为了使本发明的目的、技术方案和有益效果更加清楚,本发明提供如下附图进行说明:
图1为本发明所述的滑动窗重建方法的流程图;
具体实施方式
下面将结合附图,对本发明的优选实施例进行详细的描述。
图1为用于提高大视野全息成像质量的自适应滑动窗重建方法的流程图,如图1所示,一种用于提高大视野全息成像质量的自适应滑动窗重建方法,具体包括以下步骤:
s1:读入待重建全息图,进行预处理;
将其从rgb彩色空间映射为灰度图像,图像rgb空间向灰度转换公式如下式所示,其中i、j为行列坐标,r、g、b分别为三原色中红、绿、蓝通道信息:
gray(i,j)=0.229×r(i,j)+0.587×g(i,j)+0.114×b(i,j)。
s2:设置分块尺寸,对图像进行裁剪;
s21:设图像块分块边长为l,整幅图像行列值为m和n,单位为像素(pixel);对应整除取其行、列所能排布图像块的个数row和col:
row=m|l
col=n|l
s22:取整幅图像中间row×l行与col×l列区域为有效重建区域,完成图像裁剪。
s3:初步分块、重建,记录各块平均亮度,根据数据离出率自适应设置窗体滑动步长;
s31:在有效区域内进行图像块的无交叠初步排布,即传统分块方法,对各图像块逐一重建,并记录各重建图像块平均亮度值,记录在row×col的矩阵中;
s32:对矩阵中所记录的平均亮度数据进行利群(outlier)分析,采用目前国际通用标准:以低于箱形图(boxplot)下箱体的1.5倍个四分位间距(inter-quartilerange,iqr),或是高于箱形图上箱体的1.5倍个四分位间距作为离群值界限;
s33:统计离群值个数numout,并计算各块平均灰度值的离群率ratioout:
s34:根据利群率设置自适应的窗体滑动步长step,其中[·]表示四舍五入取整:
该滑动步长的设置理由在于,图像各块平均灰度值离出率越高说明图像块重建后块间差异越大,需要更为细致的滑动窗进行处理。
s4:按照s1中设置分块尺寸与s3中计算所得自适应步长移动窗体,对图像有效区域进行窗体遍历,并对各窗内图像块进行重建;
某一全息图像块,当其进行p次行方向滑动、q次列方向滑动后,以全息整图有效区域最左上角像素点为参考点,记为(1,1),则该图像块的行:l;列:l,对应覆盖在全息图有效区域中的范围为,行:(p-1)×l+1至p×l;列:(q-1)×l+1至q×l;
设置两个与全息图有效区域等大的空矩阵matrix1和matrix2,分别用于叠加各全息图块重建图recp,q对应位置处像素灰度值和叠加次数;
其中图像:x=(p-1)×l+i,y=(q-1)×l+j。
s5:将s4中处理所得的两个矩阵相同坐标处数据对应相除,得到最终经过自适应滑动窗分块处理后叠加平均的重建图像rec为:rec(x,y)=marix1(x,y)/matrix2(x,y)。
最后说明的是,以上优选实施例仅用以说明本发明的技术方案而非限制,尽管通过上述优选实施例已经对本发明进行了详细的描述,但本领域技术人员应当理解,可以在形式上和细节上对其作出各种各样的改变,而不偏离本发明权利要求书所限定的范围。