CN107782791A - 一种用于牙体的声显微成像方法 - Google Patents

一种用于牙体的声显微成像方法 Download PDF

Info

Publication number
CN107782791A
CN107782791A CN201710785542.8A CN201710785542A CN107782791A CN 107782791 A CN107782791 A CN 107782791A CN 201710785542 A CN201710785542 A CN 201710785542A CN 107782791 A CN107782791 A CN 107782791A
Authority
CN
China
Prior art keywords
mrow
msub
msubsup
mtr
mtd
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.)
Pending
Application number
CN201710785542.8A
Other languages
English (en)
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.)
AVIC TECHNOLOGY FOUNDATION ESTABLISHMENT
Original Assignee
AVIC TECHNOLOGY FOUNDATION ESTABLISHMENT
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 AVIC TECHNOLOGY FOUNDATION ESTABLISHMENT filed Critical AVIC TECHNOLOGY FOUNDATION ESTABLISHMENT
Priority to CN201710785542.8A priority Critical patent/CN107782791A/zh
Publication of CN107782791A publication Critical patent/CN107782791A/zh
Pending legal-status Critical Current

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
    • G01N29/04—Analysing solids
    • G01N29/06—Visualisation of the interior, e.g. acoustic microscopy
    • G01N29/0654—Imaging
    • G01N29/0681—Imaging by acoustic microscopy, e.g. scanning acoustic microscopy
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2291/00—Indexing codes associated with group G01N29/00
    • G01N2291/02—Indexing codes associated with the analysed material
    • G01N2291/024—Mixtures
    • G01N2291/02483—Other human or animal parts, e.g. bones

