一种非朗伯表面快速三维重构方法技术

本发明专利技术属于图像理解与计算机视觉技术领域，特别涉及一种非朗伯表面快速三维重构方法。提供的技术方案是：一种非朗伯表面快速三维重构方法，依次包括下述步骤：一．构建适于非朗伯表面且接近实际成像条件下的成像模型；二．由上述成像模型和摄像机获取的灰度图像，建立非朗伯表面下SFS问题的图像辐照度方程；三．利用粘性解理论计算上述图像辐照度方程的粘性解；四．根据图像辐照度方程解的结果，重构出非朗伯表面的三维形状信息。本发明专利技术的优点是：直接对非朗伯表面下建立的图像辐照度方程进行求解，同时利用弱解理论——粘性解理论计算图像辐照度方程的粘性解，使非朗伯表面的SFS问题良性化；具有快速、准确的优点。

【技术实现步骤摘要】

本专利技术属于图像理解与计算机视觉
，主要涉及到由单幅图像明暗变化实现三维重构的技术，特别涉及到。
技术介绍
为了满足工农业产品检测、文档图像复原、医学图像分析、生物特征识别以及地形、地貌测量等领域中三维建模的需求，有效重构物体表面的三维形状信息，近年来已经成为图像理解与计算机视觉
中一个非常重要的研究方向。然而，如何快速、准确重构非朗伯表面的三维形状，至今仍未得到较好解决。目前，三维形状信息的获取手段主要有两大类接触式测量和非接触式测量。尽管接触式测量方法已经比较成熟，而且还可以获得很高精度的数据，但其本身也存在较多不足①接触式测量使得测头容易磨损及划伤被测表面，测量过程需要人工干预，因而测量成本较高，自动化程度较低；②接触式测量采用逐点扫描的方法，且扫描速度受到机械运动的限制，故测量时间较长、测量效率较低等，不太适于快速测量的场合；③对大型物件，尤其是软质材料测量效果不好，对测头不能触及的表面无法测量，而且对测量环境要求较高，故其应用范围往往受到限制。总之，接触式测量难以满足快速、有效的测量需求。在申请号为03153504. 6的中国专利技术专利中，申请人提到了非接触式测量的几种方法光学传感器法、激光扫描法、立体视觉法、投影栅相位法等。随着计算机视觉和图像技术的发展，利用图像数据获取目标三维形状信息的非接触式测量技术已在社会生产的各个方面显示出越来越重要的地位和作用。由单幅图像明暗变化恢复形状(简称SFS)技术是实现物体表面三维重构的关键方法之一，其目标是利用单幅图像的明暗变化恢复物体表面各点的法向向量，进而可以获得其相对高度。许多研究者试图找出最有效的实现SFS技术的方法，并取得了一定的进展。目前应用最广泛的主要有四种方法演化和偏微分方程方法、最优化方法、局部化方法以及线性化方法等。演化和偏微分方程方法是从图像中一组已知高度值的边界点或奇异点出发，逐步确定出所有图像像素点对应的物体表面的形状信息，它是对SFS问题建立的一阶非线性偏微分方程本身直接求解。最优化方法就是将SFS问题中的图像辐照度方程表示为能量函数形式，然后通过附加约束条件将其转化为求解泛函极值的问题。局部化方法的思路就是在物体表面的每一点处，假定表面形状为平面、球面或抛物面等，并将表面反射模型与假设的局部形状结合起来，构成关于局部形状参数的线性偏微分方程组，然后利用已知的边界条件来求解该方程组得到物体表面的三维形状。线性化方法就是通过将反射图线性化，使非线性问题转化为线性问题，进而求解得到图像辐照度方程的解。在上述这些实现SFS技术的方法中，各种方法存在一定的优点，也不可避免地存在不足之处。如线性化方法速度较快，但重构精度不高，求解过程缺乏可信度；局部化方法使用范围受到限制，且精度不高；演化和偏微分方程方法虽然精度较高，但需要考虑边界条件；最优化方法不需要考虑边界条件，但需要引入各种约束条件等等。同时，上述方法通常只是针对朗伯表面设计的，如果应用到非朗伯表面，往往误差较大。另外，SFS方法中建立的图像辐照度方程通常是Hamilton-Jacobi类型的一阶非线性偏微分方程，一般情况下不存在经典意义下的光滑解，通常认为是病态问题。为了使SFS问题良性化，考虑图像辐照度方程在弱解意义下的特性。Hamilton-Jacobi方程粘性解理论是一种非常好的弱解理论，可定量、直观地描述一阶或二阶偏微分方程的弱解。法国学者Prados提出了使用粘性解理论对朗伯表面下建立的图像辐照度方程解的特性进行了分析。但遗憾的是，对于非朗伯表面建立的图像辐照度方程的研究，尚未发现国内外有相关文献报导。
技术实现思路
本专利技术的目的在于针对现有技术的不足，提出。本专利技术所述的方法不需要特定的边界条件，而且使SFS问题良性化，具有快速、准确的优点。为了实现上述目的，本专利技术采用的技术方案如下 ，依次包括下述步骤 一.构建适于非朗伯表面且接近实际成像条件下的成像模型从三个阶段入手，光源入射过程、非朗伯表面对入射光的反射过程、摄像机对反射光的记录过程，涉及到光源的分布、非朗伯表面的反射特性以及摄像机的投影方式的建模，其中所述光源为近点光源，而且光源辐照度随着物体与光源之间距离的增大而存在衰减；非朗伯表面的反射特性使用Oren-Nayer反射模型来近似；摄像机的投影方式采用透视投影方式。二.由上述成像模型和摄像机获取的灰度图像，建立非朗伯表面下SFS问题的图像辐照度方程该步骤依次包括下述步骤， (一)建立三维笛卡尔直角坐标系，设成像平面为X二 (X1, X2)位于Z = -/处，/ >0为摄像机的焦距，摄像机的光轴与Z轴重合,光源位于投影中心，即光心O处。非朗伯表面S可以表示为Z = Z(X), X=(Z11Z2)，像平面上一点(X,-/)是物体表面S上点P=在成像平面的投影点，r为光心O到P点的距离，根据透视投影成像原理有权利要求1.，依次包括下述步骤 一.构建适于非朗伯表面且接近实际成像条件下的成像模型从三个阶段入手，光源入射过程、非朗伯表面对入射光的反射过程、摄像机对反射光的记录过程，涉及到光源的分布、非朗伯表面的反射特性以及摄像机的投影方式的建模，其中所述光源为近点光源，而且光源辐照度随着物体与光源之间距离的增大而存在衰减；非朗伯表面的反射特性使用Oren-Nayer反射模型来近似；摄像机的投影方式采用透视投影方式； 二.由上述成像模型和摄像机获取的灰度图像，建立非朗伯表面下SFS问题的图像辐照度方程该步骤里依次包括下述步骤， (一)建立三维笛卡尔直角坐标系，设成像平面为X = ( , X2)位于Z = -/处，/ >0为摄像机的焦距，摄像机的光轴与Z轴重合,光源位于投影中心，即光心O处， 非朗伯表面S可以表示为z = Z(X), X=(X1J2)，像平面上一点(X,-/)是物体表面S上点P=在成像平面的投影点，r为光心O到P点的距离，根据透视投影成像原理有全文摘要本专利技术属于图像理解与计算机视觉
，特别涉及。提供的技术方案是，依次包括下述步骤一．构建适于非朗伯表面且接近实际成像条件下的成像模型；二．由上述成像模型和摄像机获取的灰度图像，建立非朗伯表面下SFS问题的图像辐照度方程；三．利用粘性解理论计算上述图像辐照度方程的粘性解；四．根据图像辐照度方程解的结果，重构出非朗伯表面的三维形状信息。本专利技术的优点是直接对非朗伯表面下建立的图像辐照度方程进行求解，同时利用弱解理论——粘性解理论计算图像辐照度方程的粘性解，使非朗伯表面的SFS问题良性化；具有快速、准确的优点。文档编号G06T15/00GK102819866SQ20121029869公开日2012年12月12日 申请日期2012年8月22日 优先权日2012年8月22日专利技术者王国珲, 程锦, 苏炜, 宋玉贵, 王建 申请人:西安工业大学本文档来自技高网...

