CN112700487B - 一种头颅侧位片中测量标尺刻度获取方法及系统 - Google Patents

Abstract

本发明提供了一种头颅侧位片中测量标尺检测方法及系统，包括步骤：在头颅侧位片中根据头颅标志点确定头颅测量标尺的图片；根据第一预设算法对头颅测量标尺的图片进行数字化处理，确定头颅测量标尺在所述头颅侧位片中所处的相对位置信息；根据第二预设算法对头颅测量标尺图进行图像处理，获取在头颅侧位片中头颅测量标尺的刻度间隙和标尺的刻度相位；根据头颅测量标尺的刻度间隙和标尺的刻度相位，并结合人体上颌长度均值计算头颅测量标尺的刻度。本发明解决了现有技术中头颅侧位片的标尺距离需要人工参与，距离数据存在误差，效率比较低的问题。

Description

一种头颅侧位片中测量标尺刻度获取方法及系统

技术领域

本发明属于牙齿矫治技术领域，更确切的说涉及一种头颅侧位片中测量标尺刻度获取方法及系统、电子设备、及计算机存储介质。

背景技术

头颅侧位片是医生用于判断患者是牙性凸还是骨性凸的判据之一，使用标准的测量软件测量相关的数据，与标准值进行比对，使得牙颌面畸形的诊断、治疗设计更加准确。在计算相关测量值的过程中，除了需要识别头颅侧位片的标志点外，还需要检测出头颅侧位片上的标尺，以便准确，定量给出测量数据。

现阶段，大部分头颅侧位片的标尺需要人工点出标尺上的两个点，并给出两点在标尺上所表示的距离，才能判定图像中像素之间对应的实际距离。这种判定方式在头颅侧位片自动化分析中有着很大的短板，需要人工参与，效率比较低。

基于以上本申请提供了解决以上技术问题的技术方案。

发明内容

本发明的主要目的是克服现有技术中存在的缺陷，提供“一种头颅侧位片中测量标尺刻度获取系统、电子设备、及计算机存储介质”，解决了现有技术中头颅侧位片的标尺距离需要人工参与，数据存在误差，效率比较低的问题。

本发明提供的技术方案如下：

一种头颅侧位片中测量标尺的刻度的获取方法，包括：

在头颅侧位片中根据头颅标志点确定头颅测量标尺的图片；根据第一预设算法对头颅测量标尺的图片进行数字化处理，确定头颅测量标尺在所述头颅侧位片中所处的相对位置信息；根据第二预设算法对头颅测量标尺图进行图像处理，获取在头颅侧位片中头颅测量标尺的刻度间隙和标尺的刻度相位；根据头颅测量标尺的刻度间隙和标尺的刻度相位，并结合人体上颌长度均值计算头颅测量标尺的刻度值。

进一步优选的，确定头颅测量标尺在所述头颅侧位片中所处的相对位置信息包括：将头颅测量标尺的图片根据预置分割算法分割成若干子图块；

将每一子图块进行傅里叶变换处理，并将傅里叶变换处理后每一子图块的变换信息中同频率的子图块进行标记整理，并提取出头颅测量标尺在所述头颅侧位片中所处的相对位置信息。

进一步优选的，将分割成若干子图块包括：

L_max＝L_K＝min(N₁,N₂)/N_B

L_min＝L₁＝min(h,L_max/K),

L_n＝(L_max/L_min)^(n-1)/(K-1)L_min)；N₁,N₂分别代表头颅测量标尺的图片的长和宽的长度(也即为像素的长度)；N_B表示长度最小的边上预设分割的子图块的数目，L_n第n可取的子图块块的大小，通过对数尺寸得到；h为子图块最小块尺寸；K是指需要给出的L值的个数，子图块的长度L_n-1值小于等于L_n；图像处理中有意义的最小的部分可以假设为h×h的平方块。

进一步优选的，确定头颅测量标尺在所述头颅侧位片中所处的相对位置信息包括：将头颅测量标尺的图片设定在预设遍历区域内，获取预设遍历区域在头颅侧位片中所处的位姿信息；所述位姿信息包括：头颅测量标尺在的头颅片中的倾斜角度和相对位置信息。

进一步优选的，包括：在预设角度范围内生成与每一角度值对应遍历线，