Landscapes

  • Physics & Mathematics (AREA)
  • Acoustics & Sound (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Medicines Containing Antibodies Or Antigens For Use As Internal Diagnostic Agents (AREA)

Abstract

本发明属于医学及无损检测技术领域，涉及一种用于牙体的声显微成像方法。本发明方法利用高分辨率宽带窄脉冲声波与被观察的牙体相互作用，产生的声学信息与牙体组织的数理联系，利用专门的声学扫描单元驱动声学透镜，使其对被检测的牙体样本进行格栅扫描，利用声学透镜实时提取被检测的牙体内部的声学信息，通过对声学透镜获取的声波信号和对应的位置信号重构，进行声显微成像显示，实现对牙体表面及其内部组织的声显微成像。试验应用结果表明，采用所构建的牙体声显微成像方法，可以快速清晰的得到牙体内部釉质、牙本质、牙髓、牙骨质等组织及其变化的特征分布，取得了较好的实际声显微成像效果。

Description

一种用于牙体的声显微成像方法
技术领域
本发明属于医学及无损检测技术领域，涉及一种用于牙体的声显微成像方法。
背景技术
牙体组织的变化是导致牙腔疾病的主要根源，通常当牙体组织发生了明显的病变时，可能会造成牙的功能失效。目前，临床上，主要是由主治医生采用目视或辅助Micro-X射线CT等方法进行牙体损坏或牙体疾病进行检查和诊断。其主要不足是：1)目视检查方法，得不到牙体组织内部信息，特别是难以掌握到一些早期牙体组织的内部变化或病变方面的信息，从而可能会失去最佳的牙病治疗时机；2)Micro-X射线CT方法，对患者辐射损害明显，也难以得到早期牙体组织的内部变化或病变方面的信息。对于离体牙，目前主要是采用切片后的光学观察的方法。其主要不足是：1)只能观察到样品表面存在光学灰度差的组织及其变化；2)不能得到牙体样品内部牙体组织和早期牙体组织的内部变化或病变方面的信息。
发明内容
本发明的目的是针对现有技术中的不足，提出一种基于来自牙体中的声学信息的声显微成像方法。
本发明的技术解决方案是，采用由声学透镜、声学激励/接收单元、声信息处理单元、声学扫描机构、声学扫描控制单元、声学成像单元、牙体样本固定平台组成的声显微成像装置进行声显微成像，其方法是，
(1)声学扫描
①将被观察的牙体样本放置在声显微成像装置的牙体样本固定平台中的声学扫描台容器内，并使牙体样本中被观察面置于声学透镜的下方；
②根据牙体样本在扫描平面内的最大投影面积S0，确定声扫描域Su的大小，
且使Su＞S0，
这里，分别为声透镜在x、y方向的最大扫描范围，且
其中，—分别为声学透镜在x方向的扫描终点和起点，
—分别为声学透镜在y方向的扫描终点和起点，
③根据被观察牙体样本的实际大小，根据式(1)，在声学扫描控制单元中设置牙体样本的(x,y)方向扫描范围；
④设置(x,y)方向声学扫描分辨率Δx、Δy，且使Δx＝Δy；
⑤设置声学透镜的扫描速度υ，
扫描分辨率按式(2)确定，
式中，Δt—为每点声学信息和位置信息采集和成像显示所需时间；
⑥扫描成像，按照所设置和选择的声学扫描参数，对牙体样本进行声学扫描；
(2)声学成像
利用声学透镜提取扫描过程中第i行第j列位置点的反射声学信息或透射声学信息对牙体样本进行声学成像,
1)反射法声学成像
①利用声学透镜提取扫描过程中第i行第j列位置点的反射声学信息
利用反射声学信息阵列表示扫描域Su内的各个检测位置的反射声学信息：
式中，
—对应牙体样本中声学反射系数，与牙体组织和牙体均匀性有关，
A0—为牙体样本中入射声波信号的幅值，
α—为入射声波在牙体样本中的声衰减系数，
ν—为入射声波在牙体样本中的声速，
—为来自第i行第j列位置牙体组织中的反射声波的传播时间，
i、j—分别表示第i行第j列位置序列，i＝1,2,...,m,j＝1,2,...,n,
m、n—分别表示反射声学信息阵列的行数和列数；
②利用声学扫描控制单元提取扫描域Su内第i行第j列位置的检测位置坐标pij(xi,yj)，利用位置阵列Pm×n表示扫描域Su内各个位置点的坐标位置：
③利用声信息处理单元对进行数字化后，按照式(6)转换为成像信号dkl(rkl,gkl,bkl)，此时：
这里，rkl,gkl,bkl—分别对应每个成像点的颜色值，
kr,kg,kb—分别为信号幅值与rkl,gkl,bkl颜色值之间的转换系数，
k、l—分别对应成像域中的位置坐标序列，且取
式(7)中，Δnx、Δny—分别为扫描域Su中每个位置点在成像域Sp中的像素点数，此时若用Δp表示Su中每个位置点在成像域Sp中的像素大小，则有：
Δp＝Δnx×Δny×Δs
这里，Δs—为声学成像单元中所用图像显示器的像素点的物理大小，
将扫描域Su中每个检测点位置坐标pij(xi,yj)映射到成像域Sp中位置坐标qkl(xk,yl)，并结合式(4)和式(5)，对每个点成像信号dkl(rkl,gkl,bkl)和成像位置坐标qkl(xk,yl)进行融合后，将对应扫描区域Su的成像结果在成像域Sp中通过矩阵QM×N表示，即有：
式中(8)，M＝m×Δnx，N＝n×Δny
2)透射法声学成像
①利用声学透镜提取扫描过程中第i行第j列位置点的透射声波信息
利用透射声学信息阵列表示扫描域Su内的各个检测位置的透射声学信息：
式中，
—对应牙体样本中声学透射系数，与牙体组织和牙体均匀性有关，
α—为入射声波在牙体样本中的声衰减系数，
ν—为入射声波在牙体样本中的声速，
—为来自第i行第j列位置牙体组织中的透射声波的传播时间，
i、j—分别表示第i行第j列位置序列，i＝1,2,...,m,j＝1,2,...,n,
m、n—分别表示透射声学信息阵列的行数和列数；
②利用声学扫描控制单元提取扫描域Su内第i行第j列位置的检测位置坐标pij(xi,yj)，利用式(5)中的位置阵列Pm×n表示扫描域Su内各个位置点的坐标位置，
③利用声信息处理单元对进行数字化后，按照式(11)转换为成像信号dkl(rkl,gkl,bkl)，此时：
将扫描域Su中每个检测点位置坐标pij(xi,yj)映射到成像域Sp中位置坐标qkl(xk,yl)，并结合式(5)和式(10)，对每个点成像信号dkl(rkl,gkl,bkl)和成像位置坐标qkl(xk,yl)进行融合后，利用式(8)将对应扫描区域Su的成像结果在成像域Sp中通过矩阵QM×N表示。
(3)扫描成像显示
根据设定的扫描成像参数，声学透镜在声学激励/接收单元作用下，通过声学扫描机构和声学扫描控制单元使声学透镜在扫描域Su内对牙体样本进行扫描，由声信息处理单元和声学扫描控制单元实时获取扫描域Su中每个位置点的反射声学信息信号和位置坐标pij(xi,yj)信号，然后由声学成像单元按照式(6、7、8)中的映射关系，转换为成像信号在声学成像单元中创建的成像颜色板中进行成像显示，或由声信息处理单元和声学扫描控制单元实时获取扫描域Su中每个位置点的反射声学信息信号和位置坐标pij(xi,yj)信号，然后由声学成像单元按照式(7、8、11)中的映射关系，转换为成像信号在声学成像单元中创建的成像颜色板中进行成像显示，成像颜色板的大小由成像域Sp按照式(12)确定，
Sp＝M×N (12)。
所述的扫描域Su的确定方法是，
Su＝(2.0～3.3)So
这里，So为被观察牙体样本在扫描平面内的最大投影面积，
—分别为被观察牙体样本最大投影面积在x方向的终点和起点，
—分别为被观察牙体样本最大投影面积在y方向的终点和起点，
使被观察的牙体样本位于最佳的声学扫描成像区。
所述声学成像中的的确定方法是，
所述的声学透镜的焦点大小，根据被观察牙体样本的特征和观察分析要求，在10微米-100微米范围内。
本发明具有的优点和有益效果，
本发明利用高分辨率脉冲声波与被观察的牙体相互作用，产生的声学信息与牙体内部不同组织部位的数理联系，通过专门的声学扫描单元驱动声学透镜，使之对被检测的牙体样本进行光栅式扫描，并实时提取每个扫描位置点的声学信息和位置信息，通过对声学透镜接收到的声学信息和对应的位置信息进行重构，实现对牙体表面及其内部组织的声显微成像显示，进而对牙体表面及牙体内部组织和牙体组织变化或牙体缺陷进行快速声显微成像，为牙体疾病和牙体早期组织变化的临床研究、诊断、牙病预防等提供一种快速环保的声显微成像分析方法。
本发明利用牙体组织及其变化与其声弹性特性，进行牙体组织声显微成像，可以得到牙体组织表面及其内部生物组织结构和缺陷的声图像信息，为研究牙体组织的内部变化或病变的诊断提供了一种环保无害的可视化方法，从而有利于及时掌握一些早期牙体组织的内部变化或病变方面的信息，及时制定最佳的牙病治疗时机；对于离体牙，可以观察到牙体表面及其内部生物组织的变化或病变等方面的信息，而对被观察牙体样本没有任何生物破坏。利用本发明，可以用于快速方便地实现对牙体表面、牙体内部组织及其变化以及牙体缺陷的声显微成像分析，同时也为临床牙体疾病的诊断和检查提供了一种快速环保的可视化成像方法。
附图说明
图1是本发明的牙体声显微成像方法原理示意图。
图2是本发明的牙体声显微成像实现步骤。
具体实施方式
下面结合附图对本发明作详细说明
采用由声学透镜1、声学激励/接收单元2、声信息处理单元3、声学扫描机构4、声学扫描控制单元5、声学成像单元6、牙体样本固定平台7组成的声显微成像装置对牙体组织进行声显微成像，参见图1所示。
进行声显微成像方法的流程如图2所示，
(1)声学扫描
①将被观察的牙体样本8放置在声显微成像装置的牙体样本固定平台7中的声学扫描台容器内，并使牙体样本8中被观察面置于声学透镜1的下方；
②根据牙体样本8在扫描平面内的最大投影面积S0，确定声扫描域Su的大小，
且使Su＞S0，
这里，分别为声透镜1在x、y方向的最大扫描范围，且
其中，—分别为声学透镜1在x方向的扫描终点和起点，
—分别为声学透镜1在y方向的扫描终点和起点，
③根据被观察牙体样本8的实际大小，根据式(1)，在声学扫描控制单元5中设置牙体样本8的(x,y)方向扫描范围；
④设置(x,y)方向声学扫描分辨率Δx、Δy，且使Δx＝Δy；
⑤设置声学透镜1的扫描速度υ，
扫描分辨率按式(2)确定，
式中，Δt—为每点声学信息和位置信息采集和成像显示所需时间；
⑥扫描成像，按照所设置和选择的声学扫描参数，对牙体样本8进行声学扫描；
(2)声学成像
利用声学透镜1提取扫描过程中第i行第j列位置点的反射声学信息或透射声学信息对牙体样本8进行声学成像,
1)反射法声学成像
①利用声学透镜1提取扫描过程中第i行第j列位置点的反射声学信息
利用反射声学信息阵列表示扫描域Su内的各个检测位置的反射声学信息：
式中，
—对应牙体样本中声学反射系数，与牙体组织和牙体均匀性有关，
A0—为牙体样本中入射声波信号的幅值，
α—为入射声波在牙体样本中的声衰减系数，
ν—为入射声波在牙体样本中的声速，
—为来自第i行第j列位置牙体组织中的反射声波的传播时间，
i、j—分别表示第i行第j列位置序列，i＝1,2,...,m,j＝1,2,...,n,
m、n—分别表示反射声学信息阵列的行数和列数；
②利用声学扫描控制单元5提取扫描域Su内第i行第j列位置的检测位置坐标pij(xi,yj)，利用位置阵列Pm×n表示扫描域Su内各个位置点的坐标位置：
③利用声信息处理单元3对进行数字化后，按照式(6)转换为成像信号dkl(rkl,gkl,bkl)，此时：
这里，rkl,gkl,bkl—分别对应每个成像点的颜色值，
kr,kg,kb—分别为信号幅值与rkl,gkl,bkl颜色值之间的转换系数，
k、l—分别对应成像域中的位置坐标序列，且取
式(7)中，Δnx、Δny—分别为扫描域Su中每个位置点在成像域Sp中的像素点数，此时若用Δp表示Su中每个位置点在成像域Sp中的像素大小，则有：
Δp＝Δnx×Δny×Δs
这里，Δs—为声学成像单元6中所用图像显示器的像素点的物理大小，
将扫描域Su中每个检测点位置坐标pij(xi,yj)映射到成像域Sp中位置坐标qkl(xk,yl)，并结合式(4)和式(5)，对每个点成像信号dkl(rkl,gkl,bkl)和成像位置坐标qkl(xk,yl)进行融合后，将对应扫描区域Su的成像结果在成像域Sp中通过矩阵QM×N表示，即有：
式中(8)，M＝m×Δnx，N＝n×Δny
2)透射法声学成像
①利用声学透镜1提取扫描过程中第i行第j列位置点的透射声波信息
利用透射声学信息阵列表示扫描域Su内的各个检测位置的透射声学信息：
式中，
—对应牙体样本中声学透射系数，与牙体组织和牙体均匀性有关，
α—为入射声波在牙体样本中的声衰减系数，
ν—为入射声波在牙体样本中的声速，
—为来自第i行第j列位置牙体组织中的透射声波的传播时间，
i、j—分别表示第i行第j列位置序列，i＝1,2,...,m,j＝1,2,...,n,
m、n—分别表示透射声学信息阵列的行数和列数；
②利用声学扫描控制单元5提取扫描域Su内第i行第j列位置的检测位置坐标pij(xi,yj)，利用式(5)中的位置阵列Pm×n表示扫描域Su内各个位置点的坐标位置，
③利用声信息处理单元3对进行数字化后，按照式(11)转换为成像信号dkl(rkl,gkl,bkl)，此时：
将扫描域Su中每个检测点位置坐标pij(xi,yj)映射到成像域Sp中位置坐标qkl(xk,yl)，并结合式(5)和式(10)，对每个点成像信号dkl(rkl,gkl,bkl)和成像位置坐标qkl(xk,yl)进行融合后，利用式(8)将对应扫描区域Su的成像结果在成像域Sp中通过矩阵QM×N表示。
(3)扫描成像显示
根据设定的扫描成像参数，声学透镜1在声学激励/接收单元2作用下，通过声学扫描机构4和声学扫描控制单元5使声学透镜1在扫描域Su内对牙体样本进行扫描，由声信息处理单元3和声学扫描控制单元5实时获取扫描域Su中每个位置点的反射声学信息信号和位置坐标pij(xi,yj)信号，然后由声学成像单元6按照式(6、7、8)中的映射关系，转换为成像信号在声学成像单元6中创建的成像颜色板中进行成像显示，或由声信息处理单元3和声学扫描控制单元5实时获取扫描域Su中每个位置点的反射声学信息信号和位置坐标pij(xi,yj)信号，然后由声学成像单元6按照式(7、8、11)中的映射关系，转换为成像信号在声学成像单元6创建的成像颜色板中进行成像显示，成像颜色板的大小由成像域Sp按照式(12)确定，
Sp＝M×N (12)。
扫描域Su的确定方法是，
Su＝(2.0～3.3)So
这里，So为被观察牙体样本8在扫描平面内的最大投影面积，
—分别为被观察牙体样本8最大投影面积在x方向的终点和起点，
—分别为被观察牙体样本8最大投影面积在y方向的终点和起点，
使被观察的牙体样本8位于最佳的声学扫描成像区。
声学成像中的的确定方法是，
声学透镜1的焦点大小，根据被观察牙体样本的特征和观察分析要求，在10微米-100微米范围内。
实施例
选择中航复合材料有限责任公司生产的声显微扫描系统，选用10微米、50微米和100微米焦点声束的声学透镜，分别对来自临床的人体实际牙样本，牙样本在扫描平面内的最大投影面积S0约为10mm×12mm，将牙体样放置在本牙体样本固定平台中心位置，并使牙体样位于声学透镜的下方，取 分别取和采用反射法和透射法分别对牙样本进行了声显微成像试验分析和应用，应用结果表明，采用所构建的牙体组织声显微成像方法，可以快速清晰的得到牙体内部釉质、牙本质、牙髓、牙骨质等组织及其变化的特征分布，取得了较好的实际声显微成像效果。

