买专利卖专利找龙图腾，真高效！ 查专利查商标用IPTOP,全免费！专利年费监控用IP管家,真方便！

申请/专利权人：南京医科大学

摘要：本发明提供一种基于三维空间面型的正畸专用微笑美学评测方法及系统，通过建立以两台摄像机的中点为世界坐标原点的三维坐标系；受试者在自然头位下发出姿势性微笑，左摄像头L、右摄像头R同时对受试者的微笑进行采集，获得面部图像；标定左、右图像上的面部特征点，确定面部特征点的三维世界坐标；确定参考平面及参考线，根据面部特征点的三维坐标计算出用于微笑评测的评测指标。该种基于三维空间面型的正畸专用微笑美学评测方法及系统，针对口腔医学治疗前、后微笑美学评测，所得到的指标能够辅助口腔医师对患者微笑作出客观、定量的评价，具有实用范围广，适应性强，操作简单，灵活性强的特点。

主权项：1.一种基于三维空间面型的正畸专用微笑美学评测方法，其特征在于：包括以下步骤，S1、利用张正友标定法对左摄像头L、右摄像头R进行标定，得到摄像头的内参、外参及畸变参数，建立以两台摄像机的中点为世界坐标原点的三维坐标系；S2、受试者在自然头位下发出姿势性微笑，左摄像头L、右摄像头R同时对受试者的微笑进行采集，获得面部图像；S3、在步骤S2所获得面部图像上标定左、右图像上的面部特征点，确定面部特征点的三维世界坐标；S31、双目标定：根据步骤S1所获得的左、右摄像头的内参和畸变参数，结合标定角点对左右摄像头进行双目标定，得到左、右摄像头相机坐标系间的变换矩阵，包括旋转矢量和平移矢量，以及对应的本征矩阵和基础矩阵；S32、立体矫正：根据步骤S1所获得的左、右摄像头的内参、畸变参数及步骤S31获得的左、右摄像头坐标系间的变换矩阵，对左、右摄像头得到的图像进行立体矫正，使得左、右摄像头的图像位于同一平面且行对准，求得到相应的重投影矩阵；S33、标定面部特征点：标定立体矫正后的图像上的面部特征点，需标定的面部特征点包括：瞳孔点、鼻翼点、鼻下点、颏前点、口角点、上唇上缘中点、左唇峰点、右唇峰点、左侧上唇中点、右侧上唇中点、上口点、下口点、下唇下缘中点、左侧下唇中点和右侧下唇中点，具体为，瞳孔点：瞳孔中心点，分为左瞳孔点Nlp和右瞳孔点Nrp；鼻翼点：鼻部宽度方向最远点，分为左鼻翼点Nln和右鼻翼点Nrn；鼻下点：鼻小柱与上唇连接点Nun；颏前点：颏部最凸点Ncf；口角点：口角两侧最外点，分为左口角点Nlm和右口角点Nrm；上唇上缘中点：上唇上缘的中点Nuluc，将上唇上缘分为左右两个部分；左唇峰点：左侧唇峰点Nllp；右唇峰点：右侧唇峰点Nrlp；左侧上唇中点：上唇上缘左侧的中点Nlulc；右侧上唇中点：上唇上缘右侧的中点Nrulc；上口点：上唇下缘中点Numc；下口点：下唇上缘中点Nlmc；下唇下缘中点：下唇下缘的中点Nllc，将下唇下缘分为左右两个部分；左侧下唇中点：下唇下缘左侧的中点Nllcl；右侧下唇中点：下唇下缘右侧的中点Nllcr；左上中切牙龈缘点：左上中切牙龈缘最上点Nltu；右上中切牙龈缘点：右上中切牙龈缘最上点Nrtu；左上中切牙切缘点：左上中切牙切缘最下点Nltl；左上尖牙远中点：左上尖牙最远中点Nltc；右上尖牙远中点：右上尖牙最远中点Nrtc；左上尖牙牙尖点：左上尖牙牙尖顶点Nltt；右上尖牙牙尖点：右上尖牙牙尖顶点Nrtt；牙列中线：左、右中切牙间的平行牙长轴作为牙列中线l；牙列中线上端点：牙列中线上方端点Ntcu；牙列中线下端点：牙列中线下方端点Ntcl；S34、计算面部特征点的三维世界坐标：由步骤S33标记得到的面部特征点N在左、右图像上的二维坐标得到面部特征点在左、右相机图像间的视差，根据步骤S32所获得的重投影矩阵，计算面部特征点在世界坐标系中的三维坐标xN,yN,zN；S4、确定参考平面及参考线，包括正中矢状面、咬合平面和牙列中线，根据面部特征点的三维坐标计算出用于微笑评测的评测指标。