在生成的所述遍历线上获取头颅测量标尺图的像素值；将获取头颅测量标尺图的像素值进行傅里叶变换得到得到不同角度下，不同的位置的傅里叶变换图；在生成的傅里叶变换图中获取幅值最高点对应的遍历线信息，并在遍历线信息提取头颅测量标尺图的所述位姿信息。

进一步优选的，还包括：通过Hough变换算法在所述头颅侧位片中获取直线位置信息，将获取的直线位置信息设置头颅测量标尺在所述头颅侧位片中所处的相对位置信息。

进一步优选的，获取标尺的刻度间隙和标尺的刻度相位包括：

利用功率密度函数获取头颅测量标尺在的头颅片中的倾斜角度；

根据所述倾斜角度获取与标尺刻度线垂直的直线，并提取直线的像素值；

将提取的直线上像素的像素值通过累积平均归一化平方差函数获取头颅测量标尺的相对位置信息的基频，通过基频获取头颅测量标尺的刻度间隙值；

对头颅测量标尺图进行最小二乘拟合得到标尺的正弦模型，通过正弦模型可以得到标尺的刻度相位。

进一步优选的，测量标尺的倾斜角度包括：

功率谱密度最大的点的坐标为(v,u)。

进一步优选的，每个子图块的功率密度包括：PSDF(ω)＝F{Y_B}F^*{Y_B}，F为Y_B的频谱。

进一步优选的，获取所述基频值包括：

自相关ACF函数：为函数值y(j)与自身平移T后的函数值y(j+T)之间的相似性；T为y平移数。

进一步优选的，包括：通过最小二乘拟合得到刻度标尺的正弦模型，

(Y_B(x)-M(x,α))²，

M--正弦模型F--为的Y_B频谱，2π/λ为在尺子方向上每个像素周期；λ为标尺刻度之间的像素数。α₁,α₂分别表示sin函数和cos函数的幅值，α＝{α₁,α₂}通过迭代求出的。

一种头颅侧位片中测量标尺刻度获取系统，可执行上述任一项所述头颅测量中的测量标尺刻度的获取方法，包括，测量标尺图像获取模块，在头颅侧位片中根据头颅标志点确定头颅测量标尺的图片；

标尺图像位置获取模块，根据第一预设算法通过所述测量标尺图像获取模块中获取的头颅测量标尺的图片进行数字化处理，确定头颅测量标尺在所述头颅侧位片中所处的相对位置信息；

标尺图像处理模块，根据第二预设算法对标尺图像位置获取模块中获取的头颅测量标尺的位置信息进行图像处理，获取在头颅侧位片中头颅测量标尺的刻度间隙和标尺的刻度相位；

标尺刻度获取模块，根据标尺图像处理模块中获取的头颅测量标尺的刻度间隙和标尺的刻度相位，并结合人体上颌长度均值计算头颅测量标尺的刻度。

一种电子设备，包括处理器和存储器，所述处理器执行所述存储器存储的计算机指令，使得所述电子设备执行权利上述任一项所述的头颅侧位片中测量标尺刻度获取方法。

一种计算机存储介质，包括计算机指令，当所述计算机指令在电子设备上运行时，使得所述电子设备执行上述任一项所述的头颅侧位片中测量标尺刻度获取方法。

本发明提供头颅侧位片中测量标尺刻度获取方法及系统、电子设备、计算机存储介质，能够带来以下至少一种有益效果：

在本申请中在头颅侧位片中粗定位出标尺的位置，标尺精确位置的确定，进一步迭代算出标尺的刻度间隙，以及标尺的刻度相位，通过计算机语言实现，以数据形式展示，解决了现有技术的人工参与效率低下的问题。

在本申请中在标尺刻度求取，是通过标准数据库中以不同换牙阶段人体的上颌长的平均值不同，并以该平均作为参考通过傅里叶变换，Hough变换等多种方式完成标尺刻度的求取，实现了自动化求取，解放了手工参与，使其数据更加精细化，进一步的为隐形正畸提供了精确的数据依据，同时提高了工作效率；特别针对疑难病例给出精确的数据信息，给有需求的患者更加完美的体验。

附图说明

下面将以明确易懂的方式，结合附图说明优选实施方式，对上述特性、技术特征、优点及其实现方式予以进一步说明。

图1为本发明头颅侧位片中测量标尺刻度获取方法实施例的流程图；

图2为本发明头颅侧位片中测量标尺刻度图；

图3为本发明头颅侧位片中测量标尺刻度获取方法实施例另一流程图；