Claims (4)

1.一种用于牙体的声显微成像方法，采用由声学透镜(1)、声学激励/接收单元(2)、声信息处理单元(3)、声学扫描机构(4)、声学扫描控制单元(5)、声学成像单元(6)、牙体样本固定平台(7)组成的声显微成像装置，其特征是，
(1)声学扫描
①将被观察的牙体样本(8)放置在声显微成像装置的牙体样本固定平台(7)中的声学扫描台容器内，并使牙体样本(8)中被观察面置于声学透镜(1)的下方；
②根据牙体样本(8)在扫描平面内的最大投影面积S0，确定声扫描域Su的大小，
<mrow> <msup> <mi>S</mi> <mi>u</mi> </msup> <mo>=</mo> <msubsup> <mi>x</mi> <mrow> <mi>m</mi> <mi>a</mi> <mi>x</mi> </mrow> <mi>u</mi> </msubsup> <mo>&amp;times;</mo> <msubsup> <mi>y</mi> <mrow> <mi>m</mi> <mi>a</mi> <mi>x</mi> </mrow> <mi>u</mi> </msubsup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow>
且使Su＞S0，
这里，分别为声透镜(1)在x、y方向的最大扫描范围，且
<mrow> <msubsup> <mi>x</mi> <mrow> <mi>m</mi> <mi>a</mi> <mi>x</mi> </mrow> <mi>u</mi> </msubsup> <mo>=</mo> <msubsup> <mi>x</mi> <mrow> <mi>e</mi> <mi>n</mi> <mi>d</mi> </mrow> <mi>u</mi> </msubsup> <mo>-</mo> <msubsup> <mi>x</mi> <mrow> <mi>s</mi> <mi>t</mi> </mrow> <mi>u</mi> </msubsup> </mrow>
<mrow> <msubsup> <mi>y</mi> <mrow> <mi>m</mi> <mi>a</mi> <mi>x</mi> </mrow> <mi>u</mi> </msubsup> <mo>=</mo> <msubsup> <mi>y</mi> <mrow> <mi>e</mi> <mi>n</mi> <mi>d</mi> </mrow> <mi>u</mi> </msubsup> <mo>-</mo> <msubsup> <mi>y</mi> <mrow> <mi>s</mi> <mi>t</mi> </mrow> <mi>u</mi> </msubsup> </mrow>
其中，—分别为声学透镜(1)在x方向的扫描终点和起点，
—分别为声学透镜(1)在y方向的扫描终点和起点，
③根据被观察牙体样本(8)的实际大小，根据式(1)，在声学扫描控制单元(5)中设置牙体样本(8)的(x,y)方向扫描范围；
④设置(x,y)方向声学扫描分辨率Δx、Δy，且使Δx＝Δy；
⑤设置声学透镜(1)的扫描速度υ，
扫描分辨率按式(2)确定，
<mrow> <mi>&amp;upsi;</mi> <mo>&amp;le;</mo> <mfrac> <mrow> <mi>&amp;Delta;</mi> <mi>x</mi> </mrow> <mrow> <mi>&amp;Delta;</mi> <mi>t</mi> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow>
式中，Δt—为每点声学信息和位置信息采集和成像显示所需时间；
⑥扫描成像，按照所设置和选择的声学扫描参数，对牙体样本(8)进行声学扫描；
(2)声学成像
利用声学透镜(1)提取扫描过程中第i行第j列位置点的反射声学信息或透射声学信息对牙体样本(8)进行声学成像,
1)反射法声学成像
①利用声学透镜(1)提取扫描过程中第i行第j列位置点的反射声学信息
<mrow> <msubsup> <mi>u</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> <mi>R</mi> </msubsup> <mrow> <mo>(</mo> <msubsup> <mi>t</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> <mi>R</mi> </msubsup> <mo>)</mo> </mrow> <mo>=</mo> <msubsup> <mi>Z</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> <mi>R</mi> </msubsup> <msub> <mi>A</mi> <mn>0</mn> </msub> <msup> <mi>e</mi> <mrow> <mo>-</mo> <msubsup> <mi>&amp;alpha;vt</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> <mi>R</mi> </msubsup> </mrow> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow>
利用反射声学信息阵列表示扫描域Su内的各个检测位置的反射声学信息：
<mrow> <msubsup> <mi>U</mi> <mrow> <mi>m</mi> <mo>&amp;times;</mo> <mi>n</mi> </mrow> <mi>R</mi> </msubsup> <mo>=</mo> <msub> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msubsup> <mi>u</mi> <mn>11</mn> <mi>R</mi> </msubsup> <mo>(</mo> <msubsup> <mi>t</mi> <mn>11</mn> <mi>R</mi> </msubsup> <mo>)</mo> <mo>,</mo> <mn>...</mn> <mo>,</mo> <msubsup> <mi>&amp;mu;</mi> <mrow> <mn>1</mn> <mi>j</mi> </mrow> <mi>R</mi> </msubsup> <mo>(</mo> <msubsup> <mi>t</mi> <mrow> <mn>1</mn> <mi>j</mi> </mrow> <mi>R</mi> </msubsup> <mo>)</mo> <mo>,</mo> <mn>..</mn> <mo>,</mo> <msubsup> <mi>&amp;mu;</mi> <mrow> <mn>1</mn> <mi>n</mi> </mrow> <mi>R</mi> </msubsup> <mo>(</mo> <msubsup> <mi>t</mi> <mrow> <mn>1</mn> <mi>n</mi> </mrow> <mi>R</mi> </msubsup> <mo>)</mo> </mtd> </mtr> <mtr> <mtd> <msubsup> <mi>u</mi> <mn>21</mn> <mi>R</mi> </msubsup> <mo>(</mo> <msubsup> <mi>t</mi> <mn>21</mn> <mi>R</mi> </msubsup> <mo>)</mo> <mo>,</mo> <mn>...</mn> <mo>,</mo> <msubsup> <mi>&amp;mu;</mi> <mrow> <mn>2</mn> <mi>j</mi> </mrow> <mi>R</mi> </msubsup> <mo>(</mo> <msubsup> <mi>t</mi> <mrow> <mn>2</mn> <mi>j</mi> </mrow> <mi>R</mi> </msubsup> <mo>)</mo> <mo>,</mo> <mn>..</mn> <mo>,</mo> <msubsup> <mi>&amp;mu;</mi> <mrow> <mn>2</mn> <mi>n</mi> </mrow> <mi>R</mi> </msubsup> <mo>(</mo> <msubsup> <mi>t</mi> <mrow> <mn>2</mn> <mi>n</mi> </mrow> <mi>R</mi> </msubsup> <mo>)</mo> </mtd> </mtr> <mtr> <mtd> <mn>....</mn> </mtd> </mtr> <mtr> <mtd> <msubsup> <mi>u</mi> <mrow> <mi>i</mi> <mn>1</mn> </mrow> <mi>R</mi> </msubsup> <mo>(</mo> <msubsup> <mi>t</mi> <mrow> <mi>i</mi> <mn>1</mn> </mrow> <mi>R</mi> </msubsup> <mo>)</mo> <mo>,</mo> <mn>...</mn> <mo>,</mo> <msubsup> <mi>&amp;mu;</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> <mi>R</mi> </msubsup> <mo>(</mo> <msubsup> <mi>t</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> <mi>R</mi> </msubsup> <mo>)</mo> <mo>,</mo> <mn>..</mn> <mo>,</mo> <msubsup> <mi>&amp;mu;</mi> <mrow> <mi>i</mi> <mi>n</mi> </mrow> <mi>R</mi> </msubsup> <mo>(</mo> <msubsup> <mi>t</mi> <mrow> <mi>i</mi> <mi>n</mi> </mrow> <mi>R</mi> </msubsup> <mo>)</mo> </mtd> </mtr> <mtr> <mtd> <mn>...</mn> </mtd> </mtr> <mtr> <mtd> <msubsup> <mi>u</mi> <mrow> <mi>m</mi> <mn>1</mn> </mrow> <mi>R</mi> </msubsup> <mo>(</mo> <msubsup> <mi>t</mi> <mrow> <mi>m</mi> <mn>1</mn> </mrow> <mi>R</mi> </msubsup> <mo>)</mo> <mo>,</mo> <mn>...</mn> <mo>,</mo> <msubsup> <mi>&amp;mu;</mi> <mrow> <mi>m</mi> <mi>j</mi> </mrow> <mi>R</mi> </msubsup> <mo>(</mo> <msubsup> <mi>t</mi> <mrow> <mi>m</mi> <mi>j</mi> </mrow> <mi>R</mi> </msubsup> <mo>)</mo> <mo>,</mo> <mn>..</mn> <mo>,</mo> <msubsup> <mi>&amp;mu;</mi> <mrow> <mi>m</mi> <mi>n</mi> </mrow> <mi>R</mi> </msubsup> <mo>(</mo> <msubsup> <mi>t</mi> <mrow> <mi>m</mi> <mi>n</mi> </mrow> <mi>R</mi> </msubsup> <mo>)</mo> </mtd> </mtr> </mtable> </mfenced> <mrow> <mi>m</mi> <mo>&amp;times;</mo> <mi>n</mi> </mrow> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>)</mo> </mrow> </mrow>
式中，
—对应牙体样本(8)中声学反射系数，与牙体组织和牙体均匀性有关，
A0—为牙体样本(8)中入射声波信号的幅值，
α—为入射声波在牙体样本(8)中的声衰减系数，
ν—为入射声波在牙体样本中的声速，
—为来自第i行第j列位置牙体组织中的反射声波的传播时间，
i、j—分别表示第i行第j列位置序列，i＝1,2,...,m,j＝1,2,...,n,
m、n—分别表示反射声学信息阵列的行数和列数；
②利用声学扫描控制单元(5)提取扫描域Su内第i行第j列位置的检测位置坐标pij(xi,yj)，利用位置阵列Pm×n表示扫描域Su内各个位置点的坐标位置：
<mrow> <msub> <mi>P</mi> <mrow> <mi>m</mi> <mo>&amp;times;</mo> <mi>n</mi> </mrow> </msub> <mo>=</mo> <msub> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mo>(</mo> <msub> <mi>x</mi> <mn>1</mn> </msub> <mo>,</mo> <msub> <mi>y</mi> <mn>1</mn> </msub> <mo>)</mo> <mo>,</mo> <mo>...</mo> <mo>,</mo> <mo>(</mo> <msub> <mi>x</mi> <mn>1</mn> </msub> <mo>,</mo> <msub> <mi>y</mi> <mi>j</mi> </msub> <mo>)</mo> <mo>,</mo> <mn>..</mn> <mo>,</mo> <mo>(</mo> <msub> <mi>x</mi> <mn>1</mn> </msub> <mo>,</mo> <msub> <mi>y</mi> <mi>n</mi> </msub> <mo>)</mo> </mtd> </mtr> <mtr> <mtd> <mo>(</mo> <msub> <mi>x</mi> <mn>2</mn> </msub> <mo>,</mo> <msub> <mi>y</mi> <mn>1</mn> </msub> <mo>)</mo> <mo>,</mo> <mo>...</mo> <mo>,</mo> <mo>(</mo> <msub> <mi>x</mi> <mn>2</mn> </msub> <mo>,</mo> <msub> <mi>y</mi> <mi>j</mi> </msub> <mo>)</mo> <mo>,</mo> <mn>..</mn> <mo>,</mo> <mo>(</mo> <msub> <mi>x</mi> <mn>2</mn> </msub> <mo>,</mo> <msub> <mi>y</mi> <mi>n</mi> </msub> <mo>)</mo> </mtd> </mtr> <mtr> <mtd> <mo>...</mo> </mtd> </mtr> <mtr> <mtd> <mo>(</mo> <msub> <mi>x</mi> <mi>i</mi> </msub> <mo>,</mo> <msub> <mi>y</mi> <mn>1</mn> </msub> <mo>)</mo> <mo>,</mo> <mo>...</mo> <mo>,</mo> <mo>(</mo> <msub> <mi>x</mi> <mi>i</mi> </msub> <mo>,</mo> <msub> <mi>y</mi> <mi>j</mi> </msub> <mo>)</mo> <mo>,</mo> <mn>..</mn> <mo>,</mo> <mo>(</mo> <msub> <mi>x</mi> <mi>i</mi> </msub> <mo>,</mo> <msub> <mi>y</mi> <mi>n</mi> </msub> <mo>)</mo> </mtd> </mtr> <mtr> <mtd> <mo>...</mo> </mtd> </mtr> <mtr> <mtd> <mo>(</mo> <msub> <mi>x</mi> <mi>m</mi> </msub> <mo>,</mo> <msub> <mi>y</mi> <mn>1</mn> </msub> <mo>)</mo> <mo>,</mo> <mo>...</mo> <mo>,</mo> <mo>(</mo> <msub> <mi>x</mi> <mi>m</mi> </msub> <mo>,</mo> <msub> <mi>y</mi> <mi>j</mi> </msub> <mo>)</mo> <mo>,</mo> <mn>..</mn> <mo>,</mo> <mo>(</mo> <msub> <mi>x</mi> <mi>m</mi> </msub> <mo>,</mo> <msub> <mi>y</mi> <mi>n</mi> </msub> <mo>)</mo> </mtd> </mtr> </mtable> </mfenced> <mrow> <mi>m</mi> <mo>&amp;times;</mo> <mi>n</mi> </mrow> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>5</mn> <mo>)</mo> </mrow> </mrow>
③利用声信息处理单元(3)对进行数字化后，按照式(6)转换为成像信号dkl(rkl,gkl,bkl)，此时：
<mrow> <mtable> <mtr> <mtd> <mrow> <msub> <mi>r</mi> <mrow> <mi>k</mi> <mi>l</mi> </mrow> </msub> <mo>=</mo> <msub> <mi>k</mi> <mi>r</mi> </msub> <msubsup> <mi>u</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> <mi>R</mi> </msubsup> <mrow> <mo>(</mo> <msubsup> <mi>t</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> <mi>R</mi> </msubsup> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>g</mi> <mrow> <mi>k</mi> <mi>l</mi> </mrow> </msub> <mo>=</mo> <msub> <mi>k</mi> <mi>g</mi> </msub> <msubsup> <mi>u</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> <mi>R</mi> </msubsup> <mrow> <mo>(</mo> <msubsup> <mi>t</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> <mi>R</mi> </msubsup> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>b</mi> <mrow> <mi>k</mi> <mi>l</mi> </mrow> </msub> <mo>=</mo> <msub> <mi>k</mi> <mi>b</mi> </msub> <msubsup> <mi>u</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> <mi>R</mi> </msubsup> <mrow> <mo>(</mo> <msubsup> <mi>t</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> <mi>R</mi> </msubsup> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>6</mn> <mo>)</mo> </mrow> </mrow>
这里，rkl,gkl,bkl—分别对应每个成像点的颜色值，
kr,kg,kb—分别为信号幅值与rkl,gkl,bkl颜色值之间的转换系数，
k、l—分别对应成像域中的位置坐标序列，且取
<mrow> <mtable> <mtr> <mtd> <mrow> <mi>k</mi> <mo>=</mo> <msub> <mi>&amp;Delta;n</mi> <mi>x</mi> </msub> <mo>&amp;times;</mo> <mi>i</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>l</mi> <mo>=</mo> <msub> <mi>&amp;Delta;n</mi> <mi>y</mi> </msub> <mo>&amp;times;</mo> <mi>j</mi> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>7</mn> <mo>)</mo> </mrow> </mrow>
式(7)中，Δnx、Δny—分别为扫描域Su中每个位置点在成像域Sp中的像素点数，此时若用Δp表示Su中每个位置点在成像域Sp中的像素大小，则有：
Δp＝Δnx×Δny×Δs
这里，Δs—为声学成像单元(6)中所用图像显示器的像素点的物理大小，
将扫描域Su中每个检测点位置坐标pij(xi,yj)映射到成像域Sp中位置坐标qkl(xk,yl)，并结合式(4)和式(5)，对每个点成像信号dkl(rkl,gkl,bkl)和成像位置坐标qkl(xk,yl)进行融合后，将对应扫描区域Su的成像结果在成像域Sp中通过矩阵QM×N表示，即有：
<mrow> <msub> <mi>Q</mi> <mrow> <mi>M</mi> <mo>&amp;times;</mo> <mi>N</mi> </mrow> </msub> <mo>=</mo> <msub> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mo>&amp;lsqb;</mo> <msub> <mi>x</mi> <mn>1</mn> </msub> <mo>,</mo> <msub> <mi>y</mi> <mn>1</mn> </msub> <mo>,</mo> <mo>(</mo> <msub> <mi>r</mi> <mn>11</mn> </msub> <mo>,</mo> <msub> <mi>g</mi> <mn>11</mn> </msub> <mo>,</mo> <msub> <mi>b</mi> <mn>11</mn> </msub> <mo>)</mo> <mo>&amp;rsqb;</mo> <mo>,</mo> <mn>...</mn> <mo>,</mo> <mo>&amp;lsqb;</mo> <msub> <mi>x</mi> <mn>1</mn> </msub> <mo>,</mo> <msub> <mi>y</mi> <mi>l</mi> </msub> <mo>,</mo> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mn>1</mn> <mi>l</mi> </mrow> </msub> <mo>,</mo> <msub> <mi>g</mi> <mrow> <mn>1</mn> <mi>l</mi> </mrow> </msub> <mo>,</mo> <msub> <mi>b</mi> <mrow> <mn>1</mn> <mi>l</mi> </mrow> </msub> <mo>)</mo> <mo>&amp;rsqb;</mo> <mo>,</mo> <mn>..</mn> <mo>,</mo> <mo>&amp;lsqb;</mo> <msub> <mi>x</mi> <mn>1</mn> </msub> <mo>,</mo> <msub> <mi>y</mi> <mi>N</mi> </msub> <mo>,</mo> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mn>1</mn> <mi>N</mi> </mrow> </msub> <mo>,</mo> <msub> <mi>g</mi> <mrow> <mn>1</mn> <mi>N</mi> </mrow> </msub> <mo>,</mo> <msub> <mi>b</mi> <mrow> <mn>1</mn> <mi>N</mi> </mrow> </msub> <mo>)</mo> <mo>&amp;rsqb;</mo> </mtd> </mtr> <mtr> <mtd> <mo>&amp;lsqb;</mo> <msub> <mi>x</mi> <mn>2</mn> </msub> <mo>,</mo> <msub> <mi>y</mi> <mn>1</mn> </msub> <mo>,</mo> <mo>(</mo> <msub> <mi>r</mi> <mn>21</mn> </msub> <mo>,</mo> <msub> <mi>g</mi> <mn>21</mn> </msub> <mo>,</mo> <msub> <mi>b</mi> <mn>21</mn> </msub> <mo>)</mo> <mo>&amp;rsqb;</mo> <mo>,</mo> <mn>...</mn> <mo>,</mo> <mo>&amp;lsqb;</mo> <msub> <mi>x</mi> <mn>2</mn> </msub> <mo>,</mo> <msub> <mi>y</mi> <mi>l</mi> </msub> <mo>,</mo> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mn>2</mn> <mi>l</mi> </mrow> </msub> <mo>,</mo> <msub> <mi>g</mi> <mrow> <mn>2</mn> <mi>l</mi> </mrow> </msub> <mo>,</mo> <msub> <mi>b</mi> <mrow> <mn>2</mn> <mi>l</mi> </mrow> </msub> <mo>)</mo> <mo>&amp;rsqb;</mo> <mo>,</mo> <mn>..</mn> <mo>,</mo> <mo>&amp;lsqb;</mo> <msub> <mi>x</mi> <mn>2</mn> </msub> <mo>,</mo> <msub> <mi>y</mi> <mi>N</mi> </msub> <mo>,</mo> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mn>2</mn> <mi>N</mi> </mrow> </msub> <mo>,</mo> <msub> <mi>g</mi> <mrow> <mn>2</mn> <mi>N</mi> </mrow> </msub> <mo>,</mo> <msub> <mi>b</mi> <mrow> <mn>2</mn> <mi>N</mi> </mrow> </msub> <mo>)</mo> <mo>&amp;rsqb;</mo> </mtd> </mtr> <mtr> <mtd> <mo>...</mo> </mtd> </mtr> <mtr> <mtd> <mo>&amp;lsqb;</mo> <msub> <mi>x</mi> <mi>k</mi> </msub> <mo>,</mo> <msub> <mi>y</mi> <mn>1</mn> </msub> <mo>,</mo> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>k</mi> <mn>1</mn> </mrow> </msub> <mo>,</mo> <msub> <mi>g</mi> <mrow> <mi>k</mi> <mn>1</mn> </mrow> </msub> <mo>,</mo> <msub> <mi>b</mi> <mrow> <mi>k</mi> <mn>1</mn> </mrow> </msub> <mo>)</mo> <mo>&amp;rsqb;</mo> <mo>,</mo> <mn>...</mn> <mo>,</mo> <mo>&amp;lsqb;</mo> <msub> <mi>x</mi> <mi>k</mi> </msub> <mo>,</mo> <msub> <mi>y</mi> <mi>l</mi> </msub> <mo>,</mo> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>k</mi> <mi>l</mi> </mrow> </msub> <mo>,</mo> <msub> <mi>g</mi> <mrow> <mi>k</mi> <mi>l</mi> </mrow> </msub> <mo>,</mo> <msub> <mi>b</mi> <mrow> <mi>k</mi> <mi>l</mi> </mrow> </msub> <mo>)</mo> <mo>&amp;rsqb;</mo> <mo>,</mo> <mn>..</mn> <mo>,</mo> <mo>&amp;lsqb;</mo> <msub> <mi>x</mi> <mi>k</mi> </msub> <mo>,</mo> <msub> <mi>y</mi> <mi>N</mi> </msub> <mo>,</mo> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>k</mi> <mi>N</mi> </mrow> </msub> <mo>,</mo> <msub> <mi>g</mi> <mrow> <mi>k</mi> <mi>N</mi> </mrow> </msub> <mo>,</mo> <msub> <mi>b</mi> <mrow> <mi>k</mi> <mi>N</mi> </mrow> </msub> <mo>)</mo> <mo>&amp;rsqb;</mo> </mtd> </mtr> <mtr> <mtd> <mo>...</mo> </mtd> </mtr> <mtr> <mtd> <mo>&amp;lsqb;</mo> <msub> <mi>x</mi> <mi>M</mi> </msub> <mo>,</mo> <msub> <mi>y</mi> <mn>1</mn> </msub> <mo>,</mo> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>M</mi> <mn>1</mn> </mrow> </msub> <mo>,</mo> <msub> <mi>g</mi> <mrow> <mi>M</mi> <mn>1</mn> </mrow> </msub> <mo>,</mo> <msub> <mi>b</mi> <mrow> <mi>M</mi> <mn>1</mn> </mrow> </msub> <mo>)</mo> <mo>&amp;rsqb;</mo> <mo>,</mo> <mn>...</mn> <mo>,</mo> <mo>&amp;lsqb;</mo> <msub> <mi>x</mi> <mi>M</mi> </msub> <mo>,</mo> <msub> <mi>y</mi> <mi>l</mi> </msub> <mo>,</mo> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>M</mi> <mi>l</mi> </mrow> </msub> <mo>,</mo> <msub> <mi>g</mi> <mrow> <mi>M</mi> <mi>l</mi> </mrow> </msub> <mo>,</mo> <msub> <mi>b</mi> <mrow> <mi>M</mi> <mi>l</mi> </mrow> </msub> <mo>)</mo> <mo>&amp;rsqb;</mo> <mo>,</mo> <mn>..</mn> <mo>,</mo> <mo>&amp;lsqb;</mo> <msub> <mi>x</mi> <mi>M</mi> </msub> <mo>,</mo> <msub> <mi>y</mi> <mi>N</mi> </msub> <mo>,</mo> <mo>(</mo> <msub> <mi>r</mi> <mrow> <mi>M</mi> <mi>N</mi> </mrow> </msub> <mo>,</mo> <msub> <mi>g</mi> <mrow> <mi>M</mi> <mi>N</mi> </mrow> </msub> <mo>,</mo> <msub> <mi>b</mi> <mrow> <mi>M</mi> <mi>N</mi> </mrow> </msub> <mo>)</mo> <mo>&amp;rsqb;</mo> </mtd> </mtr> </mtable> </mfenced> <mrow> <mi>M</mi> <mo>&amp;times;</mo> <mi>N</mi> </mrow> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>8</mn> <mo>)</mo> </mrow> </mrow>
式中(8)，M＝m×Δnx，N＝n×Δny
2)透射法声学成像
①利用声学透镜(1)提取扫描过程中第i行第j列位置点的透射声波信息
<mrow> <msubsup> <mi>u</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> <mi>T</mi> </msubsup> <mrow> <mo>(</mo> <msubsup> <mi>t</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> <mi>T</mi> </msubsup> <mo>)</mo> </mrow> <mo>=</mo> <msubsup> <mi>Z</mi> <mrow> <msub> <mi>ij</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> </msub> </mrow> <mi>T</mi> </msubsup> <msub> <mi>A</mi> <mn>0</mn> </msub> <msup> <mi>e</mi> <mrow> <mo>-</mo> <msubsup> <mi>&amp;alpha;vt</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> <mi>T</mi> </msubsup> </mrow> </msup> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>9</mn> <mo>)</mo> </mrow> </mrow>
利用透射声学信息阵列表示扫描域Su内的各个检测位置的透射声学信息：
<mrow> <msubsup> <mi>U</mi> <mrow> <mi>m</mi> <mo>&amp;times;</mo> <mi>n</mi> </mrow> <mi>T</mi> </msubsup> <mo>=</mo> <msub> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <msubsup> <mi>u</mi> <mn>11</mn> <mi>T</mi> </msubsup> <mo>(</mo> <msubsup> <mi>t</mi> <mn>11</mn> <mi>T</mi> </msubsup> <mo>)</mo> <mo>,</mo> <mo>...</mo> <mo>,</mo> <msubsup> <mi>&amp;mu;</mi> <mrow> <mn>1</mn> <mi>j</mi> </mrow> <mi>T</mi> </msubsup> <mo>(</mo> <msubsup> <mi>t</mi> <mrow> <mn>1</mn> <mi>j</mi> </mrow> <mi>T</mi> </msubsup> <mo>)</mo> <mo>,</mo> <mn>..</mn> <mo>,</mo> <msubsup> <mi>u</mi> <mrow> <mn>1</mn> <mi>n</mi> </mrow> <mi>T</mi> </msubsup> <mo>(</mo> <msubsup> <mi>t</mi> <mrow> <mn>1</mn> <mi>n</mi> </mrow> <mi>T</mi> </msubsup> <mo>)</mo> </mtd> </mtr> <mtr> <mtd> <msubsup> <mi>u</mi> <mn>21</mn> <mi>T</mi> </msubsup> <mo>(</mo> <msubsup> <mi>t</mi> <mn>21</mn> <mi>T</mi> </msubsup> <mo>)</mo> <mo>,</mo> <mo>...</mo> <mo>,</mo> <msubsup> <mi>&amp;mu;</mi> <mrow> <mn>2</mn> <mi>j</mi> </mrow> <mi>T</mi> </msubsup> <mo>(</mo> <msubsup> <mi>t</mi> <mrow> <mn>2</mn> <mi>j</mi> </mrow> <mi>T</mi> </msubsup> <mo>)</mo> <mo>,</mo> <mn>..</mn> <mo>,</mo> <msubsup> <mi>u</mi> <mrow> <mn>2</mn> <mi>n</mi> </mrow> <mi>T</mi> </msubsup> <mo>(</mo> <msubsup> <mi>t</mi> <mrow> <mn>2</mn> <mi>n</mi> </mrow> <mi>T</mi> </msubsup> <mo>)</mo> </mtd> </mtr> <mtr> <mtd> <mo>...</mo> </mtd> </mtr> <mtr> <mtd> <msubsup> <mi>u</mi> <mrow> <mi>i</mi> <mn>1</mn> </mrow> <mi>T</mi> </msubsup> <mo>(</mo> <msubsup> <mi>t</mi> <mrow> <mi>i</mi> <mn>1</mn> </mrow> <mi>T</mi> </msubsup> <mo>)</mo> <mo>,</mo> <mo>...</mo> <mo>,</mo> <msubsup> <mi>&amp;mu;</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> <mi>T</mi> </msubsup> <mo>(</mo> <msubsup> <mi>t</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> <mi>T</mi> </msubsup> <mo>)</mo> <mo>,</mo> <mn>..</mn> <mo>,</mo> <msubsup> <mi>u</mi> <mrow> <mi>i</mi> <mi>n</mi> </mrow> <mi>T</mi> </msubsup> <mo>(</mo> <msubsup> <mi>t</mi> <mrow> <mi>i</mi> <mi>n</mi> </mrow> <mi>T</mi> </msubsup> <mo>)</mo> </mtd> </mtr> <mtr> <mtd> <mo>...</mo> </mtd> </mtr> <mtr> <mtd> <msubsup> <mi>u</mi> <mrow> <mi>m</mi> <mn>1</mn> </mrow> <mi>T</mi> </msubsup> <mo>(</mo> <msubsup> <mi>t</mi> <mrow> <mi>m</mi> <mn>1</mn> </mrow> <mi>T</mi> </msubsup> <mo>)</mo> <mo>,</mo> <mo>...</mo> <mo>,</mo> <msubsup> <mi>&amp;mu;</mi> <mrow> <mi>m</mi> <mi>j</mi> </mrow> <mi>T</mi> </msubsup> <mo>(</mo> <msubsup> <mi>t</mi> <mrow> <mi>m</mi> <mi>j</mi> </mrow> <mi>T</mi> </msubsup> <mo>)</mo> <mo>,</mo> <mn>..</mn> <mo>,</mo> <msubsup> <mi>u</mi> <mrow> <mi>m</mi> <mi>n</mi> </mrow> <mi>T</mi> </msubsup> <mo>(</mo> <msubsup> <mi>t</mi> <mrow> <mi>m</mi> <mi>n</mi> </mrow> <mi>T</mi> </msubsup> <mo>)</mo> </mtd> </mtr> </mtable> </mfenced> <mrow> <mi>m</mi> <mo>&amp;times;</mo> <mi>n</mi> </mrow> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>10</mn> <mo>)</mo> </mrow> </mrow>
式中，
—对应牙体样本中声学透射系数，与牙体组织和牙体均匀性有关，
α—为入射声波在牙体样本中的声衰减系数，
ν—为入射声波在牙体样本中的声速，
—为来自第i行第j列位置牙体组织中的透射声波的传播时间，
i、j—分别表示第i行第j列位置序列，i＝1,2,...,m,j＝1,2,...,n,
m、n—分别表示透射声学信息阵列的行数和列数；
②利用声学扫描控制单元(5)提取扫描域Su内第i行第j列位置的检测位置坐标pij(xi,yj)，利用式(5)中的位置阵列Pm×n表示扫描域Su内各个位置点的坐标位置，
③利用声信息处理单元(3)对进行数字化后，按照式(11)转换为成像信号dkl(rkl,gkl,bkl)，此时：
<mrow> <mtable> <mtr> <mtd> <mrow> <msub> <mi>r</mi> <mrow> <mi>k</mi> <mi>l</mi> </mrow> </msub> <mo>=</mo> <msub> <mi>k</mi> <mi>r</mi> </msub> <msubsup> <mi>u</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> <mi>T</mi> </msubsup> <mrow> <mo>(</mo> <msubsup> <mi>t</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> <mi>T</mi> </msubsup> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>g</mi> <mrow> <mi>k</mi> <mi>l</mi> </mrow> </msub> <mo>=</mo> <msub> <mi>k</mi> <mi>g</mi> </msub> <msubsup> <mi>u</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> <mi>T</mi> </msubsup> <mrow> <mo>(</mo> <msubsup> <mi>t</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> <mi>T</mi> </msubsup> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>b</mi> <mrow> <mi>k</mi> <mi>l</mi> </mrow> </msub> <mo>=</mo> <msub> <mi>k</mi> <mi>b</mi> </msub> <msubsup> <mi>u</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> <mi>T</mi> </msubsup> <mrow> <mo>(</mo> <msubsup> <mi>t</mi> <mrow> <mi>i</mi> <mi>j</mi> </mrow> <mi>T</mi> </msubsup> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>11</mn> <mo>)</mo> </mrow> </mrow>
将扫描域Su中每个检测点位置坐标pij(xi,yj)映射到成像域Sp中位置坐标qkl(xk,yl)，并结合式(5)和式(10)，对每个点成像信号dkl(rkl,gkl,bkl)和成像位置坐标qkl(xk,yl)进行融合后，利用式(8)将对应扫描区域Su的成像结果在成像域Sp中通过矩阵QM×N表示。
(3)扫描成像显示
根据设定的扫描成像参数，声学透镜(1)在声学激励/接收单元(2)作用下，通过声学扫描机构(4)和声学扫描控制单元(5)使声学透镜(1)在扫描域Su内对牙体样本进行扫描，由声信息处理单元(3)和声学扫描控制单元(5)实时获取扫描域Su中每个位置点的反射声学信息信号和位置坐标pij(xi,yj)信号，然后由声学成像单元(6)按照式(6、7、8)中的映射关系，转换为成像信号在声学成像单元(6)中创建的成像颜色板中进行成像显示，或由声信息处理单元(3)和声学扫描控制单元(5)实时获取扫描域Su中每个位置点的反射声学信息信号和位置坐标pij(xi,yj)信号，然后由声学成像单元(6)按照式(7、8、11)中的映射关系，转换为成像信号在声学成像单元(6)中创建的成像颜色板中进行成像显示，成像颜色板的大小由成像域Sp按照式(12)确定，
Sp＝M×N (12)。
2.根据权利要求1所述的一种用于牙体的声显微成像方法，其特征是，所述的扫描域Su的确定方法是，
Su＝(2.0～3.3)So
这里，So为被观察牙体样本(8)在扫描平面内的最大投影面积，
—分别为被观察牙体样本(8)最大投影面积在x方向的终点和起点，
—分别为被观察牙体样本(8)最大投影面积在y方向的终点和起点，
使被观察的牙体样本(8)位于最佳的声学扫描成像区。
3.根据权利要求1所述的一种用于牙体的声显微成像方法，其特征是，所述声学成像中的的确定方法是，
<mrow> <msubsup> <mi>x</mi> <mrow> <mi>s</mi> <mi>t</mi> </mrow> <mi>u</mi> </msubsup> <mo>=</mo> <mrow> <mo>(</mo> <mn>0.2</mn> <mo>~</mo> <mn>0.4</mn> <mo>)</mo> </mrow> <msubsup> <mi>x</mi> <mrow> <mi>s</mi> <mi>t</mi> </mrow> <mi>t</mi> </msubsup> </mrow>
<mrow> <msubsup> <mi>x</mi> <mrow> <mi>e</mi> <mi>n</mi> <mi>d</mi> </mrow> <mi>u</mi> </msubsup> <mo>=</mo> <mrow> <mo>(</mo> <mn>1.4</mn> <mo>~</mo> <mn>1.8</mn> <mo>)</mo> </mrow> <msubsup> <mi>x</mi> <mrow> <mi>e</mi> <mi>n</mi> <mi>d</mi> </mrow> <mi>t</mi> </msubsup> </mrow>
<mrow> <mtable> <mtr> <mtd> <mrow> <msubsup> <mi>y</mi> <mrow> <mi>s</mi> <mi>t</mi> </mrow> <mi>u</mi> </msubsup> <mo>=</mo> <mrow> <mo>(</mo> <mn>0.2</mn> <mo>~</mo> <mn>0.4</mn> <mo>)</mo> </mrow> <msubsup> <mi>y</mi> <mrow> <mi>s</mi> <mi>t</mi> </mrow> <mi>t</mi> </msubsup> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msubsup> <mi>y</mi> <mrow> <mi>e</mi> <mi>n</mi> <mi>d</mi> </mrow> <mi>u</mi> </msubsup> <mo>=</mo> <mrow> <mo>(</mo> <mn>1.4</mn> <mo>~</mo> <mn>1.8</mn> <mo>)</mo> </mrow> <msubsup> <mi>y</mi> <mrow> <mi>e</mi> <mi>n</mi> <mi>d</mi> </mrow> <mi>t</mi> </msubsup> </mrow> </mtd> </mtr> </mtable> <mo>.</mo> </mrow>
4.根据权利要求1所述的一种用于牙体的声显微成像方法，其特征是，所述的声学透镜(1)的焦点大小，根据被观察牙体样本的特征和观察分析要求，在10微米-100微米范围内。
CN201710785542.8A 2017-09-04 2017-09-04 一种用于牙体的声显微成像方法 Pending CN107782791A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710785542.8A CN107782791A (zh) 2017-09-04 2017-09-04 一种用于牙体的声显微成像方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710785542.8A CN107782791A (zh) 2017-09-04 2017-09-04 一种用于牙体的声显微成像方法

Publications (1)

Publication Number Publication Date
CN107782791A true CN107782791A (zh) 2018-03-09

Family

ID=61437532

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710785542.8A Pending CN107782791A (zh) 2017-09-04 2017-09-04 一种用于牙体的声显微成像方法

Country Status (1)

Country Link
CN (1) CN107782791A (zh)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102608036A (zh) * 2012-03-20 2012-07-25 中北大学 基于声学透镜和传感器阵列的三维光声成像系统及方法
WO2012135977A1 (en) * 2011-04-08 2012-10-11 Carestream Health, Inc. Intra oral camera having a liquid lens for image stabilization
CN205538646U (zh) * 2016-01-29 2016-08-31 湖南工业大学 一种基于led的光声显微装置
KR20160106790A (ko) * 2015-03-02 2016-09-13 한국광기술원 피부 및 치아 검사 시스템 및 방법
CN106943120A (zh) * 2017-04-21 2017-07-14 厦门大学 一种光声显微镜及其监测微泡在生物组织中爆破的方法
CN107072530A (zh) * 2014-09-16 2017-08-18 卡尔斯特里姆保健公司 使用激光投影的牙齿表面成像设备

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2012135977A1 (en) * 2011-04-08 2012-10-11 Carestream Health, Inc. Intra oral camera having a liquid lens for image stabilization
CN102608036A (zh) * 2012-03-20 2012-07-25 中北大学 基于声学透镜和传感器阵列的三维光声成像系统及方法
CN107072530A (zh) * 2014-09-16 2017-08-18 卡尔斯特里姆保健公司 使用激光投影的牙齿表面成像设备
KR20160106790A (ko) * 2015-03-02 2016-09-13 한국광기술원 피부 및 치아 검사 시스템 및 방법
CN205538646U (zh) * 2016-01-29 2016-08-31 湖南工业大学 一种基于led的光声显微装置
CN106943120A (zh) * 2017-04-21 2017-07-14 厦门大学 一种光声显微镜及其监测微泡在生物组织中爆破的方法

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
BETTINA BASRANI: "《牙体牙髓放射影像学第二版》", 31 July 2016 *
BIN RAO ET AL.: "Photoacoustic microscopy of human teeth", 《PROC. OF SPIE》 *
刘松平等: "复合材料冲击损伤高分辨率超声成像检测与损伤行为分析", 《机械工程学报》 *
岳松龄: "《现代龋病学》", 30 April 2009 *
应烨军 等: "熔焊焊缝超声扫描成像检测", 《无损检测》 *
王亚非 等: "声显微成像", 《压电与声光》 *

Similar Documents

Publication Publication Date Title
CN101101277B (zh) 一种高清晰焊缝超声成像无损检测方法
Sherar et al. A 100 MHz B-scan ultrasound backscatter microscope
CN104883965B (zh) 被检体信息获得装置、显示方法、程序和处理装置
CN106680825B (zh) 一种声学阵列成像系统与方法
US20030101820A1 (en) Acoustic microscope
CN102944521B (zh) 非接触式光声和光学相干断层双成像装置及其检测方法
EP2514364B1 (en) Measurement system, and image forming method and program
US9786040B2 (en) Ultrasound diagnostic apparatus and ultrasound two-dimensional cross-section image generation method
CN104898123B (zh) 基于角域虚拟源的水浸超声合成孔径聚焦成像方法
CN101036162A (zh) 在显示的图像数据中保持一致的解剖视图的方法和系统
CN205758513U (zh) 一种皮肤疾病多模态成像检测系统
JP2000028589A (ja) ３次元超音波画像化装置
CN103040429A (zh) 一种用于口腔的光学影像检测装置和成像方法
CN104323762A (zh) 一种基于光声显微成像的鲜红斑痣血管定量化检测装置
CN105433906A (zh) 一种扫描暗场激光散斑血流成像方法及装置
Yao An ultrasonic method for 3D reconstruction of surface topography
CN106257323B (zh) 叠层成像术成像的方法
CN102512207A (zh) B超探头式光声成像乳腺扫描仪及其扫描方法
CN102496320A (zh) 一种基于ct体数据的实时超声图像模拟方法
CN107782791A (zh) 一种用于牙体的声显微成像方法
CN117731327A (zh) 三维剪切波弹性成像方法和超声成像系统
CN107874741A (zh) 生物组织光学检测探头
JP5354885B2 (ja) 超音波診断システム
EP4076208B1 (en) Systems and methods for assessing a placenta
CN107515250A (zh) 一种用于牙体的声显微成像装置

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
RJ01 Rejection of invention patent application after publication

Application publication date: 20180309

RJ01 Rejection of invention patent application after publication