全文数据：基于三维空间面型的正畸专用微笑美学评测方法及系统技术领域本发明涉及一种基于三维空间面型的正畸专用微笑美学评测方法及系统。背景技术微笑是人类共同的语言，在人际交往和情感表达中具有重要作用，在口腔医学的临床实践中，微笑的美观与和谐是口腔医师和就诊患者所共同追求的目标之一。微笑已成为口腔医学临床检查诊断、疗效预测及最终评价不可忽视的参考资料。但是，针对口腔正畸治疗前后面部美学尤其是微笑美学的评定，至今主要通过医师和患者感受等定性判断，没有一个统一、标准、客观且量化的第三方测评系统。同时，能够针对基于三维空间的面型微笑进行评测的方法，国内外目前也鲜见报道。上述问题是在正畸治疗前后的微笑美学评测过程中应当予以考虑并解决的问题。发明内容基于当今时代的需求和现有技术存在的不足，本发明的目的是提供一种基于三维空间面型的正畸专用微笑美学评测方法及系统解决现有技术中存在的通过医师和患者感受等定性判断，没有一个统一、标准、客观且量化的问题。本发明的技术解决方案是：一种基于三维空间面型的正畸专用微笑美学评测方法包括以下步骤，S1、利用张正友标定法对左摄像头L、右摄像头R进行标定，得到摄像头的内参、外参及畸变参数，建立以两台摄像机的中点为世界坐标原点的三维坐标系；S2、受试者在自然头位下发出姿势性微笑，左摄像头L、右摄像头R同时对受试者的微笑进行采集，获得面部图像；S3、在步骤S2所获得面部图像上标定左、右图像上的面部特征点，确定面部特征点的三维世界坐标；S4、确定参考平面及参考线，包括正中矢状面、咬合平面和牙列中线，根据面部特征点的三维坐标计算出用于微笑评测的评测指标。进一步地，步骤S3包括：S31、双目标定：根据步骤S1所获得的左、右摄像头的内参和畸变参数，结合标定角点对左右摄像头进行双目标定，得到左、右摄像头相机坐标系间的变换矩阵，包括旋转矢量和平移矢量，以及对应的本征矩阵和基础矩阵；S32、立体矫正：根据步骤S1所获得的左、右摄像头的内参、畸变参数及步骤S31获得的左、右摄像头坐标系间的变换矩阵，对左、右摄像头得到的图像进行立体矫正，使得左、右摄像头的图像位于同一平面且行对准，求得到相应的重投影矩阵；S33、标定面部特征点：标定立体矫正后的图像上的面部特征点，需标定的面部特征点包括：瞳孔点、鼻翼点、鼻下点、颏前点、口角点、上唇上缘中点、左唇峰点、右唇峰点、左侧上唇中点、右侧上唇中点、上口点、下口点、下唇下缘中点、左侧下唇中点和右侧下唇中点，具体为，瞳孔点：瞳孔中心点，分为左瞳孔点Nlp和右瞳孔点Nrp；鼻翼点：鼻部宽度方向最远点，分为左鼻翼点Nln和右鼻翼点Nrn；鼻下点：鼻小柱与上唇连接点Nun；颏前点：颏部最凸点Ncf；口角点：口角两侧最外点，分为左口角点Nlm和右口角点Nrm；上唇上缘中点：上唇上缘的中点Nuluc，将上唇上缘分为左右两个部分；左唇峰点：左侧唇峰点Nllp；右唇峰点：右侧唇峰点Nrlp；左侧上唇中点：上唇上缘左侧的中点Nlulc；右侧上唇中点：上唇上缘右侧的中点Nrulc；上口点：上唇下缘中点Numc；下口点：下唇上缘中点Nlmc；下唇下缘中点：下唇下缘的中点Nllc，将下唇下缘分为左右两个部分；左侧下唇中点：下唇下缘左侧的中点Nllcl；右侧下唇中点：下唇下缘右侧的中点Nllcr；左上中切牙龈缘点：左上中切牙龈缘最上点Nltu；右上中切牙龈缘点：右上中切牙龈缘最上点Nrtu；左上中切牙切缘点：左上中切牙切缘最下点Nltl；左上尖牙远中点：左上尖牙最远中点Nltc；右上尖牙远中点：右上尖牙最远中点Nrtc；左上尖牙牙尖点：左上尖牙牙尖顶点Nltt；右上尖牙牙尖点：右上尖牙牙尖顶点Nrtt；牙列中线：左、右中切牙间的平行牙长轴作为牙列中线l；牙列中线上端点：牙列中线上方端点Ntcu；牙列中线下端点：牙列中线下方端点Ntcl；S34、计算面部特征点的三维世界坐标：由步骤S33标记得到的面部特征点N在左、右图像上的二维坐标得到面部特征点在左、右相机图像间的视差，根据步骤S32所获得的重投影矩阵，计算面部特征点在世界坐标系中的三维坐标xN,yN,zN。进一步地，步骤S4具体为，S41、确定头部正中矢状面S1：以左、右瞳孔连线作为头部水平基准线，过鼻下点作垂直于瞳孔连线的平面，将此平面作为头部正中矢状面，头部正中矢状面S1为：A1x+B1y+C1z+D1＝0，其中参数A1、B1、C1及D1的求解具体为，左瞳孔点Nlp的坐标为xlp,ylp,zlp，右瞳孔点Nrp的坐标为xrp,yrp,zrp，左右瞳孔连线的方向向量为xlp-xrp,ylp-yrp,zlp-zrp，头部正中矢状面的法向量为A1,B1,C1，由于左右瞳孔连线垂直于头部正中矢状面，则左右瞳孔连线的方向向量与头部正中矢状面的法向量平行，则：头部正中矢状面经过鼻下点Nun，鼻下点坐标为xun,yun,zun，将其带入平面方程，得到：A1xun+B1yun+C1zun+D1＝0联立求解以上两个方程，可求得参数A1、B1、C1及D1；S42、确定咬合平面S2：左上尖牙牙尖顶点、左上中切牙切缘点、右上尖牙牙尖顶点所成平面为咬合平面，咬合平面S2为：A2x+B2y+C2z+D2＝0，其中，参数A2、B2、C2及D2求解具体为，由于左上尖牙牙尖顶点Nltt、左上中切牙切缘点Nltl、右上尖牙牙尖顶点Nrtt经过咬合平面，将三点坐标代入头部正中矢状面方程，求得参数A2、B2、C2及D2；S43、确定牙列中线：左、右中切牙间的平行牙长轴作为牙列中线，分为上牙列中线和下牙列中线；S44、根据步骤S3所得的三维坐标，计算出用于评价微笑的评测指标。进一步地，步骤S44，计算出用于评价微笑的评测指标，具体为，S441、对于左右口角对称性评测指标，计算过程如下：左口角点Nlm的坐标为xlm，ylm，zlm，左口角点Nlm到矢状面的距离：右口角点Nrm的坐标为xrm，yrm，zrm，右口角点Nrm到矢状面的距离：得出左右口角对称性评测指标为|d1-d2|；S442、对于上唇对称性评测指标，计算过程如下：上唇上缘中点Nuluc的坐标为xuluc，yuluc，zuluc，其到矢状面的距离：左唇峰点Nllp的坐标为xllp,yllp,zllp到正中矢状面的距离为：右唇峰点Nrlp的坐标为xrlp,yrlp,zrlp，其到正中矢状面的距离为：最后，上唇对称性评测指标为d3+|d4-d5|；S443、对于下唇对称性评测指标，计算过程如下：下唇下缘中点Nllc的坐标为xllc，yllc，zllc，其到矢状面的距离为：左侧下唇中点Nllcl的坐标为xllcl，yllcl，zllcl，其到正中矢状面的距离为：右侧下唇中点Nllcr的坐标为xllcr，yllcr，zllcr，其到正中矢状面的距离为：最后，下唇对称性评测指标为d6+|d7-d8|；S444、对于牙列对称性评测指标，计算过程如下：牙列中线上端点的坐标为xtcu，ytcu，ztcu，其到矢状面的距离为：牙列中线下端点的坐标为xtcl，ytcl，ztcl，其到矢状面的距离为：最后，牙列对称性评测指标为：S445、对于鼻部对称性评测指标，计算过程如下：左鼻翼点Nln的坐标为xln，yln，zln，其到矢状面的距离：右鼻翼点Nrn的坐标为xrn，yrn，zrn，其到矢状面的距离：最后，鼻部对称性评测指标为|d11-d12|；S446、对于下巴颏对称性评测指标，计算过程如下：颏前点Ncf的坐标为xcf,ycf,zcf，其到矢状面的距离即为下巴颏对称性评测指标为：S447、对于口角宽度评测指标，计算过程如下：左口角点Nlm的坐标为xlm，ylm，zlm，右口角点Nrm的坐标为xrm，yrm，zrm，两点间的距离即为口角宽度评测指标为：S448、对于唇间隙评测指标，计算过程如下：上口点Numc的坐标为xumc，yumc，zumc，下口点Nlmc的坐标为xlmc，ylmc，zlmc，两点间的距离即为唇间隙评测指标为：S449、对于牙龈暴露量评测指标，计算过程如下：左侧上唇中点Nlulc的坐标为xlulc，ylulc，zlulc，左上中切牙龈缘点Nltu的坐标为xltu，yltu，zltu，两点间的距离为：右侧上唇中点Nrulc的坐标为xrulc，yrulc，zrulc，右上中切牙龈缘点Nrtu的坐标为xrtu，yrtu，zrtu，两点间的距离为：最后，牙龈暴露量评测指标为d16+d17；S4410、对于笑容指数评测指标：口角宽度为d14，唇间隙为d15，最后笑容指数评测指标为d14d15；S4411、对于鼻唇沟高度评测指标，计算过程如下：鼻下点Nun的坐标为xun，yun，zun，上唇上缘中点Nuluc的坐标为xuluc，yuluc，zuluc，两点间的距离即为鼻唇沟高度评测指标为：S4412、对于颏唇沟高度评测指标，计算过程如下：颏前点Ncf的坐标为xcf,ycf,zcf，下唇下缘中点Nllc的坐标为xllc,yllc,zllc，两点间的距离即为颏唇沟高度评测指标为：S4413、对于颊廊比率评测指标，计算过程如下：左上尖牙最远中点Nltc的坐标为xltc，yltc，zltc，右上尖牙最远中点Nrtc的坐标为xrtc，yrtc，zrtc，两点间的距离为：口角宽度为d14，该距离与口角宽度的比值即为颊廊比率评测指标为d20d14；S4414、对于咬合平面的倾斜度评测指标，计算过程如下：咬合平面S2为A2x+B2y+C2z+D2＝0，其法向量m为A2,B2,C2；左瞳孔点Nlp的坐标为xlp,ylp,zlp，右瞳孔点Nrp的坐标为xrp,yrp,zrp，所以瞳孔连线的方向向量n为xlp-xrp,ylp-yrp,zlp-zrp；所以咬合平面和瞳孔连线的角度即为咬合平面的倾斜度评测指标为：进一步地，步骤S2具体为，S21、受试者坐于距左、右摄像头中点水平位置，双眼平视前方，调整采集器支架的高度使受试者头部处于两相机公共视场的区域内；S22、操作者辅助受试者摆出自然头位并引导受试者自然微笑，同时左摄像头L、右摄像机R同步采集微笑照片，得到受试者在自然头位下微笑的面部图像。一种基于三维空间面型的正畸专用微笑美学评测系统，包括面部图像采集模块、数据分析处理模块和显示模块；面部图像采集模块：通过左摄像头L、右摄像头R同时对受试者进行采集，获得面部图像，并将所得面部图像发送给数据分析处理模块；数据分析处理模块：获得面部图像采集模块发送的面部图像，对面部图像采用上述任一项所述基于三维空间面型的正畸专用微笑美学评测方法的进行分析，确定面部特征点的三维坐标，进而得出评测指标；显示模块：用于显示面部图像采集模块传来的面部图像信息，使操作者对其进行面部特征点的人工标定，并显示数据分析处理模块所得评测指标。进一步地，左摄像头L、右摄像头R分别设于图像采集支架上，图像采集支架包括横杆、云台、高度调节套管和三脚架底座，横杆的两端分别设有云台，左摄像头L、右摄像头R分别设有云台上，高度调节套管的一端连接横杆的中点，高度调节套管的另一端连接三脚架底座上。进一步地，高度调节套管包括中轴、竖管和曲臂构成，中轴采用套管结构，中轴的端部设有通孔，通孔内设有竖管，中轴的侧部设有曲臂，曲臂的端部螺纹穿入中轴并接触竖管，竖管标有刻度，竖杆的端部连接横杆的中点。本发明的有益效果是，该种基于三维空间面型的正畸专用微笑美学评测方法及系统，与现有技术相比，具有以下优点：一、该种基于三维空间面型的正畸专用微笑美学评测方法及系统，针对口腔医学治疗前、后微笑美学评测，具有实用范围广，适应性强，操作简单，灵活性强的特点。二、本发明首次提出基于三维空间的微笑评测系统及方法，能够推进微笑面像更好地应用于口腔正畸临床检查诊断、疗效预测及治疗效果评价。三、本发明设计了一个统一、标准且客观的第三方测评系统及方法，实现口腔医学治疗后面部美学尤其是微笑美学的定量分析及评测，所得到的指标能够辅助口腔医师对患者微笑作出客观、定量的评价。附图说明图1是本发明实施例基于三维空间面型的正畸专用微笑美学评测方法的流程示意图；图2是实施例中面部特征点及参考平面和参考线的说明示意图；图3是实施例中采集器支架的示意图；其中：1-云台，2-横杆，3-竖管，4-中轴，5-曲臂，6-三脚架底座。具体实施方式下面结合附图详细说明本发明的优选实施例。实施例一种基于三维空间面型的正畸专用微笑美学评测方法，如图1，包括以下步骤，S1.利用张正友标定法对左摄像头L、右摄像头R进行标定，得到摄像头的内参、外参及畸变参数，建立以两台摄像机的中点为世界坐标原点的三维坐标系。其中，左摄像头L为受试者右手侧摄像头，右摄像头R为受试者左手侧摄像头。建立以两台摄像机的中点为世界坐标原点的三维坐标系，具体为：取左摄像头镜头L的中心点O1与右摄像头镜头R的中心点Or连线中点O为坐标系原点，两台摄像机镜头中心连线为X轴，与地面平行，Y轴垂直于地面，Z轴垂直于X轴与Y轴。S2.受试者在自然头位下发出姿势性微笑，左、右摄像机同时对此进行采集，获得面部图像。步骤S2具体为，S21.受试者坐于距左、右摄像头中点水平距离1.5m处，双眼平视前方，调整采集器支架的高度使受试者头部处于两相机公共视场的区域内。S22.操作者辅助受试者摆出自然头位并引导受试者自然微笑，同时左、右摄像机同步采集微笑照片，得到受试者在自然头位下微笑的面部图像。S3.由操作者人工标定左、右图像上的面部特征点，确定面部特征点的三维世界坐标。具体为，特征点N在摄像头L与摄像头R所采集的左、右图像上的对应点为Nl、Nr，坐标为根据双目立体视觉的原理，由面部特征点的二维坐标和两摄像头的内参、外参，确定面部特征点的三维世界坐标xN,yN,zN。步骤S3具体为，S31.双目标定：根据步骤S1所获得的左、右摄像头的内参和畸变参数，结合标定角点对左右摄像头进行双目标定，得到左、右摄像头相机坐标系间的变换矩阵，包括旋转矢量和平移矢量，以及对应的本征矩阵和基础矩阵。S32.立体矫正：根据步骤S1所获得的左、右摄像头的内参、畸变参数及步骤S31获得的左、右摄像头坐标系间的变换矩阵，对左、右摄像头得到的图像进行立体矫正，使得左、右摄像头的图像位于同一平面且行对准，求得到相应的重投影矩阵。S33.标定面部特征点：操作者手工标定立体矫正后的图像上的面部特征点，需标定的面部特征点包括：瞳孔点、鼻翼点、鼻下点、颏前点、口角点、上唇上缘中点、左唇峰点、右唇峰点、左侧上唇中点、右侧上唇中点、上口点、下口点、下唇下缘中点、左侧下唇中点和右侧下唇中点，面部特征点位置及序号如图2所示，需标定的面部特征点具体包括：瞳孔点：瞳孔中心点，分为左瞳孔点Nlp和右瞳孔点Nrp；鼻翼点：鼻部宽度方向最远点，分为左鼻翼点Nln和右鼻翼点Nrn；鼻下点：鼻小柱与上唇连接点Nun；颏前点：颏部最凸点Ncf；口角点：口角两侧最外点，分为左口角点Nlm和右口角点Nrm；上唇上缘中点：上唇上缘的中点Nuluc，将上唇上缘分为左右两个部分；左唇峰点：左侧唇峰点Nllp；右唇峰点：右侧唇峰点Nrlp；左侧上唇中点：上唇上缘左侧的中点Nlulc；右侧上唇中点：上唇上缘右侧的中点Nrulc；上口点：上唇下缘中点Numc；下口点：下唇上缘中点Nlmc；下唇下缘中点：下唇下缘的中点Nllc，将下唇下缘分为左右两个部分；左侧下唇中点：下唇下缘左侧的中点Nllcl；右侧下唇中点：下唇下缘右侧的中点Nllcr；左上中切牙龈缘点：左上中切牙龈缘最上点Nltu；右上中切牙龈缘点：右上中切牙龈缘最上点Nrtu；左上中切牙切缘点：左上中切牙切缘最下点Nltl；左上尖牙远中点：左上尖牙最远中点Nltc；右上尖牙远中点：右上尖牙最远中点Nrtc；左上尖牙牙尖点：左上尖牙牙尖顶点Nltt；右上尖牙牙尖点：右上尖牙牙尖顶点Nrtt；牙列中线：左、右中切牙间的平行牙长轴作为牙列中线l；牙列中线上端点：牙列中线上方端点Ntcu；牙列中线下端点：牙列中线下方端点Ntcl。S34.计算特征点三维世界坐标：，由步骤S33标记得到的特征点N在左、右图像上的二维坐标得到特征点在左、右相机图像间的视差，根据步骤S32所获得的重投影矩阵，计算特征点在世界坐标系中的三维坐标xN,yN,zN。S4.确定参考平面及参考线，包括正中矢状面、咬合平面和牙列中线，根据面部特征点的三维坐标计算出用于微笑评测的各项评测指标。步骤S4具体为，S41.确定头部正中矢状面S1：以左、右瞳孔连线作为头部水平基准线，过鼻下点作垂直于瞳孔连线的平面，将此平面作为头部正中矢状面。步骤S41中，定义头部正中矢状面S1为：A1x+B1y+C1z+D1＝0，其中参数A1、B1、C1及C1的求解具体为，左瞳孔点Nlp的坐标为xlp,ylp,zlp，右瞳孔点Nrp的坐标为xrp,yrp,zrp，左右瞳孔连线的方向向量为xlp-xrp,ylp-yrp,zlp-zrp，头部正中矢状面的法向量为A1,B1,C1，由于左右瞳孔连线垂直于头部正中矢状面，则左右瞳孔连线的方向向量与头部正中矢状面的法向量平行，则：头部正中矢状面经过鼻下点Nun，鼻下点坐标为xun，yun，zun，将其带入平面方程，得到：A1xun+B1yun+C1zun+D1＝0联立求解以上两个方程，可求得参数A1、B1、C1及D1。S42.确定咬合平面S2：左上尖牙牙尖顶点、左上中切牙切缘点、右上尖牙牙尖顶点所成平面为咬合平面。步骤S42中，定义咬合平面S2为：A2x+B2y+C2z+D2＝0，其中，参数A2、B2、C2及D2求解具体为，由于左上尖牙牙尖顶点Nltt、左上中切牙切缘点Nltl、右上尖牙牙尖顶点Nrtt经过咬合平面，将三点坐标代入头部正中矢状面方程，求得参数A2、B2、C2及D2。S43.确定牙列中线：左、右中切牙间的平行牙长轴作为牙列中线，分为上牙列中线和下牙列中线。S44.根据步骤S3所得的三维坐标，计算出用于评价微笑的评测指标。步骤S44中，计算出用于评价微笑的评测指标，具体为，S441、对于左右口角对称性评测指标，计算过程如下：左口角点Nlm的坐标为xlm，ylm，zlm，左口角点Nlm到矢状面的距离：右口角点Nrm的坐标为xrm，yrm，zrm，右口角点Nrm到矢状面的距离：得出左右口角对称性评测指标为|d1-d2|；S442、对于上唇对称性评测指标，计算过程如下：上唇上缘中点Nuluc的坐标为xuluc，yuluc，zuluc，其到矢状面的距离：左唇峰点Nllp的坐标为xllp,yllp,zllp到正中矢状面的距离为：右唇峰点Nrlp的坐标为xrlp,yrlp,zrlp，其到正中矢状面的距离为：最后，上唇对称性评测指标为d3+|d4-d5|；S443、对于下唇对称性评测指标，计算过程如下：下唇下缘中点Nllc的坐标为xllc，yllc，zllc，其到矢状面的距离为：左侧下唇中点Nllcl的坐标为xllcl，yllcl，zllcl，其到正中矢状面的距离为：右侧下唇中点Nllcr的坐标为xllcr，yllcr，zllcr，其到正中矢状面的距离为：最后，下唇对称性评测指标为d6+|d7-d8|；S444、对于牙列对称性评测指标，计算过程如下：牙列中线上端点的坐标为xtcu，ytcu，ztcu，其到矢状面的距离为：牙列中线下端点的坐标为xtcl，ytcl，ztcl，其到矢状面的距离为：最后，牙列对称性评测指标为：S445、对于鼻部对称性评测指标，计算过程如下：左鼻翼点Nln的坐标为xln，yln，zln，其到矢状面的距离：右鼻翼点Nrn的坐标为xrn，yrn，zrn，其到矢状面的距离：最后，鼻部对称性评测指标为|d11-d12|；S446、对于下巴颏对称性评测指标，计算过程如下：颏前点Ncf的坐标为xcf,ycf,zcf，其到矢状面的距离即为下巴颏对称性评测指标为：S447、对于口角宽度评测指标，计算过程如下：左口角点Nlm的坐标为xlm，ylm，zlm，右口角点Nrm的坐标为xrm，yrm，zrm，两点间的距离即为口角宽度评测指标为：S448、对于唇间隙评测指标，计算过程如下：上口点Numc的坐标为xumc，yumc，zumc，下口点Nlmc的坐标为xlmc，ylmc，zlmc，两点间的距离即为唇间隙评测指标为：S449、对于牙龈暴露量评测指标，计算过程如下：左侧上唇中点Nlulc的坐标为xlulc，ylulc，zlulc，左上中切牙龈缘点Nltu的坐标为xltu，yltu，zltu，两点间的距离为：右侧上唇中点Nrulc的坐标为xrulc，yrulc，zrulc，右上中切牙龈缘点Nrtu的坐标为xrtu，yrtu，zrtu，两点间的距离为：最后，牙龈暴露量评测指标为d16+d17；S4410、对于笑容指数评测指标：口角宽度为d14，唇间隙为d15，最后笑容指数评测指标为d14d15；S4411、对于鼻唇沟高度评测指标，计算过程如下：鼻下点Nun的坐标为xun,yun,zun，上唇上缘中点Nuluc的坐标为xuluc，yuluc，zuluc，两点间的距离即为鼻唇沟高度评测指标为：S4412、对于颏唇沟高度评测指标，计算过程如下：颏前点Ncf的坐标为xcf,ycf,zcf，下唇下缘中点Nllc的坐标为xllc,yllc,zllc，两点间的距离即为颏唇沟高度评测指标为：S4413、对于颊廊比率评测指标，计算过程如下：左上尖牙最远中点Nltc的坐标为xltc，yltc，zltc，右上尖牙最远中点Nrtc的坐标为xrtc，yrtc，zrtc，两点间的距离为：口角宽度为d14，该距离与口角宽度的比值即为颊廊比率评测指标为d20d14；S4414、对于咬合平面的倾斜度评测指标，计算过程如下：咬合平面S2为A2x+B2y+C2z+D2＝0，其法向量m为A2,B2,C2；左瞳孔点Nlp的坐标为xlp,ylp,zlp，右瞳孔点Nrp的坐标为xrp,yrp,zrp，所以瞳孔连线的方向向量n为xlp-xrp,ylp-yrp,zlp-zrp；所以咬合平面和瞳孔连线的角度即为咬合平面的倾斜度评测指标为：实施例的步骤S44中，评测指标及其定义如下表1所示。表1实施例中评测指标及其定义实施例的基于三维空间面型的正畸专用微笑美学评测方法，能够进行统一、标准且客观的测评，实现口腔医学治疗后面部美学尤其是微笑美学的定量分析及评测，所得到的指标能够辅助口腔医师对患者微笑作出客观、定量的评价。实施例还提供一种基于三维空间面型的正畸专用微笑美学评测系统，包括面部图像采集模块、数据分析处理模块和显示模块。面部图像采集模块：通过左摄像头L、右摄像头R同时对受试者进行采集，获得面部图像，并将所得面部图像发送给数据分析处理模块。数据分析处理模块：获得面部图像采集模块发送的面部图像，对面部图像采用上述任一项所述基于三维空间面型的正畸专用微笑美学评测方法的进行分析，确定面部特征点的三维坐标，进而得出评测指标。显示模块：用于显示面部图像采集模块传来的面部图像信息，使操作者对其进行面部特征点的人工标定，并显示数据分析处理模块所得评测指标。实施例中，如图3，左摄像头L、右摄像头R分别设于图像采集支架上，图像采集支架包括横杆2、云台1、高度调节套管和三脚架底座6，横杆2的两端分别设有云台1，左摄像头L、右摄像头R分别设有云台1上，高度调节套管的一端连接横杆2的中点，高度调节套管的另一端连接三脚架底座6上。高度调节套管包括中轴4、竖管3和曲臂5构成，中轴4采用套管结构，中轴4的端部设有通孔，通孔内设有竖管3，中轴4的侧部设有曲臂5，曲臂5的端部螺纹穿入中轴4并接触竖管3，竖管3标有刻度，竖杆的端部连接横杆2的中点。实施例中，摄像头的分辨率为2048X1536，拍摄范围为0.8-1.4m。分辨率和拍摄范围可以根据实际需求及情况进行调整，由此满足人体真实面部所需的精确率需求，适应头颅尺寸的个体差异。竖管3上标有的刻度，可以标定当前左摄像头L、右摄像头R所在的高度。图像采集支架为光学图像采集器提供机械支撑，并调节拍摄的范围。该种基于三维空间面型的正畸专用微笑美学评测方法及系统，针对口腔医学治疗前、后微笑美学评测，具有实用范围广，适应性强，操作简单，灵活性强的特点。实施例首次提出基于三维空间的微笑评测方法及系统，能够推进微笑面像更好地应用于口腔正畸临床检查诊断、疗效预测及治疗效果评价。