图4为本发明头颅侧位片中测量标尺刻度获取方法实施例另一流程图；

图5为本发明头颅侧位片上颌长的位置展示图；

图6为本发明头颅侧位片中测量标尺刻度图图像处理后的分割子图块；

图7为本发明头颅侧位片中测量标尺刻度图的变换数据图谱；

图8为本发明头颅侧位片中测量标尺刻度图图像处理后标尺图；

图9为本发明头颅侧位片中测量标尺刻度获取系统实施例的结构图；

图10为本发明头颅侧位片中测量标尺图算法遍历图；

图11为本发明电子设备的结构图。

具体实施方式

为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对照附图说明本发明的具体实施方式。显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图，并获得其他的实施方式。

本申请是是求取头颅侧位片标志点之间的距离信息，通过距离信息与标准值进行比对，通过确切的数值信息判断矫治患者的口腔状态，方便医生制定患者的矫治方案，具体的通过求取侧位片中标尺刻度，根据像素值以及标尺距离进一步的求取相应的标尺刻度距离，刻度距离用来分析像素/刻度值用来确定中侧位片中两点之间的距离；现有技术中都是手工操作，或者凭借医生的经验完成测量，不满足现有大环境的数字化时代，无法满足自动诊疗。基于此参见以下的实施例，本申请提供了解决技术问题的技术方案。

参见图1至图11所示，本发明提供了一种头颅侧位片中测量标尺检测方法的实施例，包括步骤：

步骤S100在头颅侧位片中根据头颅标志点确定头颅测量标尺的图片；

根据头颅侧位片拍摄的规律，标尺在图中的位置参见图2所示的右上角的位置，在进行拍摄头颅侧位片时标尺的放置位置一般为鼻根点(N)的前方，结合图5所示，即面向右侧的侧位片其标尺的位置应为在鼻根点右侧，面向左侧的侧位片的标尺位置为鼻根点左侧，以鼻根点为参考，根据鼻根点的坐标截出只带有标尺的图片，可以以世界坐标系为基础，则以头颅侧位片为基础获取相应的标尺的图片。

在本申请中是以鼻根点作为颅部标志点之一，可以通过不限于随机森林的方法自动查找得到，在头颅中有多个标志点。本申请中针对于标志点的选取不做限制。

步骤S200根据第一预设算法对头颅测量标尺的图片进行数字化处理，确定头颅测量标尺在所述头颅侧位片中所处的相对位置信息；

参见图6所示，针对获取到的头颅测量标尺进行数字化处理，也即转化为数字信息，并进行定位其位置关系，根据位置关系求取相关的标尺信息。

本申请中针对第一预设算法获取头颅侧位片中头颅测量标尺的相对位置信息提供以下的实施例方式参见以下：

参见图3和图6所示，第一种实施例方式为：

步骤S211将头颅测量标尺的图片根据预置分割算法分割成若干子图块；

步骤S212将每一子图块进行傅里叶变换处理；

步骤S213并将傅里叶变换处理后的子图块中同频率的子图块进行标记整理；

步骤S214提取出头颅测量标尺在所述头颅侧位片中所处的相对位置信息。

具体的，为了准确定位标尺的位置，可以将参见图6所示，将获得的头颅测量标尺的图片分成不同的块，块的大小为L，对每个小块做傅里叶变-换，具有相同频率特征的小块进行标记，小块区域就是需要精准定位的标尺位置。如下图6所示，为某一粗定位图分割成一小块后傅里叶计算结果图，矩形块圈出了具有相同傅里叶特征的小块，则同频率小块的进行整理后，则为头颅测量标尺的相对位置。

具体的根据预置分割算法分割成若干子图块的实施方式为：

不同的块的大小L可以估计出不同的标尺间隔距离d(L)。对于输入的图像，{L₁,L₂,...L_K}继续自动选择包含部分标尺图像的块，K是指需要给出的L值的个数。最大的方块边L(L_max)和最小的边(L_min)都设置为输入图像最小边的分数倍；

L_max＝L_K＝min(N₁,N₂)/N_B

L_min＝L₁＝min(h,L_max/K),

L_n＝(L_max/L_min)⁽ⁿ-1)/(K-1)L_min)；

