CN116753914A - 一种基于图像修正的高精度倾角测量装置及测量方法 - Google Patents
一种基于图像修正的高精度倾角测量装置及测量方法 Download PDFInfo
- Publication number
- CN116753914A CN116753914A CN202310873057.1A CN202310873057A CN116753914A CN 116753914 A CN116753914 A CN 116753914A CN 202310873057 A CN202310873057 A CN 202310873057A CN 116753914 A CN116753914 A CN 116753914A
- Authority
- CN
- China
- Prior art keywords
- point
- pendulum
- pixel
- shadow
- measuring
- 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
- 238000003702 image correction Methods 0.000 title claims abstract description 23
- 238000000034 method Methods 0.000 title claims abstract description 20
- 238000005259 measurement Methods 0.000 claims abstract description 33
- 230000006698 induction Effects 0.000 claims abstract description 6
- 238000013016 damping Methods 0.000 claims abstract description 5
- 230000003068 static effect Effects 0.000 claims abstract description 4
- 239000003550 marker Substances 0.000 claims description 24
- 230000003287 optical effect Effects 0.000 claims description 15
- 238000003672 processing method Methods 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 5
- 238000007789 sealing Methods 0.000 claims description 4
- PCTMTFRHKVHKIS-BMFZQQSSSA-N (1s,3r,4e,6e,8e,10e,12e,14e,16e,18s,19r,20r,21s,25r,27r,30r,31r,33s,35r,37s,38r)-3-[(2r,3s,4s,5s,6r)-4-amino-3,5-dihydroxy-6-methyloxan-2-yl]oxy-19,25,27,30,31,33,35,37-octahydroxy-18,20,21-trimethyl-23-oxo-22,39-dioxabicyclo[33.3.1]nonatriaconta-4,6,8,10 Chemical compound C1C=C2C[C@@H](OS(O)(=O)=O)CC[C@]2(C)[C@@H]2[C@@H]1[C@@H]1CC[C@H]([C@H](C)CCCC(C)C)[C@@]1(C)CC2.O[C@H]1[C@@H](N)[C@H](O)[C@@H](C)O[C@H]1O[C@H]1/C=C/C=C/C=C/C=C/C=C/C=C/C=C/[C@H](C)[C@@H](O)[C@@H](C)[C@H](C)OC(=O)C[C@H](O)C[C@H](O)CC[C@@H](O)[C@H](O)C[C@H](O)C[C@](O)(C[C@H](O)[C@H]2C(O)=O)O[C@H]2C1 PCTMTFRHKVHKIS-BMFZQQSSSA-N 0.000 claims description 3
- 238000012937 correction Methods 0.000 claims description 3
- 238000003384 imaging method Methods 0.000 claims description 3
- 238000000691 measurement method Methods 0.000 claims 7
- 230000007547 defect Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 230000005484 gravity Effects 0.000 description 2
- 238000009434 installation Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 229910001374 Invar Inorganic materials 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000000919 ceramic Substances 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 208000002173 dizziness Diseases 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process 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
- 238000010248 power generation Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C9/00—Measuring inclination, e.g. by clinometers, by levels
- G01C9/12—Measuring inclination, e.g. by clinometers, by levels by using a single pendulum plumb lines G01C15/10
-
- 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
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B7/00—Measuring arrangements characterised by the use of electric or magnetic techniques
- G01B7/003—Measuring arrangements characterised by the use of electric or magnetic techniques for measuring position, not involving coordinate determination
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C9/00—Measuring inclination, e.g. by clinometers, by levels
- G01C9/02—Details
- G01C9/06—Electric or photoelectric indication or reading means
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H9/00—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves by using radiation-sensitive means, e.g. optical means
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/40—Extraction of image or video features
- G06V10/44—Local feature extraction by analysis of parts of the pattern, e.g. by detecting edges, contours, loops, corners, strokes or intersections; Connectivity analysis, e.g. of connected components
-
- 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/10016—Video; Image sequence
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Theoretical Computer Science (AREA)
- Multimedia (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明提出一种基于图像修正的高精度倾角测量装置及测量方法,采用单摆原理和磁感应原理来测量倾角动态变化,并采用特别设计的算法对测量结果实时标定,消除系统误差,提升测量精度。测量装置包括底座固定板,测量时底座固定板的底面与待测倾角面紧密接触,底座固定板上部依次设有磁场建立单元和标定板,摆锤悬挂在标定板的正上方,磁场建立单元用于产生磁场,通过磁场与摆锤产生磁感应,基于电磁阻尼原理使摆锤的摆动迅速收敛至静止,标定板用于对测量误差进行标定,高清视频摄像头设置在摆锤的正上方,高清视频摄像头用于实时拍摄摆锤的摆动幅度,基于摆锤在图像中的位置,得到摆锤的倾斜角度,待测倾角面的倾角等于摆锤的倾斜角度。
Description
技术领域
本发明涉及图像处理领域,具体是一种基于图像修正的高精度倾角测量装置及测量方法。
背景技术
随着社会经济快速发展,高精度倾角测量在工程领域变得日益重要,尤其是在工业、军工、大型工程等领域,如设备安平、高层建筑物、大型桥梁、风力发电塔筒和军工等领域。目前基于MEMS技术的倾角传感器已得到广泛应用,通过测量加速度来计算物体的倾斜角度并将其进行数字化。与传统水泡式条式测量仪相比,基于MEMS的倾角传感器具有体积小、精度高等优点,但倾角传感器容易受外界温度影响,如在环境温差较大、高振动动等场合,对测量精度影响较大,此外,现有的传感器随着测量的持续,倾角传感器的误差会出现累积,导致测量精度大大降低。
发明内容
本发明的目的是针对现有技术存在的上述问题,克服现有技术的不足之处,设计一种基于图像修正的高精度倾角测量装置,克服现有的传感器测量精度受温度影响的缺陷,消除测量过程中的累积误差。
一种基于图像修正的高精度倾角测量装置,包括上下表面平行的底座固定板,测量时所述底座固定板的底面与待测倾角面紧密接触,所述底座固定板上部依次设有磁场建立单元和标定板,所述标定板与所述底座固定板上表面相平行,摆锤悬挂在所述标定板的正上方,所述磁场建立单元用于产生磁场,通过磁场与所述摆锤产生磁感应,基于电磁阻尼原理使所述摆锤的摆动迅速收敛至静止,所述标定板用于对测量误差进行标定,高清视频摄像头设置在所述摆锤的正上方,所述高清视频摄像头用于实时拍摄所述摆锤的摆动幅度,基于所述摆锤在图像中的位置,得到所述摆锤的倾斜角度,待测倾角面的倾角等于所述摆锤的倾斜角度。
优选的,所述磁场建立单元、所述标定板、所述摆锤和所述高清视频摄像头均设置在低温度系数的密封罩内,消除外界温度变化对测量的影响,提高了倾角的测量精度。
优选的,所述密封罩的顶部设有密封盖,所述高清视频摄像头和所述摆锤固定设置在所述密封盖上,便于安装。
优选的,所述底座固定板上侧中心位置设置有凹槽,所述磁场建立单元固定设置在凹槽内,便于所述标定板与所述底座固定板连接。
优选的,所述标定板上绘制一个圆形,所述摆锤投影落在所述标定板的圆形内,所述高清视频摄像头拍摄范围涵盖所述标定板的边界,所述圆形上设有8个圆形的标志点,任意两个相邻的标志点与圆心的夹角为45°,并且,8个标志点与圆心的距离均相等,所述摆锤是球形,所述标志点的直径小于所述摆锤落在所述投影在所述标定板上投影的直径,便于倾角的计算。
本发明还提供一种基于图像修正的高精度倾角测量装置的测量方法,包括如下步骤:
步骤1:所述底座固定板的底面与待测倾角面紧密接触,
步骤2:高清视频摄像头进行拍摄,获得标定板图像;
步骤3:进行标定板图像中阴影区域的识别;
步骤4:当阴影区域的个数为9个时,判定为8个标志点和1个测量点;当阴影区域的个数为8个时,判定为8个标志点和1个测量点;当阴影区域的个数既不是8个也不是9个时,判定为无效数据,重复步骤2-4;
步骤5:根据标志点和测量点求解摆锤的摆动角度。
进一步的,步骤3包括如下步骤:
步骤3.1:将标定板图像中所有像素点存入像素集合A,对于像素集合A中任一像素点i,将像素点i的灰度值g(i)与灰度阈值δ进行比较,当g(i)>δ时,判定g(i)灰度值不满足要求,将i从像素集合A中删除,重新执行步骤3.1,当g(i)≤δ时,将i从像素集合A中删除,并将像素点i存入阴影集合Pn中,n是执行次数,n=1,2,3,……,继续进行步骤3.2;
步骤3.2:对于阴影集合Pn中像素点i,坐标为(x,y),依次遍历与像素点i相邻的8个邻域像素点:(x-1,y),(x-1,y-1),(x-1,y+1),(x,y-1),(x,y+1),(x+1,y),(x+1,y-1),(x+1,y+1),如果邻域像素点的灰度值≤δ,则将其存入集合Pn中,再依次遍历这些邻域像素点的领域像素点,重复下去直到没有邻域像素点为止;
步骤3.3:重复步骤3.1和3.2,找出标定板图像中所有阴影区域的结合(P1,P2,……,Pk),其中,Pk表示第k个阴影集合,k是阴影区域的个数。
进一步的,步骤2高清视频摄像头进行拍摄时,理论上摄影头光轴与测量平面的夹角为90°,由于外界因素存在摄像头产生倾斜的情况,使得摄像头光轴与测量平面夹角与理论值不符,此时需进行摄像头光轴倾角校正,点O为相机,OM为理论光轴,垂直于平面ABCD,AB垂直于CD,AB作为x轴,CD作为y轴,
假定OM’为实际光轴,将OM’分别投影到OAB平面和OCD平面得到点M1’和点M2’,OM’分别在OAB平面和OCD平面求解∠MOM1’和∠MOM2’,
A’、B’为标志点A和标志点B在摄像头倾斜后图片上的成像,则:
∠MOM′1=arctan(MM′1/OM)=arctan((BM-BM′1)/OM)
同理,
∠MOM′2=arctan(MM′2/OM)=arctan((DM-DM′2)/OM)
其中,BM1’满足下式:
同理,
其中,XA′、XB′分别表示标志点A、B在图片阴影区域的中心点的x轴的坐标,表示点M1’的x轴的坐标,
C’、D’分别表示标志点C、D在图片阴影区域的中心点,YC′、YD′分别表示标志点C、D在图片阴影区域的中心点的y轴的坐标,表示点M2’的y轴的坐标。
进一步的,步骤5进行标志点和测量点的区分具体操作如下:
当k=9时,
首先,分别获取各阴影区域中心点坐标,
然后,分别计算各阴影区域中心点与8个标志点的标定坐标的欧式距离,将获得的9个欧式距离按照从小到大的顺序排列,欧式距离最大的阴影区域为测量点,其余阴影区域为标志点;
当k=8时,
首先,分别获取各阴影区域中心点坐标,
然后,找出包含像素点最多的阴影区域,为测量点和标志点重合的阴影区域,其余阴影区域为标志点,
接下来,找出与重合阴影区域中的标志点关于圆心相对称的标志点,基于对称标志点的中心坐标得到重合阴影区域中的标志点的中心坐标,
提取重合阴影区域的边缘轮廓点集,其中,边缘轮廓点集通过依次遍历重合阴影区域内所有的像素点,当任一像素点的上下左右四个邻域的像素点中至少有1个像素点不在阴影区域集合内,则将该像素点放入轮廓点集中,
进一步,计算边缘轮廓点集与被重合标志点实际坐标P’的欧式距离,从轮廓点集中删除欧式距离与标志点半径差值小于距离阈值的点;
最后,基于剩余轮廓点集拟合圆心,即为测量点的中心坐标。
进一步的,摆锤摆动倾角在OAB平面投影角∠Z1OM,具体计算如下:
其中,表示Z’到M1’的像素距离,/>满足:
其中,OM1’满足:
摆锤摆动倾角在OCD平面投影角∠Z2OM,具体计算如下:
其中,表示Z’到M2’的像素距离,/>满足:
其中,OM2’满足:
与现有技术相比,本发明的有益效果:
本发明的目的是针对现有技术存在的上述问题,克服现有技术的不足之处,设计一种基于图像修正的高精度倾角测量装置,克服现有的传感器测量精度受温度影响的缺陷,消除测量过程中的累积误差。
1.本发明采用单摆原理和磁感应原理来测量倾角动态变化,通过高清摄像机测量重力摆的角度,同时通过电磁阻尼原理使得重力摆迅速收敛到中心位置,测量精度和稳定性高。
2.本发明基于高精度标定板的方式,采用特别设计的算法对测量结果进行实时标定,消除系统误差,提升测量精度。
3.本发明密封罩采用低温度系数材料,克服周围环境温度变化对测量精度的影响,测量结果更可靠。
4.本发明通过自设的算法有效解决了由于摄像头的微晃动,可能导致标志点和测量点在图片成像中出现眩晕效果的情况,模型精度高。
附图说明
图1是本发明设计的一种基于图像修正的高精度倾角测量装置示意图;
图2为本发明的标志点和测量点没有重合的标定板图像示意图;
图3为本发明的标志点和测量点重合的标定板图像示意图;
图4为本发明的摄像头没有倾斜情况下,摄像头与测量平面的理论示意图;
图5为本发明的摄像头倾斜情况下,摄像头与测量平面的示意图;
图6为本发明的摄像头在OAB平面的投影示意图;
图7为本发明的重合阴影区域中的标志点和对称标志点的投影示意图;
图8为本发明的测量点在OAB面的投影示意图。
具体实施方式
下面结合附图以及具体实施方法对本发明一种基于图像修正的高精度倾角测量装置作进一步详细说明。
本发明提出一种基于图像修正的高精度倾角测量装置,如图1所示,包括上下表面平行的底座固定板1,测量时所述底座固定板1的底面与待测倾角面紧密接触,所述底座固定板1上部依次设有磁场建立单元2和标定板3,所述标定板3与所述底座固定板1上表面相平行,摆锤4悬挂在所述标定板3的正上方,所述磁场建立单元2用于产生磁场,通过磁场与所述摆锤4产生磁感应,基于电磁阻尼原理使所述摆锤4的摆动迅速收敛至静止,所述标定板3用于对测量误差进行标定,高清视频摄像头6设置在所述摆锤4的正上方,所述高清视频摄像头6用于实时拍摄所述摆锤4的摆动幅度,基于所述摆锤4在图像中的位置,得到所述摆锤4的倾斜角度,待测倾角面的倾角等于所述摆锤4的倾斜角度。
进一步的,所述磁场建立单元2、所述标定板3、所述摆锤4和所述高清视频摄像头6均设置在低温度系数的密封罩5内,消除外界温度变化对测量的影响,提高了倾角的测量精度。密封罩5优选低温度系数的殷钢、陶瓷等材质制造,克服周围环境温度变化对测量精度的影响,测量结果更可靠。
进一步的,所述密封罩5的顶部设有密封盖7,所述高清视频摄像头6和所述摆锤4固定设置在所述密封盖7上,便于安装。
进一步的,所述底座固定板1上侧中心位置设置有凹槽,所述磁场建立单元2固定设置在凹槽内,便于所述标定板3与所述底座固定板1连接。
优选的,所述标定板3上绘制一个圆形,所述摆锤4投影落在所述标定板3的圆形内,所述高清视频摄像头6拍摄范围涵盖所述标定板3的边界,所述圆形上设有8个圆形的标志点,任意两个相邻的标志点与圆心的夹角为45°,并且,8个标志点与圆心的距离均相等,所述摆锤4是球形,所述标志点的直径小于所述摆锤4落在所述投影在所述标定板3上投影的直径,便于倾角的计算。如图2和图3所示,小圆点表示标志点,构成多组垂直的直角坐标系,图2中大圆点表示摄像头拍摄摆锤位置,本发明中称为测量点,即根据标志点求解测量点的摆动角度。本发明中为了方便计算,设定测量点的半径是标志点半径的2倍。
本发明还提供一种基于图像修正的高精度倾角测量装置的测量方法,包括如下步骤:
步骤1:所述底座固定板1的底面与待测倾角面紧密接触,
步骤2:高清视频摄像头6进行拍摄,获得标定板图像;
步骤3:进行标定板图像中阴影区域的识别;
步骤4:如图2所示,当阴影区域的个数为9个时,判定为8个标志点和1个测量点;如图3所示,当阴影区域的个数为8个时,判定为8个标志点和1个测量点;当阴影区域的个数既不是8个也不是9个时,判定为无效数据,重复步骤2-4;
步骤5:根据标志点和测量点求解摆锤4的摆动角度。
进一步的,步骤3包括如下步骤:
步骤3.1:将标定板图像中所有像素点存入像素集合A,对于像素集合A中任一像素点i,将像素点i的灰度值g(i)与灰度阈值δ进行比较,当g(i)>δ时,判定g(i)灰度值不满足要求,将g(i)从像素集合A中删除,重新执行步骤3.1,当g(i)≤δ时,将g(i)从像素集合A中删除,并将像素点i存入阴影集合P中,继续进行步骤3.2;
步骤3.2:对于阴影集合P中像素点i,坐标为(x,y),依次遍历与像素点i相邻的8个邻域像素点:(x-1,y),(x-1,y-1),(x-1,y+1),(x,y-1),(x,y+1),(x+1,y),(x+1,y-1),(x+1,y+1),去除超出标定板图像坐标范围的邻域像素点,对剩余的邻域像素点进行灰度判断,如果剩余邻域像素点中任一像素点j的灰度值g(j)≤δ,则将该邻域像素点j存入阴影集合P中,重复遍历阴影集合P,直到阴影集合P中所有邻域像素点的灰度值均>阈值δ;
步骤3.3:重复步骤3.1和3.2,找出标定板图像中所有阴影区域的结合P1,P2,……,Pk,其中,Pk表示第k个阴影集合,k是阴影区域的个数。
进一步的,步骤2高清视频摄像头6进行拍摄时,理论上摄影头光轴与测量平面的夹角为90°,但由于外界因素存在摄像头产生倾斜的情况,使得摄像头光轴与测量平面夹角与理论值不符,此时需进行摄像头光轴倾角校正,如图4所示,点O为相机,OM为理论光轴,垂直于平面ABCD,AB垂直于CD,AB记为x轴,CD记为y轴,
如图5所示,假定OM’为实际光轴,将OM’分别投影到OAB平面和OCD平面得到点M1’和点M2’,OM’分别在OAB平面和OCD平面求解∠MOM1’和∠MOM2’,
以OAB平面的投影为例,图6所示,A’、B’为标志点A和标志点B在摄像头倾斜后图片上的成像,则:
∠MOM′1=arctan(MM′1/OM)=arctan((BM-BM′1)/OM)
同理,
∠MOM′2=arctan(MM′2/OM)=arctan((DM-DM′2)/OM)
其中,BM1’满足下式:
同理,
其中,XA′、XB′分别表示标志点A、B在图片阴影区域的中心点的x轴的坐标,表示点M1’的x轴的坐标,
C’、D’分别表示标志点C、D在图片阴影区域的中心点,YC′、YD′分别表示标志点C、D在图片阴影区域的中心点的y轴的坐标,表示点M2’的y轴的坐标。
进一步的,步骤5进行标志点和测量点的区分具体操作如下:
(1)当k=9时,
首先,分别获取各阴影区域中心点坐标,
然后,分别计算各阴影区域中心点与8个标志点的标定坐标的欧式距离,将获得的9个欧式距离按照从小到大的顺序排列,欧式距离最大的阴影区域为测量点,其余阴影区域为标志点,其中,8个标志点的理论坐标分别为:(rb,0)、(-rb,0)、(0,rb)、(0,-rb)、 和/>其中rb是标志点与圆心的距离;
(2)当k=8时,表示测量点和其中一个标志点阴影区重合,
首先,分别获取各阴影区域中心点坐标,
然后,找出包含像素点最多的阴影区域,为测量点和标志点重合的阴影区域,其余阴影区域为标志点,
接下来,找出与重合阴影区域中的标志点关于圆心相对称的标志点,基于对称标志点的中心坐标得到重合阴影区域中的标志点的中心坐标,
如图7所示,在摄像头产生倾斜的情况,通过计算出的摄像头在OAB和OCD两个平面的夹角,以及与重合阴影区域中的标志点P关于圆心相对称的标志点Q,计算出被重合的标志点的中心点坐标,具体公式如下:
其中,XP′、YP′分别表示重合阴影区域中的标志点P在图片阴影区域的中心点的x轴、y轴的坐标,XQ′、YQ′分别表示标志点P关于圆心相对称的标志点Q在图片阴影区域的中心点的x轴、y轴的坐标,P’表示重合阴影区域中的标志点P对应的摄像头拍摄图片上的实际点,Q’表示标志点P关于圆心相对称的标志点Q摄像头拍摄图片上的实际点,表示M1’和P’之间的像素距离,M′2Q′表示M2’和Q’之间的像素距离,
提取重合阴影区域的边缘轮廓点集,其中,边缘轮廓点集通过依次遍历重合阴影区域内所有的像素点,当任一像素点的上下左右四个邻域的像素点中至少有1个像素点不在阴影区域集合内,则将该像素点放入轮廓点集中,
进一步,计算边缘轮廓点集与被重合标志点实际坐标P’的欧式距离,从轮廓点集中删除欧式距离与标志点半径差值小于距离阈值&的点,&=(rc+rb)/2,其中,rc表示测量点到圆心的半径;
最后,基于剩余轮廓点集拟合圆心,即为测量点的中心。
进一步的,同样将摆锤4的倾角拆分到两个平面进行计算,如图8所示,摆锤4的摆动角度∠Z1OM为摆锤4倾角在OAB平面投影角,具体计算如下:
其中,表示Z’到M1’的像素距离,/>满足:
其中,OM1’满足:
同理可求得在OCD平面的倾角∠Z2OM,具体计算如下:
其中,表示Z’到M2’的像素距离,/>满足:
其中,OM2’满足:
本发明提出一种基于图像修正的高精度倾角测量装置,上述实施例只为说明本发明的技术构思及特点,其目的在于让熟悉此项技术的人士能够了解本发明的内容并据以实施,并不能以此限制本发明的保护范围。凡根据本发明精神实质所作的等效变化或修饰,都应涵盖在本发明的保护范围之内。
Claims (10)
1.一种基于图像修正的高精度倾角测量装置,其特征在于:包括上下表面平行的底座固定板(1),测量时所述底座固定板(1)的底面与待测倾角面紧密接触,所述底座固定板(1)上部依次设有磁场建立单元(2)和标定板(3),所述标定板(3)与所述底座固定板(1)上表面相平行,摆锤(4)悬挂在所述标定板(3)的正上方,所述磁场建立单元(2)用于产生磁场,通过磁场与所述摆锤(4)产生磁感应,基于电磁阻尼原理使所述摆锤(4)的摆动迅速收敛至静止,所述标定板(3)用于对测量误差进行标定,高清视频摄像头(6)设置在所述摆锤(4)的正上方,所述高清视频摄像头(6)用于实时拍摄所述摆锤(4)的摆动幅度,基于所述摆锤(4)在图像中的位置,得到所述摆锤(4)的倾斜角度,待测倾角面的倾角等于所述摆锤(4)的倾斜角度。
2.根据权利要求1所述的一种基于图像修正的高精度倾角测量装置,其特征在于:
所述磁场建立单元(2)、所述标定板(3)、所述摆锤(4)和所述高清视频摄像头(6)均设置在低温度系数的密封罩(5)内,消除外界温度变化对测量的影响,提高了倾角的测量精度。
3.根据权利要求2所述的一种基于图像修正的高精度倾角测量装置,其特征在于:
所述密封罩(5)的顶部设有密封盖(7),所述高清视频摄像头(6)和所述摆锤(4)固定设置在所述密封盖(7)上,便于安装。
4.根据权利要求1所述的一种基于图像修正的高精度倾角测量装置,其特征在于:
所述底座固定板(1)上侧中心位置设置有凹槽,所述磁场建立单元(2)固定设置在凹槽内,便于所述标定板(3)与所述底座固定板(1)连接。
5.基于权利要求1所述的一种基于图像修正的高精度倾角测量装置的测量方法,其特征在于:
所述标定板(3)上绘制一个圆形,所述摆锤(4)投影落在所述标定板(3)的圆形内,所述高清视频摄像头(6)拍摄范围涵盖所述标定板(3)的边界,所述圆形上设有8个圆形的标志点,任意两个相邻的标志点与圆心的夹角为45°,并且,8个标志点与圆心的距离均相等,所述摆锤(4)是球形,所述标志点的直径小于所述摆锤(4)落在所述投影在所述标定板(3)上投影的直径,便于倾角的计算。
6.根据权利要求5所述的的一种基于图像修正的高精度倾角测量方法,其特征在于,包括如下步骤:
步骤1:所述底座固定板(1)的底面与待测倾角面紧密接触,
步骤2:高清视频摄像头(6)进行拍摄,获得标定板图像;
步骤3:进行标定板图像中阴影区域的识别;
步骤4:当阴影区域的个数为9个时,判定为8个标志点和1个测量点;当阴影区域的个数为8个时,判定为8个标志点和1个测量点;当阴影区域的个数既不是8个也不是9个时,判定为无效数据,重复步骤2-4;
步骤5:根据标志点和测量点求解摆锤(4)的摆动角度。
7.根据权利要求6所述的一种基于图像修正的高精度倾角测量方法,其特征在于:
步骤3包括如下步骤:
步骤3.1:将标定板图像中所有像素点存入像素集合A,对于像素集合A中任一像素点i,将像素点i的灰度值g(i)与灰度阈值δ进行比较,当g(i)>δ时,判定g(i)灰度值不满足要求,将i从像素集合A中删除,重新执行步骤3.1,当g(i)≤δ时,将i从像素集合A中删除,并将像素点i存入阴影集合Pn中,n是执行次数,n=1,2,3,……,继续进行步骤3.2;
步骤3.2:对于阴影集合Pn中像素点i,坐标为(x,y),依次遍历与像素点i相邻的8个邻域像素点:(x-1,y),(x-1,y-1),(x-1,y+1),(x,y-1),(x,y+1),(x+1,y),(x+1,y-1),(x+1,y+1),如果邻域像素点的灰度值≤δ,则将其存入集合Pn中,再依次遍历这些邻域像素点的领域像素点,重复下去直到没有邻域像素点为止;
步骤3.3:重复步骤3.1和3.2,找出标定板图像中所有阴影区域的结合(P1,P2,……,Pk),其中,Pk表示第k个阴影集合,k是阴影区域的个数。
8.根据权利要求7所述的一种基于图像修正的高精度倾角测量方法,其特征在于:
步骤2高清视频摄像头(6)进行拍摄时,理论上摄影头光轴与测量平面的夹角为90°,由于外界因素存在摄像头产生倾斜的情况,使得摄像头光轴与测量平面夹角与理论值不符,此时需进行摄像头光轴倾角校正,点O为相机,OM为理论光轴,垂直于平面ABCD,AB垂直于CD,AB作为x轴,CD作为y轴,
假定OM’为实际光轴,将OM’分别投影到OAB平面和OCD平面得到点M1’和点M2’,OM’分别在OAB平面和OCD平面求解∠MOM1’和∠MOM2’,
A’、B’为标志点A和标志点B在摄像头倾斜后图片上的成像,则:
∠MOM′1=arctan(MM′1/OM)=arctan((BM-BM′1)/OM)
同理,
∠MOM′2=arctan(MM′2/OM)=arctan((DM-DM′2)/OM)
其中,BM1’满足下式:
同理,
其中,XA′、XB′分别表示标志点A、B在图片阴影区域的中心点的x轴的坐标,表示点M1’的x轴的坐标,
C’、D’分别表示标志点C、D在图片阴影区域的中心点,YC′、YD′分别表示标志点C、D在图片阴影区域的中心点的y轴的坐标,表示点M2’的y轴的坐标。
9.根据权利要求7所述的一种基于图像修正的高精度倾角测量方法,其特征在于:
步骤5进行标志点和测量点的区分具体操作如下:
(1)当k=9时,
首先,分别获取各阴影区域中心点坐标,
然后,分别计算各阴影区域中心点与8个标志点的标定坐标的欧式距离,将获得的9个欧式距离按照从小到大的顺序排列,欧式距离最大的阴影区域为测量点,其余阴影区域为标志点;
(2)当k=8时,
首先,分别获取各阴影区域中心点坐标,
然后,找出包含像素点最多的阴影区域,为测量点和标志点重合的阴影区域,其余阴影区域为标志点,
接下来,找出与重合阴影区域中的标志点关于圆心相对称的标志点,基于对称标志点的中心坐标得到重合阴影区域中的标志点的中心坐标,
提取重合阴影区域的边缘轮廓点集,其中,边缘轮廓点集通过依次遍历重合阴影区域内所有的像素点,当任一像素点的上下左右四个邻域的像素点中至少有1个像素点不在阴影区域集合内,则将该像素点放入轮廓点集中,
进一步,计算边缘轮廓点集与被重合标志点实际坐标P’的欧式距离,从轮廓点集中删除欧式距离与标志点半径差值小于距离阈值的点;
最后,基于剩余轮廓点集拟合圆心,即为测量点的中心坐标。
10.根据权利要求9所述的一种基于图像修正的高精度倾角测量方法,其特征在于:
摆锤(4)摆动倾角在OAB平面投影角∠Z1OM,具体计算如下:
其中,Z′表示Z在图片阴影区域的中心点,表示Z’到M1’的像素距离,/>满足:
其中,XZ′表示点Z1’的x轴的坐标,OM1’满足:
摆锤(4)摆动倾角在OCD平面投影角∠Z2OM,具体计算如下:
其中,表示Z’到M2’的像素距离,/>满足:
其中,YZ′表示点Z1’的y轴的坐标,OM2’满足:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310873057.1A CN116753914B (zh) | 2023-07-17 | 2023-07-17 | 一种基于图像修正的高精度倾角测量装置及测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310873057.1A CN116753914B (zh) | 2023-07-17 | 2023-07-17 | 一种基于图像修正的高精度倾角测量装置及测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116753914A true CN116753914A (zh) | 2023-09-15 |
CN116753914B CN116753914B (zh) | 2024-03-29 |
Family
ID=87957255
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310873057.1A Active CN116753914B (zh) | 2023-07-17 | 2023-07-17 | 一种基于图像修正的高精度倾角测量装置及测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116753914B (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110093238A1 (en) * | 2008-05-26 | 2011-04-21 | Roberto Biasi | Electromechanical device for measuring the inclination of a support plane with high resolution, high accuracy and low sensitivity to outside disturbances |
CN105466394A (zh) * | 2016-01-05 | 2016-04-06 | 上海筑邦测控科技有限公司 | 一种基于吊锤位置视频识别技术的倾角传感器 |
CN109612430A (zh) * | 2018-11-29 | 2019-04-12 | 昆山睿力得软件技术有限公司 | 一种基于视觉引导的经纬仪测量方法 |
-
2023
- 2023-07-17 CN CN202310873057.1A patent/CN116753914B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110093238A1 (en) * | 2008-05-26 | 2011-04-21 | Roberto Biasi | Electromechanical device for measuring the inclination of a support plane with high resolution, high accuracy and low sensitivity to outside disturbances |
CN105466394A (zh) * | 2016-01-05 | 2016-04-06 | 上海筑邦测控科技有限公司 | 一种基于吊锤位置视频识别技术的倾角传感器 |
CN109612430A (zh) * | 2018-11-29 | 2019-04-12 | 昆山睿力得软件技术有限公司 | 一种基于视觉引导的经纬仪测量方法 |
Also Published As
Publication number | Publication date |
---|---|
CN116753914B (zh) | 2024-03-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107256568B (zh) | 一种高精度机械臂手眼相机标定方法及标定系统 | |
CN112629431B (zh) | 土木结构变形监测方法及相关设备 | |
CN106971408B (zh) | 一种基于时空转换思想的摄像机标定方法 | |
CN101539397B (zh) | 物体三维姿态的精密光学测量方法 | |
CN101694370B (zh) | 大尺寸工业摄影测量系统的空间误差场获取方法及基准装置 | |
CN110285770B (zh) | 一种桥梁挠度变化测量方法、装置及设备 | |
CN109887041B (zh) | 一种机械臂控制数字相机摄影中心位置和姿态的方法 | |
CN111189403B (zh) | 一种隧道变形的监测方法、装置及计算机可读存储介质 | |
CN109883444A (zh) | 一种姿态角耦合误差补偿方法、装置及电子设备 | |
CN109855822A (zh) | 一种基于无人机的高铁桥梁竖向动扰度测量方法 | |
CN109798874A (zh) | 一种高铁桥梁竖向动扰度测量方法 | |
CN109087360A (zh) | 一种机器人相机外参的标定方法 | |
CN105627916B (zh) | 一种建立跟踪仪地理坐标系及六自由度测量的方法 | |
CN106292135B (zh) | Tdi ccd相机焦面映射至三轴气浮转台球心的方法 | |
CN103364167B (zh) | 一种观察窗折射偏移校正方法 | |
CN115222819A (zh) | 一种基于机场大范围场景下多模信息参照的相机自标定和目标追踪方法 | |
CN116753914B (zh) | 一种基于图像修正的高精度倾角测量装置及测量方法 | |
WO2022126339A1 (zh) | 土木结构变形监测方法及相关设备 | |
CN109813510A (zh) | 基于无人机的高铁桥梁竖向动扰度测量方法 | |
CN107976146B (zh) | 一种线阵ccd相机的自标定方法及测量方法 | |
CN111750896A (zh) | 云台标定方法、装置、电子设备及存储介质 | |
CN111047552B (zh) | 三维挠度测量方法及相关产品 | |
CN108548506A (zh) | 一种利用光学标志对高精度平面进行平面度测量的方法 | |
CN115046498B (zh) | 一种单目旋转结构光三维测量系统标定方法 | |
CN111561867A (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 |