权利要求：1.一种基于三维空间面型的正畸专用微笑美学评测方法，其特征在于：包括以下步骤，S1、利用张正友标定法对左摄像头L、右摄像头R进行标定，得到摄像头的内参、外参及畸变参数，建立以两台摄像机的中点为世界坐标原点的三维坐标系；S2、受试者在自然头位下发出姿势性微笑，左摄像头L、右摄像头R同时对受试者的微笑进行采集，获得面部图像；S3、在步骤S2所获得面部图像上标定左、右图像上的面部特征点，确定面部特征点的三维世界坐标；S4、确定参考平面及参考线，包括正中矢状面、咬合平面和牙列中线，根据面部特征点的三维坐标计算出用于微笑评测的评测指标。2.如权利要求1所述的基于三维空间面型的正畸专用微笑美学评测方法，其特征在于：步骤S3包括：S31、双目标定：根据步骤S1所获得的左、右摄像头的内参和畸变参数，结合标定角点对左右摄像头进行双目标定，得到左、右摄像头相机坐标系间的变换矩阵，包括旋转矢量和平移矢量，以及对应的本征矩阵和基础矩阵；S32、立体矫正：根据步骤S1所获得的左、右摄像头的内参、畸变参数及步骤S31获得的左、右摄像头坐标系间的变换矩阵，对左、右摄像头得到的图像进行立体矫正，使得左、右摄像头的图像位于同一平面且行对准，求得到相应的重投影矩阵；S33、标定面部特征点：标定立体矫正后的图像上的面部特征点，需标定的面部特征点包括：瞳孔点、鼻翼点、鼻下点、颏前点、口角点、上唇上缘中点、左唇峰点、右唇峰点、左侧上唇中点、右侧上唇中点、上口点、下口点、下唇下缘中点、左侧下唇中点和右侧下唇中点，具体为，瞳孔点：瞳孔中心点，分为左瞳孔点Nlp和右瞳孔点Nrp；鼻翼点：鼻部宽度方向最远点，分为左鼻翼点Nln和右鼻翼点Nrn；鼻下点：鼻小柱与上唇连接点Nun；颏前点：颏部最凸点Ncf；口角点：口角两侧最外点，分为左口角点Nlm和右口角点Nrm；上唇上缘中点：上唇上缘的中点Nuluc，将上唇上缘分为左右两个部分；左唇峰点：左侧唇峰点Nllp；右唇峰点：右侧唇峰点Nrlp；左侧上唇中点：上唇上缘左侧的中点Nlulc；右侧上唇中点：上唇上缘右侧的中点Nrulc；上口点：上唇下缘中点Numc；下口点：下唇上缘中点Nlmc；下唇下缘中点：下唇下缘的中点Nllc，将下唇下缘分为左右两个部分；左侧下唇中点：下唇下缘左侧的中点Nllcl；右侧下唇中点：下唇下缘右侧的中点Nllcr；左上中切牙龈缘点：左上中切牙龈缘最上点Nltu；右上中切牙龈缘点：右上中切牙龈缘最上点Nrtu；左上中切牙切缘点：左上中切牙切缘最下点Nltl；左上尖牙远中点：左上尖牙最远中点Nltc；右上尖牙远中点：右上尖牙最远中点Nrtc；左上尖牙牙尖点：左上尖牙牙尖顶点Nltt；右上尖牙牙尖点：右上尖牙牙尖顶点Nrtt；牙列中线：左、右中切牙间的平行牙长轴作为牙列中线l；牙列中线上端点：牙列中线上方端点Ntcu；牙列中线下端点：牙列中线下方端点Ntcl；S34、计算面部特征点的三维世界坐标：由步骤S33标记得到的面部特征点N在左、右图像上的二维坐标得到面部特征点在左、右相机图像间的视差，根据步骤S32所获得的重投影矩阵，计算面部特征点在世界坐标系中的三维坐标xN,yN,zN。3.如权利要求2所述的基于三维空间面型的正畸专用微笑美学评测方法，其特征在于：步骤S4具体为，S41、确定头部正中矢状面S1：以左、右瞳孔连线作为头部水平基准线，过鼻下点作垂直于瞳孔连线的平面，将此平面作为头部正中矢状面，头部正中矢状面S1为：A1x+B1y+C1z+D1＝0，其中参数A1、B1、C1及D1的求解具体为，左瞳孔点Nlp的坐标为xlp,ylp,zlp，右瞳孔点Nrp的坐标为xrp,yrp,zrp，左右瞳孔连线的方向向量为xlp-xrp,ylp-yrp,zlp-zrp，头部正中矢状面的法向量为A1,B1,C1，由于左右瞳孔连线垂直于头部正中矢状面，则左右瞳孔连线的方向向量与头部正中矢状面的法向量平行，则：头部正中矢状面经过鼻下点Nun，鼻下点坐标为xun,yun,zun，将其带入平面方程，得到：A1xun+B1yun+C1zun+D1＝0联立求解以上两个方程，可求得参数A1、B1、C1及D1；S42、确定咬合平面S2：左上尖牙牙尖顶点、左上中切牙切缘点、右上尖牙牙尖顶点所成平面为咬合平面，咬合平面S2为：A2x+B2y+C2z+D2＝0，其中，参数A2、B2、C2及D2求解具体为，由于左上尖牙牙尖顶点Nltt、左上中切牙切缘点Nltl、右上尖牙牙尖顶点Nrtt经过咬合平面，将三点坐标代入头部正中矢状面方程，求得参数A2、B2、C2及D2；S43、确定牙列中线：左、右中切牙间的平行牙长轴作为牙列中线，分为上牙列中线和下牙列中线；S44、根据步骤S3所得的三维坐标，计算出用于评价微笑的评测指标。4.如权利要求3所述的基于三维空间面型的正畸专用微笑美学评测方法，其特征在于：步骤S44，计算出用于评价微笑的评测指标，具体为，S441、对于左右口角对称性评测指标，计算过程如下：左口角点Nlm的坐标为xlm,ylm,zlm，左口角点Nlm到矢状面的距离：右口角点Nrm的坐标为xrm,yrm,zrm，右口角点Nrm到矢状面的距离：得出左右口角对称性评测指标为|d1-d2|；S442、对于上唇对称性评测指标，计算过程如下：上唇上缘中点Nuluc的坐标为xuluc,yuluc,zuluc，其到矢状面的距离：左唇峰点Nllp的坐标为xllp,yllp,zllp到正中矢状面的距离为：右唇峰点Nrlp的坐标为xrlp,yrlp,zrlp，其到正中矢状面的距离为：最后，上唇对称性评测指标为d3+|d4-d5|；S443、对于下唇对称性评测指标，计算过程如下：下唇下缘中点Nllc的坐标为xllc,yllc,zllc，其到矢状面的距离为：左侧下唇中点Nllcl的坐标为xllcl,yllcl,zllcl，其到正中矢状面的距离为：右侧下唇中点Nllcr的坐标为xllcr,yllcr,zllcr，其到正中矢状面的距离为：最后，下唇对称性评测指标为d6+|d7-d8|；S444、对于牙列对称性评测指标，计算过程如下：牙列中线上端点的坐标为xtcu,ytcu,ztcu，其到矢状面的距离为：牙列中线下端点的坐标为xtcl,ytcl,ztcl，其到矢状面的距离为：最后，牙列对称性评测指标为：S445、对于鼻部对称性评测指标，计算过程如下：左鼻翼点Nln的坐标为xln,yln,zln，其到矢状面的距离：右鼻翼点Nrn的坐标为xrn,yrn,zrn，其到矢状面的距离：最后，鼻部对称性评测指标为|d11-d12|；S446、对于下巴颏对称性评测指标，计算过程如下：颏前点Ncf的坐标为xcf,ycf,zcf，其到矢状面的距离即为下巴颏对称性评测指标为：S447、对于口角宽度评测指标，计算过程如下：左口角点Nlm的坐标为xlm,ylm,zlm，右口角点Nrm的坐标为xrm,yrm,zrm，两点间的距离即为口角宽度评测指标为：S448、对于唇间隙评测指标，计算过程如下：上口点Numc的坐标为xumc,yumc,zumc，下口点Nlmc的坐标为xlmc,ylmc,zlmc，两点间的距离即为唇间隙评测指标为：S449、对于牙龈暴露量评测指标，计算过程如下：左侧上唇中点Nlulc的坐标为xlulc,ylulc,zlulc，左上中切牙龈缘点Nltu的坐标为xltu,yltu,zltu，两点间的距离为：右侧上唇中点Nrulc的坐标为xrulc,yrulc,zrulc，右上中切牙龈缘点Nrtu的坐标为xrtu,yrtu,zrtu，两点间的距离为：最后，牙龈暴露量评测指标为d16+d17；S4410、对于笑容指数评测指标：口角宽度为d14，唇间隙为d15，最后笑容指数评测指标为d14d15；S4411、对于鼻唇沟高度评测指标，计算过程如下：鼻下点Nun的坐标为xun,yun,zun，上唇上缘中点Nuluc的坐标为xuluc,yuluc,zuluc，两点间的距离即为鼻唇沟高度评测指标为：S4412、对于颏唇沟高度评测指标，计算过程如下：颏前点Ncf的坐标为xcf,ycf,zcf，下唇下缘中点Nllc的坐标为xllc,yllc,zllc，两点间的距离即为颏唇沟高度评测指标为：S4413、对于颊廊比率评测指标，计算过程如下：左上尖牙最远中点Nltc的坐标为xltc,yltc,zltc，右上尖牙最远中点Nrtc的坐标为xrtc,yrtc,zrtc，两点间的距离为：口角宽度为d14，该距离与口角宽度的比值即为颊廊比率评测指标为d20d14；S4414、对于咬合平面的倾斜度评测指标，计算过程如下：咬合平面S2为A2x+B2y+C2z+D2＝0，其法向量m为A2,B2,C2；左瞳孔点Nlp的坐标为xlp,ylp,zlp，右瞳孔点Nrp的坐标为xrp,yrp,zrp，所以瞳孔连线的方向向量n为xlp-xrp,ylp-yrp,zlp-zrp；所以咬合平面和瞳孔连线的角度即为咬合平面的倾斜度评测指标为：5.如权利要求1所述的基于三维空间面型的正畸专用微笑美学评测方法，其特征在于：步骤S2具体为，S21、受试者坐于距左、右摄像头中点水平位置，双眼平视前方，调整采集器支架的高度使受试者头部处于两相机公共视场的区域内；S22、操作者辅助受试者摆出自然头位并引导受试者自然微笑，同时左摄像头L、右摄像机R同步采集微笑照片，得到受试者在自然头位下微笑的面部图像。6.一种基于三维空间面型的正畸专用微笑美学评测系统，其特征在于：包括面部图像采集模块、数据分析处理模块和显示模块；面部图像采集模块：通过左摄像头L、右摄像头R同时对受试者进行采集，获得面部图像，并将所得面部图像发送给数据分析处理模块；数据分析处理模块：获得面部图像采集模块发送的面部图像，对面部图像采用权利要求1-5任一项所述基于三维空间面型的正畸专用微笑美学评测方法的进行分析，确定面部特征点的三维坐标，进而得出评测指标；显示模块：用于显示面部图像采集模块传来的面部图像信息，使操作者对其进行面部特征点的人工标定，并显示数据分析处理模块所得评测指标。7.如权利要求6所述的基于三维空间面型的正畸专用微笑美学评测系统，其特征在于：左摄像头L、右摄像头R分别设于图像采集支架上，图像采集支架包括横杆、云台、高度调节套管和三脚架底座，横杆的两端分别设有云台，左摄像头L、右摄像头R分别设有云台上，高度调节套管的一端连接横杆的中点，高度调节套管的另一端连接三脚架底座上。8.如权利要求7所述的基于三维空间面型的正畸专用微笑美学评测系统，其特征在于：高度调节套管包括中轴、竖管和曲臂构成，中轴采用套管结构，中轴的端部设有通孔，通孔内设有竖管，中轴的侧部设有曲臂，曲臂的端部螺纹穿入中轴并接触竖管，竖管标有刻度，竖杆的端部连接横杆的中点。

百度查询： 南京医科大学 基于三维空间面型的正畸专用微笑美学评测方法及系统