CN106952331A - 一种基于三维模型的纹理映射方法和装置 - Google Patents

一种基于三维模型的纹理映射方法和装置 Download PDF

Info

Publication number
CN106952331A
CN106952331A CN201710115610.XA CN201710115610A CN106952331A CN 106952331 A CN106952331 A CN 106952331A CN 201710115610 A CN201710115610 A CN 201710115610A CN 106952331 A CN106952331 A CN 106952331A
Authority
CN
China
Prior art keywords
characteristic point
value
point
free form
form surface
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
Application number
CN201710115610.XA
Other languages
English (en)
Other versions
CN106952331B (zh
Inventor
刘星明
蔡铁
王辉静
唐琪
闫立军
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shenzhen Institute of Information Technology
Original Assignee
Shenzhen Institute of Information Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shenzhen Institute of Information Technology filed Critical Shenzhen Institute of Information Technology
Priority to CN201710115610.XA priority Critical patent/CN106952331B/zh
Publication of CN106952331A publication Critical patent/CN106952331A/zh
Application granted granted Critical
Publication of CN106952331B publication Critical patent/CN106952331B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T15/003D [Three Dimensional] image rendering
    • G06T15/04Texture mapping

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Graphics (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Image Generation (AREA)
  • Image Analysis (AREA)

Abstract

本发明提供了一种基于三维模型的纹理映射方法和装置,纹理映射装置通过同时选取纹理图片和三维模型的特征点,使得通过对该特征点对进行奇异值分解所得到的自由曲面映射初始值的准确率大大提高,同时通过对二维图像特征点和三维几何特征点的数据进行归一化处理,减少了尺度因素对自由曲面映射初始值误差的影响,再将该自由曲面映射初始值作为迭代初值代入根据特征点所建立的误差距离函数中进行迭代计算优化后,得到精准的自由曲面映射最终值,再根据该自由曲面映射最终值,将纹理图片映射至根据三维模型数据采集系统所采集的点云数据所建立的三维自由曲面初始结构上,以生成精准的仿真三维模型。

Description

一种基于三维模型的纹理映射方法和装置
技术领域
本发明涉及本发明属于计算机视觉技术领域,尤其涉及一种基于三维模型的纹理映射的方法和装置。
背景技术
现实世界的三维模型通常无法用传统的规则多边形、多面体或球面进行表示,或者说无法用规则的多面体或球形精确的逼近曲面,而自由曲面具有不规则的形体(B样条曲面、隐式曲面和细分曲面等),能够实现任意曲面的表达。如何实现自由曲面的重建以及自由曲面的映射近年引起很多研究专家和学者的关注和研究,可以说,实现复杂结构的自由曲面表示以及自由曲面映射是一项十分重要且具有很大挑战的研究内容,具有十分重要的研究意义和应用价值。
现有的纹理映射方法中,是通过将采集纹理图片的相机固定在三维模型数据采集系统上,通过系统标定实现纹理图片和三维模型数据的同时采集,但在采集结构复杂的三维模型以得到自由曲面目标物体的数据时,由于纹理相机固定在三维模型数据采集系统上,而导致三维模型数据采集系统无法采集到完整的数据,使得根据所采集数据将纹理图片映射到三维模型所生成的仿真三维模型的精确度大大降低,无法实现建立精准的仿真三维模型。
发明内容
有鉴于此,本发明实施例提供了一种基于三维模型的纹理映射方法和装置,通过以解决现有纹理映射方法中因纹理图片不独立于三维模型数据采集系统,而导致不能精准地生成与复杂的三维模型相应的仿真三维模型的问题。
第一方面,本发明实施例提供了一种基于三维模型的纹理映射方法,所述方法包括:
获取所述三维模型的点云数据,识别并删除所述点云数据中散乱的点数据,建立三维自由曲面初始结构;
根据物像关系或光线映射关系交互选取纹理图片和所述三维自由曲面初始结构对应的特征点,得到预设数量的特征点对其中,为点集,点xi为所述纹理图片上的二维图像特征点,点Xi为所述三维自由曲面初始结构的三维几何特征点,k为所述特征点对的个数;
对所述预设数量的特征点对进行奇异值分解,得到所述三维模型的自由曲面映射初始值;
构建变换的误差距离函数其中,λi为尺度因子,Di为所述纹理图片中两个图像特征点之间的欧氏距离,P为待求自由曲面映射值;
将所述自由曲面映射初始值作为所述误差距离函数中P的初值,通过列文伯格-马夸尔特算法和预设的迭代终止阈值或迭代次数,对所述误差距离函数进行迭代计算得到所述三维模型的自由曲面映射最终值;
根据所述自由曲面映射最终值将所述纹理图片映射至所述三维自由曲面初始结构,生成仿真三维模型。
第二方面,本发明实施例提供了一种基于三维模型的纹理映射装置,所述装置包括三维自由曲面初始结构建立单元、特征点对选取单元、特征点对分解单元、误差距离函数构建单元、自由曲面映射最终值计算单元和纹理映射处理单元。
其中,三维自由曲面初始结构建立单元,用于获取所述三维模型的点云数据,识别并删除所述点云数据中散乱的点数据,建立三维自由曲面初始结构;
特征点对选取单元,用于根据物像关系或光线映射关系交互选取纹理图片和所述三维自由曲面初始结构对应的特征点,得到预设数量的特征点对其中,为点集,点xi为所述纹理图片上的二维图像特征点,点Xi为所述三维自由曲面初始结构的三维几何特征点,k为所述特征点对的个数;
特征点对分解单元,用于对所述预设数量的特征点对进行奇异值分解,得到所述三维模型的自由曲面映射初始值;
误差距离函数构建单元,用于构建变换的误差距离函数其中,λi为尺度因子,Di为所述纹理图片中两个图像特征点之间的欧氏距离,P为待求自由曲面映射值;
自由曲面映射最终值计算单元,用于将所述自由曲面映射初始值作为所述误差距离函数中P的初值,通过列文伯格-马夸尔特算法和预设的迭代终止阈值或迭代次数,对所述误差距离函数进行迭代计算得到所述三维模型的自由曲面映射最终值;
纹理映射处理单元,用于根据所述自由曲面映射最终值将所述纹理图片映射至所述三维自由曲面初始结构,生成仿真三维模型。
在本发明实施例中,纹理映射装置通过同时选取纹理图片和三维模型的特征点,使得通过对该特征点对进行奇异值分解所得到的自由曲面映射初始值的准确率大大提高,同时通过对二维图像特征点和三维几何特征点的数据进行归一化处理,减少了尺度因素对自由曲面映射初始值误差的影响,再将该自由曲面映射初始值作为迭代初值代入根据特征点所建立的误差距离函数中进行迭代计算优化后,得到精准的自由曲面映射最终值,再根据该自由曲面映射最终值,将纹理图片映射至根据三维模型数据采集系统所采集的点云数据所建立的三维自由曲面初始结构上,以生成精准的仿真三维模型。
附图说明
图1是本发明实施例提供的一种基于三维模型的纹理映射方法的实现流程图;
图2是本发明实施例提供的一种删除孤立的点数据的方法的具体实现流程图;
图3是本发明实施例提供的一种删除偏离的点数据的方法的具体实现流程图;
图4是本发明实施例提供的一种基于三维模型的纹理映射装置的结构框图。
具体实施方式
为使得本发明的发明目的、特征、优点能够更加的明显和易懂,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而非全部实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
应理解,在本发明各个实施例中,各过程的序号的大小并不意味着执行顺序的先后,各过程的执行顺序应以其功能和内在逻辑确定,而不应对本发明实施例的实施过程构成任何限定。
以下结合具体实施例对本发明的实现进行详细描述:
图1示出了本发明实施例提供的一种基于三维模型的纹理映射方法的实现流程,详述如下:
在步骤S101中,获取所述三维模型的点云数据,识别并删除所述点云数据中散乱的点数据,建立三维自由曲面初始结构。
在本发明实施例中,所述三维模型的点云数据为三维模型数据采集系统所采集的需要进行仿真建模的三维模型的点数据,该点云数据用于构建三维自由曲面初始结构。优先的,为了避免数据采集过程中引入的散乱的点数据造成的噪声干扰,如三维模型中的光斑投影或反射投影,在根据该点云数据建立三维自由曲面初始结构之前,还包括识别并删除该点云数据中散乱的点数据,以使得所建立的三维自由曲面初始结构更加精准。
在这里,在建立三维自由曲面初始结构之前,还包括建立第一坐标系,在所述第一坐标系上根据点云数据构建三维自由曲面初始结构。
在步骤S102中,根据物像关系或光线映射关系交互选取纹理图片和所述三维自由曲面初始结构对应的特征点,得到预设数量的特征点对其中,为点集,点xi为所述纹理图片上的二维图像特征点,点Xi为所述三维自由曲面初始结构的三维几何特征点,k为所述特征点对的个数。
在本发明实施例中,所采用的纹理图片为独立于三维模型数据采集系统的纹理图片,即纹理图片与三维模型数据不是由同一采集系统同时进行采集的,该纹理图片可以是由纹理相机采集的纹理图片,也可以对原始图片进行纹理处理后的图片,这里不做具体限定。由于所采用的纹理图片独立三维模型数据采集系统,因此使得所采用的纹理图片或三维模型各自所采集的数据更加完整,从而使得根据这些数据所生成的仿真三维模型更加准确。本发明实施例中,通过将纹理图片精准地映射到三维自由曲面初始结构上以得到精准的仿真三维模型。
在这里,所采用的纹理图片为二维图像,从纹理图片上选取的特征点为二维图像特征点xi。在这里,二维图像特征点xi为纹理图像中含图像特征的点,优选为特征圆点、角点或其他具有明显色彩的特征点。三维自由曲面初始结构为三维几何图像,从三维自由曲面初始结构上选取的特征点为三维几何特征点Xi。在这里,三维几何特征点为三维自由曲面初始结构中含几何特征的点,优选为曲线或曲面的特征点、几何角点、几何特征点。所选取的二维图像特征点xi和三维几何特征点Xi存在对应关系。在这里,在得到预设数量的特征点对后,建立第二坐标系,并将预设数量的特征点对对应标识在该第二坐标系上。
在这里,预设数量对应为特征点对的k,其中,k为用户输入的数值,一般取k≥6,以避免因特征点对所选取的个数少于小孔成像方程式未知数的个数而导致无法求解得到自由曲面映射初始值p0。如果所选取的特征点对的个数超过小孔成像方程式未知数的个数,或远远超过小孔成像方程式未知数的个数,可根据最小二乘法求解得到自由曲面映射初始值p0
在这里,点云数据中散乱的点数据为孤立和/或偏离的点数据。三维几何特征点Xi存在有m个与其有连接关系的邻近点Xj,Xi的m个邻近点Xj表示为其中,当三维几何特征点Xi与其m个邻近点Xj的平均距离大于预设的距离阈值时,该三维几何特征点Xi为孤立的点数据;当三维几何特征点Xi与其m个邻近点Xj的法向量夹角大于预设角度时,该三维几何特征点Xi为偏离的点数据。
具体地,所述点云数据中散乱的点数据为孤立的点数据时,步骤S101中的识别并删除所述点云数据中散乱的点数据的步骤中,包括如图2所示的删除孤立的点数据的具体实现流程:
在步骤S201中,识别并获取与所述点云数据中的三维几何特征点Xi存在连接关系的m个邻近点Xj,所述m为正整数。
在本发明实施例中,将点云数据中的任意点数据表示为三维几何特征点Xi,将与每个三维几何特征点Xi存在直接的连接关系的点数据作为三维几何特征点Xi的邻近点Xj
在步骤S202中,根据平均距离函数计算得到所述三维几何特征点Xi与m个所述邻近点Xj的平均距离。
在本发明实施例中,先分别计算三维几何特征点Xi与m个三邻近点Xj的距离,将计算得到的距离相加后除以m,得到三维几何特征点Xi与m个邻近点Xj的平均距离,比如,分别计算三维几何特征点Xi与m个邻近点Xj得到的距离为a1、a2…am,那么,三维几何特征点Xi与m个邻近点Xj的距离的平均距离为a1、a2…am相加后得到的和乘以的值。
在步骤S203中,确定所述平均距离函数的值是否大于预设的距离阈值。
在本发明实施例中,通过确定平均距离函数的值是否大于预设的距离阈值来判断三维几何特征Xi是否为孤立的点数据。
在步骤S204中,如果所述平均距离函数的值大于预设的距离阈值,删除所述三维几何特征点Xi
或者,当无法计算三维几何特征点Xi与m个邻近点Xj的平均距离时,通过判断三维几何特征点的法向量夹角与预设角度的关系来确定是否删除所述三维几何特征点Xi,即所述点云数据中散乱的点数据为偏离的点数据时,步骤S101中的识别并删除所述点云数据中散乱的点数据的步骤中,还包括如图3所示的删除孤立的点数据的具体实现流程:
在步骤S301中,识别并获取与所述点云数据中的三维几何特征点Xi存在连接关系的m个邻近点Xj,所述m为正整数。
在步骤S302中,确定所述三维几何特征点Xi与m个所述邻近点Xj的法向量夹角是否大于预设角度。
在本发明实施例中,通过确定三维几何特征点Xi与m个邻近点Xj的法向量夹角是否大于预设角度来确定该三维几何特征点Xi是否为偏离的点数据。
在步骤S303中,如果所述三维几何特征点Xi与m个所述邻近点Xj的法向量夹角大于预设角度,删除所述三维几何特征点Xi
在本发明实施例中,预设角度是由用户进行设置的,可以设定为任意角度值。优选的,预设角度为65°,即当三维几何特征点Xi与其邻近点的法向量夹角θ≥65°时,则认为该三维几何特征点Xi为偏离的点数据。
在这里,对孤立的点数据或者偏离的点数据的识别及删除过程不分先后,可以同时进行,也可以单独进行,这里不做具体限定。
在这里,用户根据点云数据中各个点数据周围的连接关系,选取适用该点云数据的网格模型,优选结构简单,便于表达的三维点云网格模型。在这里,自由曲面可以由样条曲面、隐式曲面表示。
具体地,在确定使用的网格模型后,纹理映射装置根据去除散乱的点数据后的点云数据,利用特征曲线和网格曲面进行逼近拟合,以建立三维自由曲面初始结构。
在步骤S103中,对所述预设数量的特征点对进行奇异值分解,得到所述三维模型的自由曲面映射初始值p0
在本发明实施例中,在对所述预设数量的特征点对进行奇异值分解之前,还包括对所述预设数量的特征点对进行归一化处理,以得到满足预设数值要求的特征点对。
具体地,步骤S103通过奇异值分解以得到三维模型的自由曲面映射初始值p0的具体实现步骤如图4所示:
在步骤S401中,对所述特征点对中的进行归一化处理,使得二维图像特征点xi距离中心坐标的平均距离满足第一预设数值要求。
在本发明实施例中,中心坐标为根据二维图像特征点xi的各个坐标之和计算出来的平均坐标,比如,分别将k个X轴的坐标和k个Y轴的坐标相加后求得的平均值作为由X轴和Y轴组成的二维平面上的中心坐标[(x1+x2+…+xk)/k,(y1+y2+…+yk)]。在这里,在对进行归一化处理后,二维图像特征点xi距离中心坐标的平均距离为(即为第一预设数据要求),或是为其他预设的小范围数值。
在步骤S402中,对所述特征点对中的进行归一化处理,使得Xi距离中心坐标的平均距离满足第二预设数值要求。
在本发明实施例中,中心坐标为根据三维几何特征点Xi的各个坐标之和计算出来的平均坐标,比如,分别将k个X轴的坐标、k个Y轴的坐标和k个Z轴的坐标相加后求得的平均值作为由X轴、Y轴和Z轴组成的三维平面上的中心坐标[(x1+x2+…+xk)/k,(y1+y2+…+yk),(z1+z2+…+zk)]。在这里,在对进行归一化处理后,三维几何特征点Xi距离中心坐标的平均距离为(即为第二预设数据要求),或是为其他预设的小范围数值。
在步骤S403中,对进行归一化处理后的所述特征点对进行奇异值分解,将奇异值矩阵中最小的奇异值对应的右奇异向量V中的奇异向量作为所述三维模型的自由曲面映射初始值P0
在本发明实施例中,通过对特征点对进行归一化处理,以减少尺度因素对数据的影响,降低自由曲面映射初始值P0的计算误差。在这里,尺度因素为影响奇异值分解结果的数据尺度。
在步骤S104中,构建变换的误差距离函数其中,λi为尺度因子,Di为所述纹理图片中两个图像特征点之间的欧氏距离,P为待求自由曲面映射值。
在本发明实施例中,λi的取值与曲面空间相关,0.8≤λi≤1,λi由用户根据曲面的结构来设置,其中,处于曲面边缘的数据点设置λi=0.8,处于数据中心的点设置的λi=1.0,λi∈[0.8,1.0]。在这里,处于数据中心的点为根据自由曲面上所有点的均值计算得到的数据点。
在这里,二维图像特征点xi的齐次坐标表示为:xi=(xi1,yi1,1)T,三维几何特征点Xi的齐次坐标表示为:Xi=(Xi1,Yi1,Zi1,1)T
在这里,D为纹理图片中两个图像特征点之间的欧氏距离,具体为
其中,x11和x12分别为x1在X轴和Y轴上的坐标点,x21和x22分别为x2在X轴和Y轴上的坐标点,D(x1,x2)为纹理图片上x1和x2之间的欧氏距离。
在步骤S105中,将所述自由曲面映射初始值P0作为所述误差距离函数中P的初值,通过列文伯格-马夸尔特算法和预设的迭代终止阈值或迭代次数,对所述误差距离函数进行迭代计算得到所述三维模型的自由曲面映射最终值。
在本发明实施例中,通过列文伯格-马夸尔特算法进行算法迭代优化,以求解得到更加精确的待求自由曲面映射值P,以得到精准的自由曲面映射最终值。
具体地,对所述误差距离函数进行迭代计算得到所述三维模型的自由曲面映射最终值的步骤,包括:
对所述误差函数进行迭代计算出所述待求自由曲面映射值P,对计算出的所述待求自由曲面映射值P进行坐标变换后,再对所述待求自由曲面映射值P进行逆归一化处理,得到所述三维模型的自由曲面映射最终值。
优选的,在对所述误差距离函数进行迭代计算得到所述三维模型的自由曲面映射最终值的步骤中,还包括:
在每一次进行迭代计算之前,识别并删除距离过大的特征点对,再对所述误差距离函数进行下一次迭代计算。
通过在每次迭代计算前删除距离过大的特征点对后,再进行下一步的迭代优化计算,以保证迭代算法的收敛,提高迭代计算的效率。
在步骤S106中,根据所述自由曲面映射最终值将所述纹理图片映射至所述三维自由曲面初始结构,生成仿真三维模型。
在本发明实施例中,通过同时选取纹理图片和三维模型的特征点,使得通过对该特征点对进行奇异值分解所得到的自由曲面映射初始值p0的准确率大大提高,同时通过对二维图像特征点和三维几何特征点的数据进行归一化处理,减少了尺度因素对自由曲面映射初始值p0误差的影响,再将该自由曲面映射初始值p0作为迭代初值代入根据特征点所建立的误差距离函数中的p进行迭代计算优化后,得到精准的自由曲面映射最终值,再根据该自由曲面映射最终值,将纹理图片映射至根据三维模型数据采集系统所采集的点云数据所建立的三维自由曲面初始结构上,以生成精准的仿真三维模型。
本领域普通技术人员可以理解实现上述各实施例方法中的全部或部分步骤是可以通过程序来指令相关的硬件来完成,相应的程序可以存储于一计算机可读取存储介质中,上述的存储介质,如ROM/RAM、磁盘或光盘等。
图4示出了本发明实施例提供的一种基于三维模型的纹理映射装置41的具体结构框图,为了便于说明,仅示出了与本发明实施例相关的部分。
本发明实施例提供了一种基于三维模型的纹理映射装置4,包括三维自由曲面初始结构建立单元41、特征点对选取单元42、特征点对分解单元43、误差距离函数构建单元44、自由曲面映射最终值计算单元45和纹理映射处理单元46;
其中,三维自由曲面初始结构建立单元41,用于获取所述三维模型的点云数据,识别并删除所述点云数据中散乱的点数据,建立三维自由曲面初始结构;
特征点对选取单元42,用于根据物像关系或光线映射关系交互选取纹理图片和所述三维自由曲面初始结构对应的特征点,得到预设数量的特征点对其中,为点集,点xi为所述纹理图片上的二维图像特征点,点Xi为所述三维自由曲面初始结构的三维几何特征点,k为所述特征点对的个数;
特征点对分解单元43,用于对所述预设数量的特征点对进行奇异值分解,得到所述三维模型的自由曲面映射初始值p0
误差距离函数构建单元44,用于构建变换的误差距离函数其中,λi为尺度因子,Di为所述纹理图片中两个图像特征点之间的欧氏距离,P为待求自由曲面映射值;
自由曲面映射最终值计算单元45,用于将所述自由曲面映射初始值p0作为所述误差距离函数中P的初值,通过列文伯格-马夸尔特算法和预设的迭代终止阈值或迭代次数,对所述误差距离函数进行迭代计算得到所述三维模型的自由曲面映射最终值;
纹理映射处理单元46,用于根据所述自由曲面映射最终值将所述纹理图片映射至所述三维自由曲面初始结构,生成仿真三维模型。
进一步地,所述散乱的点数据为孤立的点数据时;所述三维自由曲面初始结构建立单元中包括第一邻近点识别子单元、平均距离计算子单元、距离确定子单元和第一删除子单元;
其中,第一邻近点识别子单元,用于识别并获取与所述点云数据中的所述三维几何特征点Xi存在连接关系的m个邻近点Xj,所述m为正整数;
平均距离计算子单元,用于根据平均距离函数计算得到所述三维几何特征点Xi与m个所述邻近点Xj的平均距离;
距离确定子单元,用于确定所述平均距离函数的值是否大于预设的距离阈值;
第一删除子单元,用于如果所述平均距离函数的值大于预设的距离阈值,删除所述三维几何特征点Xi
或者,所述散乱的点数据为偏离的点数据;所述三维自由曲面初始结构建立单元中还包括第二邻近点识别子单元、夹角确定子单元和第二删除子单元;
其中第二邻近点识别子单元,用于识别并获取与所述点云数据中的所述三维几何特征点Xi存在连接关系的m个邻近点Xj,所述m为正整数;
夹角确定子单元,用于确定所述三维几何特征点Xi与m个所述邻近点Xj的法向量夹角是否大于预设角度;
第二删除子单元,如果所述三维几何特征点Xi与m个所述三维几何特征点Xj的法向量夹角大于预设角度,删除所述三维几何特征点Xi
进一步地,所述特征点对分解单元43,包括第一归一化处理子单元、第一归二化处理子单元和特征点对分解子单元;
其中,第一归一化处理子单元,用于对所述特征点对中的进行归一化处理,使得xi距离中心坐标的平均距离满足第一预设数值要求;
第二归一化处理子单元,用于对所述特征点对中的进行归一化处理,使得Xi距离中心坐标的平均距离满足第二预设数值要求;
特征点对分解子单元,用于对进行归一化处理后的所述特征点对进行奇异值分解,将奇异值矩阵中最小的奇异值对应的右奇异向量V中的奇异向量作为所述三维模型的自由曲面映射初始值p0
进一步地,所述自由曲面映射最终值计算单元45,具体用于:
对所述误差函数进行迭代计算出所述待求自由曲面映射值P,对计算出的所述待求自由曲面映射值P进行坐标变换后,再对所述待求自由曲面映射值P进行逆归一化处理,得到所述三维模型的自由曲面映射最终值。
进一步地,所述自由曲面映射最终值计算单元45,还包括:
特征点对删除子单元,用于在每一次进行迭代计算之前,识别并删除距离过大的特征点对,再对所述误差距离函数进行下一次迭代计算。
在本发明实施例中,纹理映射装置通过同时选取纹理图片和三维模型的特征点,使得通过对该特征点对进行奇异值分解所得到的自由曲面映射初始值p0的准确率大大提高,同时通过对二维图像特征点和三维几何特征点的数据进行归一化处理,减少了尺度因素对自由曲面映射初始值p0误差的影响,再将该自由曲面映射初始值p0作为迭代初值代入根据特征点所建立的误差距离函数中的p进行迭代计算优化后,得到精准的自由曲面映射最终值,再根据该自由曲面映射最终值,将纹理图片映射至根据三维模型数据采集系统所采集的点云数据所建立的三维自由曲面初始结构上,以生成精准的仿真三维模型。
本领域普通技术人员可以意识到,结合本文中所公开的实施例描述的各示例的单元及步骤,能够以电子硬件、或者计算机软件和电子硬件的结合来实现。这些功能究竟以硬件还是软件方式来执行,取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用来使用不同装置来实现所描述的功能,但是这种实现不应认为超出本发明的范围。
所属领域的技术人员可以清楚地了解到,为描述的方便和简洁,上述描述的装置和单元的具体工作过程,可以参考前述装置实施例中的对应过程,在此不再赘述。
在本申请所提供的几个实施例中,应该理解到,所揭露的装置和方法,可以通过其它的方式实现。例如,以上所描述的装置实施例仅仅是示意性的,例如,上述单元的划分,仅仅为一种逻辑功能划分,实际实现时可以有另外的划分方式,例如多个单元或组件可以结合或者可以集成到另一个系统,或一些特征可以忽略,或不执行。另一点,所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口,装置或单元的间接耦合或通信连接,可以是电性,机械或其它的形式。
上述作为分离部件说明的单元可以是或者也可以不是物理上分开的,作为单元显示的部件可以是或者也可以不是物理单元,即可以位于一个地方,或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。
另外,在本发明各个实施例中的各功能单元可以集成在一个处理单元中,也可以是各个单元单独物理存在,也可以两个或两个以上单元集成在一个单元中。
上述功能如果以软件功能单元的形式实现并作为独立的产品销售或使用时,可以存储在一个计算机可读取存储介质中。基于这样的理解,本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质中,包括若干指令用以使得一台计算机设备(可以是个人计算机,服务器,或者网络设备等)执行本发明各个实施例上述装置的全部或部分步骤。而前述的存储介质包括:U盘、移动硬盘、只读存储器(ROM,Read-Only Memory)、随机存取存储器(RAM,Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
以上上述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到状态或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应上述以权利要求的保护范围为准。

Claims (10)

1.一种基于三维模型的纹理映射方法,其特征在于,所述方法包括:
获取所述三维模型的点云数据,识别并删除所述点云数据中散乱的点数据,建立三维自由曲面初始结构;
根据物像关系或光线映射关系交互选取纹理图片和所述三维自由曲面初始结构对应的特征点,得到预设数量的特征点对其中,为点集,点xi为所述纹理图片上的二维图像特征点,点Xi为所述三维自由曲面初始结构的三维几何特征点,k为所述特征点对的个数;
对所述预设数量的特征点对进行奇异值分解,得到所述三维模型的自由曲面映射初始值;
构建变换的误差距离函数其中,λi为尺度因子,Di为所述纹理图片中两个图像特征点之间的欧氏距离,P为待求自由曲面映射值;
将所述自由曲面映射初始值作为所述误差距离函数中P的初值,通过列文伯格-马夸尔特算法和预设的迭代终止阈值或迭代次数,对所述误差距离函数进行迭代计算得到所述三维模型的自由曲面映射最终值;
根据所述自由曲面映射最终值将所述纹理图片映射至所述三维自由曲面初始结构,生成仿真三维模型。
2.如权利要求1所述的方法,其特征在于,所述散乱的点数据为孤立的点数据;
所述识别并删除所述点云数据中散乱的点数据的步骤,具体包括:
识别并获取与所述点云数据中的所述三维几何特征点Xi存在连接关系的m个邻近点Xj,所述m为正整数;
根据平均距离函数计算得到所述三维几何特征点Xi与m个所述邻近点Xj的平均距离;
确定所述平均距离函数的值是否大于预设的距离阈值;
如果所述平均距离函数的值大于预设的距离阈值,删除所述三维几何特征点Xi
或者,所述散乱的点数据为偏离的点数据;
所述识别并删除所述点云数据中散乱的点数据的步骤,具体包括:
识别并获取与所述点云数据中的所述三维几何特征点Xi存在连接关系的m个邻近点Xj,所述m为正整数;
确定所述三维几何特征点Xi与m个所述邻近点Xj的法向量夹角是否大于预设角度;
如果所述三维几何特征点Xi与m个所述三维几何特征点Xj的法向量夹角大于预设角度,删除所述三维几何特征点Xi
3.如权利要求1或2所述的方法,其特征在于,在所述对所述预设数量的特征点对进行奇异值分解,得到所述三维模型的自由曲面映射初始值的步骤,包括:
对所述特征点对中的进行归一化处理,使得xi距离中心坐标的平均距离满足第一预设数值要求;
对所述特征点对中的进行归一化处理,使得Xi距离中心坐标的平均距离满足第二预设数值要求;
对进行归一化处理后的所述特征点对进行奇异值分解,将奇异值矩阵中最小的奇异值对应的右奇异向量V中的奇异向量作为所述三维模型的自由曲面映射初始值。
4.如权利要求1所述的方法,其特征在于,所述对所述误差距离函数进行迭代计算得到所述三维模型的自由曲面映射最终值的步骤,包括:
对所述误差函数进行迭代计算出所述待求自由曲面映射值P,对计算出的所述待求自由曲面映射值P进行坐标变换后,再对所述待求自由曲面映射值P进行逆归一化处理,得到所述三维模型的自由曲面映射最终值。
5.如权利要求1所述的方法,其特征在于,所述对所述误差距离函数进行迭代计算得到所述三维模型的自由曲面映射最终值的步骤中,还包括:
在每一次进行迭代计算之前,识别并删除距离过大的特征点对,再对所述误差距离函数进行下一次迭代计算。
6.一种基于三维模型的纹理映射装置,其特征在于,所述装置包括:
三维自由曲面初始结构建立单元,用于获取所述三维模型的点云数据,识别并删除所述点云数据中散乱的点数据,建立三维自由曲面初始结构;
特征点对选取单元,用于根据物像关系或光线映射关系交互选取纹理图片和所述三维自由曲面初始结构对应的特征点,得到预设数量的特征点对其中,为点集,点xi为所述纹理图片上的二维图像特征点,点Xi为所述三维自由曲面初始结构的三维几何特征点,k为所述特征点对的个数;
特征点对分解单元,用于对所述预设数量的特征点对进行奇异值分解,得到所述三维模型的自由曲面映射初始值;
误差距离函数构建单元,用于构建变换的误差距离函数其中,λi为尺度因子,Di为所述纹理图片中两个图像特征点之间的欧氏距离,P为待求自由曲面映射值;
自由曲面映射最终值计算单元,用于将所述自由曲面映射初始值作为所述误差距离函数中P的初值,通过列文伯格-马夸尔特算法和预设的迭代终止阈值或迭代次数,对所述误差距离函数进行迭代计算得到所述三维模型的自由曲面映射最终值;
纹理映射处理单元,用于根据所述自由曲面映射最终值将所述纹理图片映射至所述三维自由曲面初始结构,生成仿真三维模型。
7.如权利要求6所述的装置,其特征在于,所述散乱的点数据为孤立的点数据;所述三维自由曲面初始结构建立单元中包括:
第一邻近点识别子单元,用于识别并获取与所述点云数据中的所述三维几何特征点Xi存在连接关系的m个邻近点Xj,所述m为正整数;
平均距离计算子单元,用于根据平均距离函数计算得到所述三维几何特征点Xi与m个所述邻近点Xj的平均距离;
距离确定子单元,用于确定所述平均距离函数的值是否大于预设的距离阈值;
第一删除子单元,用于如果所述平均距离函数的值大于预设的距离阈值,删除所述三维几何特征点Xi
或者,所述散乱的点数据为偏离的点数据;所述三维自由曲面初始结构建立单元中还包括:
第二邻近点识别子单元,用于识别并获取与所述点云数据中的所述三维几何特征点Xi存在连接关系的m个邻近点Xj,所述m为正整数;
夹角确定子单元,用于确定所述三维几何特征点Xi与m个所述邻近点Xj的法向量夹角是否大于预设角度;
第二删除子单元,如果所述三维几何特征点Xi与m个所述三维几何特征点Xj的法向量夹角大于预设角度,删除所述三维几何特征点Xi
8.如权利要求6或7所述的装置,其特征在于,所述特征点对分解单元包括:
第一归一化处理子单元,用于对所述特征点对中的进行归一化处理,使得xi距离中心坐标的平均距离满足第一预设数值要求;
第二归一化处理子单元,用于对所述特征点对中的进行归一化处理,使得Xi距离中心坐标的平均距离满足第二预设数值要求;
特征点对分解子单元,用于对进行归一化处理后的所述特征点对进行奇异值分解,将奇异值矩阵中最小的奇异值对应的右奇异向量V中的奇异向量作为所述三维模型的自由曲面映射初始值。
9.如权利要求6所述的装置,其特征在于,所述自由曲面映射最终值计算单元具体用于:
对所述误差函数进行迭代计算出所述待求自由曲面映射值P,对计算出的所述待求自由曲面映射值P进行坐标变换后,再对所述待求自由曲面映射值P进行逆归一化处理,得到所述三维模型的自由曲面映射最终值。
10.如权利要求6所述的装置,其特征在于,所述自由曲面映射最终值计算单元还包括:
特征点对删除子单元,用于在每一次进行迭代计算之前,识别并删除距离过大的特征点对,再对所述误差距离函数进行下一次迭代计算。
CN201710115610.XA 2017-02-28 2017-02-28 一种基于三维模型的纹理映射方法和装置 Expired - Fee Related CN106952331B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710115610.XA CN106952331B (zh) 2017-02-28 2017-02-28 一种基于三维模型的纹理映射方法和装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710115610.XA CN106952331B (zh) 2017-02-28 2017-02-28 一种基于三维模型的纹理映射方法和装置

Publications (2)

Publication Number Publication Date
CN106952331A true CN106952331A (zh) 2017-07-14
CN106952331B CN106952331B (zh) 2020-12-08

Family

ID=59466610

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710115610.XA Expired - Fee Related CN106952331B (zh) 2017-02-28 2017-02-28 一种基于三维模型的纹理映射方法和装置

Country Status (1)

Country Link
CN (1) CN106952331B (zh)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109472741A (zh) * 2018-09-30 2019-03-15 先临三维科技股份有限公司 三维拼接方法及装置
CN111369659A (zh) * 2018-12-26 2020-07-03 杭州海康威视数字技术股份有限公司 一种基于三维模型的纹理映射方法、装置及设备
CN111657858A (zh) * 2019-03-07 2020-09-15 株式会社日立制作所 图像诊断装置、图像处理方法和程序
CN112381823A (zh) * 2021-01-13 2021-02-19 北京智拓视界科技有限责任公司 用于图像的几何特征的提取方法和相关产品
CN115205562A (zh) * 2022-07-22 2022-10-18 四川云数赋智教育科技有限公司 一种基于特征点的任意试卷配准方法
WO2023029706A1 (zh) * 2021-09-02 2023-03-09 山东建筑大学 一种塔机塔身损伤状态特征点映射模型的建立方法与快速判别损伤的方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100245374A1 (en) * 2009-03-24 2010-09-30 Advanced Micro Devices, Inc. Method and apparatus for angular invariant texture level of detail generation
CN102944174A (zh) * 2012-11-28 2013-02-27 北京矿冶研究总院 一种点云数据的处理方法及系统
CN103824270A (zh) * 2013-09-25 2014-05-28 浙江树人大学 一种快速的离散三维点云滤波方法
CN104200507A (zh) * 2014-08-12 2014-12-10 南京理工大学 一种三维点云点法向量估计方法
US9280848B1 (en) * 2011-10-24 2016-03-08 Disney Enterprises Inc. Rendering images with volumetric shadows using rectified height maps for independence in processing camera rays

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100245374A1 (en) * 2009-03-24 2010-09-30 Advanced Micro Devices, Inc. Method and apparatus for angular invariant texture level of detail generation
US9280848B1 (en) * 2011-10-24 2016-03-08 Disney Enterprises Inc. Rendering images with volumetric shadows using rectified height maps for independence in processing camera rays
CN102944174A (zh) * 2012-11-28 2013-02-27 北京矿冶研究总院 一种点云数据的处理方法及系统
CN103824270A (zh) * 2013-09-25 2014-05-28 浙江树人大学 一种快速的离散三维点云滤波方法
CN104200507A (zh) * 2014-08-12 2014-12-10 南京理工大学 一种三维点云点法向量估计方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
刘星明 等: "《基于相位恢复的三维可视化应用研究》", 《计算机科学》 *
薛耀红: "《点云数据配准及曲面细分技术研究》", 《中国博士学位论文全文数据库 信息科技辑》 *

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109472741A (zh) * 2018-09-30 2019-03-15 先临三维科技股份有限公司 三维拼接方法及装置
CN111369659A (zh) * 2018-12-26 2020-07-03 杭州海康威视数字技术股份有限公司 一种基于三维模型的纹理映射方法、装置及设备
CN111369659B (zh) * 2018-12-26 2023-05-02 杭州海康威视数字技术股份有限公司 一种基于三维模型的纹理映射方法、装置及设备
CN111657858A (zh) * 2019-03-07 2020-09-15 株式会社日立制作所 图像诊断装置、图像处理方法和程序
CN111657858B (zh) * 2019-03-07 2023-08-01 株式会社日立制作所 图像诊断装置、图像处理方法和存储介质
CN112381823A (zh) * 2021-01-13 2021-02-19 北京智拓视界科技有限责任公司 用于图像的几何特征的提取方法和相关产品
CN112381823B (zh) * 2021-01-13 2021-04-16 北京智拓视界科技有限责任公司 用于图像的几何特征的提取方法和相关产品
WO2023029706A1 (zh) * 2021-09-02 2023-03-09 山东建筑大学 一种塔机塔身损伤状态特征点映射模型的建立方法与快速判别损伤的方法
CN115205562A (zh) * 2022-07-22 2022-10-18 四川云数赋智教育科技有限公司 一种基于特征点的任意试卷配准方法
CN115205562B (zh) * 2022-07-22 2023-03-14 四川云数赋智教育科技有限公司 一种基于特征点的任意试卷配准方法

Also Published As

Publication number Publication date
CN106952331B (zh) 2020-12-08

Similar Documents

Publication Publication Date Title
CN106952331A (zh) 一种基于三维模型的纹理映射方法和装置
CN100562895C (zh) 一种基于区域分割和分段学习的三维人脸动画制作的方法
CN103247041B (zh) 一种基于局部采样的多几何特征点云数据的分割方法
CN109934065A (zh) 一种用于手势识别的方法和装置
CN108171780A (zh) 一种基于激光雷达构建室内真实三维地图的方法
CN103701466A (zh) 基于特征保留的散乱点云压缩算法
CN108151679A (zh) 一种土地面积测量的方法及系统
CN108763827B (zh) 一种输电塔有限元模型建立方法及装置
CN106846485A (zh) 一种室内三维建模方法及装置
CN105261062B (zh) 一种人物分段建模方法
EP3846135A1 (en) Method and device for generating body leaf-shaped structure, and computer readable storage medium
CN106981097A (zh) 一种基于分区局部光顺权因子的t样条曲面拟合方法
CN104952086B (zh) 用于管理结构数据的设备及方法
CN110084894A (zh) 三维模型的局部放大展示方法、装置与电子设备
CN112562085A (zh) 一种bim模型简化方法及简化装置
CN106802958B (zh) Cad数据到gis数据的转换方法及系统
Bertault et al. Drawing hypergraphs in the subset standard (short demo paper)
Danilov Unstructured tetrahedral mesh generation technology
CN109684517A (zh) 一种历史数据存储方法、读写方法、存储装置和设备
CN107506572B (zh) 获取目标点的高度的方法和装置
Zhou 3D urban modeling from city-scale aerial LiDAR data
Cascón et al. A new meccano technique for adaptive 3-D triangulations
CN107895398B (zh) 结合视觉显著性的浮雕编辑方法
CN116258840A (zh) 层级细节表示树的生成方法、装置、设备及存储介质
Strodthoff et al. Layered Reeb graphs for three-dimensional manifolds in boundary representation

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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20201208

CF01 Termination of patent right due to non-payment of annual fee