N₁,N₂分别代表头颅测量标尺的图片的长和宽的长度(也即为像素的长度)；N_B表示长度最小的边上预设分割的子图块的数目，L_n第n可取的子图块块的大小，通过对数尺寸得到；h为子图块最小块尺寸；K是指需要给出的L值的个数，子图块的长度L_n-1值小于等于L_n；图像处理中有意义的最小的部分可以假设为h×h的平方块。

从最小的L₁开始，假设第i个L，即，L_i那么将图像分成(N₁/L_i)*(N₂/L_i)块，对每一个小方块进行信号处理，直到发现至少个B_min连续的小方块有相同的频率特征，B_min为连续的具有相同频率方块的最少个数，则停止选择，即，选择的块长为L_i。如果第i个不满足条件特征，那么选取第i+1个L。通过数值实验发现，B_min是最小的，足以消除大多数的二元关系(单纯地由于噪声信号而选出来的候选点)的数目。自然，更大的将会减少伪正数据(背景环境)。可是，同时，如果B_min太大，将会同样减少伪负数目。

第二种实施例方式为：

步骤S221将头颅测量标尺的图片设定在预设遍历区域内；

步骤S222获取预设遍历区域在头颅侧位片中所处的位姿信息；所述位姿信息包括：头颅测量标尺在的头颅片中的倾斜角度和位置信息(也即为相对于位置信息)。

参见图10所示，遍历过程为：在预设角度范围内生成与每一角度值对应遍历线，在生成的所述遍历线上获取头颅测量标尺图的像素值；将获取头颅测量标尺图的像素值进行傅里叶变换得到不同幅值的傅里叶变换标尺图；在生成的不同幅值的傅里叶变换标尺图中获取幅值最高点对应的遍历线信息，并在遍历线信息提取头颅测量标尺图的所述位姿信息。

同样是利用标尺刻度存在周期性的规律，采用遍历的方法寻找标尺的位置以及角度。图10中θ为遍历直线与竖直方向的夹角，根据经验规律θ的取值范围假设为为[-10°,10°]，每一根遍历线上取对应的像素值，做一次傅里叶变换，于是每个θ下都会生成一个傅里叶变换图如下图7所示，图7中幅值最高的点对应的遍历线即可为标尺刻度所在的线。并在标尺刻度所在的线的函数中提取出刻度标尺的相对位置信息。

第三种实施例方式为：

结参见图8所示，通过Hough变换算法在所述头颅侧位片中获取直线的位置信息，将获取的直线位置信息设置头颅测量标尺在所述头颅侧位片中所处的相对位置信息。

步骤S300根据第二预设算法对头颅测量标尺图进行图像处理，获取在头颅侧位片中头颅测量标尺的刻度间隙和标尺的刻度相位；

具体的，在精确定位标尺相对位置信息，可以利用自相关函数受噪声影响比较低的特性，计算累积平均归一化平方差函数(CSDF Cumulative Mean NormalisedDifference Functions)算出基频，基频算出后，通过最小二乘拟合得到标尺的正弦模型。通过基频可以得到标尺的刻度间隙。

通过最小二乘拟合得到标尺的正弦模型包括：

(Y_B(x)-M(x；α))²，--(1)，Y_B(x)为切割出来的标尺栅格函数；

M--正弦模型F--为YB的频谱，2π/λ为在尺子方向上每个像素周期；λ为图像像素。是需要迭代求出的，α₁α₂分别表示sin函数和cos函数的幅值，α＝{α₁,α₂}。

步骤S310利用功率密度函数获取头颅测量标尺在的头颅片中的倾斜角度；在本申请的第一种位置获取方式中，可以将标尺划分为一个个小方块，对小方块求功率密度(自相关函数的频谱)，在功率谱密度图中求出功率谱密度最大的点的坐标，则标尺的方向：

测量标尺的倾斜角度包括：功率谱密度最大的点的坐标为(v,u)；

每个子图块的功率密度：PSDF(ω)＝F{Y_B}F^*{Y_B}--(4)

步骤S320根据所述倾斜角度方向获取与标尺刻度线垂直的直线，并提取直线的像素值；直线方程包括：(x-a)/(b-y)＝tanθ--(5)假设确定标尺的相对位置为(a,b)，(如果是截出的标尺的图，坐标为截图的中心点位置)，标尺的倾斜角为θ。

步骤S330将提取的直线上像素的像素值通过累积平均归一化平方差函数获取头颅测量标尺的相对位置信息的基频，通过基频获取头颅测量标尺的刻度间隙值；获取所述基频值包括：

