CN104778720A - 一种基于空间不变特性的快速体积测量方法 - Google Patents
一种基于空间不变特性的快速体积测量方法 Download PDFInfo
- Publication number
- CN104778720A CN104778720A CN201510231222.9A CN201510231222A CN104778720A CN 104778720 A CN104778720 A CN 104778720A CN 201510231222 A CN201510231222 A CN 201510231222A CN 104778720 A CN104778720 A CN 104778720A
- Authority
- CN
- China
- Prior art keywords
- prime
- point
- volume
- dimensional
- coordinate system
- 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
- 238000000691 measurement method Methods 0.000 title abstract description 3
- 238000000034 method Methods 0.000 claims abstract description 41
- 238000004364 calculation method Methods 0.000 claims abstract description 10
- 239000011159 matrix material Substances 0.000 claims description 38
- 238000005316 response function Methods 0.000 claims description 6
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 4
- 238000012360 testing method Methods 0.000 claims description 4
- 230000001186 cumulative effect Effects 0.000 claims description 3
- 238000011156 evaluation Methods 0.000 claims description 3
- 238000000605 extraction Methods 0.000 claims description 3
- 238000013519 translation Methods 0.000 claims description 3
- 238000005259 measurement Methods 0.000 abstract description 11
- 230000001788 irregular Effects 0.000 abstract description 8
- 230000002093 peripheral effect Effects 0.000 abstract description 2
- 238000003708 edge detection Methods 0.000 abstract 1
- 230000008569 process Effects 0.000 description 8
- 238000005516 engineering process Methods 0.000 description 6
- 239000002689 soil Substances 0.000 description 5
- 230000010354 integration Effects 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 238000011161 development Methods 0.000 description 3
- 238000005457 optimization Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 238000011835 investigation Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 238000004064 recycling Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 231100001261 hazardous Toxicity 0.000 description 1
- 238000004321 preservation Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- VMXUWOKSQNHOCA-UKTHLTGXSA-N ranitidine Chemical compound [O-][N+](=O)\C=C(/NC)NCCSCC1=CC=C(CN(C)C)O1 VMXUWOKSQNHOCA-UKTHLTGXSA-N 0.000 description 1
- 230000000452 restraining effect Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 239000004575 stone Substances 0.000 description 1
- 238000012876 topography Methods 0.000 description 1
Landscapes
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明提供了一种基于空间不变特性的快速体积测量方法,所述方法包括:对所获取的图像进行特征点提取与匹配、三维重建;在三维点云中加入特征点间实际距离与在相对坐标系内与真实世界的距离比例,作为尺度信息;利用特征点三维坐标,构建Delaunay三角网格,加入边界边检测元素,拟合最优基准面;根据离散积分的思想,将三角网格投影到基准面上,求出三角网所包围的不规则物体的体积。即可利用“一个相机,一把尺”完成物体所有特征点间实际距离的测量,省去了实地测量的麻烦;同时在构建Delaunay三角网格时添加了检测边界边的数据结构,可以有效的得到所构成边界曲面的外围基准面,进一步用于体积计算。
Description
技术领域
本发明涉及图像处理领域,特别涉及土地调查中利用近景摄影技术实现不规则体体积自动求取的方法。
背景技术
土地是人类赖以生存和发展的宝贵自然资源,土地的持续利用是人类社会持续发展的基石。在人口日益膨胀的今天,尤其是在我国,土地承受着巨大的压力,人地关系日趋紧张。经过“十五”、“十一五”的努力,“天上看、地上查、网上管”等现代高科技手段在土地监管以及土地测量领域已经全方位应用,航空、航天遥感技术和GPS全球定位系统等空间技术在土地管理中也已经广泛应用,大大提高了土地信息获取的效率和准确性,基本实现土地基础数据的充分共享,以及目标的精确测量,为社会经济发展决策提供服务。
在现有的土地调查技术中,“天上看”主要是以卫星、无人机遥感为主,尽管在大范围测量上有不可替代的作用,但是也存在先天的不足。其中,卫星遥感测量,周期长,时效性差,精度不能达到较高的标准。无人机遥感测绘,成本高,在高程变化复杂的区域,精度难以得到保障。“地上查”主要以GNSS/全站仪或是皮尺测量为主,这两种测量方法都是属于与物体实际坐标有关的绝对测量。基于GNSS的测量方法流程复杂,对工作人员要求较高,布点受地形影响较大,过度较慢,非常不适用于自动化测量及监测;而使用皮尺测量,有大量的外业工作,只适用于小区域测量,并且对于不规则体测量、人类不能到达或较危险区域,它的适用性较差。因此利用近景摄影技术进行土地测量,能够为土地测量提供一条新的快速采集土地信息的新途径,能够广泛应用于三维山体测量、滑坡监测、沙坑容积测量等许多不规则体测量计算的场合。它能构实现目标的精确测量,具有速度较快、精度较高而均匀、成本低、不损坏原地物、不受气候及季节的限制等优点,在土地测量中的应用具有积极的现实意义。
客观世界在空间上是三维的,而现有的图像采集装置所获取的图像是二维的。尽管图像中含有某些形式的三维空间信息,这些信息包括物体边与边之间的几何关系、两幅图像的视差关系、两幅图像中特征点的对应关系以及物体轮廓信息等等,但要真正在计算机中使用这些信息进行进一步的应用处理,就必须采用三维重建技术从二维图像中合理地提取并表达这些三维信息。结合基于图像的三维重建技术在土地测量中的应用具有积极的现实意义,它有着传统手段无法比拟的优势,有着诸多便利之处。它可将大量的外业测量工作移到室内完成,它具有速度较快、精度较高而均匀、成本低、不损坏原地物、不受气候及季节的限制等优点。新型体积测量技术在土地测量中的应用广泛,能实现精确测量的目标。由于基于图像的三维重建技术所需要的成本低,具有很大的灵活性,又能达到简单的三维建模功能,因此,在需要真实感建模的场合,基于图像的建模无疑具有很高的实用价值。
发明内容
发明目的:针对上述现有技术,提出一种基于空间不变特性的快速体积测量方法,实现快速获取不规则物体体积。
技术方案:一种基于空间不变特性的快速体积测量方法,包括如下步骤:
步骤(1),采用单台摄像机获取物体的多幅图像序列,并对所述摄像机进行标定;
步骤(2),基于双目立体视觉技术,提取所述多幅图像上的特征点后,对获取的多张图像序列进行立体匹配,找出任意两幅图像中同名像素点的对应坐标关系;
步骤(3),任意选取所述图像序列中的两幅图像,以拍摄其中一幅图像的摄像机为坐标原点建立相对坐标系,在所述相对坐标系中对所述步骤(2)提取的多幅图像序列上的所有特征点进行三维重建,得到整个物体的三维点云结构;
步骤(4),在所述整个物体的三维点云结构中,选择并测量两个特征点所对应物体上的实际距离,计算得到所述实际距离与这两个特征点在所述相对坐标系中的距离的比值,利用所述比值计算得到三维点云结构中所有特征点间的实际距离;
步骤(5),依据重建的三维点云结构建立物体的Delaunay三角网模型,根据所述Delaunay三角网模型获取计算物体体积的基准面,根据所述基准面以及所有特征点间的实际距离计算物体的体积。
作为本发明的优选方案,依据三角网生长算法建立所述Delaunay三角网模型,并在构建所述Delaunay三角网模型时加入检测三角网边界边的数据成员useCount,所述数据成员useCount用于记录构建的每个三角形的每条边作为基础边来构建新三角形的次数;完成所有三角形建立步骤后,选择数据成员useCount为1的所有边作为Delaunay三角网模型的边界边,再将由所有边界边构成的封闭图形拟合成一个光滑平面,所述光滑平面作为计算物体体积的基准面。
作为本发明的优选方案,所述步骤(2)具体为:
(21),对获取的图像进行预处理,所述预处理包括图片平滑、锐化处理;
(22),特征点提取:对任意两幅图像,先通过一个角点响应函数CRF去判断各个待测点是否为FAST角点,所述CRF函数如下:
N=ΣfCRF(I(p),I(x)) (1)
式(1)中,N为响应函数CRF值的和,即FAST角点的个数;I(x)为待测点周围的任意一点的灰度值;I(p)表示当前待测点的灰度值;p表示当前待测点;式(2)中,εd为一个中心待测点与圆周任一点像素灰度差值的阈值;
设定阈值εd,提取出超过阈值Nf数目的FAST角点;然后利用Harris角点的评价函数找到前Nf个较好的FAST角点;最后使用金字塔算法得到多尺度的待测图像,确定出最终的FAST角点;
(23),采用BRIEF描述子对提取的FAST角点进行描述;
(24),利用ORB特征点匹配算法对选取的FAST角点进行立体匹配与解算,得到任意两幅图片特征点匹配的基础矩阵F。
作为本发明的优选方案,所述步骤(3)具体为:
(31),任意选取两幅图像,以拍摄其中一幅图像的摄像机为坐标原点,以它的摄像机坐标系作为相对坐标系,计算得到拍摄这两幅图像的摄像机的位置关系,即摄像机的外部参数;
(32),结合摄像机标定时得到的摄像机内部参数,根据步骤(2)所得的这两幅图像的对应基础矩阵F,计算得到其本质矩阵E;
(33),对本质矩阵E进行奇异值分解,得到摄像机外部参数(R|t)的候选值,建立拍摄这两幅图像的两个不同位置的相机的投影矩阵P1和P2:
P1=K[I O] (3)
P2=K[R t] (4)
其中I是3x3的单位矩阵;O是3x1的全零矩阵;K为摄像机内部参数矩阵;R为3x3的旋转矩阵;t为三维平移向量;R、t均为相机外部参数;
(34),根据恢复出的投影矩阵P1和P2,在相对坐标系中利用SFM算法进行迭代计算,重建出这两幅图像上特征点对应的三维点云结构,即得到特征点在相对坐标系内的坐标信息;
(35),在所述步骤(31)建立的相对坐标系下,按照步骤(31)至步骤(34),对其任意两幅图像做三维重建,得到整个物体在该相对坐标系下的三维点云结构。
作为本发明的优选方案,所述步骤(5)中,根据所述基准面以及所有特征点间的实际距离计算物体的体积具体步骤为:
(1),设Delaunay三角网模型中任意三角形的三个顶点坐标为:A(x1,y1,z1),B(x2,y2,z2),C(x3,y3,z3),使用向量法进行投影计算任意三角形的三个顶点投射到所述基准面上的投影点坐标A1(x1',y1',z1')、B1(x'2,y'2,z'2)、C1(x'3,y'3,z'3);
(2),计算Delaunay三角网模型中任意三角形与所述基准面形成的单个三棱柱整体体积:
21),计算底面△A1B1C1的面积:
设
底面△A1B1C1的面积
22),采用任意三角形的三个顶点的高程平均值来计算单个三棱柱的高h:
23),计算单个三棱柱的整体体积V为:
(3),将Delaunay三角网模型中所有的n个三角形对应生成的三棱柱的体积相加合,即Vi为式(10)得到的单个三棱柱的整体体积;
(4),将计算出体积V乘以步骤(4)得到的尺度信息的立方,得到物体的总体积。
有益效果:本发明的一种基于空间不变特性的快速体积测量技术,在双目立体视觉技术的基础上,仅利用“一个相机、一把尺”,借助空间三维重建,完成物体所有特征点间实际距离的测量;利用一台相机获取在不同方位拍摄的多幅图像序列,算出拍摄每两张相片的相机的相对位置,建立相对坐标系,在相对坐标系内,利用立体匹配算法和三维重建算法对图像的特征点进行立体匹配与物体三维重建;对重建后的某两个特征点加入相对约束(两个特征点之间实际的相对距离),获得物体在相对坐标系内与物体实际尺寸的比例,作为尺度信息,从而完成利用物体三维重建达到摄影测量的目的。再利用三维重建之后特征点的三维坐标信息,构建一种自动基准面检测的Delaunay三角网格算法,根据网格信息求取合适基准面,用于进行不规则体体积快速计算。相对于以往摄影测量中的绝对位置测量修改为相对位置测量,本发明将通过对空间地物在相对坐标系中进行三维重建,再利用尺度信息获取物体各点之间实际的距离,从而实现体积的快速计算技术。技术人员只需要照相机、皮尺等测量工具,通过图像处理及三维重建技术,便可以实现快速体积测量。
在本发明中,利用特征点的二维坐标信息,改进现有的Delaunay构网算法,在三角网生长算法的基础上构建二维的Delaunay三角网格。在构网过程中,为了弥补Delaunay三角网格算法在体积计算领域的局限性,本方法加入了检测三角网边界边的数据成员useCount。若该数据成员值为0,说明对应的线段还没有在构建三角形的过程中被使用过;若该数据成员值为1,则说明对应的线段仅仅在一个构建的三角形中被使用过,并不涉及其他三角形,这样的三角边即为整个三角网格的边界边;若该数据成员值为2,则说明这条线段在构建三角形的过程中被使用了两次,即不是三角网格的边界边。构造了上述数据成员,则通过该数据成员的值即可将构成三角网边界的线段全部检测出来,再将其位置记录下来之后,便可获得这些折线段所构成的一个面。但由于其所构成的这个面并不是一个真正意义上的面,仅仅是一系列相连的折线段首尾接合所围成的一个封闭图形,因此,本方法采用最小二乘法,将其拟合成一个光滑的平面,作为不规则体体积计算的一个基准面,进一步用于Delaunay三角网所包围部分的体积计算。利用三维重建后物体在相对坐标系内的三维坐标信息、尺度信息,以及通过三角网格拟合出的基准面位置,将外围物体投影在基准面上,利用积分求体积的方法对三角网格所包围的不规则体进行体积计算。
附图说明
图1是本发明的系统流程示意图。
具体实施方式
下面结合附图对本发明做更进一步的解释。
如图1所示,一种基于空间不变特性的快速体积测量方法,包括如下步骤:
步骤(1),采用单台摄像机获取物体的多幅图像序列,并对所述摄像机进行标定:
首先,利用一台摄像机获取物体的图像序列,建立场景中物体在摄像机坐标系中与物体在图像平面上对应像点之间的关系,具体为:利用一台摄像机在不同的位置,从不同角度获取同一物体的多幅图像,所获取的每一幅图像均为物体的整体图像,每次小角度的改变相机位置进行拍摄,直至所摄的图像序列包含物体的所有信息,即所获取的多幅图像需要覆盖整个被测物体;然后进行摄像机标定与解算,获取摄像机的内部参数。
步骤(2),基于双目立体视觉技术,提取特征点后,对获取的多张图像序列进行立体匹配,以此找出任意两幅图像中同名像素点的对应坐标关系,即通过计算物体的空间点在两幅图像中的视差来获得其三维坐标值,具体步骤如下:
(21),对获取的图像进行预处理来提高图像质量,预处理包括图片平滑、锐化处理;
(22),特征点提取:对任意两幅图像,依据物体的物理特征来获取合适的特征点,即将角点、边缘点等易于分辨的特征信息作为后续特征点匹配的基元。本实施例采用FAST(Features From Accelerated Segment Test)特征点检测方法提取特征点:先通过一个角点响应函数CRF去判断各个待测点是否为FAST角点,其待测点即易于分辨的特征信息点,CRF函数如下:
N=ΣfCRF(I(p),I(x)) (1)
式(1)中,N为响应函数CRF值的和,即FAST角点的个数;I(x)为待测点周围的任意一点的灰度值;I(p)表示当前待测点的灰度值;p表示当前待测点;式(2)中,εd为一个中心待测点与圆周任一点像素灰度差值的阈值。
通过设定阈值εd,提取出超过阈值Nf数目的FAST角点;然后利用Harris角点的评价函数找到前Nf个较好的FAST角点;最后使用金字塔算法得到多尺度的待测图像,从而确定出最终要得到的FAST角点。
(23),采用BRIEF(Binary Robust Independent Element Feature)描述子对提取的FAST角点进行描述。
(24),利用ORB(Oriented FAST and rotated BRIEF)特征点匹配算法对选取的FAST角点进行立体匹配与解算,得到任意两幅图片特征点匹配的基础矩阵F;即多幅图像中,任意每两幅图像组合均得到其对应的一个基础矩阵F。
步骤(3),任意选取两幅图像,以拍摄其中一幅图像的摄像机为坐标原点建立相对坐标系,在该相对坐标系中对步骤(22)最终确定的所有FAST角点进行三维重建,得到整个物体的三维点云结构,具体为:
(31),任意选取两幅图像,以拍摄其中一幅图像的摄像机为坐标原点,以它的摄像机坐标系作为相对坐标系,计算得到拍摄这两幅图像的摄像机的位置关系,即摄像机的外部参数;
(32),结合摄像机标定时计算得到的摄像机内部参数,根据步骤(2)所得的这两幅图像的对应基础矩阵F,计算得到其本质矩阵E;
(33),对本质矩阵E进行奇异值分解,得到摄像机外部参数(R|t)的候选值,建立拍摄这两幅图像的两个不同位置的相机的投影矩阵P1和P2:
P1=K[I O] (3)
P2=K[R t] (4)
其中I是3x3的单位矩阵;O是3x1的全零矩阵;K为摄像机内部参数矩阵;R为3x3的旋转矩阵;t为三维平移向量;R、t均为相机外部参数;
(34),根据恢复出的投影矩阵P1和P2,在相对坐标系中利用SFM(Structure FromMotion)算法进行迭代计算,重建出这两幅图像上特征点对应的三维点云结构,即得到特征点在相对坐标系内的坐标信息。
由于本方案只使用了一台相机,故相机的内部参数是不变的;因此,在上述步骤(31)建立的相对坐标系下,按照步骤(31)至步骤(34)的方法,对其任意两幅图像做三维重建,得到整个物体在该相对坐标系下的三维点云结构。
步骤(4),在已知了所有选取的特征点在相对坐标系内的坐标信息后,计算出各个特征点之间的相对距离信息;然后在特征点中选取两个特征点,通过用皮尺丈量真实物体上着两个特征点的实际距离,实际距离与这两个特征点在相对坐标系内的相对距离的比值,即作为尺度信息。由此尺度信息可从三维点云结构中的任意两个特征点之间的相对距离计算出其在真实世界中的实际距离。
步骤(5),依据重建的三维点云结构建立物体的Delaunay三角网模型,基于该Delaunay三角网模型和所有特征点间的实际距离计算物体的体积,具体步骤如下:
(51),将重建后的三维点云结构投影到二维平面上,得到重建出来的特征点的二维坐标信息,然后依据三角网生长算法构建Delaunay三角网模型。构建Delaunay三角网格,即构建三种相关的数据结构:点、线、面,拓扑数据结构如下:
(1)点(T_point)类封装了离散采样点的三维坐标,构造出点的数据结构;
(2)线(T_Line)类封装了三角形的边的信息,其中每条边的两个端点的数据可以存储采样点在容器中的序号。
(3)面(T_tri)类封装了生成的每个三角形的信息,三角形的边信息以序号类型放入类中,可以随时调用,数据结构由边扩展到了面,用于最终三角网格的生成。
(4)三角网(T_TIN)类包含了构建Delaunay三角网格的核心函数,有了前面点、线、面数据结构的构建基础,可以方便构建出Delaunay三角网格。
构建上述三种具有拓扑关系的数据结构,能够便于对重建点云的整合和管理,利用三角网生长算法的思想进行三角网格的构建。
在构建Delaunay三角网模型时加入检测三角网边界边的数据成员useCount,数据成员useCount用于记录构建的每个三角形的每条边作为基础边来构建新三角形的次数;完成所有三角形建立步骤后,选择数据成员useCount为1的所有边作为Delaunay三角网模型的边界边。具体步骤如下:
(a),首先将重建出来的特征点的二维坐标信息放入初始数组容器中,并进行编号保存;
(b),以编号为0和1的点作为起始点,按照三角网判别法则构建一个寻找最大最小角的函数来找到满足Delaunay构网条件的第三个点,作为2号点;
(c),以步骤(b)中确立的3点作为初始三角形,并判断该三角形的三条边是否为可拓展边:若为可拓展边,则以可拓展边的两顶点作为基准再次运用寻找最大最小角的函数找寻第三个点,并与第三个点构成新的三角形,从而完成新三角网的拓展;若不是可拓展边,则停止拓展并给予该边标记,该不可拓展的边即是边界边;
(d),每遍历一条边,该边的数据成员useCount的值就相应的增加1,直到Delaunay三角网构网结束后,统计useCount的值,若该数据成员值为0,说明这条线段还没有在构建三角形中被使用过;若该数据成员为1,则说明该线段仅仅在一个三角形中被使用,即为整个三角网的边界边;若该数据成员为2,则说明这条线段在构建三角形的过程中被使用了两次,即非网格的边界边。
(52),完成所有三角形构建步骤后,根据数据成员useCount得到Delaunay三角网模型的所有边界边,再将由所有边界边构成的封闭图形拟合成一个光滑平面,该光滑平面作为计算物体体积的基准面,具体步骤如下:
在获取了所有的边界边信息后,找到所有边界边的两个端点的编号,并记录其端点的三维坐标信息;设待拟合的最优基准面的空间平面方程的表达式为:
ax+by+cz-d=0 (5)
其中,a、b、c、d均为方程系数,x、y、z为基准面坐标值;
要获得最佳拟合平面,应在条件a2+b2+c2=1的约束下,获取式(6)中e的最小值,即误差最小值:
其中,xi、yi、zi分别为空间平面上任意一点的坐标值;n为空间平面上点的个数;式(6)即等式约束的优化问题,利用拉格朗日乘数法求函数极值,令构造优化模型f(a,b,c,d,λ):
其中,λ为拉格朗日乘数;
令f分别对a,b,c,d求偏导,并令偏导数为零,整理如下:
由 得:
令 带入得:
令 由 得:
同理,由 和 得:
将式(10),式(11)和式(12)整理成矩阵形式,得:
由上述矩阵形式看出,可以将等式约束的优化问题转化为求解矩阵特征值和特征向量的问题,令 则上式化简为:
所以求出矩阵A的最小特征值λmin对应的特征向量α,就解算出了系数a,b,c的值,带入式(8)即可求出d的值。
由此,利用带约束的最小二乘算法将边界边拟合成了一个最优基准面,将用于后面的不规则体体积计算。
(53),利用拟合得到的最优基准面以及步骤(4)得到的尺度信息,利用积分求体积的方法求取物体的体积,即采用基于三棱柱体积计算的离散积分算法计算物体的体积:计算由Delaunay三角网格剖分后的每一个三角形与最优基准面投影构成的小三棱柱的体积,最后将所有三棱柱的体积进行加合,就得到了整个三角网所包围的物体的体积,具体的步骤如下:
设Delaunay三角网模型中任意三角形的三个顶点坐标为:A(x1,y1,z1),B(x2,y2,z2),C(x3,y3,z3);三个顶点投射到最优基准面的空间平面上的投影点为A1、B1、C1;已知最优基准面的空间平面的方程为:ax+by+cz=d;
(1)计算任意三角形的三个顶点投射到最优基准面上的投影点坐标A1(x1',y1',z1')、B1(x'2,y'2,z'2)、C1(x'3,y'3,z'3):
本实施例中使用向量法进行投影计算:
其中,G点为重建出三维坐标的每个特征点;向量为单位化的空间平面法向量;Q点为空间平面上的任意一点;G′点为投影结果。
(2)计算Delaunay三角网模型中任意三角形与最优基准面形成的单个三棱柱整体体积:
1)计算底面△A1B1C1的面积:
设
底面△A1B1C1的面积
2)计算单个三棱柱的高:
本方法中三棱柱的高采用三个点的高程平均值来估计,能够更加接近于真实体积,计算三棱柱的高步骤如下:
3)计算单个三棱柱的整体体积V为:
以上过程即为每个小三角网格与基准面构成的小三棱柱的体积,将所有三角形生成的三棱柱的体积相加合,即Vi即式(10)得到的单个三棱柱的整体体积;最后,将计算出体积V乘以步骤(4)得到的尺度信息的立方,即可以得到物体的总体积。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。
Claims (5)
1.一种基于空间不变特性的快速体积测量方法,其特征在于,包括如下步骤:
步骤(1),采用单台摄像机获取物体的多幅图像序列,并对所述摄像机进行标定;
步骤(2),基于双目立体视觉技术,提取所述多幅图像上的特征点后,对获取的多张图像序列进行立体匹配,找出任意两幅图像中同名像素点的对应坐标关系;
步骤(3),任意选取所述图像序列中的两幅图像,以拍摄其中一幅图像的摄像机为坐标原点建立相对坐标系,在所述相对坐标系中对所述步骤(2)提取的多幅图像序列上的所有特征点进行三维重建,得到整个物体的三维点云结构;
步骤(4),在所述整个物体的三维点云结构中,选择并测量两个特征点所对应物体上的实际距离,计算得到所述实际距离与这两个特征点在所述相对坐标系中的距离的比值,利用所述比值计算得到三维点云结构中所有特征点间的实际距离;
步骤(5),依据重建的三维点云结构建立物体的Delaunay三角网模型,根据所述Delaunay三角网模型获取计算物体体积的基准面,根据所述基准面以及所有特征点间的实际距离计算物体的体积。
2.根据权利要求1所述的一种基于空间不变特性的快速体积测量方法,其特征在于:依据三角网生长算法建立所述Delaunay三角网模型,并在构建所述Delaunay三角网模型时加入检测三角网边界边的数据成员useCount,所述数据成员useCount用于记录构建的每个三角形的每条边作为基础边来构建新三角形的次数;完成所有三角形建立步骤后,选择数据成员useCount为1的所有边作为Delaunay三角网模型的边界边,再将由所有边界边构成的封闭图形拟合成一个光滑平面,所述光滑平面作为计算物体体积的基准面。
3.根据权利要求1所述的一种基于空间不变特性的快速体积测量方法,其特征在于:所述步骤(2)具体为:
(21),对获取的图像进行预处理,所述预处理包括图片平滑、锐化处理;
(22),特征点提取:对任意两幅图像,先通过一个角点响应函数CRF去判断各个待测点是否为FAST角点,所述CRF函数如下:
N=ΣfCRF(l(p),l(x)) (1)
式(1)中,N为响应函数CRF值的和,即FAST角点的个数;l(x)为待测点周围的任意一点的灰度值;l(p)表示当前待测点的灰度值;p表示当前待测点;式(2)中,εd为一个中心待测点与圆周任一点像素灰度差值的阈值;
设定阈值εd,提取出超过阈值Nf数目的FAST角点;然后利用Harris角点的评价函数找到前Nf个较好的FAST角点;最后使用金字塔算法得到多尺度的待测图像,确定出最终的FAST角点;
(23),采用BRIEF描述子对提取的FAST角点进行描述;
(24),利用ORB特征点匹配算法对选取的FAST角点进行立体匹配与解算,得到任意两幅图片特征点匹配的基础矩阵F。
4.根据权利要求3所述的一种基于空间不变特性的快速体积测量方法,其特征在于:所述步骤(3)具体为:
(31),任意选取两幅图像,以拍摄其中一幅图像的摄像机为坐标原点,以它的摄像机坐标系作为相对坐标系,计算得到拍摄这两幅图像的摄像机的位置关系,即摄像机的外部参数;
(32),结合摄像机标定时得到的摄像机内部参数,根据步骤(2)所得的这两幅图像的对应基础矩阵F,计算得到其本质矩阵E;
(33),对本质矩阵E进行奇异值分解,得到摄像机外部参数(R|t)的候选值,建立拍摄这两幅图像的两个不同位置的相机的投影矩阵P1和P2:
P1=K[I O] (3)
P2=K[R t] (4)
其中I是3x3的单位矩阵;O是3x1的全零矩阵;K为摄像机内部参数矩阵;R为3x3的旋转矩阵;t为三维平移向量;R、t均为相机外部参数;
(34),根据恢复出的投影矩阵P1和P2,在相对坐标系中利用SFM算法进行迭代计算,重建出这两幅图像上特征点对应的三维点云结构,即得到特征点在相对坐标系内的坐标信息;
(35),在所述步骤(31)建立的相对坐标系下,按照步骤(31)至步骤(34),对其任意两幅图像做三维重建,得到整个物体在该相对坐标系下的三维点云结构。
5.根据权利要求4所述的一种基于空间不变特性的快速体积测量方法,其特征在于:所述步骤(5)中,根据所述基准面以及所有特征点间的实际距离计算物体的体积具体步骤为:
(1),设Delaunay三角网模型中任意三角形的三个顶点坐标为:A(x1,y1,z1),B(x2,y2,z2),C(x3,y3,z3),使用向量法进行投影计算任意三角形的三个顶点投射到所述基准面上的投影点坐标A1(x′1,y′1,z′1)、B1(x′2,y′2,z′2)、C1(x′3,y′3,z′3);
(2),计算Delaunay三角网模型中任意三角形与所述基准面形成的单个三棱柱整体体积:
21),计算底面△A1B1C1的面积:
设
底面△A1B1C1的面积
22),采用任意三角形的三个顶点的高程平均值来计算单个三棱柱的高h:
23),计算单个三棱柱的整体体积V为:
(3),将Delaunay三角网模型中所有的n个三角形对应生成的三棱柱的体积相加合,即Vt为式(10)得到的单个三棱柱的整体体积;
(4),将计算出体积V乘以步骤(4)得到的尺度信息的立方,得到物体的总体积。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510231222.9A CN104778720B (zh) | 2015-05-07 | 2015-05-07 | 一种基于空间不变特性的快速体积测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510231222.9A CN104778720B (zh) | 2015-05-07 | 2015-05-07 | 一种基于空间不变特性的快速体积测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104778720A true CN104778720A (zh) | 2015-07-15 |
CN104778720B CN104778720B (zh) | 2018-01-16 |
Family
ID=53620166
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510231222.9A Expired - Fee Related CN104778720B (zh) | 2015-05-07 | 2015-05-07 | 一种基于空间不变特性的快速体积测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104778720B (zh) |
Cited By (30)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106408650A (zh) * | 2016-08-26 | 2017-02-15 | 中国人民解放军国防科学技术大学 | 在轨掠飞成像对空间目标三维重建与测量方法 |
CN106643494A (zh) * | 2016-12-22 | 2017-05-10 | 上海华测导航技术股份有限公司 | 一种矿山堆料体积测量方法及系统 |
CN106652020A (zh) * | 2016-12-05 | 2017-05-10 | 成都通甲优博科技有限责任公司 | 一种基于模型的塔杆三维重建方法 |
CN106954012A (zh) * | 2017-03-29 | 2017-07-14 | 武汉嫦娥医学抗衰机器人股份有限公司 | 一种高清多相机全景立体成像系统及方法 |
CN107481282A (zh) * | 2017-08-18 | 2017-12-15 | 成都通甲优博科技有限责任公司 | 体积测算方法、装置及用户终端 |
CN108182722A (zh) * | 2017-07-27 | 2018-06-19 | 桂林航天工业学院 | 一种三维物体边缘优化的真正射影像生成方法 |
CN108446597A (zh) * | 2018-02-14 | 2018-08-24 | 天目爱视(北京)科技有限公司 | 一种基于可见光相机的生物特征3d数据采集方法及装置 |
CN108470151A (zh) * | 2018-02-14 | 2018-08-31 | 天目爱视(北京)科技有限公司 | 一种生物特征模型合成方法及装置 |
CN108492330A (zh) * | 2018-02-14 | 2018-09-04 | 天目爱视(北京)科技有限公司 | 一种多目视觉深度计算方法及装置 |
CN108520230A (zh) * | 2018-04-04 | 2018-09-11 | 北京天目智联科技有限公司 | 一种3d四维手部图像数据识别方法及设备 |
CN109003307A (zh) * | 2018-06-11 | 2018-12-14 | 西北工业大学 | 基于水下双目视觉测量的捕鱼网目尺寸设计方法 |
CN109035321A (zh) * | 2017-06-09 | 2018-12-18 | 河北卓达建材研究院有限公司 | 一种建筑物的体积估算方法 |
CN109089100A (zh) * | 2018-08-13 | 2018-12-25 | 西安理工大学 | 一种双目立体视频的合成方法 |
CN109163775A (zh) * | 2018-08-30 | 2019-01-08 | 北京广天夏科技有限公司 | 一种基于带式运输机的质量测量方法及装置 |
CN109272542A (zh) * | 2018-08-24 | 2019-01-25 | 华南理工大学 | 一种三维物体体积的确定方法 |
CN109509535A (zh) * | 2018-10-08 | 2019-03-22 | 北京健康有益科技有限公司 | 食物体积的获取方法、食物热量的获取方法、电子设备 |
CN109870106A (zh) * | 2019-02-14 | 2019-06-11 | 南京邮电大学 | 一种基于无人机图片的建筑物体积测量方法 |
CN109903327A (zh) * | 2019-03-04 | 2019-06-18 | 西安电子科技大学 | 一种稀疏点云的目标物尺寸测量方法 |
CN109990703A (zh) * | 2019-03-18 | 2019-07-09 | 桂林电子科技大学 | 一种预制构件的尺寸检测方法及系统 |
CN110095062A (zh) * | 2019-04-17 | 2019-08-06 | 北京华捷艾米科技有限公司 | 一种物体体积参数测量方法、装置及设备 |
CN110260801A (zh) * | 2019-05-13 | 2019-09-20 | 平安科技(深圳)有限公司 | 用于测量物料体积的方法和装置 |
CN110276732A (zh) * | 2019-06-18 | 2019-09-24 | 昆明理工大学 | 一种顾及地形特征线要素的山区点云空洞修复方法 |
CN110532719A (zh) * | 2019-09-04 | 2019-12-03 | 上海宝冶集团有限公司 | 一种施工场地土方量计算方法及装置 |
CN111156922A (zh) * | 2019-12-23 | 2020-05-15 | 苏州迭慧智能科技有限公司 | 一种利用外形轮廓测量的方法 |
WO2020124976A1 (zh) * | 2018-12-21 | 2020-06-25 | 北京市商汤科技开发有限公司 | 图像处理方法及装置、电子设备及存储介质 |
CN111724436A (zh) * | 2020-06-23 | 2020-09-29 | 哈尔滨工业大学 | 一种基于点云数据处理的隧道喷射方量估算方法 |
CN111811655A (zh) * | 2020-06-03 | 2020-10-23 | 浙江大学 | 光的偏振状态评测方法、应用、装置、电子设备及介质 |
CN112365534A (zh) * | 2020-11-12 | 2021-02-12 | 上海明华电力科技有限公司 | 一种基于单目相机三维重建的大型煤堆体积测量方法 |
CN112950702A (zh) * | 2021-02-01 | 2021-06-11 | 华电淄博热电有限公司 | 一种基于三维点云的煤堆体积计算方法 |
CN114220251A (zh) * | 2021-12-23 | 2022-03-22 | 西安芯海微电子科技有限公司 | 一种水位监测预警电路、方法和电子设备 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101281662A (zh) * | 2008-05-16 | 2008-10-08 | 中山大学 | 一种三棱柱体三维地层剖分方法及开挖土方估算方法 |
US20120229595A1 (en) * | 2011-03-11 | 2012-09-13 | Miller Michael L | Synthesized spatial panoramic multi-view imaging |
CN102692190A (zh) * | 2011-03-24 | 2012-09-26 | 北京林业大学 | 一种基于数字高程模型的树木体积计测方法 |
US20120275667A1 (en) * | 2011-04-29 | 2012-11-01 | Aptina Imaging Corporation | Calibration for stereoscopic capture system |
-
2015
- 2015-05-07 CN CN201510231222.9A patent/CN104778720B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101281662A (zh) * | 2008-05-16 | 2008-10-08 | 中山大学 | 一种三棱柱体三维地层剖分方法及开挖土方估算方法 |
US20120229595A1 (en) * | 2011-03-11 | 2012-09-13 | Miller Michael L | Synthesized spatial panoramic multi-view imaging |
CN102692190A (zh) * | 2011-03-24 | 2012-09-26 | 北京林业大学 | 一种基于数字高程模型的树木体积计测方法 |
US20120275667A1 (en) * | 2011-04-29 | 2012-11-01 | Aptina Imaging Corporation | Calibration for stereoscopic capture system |
Cited By (41)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106408650B (zh) * | 2016-08-26 | 2018-12-11 | 中国人民解放军国防科学技术大学 | 在轨掠飞成像对空间目标三维重建与测量方法 |
CN106408650A (zh) * | 2016-08-26 | 2017-02-15 | 中国人民解放军国防科学技术大学 | 在轨掠飞成像对空间目标三维重建与测量方法 |
CN106652020A (zh) * | 2016-12-05 | 2017-05-10 | 成都通甲优博科技有限责任公司 | 一种基于模型的塔杆三维重建方法 |
CN106652020B (zh) * | 2016-12-05 | 2019-12-24 | 成都通甲优博科技有限责任公司 | 一种基于模型的塔杆三维重建方法 |
CN106643494A (zh) * | 2016-12-22 | 2017-05-10 | 上海华测导航技术股份有限公司 | 一种矿山堆料体积测量方法及系统 |
CN106954012A (zh) * | 2017-03-29 | 2017-07-14 | 武汉嫦娥医学抗衰机器人股份有限公司 | 一种高清多相机全景立体成像系统及方法 |
CN109035321A (zh) * | 2017-06-09 | 2018-12-18 | 河北卓达建材研究院有限公司 | 一种建筑物的体积估算方法 |
CN108182722A (zh) * | 2017-07-27 | 2018-06-19 | 桂林航天工业学院 | 一种三维物体边缘优化的真正射影像生成方法 |
CN108182722B (zh) * | 2017-07-27 | 2021-08-06 | 桂林航天工业学院 | 一种三维物体边缘优化的真正射影像生成方法 |
CN107481282A (zh) * | 2017-08-18 | 2017-12-15 | 成都通甲优博科技有限责任公司 | 体积测算方法、装置及用户终端 |
CN107481282B (zh) * | 2017-08-18 | 2020-04-28 | 成都通甲优博科技有限责任公司 | 体积测算方法、装置及用户终端 |
CN108492330A (zh) * | 2018-02-14 | 2018-09-04 | 天目爱视(北京)科技有限公司 | 一种多目视觉深度计算方法及装置 |
CN108470151A (zh) * | 2018-02-14 | 2018-08-31 | 天目爱视(北京)科技有限公司 | 一种生物特征模型合成方法及装置 |
CN108446597A (zh) * | 2018-02-14 | 2018-08-24 | 天目爱视(北京)科技有限公司 | 一种基于可见光相机的生物特征3d数据采集方法及装置 |
CN108446597B (zh) * | 2018-02-14 | 2019-06-25 | 天目爱视(北京)科技有限公司 | 一种基于可见光相机的生物特征3d数据采集方法及装置 |
CN108520230A (zh) * | 2018-04-04 | 2018-09-11 | 北京天目智联科技有限公司 | 一种3d四维手部图像数据识别方法及设备 |
CN109003307A (zh) * | 2018-06-11 | 2018-12-14 | 西北工业大学 | 基于水下双目视觉测量的捕鱼网目尺寸设计方法 |
CN109003307B (zh) * | 2018-06-11 | 2021-10-22 | 西北工业大学 | 基于水下双目视觉测量的捕鱼网目尺寸设计方法 |
CN109089100B (zh) * | 2018-08-13 | 2020-10-23 | 西安理工大学 | 一种双目立体视频的合成方法 |
CN109089100A (zh) * | 2018-08-13 | 2018-12-25 | 西安理工大学 | 一种双目立体视频的合成方法 |
CN109272542A (zh) * | 2018-08-24 | 2019-01-25 | 华南理工大学 | 一种三维物体体积的确定方法 |
CN109163775A (zh) * | 2018-08-30 | 2019-01-08 | 北京广天夏科技有限公司 | 一种基于带式运输机的质量测量方法及装置 |
CN109509535A (zh) * | 2018-10-08 | 2019-03-22 | 北京健康有益科技有限公司 | 食物体积的获取方法、食物热量的获取方法、电子设备 |
JP2021520577A (ja) * | 2018-12-21 | 2021-08-19 | 北京市商▲湯▼科技▲開▼▲發▼有限公司Beijing Sensetime Technology Development Co.,Ltd. | 画像処理方法及び装置、電子機器並びに記憶媒体 |
JP7026825B2 (ja) | 2018-12-21 | 2022-02-28 | 北京市商▲湯▼科技▲開▼▲發▼有限公司 | 画像処理方法及び装置、電子機器並びに記憶媒体 |
WO2020124976A1 (zh) * | 2018-12-21 | 2020-06-25 | 北京市商汤科技开发有限公司 | 图像处理方法及装置、电子设备及存储介质 |
CN109870106A (zh) * | 2019-02-14 | 2019-06-11 | 南京邮电大学 | 一种基于无人机图片的建筑物体积测量方法 |
CN109903327A (zh) * | 2019-03-04 | 2019-06-18 | 西安电子科技大学 | 一种稀疏点云的目标物尺寸测量方法 |
CN109990703A (zh) * | 2019-03-18 | 2019-07-09 | 桂林电子科技大学 | 一种预制构件的尺寸检测方法及系统 |
CN110095062A (zh) * | 2019-04-17 | 2019-08-06 | 北京华捷艾米科技有限公司 | 一种物体体积参数测量方法、装置及设备 |
CN110260801A (zh) * | 2019-05-13 | 2019-09-20 | 平安科技(深圳)有限公司 | 用于测量物料体积的方法和装置 |
CN110276732B (zh) * | 2019-06-18 | 2022-06-07 | 昆明理工大学 | 一种顾及地形特征线要素的山区点云空洞修复方法 |
CN110276732A (zh) * | 2019-06-18 | 2019-09-24 | 昆明理工大学 | 一种顾及地形特征线要素的山区点云空洞修复方法 |
CN110532719A (zh) * | 2019-09-04 | 2019-12-03 | 上海宝冶集团有限公司 | 一种施工场地土方量计算方法及装置 |
CN111156922A (zh) * | 2019-12-23 | 2020-05-15 | 苏州迭慧智能科技有限公司 | 一种利用外形轮廓测量的方法 |
CN111811655A (zh) * | 2020-06-03 | 2020-10-23 | 浙江大学 | 光的偏振状态评测方法、应用、装置、电子设备及介质 |
CN111724436A (zh) * | 2020-06-23 | 2020-09-29 | 哈尔滨工业大学 | 一种基于点云数据处理的隧道喷射方量估算方法 |
CN112365534A (zh) * | 2020-11-12 | 2021-02-12 | 上海明华电力科技有限公司 | 一种基于单目相机三维重建的大型煤堆体积测量方法 |
CN112365534B (zh) * | 2020-11-12 | 2023-02-03 | 上海明华电力科技有限公司 | 一种基于单目相机三维重建的大型煤堆体积测量方法 |
CN112950702A (zh) * | 2021-02-01 | 2021-06-11 | 华电淄博热电有限公司 | 一种基于三维点云的煤堆体积计算方法 |
CN114220251A (zh) * | 2021-12-23 | 2022-03-22 | 西安芯海微电子科技有限公司 | 一种水位监测预警电路、方法和电子设备 |
Also Published As
Publication number | Publication date |
---|---|
CN104778720B (zh) | 2018-01-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104778720B (zh) | 一种基于空间不变特性的快速体积测量方法 | |
CN103363962B (zh) | 一种基于多光谱影像的湖泊水储量遥感估算方法 | |
CN106896213B (zh) | 一种基于点云数据的岩体结构面智能识别与信息提取方法 | |
CN109472802B (zh) | 一种基于边缘特征自约束的表面网格模型构建方法 | |
CN103604417B (zh) | 物方信息约束的多视影像双向匹配策略 | |
CN104952107A (zh) | 基于车载LiDAR点云数据的桥梁三维重建方法 | |
CN106886794B (zh) | 顾及高阶结构特征的异源遥感影像同名点匹配方法 | |
CN104268935A (zh) | 一种基于特征的机载激光点云与影像数据融合系统及方法 | |
CN109342337A (zh) | 一种重度盐渍土盐分获取方法、系统及装置 | |
CN102788979B (zh) | 一种基于后向投影InSAR成像配准的GPU实现方法 | |
CN111667569B (zh) | 一种基于Rhino及Grasshopper的三维实景土方可视化精准测算方法 | |
CN112698328B (zh) | 一种用于大坝及滑坡变形gb-sar监测的相位解缠方法及系统 | |
CN103454636A (zh) | 基于多像素协方差矩阵的差分干涉相位估计方法 | |
CN112801983A (zh) | 基于摄影测量的边坡全域沉降检测方法及系统 | |
CN115877421A (zh) | 一种输电通道地质敏感区的形变检测方法及装置 | |
Yang et al. | Linking persistent scatterers to the built environment using ray tracing on urban models | |
CN108447126B (zh) | 基于参考平面的移动测量系统激光点云精度评估方法 | |
CN110009740A (zh) | 基于运动恢复结构的地质露头快速三维重建方法 | |
Yaseen et al. | Local interpolation of coseismic displacements measured by InSAR | |
CN106846477B (zh) | 一种编录野外地质影像的地质标记解译建模方法 | |
Zhang et al. | Extraction of three-dimensional architectural data from QuickBird images | |
CN116994012A (zh) | 一种基于生态修复的图斑匹配系统及匹配方法 | |
Zhang et al. | 3D Visualization of Landslide Based on Close-Range Photogrammetry. | |
RU2570334C1 (ru) | Способ картографирования ледниковой геоморфологии | |
CN114239379A (zh) | 一种基于形变检测的输电线路地质灾害分析方法及系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
EXSB | Decision made by sipo to initiate substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180116 |
|
CF01 | Termination of patent right due to non-payment of annual fee |