本发明专利技术公开了一种利用单摄像头传回的图片，自动测量前景物体高度的新方法。此方法采用数字图像处理技术，在计算前景高度的过程中，首先通过标定，取得图像平面和世界坐标系之间像素点的映射关系，即计算出：单应矩阵ＨｇＭａｔｉｘ，竖直方向的消影点ｍｖ，横向的消影点ｍｘ，纵向的消影点ｍｙ，进而计算出世界坐标系和图像平面坐标系之间的投影矩阵ＰｒｏＭａｔｉｘ；然后，取得前景区域的高点ｈｉｇｈｎｏｄ和低点ｌｏｗｎｏｄ，再计算绝对高点ｔｒｕｅｈｉｇｈ及高点ｈｉｇｈｎｏｄ在地平面上的投影点即绝对低点ｔｒｕｅｌｏｗ；最后计算出绝对高点ｔｒｕｅｈｉｇｈ的Ｚ值和ｌｏｗｎｏｄ与ｔｒｕｅｌｏｗ所对应的时间坐标点之间的距离值，通过勾股定理，就可取得前景物体高度的较精确的估计值。此方法解决了倾斜前景物体的高度测量问题，提高前景物体高度估计测量的精确度；又由于本方法采用的是单摄像头的测量方法，降低了实际应用时的成本。
【技术实现步骤摘要】
本专利技术方法涉及一种通过单摄像头传回的图片，自动测量前景物体高度的新方 法，广泛适用于单摄像头估测前景高度的问题，如草的高度估测，牛羊高度的估测等 等；同时，本方法考虑到现实中前景不一定都垂直的情况，对于倾斜前景也给出了比较 好的估测精度。
技术介绍
随着计算软硬件的发展，利用计算机视觉方法对场景进行测量已发展成为计算 机视觉中的一个重要课题，计算机视觉测量已弓I起越来越多的学者兴趣和关注。目前计 算机视觉测量方法主要有下面几种一是通过同一场景的两幅或多幅图像，重建场景 的三维结构，得到前景物体的三维数据从而达到测量目的。这种方法的估测精度对图像 间匹配点的定位精度要求很高，这样参数误差都会对测量的结果的影响就比较大，但是 图像匹配点的精确定位又是计算机视觉领域的经典难题；二是采用单幅图像测量方法， 通过图像平面和世界坐标系间的投影矩阵计算出物体的三维结构，这种方法避开了图像 间匹配点的定位问题，从而有比较高的测量精度。
技术实现思路
本专利技术的目的是为了通过单幅图像测量前景物体的高度，实现了对倾斜物体高 度的较高精度测量。扩展了单幅图像测量的应用范围。本专利技术公开了一种利用单摄像头传回的图片，自动测量前景物体高度的新方 法。此方法采用数字图像处理技术，在计算前景高度的过程中，首先通过标定，取得图 像平面和世界坐标系之间像素点的映射关系，即计算出单应矩阵HgMatix，竖直方向 的消影点mv，横向的消影点mx，纵向的消影点my，进而计算出世界坐标系和图像平面 坐标系之间的投影矩阵ProMatix;然后，取得前景区域的高点highnod和低点lownod，再 计算绝对高点truehigh和highnod点对应的投影点truelow ；最后计算出绝对高点truehigh 的Z值和lownod与truelow所对应的时间坐标点之间的距离值，通过勾股定理，就可取得 前景物体高度的较精确的估计值。此方法解决了倾斜前景物体的高度测量问题，提高前 景物体高度估计测量的精确度；又由于本方法采用的是单摄像头的测量方法，降低了实 际应用时的成本。本专利技术方法的条件和步骤如下1.所需设备待校正的成像系统，的普通PC机，标定模版；其中成像系统可以 是相机或摄像机等(本文采用的是摄像头)；标定模版为自制的黑白起棋盘板。2.实现步骤本专利技术方法(010)部分，得到世界坐标系的地平面与图像平面之间的单应矩阵 Hg，以及垂直方向的消影点mvPoint，横向消影点mxPoint，纵向消影点myPoint，具体步骤如下步骤C011:开始，利用摄像头拍摄标定模板包含的棋盘纸，这里要拍摄的包含 两个方向的标定模板图像，水平放置即放置在地平面上的标定模板图像plane—Img，和垂 直地平面放置的垂直模板图像vertical—Img。通过普通PC机读取包含这两个方向的标定 模板的图片——之后摄像头固定不动。步骤C012:对于水平放置的标定模板图像plane—Img，找到模板图像plane— Img上的角点及其在图像平面上的坐标COmer，并显示出来；对于垂直放置的标定 模板图像vertical—Img，找到模板图像vertical—Img上的角点及其在图像平面上的坐标 corner—l，并显示出来；其中i = 1，2，…，m，这里m表示网格上横向角点个数，j = 1, 2，…，n,这里n表示网格上纵向角点个数；步骤C013 接着上步C012得出的COmer，假定找到的第一个角点坐标对应 的地平面齐次坐标为(0，0，0，1)，第二个角点坐标对应的地平面齐次坐标为(0，1， 0，1)，依次类推，找出每个图像平面上角点对应的地平面坐标。根据这些相对应的点， 计算出图像平面与地平面之间的单应矩阵HgMatix。步骤C014:接着C012得出的COrner，分别拟合横向直线和纵向直线，由拟 合出的一系列横向直线，应用最小二乘，算出它们的交点即横向消影点mxPoint;由拟合 出的一系列纵向直线，应用最小二乘，算出它们的交点即纵向消影点myPoint;步骤C015 接着C012得出的corner—l，拟合出纵向直线，由拟合出的一系 列纵向直线，应用最小二乘，算出它们的交点即垂直方向的消影点mvPoint;本专利技术方法(020)部分，计算摄像头的内部参数Intri和世界坐标系与图像平面坐 标系之间的投影矩阵ProMatix，具体步骤如下步骤C021:根据步骤C013得到的单应矩阵Hg和步骤C015得到的mvPoint，由 单应矩阵与摄像头内参数之间的关系、旋转矩阵的特性，列出三组线性方程，计算出内 部参数Intri。步骤C022 根据步骤CO 13得到的单应矩阵HgMatix和步骤CO 15得到的 mvPoint,以及步骤C021得出的Intri，由单应矩阵HgMatix与投影矩阵ProMatix、垂直 方向消影点mvPoint之间的关系，计算出世界坐标系与图像平面坐标系之间的投影矩阵 ProMatix o本专利技术方法(030)部分，计算前景的高度high具体步骤如下步骤C031 对于前景的提取，不同的应用环境有不同的方法。本文采用的是高 斯混合模型的前景提取法，提取前景区域front—Img。步骤C032 计算垂直消影点mvPoint和前景区域front—Img所有点之间的距离，取最大值hH和最小值1H对应的图像点。如果消影点的纵向值大于零，则取最大值对 应的点为最高点highnod，最小值对应的图像点为最低点lownod ；如消影点的纵向值小于 零，则相反。步骤C033:计算C032中最大值hH和最小值1H的差值length，在最低点lownod和垂直消影点mvPoint连线上，从最低点lownod向最高点所在方向取长度为length的线 段，交点为绝对高点truehigh ； 步骤C034 垂直消影点mvPoint与最高点highnod的连线L1，横向消影点 mxPoint与最低点lownod的连线L2，纵向消影点myPoint与最低点lownod的连线L3 ；分别计算L1与L2的交点P1和L1与L3的交点P2 ；再计算lownod与P1的距离和lownod 与P2的距离，取其距离值小者对应的交点(P1或P2)作为绝对低点truelow ；步骤C035 根据单应矩阵HgMatix计算最低点lownod和绝对低点truelow在 地平面上的坐标，并计算出在地平面上的距离lengthL ；由算出的最低点lownod对应在 地平面上的坐标和绝对高点truehigh，计算出绝对高点truehigh在世界坐标系上的Z值 truehighZ ；步骤C036:上步的lengthL和truehighZ，根据勾股定理，计算出图像上最高点highnod和最低点lownod对应的世界坐标系中的距离，即前景的高度high。有益效果本专利技术的目的是提供一种能够计算任意放置(包括倾斜放置)物体的高度。其 优点在于采用数字图像处理技术，方法灵活，易于实现，成本较低；同时使用该方法简 单通用，对于任意相机或摄像机，任意分辨率的图像都可使用。附图说明图1为水平标定模板图像plane—Img，和垂直模板图像vertical—Img ；图2为地平面上的横向消影点mxPoint和纵向消影点my本文档来自技高网...
【技术保护点】
单副图像测量前景高度的方法，设备由待校正的成像系统，普通ＰＣ机，标定模版组成；其中成像系统可以是相机或摄像机等（本文采用的是摄像头）；标定模版为自制的黑白起棋盘板。其步骤如下：１．实现步骤：本专利技术方法（０１０）部分，得到世界坐标系的地平面与图像平面之间的单应矩阵Ｈｇ，以及垂直方向的消影点ｍｖＰｏｉｎｔ，横向消影点ｍｘＰｏｉｎｔ，纵向消影点ｍｙＰｏｉｎｔ，具体步骤如下：步骤Ｃ０１１：开始，利用摄像头拍摄标定模板包含的棋盘纸，这里要拍摄的包含两个方向的标定模板图像，水平放置即放置在地平面上的标定模板图像ｐｌａｎｅ＿Ｉｍｇ，和垂直地平面放置的垂直模板图像ｖｅｒｔｉｃａｌ＿Ｉｍｇ。通过普通ＰＣ机读取包含这两个方向的标定模板的图片－－之后摄像头固定不动。步骤Ｃ０１２：对于水平放置的标定模板图像ｐｌａｎｅ＿Ｉｍｇ，找到模板图像ｐｌａｎｅ＿Ｉｍｇ上的角点及其在图像平面上的坐标ｃｏｒｎｅｒ［ｉ↑［＊］ｊ］，并显示出来；对于垂直放置的标定模板图像ｖｅｒｔｉｃａｌ＿Ｉｍｇ，找到模板图像ｖｅｒｔｉｃａｌ＿Ｉｍｇ上的角点及其在图像平面上的坐标ｃｏｒｎｅｒ＿１［ｉ↑［＊］ｊ］，并显示出来；其中ｉ＝１，２，…，ｍ，这里ｍ表示网格上横向角点个数，ｊ＝１，２，…，ｎ，这里ｎ表示网格上纵向角点个数；步骤Ｃ０１３：接着上步Ｃ０１２得出的ｃｏｒｎｅｒ［ｉ↑［＊］ｊ］，假定找到的第一个角点坐标对应的地平面齐次坐标为（０，０，０，１），第二个角点坐标对应的地平面齐次坐标为（０，１，０，１），依次类推，找出每个图像平面上角点对应的地平面坐标。根据这些相对应的点，计算出图像平面与地平面之间的单应矩阵ＨｇＭａｔｉｘ。步骤Ｃ０１４：接着Ｃ０１２得出的ｃｏｒｎｅｒ［ｉ↑［＊］ｊ］，分别拟合横向直线和纵向直线，由拟合出的一系列横向直线，应用最小二乘，算出它们的交点即横向消影点ｍｘＰｏｉｎｔ；由拟合出的一系列纵向直线，应用最小二乘，算出它们的交点即纵向消影点ｍｙＰｏｉｎｔ；步骤Ｃ０１５：接着Ｃ０１２得出的ｃｏｒｎｅｒ＿１［ｉ↑［＊］ｊ］，拟合出纵向直线，由拟合出的一系列纵向直线，应用最小二乘，算出它们的交点即垂直方向的消影点ｍｖＰｏｉｎｔ；本专利技术方法（０２０）部分，计算摄像头的内部参数Ｉｎｔｒｉ和世界坐标系与图像平面坐标系之间的投影矩阵ＰｒｏＭａｔｉｘ，具体步骤如下：步骤Ｃ０２１：根据步骤Ｃ０１３得到的单应矩阵Ｈｇ和步骤Ｃ０１５得到的ｍｖＰｏｉｎｔ，由单应矩阵与摄像头内参数之间的关系、旋转矩阵的特性，列出三组线性方程，计算出内部参数Ｉｎｔｒｉ。步骤Ｃ０２２：根据步骤Ｃ０１３得到的单应矩阵ＨｇＭａｔｉｘ和步骤Ｃ０１５得到的ｍｖＰｏｉｎｔ，以及步骤Ｃ０２１得出的Ｉｎｔｒｉ，由单应矩阵ＨｇＭａｔｉｘ与投影矩阵ＰｒｏＭａｔｉｘ、垂直方向消影点ｍｖＰｏｉｎｔ之间的关系，计算出世界坐标系与图像平面坐标系之间的投影矩阵ＰｒｏＭａｔｉｘ。本专利技术方法（０３０）部分，计算前景的高度ｈｉｇｈ具体步骤如下：步骤Ｃ０３１：对于前景的提取，不同的应用环境有不同的方法。本文采用的是高斯混合模型的前景提取法，提取前景区域ｆｒｏｎｔ＿Ｉｍｇ。步骤Ｃ０３２：计算垂直消影点ｍｖＰｏｉｎｔ和前景区域ｆｒｏｎｔ＿Ｉｍｇ所有点之间的距离，取最大值ｈＨ和最小值１Ｈ对应的图像点。如果消影点的纵向值大于零，则取最大值对应的点为最高点ｈｉｇｈｎｏｄ，最小值对应的图像点为最低点ｌｏｗｎｏｄ；如消影点的纵向值小于零，则相反。步骤Ｃ０３３：计算Ｃ０３２中最大值ｈＨ和最小值１Ｈ的差值ｌｅｎｇｔｈ，在最低点ｌｏｗｎｏｄ和垂直消影点ｍｖＰｏｉｎｔ连线上，从最低点ｌｏｗｎｏｄ向最高点所在方向取长度为ｌｅｎｇｔｈ的线段，交点为绝对高点ｔｒｕｅｈｉｇｈ；步骤Ｃ０３４：垂直消影点ｍｖＰｏｉｎｔ与最高点ｈｉｇｈｎｏｄ的连线Ｌ１，横向消影点ｍｘＰｏｉｎｔ与最低点ｌｏｗｎｏｄ的连线Ｌ２，纵向消影点ｍｙＰｏｉｎｔ与最低点ｌｏｗｎｏｄ的连线Ｌ３；分别计算Ｌ１与Ｌ２的交点Ｐ１和Ｌ１与Ｌ３的交点Ｐ２；再计算ｌｏｗｎｏｄ与Ｐ１的距离和ｌｏｗｎｏｄ与Ｐ２的距离，取其距离值小者对应的交点（Ｐ１或Ｐ２）作为绝对低点ｔｒｕｅｌｏｗ；步骤Ｃ０３５：根据单应矩阵ＨｇＭａｔｉｘ计算最低点ｌｏｗｎｏｄ和绝对低点ｔｒｕｅｌｏｗ在地平面上的坐标，并计算出在地平面上的距离ｌｅｎｇｔｈＬ；由算出的最低点ｌｏｗｎｏｄ对应在地平面上的坐标和绝对高点ｔｒｕｅｈｉｇｈ，计算出绝对高点ｔｒｕｅｈｉｇｈ在世界坐标系上的Ｚ值ｔｒｕｅｈｉｇｈＺ；步骤Ｃ０３６：上步的ｌｅｎｇｔｈＬ和ｔｒｕｅｈｉｇｈＺ，根据勾股定理，计算出图像上最高点ｈｉｇｈｎｏｄ和最低点ｌｏｗｎｏｄ对应的世界坐标系中的距离，即前景的高度ｈｉｇｈ。...
【技术特征摘要】
【专利技术属性】
技术研发人员：李文辉，齐刚，傅博，
申请(专利权)人：吉林大学，
类型：发明
国别省市：82
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。