自相关ACF函数：为函数值y(j)与自身平移T后的函数值y(j+T)之间的相似性；T为y平移数；本申请中自相关函数ACF：反应的是同一序列在不同时刻的取值之间的相关程度。y只获取的一系列的数值，要对这些数值做自相关处理，而T值就是相当于将这个数列延迟了T以后的数值。每一个j下，得到一个ACF值，进而得到一个CSDF值。W是一个窗，一般做信号处理的时候都要做一个加窗处理，这里可以认为是截取的数据的长度。j为第j个数据点。

步骤S340对头颅测量标尺进行最小二乘拟合得到标尺的正弦模型，通过正弦模型可以得到标尺的刻度相位。

(Y_B(x)-M(x,α))²，

M--正弦模型F--为Y_B的频谱，为在尺子方向上每个像素周期；λ为标尺刻度之间的像素数；α₁,α₂分别表示sin函数和cos函数的幅值，α＝{α₁,α₂}通过迭代求出的。

步骤S400根据头颅测量标尺的刻度间隙和标尺的刻度相位结合人体上颌长度均值获取头颅测量标尺的刻度相位，计算头颅测量标尺的刻度。

S＝[(L_ANP/g)*λ]--(7)

S--头颅测量标尺的刻度值,L_ANP--人体上颌长度均值，g人体上颌长度均值长度对应的像素数，λ--头颅测量标尺的刻度间隙之间的像素数，[]为取整。

参见图5所示，上颌长(ANS-Ptm)：翼上颌裂点与前鼻棘点与在FH平面上垂足之间距离，代表上颌的长度。

参见图9所示，本申请还提供了一种头颅侧位片中测量标尺刻度获取系统的实施例，可执行上述任一项所述头颅测量中的测量标尺刻度的获取方法的实施例，包括：测量标尺图像获取模块100，在头颅侧位片中根据头颅标志点确定头颅测量标尺的图片；标尺图像位置获取模块200，根据第一预设算法通过所述测量标尺图像获取模块中获取的头颅测量标尺的图片进行数字化处理，确定头颅测量标尺在所述头颅侧位片中所处的相对位置信息；标尺图像处理模块300，根据第二预设算法对标尺图像位置获取模块中获取的头颅测量标尺的位置信息进行图像处理，获取在头颅侧位片中头颅测量标尺的刻度间隙和标尺的刻度相位；标尺刻度获取模块400，根据标尺图像处理模块中获取的头颅测量标尺的刻度间隙和标尺的刻度相位，并结合人体上颌长度均值计算头颅测量标尺的刻度。

具体的，本申请的一种头颅侧位片中测量标尺刻度获取系统的实施例可执行上述一种头颅侧位片中测量标尺刻度获取方法实施例，在此不再赘述。

在本申请中在头颅侧位片中粗定位出标尺的位置，标尺精确位置的确定，进一步迭代算出标尺的刻度间隙，以及标尺的刻度相位，通过计算机语言实现，以数据形式展示，解决了现有技术的人工参与效率低下的问题。

在本申请中在标尺刻度求取，是通过标准数据库中以不同换牙阶段人体的上颌长的平均值不同，并以该平均作为参考通过傅里叶变换，Hough变换等多种方式完成标尺刻度的求取，实现了自动化求取，解放了手工参与，使其数据更加精细化，进一步的为隐形正畸提供了精确的数据依据，同时提高了工作效率；特别针对疑难病例给出精确的数据信息，给有需求的患者更加完美的体验。

本实施例提供一种电子设备，其结构框图如图11所示，该电子设备1000可以是平板电脑、笔记本电脑或台式电脑。电子设备1000还可能被称为便携式终端、膝上型终端、台式终端等其他名称。

电子设备1000内置有处理器1001和存储器1002，其中，存储1002上存储有计算机程序，处理器1001运行存储器1002中的计算机程序时实现实施例一种头颅侧位片中测量标尺检测方法。

处理器1001可以包括一个或多个处理核心，比如4核心处理器、8核心处理器等。处理器1001可以采用DSP(Digital Signal Processing，数字信号处理)、FPGA(Field—Programmable Gate Array，现场可编程门阵列)、PLA(Programmable Logic Array，可编程逻辑阵列)中的至少一种硬件形式来实现。处理器1001也可以包括主处理器和协处理器，主处理器是用于对在唤醒状态下的数据进行处理的处理器，也称CPU(Central ProcessingUnit，中央处理器)；协处理器是用于对在待机状态下的数据进行处理的低功耗处理器。

