CN107845077B - 一种校正dmd微镜元成像非对称形变的方法和存储介质 - Google Patents
一种校正dmd微镜元成像非对称形变的方法和存储介质 Download PDFInfo
- Publication number
- CN107845077B CN107845077B CN201711057691.9A CN201711057691A CN107845077B CN 107845077 B CN107845077 B CN 107845077B CN 201711057691 A CN201711057691 A CN 201711057691A CN 107845077 B CN107845077 B CN 107845077B
- Authority
- CN
- China
- Prior art keywords
- image
- dmd
- detector
- coded
- stripes
- 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
- 238000003384 imaging method Methods 0.000 title claims abstract description 39
- 238000000034 method Methods 0.000 title claims abstract description 38
- 238000012937 correction Methods 0.000 claims abstract description 12
- 230000003287 optical effect Effects 0.000 claims description 21
- 238000012545 processing Methods 0.000 claims description 20
- 238000005070 sampling Methods 0.000 claims description 13
- 238000011156 evaluation Methods 0.000 claims description 12
- 238000004590 computer program Methods 0.000 claims description 11
- 238000001914 filtration Methods 0.000 claims description 9
- 230000007306 turnover Effects 0.000 claims description 4
- 230000000903 blocking effect Effects 0.000 claims 2
- 238000012935 Averaging Methods 0.000 claims 1
- 238000010586 diagram Methods 0.000 description 8
- 238000004364 calculation method Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000000701 chemical imaging Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000009499 grossing Methods 0.000 description 2
- 230000003595 spectral effect Effects 0.000 description 2
- 239000000470 constituent Substances 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000003708 edge detection Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000012163 sequencing technique Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/80—Geometric correction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/12—Edge-based segmentation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/60—Analysis of geometric attributes
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/80—Analysis of captured images to determine intrinsic or extrinsic camera parameters, i.e. camera calibration
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T9/00—Image coding
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
- H04N19/00—Methods or arrangements for coding, decoding, compressing or decompressing digital video signals
- H04N19/85—Methods or arrangements for coding, decoding, compressing or decompressing digital video signals using pre-processing or post-processing specially adapted for video compression
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20024—Filtering details
- G06T2207/20032—Median filtering
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Multimedia (AREA)
- Geometry (AREA)
- Signal Processing (AREA)
- Image Processing (AREA)
Abstract
本发明公开了一种校正DMD微镜元成像非对称形变的方法和存储介质,所述方法利用DMD可加载标准回字形编码条纹这一特点,通过对加载后的编码图像中的回字形编码条纹的斜率、偏移数以及校正精度的计算,实现DMD镜元和探测器像元位置关系的精确标定,解决了传统通过调节探测器角度解决这一问题时导致成像***不能全视场清晰成像这一难题。
Description
技术领域
本发明涉及编码成像领域,尤其涉及一种校正DMD微镜元成像非对称形变的方法和存储介质。
背景技术
DMD(全称为Digital Micromirror Device,中文意思为“数字微镜器件”,组成单元为微镜元)以其灵活快捷的编码方式得到了广泛应用,并推动了编码孔径光谱成像技术的快速发展,应用领域包括编码孔径光谱成像仪、超分辨率编码孔径成像等。编码孔径光谱成像技术,具体实现如下:在一次像面采用特定的编码模板对目标进行编码,结合特殊的采样成像方式,获得满足景物重构的采样数据量,实现空间信息和光谱信息的高精度重构,相关知识可参考孙萍等人的论文“编码孔径和解码孔径的设计”以及谢存等人的论文“光电成像中的一种编码解码方法”等。
编码成像***中DMD处于倾斜的反射光路中,导致每个正方形的微镜元经后置光路在探测器上的成像产生非对称形变,所有微镜元形变叠加后,将会使编码微镜元和探测器像元原有的对应关系发生偏移无法精确配准,由此产生的采样误差会导致图谱重构的精度降低,甚至无法实现图谱重构。
对于上述问题,目前常用的解决方法是将探测器倾斜放置,与DMD处于镜像对称位置,这一方法通过简单的机械调节即可实现,但是探测器倾斜会引起***离焦,且离焦量随着视场的增大而增大,而成像***的焦深范围非常有限(通常为微米量级),因此,使用探测器倾斜补偿方法会导致清晰视场范围大幅降低,实际应用中会受到很大的限制。
发明内容
为此,需要提供一种校正DMD微镜元成像非对称形变的的技术方案,用以解决目前的利用探测器倾斜补偿来调整采样误差的方法会导致清晰视场范围大幅降低的问题。
为实现上述目的,发明人提供了一种校正DMD微镜元成像非对称形变的方法,所述方法包括以下步骤:
(1)在编码成像***的DMD上加载规则回字形编码图像,所述编码图像包括n组回字形编码条纹,所述回字形编码条纹中心位置相同;
(2)控制带成像***的探测器对编码图像进行采集,并对同一编码图像采样S次;
(3)对同一编码图像采集的S幅图像的每一像素点实施中值滤波运算,并进行二值化处理和prewitt算子处理;
(4)获取经过步骤(3)处理的编码图像中的某一组回字形编码条纹,以该回字形编码条纹的一边为基准,计算与该边相交的另一边的斜率并记录;具体包括:选取与基准边相交的另一边上的两个点(an,bn)和(cn,dn),相应的斜率为:
(5)对编码图像中每一组回字形编码条纹重复步骤(4),依次计算编码图像中各个回字形编码条纹中的斜率,建立与各个回字形编码条纹的编号顺序相对应的斜率数组,对应n组回字形编码条纹的斜率数组为(η1,η2,...,ηn);
(6)根据斜率数组依次计算出每行DMD编码元对应探测器像元的偏离数,对每一行像元的偏离数做取整运算;偏离数计算方式如下:其中,符号“[]”表示向下取整;对应n组回字形编码条纹的探测器像素点偏离数为(k1,k2,...,kn);
(7)以编码图像的中心点所在的水平轴线为基准,将各组回字形编码条纹的上半部分向靠近垂直轴线方向移动相对应的探测器像素点偏离数,并将各个回字形编码条纹的下半部分沿着上半部分移动的反方向移动相对应的探测器像素点偏离数,从而恢复出与DMD所加载的规则回字形编码图像形状相一致的图像;
(8)确定DMD编码元与探测器像素点的位置对应关系,并使用结构相似度评价函数对校正精度进行评价,结构相似性评价函数如下:
进一步地,所述方法还包括以下步骤:
将DMD每个正方形DMD微镜元沿其对角线翻转±α°,使得入射光信号得以反射到光路中两个不同方向,从而实现对光的选通。
进一步地,所述方法应用于调整处于倾斜的反射光路中的DMD每个正方形微镜元经后置光路在探测器上产生的非对称形变。
进一步地,所述S的数值为10以上。
发明人还提供了一种存储介质,所述存储介质存储有计算机程序,所述计算机程序被处理器执行时实现以下步骤:
(1)控制编码成像***的DMD加载规则回字形编码图像,所述编码图像包括n组回字形编码条纹,所述回字形编码条纹中心位置相同;
(2)控制带成像***的探测器对编码图像进行采集,并对同一编码图像采样S次;
(3)对同一编码图像采集的S幅图像的每一像素点实施中值滤波运算,并进行二值化处理和prewitt算子处理;
(4)获取经过步骤(3)处理的编码图像中的某一组回字形编码条纹,以该回字形编码条纹的一边为基准,计算与该边相交的另一边的斜率并记录;具体包括:选取与基准边相交的另一边上的两个点(an,bn)和(cn,dn),相应的斜率为:
(5)对编码图像中每一组回字形编码条纹重复步骤(4),依次计算编码图像中各个回字形编码条纹中的斜率,建立与各个回字形编码条纹的编号顺序相对应的斜率数组,对应n组回字形编码条纹的斜率数组为(η1,η2,...,ηn);
(6)根据斜率数组依次计算出每行DMD编码元对应探测器像元的偏离数,对每一行像元的偏离数做取整运算;偏离数计算方式如下:其中,符号“[]”表示向下取整;对应n组回字形编码条纹的探测器像素点偏离数为(k1,k2,...,kn);
(7)以编码图像的中心点所在的水平轴线为基准,将各组回字形编码条纹的上半部分向靠近垂直轴线方向移动相对应的探测器像素点偏离数,并将各个回字形编码条纹的下半部分沿着上半部分移动的反方向移动相对应的探测器像素点偏离数,从而恢复出与DMD所加载的规则回字形编码图像形状相一致的图像;
(8)确定DMD编码元与探测器像素点的位置对应关系,并使用结构相似度评价函数对校正精度进行评价,结构相似性评价函数如下:
进一步地,所述计算机程序被处理器执行时还实现以下步骤:
控制DMD每个正方形DMD微镜元沿其对角线翻转±α°,使得入射光信号得以反射到光路中不同方向,从而实现对光的选通。
进一步地,所述计算机程序被处理器执行时,应用于调整处于倾斜的反射光路中的DMD每个正方形微镜元经后置光路在探测器上产生的非对称形变。
进一步地,所述计算机程序被处理器执行时:对同一编码图像采样的次数S为10以上。
由于编码成像***属于静态采样***,***一旦集成固化,由倾斜光路引起的DMD编码模板非对称形变导致的采样误差将固定下来。针对这一特点,本发明提出基于规则回字形编码条纹对编码图像非对称形变进行定标校正,通过对加载后的编码图像中的回字形编码条纹的斜率、偏移数以及校正精度的计算,实现DMD镜元和探测器像元位置关系的精确标定,解决了传统调节探测器角度导致成像***不能全视场清晰成像这一难题。
附图说明
图1为本发明一实施例涉及的校正DMD微镜元成像非对称形变的方法的流程图;
图2为DMD微镜元与探测器像元理想对应关系的示意图;
图3为DMD工作时处于倾斜的反射光路中的示意图;
图4为DMD加载的规则回字形编码图像的示意图;
图5为探测器采集到的非对称形变编码条纹的示意图;
图6为对探测器重复采样的S幅图像取中值后的编码条纹图像;
图7为对取中值后的编码图像进行二值化处理后得到的图像;
图8为对二值化处理后的编码图像再使用prewitt算子进行处理后得到的图像;
图9为编码图像中各组回字形编码条纹的斜率计算的示意图;
图10为完成定标校正后的实际编码图像的示意图。
附图标记说明:
1、DMD;
2、成像***;
3、探测器;
4、控制计算机;
5、前置成像***;
6、DMD等效位置;
7、DMD实际位置;
8、中继成像***;
9、探测器。
具体实施方式
为详细说明技术方案的技术内容、构造特征、所实现目的及效果,以下结合具体实施例并配合附图详予说明。
如图2所示,为DMD微镜元与探测器像元理想对应关系的示意图。DMD(1)与探测器(3)分别与控制计算机(4)连接。控制计算机(4)发送编码信号至DMD(1),使其记载规则的回形条纹编码图像;以及发送图像采集信号至探测器(3),控制探测器(3)采集图像。DMD(1)与探测器(3)之间设置有成像***(2)。由图中可以看出,此时,DMD微镜元和探测器像元位置关系严格对应。
如图3所示,DMD正常工作状态下处于倾斜入射的反射光路中,由于DMD微镜元工作时沿对角线旋转±α°,因此DMD实际位置(7)和DMD等效位置(6)之间也有一个α°的夹角。此时,DMD的微镜元会在一个方向上受到压缩而发生非对称形变。即光路依次经过前置成像***(5)、DMD等效位置(6)、中继成像***(8)相较于理想光路将发生偏移。当探测器(9)接收到的所有镜元形变经过叠加后,将会使编码微镜元和探测器像元原有的对应关系发生偏移而无法精确配准。
请参阅图1,为本发明一实施例涉及的校正DMD微镜元成像非对称形变的方法的流程图。所述方法可用于对倾斜光路导致的DMD微镜元成像非对称形变进行校正,从而实现探测器采集的编码图像像素和DMD编码孔径编码元之间的精密匹配。所述方法包括以下步骤:
首先进入步骤S101在编码成像***的DMD上加载规则回字形编码图像,所述编码图像包括n组回字形编码条纹,所述回字形编码条纹中心位置相同。如图4所示,为DMD加载的规则回字形编码图像的示意图。
而后进入步骤S102控制带成像***的探测器对编码图像进行采集,并对同一编码图像采样S次。为了降低由环境振动等外界干扰带来的误差,在本实施方式中,所述S的数值为10以上。在探测器上采集到的编码图像如图5所示,对比图4和图5可知,倾斜的反射光路使得探测器采集到的条纹在一维方向上发生了形变。
而后进入步骤S103对同一编码图像采集的S幅图像的每一像素点实施中值滤波运算,并进行二值化处理和prewitt算子处理。
中值滤波法是一种非线性平滑技术,它将每一像素点的灰度值设置为该点某邻域窗口内的所有像素点灰度值的中值。中值滤波是基于排序统计理论的一种能有效抑制噪声的非线性信号处理技术,中值滤波的基本原理是把数字图像或数字序列中一点的值用该点的一个邻域中各点值的中值代替,让周围的像素值接近真实值,从而消除孤立的噪声点。方法是用某种结构的二维滑动模板,将板内像素按照像素值的大小进行排序,生成单调上升(或下降)的二维数据序列。二维中值滤波输出为g(x,y)=med{f(x-k,y-l),(k,l∈W)},其中,f(x,y),g(x,y)分别为原始图像和处理后图像。W为二维模板,通常为3*3,5*5区域,也可以是不同的的形状,如线状、圆形、十字形、圆环形等。对图5中探测器重复采样的S幅图像取中值后的编码图像如图6所示。
图像的二值化,就是将图像上的像素点的灰度值设置为0或255,也就是将整个图像呈现出明显的只有黑和白的视觉效果。一幅图像包括目标物体、背景还有噪声,要想从多值的数字图像中直接提取出目标物体,最常用的方法就是设定一个阈值T,用T将图像的数据分成两部分:大于T的像素群和小于T的像素群。这是研究灰度变换的最特殊的方法,称为图像的二值化(BINARIZATION)。通过二值化处理,可以有效简化图像数据,以便后续进行计算,对图6进行二值化处理后的图像如图7所示。
Prewitt算子是一种一阶微分算子的边缘检测,利用像素点上下、左右邻点的灰度差,在边缘处达到极值检测边缘,去掉部分伪边缘,对噪声具有平滑作用。其原理是在图像空间利用两个方向模板与图像进行邻域卷积来完成的,这两个方向模板一个检测水平边缘,一个检测垂直边缘。
对数字图像f(x,y),Prewitt算子的定义如下:
G(i)=|[f(i-1,j-1)+f(i-1,j)+f(i-1,j+1)]-[f(i+1,j-1)+f(i+1,j)+f(i+1,j+1)]|
G(j)=|[f(i-1,j+1)+f(i,j+1)+f(i+1,j+1)]-[f(i-1,j-1)+f(i,j-1)+f(i+1,j-1)]|
则P(i,j)=max[G(i),G(j)]或P(i,j)=G(i)+G(j);
经典Prewitt算子认为:凡灰度新值大于或等于阈值的像素点都是边缘点。即选择适当的阈值T,若P(i,j)≥T,则(i,j)为边缘点,P(i,j)为边缘图像。对图7进行Prewitt算子处理后的图像如图8所示。
而后进入步骤S104获取经过处理的编码图像中的某一组回字形编码条纹,以该回字形编码条纹的一边为基准,计算与该边相交的另一边的斜率并记录。具体包括:选取与基准边相交的另一边上的两个点(an,bn)和(cn,dn),相应的斜率为:
而后进入步骤S105对编码图像中每一组回字形编码条纹,依次计算编码图像中各个回字形编码条纹中的斜率,建立与各个回字形编码条纹的编号顺序相对应的斜率数组,对应n组回字形编码条纹的斜率数组为(η1,η2,...,ηn)。请参阅图9,为以每个编码条纹的横条纹为基准,计算与横条纹相交的竖条纹的斜率实施例。为了便于计算,计算斜率的竖条纹在各个编码条纹中的位置一致,例如图9各个编码条纹中计算斜率的边均为条纹的左侧边或者右侧边。当然,在另一些实施例中,也可以是以每个编码条纹的竖条纹为基准,计算与横条纹相交的横条纹的斜率。
而后进入步骤S106根据斜率数组依次计算出每行DMD编码元对应探测器像元的偏离数。对每一行像元的偏离数做取整运算;偏离数计算方式如下:其中,符号“[]”表示向下取整;对应n组回字形编码条纹的探测器像素点偏离数为(k1,k2,...,kn)。
而后进入步骤S107以编码图像的中心点所在的水平轴线为基准,将各组回字形编码条纹的上半部分向靠近垂直轴线方向移动相对应的探测器像素点偏离数,并将各个回字形编码条纹的下半部分沿着上半部分移动的反方向移动相对应的探测器像素点偏离数,从而恢复出与DMD所加载的规则回字形编码图像形状相一致的图像。例如图9中,调整前的编码图像整体向左倾斜(计算出的斜率为负数),以当前编码图像的中心点所在的水平轴线为基准,上半部分每行向右移动对应的探测器像元的偏离数,下半部分每行向左的探测器像元的偏离数。反之,如果在另一些实施例中,调整前的编码图像整体向右倾斜(计算出的斜率为正数),在进行调整过程中,以当前编码图像的中心点所在的水平轴线为基准,上半部分每行向左移动对应的探测器像元的偏离数,下半部分每行向右的探测器像元的偏离数。
而后进入步骤S108确定DMD编码元与探测器像素点的位置对应关系,并使用结构相似度评价函数对校正精度进行评价。结构相似性评价函数如下:
其中,Amn和Bmn分别为理想图像A和校正图像B中在坐标为(m,n)处的像素点的像素值;和分别为理想图像A和校正图像B上所有像素点的平均像素值。所述校正图像B为经过步骤S107调整后的图像,所述理想图像A为标准规则的回字形条纹图像,具体如图4所示。结构相似度是指校正后的回字形(校正图像B)和标准回字形条纹(理想图像A)的相似度,可以用于评价校正效果。DMD编码元与探测器像素点的位置对应关系是指步骤S107中调整过程中移动像素点的数量,即偏移数。
在某些实施例中,将DMD每个正方形微镜元沿其对角线翻转±α°,使得入射光信号得以反射到光路中两个不同方向,从而实现对光的选通。通过控制微镜元的翻转,可实现对光信号的编码。
在某些实施例中,所述方法应用于调整处于倾斜的反射光路中的DMD每个正方形微镜元经后置光路在探测器上产生的非对称形变,从而避免了传统使用倾斜探测器校正编码微镜元非对称形变时视场离焦模糊的缺陷。
发明人还提供了一种存储介质,所述存储介质存储有计算机程序,所述计算机程序被处理器执行时实现以下步骤:
(1)控制编码成像***的DMD加载规则回字形编码图像,所述编码图像包括n组回字形编码条纹,所述回字形编码条纹中心位置相同;
(2)控制带成像***的探测器对编码图像进行采集,并对同一编码图像采样S次;优选的,所述计算机程序被处理器执行时:对同一编码图像采样的次数S为10以上;
(3)对同一编码图像采集的S幅图像的每一像素点实施中值滤波运算,并进行二值化处理和prewitt算子处理;
(4)获取经过步骤(3)处理的编码图像中的某一组回字形编码条纹,以该回字形编码条纹的一边为基准,计算与该边相交的另一边的斜率并记录;具体包括:选取与基准边相交的另一边上的两个点(an,bn)和(cn,dn),相应的斜率为:
(5)对编码图像中每一组回字形编码条纹重复步骤(4),依次计算编码图像中各个回字形编码条纹中的斜率,建立与各个回字形编码条纹的编号顺序相对应的斜率数组,对应n组回字形编码条纹的斜率数组为(η1,η2,...,ηn);
(6)根据斜率数组依次计算出每行DMD编码元对应探测器像元的偏离数,对每一行像元的偏离数做取整运算;偏离数计算方式如下:其中,符号“[]”表示向下取整;对应n组回字形编码条纹的探测器像素点偏离数为(k1,k2,...,kn);
(7)以编码图像的中心点所在的水平轴线为基准,将各组回字形编码条纹的上半部分向靠近垂直轴线方向移动相对应的探测器像素点偏离数,并将各个回字形编码条纹的下半部分沿着上半部分移动的反方向移动相对应的探测器像素点偏离数,从而恢复出与DMD所加载的规则回字形编码图像形状相一致的图像;
(8)确定DMD编码元与探测器像素点的位置对应关系,并使用结构相似度评价函数对校正精度进行评价,结构相似性评价函数如下:
在某些实施例中,所述计算机程序被处理器执行时还实现以下步骤:控制DMD每个正方形DMD微镜元沿其对角线翻转±α°,使得入射光信号得以反射到光路中不同方向,从而实现对光的选通。
在某些实施例中,所述计算机程序被处理器执行时,应用于调整处于倾斜的反射光路中的DMD每个正方形微镜元经后置光路在探测器上产生的非对称形变,从而避免了传统使用倾斜探测器校正编码微镜元非对称形变时视场离焦模糊的缺陷。
本发明提供了一种校正DMD微镜元成像非对称形变的方法和存储介质,所述方法利用DMD可加载标准回字形编码条纹这一特点,通过对加载后的编码图像中的回字形编码条纹的斜率、偏移数以及校正精度的计算,实现DMD镜元和探测器像元位置关系的精确标定,解决了传统调节探测器角度导致成像***不能全视场清晰成像这一难题。
需要说明的是,尽管在本文中已经对上述各实施例进行了描述,但并非因此限制本发明的专利保护范围。因此,基于本发明的创新理念,对本文所述实施例进行的变更和修改,或利用本发明说明书及附图内容所作的等效结构或等效流程变换,直接或间接地将以上技术方案运用在其他相关的技术领域,均包括在本发明的专利保护范围之内。
Claims (6)
1.一种校正DMD微镜元成像非对称形变的方法,其特征在于,所述方法包括以下步骤:
1)在编码成像***的DMD上加载规则回字形编码图像,所述编码图像包括n组回字形编码条纹,所述回字形编码条纹中心位置相同;
2)控制带成像***的探测器对所述编码图像进行采集,并对同一所述编码图像采样S次,所述S的数值为10以上;
3)对同一所述编码图像采集的S幅图像的每一像素点实施中值滤波运算,并进行二值化处理和prewitt算子处理;
4)获取经过步骤3)处理的编码图像中的某一组回字形编码条纹,以该回字形编码条纹的一边为基准,计算与该回字形编码条纹的一边相交的另一边的斜率并记录;具体包括:选取与基准边相交的另一边上的两个点(an,bn)和(cn,dn),相应的斜率为:
5)对编码图像中每一组回字形编码条纹重复步骤4),依次计算所述编码图像中各个回字形编码条纹中的斜率,建立与各个回字形编码条纹的编号顺序相对应的斜率数组,对应n组回字形编码条纹的斜率数组为(η1,η2,...,ηn);
6)根据斜率数组依次计算出每行DMD编码元对应探测器像素点的偏离数,对每一行像元的偏离数做取整运算;偏离数计算方式如下:其中,符号[]表示向下取整;对应n组回字形编码条纹的探测器像素点偏离数为(k1,k2,...,kn);
7)以所述编码图像的中心点所在的水平轴线为基准,将各组回字形编码条纹的上半部分向靠近垂直轴线方向移动相对应的探测器像素点偏离数,并将各个回字形编码条纹的下半部分沿着上半部分移动的反方向移动相对应的探测器像素点偏离数,从而恢复出与DMD所加载的规则回字形编码图像形状相一致的图像;
8)确定DMD编码元与探测器像素点的位置对应关系,并使用结构相似度评价函数对校正精度进行评价,结构相似性评价函数如下:
其中:所述prewitt算子处理采用如下方式进行:
G(i)=|[f(i-1,j-1)+f(i-1,j)+f(i-1,j+1)]-[f(i+1,j-1)+f(i+1,j)+f(i+1,j+1)]|
G(j)=|[f(i-1,j+1)+f(i,j+1)+f(i+1,j+1)]-[f(i-1,j-1)+f(i,j-1)+f(i+1,j-1)]|
则P(i,j)=max[G(i),G(j)]或P(i,j)=G(i)+G(j);
选择适当的阈值T,若P(i,j)≥T,则(i,j)为边缘点,P(i,j)为边缘图像。
2.根据权利要求1所述的校正DMD微镜元成像非对称形变的方法,其特征在于,所述方法还包括以下步骤:
将DMD每个正方形微镜元沿其对角线翻转±α°,使得入射光信号得以反射到光路中两个不同方向,从而实现光束的通过与阻断。
3.根据权利要求1所述的校正DMD微镜元成像非对称形变的方法,其特征在于,所述方法应用于调整处于倾斜的反射光路中的DMD每个正方形微镜元经后置光路在探测器上产生的非对称形变。
4.存储一种校正DMD微镜元成像非对称形变的方法的程序的存储介质,其特征在于,所述存储介质存储有计算机程序,所述计算机程序被处理器执行时实现以下步骤:
1)控制编码成像***的DMD加载规则回字形编码图像,所述编码图像包括n组回字形编码条纹,所述回字形编码条纹中心位置相同;
2)控制带成像***的探测器对所述编码图像进行采集,并对同一所述编码图像采样S次,所述S的数值为10以上;
3)对同一所述编码图像采集的S幅图像的每一像素点实施中值滤波运算,并进行二值化处理和prewitt算子处理;
4)获取经过步骤3)处理的编码图像中的某一组回字形编码条纹,以该回字形编码条纹的一边为基准,计算与该回字形编码条纹的一边相交的另一边的斜率并记录;具体包括:选取与基准边相交的另一边上的两个点(an,bn)和(cn,dn),相应的斜率为:
5)对编码图像中每一组回字形编码条纹重复步骤4),依次计算所述编码图像中各个回字形编码条纹中的斜率,建立与各个回字形编码条纹的编号顺序相对应的斜率数组,对应n组回字形编码条纹的斜率数组为(η1,η2,...,ηn);
6)根据斜率数组依次计算出每行DMD编码元对应探测器像素点的偏离数,对每一行像元的偏离数做取整运算;偏离数计算方式如下:其中,符号[]表示向下取整;对应n组回字形编码条纹的探测器像素点偏离数为(k1,k2,...,kn);
7)以所述编码图像的中心点所在的水平轴线为基准,将各组回字形编码条纹的上半部分向靠近垂直轴线方向移动相对应的探测器像素点偏离数,并将各个回字形编码条纹的下半部分沿着上半部分移动的反方向移动相对应的探测器像素点偏离数,从而恢复出与DMD所加载的规则回字形编码图像形状相一致的图像;
8)确定DMD编码元与探测器像素点的位置对应关系,并使用结构相似度评价函数对校正精度进行评价,结构相似性评价函数如下:
其中:所述prewitt算子处理采用如下方式进行:
G(i)=|[f(i-1,j-1)+f(i-1,j)+f(i-1,j+1)]-[f(i+1,j-1)+f(i+1,j)+f(i+1,j+1)]|
G(j)=|[f(i-1,j+1)+f(i,j+1)+f(i+1,j+1)]-[f(i-1,j-1)+f(i,j-1)+f(i+1,j-1)]|
则P(i,j)=max[G(i),G(j)]或P(i,j)=G(i)+G(j);
选择适当的阈值T,若P(i,j)≥T,则(i,j)为边缘点,P(i,j)为边缘图像。
5.根据权利要求4所述的存储介质,其特征在于,所述计算机程序被处理器执行时还实现以下步骤:
控制DMD每个正方形DMD微镜元沿其对角线翻转±α°,使得入射光信号得以反射到光路中两个不同方向,从而实现光束的通过与阻断。
6.根据权利要求4所述的存储介质,其特征在于,所述计算机程序被处理器执行时,应用于调整处于倾斜的反射光路中的DMD每个正方形微镜元经后置光路在探测器上产生的非对称形变。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711057691.9A CN107845077B (zh) | 2017-11-01 | 2017-11-01 | 一种校正dmd微镜元成像非对称形变的方法和存储介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711057691.9A CN107845077B (zh) | 2017-11-01 | 2017-11-01 | 一种校正dmd微镜元成像非对称形变的方法和存储介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107845077A CN107845077A (zh) | 2018-03-27 |
CN107845077B true CN107845077B (zh) | 2021-11-16 |
Family
ID=61681141
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711057691.9A Active CN107845077B (zh) | 2017-11-01 | 2017-11-01 | 一种校正dmd微镜元成像非对称形变的方法和存储介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107845077B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110146032B (zh) * | 2019-05-23 | 2020-09-25 | 复旦大学 | 基于光场分布的合成孔径相机标定方法 |
CN112712478B (zh) * | 2020-12-22 | 2022-11-08 | 安徽地势坤光电科技有限公司 | 一种数字微镜工作角度误差的修正方法及设备 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1174825A1 (en) * | 2000-07-12 | 2002-01-23 | Deutsche Thomson-Brandt Gmbh | Method for processing video pictures and apparatus for processing video pictures |
CN1816456A (zh) * | 2003-07-07 | 2006-08-09 | 联邦科学和工业研究机构 | 形成反射图样的方法 |
CN101656857A (zh) * | 2008-08-19 | 2010-02-24 | 精工爱普生株式会社 | 投影型显示装置和显示方法 |
CN102799080A (zh) * | 2012-09-03 | 2012-11-28 | 中国科学院光电技术研究所 | 一种掩模光栅对准标记倾斜成像的消除方法 |
CN104730820A (zh) * | 2013-12-19 | 2015-06-24 | 卡西欧计算机株式会社 | 投影装置及几何学校正调整方法 |
CN106643562A (zh) * | 2016-10-27 | 2017-05-10 | 天津大学 | 基于时域空域混合编码的结构光条纹投射方法 |
-
2017
- 2017-11-01 CN CN201711057691.9A patent/CN107845077B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1174825A1 (en) * | 2000-07-12 | 2002-01-23 | Deutsche Thomson-Brandt Gmbh | Method for processing video pictures and apparatus for processing video pictures |
CN1816456A (zh) * | 2003-07-07 | 2006-08-09 | 联邦科学和工业研究机构 | 形成反射图样的方法 |
CN101656857A (zh) * | 2008-08-19 | 2010-02-24 | 精工爱普生株式会社 | 投影型显示装置和显示方法 |
CN102799080A (zh) * | 2012-09-03 | 2012-11-28 | 中国科学院光电技术研究所 | 一种掩模光栅对准标记倾斜成像的消除方法 |
CN104730820A (zh) * | 2013-12-19 | 2015-06-24 | 卡西欧计算机株式会社 | 投影装置及几何学校正调整方法 |
CN106643562A (zh) * | 2016-10-27 | 2017-05-10 | 天津大学 | 基于时域空域混合编码的结构光条纹投射方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107845077A (zh) | 2018-03-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9451132B2 (en) | System for capturing a document in an image signal | |
US9420276B2 (en) | Calibration of light-field camera geometry via robust fitting | |
CN107607040B (zh) | 一种适用于强反射表面的三维扫描测量装置及方法 | |
Joshi et al. | PSF estimation using sharp edge prediction | |
CN108776980B (zh) | 一种面向微透镜光场相机的标定方法 | |
Placht et al. | Rochade: Robust checkerboard advanced detection for camera calibration | |
US10353191B2 (en) | Circular scanning technique for large area inspection | |
CN113205593B (zh) | 一种基于点云自适应修复的高反光表面结构光场三维重建方法 | |
CN106504290B (zh) | 一种高精度的摄像机动态标定方法 | |
US20040056966A1 (en) | Method and apparatus for image mosaicing | |
Hasinoff et al. | Confocal stereo | |
TWI480578B (zh) | 偵測廣角鏡頭之光學中心的方法以及光學中心偵測裝置 | |
CN107248139B (zh) | 基于显著视觉和dmd阵列分区控制的压缩感知成像方法 | |
Dufour et al. | Integrated digital image correlation for the evaluation and correction of optical distortions | |
WO2022126870A1 (en) | Three-dimensional imaging method and method based on light field camera and three-dimensional imaging measuring production line | |
CN105894521A (zh) | 基于高斯拟合的亚像素边缘检测方法 | |
CN107845077B (zh) | 一种校正dmd微镜元成像非对称形变的方法和存储介质 | |
CN112801870B (zh) | 一种基于网格优化的图像拼接方法,拼接***及可读存储介质 | |
CN104200456B (zh) | 一种用于线结构光三维测量的解码方法 | |
CN113447241B (zh) | 分割投影成像***的快速定标方法及装置 | |
CN107220945A (zh) | 多重退化的极模糊图像的复原方法 | |
Prakash et al. | Camera calibration using adaptive segmentation and ellipse fitting for localizing control points | |
Rosenbush et al. | Super-resolution enhancement of flash LADAR range data | |
CN114998571B (zh) | 基于固定尺寸标志物的图像处理及颜色检测方法 | |
Furukawa et al. | Simultaneous shape registration and active stereo shape reconstruction using modified bundle adjustment |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |