CN114820761A - 基于图像显微扫描平台xy方向夹角测量与运动补偿方法 - Google Patents
基于图像显微扫描平台xy方向夹角测量与运动补偿方法 Download PDFInfo
- Publication number
- CN114820761A CN114820761A CN202210493682.9A CN202210493682A CN114820761A CN 114820761 A CN114820761 A CN 114820761A CN 202210493682 A CN202210493682 A CN 202210493682A CN 114820761 A CN114820761 A CN 114820761A
- Authority
- CN
- China
- Prior art keywords
- image
- included angle
- scanning platform
- calibration line
- scanning
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 46
- 230000000007 visual effect Effects 0.000 claims abstract description 40
- 238000005259 measurement Methods 0.000 claims abstract description 7
- 230000002159 abnormal effect Effects 0.000 claims description 13
- 230000004927 fusion Effects 0.000 claims description 6
- 238000007526 fusion splicing Methods 0.000 claims description 5
- 150000001875 compounds Chemical class 0.000 claims 1
- 238000009434 installation Methods 0.000 abstract description 11
- 238000001514 detection method Methods 0.000 abstract description 5
- 230000001575 pathological effect Effects 0.000 abstract description 5
- 238000012545 processing Methods 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 6
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 238000007689 inspection Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 208000003464 asthenopia Diseases 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 238000011900 installation process Methods 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012549 training Methods 0.000 description 1
Images
Classifications
-
- 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
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/26—Measuring arrangements characterised by the use of optical techniques for measuring angles or tapers; for testing the alignment of axes
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T3/00—Geometric image transformations in the plane of the image
- G06T3/40—Scaling of whole images or parts thereof, e.g. expanding or contracting
- G06T3/4038—Image mosaicing, e.g. composing plane images from plane sub-images
-
- 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/13—Edge detection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/70—Determining position or orientation of objects or cameras
- G06T7/77—Determining position or orientation of objects or cameras using statistical methods
-
- 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/10—Image acquisition modality
- G06T2207/10056—Microscopic image
- G06T2207/10061—Microscopic image from scanning electron microscope
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Geometry (AREA)
- Life Sciences & Earth Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Evolutionary Biology (AREA)
- Probability & Statistics with Applications (AREA)
- Image Analysis (AREA)
Abstract
本发明公开了一种基于图像显微扫描平台XY方向夹角测量与运动补偿方法,能够帮助仪器装配人员方便、直观地判定测量扫描平台运动方向安装的夹角,根据所测夹角重新调整平台结构装配,使扫描平台沿着X轴和Y轴两个方向运动时保持运行方向垂直或者根据运动方向的夹角对平台沿着扫描轨迹运动采图做出对应的运动补偿,解决两运动方向安装垂直度不达标导致扫描平台无法准确完成所采集视野图像的高效拼接的问题,有效提高自动扫描过程图像处理效率,降低后期图像拼接的流程实现难度,降低技术装配人员的工作难度与强度的同时保证了调试的效率和精度,对维持自动显微扫描仪的稳定采图有重要意义,在医疗病理类玻片显微视觉自动检测行业有着广泛应用前景。
Description
技术领域
本发明涉及医疗信息采集检测和临床医学技术领域,尤其涉及一种基于图像显微扫描平台XY方向夹角测量与运动补偿方法。
背景技术
伴随着检验水平与检验需求的不断上涨,自动显微扫描平台相比传统显微镜展现出了更为便捷准确高效的特性,得到了各方面的广泛应用,而传统的人工镜检需要病理专家或经过专业培训的医师进行操作,已经不再适用于待检样本数量大、镜检结果得出时间紧的情况,且人工长时间镜检操作易导致视觉疲劳等情况影响镜检结果,加之人工镜检的重复聚焦操作导致耗时耗力、效率低下。在数字图像技术和计算机技术的协同下,自动显微扫描仪得到充分发展并应用,其能够完成自动聚焦,使镜检更高效准确,大大降低了使用成本,缩短了检验时间。
自动显微扫描仪所采集到的图像质量主要受到其扫描平台安装精度的影响,安装精度的高低对扫描仪的工作性能、效率、准确性等起到关键作用,加之成像原理要求扫描平台的运动X、Y方向有较高的垂直度要求,安装垂直度是扫描平台内部结构装配问题且对安装调试要求较高,如图1所示。
然而,目前测量其运动方向安装垂直度的方法主要依靠激光粗略测量,调整安装精度较差。如图2所示的情况,图中粗实线为标定玻片上的标定线,粗虚线为标定线在相机中实际采集所得。对实际采样的带有特制标定线的标定玻片,如图3所示,进行采集图像获得若干个视野图像。常用的图像拼接方法有几何拼接、图像识别融合拼接。其中图像识别融合拼接如图4所示,将采集到的所有视野进行图像内容识别拼接,并以此得到完整的图像内容,但是仍需要进行对拼接后的图像进行裁剪、旋转等操作才能够能到完整且正确的视野图像。
如此,现有测量方法很难保障扫描平台运动过程中在X、Y方向能够满足垂直度要求。因而在实际扫描过程中,由于安装精度等问题,往往所采集到的视野图像在拼接过程中都会出现图像内容不符、图像区域过大等问题,从而导致拼接图像失败或者拼接后的图像难以满足后续工作所需。
发明内容
本发明针对上述技术问题,提出一种基于图像显微扫描平台XY方向夹角测量与运动补偿方法,基于特制标定玻片视野标定线图像识别拼接后再拟合的方法来测量当前扫描平台运动方向的夹角,在理想状态下扫描平台沿X方向与其沿Y方向步进运动夹角为90°,所拍摄的不同视野图像经过图像拼接之后的图像边缘为平齐的,且拼接后图像内容与实际玻片标本一致。
为了实现上述目的,本发明提供如下技术方案:
本发明提供了一种基于图像显微扫描平台XY方向夹角测量与运动补偿方法,包括以下步骤:
步骤一:将带有特制标定线的标定玻片放置在显微扫描平台上,通过自动聚焦找到玻片当前视野的焦平面图像,确保采集到的视野图像中标定线轮廓清晰可见;
步骤二:设定运动方向沿玻片长边方向为X方向、短边方向为Y方向、玻片左上角为原点O,其中标定玻片的标定线呈栅格式分布,标定线分别平行于X、Y方向;
步骤三:在相机中采集玻片视野,视野中包含标定线;
步骤四:以步骤三中采集到的视野作为首个视野,扫描平台沿X方向运动采集多个视野,所采集的视野个数大于3个;
步骤五:将步骤四中采集到的除首视野以外的其余所有视野进行图像识别融合拼接,利用图像拟合直线算法拟合得到拼接图像中X方向标定线的两端点坐标(x0,y0)和(x1,y1),得到首视野标定线的方程:AXx+BXyCX=0,式中,AX、BX、CX均为用X方向标定线两端点坐标表示的常数,AX=y1-y0,BX=x0-x1,CX=(y0-y1)x0+(x1-x0)y0;
步骤六:回到首视野,扫描平台沿Y方向运动采集多个视野,所采集的视野个数大于3个;将所采集到的除首视野以外的其余所有视野进行图像识别融合拼接,利用图像拟合直线算法拟合得到拼接图像中Y方向标定线的两端点坐标(x′0,y′0)和(x′1,y′1),得到该标定线的方程:AYx+BYy+CY=0,式中,AY、BY、CY均为用Y方向标定线两端点坐标表示的常数,AY=y′1-y′0,BY=x′0-x′1,CY=(y′0-y′1)x′0+(x′1-x′0)y′0;
步骤七:依靠X、Y方向标定线存在的斜率数值得到X、Y方向标定线的夹角θ,即扫描平台X、Y运动方向的夹角;
步骤八:根据夹角θ进行判定是否进行调整,当0=90°时,扫描平台XY方向相互垂直,直接进行后续采集图像操作;当θ≠90°时,假定扫描平台XY方向的夹角所给定的误差范围为90°±5°,若角度θ超出该误差范围,则安装人员需要根据角度θ作为调整参考对扫描平台的机械结构方面做出调整或重新安装;若角度θ在该误差范围内,则采用运动补偿的方式对实际扫描轨迹做出调整;运动补偿的具体措施为:调整扫描平台X方向与相机坐标系X轴平行,在Y方向采集图像时,扫描平台Y方向步长由原有步长d调整为d·sinθ,同时根据θ是否大于90°而向X方向进行正补偿或负补偿,所补偿距离为Δl=λ·d·cosθ。
与现有技术相比,本发明的有益效果为:
本发明提出一种基于图像显微扫描平台XY方向夹角测量与运动补偿方法,通过测量平台运动X和Y两个方向的夹角是否满足要求,如不符合要求则需要做出相应调整以保证后期图像拼接的准确性和运动多视野采集图像的稳定性。本发明为显微自动扫描平台安装调整的精度与采图的稳定提供了新的思路,能够帮助仪器装配人员方便、直观地判定测量扫描平台运动方向安装的夹角,最终目的在于根据所测夹角重新调整平台结构装配,使扫描平台沿着X轴和Y轴两个方向运动时保持运行方向垂直或者根据运动方向的夹角对平台沿着扫描轨迹运动采图做出对应的运动补偿,以此可以避免显微扫描平台沿扫描轨迹采集图像时的图像偏移,尽可能确保扫描平台沿着X轴和Y轴两个方向运动时保持运行方向垂直,解决了两运动方向安装垂直度不达标导致扫描平台无法准确完成所采集视野图像的高效拼接的问题,可有效提高自动扫描过程图像处理效率,降低后期图像拼接的流程实现难度,降低技术装配人员的工作难度与强度的同时保证了调试的效率和精度,对维持自动显微扫描仪的稳定采图有重要意义,在医疗病理类玻片显微视觉自动检测行业有着广泛应用前景。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明中记载的一些实施例,对于本领域普通技术人员来讲,还可以根据这些附图获得其他的附图。
图1为扫描平台结构装配图。
图2为XY两运动方向不垂直时玻片标定线与实际采集标定线对比示意图。
图3为带有特制标定线的标定玻片。
图4为实际视野采集经图像识别融合拼接示意图。
图5为本发明实施例提供的带有标定线的特制玻片首视野示意图。
图6为本发明实施例提供的玻片视野图像X方向识别拼接对比示意。
图7为本发明实施例提供的标定线加权最小二乘法拟合方法流程图。
图8为本发明实施例提供的玻片视野图像Y方向识别拼接对比示意图。
图9为本发明实施例提供的不同斜率对应两方向拼接标定线夹角情况示意图。
图10为本发明实施例提供的玻片在平台上的扫描轨迹示意图。
图11为本发明实施例提供的扫描平台Y向轨迹偏移及对应的运动补偿示意图。
图12为本发明实施例提供的基于图像显微扫描平台XY方向夹角测量与运动补偿方法流程图。
具体实施方式
下面结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。本发明中的实施例,本领域技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明提出一种基于图像显微扫描平台XY方向夹角测量与运动补偿方法,基于特制标定玻片视野标定线图像识别拼接后再拟合的方法来测量当前扫描平台运动方向的夹角,在理想状态下扫描平台沿X方向与其沿Y方向步进运动夹角为90°,所拍摄的不同视野图像经过图像拼接之后的图像边缘为平齐的,且拼接后图像内容与实际玻片标本是一致的。
如图5所示,假定运动方向沿玻片长边方向为X方向、短边方向为Y方向、玻片左上角为原点O,对带有标定线的玻片进行视野图像采集。
其中标定线呈栅格式分布,分别平行于X、Y方向。理想状态下标定玻片上的X方向标定线图像在该坐标系下的斜率k=0,X方向第一条标定线可表示为:
By+C=0
B和C为该标定线在坐标系下表示的常数且B≠0。
其余X方向均可用类似上式的方程表达。同理Y方向标定线也能够类似Y方向第一条标定线的方程表达,Y方向第一条标定线方程表示为:
A′x+C′=0
A′和C′为该标定线在坐标系下表示的常数且A′≠0。
在相机中采集玻片视野,以图中视野1作为第一个视野。以首视野为参考沿X方向采集若干连续视野图像,视野个数由调整人员指定,一般以大于3个视野为宜。并将采集到的视野图像经图像识别拼接在一起,如图6所示。
由图6可知,X方向是否水平对应的视野图像拼接虽能使标定线共线,但是会影响视野图像的上下边缘平齐与否,以及影响水平标定线与实际标定线是否重合。
针对拼接图像的标定线进行直线拟合能够确定X方向拼接图像标定线的两端点坐标,从而得到该方向下拼接图像标定线斜率,利用同一方向的多个视野图像的拼接来求取拼接图像的标定线斜率能够有效解决因采集过程或图像识别过程中出现的误差问题。以获取图5中视野1的X方向的拼接图像标定线为例,通过沿该视野X方向采集多张连续图像经图像拼接得到拼接图像,拟合拼接图像中标定线获取端点坐标为(x0,y0)和(x1,y1),若y1≠y0,则水平方向的标定线在坐标系下的方程可表示为:
化简得到:
AXx+BXy+CX=0
式中,AX=y1-y0,BX=x0-x1,CX=(y0-y1)x0+(x1-x0)y0。若y1=y0使用上式方程表达同样成立。
本发明中图像视野中标志线拟合采用加权最小二乘法拟合实现,首先利用Sobel算法提取边缘轮廓,确定参与拟合的轮廓点数目,目的是提高拟合效率,根据选定的轮廓点数目从所有轮廓点集合中随机选取点进行近似拟合一条直线,然后计算原始轮廓上每一点到拟合直线的距离,再根据统计学中标准差方式和所设定的参数阈值确定异常距离范围,若某一点到拟合直线的距离在此异常距离范围内则将该点视为异常点并进行去除。假设有100个点参与拟合直线,其中有20个点是偏离直线的异常点,那么通过标准差的方式100个点对应100个统计值,按照距离直线远近从大到小排序后将异常点剔除。如果都参与拟合直线,拟合的直线和实际会有偏差,原因是其中混杂了异常值,所以为了确保拟合的准确度,该方法可以通过设定参数阈值控制要排除的异常点数量,值越小,参与拟合时出现的异常点越多,拟合的准确度越低。由于各轮廓点到拟合直线的距离数值分布近似服从正态分布,根据正态分布的3σ准则,参数阈值一般设定为1~3。为了有效提高拟合精度,还需要合理控制迭代次数,迭代次数默认选定为5次。迭代次数越多,参与拟合的点越少,拟合的准确度越低;迭代次数越少,参与拟合的异常点越多,拟合的准确度也越低。标定线拟合方法流程如图7所示。
而后回到视野1,以扫描平台Y方向进行步进运动,同样采集多张视野图像,并进行图像识别融合拼接。拼接示意如图8所示,Y方向是否竖直对应的视野图像拼接虽能够使标定线共线,但是会影响视野图像的左右边缘平齐与否,以及影响竖直方向的标定线与实际标定线是否重合。
对Y方向拼接的视野图像进行拟合得到拼接图像中的标定线端点坐标为(x′0,y′0)和(x′1,y′1),若x′0≠x′1,则竖直方向的标定线在坐标系下的方程可表示为:
化简得到:
AYx+BYy+CY=0
式中,AY=y′1-y′0,BY=x′0-x′1,CY=(y′0-y′1)x′0+(x′1-x′0)y′0。若x′0=x′1使用上式方程表达同样成立。
因此,以视野1位首视野沿X、Y方向所采集到的视野图像并进行拼接得到的标定线方程可以用如下方程组来表示:
式中,AX、BX、CX均为用X方向标定线两端点坐标表示的常数,AY、BY、CY均为用Y方向标定线两端点坐标表示的常数。
在图像拼接过程中,由图像识别融合拼接原理可知,以视野1分别作为X方向、Y方向的首视野并参与到拼接所得到的X方向、Y方向的标定线必定垂直,因此在实际拼接过程中,X、Y方向的拼接应当在剔除视野1后所在方向上所采集的所有视野进行拼接。因视野中包含两条自相垂直的标定线,所以X方向的标定线斜率必存在,即BX≠0。若BY≠0,则方程组用斜截式直线方程表达为:
从而得到根据kX、kY可得到X、Y方向标定线的夹角,即扫描平台X、Y运动方向的夹角角度,根据角度可以选择平台调整的方式。
以kX、kY为斜率的两直线,求取其所成夹角θ可得到如下四种情况,如图9所示。
(1)当kX>0,kY>0时:
θ=arctan(kX)+(180°-arctan(kY))
(2)当kX>0,kY<0时:
θ=arctan(kX)+(180°-arctan(kY))
(3)当kX<0,kY>0时:
θ=arctan(kX)-arctan(kY)
(4)当kX<0,kY<0时:
θ=arctan(kX)-arctan(kY)
从而得到X、Y方向标定线的夹角θ:
而若BY=0,即视野图像所拟合的Y方向标定线斜率不存在,但其对应的角度为90°,所以此时X、Y方向标定线的夹角θ:
由此得到扫描平台X、Y方向的所成夹角θ:
当θ=90°时,即扫描平台XY方向垂直,此时则无需进行其他调整操作,可直接进行后续图像采集工作。当θ≠90°时,假定扫描平台XY运动方向的夹角所给定的误差范围为90°±5°,若角度θ超出该误差范围,则安装人员需要根据角度θ作为调整参考对扫描平台的机械结构方面做出调整或重新安装;若角度θ在该误差范围内,则可采用运动补偿的方式来消除该角度在扫描平台运动过程中所带来的偏差。
理想情况下,即θ=90°且扫描平台运动坐标系与相机坐标系平行,采用S形轨迹采集图像,扫描平台的运动轨迹如图10所示。而实际情况下扫描平台XY方向夹角θ往往不等于90°,在平台运动坐标系X方向与相机坐标系的X方向平行的前提下,当平台沿X方向轨迹采集图像完毕后再向Y方向步进时,运动会发现偏移,如图11所示。
图11中,角θ为扫描平台XY方向夹角,d为Y向运动步长,α和Δl分别为Y方向在扫描平台XY方向不垂直情况(图中实线箭头)与垂直情况(图中虚线箭头)下运动偏移的角度和水平距离。
以图11为例,此时θ>90°,欲得到理想情况下的图像,即扫描平台XY方向垂直情况,则扫描平台在Y方向运动过程中所需要的运动补偿为:Y向步进运动时,所需的步长由原来的d调整为:
d·cosα=d·cos(θ-90°)=d·sinθ
同时向X正方向步进运动调整距离:
Δl=d·sinα=d·sin(θ-90°)=-d·cosθ
而当θ<90°时,Y向步进运动的步长为:
d·cosα==d·cos(90°-θ)==d·sinθ
同时向X正方向步进运动调整距离:
Δl=d·sinα=d·sin(90°-θ)=d·cosθ
因此,规定Δl朝X正方向步进运动调整为正补偿(θ>90°),反之为负补偿(θ<90°),并根据运动方向夹角θ(角度θ在所给定的误差范围内且θ≠90°)可得运动补偿为:
Y向步进距离为:d·sinθ
X向步进距离为:Δl=λ·d·cosθ。式中,λ是与θ相关的常数:θ>90°时,λ=-1;θ<90°时,λ=1。
综上,本发明通过上述分析提供一种基于图像显微扫描平台XY方向夹角测量与运动补偿方法,实现对扫描平台X和Y运动方向的夹角测量,然后对扫描平台做出相应调整或对扫描轨迹做运动补偿,利于扫描平台各个视野图像的稳定采集,具体流程如图12所示,包含如下步骤:
步骤一:将带有特制标定线的标定玻片放置在显微扫描平台上,通过自动聚焦找到玻片当前视野的焦平面图像,确保采集到的视野图像中标定线轮廓清晰可见;
步骤二:设定运动方向沿玻片长边方向为X方向、短边方向为Y方向、玻片左上角为原点O,其中标定玻片的标定线呈栅格式分布,标定线分别平行于X、Y方向;
步骤三:在相机中采集玻片视野,视野中包含标定线;
步骤四:以步骤三中采集到的视野作为首个视野,扫描平台沿X方向运动采集多个视野,所采集的视野个数大于3个;
步骤五:将步骤四中采集到的除首视野以外的其余所有视野进行图像识别融合拼接,利用图像拟合直线算法拟合得到拼接图像中X方向标定线的两端点坐标(x0,y0)和(x1,y1),得到首视野标定线的方程:AXx+BXy+CX=0,式中,AX、BX、CX均为用X方向标定线两端点坐标表示的常数,AX=y1-y0,BX=x0-x1,CX=(y0-y1)·x0+(x1-x0)·y0;
步骤六:回到首视野,扫描平台沿Y方向运动采集多个视野,所采集的视野个数大于3个;将所采集到的除首视野以外的其余所有视野进行图像识别融合拼接,利用图像拟合直线算法拟合得到拼接图像中Y方向标定线的两端点坐标(x′0,y′0)和(x′1,y′1),得到该标定线的方程:AYx+BYy+CY=0,式中,AY、BY、CY均为用Y方向标定线两端点坐标表示的常数,AY=y′1-y′0,BY=x′0-x′1,CY=(y′0-y′1)x′0+(x′1-x′0)y′0;
步骤七:依靠X、Y方向标定线存在的斜率数值得到X、Y方向标定线的夹角θ,即扫描平台X、Y运动方向的夹角;
步骤八:根据夹角θ进行判定是否进行调整,当θ=90°时,扫描平台XY方向相互垂直,直接进行后续采集图像操作;当θ≠90°时,假定扫描平台XY方向的夹角所给定的误差范围为90°±5°,若角度θ超出该误差范围,则安装人员需要根据角度θ作为调整参考对扫描平台的机械结构方面做出调整或重新安装;若角度θ在该误差范围内,则采用运动补偿的方式对实际扫描轨迹做出调整;运动补偿的具体措施为:调整扫描平台X方向与相机坐标系X轴平行,在Y方向采集图像时,扫描平台Y方向步长由原有步长d调整为d·sinθ,同时根据θ是否大于90°而向X方向进行正补偿或负补偿,所补偿距离为Δl=λ·d·cosθ。
本发明所述的一种基于图像的显微扫描平台XY方向夹角测量与运动补偿方法,提出了一种对医学病理玻片显微自动扫描平台安装过程的XY运动方向夹角测量以及根据夹角做出对应的运动补偿思路,为显微自动扫描平台安装调整的精度与采图的稳定提供了新的思路,避免了两运动方向安装垂直度不达标导致扫描平台无法准确完成所采集视野图像的高效拼接的问题,可有效提高自动扫描过程图像处理效率,降低后期图像拼接的流程实现难度,在医疗病理类玻片显微视觉自动检测行业有着广泛应用前景。
结果表明:该基于图像的显微扫描平台XY方向夹角测量与运动补偿方法为后续扫描平台内部结构实际安装精度提供保障,也为扫描平台运动轨迹的精准提供方法支持,使图像稳定采集与拼接的难度得到简化,速率更为高效,可满足玻片样本显微视觉检测的快速采图拼接显示需求,基本可以满足多数医学玻片扫描平台的安装调整使用场景。
以上所述实施例,仅为本申请的具体实施方式,用以说明本申请的技术方案,而非对其限制,本申请的保护范围并不局限于此,尽管参照前述实施例对本申请进行了详细的说明,本领域的普通技术人员应当理解:任何熟悉本技术领域的技术人员在本申请揭露的技术范围内,其依然可以对前述实施例所记载的技术方案进行修改或可轻易想到变化,或者对其中部分技术特殊进行等同替换;而这些修改、变化或者替换,并不使相应技术方案的本质脱离本申请实施例技术方案的精神和范围。都应涵盖在本申请的保护范围之内。因此,本申请的保护范围应所述以权利要求的保护范围为准。
Claims (4)
1.一种基于图像显微扫描平台XY方向夹角测量与运动补偿方法,其特征在于,包括以下步骤:
步骤一:将带有特制标定线的标定玻片放置在显微扫描平台上,通过自动聚焦找到玻片当前视野的焦平面图像,确保采集到的视野图像中标定线轮廓清晰可见;
步骤二:设定运动方向沿玻片长边方向为X方向、短边方向为Y方向、玻片左上角为原点O,其中标定玻片的标定线呈栅格式分布,标定线分别平行于X、Y方向;
步骤三:在相机中采集玻片视野,视野中包含标定线;
步骤四:以步骤三中采集到的视野作为首个视野,扫描平台沿X方向运动采集多个视野,所采集的视野个数大于3个;
步骤五:将步骤四中采集到的除首视野以外的其余所有视野进行图像识别融合拼接,利用图像拟合直线算法拟合得到拼接图像中X方向标定线的两端点坐标(x0,y0)和(x1,y1),得到首视野标定线的方程:AXx+BXy+CX=0,式中,AX、BX、CX均为用X方向标定线两端点坐标表示的常数,AX=y1-y0,BX=x0-x1,CX=(y0-y1)x0+(x1-x0)y0;
步骤六:回到首视野,扫描平台沿Y方向运动采集多个视野,所采集的视野个数大于3个;将所采集到的除首视野以外的其余所有视野进行图像识别融合拼接,利用图像拟合直线算法拟合得到拼接图像中Y方向标定线的两端点坐标(x′0,y′0)和(x′1,y′1),得到该标定线的方程:AYx+BYy+CY=0,式中,AY、BY、CY均为用Y方向标定线两端点坐标表示的常数,AY=y′1-y′0,BY=x′0-x′1,CY=(y′0-y′1)x′0+(x′1-x′0)y′0;
步骤七:依靠X、Y方向标定线存在的斜率数值得到X、Y方向标定线的夹角θ,即扫描平台X、Y运动方向的夹角;
步骤八:根据夹角θ进行判定是否进行调整,当θ=90°时,扫描平台XY方向相互垂直,直接进行后续采集图像操作;当θ≠90°时,假定扫描平台XY方向的夹角所给定的误差范围为90°±5°,若角度θ超出该误差范围,则安装人员根据角度θ作为调整参考对扫描平台的机械结构方面做出调整或重新安装;若角度θ在该误差范围内,则采用运动补偿的方式对实际扫描轨迹做出调整;运动补偿的具体措施为:调整扫描平台X方向与相机坐标系X轴平行,在Y方向采集图像时,扫描平台Y方向步长由原有步长d调整为d·sinθ,同时根据θ是否大于90°而向X方向进行正补偿或负补偿,所补偿距离为Δl=λ·d·cosθ。
2.根据权利要求1所述的基于图像显微扫描平台XY方向夹角测量与运动补偿方法,其特征在于,步骤五中图像视野中标志线拟合采用加权最小二乘法拟合实现,首先利用Sobel算法提取边缘轮廓,确定参与拟合的轮廓点数目,确定拟合迭代次数,根据选定的轮廓点数目从轮廓点集合中随机选取点近似拟合一条直线,然后计算原始轮廓上每一点到拟合直线的距离,再根据统计学标准差和参数阈值方式确定异常距离范围去除异常点。
4.根据权利要求1所述的基于图像显微扫描平台XY方向夹角测量与运动补偿方法,其特征在于,步骤八中,θ>90°时,λ=-1;θ<90°时,λ=1。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210493682.9A CN114820761B (zh) | 2022-05-07 | 2022-05-07 | 基于图像显微扫描平台xy方向夹角测量与运动补偿方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210493682.9A CN114820761B (zh) | 2022-05-07 | 2022-05-07 | 基于图像显微扫描平台xy方向夹角测量与运动补偿方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114820761A true CN114820761A (zh) | 2022-07-29 |
CN114820761B CN114820761B (zh) | 2024-05-10 |
Family
ID=82511676
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210493682.9A Active CN114820761B (zh) | 2022-05-07 | 2022-05-07 | 基于图像显微扫描平台xy方向夹角测量与运动补偿方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114820761B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115876148A (zh) * | 2023-02-07 | 2023-03-31 | 昂坤视觉(北京)科技有限公司 | 一种晶圆平面度测量平台运行误差补偿方法及系统 |
CN116643393A (zh) * | 2023-07-27 | 2023-08-25 | 南京木木西里科技有限公司 | 基于显微图像偏转的处理方法及系统 |
CN117764940A (zh) * | 2023-12-19 | 2024-03-26 | 珠海圣美生物诊断技术有限公司 | 一种显微镜状态检测方法、装置、计算机设备及存储介质 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20190158813A1 (en) * | 2016-06-10 | 2019-05-23 | Lucid VR, Inc. | Real Time Re-Calibration of Stereo Cameras |
CN111402119A (zh) * | 2020-03-19 | 2020-07-10 | 深圳市生强科技有限公司 | 全景扫描成像方法、装置、玻片扫描仪及存储介质 |
CN111610621A (zh) * | 2020-01-19 | 2020-09-01 | 北京大学 | 一种双模态显微成像系统和方法 |
CN113362276A (zh) * | 2021-04-26 | 2021-09-07 | 广东大自然家居科技研究有限公司 | 板材视觉检测方法及系统 |
CN113643436A (zh) * | 2021-08-24 | 2021-11-12 | 凌云光技术股份有限公司 | 一种深度数据拼接融合方法及装置 |
-
2022
- 2022-05-07 CN CN202210493682.9A patent/CN114820761B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20190158813A1 (en) * | 2016-06-10 | 2019-05-23 | Lucid VR, Inc. | Real Time Re-Calibration of Stereo Cameras |
CN111610621A (zh) * | 2020-01-19 | 2020-09-01 | 北京大学 | 一种双模态显微成像系统和方法 |
CN111402119A (zh) * | 2020-03-19 | 2020-07-10 | 深圳市生强科技有限公司 | 全景扫描成像方法、装置、玻片扫描仪及存储介质 |
CN113362276A (zh) * | 2021-04-26 | 2021-09-07 | 广东大自然家居科技研究有限公司 | 板材视觉检测方法及系统 |
CN113643436A (zh) * | 2021-08-24 | 2021-11-12 | 凌云光技术股份有限公司 | 一种深度数据拼接融合方法及装置 |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115876148A (zh) * | 2023-02-07 | 2023-03-31 | 昂坤视觉(北京)科技有限公司 | 一种晶圆平面度测量平台运行误差补偿方法及系统 |
CN115876148B (zh) * | 2023-02-07 | 2023-05-19 | 昂坤视觉(北京)科技有限公司 | 一种晶圆平面度测量平台运行误差补偿方法及系统 |
CN116643393A (zh) * | 2023-07-27 | 2023-08-25 | 南京木木西里科技有限公司 | 基于显微图像偏转的处理方法及系统 |
CN116643393B (zh) * | 2023-07-27 | 2023-10-27 | 南京木木西里科技有限公司 | 基于显微图像偏转的处理方法及系统 |
CN117764940A (zh) * | 2023-12-19 | 2024-03-26 | 珠海圣美生物诊断技术有限公司 | 一种显微镜状态检测方法、装置、计算机设备及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN114820761B (zh) | 2024-05-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN114820761A (zh) | 基于图像显微扫描平台xy方向夹角测量与运动补偿方法 | |
CN108613628B (zh) | 一种基于双目视觉的架空输电线路弧垂测量方法 | |
CN110006905B (zh) | 一种线面阵相机结合的大口径超净光滑表面缺陷检测装置 | |
CN112629441B (zh) | 3d曲面玻璃轮廓扫描检测方法及系统 | |
EP2666050B1 (en) | Microscope slide coordinate system registration | |
JP5997989B2 (ja) | 画像測定装置、その制御方法及び画像測定装置用のプログラム | |
JPH07193109A (ja) | ウェーハ粒子解析に関する多重走査法 | |
CN111640154B (zh) | 一种基于显微视觉的立针微平面亚像素级定位方法 | |
CN114578540B (zh) | 基于图像技术的显微扫描载物台与物镜垂直度的调整方法 | |
CN113516716B (zh) | 一种单目视觉位姿测量及调整方法和系统 | |
CN111091121A (zh) | 一种基于图像处理的椭圆表盘检测矫正的方法 | |
CN110657750A (zh) | 一种用于刀具刃口钝化的检测系统和方法 | |
US11536943B2 (en) | Microscope and method for generating an image pieced together from a plurality of individual microscope images | |
CN112468800B (zh) | 广角摄像模组的测试方法及其测试系统 | |
CN112468801B (zh) | 广角摄像模组的光心测试方法及其测试系统和测试标板 | |
CN115131286A (zh) | 一种显微扫描平台运动坐标系与相机坐标系平行度调整方法 | |
TWI573098B (zh) | Wafer defect detection method | |
CN110689582A (zh) | 全站仪摄像机标定方法 | |
CN110349169A (zh) | 一种直线测量方法 | |
JP2565585B2 (ja) | スクライブライン交差領域の位置検出方法 | |
CN111578967B (zh) | 合像式一副水准标尺零点差之差的检测装置及其检测方法 | |
CN111550480B (zh) | 一种空间精密位姿对准方法 | |
CN116978005B (zh) | 基于姿态变换的显微镜图像处理系统 | |
CN113808188B (zh) | 元件的拍摄方法 | |
CN110044267A (zh) | 基于机器视觉检测技术的带钢宽度在线测量方法 |
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 |