在一些实施例中，处理器1001可以在集成有GPU(Graphics Processing Unit，图像处理器)，GPU用于负责显示屏所需要显示的内容的渲染和绘制。一些实施例中，处理器1001还可以包括AI(Artificial Intelligence，人工智能)处理器，该AI处理器用于处理有关机器学习的计算操作。

存储器1002可以包括一个或多个计算机可读存储介质，该计算机可读存储介质可以是非暂态的。存储器1002还可包括高速随机存取存储器，以及非易失性存储器，比如一个或多个磁盘存储设备、闪存存储设备。在一些实施例中，存储器1002中的非暂态的计算机可读存储介质用于存储至少一个指令、至少一段程序、代码集或指令集，该至少一条指令、至少一段程序、代码集或指令集用于被处理器1001所执行以实现本申请中实施例一种头颅侧位片中测量标尺检测方法。

在一些实施例中，电子设备1000还包括有：外围设备接口1003和外围设备。处理器1001、存储器1002和外围设备接口1003之间可以通过总线或信号线相连。外围设备可以通过总线、信号线或电路板与外围设备接口1003相连。

具体到本实施例中，外围设备可以包括口内扫描仪1004和3D打印设备1005。处理器1001通过口内扫描仪1004获取患者口内的数字化牙颌模型，处理器1001在执行计算机程序的过程中通过程序命令获取口内扫描仪1004采集的数字化牙颌模型，再通过执行实施例一种头颅侧位片中测量标尺检测方法，以获取牙龈形变参数，然后再根据获取的牙龈形变参数设计壳状牙齿矫治器，将设计好的数字化壳状牙齿矫治器模型对应的数据信息传输至3D打印设备1005，通过3D打印设备1005直接打印制备壳状牙齿矫治器。

本实施例还提供一种计算机可读存储介质，该计算机可读存储介质可以为非易失性计算机可读存储介质，该计算机可读存储介质也可以为易失性计算机可读存储介质。该计算机可读存储介质中存储有指令，当该指令在计算机上运行时，使得计算机执行实施例一种头颅侧位片中测量标尺检测方法。

实施例中的模块如果以软件功能模块的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的全部或部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(Read-only memory，ROM)、随机存取存储器(Random access memory，RAM)、磁碟或者光盘等各种可以存储程序代码的介质。

应当说明的是，上述实施例均可根据需要自由组合。以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。

Claims (13)

1.一种头颅侧位片中测量标尺的刻度的获取方法，包括：

在头颅侧位片中根据头颅标志点确定头颅测量标尺的图片；

根据第一预设算法对头颅测量标尺的图片进行数字化处理，确定头颅测量标尺在所述头颅侧位片中所处的相对位置信息；

根据第二预设算法对头颅测量标尺图进行图像处理，获取在头颅侧位片中头颅测量标尺的刻度间隙和标尺的刻度相位；

根据头颅测量标尺的刻度间隙和标尺的刻度相位，并结合人体上颌长度均值计算头颅测量标尺的刻度值，

所述获取在头颅侧位片中头颅测量标尺的刻度间隙和标尺的刻度相位，包括，

利用功率密度函数获取头颅测量标尺的头颅片中的倾斜角度，利用所述功率密度函数计算出功率谱密度最大的点的坐标，基于功率谱密度最大的点的坐标获得所述倾斜角度；

根据所述倾斜角度获取与标尺刻度线垂直的直线，并提取直线上所有像素的像素值；

将提取的直线上像素的像素值通过累积平均归一化平方差函数获取头颅测量标尺的相对位置信息的基频，通过基频获取头颅测量标尺的刻度间隙值；

对头颅测量标尺图进行最小二乘拟合得到标尺的正弦模型，通过正弦模型得到的标尺刻度之间的像素数从而得到标尺的刻度相位。

2.根据权利要求1头颅侧位片中测量标尺刻度的获取方法，确定头颅测量标尺在所述头颅侧位片中所处的相对位置信息包括：

将头颅测量标尺的图片根据预置分割算法分割成若干子图块；

将每一子图块进行傅里叶变换处理，并将傅里叶变换处理后每一子图块的变换信息中同频率的子图块进行标记整理，并提取出头颅测量标尺在所述头颅侧位片中所处的相对位置信息。

