CN107230227A - 一种基于切线法的椭圆构像偏差估计方法 - Google Patents

一种基于切线法的椭圆构像偏差估计方法 Download PDF

Info

Publication number
CN107230227A
CN107230227A CN201710408118.1A CN201710408118A CN107230227A CN 107230227 A CN107230227 A CN 107230227A CN 201710408118 A CN201710408118 A CN 201710408118A CN 107230227 A CN107230227 A CN 107230227A
Authority
CN
China
Prior art keywords
mrow
msub
mtd
plane
rotation
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
Application number
CN201710408118.1A
Other languages
English (en)
Other versions
CN107230227B (zh
Inventor
张旭
沈毅君
朱利民
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
HUST Wuxi Research Institute
Original Assignee
HUST Wuxi Research Institute
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by HUST Wuxi Research Institute filed Critical HUST Wuxi Research Institute
Priority to CN201710408118.1A priority Critical patent/CN107230227B/zh
Publication of CN107230227A publication Critical patent/CN107230227A/zh
Application granted granted Critical
Publication of CN107230227B publication Critical patent/CN107230227B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/60Analysis of geometric attributes
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/80Analysis of captured images to determine intrinsic or extrinsic camera parameters, i.e. camera calibration

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • Image Analysis (AREA)
  • Studio Devices (AREA)
  • Length Measuring Devices By Optical Means (AREA)

Abstract

本发明公开一种基于切线法的椭圆构像偏差估计方法,包括如下步骤:S101、构造圆平面与像平面的公垂面与旋转轴;S102、获取所述公垂面绕所述旋转轴的两个极限旋转角及其对应的圆上两点;S103、根据所述圆上两点计算出椭圆的两切点,所述两切点的中点为椭圆中心,计算所述椭圆中心与圆心在像平面上的投影之差即椭圆构像偏差。本发明优点如下:一、仅使用世界坐标系与摄像机坐标系,避免了涉及多个坐标系之间的坐标转换,提高了方法的简洁性与精确性。二、计算的速度与精度取决于旋转角度函数的极值计算速度与精度,根据实际应用场合的不同,可以通过调整函数极值的计算方法满足不同的精度、速度要求。

Description

一种基于切线法的椭圆构像偏差估计方法
技术领域
本发明涉及计算机视觉领域,尤其涉及一种基于切线法的椭圆构像偏差估计方法。
背景技术
圆经过透视投影,当成像平面与圆平面不平行时投影为椭圆,然而投影椭圆的中心并非圆心在像平面上的投影,两者的偏差称为椭圆构像偏差。由于圆形标志在摄像机标定中常被用到,而现有的标志点识别算法常将椭圆中心近似作为圆的中心进行计算,在进行高精密度的光学测量时,该近似过程所引入的误差不能被忽略。文献“On the accuracyof videometry”(Lenz R and Fritsch D.International Archives of Photogrammetryand Remote Sensing,vol.27(B5),pp.334-345,1984.)第一次提出当圆形特征尺寸较大,与摄像机距离过近时,椭圆构像的偏心误差会对测量的精确度造成影响。文献“SystematicGeometric Image Measurement Errors of Circular Object Targets:MathematicalFormulation and Correction”(Ahn S J,Warnecke H J,Kotowski R.PhotogrammetricRecord,vol.16(93),pp.485-502,1999.)提出了一种椭圆构像偏差的代数计算方法,但其方法限制了圆形特征的位置,且引入多个中间坐标系,偏心误差数值以近似结果给出。该方法在简便性,通用性和精确性方面均存在不足。
发明内容
本发明的目的在于通过一种基于切线法的椭圆构像偏差估计方法,来解决以上背景技术部分提到的问题。
为达此目的,本发明采用以下技术方案:
一种基于切线法的椭圆构像偏差估计方法,其包括如下步骤:
S101、构造圆平面与像平面的公垂面与旋转轴;
S102、获取所述公垂面绕所述旋转轴的两个极限旋转角即最大旋转角和最小旋转角以及所述最大旋转角和最小旋转角对应的圆上两点;
S103、根据所述圆上两点计算出椭圆的两切点,所述两切点的中点为椭圆中心,计算所述椭圆中心与圆心在像平面上的投影之差即椭圆构像偏差。
特别地,所述步骤S101包括:S1011、构造世界坐标系OXYZ与摄像机坐标系O'X'Y'Z':以圆平面为XY平面建立世界坐标系OXYZ,以光心为摄像机坐标系原点O',光轴为z轴建立摄像机坐标系O'X'Y'Z';S1012、构造与圆平面Po和像平面PI垂直且通过光心O'的平面π即公垂面π,在摄像机坐标系O'X'Y'Z'下,垂直圆平面Po的光线方向是旋转矩阵R的第三列r3,像平面PI的垂直矢量为z轴方向[0 0 1]T,由公式(10)获得圆平面与像平面的公垂面π的法线n的方向:
S1013、构造旋转轴Iimg:旋转轴Iimg的特点是通过光心O'且方向为平面π与像平面PI的交线方向,由公式(11)获得交线方向:
当公垂面π绕旋转轴Iimg旋转θ时,由公式(12)得出公垂面π与像平面PI的交线移动距离d:
d=Io·tanθ (12)
其中,Io为光心到像平面的距离,当公垂面π绕旋转轴旋转形成旋转平面族时,其与像平面PI的交线为平行线族。
特别地,所述步骤S102包括:椭圆中心Oe是两条平行切线与椭圆切点M1、M2的中点,该两条平行切线对应平行线族的两个极限位置,同时对应旋转平面族的两个极限旋转角即最大旋转角和最小旋转角,构造垂直于旋转轴Iimg且通过光心O'的平面,由公式(13)得到圆上任意一点Q在该平面的投影Q':
任意一条垂直于旋转轴Iimg的直线均可作为旋转平面族旋转角的计算基准,当选用公式(1)获得的公垂面π的法线n(即公垂线n)的方向作为角度计算的基准时,由公式(14)得出旋转角:
由公式(15)得出对应椭圆与平行线族的两个切点M1、M2的圆上的两点Q1、Q2
其中,Oc为圆心在摄像机坐标系下的坐标。
特别地,所述步骤S103包括:摄像机透视投影的关系式(16)建立了物点Q与像点M的坐标关系,将Q1、Q2、Oc分别代入式(16)可以得到椭圆的两个切点M1、M2与圆心的像点Mc即圆心的投影:
M=K[r1 r2 r3 T]Q (16)
由公式(17)获得椭圆构像中心Oe即椭圆中心的投影:
根据公式(18)计算出椭圆构像偏差Δ:
Δ=Oe-Mc (18)。
本发明提出的基于切线法的椭圆构像偏差估计方法具有如下优点:一、仅使用世界坐标系与摄像机坐标系,避免了涉及多个坐标系之间的坐标转换,提高了方法的简洁性与精确性。二、计算的速度与精度取决于旋转角度函数的极值计算速度与精度,根据实际应用场合的不同,可以通过调整函数极值的计算方法满足不同的精度、速度要求。
附图说明
图1为本发明实施例提供的基于切线法的椭圆构像偏差估计方法流程图;
图2为本发明实施例提供的世界坐标系与摄像机坐标系示意图;
图3为本发明实施例提供的公垂面与旋转轴示意图;
图4为本发明实施例提供的极限旋转角示意图。
具体实施方式
下面结合附图和实施例对本发明作进一步说明。可以理解的是,此处所描述的具体实施例仅仅用于解释本发明,而非对本发明的限定。另外还需要说明的是,为了便于描述,附图中仅示出了与本发明相关的部分而非全部内容,除非另有定义,本文所使用的所有技术和科学术语与属于本发明的技术领域的技术人员通常理解的含义相同。本文中所使用的术语只是为了描述具体的实施例,不是旨在于限制本发明。
请参照图1所示,图1为本发明实施例提供的基于切线法的椭圆构像偏差估计方法流程图。
本实施例中基于切线法的椭圆构像偏差估计方法具体包括如下步骤:
S101、构造圆平面与像平面的公垂面与旋转轴。具体过程如下:
S1011、构造世界坐标系OXYZ与摄像机坐标系O'X'Y'Z':以圆平面为XY平面建立世界坐标系OXYZ,以光心为摄像机坐标系原点O',光轴为z轴建立摄像机坐标系O'X'Y'Z',如图2所示,201指圆平面Po,202指像平面PI,203指世界坐标系OXYZ,204指摄像机坐标系O'X'Y'Z'。
S1012、由圆平面Po的法线与像平面PI的法线做叉乘运算得到其公垂面π的法线n,由法线n和像平面PI的法线做叉乘运算,得到圆平面Po与像平面PI的交线即旋转轴Iimg。如图3所示,图中301指平面π(同公垂面π)。具体过程如下:
构造与圆平面Po和像平面PI垂直且通过光心O'的平面π即公垂面π,在摄像机坐标系O'X'Y'Z'下,垂直圆平面Po的光线方向是旋转矩阵R的第三列r3,像平面PI的垂直矢量为z轴方向[0 0 1]T,由公式(10)获得圆平面与像平面的公垂面π的法线n的方向:
S1013、构造旋转轴Iimg:旋转轴Iimg的特点是通过光心O'且方向为平面π与像平面PI的交线方向,由公式(11)获得交线方向:
当公垂面π绕旋转轴Iimg旋转θ时,由公式(12)得出公垂面π与像平面PI的交线移动距离d:
d=Io·tanθ (21)
其中,Io为光心到像平面的距离,当公垂面π绕旋转轴旋转形成旋转平面族时,其与像平面PI的交线为平行线族。
S102、获取所述公垂面绕所述旋转轴的两个极限旋转角即最大旋转角和最小旋转角以及所述最大旋转角和最小旋转角对应的圆上两点。
如图4所示,求出圆上任意一点Q在垂直于旋转轴的平面上的投影Q',向量O'Q'与向量n(法线n的方向)之间的夹角即为旋转平面族各平面的旋转角度,求出该函数的最大值与最小值所对应的圆上的点Q。具体过程如下:
椭圆中心Oe是两条平行切线与椭圆切点M1、M2的中点,该两条平行切线对应平行线族的两个极限位置,同时对应旋转平面族的两个极限旋转角即最大旋转角和最小旋转角,构造垂直于旋转轴Iimg且通过光心O'的平面,由公式(13)得到圆上任意一点Q在该平面的投影Q':
任意一条垂直于旋转轴Iimg的直线均可作为旋转平面族旋转角的计算基准,当选用公式(1)获得的公垂面π的法线n(即公垂线n)的方向作为角度计算的基准时,由公式(14)得出旋转角:
由公式(15)得出对应椭圆与平行线族的两个切点M1、M2的圆上的两点Q1、Q2
其中,Oc为圆心在摄像机坐标系下的坐标。
S103、根据所述圆上两点计算出椭圆的两切点,所述两切点的中点为椭圆中心,计算所述椭圆中心与圆心在像平面上的投影之差即椭圆构像偏差。具体过程如下:
摄像机透视投影的关系式(16)建立了物点Q与像点M的坐标关系,将Q1、Q2、Oc分别代入式(16)可以得到椭圆的两个切点M1、M2与圆心的像点Mc即圆心的投影:
M=K[r1 r2 r3 T]Q (25)
由公式(17)获得椭圆构像中心Oe即椭圆中心的投影:
根据公式(18)计算出椭圆构像偏差Δ:
Δ=Oe-Mc (27)。
本发明的技术方案具有如下优点:一、仅使用世界坐标系与摄像机坐标系,避免了涉及多个坐标系之间的坐标转换,提高了方法的简洁性与精确性。二、计算的速度与精度取决于旋转角度函数的极值计算速度与精度,根据实际应用场合的不同,可以通过调整函数极值的计算方法满足不同的精度、速度要求。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程,是可以通过计算机程序来指令相关的硬件来完成,所述的程序可存储于一计算机可读取存储介质中,该程序在执行时,可包括如上述各方法的实施例的流程。其中,所述的存储介质可为磁碟、光盘、只读存储记忆体(Read-Only Memory,ROM)或随机存储记忆体(Random AccessMemory,RAM)等。
以上结合具体实施例描述了本发明的技术原理。这些描述只是为了解释本发明的原理,而不能以任何方式解释为对本发明保护范围的限制。基于此处的解释,本领域的技术人员不需要付出创造性的劳动即可联想到本发明的其它具体实施方式,这些方式都将落入本发明的保护范围之内。

Claims (4)

1.一种基于切线法的椭圆构像偏差估计方法,其特征在于,包括如下步骤:
S101、构造圆平面与像平面的公垂面与旋转轴;
S102、获取所述公垂面绕所述旋转轴的两个极限旋转角即最大旋转角和最小旋转角以及所述最大旋转角和最小旋转角对应的圆上两点;
S103、根据所述圆上两点计算出椭圆的两切点,所述两切点的中点为椭圆中心,计算所述椭圆中心与圆心在像平面上的投影之差即椭圆构像偏差。
2.根据权利要求1所述的基于切线法的椭圆构像偏差估计方法,其特征在于,所述步骤S101包括:S1011、构造世界坐标系OXYZ与摄像机坐标系O'X'Y'Z':以圆平面为XY平面建立世界坐标系OXYZ,以光心为摄像机坐标系原点O',光轴为z轴建立摄像机坐标系O'X'Y'Z';S1012、构造与圆平面Po和像平面PI垂直且通过光心O'的平面π即公垂面π,在摄像机坐标系O'X'Y'Z'下,垂直圆平面Po的光线方向是旋转矩阵R的第三列r3,像平面PI的垂直矢量为z轴方向[0 0 1]T,由公式(10)获得圆平面与像平面的公垂面π的法线n的方向:
<mrow> <mi>n</mi> <mo>=</mo> <msup> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>1</mn> </mtd> </mtr> </mtable> </mfenced> <mi>T</mi> </msup> <mo>&amp;CircleTimes;</mo> <msub> <mi>r</mi> <mn>3</mn> </msub> <mo>=</mo> <msup> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <mo>-</mo> <msub> <mi>r</mi> <mn>32</mn> </msub> </mrow> </mtd> <mtd> <msub> <mi>r</mi> <mn>31</mn> </msub> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> </mtable> </mfenced> <mi>T</mi> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow>
S1013、构造旋转轴Iimg:旋转轴Iimg的特点是通过光心O'且方向为平面π与像平面PI的交线方向,由公式(11)获得交线方向:
<mrow> <msub> <mi>l</mi> <mrow> <mi>i</mi> <mi>m</mi> <mi>g</mi> </mrow> </msub> <mo>=</mo> <mi>n</mi> <mo>&amp;CircleTimes;</mo> <msup> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>1</mn> </mtd> </mtr> </mtable> </mfenced> <mi>T</mi> </msup> <mo>=</mo> <msup> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msub> <mi>r</mi> <mn>31</mn> </msub> </mtd> <mtd> <msub> <mi>r</mi> <mn>32</mn> </msub> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> </mtable> </mfenced> <mi>T</mi> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow>
当公垂面π绕旋转轴Iimg旋转θ时,由公式(12)得出公垂面π与像平面PI的交线移动距离d:
d=Io·tanθ (3)
其中,Io为光心到像平面的距离,当公垂面π绕旋转轴旋转形成旋转平面族时,其与像平面PI的交线为平行线族。
3.根据权利要求2所述的基于切线法的椭圆构像偏差估计方法,其特征在于,所述步骤S102包括:椭圆中心Oe是两条平行切线与椭圆切点M1、M2的中点,该两条平行切线对应平行线族的两个极限位置,同时对应旋转平面族的两个极限旋转角即最大旋转角和最小旋转角,构造垂直于旋转轴Iimg且通过光心O'的平面,由公式(13)得到圆上任意一点Q在该平面的投影Q':
<mrow> <msup> <mi>Q</mi> <mo>&amp;prime;</mo> </msup> <mo>=</mo> <mi>Q</mi> <mo>-</mo> <mi>d</mi> <mi>o</mi> <mi>t</mi> <mrow> <mo>(</mo> <mi>Q</mi> <mo>,</mo> <msub> <mi>l</mi> <mrow> <mi>i</mi> <mi>m</mi> <mi>g</mi> </mrow> </msub> <mo>)</mo> </mrow> <mfrac> <msub> <mi>l</mi> <mrow> <mi>i</mi> <mi>m</mi> <mi>g</mi> </mrow> </msub> <mrow> <mo>|</mo> <mo>|</mo> <msub> <mi>l</mi> <mrow> <mi>i</mi> <mi>m</mi> <mi>g</mi> </mrow> </msub> <mo>|</mo> <msup> <mo>|</mo> <mn>2</mn> </msup> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mrow>
任意一条垂直于旋转轴Iimg的直线均可作为旋转平面族旋转角的计算基准,当选用公式(1)获得的公垂面π的法线n即公垂线n的方向作为角度计算的基准时,由公式(14)得出旋转角:
<mrow> <mi>&amp;theta;</mi> <mo>=</mo> <mi>arccos</mi> <mrow> <mo>(</mo> <mfrac> <mrow> <mi>d</mi> <mi>o</mi> <mi>t</mi> <mrow> <mo>(</mo> <msup> <mi>Q</mi> <mo>&amp;prime;</mo> </msup> <mo>,</mo> <mi>n</mi> <mo>)</mo> </mrow> </mrow> <mrow> <mo>|</mo> <mo>|</mo> <msup> <mi>Q</mi> <mo>&amp;prime;</mo> </msup> <mo>|</mo> <mo>|</mo> </mrow> </mfrac> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>5</mn> <mo>)</mo> </mrow> </mrow>
由公式(15)得出对应椭圆与平行线族的两个切点M1、M2的圆上的两点Q1、Q2
<mrow> <mtable> <mtr> <mtd> <mrow> <msub> <mi>Q</mi> <mn>1</mn> </msub> <mo>=</mo> <munder> <mrow> <mi>arg</mi> <mi>min</mi> </mrow> <mrow> <mo>|</mo> <mo>|</mo> <mi>Q</mi> <mo>-</mo> <msub> <mi>O</mi> <mi>c</mi> </msub> <mo>|</mo> <mo>|</mo> <mo>=</mo> <mi>r</mi> </mrow> </munder> <mrow> <mo>{</mo> <mrow> <mi>&amp;theta;</mi> <mrow> <mo>(</mo> <mi>Q</mi> <mo>)</mo> </mrow> </mrow> <mo>}</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>Q</mi> <mn>2</mn> </msub> <mo>=</mo> <munder> <mrow> <mi>arg</mi> <mi>min</mi> </mrow> <mrow> <mo>|</mo> <mo>|</mo> <mi>Q</mi> <mo>-</mo> <msub> <mi>O</mi> <mi>c</mi> </msub> <mo>|</mo> <mo>|</mo> <mo>=</mo> <mi>r</mi> </mrow> </munder> <mrow> <mo>{</mo> <mrow> <mi>&amp;theta;</mi> <mrow> <mo>(</mo> <mi>Q</mi> <mo>)</mo> </mrow> </mrow> <mo>}</mo> </mrow> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo> </mrow> </mrow> 1
其中,Oc为圆心在摄像机坐标系下的坐标。
4.根据权利要求3所述的基于切线法的椭圆构像偏差估计方法,其特征在于,所述步骤S103包括:摄像机透视投影的关系式(16)建立了物点Q与像点M的坐标关系,将Q1、Q2、Oc分别代入式(16)可以得到椭圆的两个切点M1、M2与圆心的像点Mc即圆心的投影:
M=K[r1 r2 r3 T]Q (7)
由公式(17)获得椭圆构像中心Oe即椭圆中心的投影:
<mrow> <msub> <mi>O</mi> <mi>e</mi> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <mrow> <mo>(</mo> <msub> <mi>M</mi> <mn>1</mn> </msub> <mo>+</mo> <msub> <mi>M</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>8</mn> <mo>)</mo> </mrow> </mrow>
根据公式(18)计算出椭圆构像偏差Δ:
Δ=Oe-Mc (9)。
CN201710408118.1A 2017-06-02 2017-06-02 一种基于切线法的椭圆构像偏差估计方法 Active CN107230227B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710408118.1A CN107230227B (zh) 2017-06-02 2017-06-02 一种基于切线法的椭圆构像偏差估计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710408118.1A CN107230227B (zh) 2017-06-02 2017-06-02 一种基于切线法的椭圆构像偏差估计方法

Publications (2)

Publication Number Publication Date
CN107230227A true CN107230227A (zh) 2017-10-03
CN107230227B CN107230227B (zh) 2020-11-03

Family

ID=59934329

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710408118.1A Active CN107230227B (zh) 2017-06-02 2017-06-02 一种基于切线法的椭圆构像偏差估计方法

Country Status (1)

Country Link
CN (1) CN107230227B (zh)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1566900A (zh) * 2003-06-11 2005-01-19 北京航空航天大学 一种空间圆几何参数的视觉测量方法
CN101706960A (zh) * 2009-10-30 2010-05-12 北京航空航天大学 一种共面圆圆心投影点的定位方法
US7756308B2 (en) * 2005-02-07 2010-07-13 Stereotaxis, Inc. Registration of three dimensional image data to 2D-image-derived data
CN102915535A (zh) * 2012-08-23 2013-02-06 深圳大学 摄像机投影变换中圆形标志点的圆心偏差修正方法及***
CN103116892A (zh) * 2013-03-15 2013-05-22 云南大学 两个相交相同圆及公切线求解摄像机内参数

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1566900A (zh) * 2003-06-11 2005-01-19 北京航空航天大学 一种空间圆几何参数的视觉测量方法
US7756308B2 (en) * 2005-02-07 2010-07-13 Stereotaxis, Inc. Registration of three dimensional image data to 2D-image-derived data
CN101706960A (zh) * 2009-10-30 2010-05-12 北京航空航天大学 一种共面圆圆心投影点的定位方法
CN102915535A (zh) * 2012-08-23 2013-02-06 深圳大学 摄像机投影变换中圆形标志点的圆心偏差修正方法及***
CN103116892A (zh) * 2013-03-15 2013-05-22 云南大学 两个相交相同圆及公切线求解摄像机内参数

Non-Patent Citations (7)

* Cited by examiner, † Cited by third party
Title
S.J.AHN: "SYSTEMATIC GEOMETRIC IMAGE MEASUREMENT ERRORS OF CIRCULAR OBJECT TARGETS: MATHEMATICAL FORMULATION AND CORRECTION", 《PHOTOGRAMMETRIC RECORD》 *
廖祥春 等: "圆形标志及其椭圆构像间中心偏差的确定", 《武汉测绘科技大学学报》 *
苗锡奎 等: "基于星箭对接环部的飞行器单目视觉位姿测量方法", 《光学学报》 *
谭燕 等: "近景摄影测量圆形标志与其椭圆构像的中心偏差估计", 《测绘科学》 *
韩建栋等: "三维视觉测量中圆中心投影误差分析方法", 《计算机科学》 *
魏振忠等: "透视投影变换中椭圆中心畸变误差模型及其仿真研究", 《仪器仪表学报》 *
黄道明: "近景摄影测量中圆形目标的投影偏心差及模拟分析", 《工程建设》 *

Also Published As

Publication number Publication date
CN107230227B (zh) 2020-11-03

Similar Documents

Publication Publication Date Title
US9928595B2 (en) Devices, systems, and methods for high-resolution multi-view camera calibration
CN109752003B (zh) 一种机器人视觉惯性点线特征定位方法及装置
CN106056664A (zh) 一种基于惯性和深度视觉的实时三维场景重构***及方法
CN110174899A (zh) 一种基于敏捷卫星的高精度成像姿态指向控制方法
CN101621711B (zh) 采用两个相同圆进行摄像机标定的方法
US9549712B2 (en) Method and apparatus for correcting focus of CT device
CN104737205B (zh) 自适应尺度及/或重力估计
CN110197510A (zh) 双目相机的标定方法、装置、无人机和存储介质
CN111123280B (zh) 激光雷达的定位方法、装置、***、电子设备及存储介质
CN103411589B (zh) 一种基于四维实数矩阵的三维图像匹配导航方法
CN104715469A (zh) 一种数据处理方法及电子设备
CN108520543B (zh) 一种对相对精度地图进行优化的方法、设备及存储介质
Wen et al. PLS-VIO: Stereo vision-inertial odometry based on point and line features
Huo et al. Corrected calibration algorithm with a fixed constraint relationship and an error compensation technique for a binocular vision measurement system
WO2019237433A1 (zh) 扫地机器人的摄像头与里程计坐标标定方法及***
WO2022078301A1 (zh) 一种数据处理的方法和装置
US20150260545A1 (en) Gyroscope calibration
CN107230227A (zh) 一种基于切线法的椭圆构像偏差估计方法
WO2023160445A1 (zh) 即时定位与地图构建方法、装置、电子设备及可读存储介质
CN114252073B (zh) 一种机器人姿态数据融合方法
CN103412337A (zh) 基于双独立旋转磁棒电磁跟踪的位置跟踪方法及***
CN110084855A (zh) 一种改进cbct几何参数标定算法
CN116124094A (zh) 基于无人机侦察图像与组合导航信息的多目标协同定位方法
WO2022073172A1 (zh) 一种基于点线特征的机器人全局最优视觉定位方法及装置
Wang et al. High precision and fast estimation of position and attitude measurement for space targets

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
EE01 Entry into force of recordation of patent licensing contract
EE01 Entry into force of recordation of patent licensing contract

Application publication date: 20171003

Assignee: Jiangsu Jihui Huake Intelligent Equipment Technology Co., Ltd.

Assignor: Wuxi research institute of the Central China University of Science and Technology

Contract record no.: X2019980000288

Denomination of invention: Ellipse imaging deviation estimation method based on tangent method

License type: Common License

Record date: 20191023

GR01 Patent grant
GR01 Patent grant