CN108830915A - 口腔影像3d模拟动画实现方法及装置 - Google Patents
口腔影像3d模拟动画实现方法及装置 Download PDFInfo
- Publication number
- CN108830915A CN108830915A CN201810520990.XA CN201810520990A CN108830915A CN 108830915 A CN108830915 A CN 108830915A CN 201810520990 A CN201810520990 A CN 201810520990A CN 108830915 A CN108830915 A CN 108830915A
- Authority
- CN
- China
- Prior art keywords
- image
- oral cavity
- layer
- reconstruction
- obtains
- 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
Links
- 210000000214 mouth Anatomy 0.000 title claims abstract description 105
- 238000000034 method Methods 0.000 title claims abstract description 41
- 238000004088 simulation Methods 0.000 title claims abstract description 36
- 230000037182 bone density Effects 0.000 claims abstract description 46
- 238000003384 imaging method Methods 0.000 claims abstract description 29
- 238000004458 analytical method Methods 0.000 claims abstract description 26
- 230000010287 polarization Effects 0.000 claims abstract description 26
- 238000006243 chemical reaction Methods 0.000 claims abstract description 8
- 238000000354 decomposition reaction Methods 0.000 claims description 56
- 238000001914 filtration Methods 0.000 claims description 24
- 230000001815 facial effect Effects 0.000 claims description 20
- 238000006073 displacement reaction Methods 0.000 claims description 12
- 238000009877 rendering Methods 0.000 claims description 11
- 230000015556 catabolic process Effects 0.000 claims description 9
- 230000004927 fusion Effects 0.000 claims description 9
- 238000005070 sampling Methods 0.000 claims description 9
- 230000002708 enhancing effect Effects 0.000 claims description 7
- 210000000887 face Anatomy 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 claims description 3
- 210000003128 head Anatomy 0.000 claims description 3
- 108010048734 sclerotin Proteins 0.000 claims description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 3
- 238000005516 engineering process Methods 0.000 abstract description 5
- 238000001739 density measurement Methods 0.000 abstract description 2
- 238000003745 diagnosis Methods 0.000 abstract description 2
- 230000000694 effects Effects 0.000 description 5
- 230000006870 function Effects 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 210000001747 pupil Anatomy 0.000 description 4
- 210000000988 bone and bone Anatomy 0.000 description 2
- 235000013399 edible fruits Nutrition 0.000 description 2
- 230000001965 increasing effect Effects 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 238000010191 image analysis Methods 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 210000003739 neck Anatomy 0.000 description 1
- 230000001225 therapeutic effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T13/00—Animation
- G06T13/20—3D [Three Dimensional] animation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T15/00—3D [Three Dimensional] image rendering
- G06T15/005—General purpose rendering architectures
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T19/00—Manipulating 3D models or images for computer graphics
- G06T19/20—Editing of 3D images, e.g. changing shapes or colours, aligning objects or positioning parts
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T3/00—Geometric image transformations in the plane of the image
- G06T3/40—Scaling of whole images or parts thereof, e.g. expanding or contracting
- G06T3/4007—Scaling of whole images or parts thereof, e.g. expanding or contracting based on interpolation, e.g. bilinear interpolation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/20—Image enhancement or restoration using local operators
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/20—Analysis of motion
- G06T7/246—Analysis of motion using feature-based methods, e.g. the tracking of corners or segments
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2200/00—Indexing scheme for image data processing or generation, in general
- G06T2200/08—Indexing scheme for image data processing or generation, in general involving all processing steps from image acquisition to 3D model generation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10081—Computed x-ray tomography [CT]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20024—Filtering details
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30008—Bone
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30036—Dental; Teeth
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2210/00—Indexing scheme for image generation or computer graphics
- G06T2210/41—Medical
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2211/00—Image generation
- G06T2211/40—Computed tomography
- G06T2211/416—Exact reconstruction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2219/00—Indexing scheme for manipulating 3D models or images for computer graphics
- G06T2219/20—Indexing scheme for editing of 3D models
- G06T2219/2012—Colour editing, changing, or manipulating; Use of colour codes
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Graphics (AREA)
- Software Systems (AREA)
- Geometry (AREA)
- Architecture (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Multimedia (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Image Generation (AREA)
- Processing Or Creating Images (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本发明公开了一种口腔影像3D模拟动画实现方法及装置，实现方法包括：S1、获取待模拟口腔的不同角度的多个口腔影像信息，口腔影像为二维图像；S2、对不同角度的多个口腔影像信息进行三维容积重建，得到三维动画，三维动画信息包括三维动画影像集及其影像序列信息；S3、对口腔影像进行灰度值分析，得到影像的骨密度结果；S4、根据S2得到的三维动画和S3得到的灰度值分析结果，对三维动画影像集中的影像进行色彩转换处理，得到彩色图像集；S5、对彩色图像集中的每一个图像进行偏振成像处理。本发明基于数字化牙片的骨密度测量技术及计算机辅助图像密度分析技术，将口腔牙片模拟转换为有色3D动画，对口腔诊疗提供可靠的治疗依据。
Description
技术领域
本发明涉及口腔影像领域，尤其涉及一种口腔影像3D模拟动画实现方法及装置。
背景技术
随着人们口腔健康意识的增强和牙齿美观要求的提高，口腔治疗/牙齿矫正需求量越来越大，长期以来口腔临床医师依靠肉眼观察患者牙片(CT片)，来估算口腔骨密度的高低，特别是种植类治疗，必须对牙槽骨的骨密度有清晰了解，医师个人的主观因素占很大比例，以估算的口腔骨密度分布作为治疗依据，是非常粗糙的治疗方式，而牙片对于没有牙科医学背景的普通群众来讲是非常难读懂的，即现有的诊疗模式为患者无法得知自己的牙病情况和接收治疗后的预期效果，不得不被动接受治疗。
发明内容
鉴于以上内容，本发明提供一种口腔影像3D模拟动画实现方法及装置，能够将牙片转换成普通群众所能理解的3D动画，本发明的技术方案如下：
一方面，本发明提供了一种口腔影像3D模拟动画实现方法，所述口腔影像3D模拟动画实现方法包括以下步骤：
S1、获取待模拟口腔的不同角度的多个口腔影像信息，所述口腔影像为二维图像；
S2、对所述不同角度的多个口腔影像信息进行三维容积重建，得到三维动画，所述三维动画信息包括三维动画影像集及其影像序列信息；
S3、对所述口腔影像进行灰度值分析，得到影像的骨密度结果；
S4、根据S2得到的三维动画和S3得到的灰度值分析结果，对所述三维动画影像集中的影像进行色彩转换处理，得到彩色图像集；
S5、对所述彩色图像集中的每一个图像进行偏振成像处理。
进一步地，所述S1包括：利用锥束CT系统按顺序旋转地锥束生成头部区域的二维图像；
所述S2中的所述三维容积重建包括锥束重建。
进一步地，所述S3中的灰度值分析包括图像分解操作、图像重建操作及骨密度标识操作，
其中，所述图像分解操作包括：以所述口腔影像为第一层分解图像，对所述口腔影像进行图像分解，得到多层分解图像，在分解过程中，每一层分解图像对应一个下采样图像和一个细节图像；
所述图像重建操作包括：以最后一层分解图像所对应的下采样图像为起始重建图像，根据上一层分解图像对应的细节图像，对所述起始重建图像进行图像重建，得到重建中间图像；将所述起始重建图像更新为重建中间图像，将最后一层分解图像更新为上一层分解图像，并循环执行中间图像的重建，直至根据第一层分解图像对应的细节图像，对所述起始重建图像进行图像重建，得到最终重建图像；
所述骨密度标识操作包括：根据所述最终重建图像，得到并标识所述口腔影像对应的口腔区域的骨密度信息。
进一步地，所述S5中对一个图像进行偏振成像处理包括：
S51、对观看者的人眼进行追踪，得到中间视点的水平坐标、瞳孔间距、当前像素的深度和焦距信息；
S52、利用以下公式得到对当前图像的位移信息：
其中，xleft是向左位移量，xc是中间视点的水平坐标，tx是瞳孔间距，Z是当前像素的深度，f是焦距；
其中，xright是向右位移量，xc是中间视点的水平坐标，tx是瞳孔间距，Z是当前像素的深度，f是焦距；
S53、根据位移信息，将当前图像生成两个偏振图像，实现3D呈现；
循环执行S51-S53，直至结束3D模拟。
进一步地，在S5之后还包括：
S61、获取口腔影像对应的患者的面部扫描结果，所述面部扫描结果包括不同角度的多个面扫影像；
S62、将所述面部扫描结果与S5得到的偏振成像处理结果进行渲染；
S63、展示渲染后的成像结果。
进一步地，所述图像分解操作包括：
S311、对第一层分解图像进行高斯滤波操作，得到第一滤波图像；
S312、对所述第一滤波图像进行下采样操作，得到所述第一层分解图像对应的下采样图像；
S313、对所述下采样图像进行高斯滤波操作，得到第二滤波图像；
S314、对所述第二滤波图像进行插值操作，使图像尺寸放大到S22操作之前的尺寸，得到插值图像；
S315、将第一层分解图像减去插值图像，得到所述第一层分解图像对应的细节图像；
S316、将所述第一层分解图像对应的下采样图像作为第二层分解图像，同理于第一层分解图像的分解步骤，得到所述第二层分解图像对应的下采样图像和细节图像，循环执行分解步骤，直至得到最后一层分解图像及其对应的下采样图像。
进一步地，所述图像重建操作包括：
S321、对最后一层分解图像所对应的下采样图像进行高斯滤波并上采样操作，得到上采样图像；
S322、获取所述最后一层分解图像的上一层分解图像所对应的细节图像；
S323、将S321得到的上采样图像与S32得到的细节图像相加，得到所述重建中间图像；
S324、判断是否存在上一层分解图像对应的细节图像，如果存在，则执行S325和S326，否则结束重建；
S325、将S321中最后一层分解图像所对应的下采样图像更新为S323得到的重建中间图像，将最后一层分解图像更新为上一层分解图像；
S326、执行S321-S324。
进一步地，所述骨密度标识操作包括：
S331、根据以下计算公式计算口腔影像对应的口腔区域的骨密度值：
Db＝(Hb-Hw)*Ck/(Hk-Hw)，其中，Db为被测量骨质骨密度，Hb为被测骨感兴趣区的CT值，Hk是体模的CT值，Hw为水的CT值，Ck为体模的密度，其中，Hb为根据S4得到的最终重建图像而获得；
S332、根据S51得到的骨密度值，对最终重建图像进行灰度标识，以表示骨密度信息。
进一步地，S316中的所述循环执行分解步骤，直至得到最后一层分解图像及其对应的下采样图像包括：直至某一层分解图像对应的第一滤波图像尺寸为1*1，则停止分解操作，并将所述某一层分解图像的上一层分解图像作为最后一层分解图像；
S315还包括：对第一层分解图像减去插值图像得到的细节图像进行指数函数增强处理，将增强处理结果作为所述第一层分解图像对应的细节图像。
另一方面，本发明一种口腔影像3D模拟动画实现装置，包括以下模块：
口腔影像获取模块，用于获取待模拟口腔的不同角度的多个口腔影像信息，所述口腔影像为二维图像；
三维容积重建模块，用于对所述不同角度的多个口腔影像信息进行三维容积重建，得到三维动画，所述三维动画信息包括三维动画影像集及其影像序列信息；
灰度值分析模块，用于对所述口腔影像进行灰度值分析，得到影像的骨密度结果；
三维色彩模块，用于根据三维容积重建模块得到的三维动画和灰度值分析模块得到的灰度值分析结果，对所述三维动画影像集中的影像进行色彩转换处理，得到彩色图像集；
偏振成像模块，用于对所述彩色图像集中的每一个图像进行偏振成像处理。
进一步地，所述灰度值分析模块包括：
图像分解单元，用于以所述口腔影像为第一层分解图像，对所述口腔影像进行图像分解，得到多层分解图像，在分解过程中，每一层分解图像对应一个下采样图像和一个细节图像；
图像重建单元，用于以最后一层分解图像所对应的下采样图像为起始重建图像，根据上一层分解图像对应的细节图像，对所述起始重建图像进行图像重建，得到重建中间图像；将所述起始重建图像更新为重建中间图像，将最后一层分解图像更新为上一层分解图像，并循环执行中间图像的重建，直至根据第一层分解图像对应的细节图像，对所述起始重建图像进行图像重建，得到最终重建图像；
骨密度标识单元，用于根据所述最终重建图像，得到并标识所述口腔影像对应的口腔区域的骨密度信息。
进一步地，所述口腔影像3D模拟动画实现装置还包括融合模块，所述融合模块包括以下单元：
面扫单元，用于获取口腔影像对应的患者的面部扫描结果，所述面部扫描结果包括不同角度的多个面扫影像；
渲染单元，用于将所述面部扫描结果与所述偏振成像模块得到的偏振成像处理结果进行渲染；
展示单元，用于展示渲染后的成像结果。
本发明具有下列优点：
a.通过图像密度分析，计算并标识图像骨密度值，完成对牙片信息的转换，使牙片信息一目了然；
b.对口腔影像三维容积重建，模拟出3D效果；
c.结合三维容积重建和图像骨密度分析结果，对三维动画中的影像进行颜色处理，使牙片更直观；
d.对三维动画中的图像进行偏振，实现裸眼3D效果；
e.最后加入人脸扫描作渲染融合，使模拟动画更逼真更清楚。
附图说明
图1是本发明实施例提供的口腔影像3D模拟动画实现方法的流程图；
图2是本发明实施例提供的偏振成像处理的方法流程图；
图3是本发明实施例提供的图像分解操作的方法流程图；
图4是本发明实施例提供的图像重建操作的方法流程图；
图5是本发明实施例提供的融合面部信息的方法流程图；
图6是本发明实施例提供的口腔影像3D模拟动画实现装置的模块框图。
具体实施方式
以下结合说明书附图及具体实施例进一步说明本发明的技术方案。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。
实施例
在本发明的一个实施例中，提供了一种口腔影像3D模拟动画实现方法，如图1所示，所述实现方法包括以下流程：
S1、获取待模拟口腔的不同角度的多个口腔影像信息，所述口腔影像为二维图像。
具体地，利用锥束CT系统按顺序旋转地锥束生成头部区域的二维图像。
S2、对所述不同角度的多个口腔影像信息进行三维容积重建，得到三维动画，所述三维动画信息包括三维动画影像集及其影像序列信息。
具体地，所述三维容积重建包括锥束重建。
S3、对所述口腔影像进行灰度值分析，得到影像的骨密度结果。
具体地，所述灰度值分析包括图像分解操作、图像重建操作及骨密度标识操作：
所述图像分解操作包括：以所述口腔影像为第一层分解图像，对所述口腔影像进行图像分解，得到多层分解图像，在分解过程中，每一层分解图像对应一个下采样图像和一个细节图像；
所述图像重建操作包括：以最后一层分解图像所对应的下采样图像为起始重建图像，根据上一层分解图像对应的细节图像，对所述起始重建图像进行图像重建，得到重建中间图像；将所述起始重建图像更新为重建中间图像，将最后一层分解图像更新为上一层分解图像，并循环执行中间图像的重建，直至根据第一层分解图像对应的细节图像，对所述起始重建图像进行图像重建，得到最终重建图像；
所述骨密度标识操作包括：根据所述最终重建图像，得到并标识所述口腔影像对应的口腔区域的骨密度信息。
S4、根据S2得到的三维动画和S3得到的灰度值分析结果，对所述三维动画影像集中的影像进行色彩转换处理，得到彩色图像集。
预先在软件中，写入X光片的灰度级和骨密度的对应值，即根据骨密度可以匹配相应的灰度值，具体地，S4是将S2的三维重建和S4的灰度值融合起来，得到一个带有色彩信息的三维动画。
S5、对所述彩色图像集中的每一个图像进行偏振成像处理，如图2所示，具体包括以下流程：
S51、对观看者的人眼进行追踪，得到中间视点的水平坐标、瞳孔间距、当前像素的深度和焦距信息；
S52、利用以下公式得到对当前图像的位移信息：
其中，xleft是向左位移量，xc是中间视点的水平坐标，tx是瞳孔间距，Z是当前像素的深度，f是焦距；
其中，xright是向右位移量，xc是中间视点的水平坐标，tx是瞳孔间距，Z是当前像素的深度，f是焦距；
S53、根据位移信息，将当前图像生成两个偏振图像，实现3D呈现；
循环执行S51-S53，直至结束3D模拟。
以下针对灰度值分析包括的图像分解操作、图像重建操作及骨密度标识操作进行具体的介绍：
参见图3，所述图像分解操作包括以下流程：
S311、对第一层分解图像进行5*5高斯滤波操作，得到第一滤波图像；
S312、对所述第一滤波图像进行下采样操作，得到所述第一层分解图像对应的下采样图像；
S313、对所述下采样图像进行5*5高斯滤波操作，得到第二滤波图像；
S314、对所述第二滤波图像进行插值操作，使图像尺寸放大到S22操作之前的尺寸，得到插值图像；
S315、将第一层分解图像减去插值图像，得到所述第一层分解图像对应的细节图像，优选地，对第一层分解图像减去插值图像得到的细节图像进行指数函数增强处理，将增强处理结果作为所述第一层分解图像对应的细节图像，这样，后续重建得到的就是增强图像；
S316、将所述第一层分解图像对应的下采样图像作为第二层分解图像，同理于第一层分解图像的分解步骤，得到所述第二层分解图像对应的下采样图像和细节图像，循环执行分解步骤，直至得到最后一层分解图像及其对应的下采样图像，即直至某一层分解图像对应的第一滤波图像尺寸为1*1(不满足继续分解条件的情况)，则停止分解操作，并将所述某一层分解图像的上一层分解图像作为最后一层分解图像。
参见图4，所述图像重建操作包括以下流程：
S321、对最后一层分解图像所对应的下采样图像进行5*5高斯滤波并上采样操作，得到上采样图像；
S322、获取所述最后一层分解图像的上一层分解图像所对应的细节图像；
S323、将S321得到的上采样图像与S322得到的细节图像相加，得到所述重建中间图像；
S324、判断是否存在上一层分解图像对应的细节图像，如果存在，则执行S325，否则结束重建；
S325、将S321中最后一层分解图像所对应的下采样图像更新为S323得到的重建中间图像，将最后一层分解图像更新为上一层分解图像并重复执行S321-S324。
另外，所述骨密度标识操作包括以下流程：
S331、根据以下计算公式计算口腔影像对应的口腔区域的骨密度值：
Db＝(Hb-Hw)*Ck/(Hk-Hw)，其中，Db为被测量骨质骨密度，Hb为被测骨感兴趣区的CT值，Hk是体模的CT值，Hw为水的CT值，Ck为体模的密度，其中，Hb为根据上述操作得到的最终重建图像而获得；
S332、根据S51得到的骨密度值，对最终重建图像进行灰度标识，以表示骨密度信息。
在本发明的一个优选实施例中，在偏振模拟3D成像之后还包括融合面部信息，参见图5，包括以下流程：
S61、获取口腔影像对应的患者的面部扫描结果，所述面部扫描结果包括不同角度的多个面扫影像；
S62、将所述面部扫描结果与S5得到的偏振成像处理结果进行渲染；
S63、展示渲染后的成像结果。
以上渲染的效果即为能够更高像素地向用户展示人脸及口腔，效果更形象更逼真。
另一方面，本发明一种口腔影像3D模拟动画实现装置，参见图6，所述装置包括以下模块：
口腔影像获取模块610，用于获取待模拟口腔的不同角度的多个口腔影像信息，所述口腔影像为二维图像；
三维容积重建模块620，用于对所述不同角度的多个口腔影像信息进行三维容积重建，得到三维动画，所述三维动画信息包括三维动画影像集及其影像序列信息；
灰度值分析模块630，用于对所述口腔影像进行灰度值分析，得到影像的骨密度结果；
三维色彩模块640，用于根据三维容积重建模块620得到的三维动画和灰度值分析模块630得到的灰度值分析结果，对所述三维动画影像集中的影像进行色彩转换处理，得到彩色图像集；
偏振成像模块650，用于对所述彩色图像集中的每一个图像进行偏振成像处理。
进一步地，所述灰度值分析模块630包括：
图像分解单元631，用于以所述口腔影像为第一层分解图像，对所述口腔影像进行图像分解，得到多层分解图像，在分解过程中，每一层分解图像对应一个下采样图像和一个细节图像；
图像重建单元632，用于以最后一层分解图像所对应的下采样图像为起始重建图像，根据上一层分解图像对应的细节图像，对所述起始重建图像进行图像重建，得到重建中间图像；将所述起始重建图像更新为重建中间图像，将最后一层分解图像更新为上一层分解图像，并循环执行中间图像的重建，直至根据第一层分解图像对应的细节图像，对所述起始重建图像进行图像重建，得到最终重建图像；
骨密度标识单元633，用于根据所述最终重建图像，得到并标识所述口腔影像对应的口腔区域的骨密度信息。
进一步地，所述口腔影像3D模拟动画实现装置还包括融合模块660，所述融合模块660包括以下单元：
面扫单元661，用于获取口腔影像对应的患者的面部扫描结果，所述面部扫描结果包括不同角度的多个面扫影像；
渲染单元662，用于将所述面部扫描结果与所述偏振成像模块650得到的偏振成像处理结果进行渲染；
展示单元663，用于展示渲染后的成像结果。
需要说明的是：上述实施例提供的口腔影像3D模拟动画实现装置在进行动画模拟时，仅以上述各功能模块的划分进行举例说明，实际应用中，可以根据需要而将上述功能分配由不同的功能模块完成，即将口腔影像3D模拟动画实现装置的内部结构划分成不同的功能模块，以完成以上描述的全部或者部分功能。另外，本实施例提供的口腔影像3D模拟动画实现装置实施例与上述实施例提供的口腔影像3D模拟动画实现方法属于同一构思，其具体实现过程详见方法实施例，这里不再赘述。
本发明基于数字化牙片的骨密度测量技术及计算机辅助图像密度分析技术，将口腔牙片模拟转换为有色3D动画，对口腔诊疗提供可靠的治疗依据。
以上所述仅为本发明的优选实施例，并非因此限制其专利范围，凡是利用本发明说明书及附图内容所作的等效结构或等效流程变换，直接或间接运用在其他相关的技术领域，均同理包括在本发明的专利保护范围内。
Claims (12)
1.一种口腔影像3D模拟动画实现方法，其特征在于，包括以下步骤：
S1、获取待模拟口腔的不同角度的多个口腔影像信息，所述口腔影像为二维图像；
S2、对所述不同角度的多个口腔影像信息进行三维容积重建，得到三维动画，所述三维动画信息包括三维动画影像集及其影像序列信息；
S3、对所述口腔影像进行灰度值分析，得到影像的骨密度结果；
S4、根据S2得到的三维动画和S3得到的灰度值分析结果，对所述三维动画影像集中的影像进行色彩转换处理，得到彩色图像集；
S5、对所述彩色图像集中的每一个图像进行偏振成像处理。
2.根据权利要求1所述的口腔影像3D模拟动画实现方法，其特征在于，所述S1包括：利用锥束CT系统按顺序旋转地锥束生成头部区域的二维图像；
所述S2中的所述三维容积重建包括锥束重建。
3.根据权利要求1所述的口腔影像3D模拟动画实现方法，其特征在于，所述S3中的灰度值分析包括图像分解操作、图像重建操作及骨密度标识操作，
其中，所述图像分解操作包括：以所述口腔影像为第一层分解图像，对所述口腔影像进行图像分解，得到多层分解图像，在分解过程中，每一层分解图像对应一个下采样图像和一个细节图像；
所述图像重建操作包括：以最后一层分解图像所对应的下采样图像为起始重建图像，根据上一层分解图像对应的细节图像，对所述起始重建图像进行图像重建，得到重建中间图像；将所述起始重建图像更新为重建中间图像，将最后一层分解图像更新为上一层分解图像，并循环执行中间图像的重建，直至根据第一层分解图像对应的细节图像，对所述起始重建图像进行图像重建，得到最终重建图像；
所述骨密度标识操作包括：根据所述最终重建图像，得到并标识所述口腔影像对应的口腔区域的骨密度信息。
4.根据权利要求1所述的口腔影像3D模拟动画实现方法，其特征在于，所述S5中对一个图像进行偏振成像处理包括：
S51、对观看者的人眼进行追踪，得到中间视点的水平坐标、瞳孔间距、当前像素的深度和焦距信息；
S52、利用以下公式得到对当前图像的位移信息：
其中，xleft是向左位移量，xc是中间视点的水平坐标，tx是瞳孔间距，Z是当前像素的深度，f是焦距；
其中，xright是向右位移量，xc是中间视点的水平坐标，tx是瞳孔间距，Z是当前像素的深度，f是焦距；
S53、根据位移信息，将当前图像生成两个偏振图像，实现3D呈现；
循环执行S51-S53，直至结束3D模拟。
5.根据权利要求1所述的口腔影像3D模拟动画实现方法，其特征在于，在S5之后还包括：
S61、获取口腔影像对应的患者的面部扫描结果，所述面部扫描结果包括不同角度的多个面扫影像；
S62、将所述面部扫描结果与S5得到的偏振成像处理结果进行渲染；
S63、展示渲染后的成像结果。
6.根据权利要求3所述的口腔影像3D模拟动画实现方法，其特征在于，所述图像分解操作包括：
S311、对第一层分解图像进行高斯滤波操作，得到第一滤波图像；
S312、对所述第一滤波图像进行下采样操作，得到所述第一层分解图像对应的下采样图像；
S313、对所述下采样图像进行高斯滤波操作，得到第二滤波图像；
S314、对所述第二滤波图像进行插值操作，使图像尺寸放大到S22操作之前的尺寸，得到插值图像；
S315、将第一层分解图像减去插值图像，得到所述第一层分解图像对应的细节图像；
S316、将所述第一层分解图像对应的下采样图像作为第二层分解图像，同理于第一层分解图像的分解步骤，得到所述第二层分解图像对应的下采样图像和细节图像，循环执行分解步骤，直至得到最后一层分解图像及其对应的下采样图像。
7.根据权利要求3所述的口腔影像3D模拟动画实现方法，其特征在于，所述图像重建操作包括：
S321、对最后一层分解图像所对应的下采样图像进行高斯滤波并上采样操作，得到上采样图像；
S322、获取所述最后一层分解图像的上一层分解图像所对应的细节图像；
S323、将S321得到的上采样图像与S32得到的细节图像相加，得到所述重建中间图像；
S324、判断是否存在上一层分解图像对应的细节图像，如果存在，则执行S325和S326，否则结束重建；
S325、将S321中最后一层分解图像所对应的下采样图像更新为S323得到的重建中间图像，将最后一层分解图像更新为上一层分解图像；
S326、执行S321-S324。
8.根据权利要求3所述的口腔影像3D模拟动画实现方法，其特征在于，所述骨密度标识操作包括：
S331、根据以下计算公式计算口腔影像对应的口腔区域的骨密度值：
Db＝(Hb-Hw)*Ck/(Hk-Hw)，其中，Db为被测量骨质骨密度，Hb为被测骨感兴趣区的CT值，Hk是体模的CT值，Hw为水的CT值，Ck为体模的密度，其中，Hb为根据S4得到的最终重建图像而获得；
S332、根据S51得到的骨密度值，对最终重建图像进行灰度标识，以表示骨密度信息。
9.根据权利要求6所述的口腔影像3D模拟动画实现方法，其特征在于，S316中的所述循环执行分解步骤，直至得到最后一层分解图像及其对应的下采样图像包括：直至某一层分解图像对应的第一滤波图像尺寸为1*1，则停止分解操作，并将所述某一层分解图像的上一层分解图像作为最后一层分解图像；
S315还包括：对第一层分解图像减去插值图像得到的细节图像进行指数函数增强处理，将增强处理结果作为所述第一层分解图像对应的细节图像。
10.一种口腔影像3D模拟动画实现装置，其特征在于，包括以下模块：
口腔影像获取模块，用于获取待模拟口腔的不同角度的多个口腔影像信息，所述口腔影像为二维图像；
三维容积重建模块，用于对所述不同角度的多个口腔影像信息进行三维容积重建，得到三维动画，所述三维动画信息包括三维动画影像集及其影像序列信息；
灰度值分析模块，用于对所述口腔影像进行灰度值分析，得到影像的骨密度结果；
三维色彩模块，用于根据三维容积重建模块得到的三维动画和灰度值分析模块得到的灰度值分析结果，对所述三维动画影像集中的影像进行色彩转换处理，得到彩色图像集；
偏振成像模块，用于对所述彩色图像集中的每一个图像进行偏振成像处理。
11.根据权利要求10所述的口腔影像3D模拟动画实现装置，其特征在于，所述灰度值分析模块包括：
图像分解单元，用于以所述口腔影像为第一层分解图像，对所述口腔影像进行图像分解，得到多层分解图像，在分解过程中，每一层分解图像对应一个下采样图像和一个细节图像；
图像重建单元，用于以最后一层分解图像所对应的下采样图像为起始重建图像，根据上一层分解图像对应的细节图像，对所述起始重建图像进行图像重建，得到重建中间图像；将所述起始重建图像更新为重建中间图像，将最后一层分解图像更新为上一层分解图像，并循环执行中间图像的重建，直至根据第一层分解图像对应的细节图像，对所述起始重建图像进行图像重建，得到最终重建图像；
骨密度标识单元，用于根据所述最终重建图像，得到并标识所述口腔影像对应的口腔区域的骨密度信息。
12.根据权利要求10所述的口腔影像3D模拟动画实现装置，其特征在于，还包括融合模块，所述融合模块包括以下单元：
面扫单元，用于获取口腔影像对应的患者的面部扫描结果，所述面部扫描结果包括不同角度的多个面扫影像；
渲染单元，用于将所述面部扫描结果与所述偏振成像模块得到的偏振成像处理结果进行渲染；
展示单元，用于展示渲染后的成像结果。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810520990.XA CN108830915B (zh) | 2018-05-28 | 2018-05-28 | 口腔影像3d模拟动画实现方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810520990.XA CN108830915B (zh) | 2018-05-28 | 2018-05-28 | 口腔影像3d模拟动画实现方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108830915A true CN108830915A (zh) | 2018-11-16 |
| CN108830915B CN108830915B (zh) | 2023-03-31 |
Family
ID=64145832
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810520990.XA Active CN108830915B (zh) | 2018-05-28 | 2018-05-28 | 口腔影像3d模拟动画实现方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108830915B (zh) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110097944A (zh) * | 2019-04-11 | 2019-08-06 | 山东大学齐鲁医院 | 一种人体器官模型的显示调控方法及系统 |
| CN110570508A (zh) * | 2019-07-23 | 2019-12-13 | 博志生物科技(深圳)有限公司 | 骨质疏松状况三维可视化渲染方法 |
| CN113554576A (zh) * | 2020-04-24 | 2021-10-26 | 上海联影医疗科技股份有限公司 | 一种多期相数据的减影方法、装置、设备及存储介质 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101030298A (zh) * | 2007-03-29 | 2007-09-05 | 杭州电子科技大学 | 多尺度自适应对比度变换的医学图像增强方法 |
| CN102215779A (zh) * | 2008-11-18 | 2011-10-12 | Ibur有限责任公司 | 用于链接物理和数字数据以用于诊断、治疗计划、患者教育、通信、制造和数据传送目的的牙科装置和方法 |
| TW201813596A (zh) * | 2016-10-06 | 2018-04-16 | 美佳境植牙股份有限公司 | 用於植體診斷的圖像生成系統及其生成方法 |
| TW201813587A (zh) * | 2016-10-06 | 2018-04-16 | 美佳境植牙股份有限公司 | 移動鏈接植體診斷系統 |
-
2018
- 2018-05-28 CN CN201810520990.XA patent/CN108830915B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101030298A (zh) * | 2007-03-29 | 2007-09-05 | 杭州电子科技大学 | 多尺度自适应对比度变换的医学图像增强方法 |
| CN102215779A (zh) * | 2008-11-18 | 2011-10-12 | Ibur有限责任公司 | 用于链接物理和数字数据以用于诊断、治疗计划、患者教育、通信、制造和数据传送目的的牙科装置和方法 |
| TW201813596A (zh) * | 2016-10-06 | 2018-04-16 | 美佳境植牙股份有限公司 | 用於植體診斷的圖像生成系統及其生成方法 |
| TW201813587A (zh) * | 2016-10-06 | 2018-04-16 | 美佳境植牙股份有限公司 | 移動鏈接植體診斷系統 |
Non-Patent Citations (1)
| Title |
|---|
| 杨登峰 等: ""重组人骨保护素对激素性股骨头坏死患者骨密度及髋关节Harris评分的影响"", 《现代生物医学进展》 * |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110097944A (zh) * | 2019-04-11 | 2019-08-06 | 山东大学齐鲁医院 | 一种人体器官模型的显示调控方法及系统 |
| CN110097944B (zh) * | 2019-04-11 | 2021-05-07 | 山东大学齐鲁医院 | 一种人体器官模型的显示调控方法及系统 |
| CN110570508A (zh) * | 2019-07-23 | 2019-12-13 | 博志生物科技(深圳)有限公司 | 骨质疏松状况三维可视化渲染方法 |
| CN110570508B (zh) * | 2019-07-23 | 2023-09-08 | 博志生物科技(深圳)有限公司 | 骨质疏松状况三维可视化渲染方法 |
| CN113554576A (zh) * | 2020-04-24 | 2021-10-26 | 上海联影医疗科技股份有限公司 | 一种多期相数据的减影方法、装置、设备及存储介质 |
| CN113554576B (zh) * | 2020-04-24 | 2024-01-19 | 上海联影医疗科技股份有限公司 | 一种多期相数据的减影方法、装置、设备及存储介质 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108830915B (zh) | 2023-03-31 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Tzou et al. | Comparison of three-dimensional surface-imaging systems | |
| EP3490489B1 (en) | Method and system for removing orthodontic brackets from the digital mesh of the dentition | |
| EP3276575A1 (en) | Method, apparatus and system for reconstructing images of 3d surface | |
| CN106901834A (zh) | 微创心脏外科手术的术前规划及手术虚拟现实模拟方法 | |
| CN106874700A (zh) | 基于Web的手术模拟方法、手术模拟装置及电子设备 | |
| Schendel et al. | 3D orthognathic surgery simulation using image fusion | |
| Sforza et al. | Soft-tissue facial anthropometry in three dimensions: from anatomical landmarks to digital morphology in research, clinics and forensic anthropology | |
| CN109157284A (zh) | 一种脑肿瘤医学影像三维重建显示交互方法及系统 | |
| CN108830915A (zh) | 口腔影像3d模拟动画实现方法及装置 | |
| CN110021445A (zh) | 一种基于vr模型的医疗系统 | |
| CN105718730A (zh) | 对受试者的疼痛量化评估方法以及实施该方法的系统 | |
| Buchaillard et al. | 3D statistical models for tooth surface reconstruction | |
| Robles et al. | A step-by-step method for producing 3D crania models from CT data | |
| Advincula et al. | Development and future trends in the application of visualization toolkit (VTK): the case for medical image 3D reconstruction | |
| WO2023274690A1 (en) | Non-invasive periodontal examination | |
| CN108805950A (zh) | 口腔牙片信息转换方法及装置 | |
| JP2007289569A (ja) | 医用画像処理装置及び医用画像処理方法 | |
| JP2015154914A (ja) | 画像処理装置、画像処理装置の制御方法およびプログラム | |
| Yang et al. | ImplantFormer: vision transformer-based implant position regression using dental CBCT data | |
| CN108877897A (zh) | 牙科诊疗方案生成方法、装置及诊疗系统 | |
| Enciso et al. | Three-dimensional head anthropometric analysis | |
| KR101657285B1 (ko) | 초음파 치료 시뮬레이션시스템 | |
| CN101176683A (zh) | 人体立体解剖图象的制作方法与应用 | |
| Khatri et al. | Unfolding the mysterious path of forensic facial reconstruction: Review of different imaging modalities | |
| Yu et al. | Deep learning single view computed tomography guided by FBP algorithm |
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 | ||
| CP03 | Change of name, title or address |
Address after: Floor 12, No. 609, Xiyuan Road, Suzhou City, Jiangsu Province, 215000 Patentee after: Dental Doctor Medical Holding Group Co.,Ltd. Country or region after: China Address before: Floor 12, No. 609, Xiyuan Road, Suzhou City, Jiangsu Province, 215008 Patentee before: DENTAL DOCTOR MEDICAL HOLDING GROUP CO.,LTD. Country or region before: China |
|
| CP03 | Change of name, title or address |