3.根据权利要求2头颅侧位片中测量标尺刻度的获取方法，将分割成若干子图块包括：

L_max＝L_K＝min(N₁,N₂)/N_B

L_min＝L₁＝min(h,L_max/K)

L_n＝(L_max/L_min)^(n-1)/(K-1)L_min

N₁,N₂分别代表头颅测量标尺的图片的长和宽的长度；N_B表示长度最小的边上预设分割的子图块的数目，L_n为第n可取的子图块的大小；h为子图块最小块尺寸；K是指需要给出的L值的个数，子图块的长度L_n-1值小于等于L_n；图像处理中有意义的最小的部分为h×h的块。

4.根据权利要求1头颅侧位片中测量标尺刻度的获取方法，确定头颅测量标尺在所述头颅侧位片中所处的相对位置信息包括：

将头颅测量标尺的图片设定在预设遍历区域内，获取预设遍历区域在头颅侧位片中所处的位姿信息；所述位姿信息包括：头颅测量标尺在的头颅片中的倾斜角度和相对位置信息。

5.根据权利要求4头颅侧位片中测量标尺刻度的获取方法，包括：在预设角度范围内生成与每一角度值对应遍历线，

在生成的所述遍历线上获取头颅测量标尺图的像素值；

将获取头颅测量标尺图的像素值进行傅里叶变换得到不同角度下，不同的位置的傅里叶变换图；

在生成的傅里叶变换图中获取幅值最高点对应的遍历线信息，并在遍历线信息提取头颅测量标尺图的所述位姿信息。

6.根据权利要求1头颅侧位片中测量标尺刻度的获取方法，还包括：通过Hough变换算法在所述头颅侧位片中获取直线的位置信息，将获取的直线位置信息设置头颅测量标尺在所述头颅侧位片中所处的相对位置信息。

7.根据权利要求1头颅侧位片中测量标尺刻度的获取方法，测量标尺的倾斜角度包括：

功率谱密度最大的点的坐标为(v,u)。

8.根据权利要求7头颅侧位片中测量标尺刻度的获取方法，每个子图块的功率密度包括：PSDF(w)＝F{Y_B}F^*{Y_B}，F为Y_B的频谱。

9.根据权利要求8头颅侧位片中测量标尺刻度的获取方法，获取所述基频值包括：

其中，相关ACF函数：为函数值y(j)与自身平移T后的函数值y(j+T)之间的相似性；T为y平移数。

10.根据权利要求9头颅侧位片中测量标尺刻度的获取方法，包括：通过最小二乘拟合(Y_B(x)-M(x,α))²得到刻度标尺的正弦模型，

M(x,y；θ_R,λ)＝α₁cos[F(x,y；θ_R,λ)]+α₂sin[F(x,y；θ_R,λ)]

其中，M表示正弦模型F表示为Y_B的频谱，2π/λ为在尺子方向上每个像素周期；λ为标尺刻度之间的像素数；α₁，α₂分别表示sin函数和cos函数的幅值，α＝{α₁，α₂}通过迭代求出的。

11.一种头颅侧位片中测量标尺刻度获取系统，可执行权利要求1-10任一项所述头颅测量中的测量标尺刻度的获取方法，包括，

测量标尺图像获取模块，在头颅侧位片中根据头颅标志点确定头颅测量标尺的图片；

标尺图像位置获取模块，根据第一预设算法通过所述测量标尺图像获取模块中获取的头颅测量标尺的图片进行数字化处理，确定头颅测量标尺在所述头颅侧位片中所处的相对位置信息；

标尺图像处理模块，根据第二预设算法对标尺图像位置获取模块中获取的头颅测量标尺的位置信息进行图像处理，获取在头颅侧位片中头颅测量标尺的刻度间隙和标尺的刻度相位；

标尺刻度获取模块，根据标尺图像处理模块中获取的头颅测量标尺的刻度间隙和标尺的刻度相位，并结合人体上颌长度均值计算头颅测量标尺的刻度。

12.一种电子设备，包括处理器和存储器，所述处理器执行所述存储器存储的计算机指令，使得所述电子设备执行权利要求1至10任一项所述的头颅侧位片中测量标尺刻度获取方法。

13.一种计算机存储介质，包括计算机指令，当所述计算机指令在电子设备上运行时，使得所述电子设备执行权利要求1至10任一项所述的头颅侧位片中测量标尺刻度获取方法。