CN100538745C - 颅骨三维模型构建方法 - Google Patents
颅骨三维模型构建方法 Download PDFInfo
- Publication number
- CN100538745C CN100538745C CNB2008100172233A CN200810017223A CN100538745C CN 100538745 C CN100538745 C CN 100538745C CN B2008100172233 A CNB2008100172233 A CN B2008100172233A CN 200810017223 A CN200810017223 A CN 200810017223A CN 100538745 C CN100538745 C CN 100538745C
- Authority
- CN
- China
- Prior art keywords
- skull
- point
- experimenter
- image
- eye socket
- 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
Images
Landscapes
- Image Processing (AREA)
Abstract
本发明公开了一种颅骨三维模型构建方法,通过采集受试者正位、左侧位和颅顶颏底位相互正交的三张X光头颅定位片图像,利用X线片图像上耳眼眶平面的四个点,对X光头颅定位片图像进行标定;根据两组耳眼眶平面点对颅骨标准模板进行缩放来获得受试者颅骨初始三维数据模型和初始特征点三维坐标;再根据初始特征点和受试者实际特征点采用基于径向基函数的形变算法进行受试者颅骨模型的构建,得到了完整的颅骨三维模型。
Description
技术领域
本发明涉及一种颅骨三维模型构建方法,特别是人类颅骨三维模型的构建方法。
背景技术
人类颅骨三维模型的构建,是口腔医学中一种重要的研究手段。
文献“X线片的颅颌骨特征点的三维确定,计算机辅助设计与图形学学报,2006,Vol.18(1),p89-93”公开了一种利用X线片确定颅颌骨特征点的方法。该方法根据立体头影测量的原理,针对颅颌骨特征点,在头颅正位、侧位定位片的基础上增加2张斜位定位片,以侧位定位片为0°,拍摄0°,45°,90°和135°4张X线头颅定位片。每次选取特征点显影最佳的2张,分别确定特征点的二维坐标并且对每对数据进行了校正,然后从中重构出颅颌骨的特征点。但是此方法使用图像处理的方法来获取4张X线片上的特征点,无法保证每个点在所有X线片上的投影是一一对应的;由于4张X线片并非正交关系,难以保证它们中心转轴一致,而特征点重建质量受拍摄中心转轴的偏移影响较大;仅能重构出少量特征点的三维坐标,无法得到完整颅骨的三维模型。
发明内容
为了克服现有技术无法得到完整的颅骨三维模型的不足,本发明提供一种颅骨三维模型构建方法,采用三张正交X线片图像,通过耳眼眶平面四个点对正交X片图像进行标定,根据两组耳眼眶平面点对颅骨标准模板进行缩放来获得颅骨初始模型,再根据初始模型特征点和实际特征点采用基于径向基函数的形变算法进行颅骨模型的构建,可以得到完整的颅骨三维模型。
本发明解决其技术问题所采用的技术方案:一种颅骨三维模型构建方法,其特点是包括以下步骤:
(a)用激光三维扫描仪对一个标准颅骨进行三维扫描,分别扫描颅骨的正位和左右侧45°三张三维图像,将图像数据读入计算机,通过对图像进行处理、定标,合成为完整的颅骨三维数据模型,在该颅骨的三维数据模型上,依次拾取所需要的解剖结构特征点,得到标准模型的解剖结构特征点,完成颅骨标准模板的建立;
(b)在受试者面部鼻尖处粘贴一个铅点,分别采集受试者正位、左侧位和颅顶颏底位相互正交的三张X光头颅定位片图像,利用X线片图像上耳眼眶平面的四个点,对X光头颅定位片图像进行标定,分别拾取耳眼眶平面的四个点在三个方向上的投影,得到耳眼眶平面四个点的三维坐标,根据两点间的距离公式,计算出左侧位片中标尺投影上两个点图像空间的像素距离,确定颅骨的实际大小;
(c)在步骤(a)建立的颅骨标准模板上,拾取其耳眼眶平面的四个点,以步骤(b)采集的正交X光头颅定位片图像上耳眼眶平面四个点的三维坐标为基准点,计算颅骨标准模板耳眼眶平面四个点与基准点一一对应重合时的变换矩阵,并利用变换矩阵对颅骨标准模板三维数据进行变换,得到受试者的颅骨初始三维数据模型和初始特征点三维坐标;
(d)将受试者的初始特征点三维坐标分别在X、Y和Z方向投影,将三个方向的投影点分别映射到标定后的三张X光头颅定位片图像上,构建出受试者实际的特征点,根据公式
di=Pi-Qi
计算每个特征点的位移di,式中Pi为实际特征点,Qi为初始特征点;通过计算N个已知特征点的位移,构造一个插值函数di=f(Pi),i=1,...,N,采用径向基函数来构造插值函数,通过解方程组: k=1,2,3;以矩阵形式可表示为:λikAij=fk,其中Aij=φ(‖Pi-Pj‖),fk(Pi)=Pi-Qj,通过解线性方程组求得λik;在计算出λik后,利用插值函数,将初始特征点的每个顶点代入径向基函数中,求出其它非特征点的位移,计算出实际模型中每个网格顶点的新坐标,构建出颅骨三维数据模型。
本发明的有益效果是:由于采用三张正交X线片图像,通过调整特征点,得到了精确的颅骨实际的三维特征点;根据两组耳眼眶平面点对颅骨标准模板进行缩放来获得颅骨初始模型,使用初始模型特征点和实际特征点,采用基于径向基函数的形变算法进行颅骨模型的构建,得到了完整的颅骨三维模型。
下面结合附图和实施例对本发明作详细说明。
附图说明
附图1是本发明颅骨三维模型构建方法的流程图。
具体实施方式
参照附图1,本实施例针对某男性青年,年龄22岁,软组织直面型,左右对称,面部形态无异常且无正畸、正颌及外伤史进行颅部硬组织三维构建。具体步骤如下:
第一步,用激光三维扫描仪,采用高保质扫描方式,对一个标准颅骨进行三维扫描,分别扫描颅骨的正位和左右侧45°三张三维图像,扫描过程中保持颅骨水平端正,这样可以保证获得的颅骨图像是端正的。将这三张图像数据读入计算机,通过运行扫描仪自带软件对图像进行处理合成为一个完整的颅骨三维数据模型。
对获得的颅骨三维数据模型进行定标。在颅骨三维数据模型上拾取耳眼眶平面的四个点,包括左右耳点和左右眼眶点,计算左右耳点的中点坐标,由这个坐标值构造一个平移矩阵,对颅骨三维数据模型进行平移,再进行旋转调整,使得耳眼眶平面位于坐标轴XOZ平面上,且左右耳点在X轴上,且关于坐标原点对称,平面四个点的Y轴坐标近似为0,这样就完成了对颅骨三维数据模型的定标,得到颅骨的标准三维数据模型。
预先确定所需要的解剖结构特征点和特征点选择的顺序,在完成定标的颅骨的标准三维数据模型上,依次拾取这些特征点,得到标准模型的解剖结构特征点,完成颅骨标准模板的建立。
第二步,在受试者面部鼻尖处粘贴一个铅点,然后将受试者面部固定在头颅定位仪中进行X光摄影,分别采集正位、左侧位和颅顶颏底位三张X光头颅定位片,这三张图像相互正交。
对采集的X线片图像,进行基本图像增强处理。根据X线片图像的视觉效果,采用直方图均衡化、对比度和亮度调节的方法对图像进行增强。还可对图像进行反色处理,使图像的特征清晰、易于辨识,达到最佳的视觉效果。
利用X线片图像上耳眼眶平面的四个点,在坐标系中对每幅图像进行平移、旋转调整,使得正位片图像中的犁状孔投影位于Y轴、耳眼眶平面四个点的投影位于X轴,使得侧位片中左右耳点投影位于原点,左右眼眶点的投影位于Z轴,颅顶颏底位片图像中犁状孔投影位于Z轴。然后利用三张X线片的正交性质,确定铅点在颅顶颏底位片图像上投影坐标是惟一的,完成三张X线片图像的标定。
在完成标定的三张正交X线片图像上,分别拾取耳眼眶平面的四个点,根据这三幅X线片图像是正交的关系,而拾取的三组点为三维空间中耳眼眶平面四个点在三个方向上的投影,根据这个原理可直接得到耳眼眶平面四个点的三维坐标。
利用侧位X光片图像存在的标尺,拾取标尺上表示1cm的两个点P1、P2,设这两个点的图像坐标分别为X1、Y1和X2、Y2,根据两点间的距离公式 计算出这两个点图像空间的像素距离,即得到了实际空间中1cm代表的图像空间的像素距离,即确定了颅骨的实际大小。
第三步,导入颅骨标准模板的三维数据模型和该标准模板所对应的解剖结构的特征点。在这个颅骨三维数据模型上,拾取其耳眼眶平面的四个点,以正交X线片图像耳眼眶平面四个点的三维坐标为基准点,计算标准模板耳眼眶平面四个点与基准点一一对应重合时的变换矩阵。利用这个变换矩阵对标准模板三维数据进行变换,得到受试者的颅骨的初始三维数据模型;同理利用这个变换矩阵对标准模板的解剖结构特征点进行变换,得到受试者的初始特征点三维坐标。
第四步,将受试者的初始特征点三维坐标分别在X、Y和Z方向投影,将三个方向的投影点分别映射到标定后的三张X线片图像上。观察每幅图像上映射的特征点是否与图像上该特征点的位置重合,对于不重合的特征点进行调整,使得映射特征点与图像实际特征点的位置重合。基于正交三视图的原理,对调整过后的所有图像特征点构建出受试者实际的特征点。
根据这两组初始特征点和实际的特征点,计算每个特征点的位移di=Pi-Qi,其中Pi为实际的特征点,Qi为初始模板的特征点。通过计算N个已知特征点的位移,构造一个插值函数di=f(Pi),i=1,...,N,采用径向基函数来构造插值函数。给定N个特征点,设一点P的三维坐标为x、y、z,则径向基函数的插值函数形式可表示为: 式中fk为径向基函数,P为对应网格顶点的三维坐标,Pi对应于模型特征点的三维坐标,λik为特征点Pi的权值,φ(r)为基函数,选择高斯函数作为基函数,N为特征点的个数,‖P-Pi‖为计算P-Pi的欧氏范数。
由实际特征点和初始模型特征点的三维坐标已知,为了求得插值函数中特征点的权值系数λik,通过解方程组: i,j=1,...,N,k=1,2,3。以矩阵形式可表示为:λikAij=fk,其中Aij=φ(‖Pi-Pj‖),fk(Pi)=Pi-Qj,通过解线性方程组求得λik。
在计算出λik后,利用这个插值函数,将初始模型上的每个顶点代入径向基函数中,求出其它非特征点的位移,进而计算出实际模型中每个网格顶点的新坐标,构建出颅骨三维数据模型。
Claims (1)
1、一种颅骨三维模型构建方法,其特征在于包括以下步骤:
(a)用激光三维扫描仪对一个标准颅骨进行三维扫描,分别扫描颅骨的正位和左右侧45°三张三维图像,将图像数据读入计算机,通过对图像进行处理、定标,合成为完整的颅骨三维数据模型,在该颅骨的三维数据模型上,依次拾取所需要的解剖结构特征点,得到标准模型的解剖结构特征点,完成颅骨标准模板的建立;
(b)在受试者面部鼻尖处粘贴一个铅点,分别采集受试者正位、左侧位和颅顶颏底位相互正交的三张X光头颅定位片图像,利用X线片图像上耳眼眶平面的四个点,对X光头颅定位片图像进行标定,分别拾取耳眼眶平面的四个点在三个方向上的投影,得到耳眼眶平面四个点的三维坐标,根据两点间的距离公式,计算出左侧位片中标尺投影上两个点图像空间的像素距离,确定颅骨的实际大小;
(c)在步骤(a)建立的颅骨标准模板上,拾取其耳眼眶平面的四个点,以步骤(b)采集的正交X光头颅定位片图像上耳眼眶平面四个点的三维坐标为基准点,计算颅骨标准模板耳眼眶平面四个点与基准点一一对应重合时的变换矩阵,并利用变换矩阵对颅骨标准模板三维数据进行变换,得到受试者的颅骨初始三维数据模型和初始特征点三维坐标;
(d)将受试者的初始特征点三维坐标分别在X、Y和Z方向投影,将三个方向的投影点分别映射到标定后的三张X光头颅定位片图像上,构建出受试者实际的特征点,根据公式
di=Pi-Qi
计算每个特征点的位移di,式中Pi为实际特征点,Qi为初始特征点;通过计算N个已知特征点的位移,构造一个插值函数di=f(Pi),i=1,...,N,采用径向基函数来构造插值函数,通过解方程组: i,j=1,...,N,k=1,2,3;以矩阵形式可表示为:λikAij=fk,其中Aij=φ(‖Pi-Pj‖),fk(Pi)=Pi-Qj,通过解线性方程组求得λik;在计算出λik后,利用插值函数,将初始特征点的每个顶点代入径向基函数中,求出其它非特征点的位移,计算出实际模型中每个网格顶点的新坐标,构建出颅骨三维数据模型。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2008100172233A CN100538745C (zh) | 2008-01-04 | 2008-01-04 | 颅骨三维模型构建方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2008100172233A CN100538745C (zh) | 2008-01-04 | 2008-01-04 | 颅骨三维模型构建方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101216953A CN101216953A (zh) | 2008-07-09 |
CN100538745C true CN100538745C (zh) | 2009-09-09 |
Family
ID=39623381
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB2008100172233A Expired - Fee Related CN100538745C (zh) | 2008-01-04 | 2008-01-04 | 颅骨三维模型构建方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN100538745C (zh) |
Families Citing this family (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101339670B (zh) * | 2008-08-07 | 2010-06-09 | 浙江工业大学 | 一种计算机辅助的三维颅面复原方法 |
JP5463866B2 (ja) * | 2009-11-16 | 2014-04-09 | ソニー株式会社 | 画像処理装置および画像処理方法、並びにプログラム |
CN101969523B (zh) * | 2010-10-21 | 2012-10-03 | 西北农林科技大学 | 一种三维扫描装置及三维扫描方法 |
CN103978789B (zh) * | 2014-05-22 | 2016-05-11 | 中国科学院苏州生物医学工程技术研究所 | 基于3d打印的头部医学模型快速成型方法 |
CN104268930B (zh) * | 2014-09-10 | 2018-05-01 | 芜湖林一电子科技有限公司 | 一种坐标对比三维扫描方法 |
CN104524695B (zh) * | 2014-12-29 | 2016-11-02 | 中国医学科学院生物医学工程研究所 | 一种基于头部厚度测量建立标准头的缩放方法 |
US10762632B2 (en) * | 2018-09-12 | 2020-09-01 | Siemens Healthcare Gmbh | Analysis of skeletal trauma using deep learning |
CN109472786B (zh) * | 2018-11-05 | 2024-05-14 | 平安科技(深圳)有限公司 | 脑出血图像处理方法、装置、计算机设备及存储介质 |
CN110169782B (zh) * | 2019-06-27 | 2023-06-02 | 北京大学第三医院(北京大学第三临床医学院) | 一种颅面骨骼结构的头影测量方法 |
CN110570430B (zh) * | 2019-09-16 | 2023-01-03 | 中南大学 | 基于体配准的眼眶骨组织分割方法 |
CN112545537B (zh) * | 2019-09-26 | 2024-03-22 | 合肥美亚光电技术股份有限公司 | 头影测量描迹图生成方法及系统 |
CN111583219B (zh) * | 2020-04-30 | 2021-05-18 | 赤峰学院附属医院 | 颅颌面软硬组织的分析方法及装置、电子设备 |
CN112017278B (zh) * | 2020-08-28 | 2024-04-26 | 西安理工大学 | 一种基于Grasshopper的泳镜造型定制设计方法 |
CN112258492B (zh) * | 2020-10-30 | 2021-09-21 | 李艳 | 颅骨不对称信息的获取方法、存储介质及电子设备 |
CN112258639B (zh) * | 2020-10-30 | 2021-09-28 | 李艳 | 颅骨矫正信息的获取方法、存储介质及电子设备 |
CN112270743B (zh) * | 2020-10-30 | 2021-10-15 | 李艳 | 颅骨六视基准面的确定方法、装置、存储介质及电子设备 |
CN112288858B (zh) * | 2020-10-30 | 2022-11-11 | 李艳 | 人体骨骼的建模方法、存储介质及电子设备 |
CN112270705B (zh) * | 2020-10-30 | 2021-08-31 | 李艳 | 颅骨的切片方法、存储介质及电子设备 |
CN116650022B (zh) * | 2023-08-01 | 2023-11-24 | 湖南科迈森医疗科技有限公司 | 超声与内窥镜图像融合辅助定位子宫病灶的方法及系统 |
CN117274418B (zh) * | 2023-10-08 | 2024-04-02 | 北京长木谷医疗科技股份有限公司 | 基于正侧位x线图像的ct图像生成方法、装置及设备 |
-
2008
- 2008-01-04 CN CNB2008100172233A patent/CN100538745C/zh not_active Expired - Fee Related
Non-Patent Citations (6)
Title |
---|
X线片的颅颌骨特征点的三维确定. 李光耀等.计算机辅助设计与图形学学报,第18卷第1卷期. 2006 |
X线片的颅颌骨特征点的三维确定. 李光耀等.计算机辅助设计与图形学学报,第18卷第1卷期. 2006 * |
基于X光图像重建的颅面三维可视化系统. 潘俊君.第五届中国计算机图形学大会. 2004 |
基于X光图像重建的颅面三维可视化系统. 潘俊君.第五届中国计算机图形学大会. 2004 * |
基于颅骨的人脸建模技术研究及在法医面貌复原中的应用. 李康.中国优秀硕士学位论文全文数据库. 2006 |
基于颅骨的人脸建模技术研究及在法医面貌复原中的应用. 李康.中国优秀硕士学位论文全文数据库. 2006 * |
Also Published As
Publication number | Publication date |
---|---|
CN101216953A (zh) | 2008-07-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN100538745C (zh) | 颅骨三维模型构建方法 | |
CN103208133B (zh) | 一种图像中人脸胖瘦的调整方法 | |
US9330206B2 (en) | Producing a three dimensional model of an implant | |
US7936911B2 (en) | 3D planning and prediction method for optimizing facial skeleton symmetry in orthognathic surgery | |
CN105250062B (zh) | 一种基于医学影像的3d打印骨骼矫形支具制备方法 | |
Tzou et al. | Comparison of three-dimensional surface-imaging systems | |
Kähler et al. | Reanimating the dead: reconstruction of expressive faces from skull data | |
CN106295170B (zh) | 基于牙冠与牙根特征的牙齿建模方法 | |
KR100327541B1 (ko) | 3차원 얼굴 모델링 시스템 및 모델링 방법 | |
AU598466B2 (en) | System and method for the display of surface structures contained within the interior region of a solid body | |
CN102663818A (zh) | 颅颌面三维形貌模型的构建方法及其装置 | |
US20200273248A1 (en) | Method for manipulating 3d objects by flattened mesh | |
CN104299250A (zh) | 基于先验模型的正面人脸图像合成方法及系统 | |
CN103325143A (zh) | 基于模型匹配的标记点自动注册方法 | |
WO2017156905A1 (zh) | 一种将二维图像转化为多视点图像的显示方法及系统 | |
CN110459301A (zh) | 基于热力图和面部关键点的脑部神经外科导航配准方法 | |
CN110458932A (zh) | 图像处理方法、装置、系统、存储介质和图像扫描设备 | |
CN115619773B (zh) | 一种三维牙齿多模态数据配准方法及系统 | |
CN112489135B (zh) | 一种虚拟三维人脸面部重建系统的标定方法 | |
CN109118455B (zh) | 一种基于现代人软组织分布的古人类头骨颅面交互复原方法 | |
CN101702056B (zh) | 基于立体图像对的立体图像显示方法 | |
CN105389444A (zh) | 一种个性化种植牙的龈缘曲线设计方法 | |
CN107578465A (zh) | 一种模拟展示人体整形效果的方法及系统 | |
JP5391293B2 (ja) | 身体再建モデル形状データ作成装置及び方法並びにプログラム | |
CN103700130A (zh) | 发音器官的网格动画生成方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
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: 20090909 Termination date: 20120104 |