CN1224938C - 由单幅图像测量场景及场景内物体几何信息的方法 - Google Patents
由单幅图像测量场景及场景内物体几何信息的方法 Download PDFInfo
- Publication number
- CN1224938C CN1224938C CN 02118167 CN02118167A CN1224938C CN 1224938 C CN1224938 C CN 1224938C CN 02118167 CN02118167 CN 02118167 CN 02118167 A CN02118167 A CN 02118167A CN 1224938 C CN1224938 C CN 1224938C
- Authority
- CN
- China
- Prior art keywords
- scene
- point
- image
- information
- straight line
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 43
- 239000011159 matrix material Substances 0.000 claims abstract description 20
- 238000013507 mapping Methods 0.000 claims abstract description 11
- 238000005259 measurement Methods 0.000 claims description 25
- 239000000284 extract Substances 0.000 claims description 18
- 230000009466 transformation Effects 0.000 claims description 11
- 238000004364 calculation method Methods 0.000 claims description 6
- 230000005540 biological transmission Effects 0.000 claims description 4
- 230000008676 import Effects 0.000 claims description 2
- 238000003672 processing method Methods 0.000 claims description 2
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 description 4
- 238000000605 extraction Methods 0.000 description 4
- 238000006243 chemical reaction Methods 0.000 description 2
- 239000004744 fabric Substances 0.000 description 2
- 230000003993 interaction Effects 0.000 description 2
- 229910052742 iron Inorganic materials 0.000 description 2
- 238000013076 uncertainty analysis Methods 0.000 description 2
- 238000007476 Maximum Likelihood Methods 0.000 description 1
- 239000003795 chemical substances by application Substances 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000007598 dipping method Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 230000008569 process Effects 0.000 description 1
Images
Landscapes
- Image Analysis (AREA)
- Image Processing (AREA)
Abstract
一种由单幅图像测量场景及场景内物体几何信息的方法,包括步骤:拍摄一幅或多幅被测场景的图像并输入计算机;从图像中提取出参考平面内的点或直线的信息,垂直平行线的信息及参考距离的端点;利用参考平面内的点或直线的信息及其与场景的对应关系,求出参考平面到图像的映射关系矩阵;利用垂直平行线的信息及映射关系矩阵计算出图像上的消隐点及参考平面的消隐线;然后利用这些信息进行测量。本发明可以进行下述测量:点在场景坐标系中的坐标及两点之间的距离、物体的高度、点到直线的距离、两直线间的夹角、平面图型的面积以及规则物体的体积等。本发明提供的测量方法测量精度高、鲁棒性好,使基于视觉的测量方法走向实用。
Description
技术领域
本发明涉及计算机辅助自动测量。
背景技术
在我们日常生活和许多领域中,测量是必不可少的,比如交通事故、刑事犯罪现场,建筑物、施工现场等,都需要进行大量的测量工作。传统方法都是采用手工测量,这种方法十分耗时,容易出现较大的测量误差,而且这种测量也容易因测量人员的主观因素产生不公正的结果。
近年来,随着计算机视觉的发展和应用,人们开始探索借助计算机视觉方法,通过图像进行测量。在这类方法中,较为典型的是通过两幅或多幅图像利用立体视觉的方法重建出被测场景的三维结构,但该方法涉及到计算机视觉领域的经典难题,即多幅图像间的匹配,因而很难投入实际应用。另一种方法是事先在场景平面中人工测量4个以上控制点的信息并以此来计算空间平面与像平面间的投影矩阵,然后通过单幅图像进行测量,这种方法的测量精度与控制点的选取精度、位置和数量直接相关,但在实际应用中,有时控制点很难选取,而且也无法保证控制点的测量精度。
发明内容
本发明的目的在于提供一种较为简便、实用并具有较高测量精度和鲁棒性的基于单幅图像的测量方法。
为实现上述目的,由单幅图像测量场景及场景内物体几何信息的方法,包括步骤:
1).拍摄一幅或多幅被测场景的图像并输入计算机;
2).从图像中提取出场景中参考平面内的点或直线的信息,与参考平面垂直的平行线信息以及参考距离的两个端点;
3).利用步骤2提取出的参考平面内点或直线的信息及其与场景的对应关系,求出参考平面到图像平面的映射关系矩阵;
4).利用步骤2提取出的点或线的信息,在图像上计算出垂直于参考平面方向上的消隐点及参考平面的消隐线。
5).通过图像处理方法从图像中提取出需测量的物体的边界信息并将其分段拟合成直线或二次曲线,对于物体的上下表面,利用图像上的对极关系或物体表面自身的性质求出该表面中心点的坐标;
6).利用上述步骤求出的参数测量场景及场景内物体的几何信息。这些信息包括:空间点的坐标、距离、直线夹角、点到直线的距离、平面图形面积及物体体积。
本发明可以进行以下测量:参考平面内任意一点在场景坐标系中的坐标及任意两点之间的距离;空间中任意一点相对于参考平面的高度、在场景坐标系中的坐标及任意两点之间的距离;点到参考平面内某一直线的距离;场景中参考平面及与参考平面相平行的平面内任意两直线间的夹角、平面图形的面积以及一些相对规则和对称的物体的体积,如方体、圆柱体、圆锥体、圆台等。同时,可以根据对各步骤的不确定性分析给出每个测量结果的误差范围。本发明提供的测量方法不需要对摄像机参数进行标定,具有较为简便、实用、测量精度高、鲁棒性好等特点,使基于视觉的测量方法走向实用。
附图说明
图1是本测量方法的流程图;
图2参考平面外点的坐标及高度的测量原理图
图3参考平面内点到直线距离的测量原理图
图4是实施例的测量场景图。
具体实施方式
本发明方法要求场景中含有一张参考平面及与该参考平面垂直的两条或多条平行线及一个参考距离,在参考平面内应包含一些已知的几何信息,如一些点或直线的信息,这些信息可以通过场景中的窗框、门框、地板的纹理、规则物体的轮廓等获得,也可以通过在场景中放置一个已知几何信息的模板获得,同时,场景中应具有一个参考距离。用普通照相机或数码相机拍摄一幅图像即可进行测量,而且不需要对摄像机参数进行标定,若能从不同角度拍摄多幅图像,则可从中选取一幅较好的进行测量或综合利用同一被测量在多幅图像中的测量值以便提高测量精度。
对于一幅图像,该方法主要包括图像获取、关键特征提取、变换矩阵计算、消隐线及消隐点计算、物体特征的提取和计算、测量等步骤组成,如图1所示,各步的具体说明如下:
1.获取图像
用普通照相机或数码相机拍摄一幅或多幅被测场景的图像,并借助扫描仪或某些专用接口将图像输入计算机。
2.关键特征提取
通过自动或人机交互方式从图像中提取出场景中参考平面内的点或直线的信息,与参考平面垂直的平行线信息,参考距离的两个端点。如果在场景中放置有模板,则从图像中提取模板中所包含的已知点或直线的几何信息。其中,可以用Canny或Robert、Sobel等边缘检测算子提取边缘信息,然后用Hough变换或最小二乘、极大似然估计等方法求取直线参数。
3.变换矩阵计算
利用第2步提取出的参考平面内点或直线的信息及其与场景的对应关系,求出参考平面到图像平面的映射关系矩阵H,其中,H为一3×3的矩阵。以点的对应为例,设场景中点的齐次坐标为Xi,图像中的对应点为mi,则由mi=sHXi即可求出映射关系矩阵H。
4.消隐线及消隐点计算
利用第2步提取出的与参考平面垂直的平行线信息,计算出图像在该方向上的消隐点Vr;利用第3步所得到的变换矩阵计算出图像中参考平面的消隐线Lv,或直接利用第2步提取出的参考平面内的信息计算消隐线。
5.物体特征的提取和计算
采用与第2步相类似的方法,通过自动或人机交互方法从图像中提取出需测量的物体的边界信息并将其分段拟合成直线或二次曲线。对于物体的上下表面,利用图像上的对极关系或物体表面自身的性质求出该表面中心点的坐标。
6.测量
本发明可以进行以下测量:参考平面内任意一点在场景坐标系中的坐标及任意两点之间的距离;空间中任意一点相对于参考平面的高度、在场景坐标系中的坐标及任意两点之间的距离,测量的前提是该点在参考平面内的垂点是已知的;场景中参考平面及与参考平面相平行的平面内任意两直线间的夹角;点到参考平面内直线距离;平面图形的面积以及一些相对规则和对称的物体的体积等。同时,可以根据对各步骤的不确定性分析给出每个测量结果的误差范围。具体测量方法如下:
(1).参考平面内点的坐标及距离的测量:对于参考平面图像上的任意一点m,利用第3步求出的映射关系矩阵的逆变换就可以求出场景中对应点的坐标值X=H-1m,进而可以算出参考平面内任意两点间的距离。
(2).参考平面外点的坐标及高度的测量:如图2所示,Lv和Vr分别为第4步求出的消隐线和消隐点;S0和E0为参考距离的两个端点,则S0、E0和Vr三点共线Lr;S1和E1为待测高度的两个端点,则借助消隐线Lv可将点S1映射到直线Lr上与其具有等高度的点H1上。这样,根据射影变换保持交比不变的性质即可求出S1和E1对应的空间距离。另一方面,由于点E1在参考平面内,其坐标可由(1)得到,由此可以得到点S1对应的空间点的坐标。
(3).场景中任意两点之间距离的测量:由(1)和(2)可以得到参考平面内任意两点间的距离以及空间点的高度,通过简单的几何代数运算即可求得场景中任意两点之间的距离。
(4).参考平面内任意两条直线间夹角的测量:设两条直线与消隐线Lv的两个交点为P1、P2,由第3步求出的映射关系矩阵H可以得到两个圆环点的图像I、J,则P1、P2和I、J四点共线,根据这四点的交比及射影变换保持交比不变的性质即可求出两条直线间的夹角。例如:设此交比为Cr,则夹角θ=ln(Cr)/(2*i)。也可以通过变换矩阵将此直线映射到空间平面上,根据两条直线的方程直接计算其夹角。
(5).参考平面及与参考平面相平行的平面内任意平面图形面积的测量:对于一般的平面图形,在图像上检测出其边缘信息并通过映射关系矩阵H将其映射到空间平面上,然后通过积分法求出其面积的近似值;对于矩形和三角形,可由(1)直接求出其边长并利用矩形和三角形面积公式计算面积;对于圆,在图像上提取其边缘并用二次曲线拟合出其方程,利用二次曲线与消隐线间的对极关系求出其中心坐标,由(1)计算该中心与任一边缘点之间的距离即为半径,从而可以利用圆的面积公式得到其面积。
(6).物体体积的测量:由(5)可以测出物体的底面积及截面积,由(2)可以测出物体的高度,这样,通过简单的代数或几何运算即可以由体积计算公式求出物体的体积。例如:对于圆锥体,由第5步提取出底部的边界信息并拟合成二次曲线并求出其中心点的坐标,这样就可以计算出底面积和高度从而求出体积;对于圆台,分别提取出上下表面的边界信息并拟合成二次曲线,计算出上下表面的面积,求出其中心点的坐标中心,则两中心点之间的距离即为高度,从而可以利用体积计算公式求出体积。
(7).点到参考平面内直线距离的测量:测量方法如图3所示,求参考平面内的点A到直线BC的距离。设BC与消隐线Lv的交点为D,LV上的两个圆环点为I、J,在Lv上计算与D、I、J交比为-1的点E,连接A、E与BC交于点H,由(1)计算A、H间的距离即为所求的点A到直线BC的距离。若要计算空间一点S到直线BC的距离,可先由(2)测出S点的高度,再计算S点在参考平面内的垂点到直线BC的距离,则通过简单的代数或几何运算即可求出S点到直线BC的距离。也可以通过变换矩阵将此直线映射到空间平面上,由(2)计算出点在场景坐标系中的坐标,然后根据直线的方程和空间点的坐标直接计算点到直线的距离。
实施例
如图4所示为一被测场景的图像,我们利用边缘检测和直线拟合算法提取出图中已知的地板信息和垂直于地板的平行线的信息,通过最小二乘法计算出映射关系矩阵、消隐线及消隐点,参考距离为铁柜的高度H0。按照本发明所述的方法,在图像中用鼠标任选两点,就可以直接从图像求出场景中这两点间的距离,测量结果如图4所示。其中,V为圆铁桶的体积,可以看出,距离D1和D2在现实空间中是很难测量的。通过本发明方法所得到的测量值如图4所示。图中各测量值的真实值如下:人的高度H1=171.5cm;控制柜高度H2=170.4cm;距离D1=302.7cm;D2=389.2cm;圆铁桶的直径为24cm,高为59.3cm,体积为V=26826.7cm3。本发明的测量相对误差一般在4%以内,这是完全可以达到测量精度要求的。
Claims (12)
1.一种由单幅图像测量场景及场景内物体几何信息的方法,包括步骤:
1)拍摄一幅或多幅被测场景的图像并输入计算机;
2)从图像中提取出场景中参考平面内的点或直线的信息,与参考平面垂直的平行线信息以及参考距离的两个端点;
3)利用步骤2提取出的参考平面内点或直线的信息及其与场景的对应关系,求出参考平面到图像平面的映射关系矩阵;
4)利用步骤2提取出的点或线的信息,在图像上计算出垂直于参考平面方向上的消隐点及参考平面的消隐线;
5)通过图像处理方法从图像中提取出需测量的物体的边界信息并将其分段拟合成直线或二次曲线,对于物体的上下表面,利用图像上的对极关系或物体表面自身的性质求出该表面中心点的坐标;
6)利用上述步骤求出的参数测量场景及场景内物体的几何信息,这些信息包括:空间点的坐标、距离、直线夹角、点到直线的距离、平面图形面积及物体体积。
2.按权利要求1所述的由单幅图像测量场景及场景内物体几何信息的方法,其特征在于在场景中放置模板,利用模板上已知的点、线或图案的几何信息计算映射关系矩阵、消隐点及消隐线。
3.按权利要求1所述的由单幅图像测量场景及场景内物体几何信息的方法,其特征在于利用待测场景中参考平面及垂直于
该平面的直线的
消隐 点进行测量。
4.按权利要求1所述的由单幅图像测量场景及场景内物体几何信息的方法,其特征在于将待测高度的点映射到参考距离所在的直线上,根据射影变换保持交比不变的性质即可求出其对应的空间高度。
5.按权利要求1所述的由单幅图像测量场景及场景内物体几何信息的方法,其特征在于利用参考平面内任意两点间的距离以及空间点的高度,求得场景中任意两点之间的距离。
6.按权利要求1所述的由单幅图像测量场景及场景内物体几何信息的方法,其特征在于根据射影变换保持交比不变的性质计算两条直线间的夹角。
7.按权利要求1所述的由单幅图像测量场景及场景内物体几何信息的方法,其特征在于通过变换矩阵将直线映射到空间平面上,根据两条直线的方程直接计算其夹角。
8.按权利要求1所述的由单幅图像测量场景及场景内物体几何信息的方法,其特征在于对于位于参考平面或平行于参考平面的平面内的平面图形,在图像上检测出其边缘信息并通过映射关系矩阵将其映射到空间平面上,然后通过积分法求出其面积的近似值。
9.按权利要求1所述的由单幅图像测量场景及场景内物体几何信息的方法,其特征在于利用与两个圆环点交比为-1的方法求取点到参考平面内直线的垂足从而计算点到直线的距离。
10.按权利要求1或4所述的由单幅图像测量场景及场景内物体几何信息的方法,其特征在于通过变换矩阵将直线映射到空间平面上,根据直线的方程和空间点的坐标直接计算点到直线的距离。
11.按权利要求1所述的由单幅图像测量场景及场景内物体几何信息的方法,其特征在于利用二次曲线与消隐线间的对极关系求出圆的中心坐标。
12.按权利要求1或11所述的由单幅图像测量场景及场景内物体几何信息的方法,其特征在于利用圆所对应的图像二次曲线的中心坐标计算圆柱、圆锥、圆台等物体的高度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 02118167 CN1224938C (zh) | 2002-04-23 | 2002-04-23 | 由单幅图像测量场景及场景内物体几何信息的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 02118167 CN1224938C (zh) | 2002-04-23 | 2002-04-23 | 由单幅图像测量场景及场景内物体几何信息的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1453740A CN1453740A (zh) | 2003-11-05 |
CN1224938C true CN1224938C (zh) | 2005-10-26 |
Family
ID=29257296
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 02118167 Expired - Fee Related CN1224938C (zh) | 2002-04-23 | 2002-04-23 | 由单幅图像测量场景及场景内物体几何信息的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN1224938C (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101226637B (zh) * | 2007-01-18 | 2010-05-19 | 中国科学院自动化研究所 | 一种车轮与地面接触点的自动检测方法 |
Families Citing this family (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN100344936C (zh) * | 2004-12-25 | 2007-10-24 | 鸿富锦精密工业(深圳)有限公司 | 影像几何元素自动判读***及方法 |
CN101625723B (zh) * | 2009-07-02 | 2012-03-28 | 浙江省电力公司 | 电力线轮廓的快速图像识别方法 |
CN102853820A (zh) * | 2011-07-01 | 2013-01-02 | 中国钢铁股份有限公司 | 高炉的落料轨迹的量测方法 |
CN103177254A (zh) * | 2011-12-26 | 2013-06-26 | 鸿富锦精密工业(深圳)有限公司 | 量测元素提取***及方法 |
CN102565303B (zh) * | 2011-12-28 | 2014-10-15 | 中国科学院水利部成都山地灾害与环境研究所 | 一种冲沟沟头溯源侵蚀速率的快速监测方法 |
MX2018005143A (es) * | 2015-10-30 | 2018-06-06 | Unilever Nv | Medicion de rizo del cabello. |
CN106871764A (zh) * | 2017-02-28 | 2017-06-20 | 成都金盘电子科大多媒体技术有限公司 | 一种基于医学影像的空间长度测量方法 |
CN109945831B (zh) * | 2019-03-28 | 2021-03-30 | 东南大学 | 一种简易的建筑物高度测量方法及装置 |
CN110070097B (zh) * | 2019-04-19 | 2023-07-25 | 戴文跃 | 一种图形对象比对方法 |
CN110176003A (zh) * | 2019-06-06 | 2019-08-27 | 三峡大学 | 基于交比图像量测技术检测建筑物表面损伤面积的方法 |
CN112950740B (zh) * | 2019-12-10 | 2024-07-12 | 中交宇科(北京)空间信息技术有限公司 | 高精地图道路中心线的生成方法、装置、设备及存储介质 |
CN113129363A (zh) * | 2019-12-30 | 2021-07-16 | 南京智能创研科技有限公司 | 一种基于特征物体和透视变换的图像距离信息提取方法 |
CN113781422A (zh) * | 2021-09-01 | 2021-12-10 | 廊坊中油朗威工程项目管理有限公司 | 基于单幅图像几何测量算法的管道施工违章识别方法 |
-
2002
- 2002-04-23 CN CN 02118167 patent/CN1224938C/zh not_active Expired - Fee Related
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101226637B (zh) * | 2007-01-18 | 2010-05-19 | 中国科学院自动化研究所 | 一种车轮与地面接触点的自动检测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN1453740A (zh) | 2003-11-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111473739B (zh) | 一种基于视频监控的隧道塌方区围岩变形实时监测方法 | |
CN1224938C (zh) | 由单幅图像测量场景及场景内物体几何信息的方法 | |
CN102103747B (zh) | 采用参照物高度的监控摄像机外部参数标定方法 | |
CN114417489B (zh) | 一种基于实景三维模型的建筑基底轮廓精细化提取方法 | |
CN113587835A (zh) | 一种利用三维激光扫描技术进行桥梁工程质量验收的方法 | |
CN110221603A (zh) | 一种基于激光雷达多帧点云融合的远距离障碍物检测方法 | |
CN111553292A (zh) | 一种基于点云数据的岩体结构面识别与产状分类方法 | |
CN103714571A (zh) | 一种基于摄影测量的单像机三维重建方法 | |
ITMI20090941A1 (it) | Metodo di restituzione fotogrammetrica assistita da nuvola di punti e relativo apparato. | |
CN103322911A (zh) | 一种基于图像的特色建筑立面图测绘方法 | |
Gruen | TOBAGO—a semi-automated approach for the generation of 3-D building models | |
KR101078238B1 (ko) | 항공 레이저 측량의 점군 집합을 활용한 하천 지역 내 3차원 인공 제방선 추출 방법 | |
CN109146990B (zh) | 一种建筑轮廓的计算方法 | |
CN102012213A (zh) | 单副图像测量前景高度的方法 | |
Jiang et al. | Determination of construction site elevations using drone technology | |
Manandhar et al. | Feature extraction from range data | |
Zhang et al. | Automatic terrain extraction using multiple image pair and back matching | |
Ioannidis et al. | 5D Multi-Purpose Land Information System. | |
CN113744393B (zh) | 一种多层级边坡滑坡变化监测方法 | |
CN113469949B (zh) | 一种基于地理信息的输电通道内隐患测距的方法 | |
Grün et al. | Optical 3-d measurement techniques | |
CN113375644A (zh) | 一种边坡产状高精度测量方法 | |
Henricsson et al. | Automated 3-D reconstruction of buildings and visualization of city models | |
CN1284961C (zh) | 基于单幅图像的平面测量方法 | |
CN1288420C (zh) | 基于平行线的单视平面测量方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20051026 Termination date: 20120423 |