【技术保护点】
一种非朗伯表面快速三维重构方法，依次包括下述步骤：一．构建适于非朗伯表面且接近实际成像条件下的成像模型：从三个阶段入手，光源入射过程、非朗伯表面对入射光的反射过程、摄像机对反射光的记录过程，涉及到光源的分布、非朗伯表面的反射特性以及摄像机的投影方式的建模，其中：所述光源为近点光源，而且光源辐照度随着物体与光源之间距离的增大而存在衰减；非朗伯表面的反射特性使用Oren？Nayer反射模型来近似；摄像机的投影方式采用透视投影方式；二．由上述成像模型和摄像机获取的灰度图像，建立非朗伯表面下SFS问题的图像辐照度方程：该步骤里依次包括下述步骤，(一)？建立三维笛卡尔直角坐标系，设成像平面为？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？位于处，为摄像机的焦距，摄像机的光轴与Z轴重合，光源位于投影中心，即光心O处,非朗伯表面S可以表示为，像平面上一点是物体表面S上点在成像平面的投影点，r为光心O到P点的距离，根据透视投影成像原理有？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？(1)令，于是物体表面S可以用函数表示？？？？？？？？？？？？？？？？？？？？？？？？？？(2)其中：是定义在实数集上的一个开集，代表图像的区域,(二)？计算光源入射到点的单位方向向量为？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？(3)(三)？利用Oren？Nayer反射模型计算非朗伯表面辐射亮度？？？？？？？？？？？？？？(4)其中：和是非朗伯表面某一点的局部坐标系，分别为光源方向向量L和摄像机方向向量V的天顶角和方位角；为Gauss分布函数的标准差，代表了表面的粗糙程度；为漫反射率；为点光源的辐射强度；,本专利技术中由于光源位于投影中心，故式（4）简化为？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？(5)其中：为摄像机获取的灰度图像；,(四)？计算摄像机遵循透视投影方式时，非朗伯表面点处的法向向量为？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？(6)(五)？由于为与之间的夹角，且令，故有？？？？？？(7)其中：将式（7）代入式（5），得到非朗伯表面下SFS问题的图像辐照度方程？？？？？？？？？？？？？？？？？？？？？？？？？？？(8)其中：为摄像机获取的灰度图像；三．利用粘性解理论计算上述图像辐照度方程的粘性解：该步骤依次包括下述步骤(一)？将图像辐照度方程（8）看作是关于的一元二次方程？？？？？？？？？？？？？？？？？？？？？？？？？？(9)(二)？求解一元二次方程（9）得？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？(10)(三)？于是，图像辐照度方程（8）等价为？？(11)(四)？一般情况下，方程（11）不存在经典解，可以利用粘性解理论计算上述图像辐照度方程的解,显然，上式是一个Hamilton？Jacobi类型的一阶非线性PDE，将简记为p，可以得到相应的Hamilton函数？(12)(五)？应用最优控制理论，式（12）可以转化为如下的控制形式？？？？？？(13)其中：；；是定义在上的单位圆面；D和R满足(六)？图像区域的均匀离散网格点：，，，其中：定义了数值算法中离散网格的大小,定义，对其在时域应用前向Euler公式得到图像辐照度方程求解算法？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？？(14)(七)？本专利技术中使用如下方法逼近？？？？(15)其中：为的分量；分别为p的后向、前向差分,(八)？对于逼近方法（15），计算以下最优问题：？？？？？？？？？？？？？？？(16)(九)？求解算法（14）转化为？？？？？？？？？？？？？？？？？？？？？？？(17)使用Newton法来进行求解，得到图像辐照度方程的粘性解,(十)？使用迭代Fast？Marching对上述求解算法（17）进行加速收敛,四．根据图像辐照度方程解的结果，重构出非朗伯表面的三维形状信息。2012102986957100001dest_path_image001.jpg,2012102986957100001dest_path_image002.jpg,2012102986957100001dest_path_image003.jpg,2012102986957100001dest_path_image004.j...

【技术特征摘要】

【专利技术属性】
技术研发人员：王国珲，程锦，苏炜，宋玉贵，王建，
申请(专利权)人：西安工业大学，
类型：发明
国别省市：