CN102655815B - 测量系统和图像形成方法 - Google Patents
测量系统和图像形成方法 Download PDFInfo
- Publication number
- CN102655815B CN102655815B CN200980162927.2A CN200980162927A CN102655815B CN 102655815 B CN102655815 B CN 102655815B CN 200980162927 A CN200980162927 A CN 200980162927A CN 102655815 B CN102655815 B CN 102655815B
- Authority
- CN
- China
- Prior art keywords
- data
- sound wave
- signal
- telecommunication
- image
- 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
- 238000000034 method Methods 0.000 title claims description 47
- 238000005259 measurement Methods 0.000 title description 2
- 238000012545 processing Methods 0.000 claims description 37
- 238000006243 chemical reaction Methods 0.000 claims description 19
- 238000005286 illumination Methods 0.000 claims description 9
- 238000001514 detection method Methods 0.000 claims description 8
- 239000006096 absorbing agent Substances 0.000 abstract description 16
- 230000003287 optical effect Effects 0.000 abstract description 10
- 238000003384 imaging method Methods 0.000 abstract description 4
- 230000000875 corresponding effect Effects 0.000 description 13
- 238000010586 diagram Methods 0.000 description 7
- 241001269238 Data Species 0.000 description 5
- 239000000203 mixture Substances 0.000 description 5
- 230000005540 biological transmission Effects 0.000 description 4
- 238000009826 distribution Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 239000000463 material Substances 0.000 description 4
- 238000003825 pressing Methods 0.000 description 4
- 238000003325 tomography Methods 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 230000002596 correlated effect Effects 0.000 description 3
- IIZPXYDJLKNOIY-JXPKJXOSSA-N 1-palmitoyl-2-arachidonoyl-sn-glycero-3-phosphocholine Chemical compound CCCCCCCCCCCCCCCC(=O)OC[C@H](COP([O-])(=O)OCC[N+](C)(C)C)OC(=O)CCC\C=C/C\C=C/C\C=C/C\C=C/CCCCC IIZPXYDJLKNOIY-JXPKJXOSSA-N 0.000 description 2
- 238000010521 absorption reaction Methods 0.000 description 2
- 210000004027 cell Anatomy 0.000 description 2
- 239000003795 chemical substances by application Substances 0.000 description 2
- 230000000052 comparative effect Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 239000000499 gel Substances 0.000 description 2
- 150000002333 glycines Chemical class 0.000 description 2
- 239000000976 ink Substances 0.000 description 2
- 239000000787 lecithin Substances 0.000 description 2
- 229940067606 lecithin Drugs 0.000 description 2
- 235000010445 lecithin Nutrition 0.000 description 2
- 230000001629 suppression Effects 0.000 description 2
- 238000013519 translation Methods 0.000 description 2
- 238000002604 ultrasonography Methods 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 229920001817 Agar Polymers 0.000 description 1
- 239000008272 agar Substances 0.000 description 1
- 230000003321 amplification Effects 0.000 description 1
- 230000033228 biological regulation Effects 0.000 description 1
- 201000011510 cancer Diseases 0.000 description 1
- 210000005056 cell body Anatomy 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- JJTNWUKHGJGCNS-UHFFFAOYSA-N hex-2-ene Chemical compound CCCC=CC.CCCC=CC JJTNWUKHGJGCNS-UHFFFAOYSA-N 0.000 description 1
- 229940028435 intralipid Drugs 0.000 description 1
- 230000001678 irradiating effect Effects 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 238000012634 optical imaging Methods 0.000 description 1
- 238000005293 physical law Methods 0.000 description 1
- 229920000642 polymer Polymers 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 238000005728 strengthening Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/56—Details of data transmission or power supply, e.g. use of slip rings
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/0093—Detecting, measuring or recording by applying one single type of energy and measuring its conversion into another type of energy
- A61B5/0095—Detecting, measuring or recording by applying one single type of energy and measuring its conversion into another type of energy by applying light and detecting acoustic waves, i.e. photoacoustic measurements
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/1702—Systems in which incident light is modified in accordance with the properties of the material investigated with opto-acoustic detection, e.g. for gases or analysing solids
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/1702—Systems in which incident light is modified in accordance with the properties of the material investigated with opto-acoustic detection, e.g. for gases or analysing solids
- G01N2021/1706—Systems in which incident light is modified in accordance with the properties of the material investigated with opto-acoustic detection, e.g. for gases or analysing solids in solids
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Pathology (AREA)
- Medical Informatics (AREA)
- General Health & Medical Sciences (AREA)
- Surgery (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Animal Behavior & Ethology (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- Heart & Thoracic Surgery (AREA)
- Biomedical Technology (AREA)
- Chemical & Material Sciences (AREA)
- Optics & Photonics (AREA)
- Radiology & Medical Imaging (AREA)
- Analytical Chemistry (AREA)
- High Energy & Nuclear Physics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Biochemistry (AREA)
- Immunology (AREA)
- General Physics & Mathematics (AREA)
- Acoustics & Sound (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
在光声成像中,存在不能从实际存在的光吸收体的图像即实际图像区分伪像以及虽然事实上什么都不存在但确定存在光吸收体的可能性。通过利用设置在样本中的声波反射表面,基于在图像数据上观察时代表光吸收体的实际图像和反射图像以及相关的伪像的相应的体素数据之间的位置关系,区分或减少可归因于声波接收器的带域的伪像。
Description
技术领域
本发明涉及用于接收声波并生成图像数据的测量系统、用于从声波生成图像数据的图像形成方法、以及用于使图像形成方法被执行的程序。
背景技术
用于使从诸如激光器的光源照射到样本的光在样本内传播并获得关于样本的内部的信息的光学成像装置的研究主要在医疗领域中一直在积极进展。作为基于这些研究的光学成像技术中的一种,提出了光声断层成像(PAT)(参见专利文献1)。
在光声断层成像中,脉冲光首先从光源产生并被照射到样本。使照射的光在样本内传播和扩散,同时体组织吸收所述光的能量并从其产生声波(一般是超声波)。所产生的声波在位于关于样本的周围的多个点处被检测,并且,所获得的信号被数学分析并被处理成图像数据。以下,这种处理被称为成像。可通过成像获得样本内的初始压力产生分布或光学特性值分布特别是光能量吸收密度分布,并且,可获得诸如恶性肿瘤的位置的关于样本的内部的信息。
引文列表
专利文献
PTL1:美国专利No.5840023
发明内容
技术问题
有时会在通过成像获得的图像中出现伪像(artifact)。伪像意指看起来如同存在某物(虽然实际上它不存在)的图像。伪像也被称为幻影(ghost)。特别地,当视角受限时,例如,当获得声波信号的位置沿一定的方向位于一平面(plane)上而不是绕样本跨过360度时,可归因于声波接收器的带域(band)的伪像表现得显著。到目前为止,不能从实际存在的光吸收体的图像(即,从实际图像)区分伪像。
因此,本发明的目的是,区分或减少(抑制)可归因于声波接收器的带域的伪像。
问题的解决方案
根据本发明,提供一种测量系统,该测量系统包括:用于接收通过向样本照射光产生的声波并用于将接收的声波转换成电信号的声波接收器,以及,用于通过使用所述电信号生成图像数据的图像形成装置,其中,声波接收器至少接收从样本内的检测目标产生的声波中的直接到达声波接收器的直接波和从检测目标产生的声波中的被设置在样本中的声波反射表面反射的反射波,并且将接收的波转换成相应的电信号,所述图像形成装置包括:用于将从直接波转换的电信号和从反射波转换的电信号转换成体素(voxel)数据或像素数据由此生成图像数据的数据生成单元;用于确定当在所述图像数据上观察时从直接波转换的体素数据或像素数据的位置与从反射波转换的体素数据或像素数据的位置关于图像数据上的与声波反射表面对应的位置是否对称的确定单元;和用于执行处理以区分或减少被确定单元确定为非对称的体素数据或像素数据的图像的处理单元。
并且,根据本发明,提供一种测量系统,该测量系统包括:用于接收通过向样本照射光产生的声波并用于将接收的声波转换成电信号的声波接收器和用于通过使用所述电信号生成图像数据的图像形成装置,其中,声波接收器至少接收从样本内的检测目标产生的声波中的直接到达声波接收器的直接波和从检测目标产生的声波中的被设置在样本中的声波反射表面反射的反射波,并且将接收的波转换成相应的电信号,所述图像形成装置包括:用于将从直接波转换的电信号和从反射波转换的电信号转换成体素数据或像素数据由此生成图像数据的数据生成单元;用于在与声波反射表面对应的图像数据上的位置处折 叠(fold)图像数据由此生成折叠的数据的折叠数据生成单元;和用于叠加(superimpose)所述折叠的数据与被折叠之前的图像数据的叠加单元。
本发明的有利效果
根据本发明,可通过执行处理以区分或减少伪像,获得高度可靠的图像数据。
附图说明
图1是示出可应用本发明的测量系统的配置的一个例子的框图。
图2是用于解释在可应用本发明的测量系统中获得的实际图像、反射图像和它们的伪像之间的位置关系的一个例子的示意图。
图3是示出根据本发明的实施例1的图像形成装置的内部配置的一个例子的框图。
图4示出本发明的实施例1的处理流程。
图5是示出根据本发明的实施例2的图像形成装置的内部配置的一个例子的框图。
图6示出本发明的实施例2的处理流程。
图7示出通过应用本发明的例子获得的结果。
图8示出通过使用基于已知的技术的图像数据生成方法获得的结果。
图9示出图7和图8中沿Y=1.0cm处的截面的结果。
具体实施方式
以下将借助于附图详细描述本发明。注意,在本发明中,声波意指包含所谓的声波、超声波和光声波并且在诸如近红外线的光(电磁波)被照射到样本时在样本内产生的弹性波。
[实施例1]
参照附图描述本发明的实施例1。在本实施例中描述用于规定、区分或减少伪像(包含去除它的情况)的系统配置。图1是可应用本 发明的测量系统的框图。在利用光声断层成像的测量系统中,从光源2产生光2(脉冲光),并且,光2通过诸如透镜和反射镜的光学装置4被照射到样本3。样本3内部的光吸收体5(即,要被检测的对象)吸收光的能量并产生声波6。产生的声波6的一部分直接到达声波接收器8,并且,产生的声波的另一部分在被用作声波反射表面的声波反射板7反射之后到达声波接收器8。在以下的描述中,在不被反射的情况下直接到达声波接收器8的声波被称为直接波,并且,在被声波反射板7反射之后到达声波接收器8的声波被称为反射波。声波接收器8接收声波以转换成电信号,并且向电信号处理电路9输出所述电信号。电信号处理电路9执行电信号的放大、数字转换等,并且将数字转换之后的电信号(即,数字信号)输出到图像形成装置10。图像形成装置10通过使用所述数字信号生成图像数据,并且,生成的图像数据通过显示装置11被显示为图像。并且,声波接收器8通过扫描控制器12被机械扫描并构成为能够在各种位置处执行测量。
在本发明中,用作声波反射表面的声波反射板7被设置在样本3上。优选地,声波反射板7在其面向样本3的侧具有平坦表面,并具有与样本的声学阻抗大大不同的声学阻抗(例如,约3×106Kg/m2·s)。并且,声波反射板7优选由透明并允许光2充分地通过它的材料制成。这是由于,这种材料使得光能够从设置声波反射板7的侧被照射到样本。当设置用于从两侧按压样本的按压板时,可使用设置在声波接收器8的相对侧的按压板作为声波反射板,或者,可以在该按压板上设置声波反射表面。
光源1意图照射可被活体的成分之中的特定的成分(例如,血色素)吸收的特定波长的光。作为光源1,测量系统包含能够产生具有5纳秒~50纳秒的持续期的脉冲光的至少一个脉冲光源。虽然优选使用具有大的输出的激光器作为光源,但是,作为激光器的替代,也可使用发光二极管等。可以使用诸如固态激光器、气体激光器、染料激光器和半导体激光器的各种激光器中的适当的一个作为激光器。另外,可从设置声波接收器8的侧照射所述光,或者,可从声波接收器8的 相对侧照射所述光。作为替代方案,可从样本的两侧照射所述光。
光学装置4包括例如用于反射光的反射镜和用于会聚或扩展光以改变光束的形状的透镜。除了反射镜和透镜以外,也可使用光学波导作为光学部件。所述光学部件可以是任何类型的,只要它使得从光源1发射的光能够以希望的形式被照射到样本3。光2优选通过透镜发散(diverge),而扩展到一定的区域中。并且,光2照射到样本3的区域优选可在样本上移动。换句话说,测量系统优选被构成为使得从光源1产生的光可在样本上移动。光的这种可动性使得能够在更大的范围上照射光。更优选地,光2照射到样本3的区域(即,照射到样本的光)与声波接收器8同步地移动。可例如通过使用诸如可移动反射镜的方式或通过机械移动光源自身的方式,实施移动光被照射到样本的区域的方法。
声波接收器8包含分别接收声波并将声波转换成电信号的一个或更多个装置,并且,它由例如利用压电现象的转换器(transducer)、利用光学共振的转换器或利用电容的变化的转换器构成。声波接收器8为任意类型的,只要它可接收声波并将声波转换成电信号。通过一维或二维布置多个用于接收声波的装置,能够同时在多个位置处接收声波,以缩短接收时间,并减少样本的振动等的影响。作为替代方案,也可通过扫描一个装置获得与通过二维或一维布置装置获得的信号类似的信号。可在样本的整个表面上设置多个装置。另外,优选在声波接收器8与样本之间施加诸如用于实现声学匹配的凝胶的声学匹配材料。
下面,将参照图2详细说明通过分析处理在图像形成装置10中获得的图像数据。图像数据意指代表关于样本的内部的信息的数据,不管所述信息是二维的还是三维的。通过排列多个像素数据构成二维信息,并且,通过排列多个体素数据构成三维信息。通过使用基于例如时域或傅立叶域处理的图像数据生成方法(图像重构方法),对在多个位置处获得的声学信号执行分析处理,获得像素数据和体素数据。虽然以下描述三维图像数据,但是,本发明也可被应用于二维图像数 据。
图2示出比较由获得的图像数据代表的图像与样本和声波接收器的相应位置之间的位置关系的概念图。通过光激励的声波从光吸收体(即,产生声波的位置)以球状形式扩展。声波的一部分作为直接波直接到达声波接收器8。声波的另一部分被具有不同的声学阻抗的声波反射表面16(由声波反射板7提供)反射,并且,所反射的波到达声波接收器8。声波接收器8接收直接波和反射波并将它们转换成相应的电信号。通过对包含直接波和反射波两者的电信号的电信号执行分析处理,获得提供光吸收体的实际图像14的体素数据和提供其反射图像17的体素数据。实际图像14是通过直接波形成的图像,并且,反射图像17是通过反射波形成的图像。反射图像17在根据关于反射的物理定律的位置处出现。如果声波反射表面16的形状和位置是已知的,那么可以计算这样的位置。当声波反射表面16是平坦表面时,实际图像14与反射图像17之间的关系同于在反射镜中所看到的实际图像和反射图像之间的关系。更具体而言,在图2中,反射图像17出现的位置被定义为以与从声波反射表面16到实际图像14的距离相同的距离沿声波反射表面16的法线方向距离实际图像14的位置,即a1=a2成立的位置。并且,由于伪像可归因于声波接收器的带域,因此,不仅出现源自实际图像14的伪像15,而且出现以类似的方式源自反射图像17的伪像18。
现在描述可归因于声波接收器的带域并且是本发明要解决的问题的伪像的产生机制。在光声断层成像中,已知由光激励的声波由跨着宽的带域的频率分量构成。一般地,由于用于接收这些声波的声波接收器仅能够接收带域的一部分,因此,从声波接收器输出的电信号采取仅由部分频带构成的信号波形。此时,由于带域的不足,信号在它返回基线时超出(overshoot beyond)基线,由此产生鸣响(ringing)(振动)。这样的产生的可归因于声波接收器的带域的鸣响分量导致伪像。据推测,其原因在于,样本的结构基本上仅通过相互强化的来自光吸收器的信号分量被再现,但是,如果出现鸣响,那么鸣响分量 也相互强化,从而导致伪像。
导致伪像的信号振动(鸣响分量)出现于当从声波接收器观看时在时间上比纯正的(genuine)信号晚的位置处。因此,当电信号通过分析处理被转换成体积数据(volume data)时,如图2所示,伪像总是出现于比实际图像和反射图像中的每一个更远离声波接收器8的位置处。换句话说,从源自实际图像的伪像15到声波反射表面16的距离b1与从反射图像18到声波反射表面16的距离b2之间的关系由b1≠b2给出。因此,当在作为边界的声波反射表面16的两侧比较这些图像时,实际存在的光吸收体的图像出现于面对称的位置处,并且,伪像出现于非面对称的位置处。
下面,参照图3和图4描述可应用本发明的用于规定与区分或减少伪像的具体方法。图3是示出根据实施例1的图像形成装置的内部配置的框图,以及图4是用于区分伪像的处理的流程图。
首先,脉冲光被照射到样本(S1-1),并且,通过脉冲光激励的声波被声波接收器接收(S1-2)。从声波接收器输出的电信号在电信号处理电路9中被放大并经受模拟-数字转换,并然后被存储为数字信号。在这些步骤中,声波接收器至少在从向样本照射光到反射波到达声波接收器的时段中保持在可接收状态。并且,电信号处理电路9需要将在这样的时段中接收的声波存储为数字信号。换句话说,电信号处理电路9至少存储直接波和反射波的相应的数字信号。另外,在样本内的多个位置处获得存储于电信号处理电路9中的数字信号。图像形成装置10中的体积数据生成单元101将在多个位置处获得的数字信号转换成代表样本内的空间信息的相应的体素数据,由此生成图像数据(S1-3)。
然后,在对称确定单元102中,确定当在生成的图像数据上观察时代表出现的图像的体素数据的位置关于与声波反射表面对应的位置是否是面对称的(S1-4)。这里,由于从声波反射表面到声波接收器的距离是已知的,因此,当在生成的图像数据上观察时与声波反射表面对应的位置也是已知的。因此,能够确定当在生成的图像数据上观 察时代表出现的图像的体素数据的位置关于与声波反射表面对应的位置是否面对称。另一可设想的方法是,利用在声波反射表面处产生的声波。在许多的情况下,样本3和声波反射板7在光吸收特性上相互不同,并且,在样本与声波反射板之间的界面(声波反射表面)处产生声波。通过接收在声波反射表面处产生的声波并将接收的声波转换成体素数据,形成基于在声波反射表面处产生的声波的图像(以下,称为“声波反射表面的图像”)。因此,通过使用代表声波反射表面的图像的体素数据(即,基于在声波反射表面处产生的声波的体素数据)的位置作为基准,可以确定其它的图像是否是面对称的。通过上述的步骤,当出现的图像不面对称时,它们可被规定为伪像。
标签附加单元103将标签附加到代表非面对称并且被规定为伪像的图像的体素数据(S1-5)。
然后,处理单元104通过使用例如着色技术对附加有标签的代表图像的体素数据(即,代表非面对称并且被规定为伪像的图像的体素数据)执行区分处理(S1-6),使得相关的图像当在显示装置11上被显示时可被区分为伪像。除了诸如简单的着色的区分处理,可以执行减少(抑制)处理从而以相对浅的颜色显示被规定为伪像的体素数据,或者以与背景水平对应的颜色示出这些体素数据,使得这些体素数据不被显示为明显的图像。
最后,显示装置11显示例如体素数据的截面图像或3D呈现图像(S1-7)。通过上述的步骤,伪像可被规定与区分或减少。
[实施例2]
虽然伪像在实施例1中被规定与区分或减少,但是,关于用于减少伪像而无需在图像数据中规定它们的系统配置描述实施例2。
图5是示出根据实施例2的图像形成装置的内部配置的框图,以及图6是用于减少伪像的处理的流程图。在图2中,生成图像数据之前的步骤与实施例1中的那些相同。换句话说,直到体积数据生成单元101生成图像数据的S2-3的处理与实施例1中的直到S1-3的处理相同。
然后,折叠数据生成单元105折叠生成的图像数据中的与从声波检测器观察时离声波反射表面更远的那侧对应的那些体素数据,即,存在提供反射图像的体素数据的区域中的那些体素数据,由此生成折叠数据(S2-4)。
并且,叠加单元106执行使折叠数据与原始图像数据叠加的处理(S2-5)。在经受叠加处理之后的图像数据中,代表通过实际存在的光吸收体形成的实际图像的体素数据附加地与代表反射图像的体素数据叠加,由此,两个体素数据的信号相互强化。但是,不附加地叠加伪像的相应的体素数据,得到的信号的强度相对地弱。假定x、y和z的位置处的原始图像数据的体素值(强度)是A1(x,y,z)并且x、y和z的位置处的折叠数据的体素值(强度)是A2(x,y,z),那么通过使用下式(1)执行具体的叠加处理。
[数值式1]
在所有的体素上执行叠加处理。即使当通过简单的乘法执行叠加处理时,也可减少伪像。但是,实际图像的体素值(强度)通过乘法被取平方。因此,为了保持体素值的定量特性即线性,优选取平方值的根。通过叠加处理,当存在背景时减少伪像。当不存在背景时,由于伪像乘以零,因此伪像可被去除。因此,背景越低,则减少伪像的效果越大。另外,可以以其它的各种方式,诸如基于下式(2)的方式,执行叠加处理。
[数值式2]
A1(x,y,z)+A2(x,y,z) (2)
最后,通过叠加处理获得的数据通过显示装置11被显示为例如截面图像或3D呈现图像(S2-6)。
通过上述的步骤,可以减少伪像而无需规定伪像的步骤。
[实施例3]
关于图1中的声波反射板7具有弯曲表面的情况描述实施例3。在本实施例中,以与实施例1和实施例2中的方式类似的方式生成图像数据。但是,当声波反射板具有弯曲表面时,得到的反射图像根据弯曲表面的形状畸变。鉴于以上的点,从反射界面的形状计算在反射图像中导致的畸变,并且,执行代表相关的反射图像的体素数据的位置的坐标变换,如同声波是被平坦的反射板反射的。换句话说,代表反射图像的体素数据的位置被坐标变换为在声波反射表面是平坦表面时要获得的体素数据的位置。当事先已知所导致的畸变时,可准备用于执行坐标变换的式子或表,并且,可向各体素数据应用准备的式子或表。
然后,能够以与实施例1和实施例2中的方式类似的方式区分和减少伪像。
[实施例4]
关于在物理上不存在提供声波反射表面的声波反射板的情况描述实施例4。通过扩展在实施例3中使用的方法以使其适于具有弯曲表面的反射板,本发明可进一步被应用于不存在反射板的情况。当样本处于空气中并且样本的声学阻抗与空气的声学阻抗不同时,在样本与空气之间的界面处反射声波。样本与空气之间的界面的形状通过3D扫描装置被读取,并且,当如实施例3那样根据所读取的界面的形状考虑在反射图像中导致的畸变时,执行代表相关的反射图像的体素数据的位置的坐标变换,如同声波被平坦的反射板反射。并且,由于在样本与空气之间的界面处(即,声波反射表面处)产生声波,因此,在样本与空气之间的界面处产生的声波的体素数据本身可被利用。换句话说,通过利用在界面处产生的声波的体素数据的位置,可检查被导致的畸变而无需通过3D扫描装置读取界面的形状。
然后,能够以与实施例1和实施例2中的方式类似的方式区分或减少伪像。
另外,也可通过执行以下的处理实现本发明。更具体而言,用于实现上述的实施例1~4的功能的软件(程序)通过网络或各种存储介质中的一些被提供给系统或装置,并且,系统或装置中的计算机(或CPU、MPU等)读取并执行所述程序。
例子
描述被实施为通过使用本发明的实施例2而确认有利的效果的实验的结果。实验系统的细节如下。样本的基材是豆油(soy-oil)注入卵磷脂(Intralipid)和水的混合物。所述混合物通过使用琼脂被成形为长方体。通过混合豆油注入卵磷脂、水和黑色墨汁(India墨)并将混合物成形为球体而制备的光吸收体被放置在样本内。样本被置于空气中,并且,通过使用Nd:YAG激光器,具有1064nm的波长的纳秒量级的脉冲光以所述光被扩展为入射到样本的整个表面上的方式从一侧被重复照射到样本。并且,由甲基戊烯(methylpentene)聚合物制成的声波透射板在其被照射脉冲光的入射面的相对侧保持与样本的表面接触。2D阵列声波接收器安装有被插入在样本与2D阵列声波接收器之间的声波透射板。声学匹配凝胶被施加到样本与声波透射板之间的界面以及声波透射板与声波接收器之间的界面中的每一个。这里使用的2D阵列声波接收器由沿X方向和Y方向中的每一个以2mm的节距排列的总共324个器件(即,X方向18个器件和Y方向18个器件)构成。这些器件中的每一个具有1MHz±40%的频带。在本例子中,没有使用声波接收板,并且,由于脉冲光入射表面保持与具有不同的声学阻抗的空气接触,因此,脉冲光入射表面被视为声波反射表面。
图7示出通过使用本发明的实施例2获得的MIP(最大强度投影)结果。这里使用的叠加方法基于式(1)。并且,图8作为比较例示出通过使用与上述的实验系统相同的实验系统借助于没有应用本发明的已知的图像数据生成方法获得的图像数据的MIP结果。在图7和图8中的每一个中,横轴表示Y方向的距离,纵轴表示Z方向的距离。在Z=4.5cm附近出现的线是光吸收体的实际图像14,在Z=5.4cm附近出现的线是声波反射表面的图像19,并且,在Z=6.1cm附近出现的 线是光吸收体的反射图像17。并且,在图8中,在Z=4.7cm附近出现的线是实际图像的伪像15,并且,在Z=6.3cm附近出现的线是反射图像的伪像18。
从图7可以看出,与图8中的相比,伪像减少。
并且,图9示出图7和图8所示的结果中的每一个中Y=1.0cm处的截面的强度变化。在图9中,横轴表示Z方向的距离,纵轴表示强度。并且,虚线表示比较例的结果,实线表示例子的结果。图9中的那些结果分别基于在Z=4.5cm附近出现的峰值被归一化。从这些结果可以看出,在分别代表实际图像和反射图像的大的峰值后面出现的伪像的强度降低到背景水平。由此,可以有效地减少伪像。
附图标记列表
1 光源
2 光
3 样本
4 光学部件
5 光吸收体或初始压力分布
6 声波
7 声波反射板
8 声波接收器
9 电信号处理电路
10 图像形成装置
11 显示装置
12 扫描控制器
14 光吸收体的实际图像
15 实际图像的伪像
16 声波反射表面
17 光吸收体的反射图像
18 反射图像的伪像
19 声波反射表面的图像
101 体积数据生成单元
102 对称确定单元
103 标签附加单元
104 处理单元
105 折叠数据生成单元
106 叠加单元
Claims (12)
1.一种测量系统,包括用于接收通过向样本照射光而产生的声波并将所接收的声波转换成电信号的声波接收器和用于通过使用所述电信号来生成图像数据的图像形成装置,
其中,所述声波接收器接收从样本内的检测目标产生的声波中的直接到达声波接收器的直接波和从检测目标产生的声波中的被设置在样本上的声波反射表面反射的反射波,并且将接收的波转换成相应的电信号,
所述图像形成装置包括:
数据生成单元,用于将从所述直接波转换的电信号和从所述反射波转换的电信号转换成体素数据或像素数据,由此生成图像数据;
确定单元,用于确定当在所述图像数据上观察时从所述直接波转换的体素数据或像素数据的位置与从所述反射波转换的体素数据或像素数据的位置关于所述图像数据上的与声波反射表面对应的位置是否是对称的;和
处理单元,用于执行处理以区分被所述确定单元确定为非对称的体素数据或像素数据的图像。
2.一种测量系统,包括用于接收通过向样本照射光而产生的声波并将所接收的声波转换成电信号的声波接收器和用于通过使用所述电信号来生成图像数据的图像形成装置,
其中,所述声波接收器接收从样本内的检测目标产生的声波中的直接到达声波接收器的直接波和从检测目标产生的声波中的被设置在样本上的声波反射表面反射的反射波,并且将接收的波转换成相应的电信号,
所述图像形成装置包括:
数据生成单元,用于将从所述直接波转换的电信号和从所述反射波转换的电信号转换成体素数据或像素数据,由此生成图像数据;
确定单元,用于确定当在所述图像数据上观察时从所述直接波转换的体素数据或像素数据的位置与从所述反射波转换的体素数据或像素数据的位置关于所述图像数据上的与声波反射表面对应的位置是否是对称的;和
处理单元,用于执行处理以减少被所述确定单元确定为非对称的体素数据或像素数据的图像。
3.一种测量系统,包括用于接收通过向样本照射光产生的声波并将接收的声波转换成电信号的声波接收器和用于通过使用所述电信号生成图像数据的图像形成装置,
其中,所述声波接收器接收从样本内的检测目标产生的声波中的直接到达所述声波接收器的直接波和从检测目标产生的声波中的被设置在样本上的声波反射表面反射的反射波,并且将接收的波转换成相应的电信号,
所述图像形成装置包括:
数据生成单元,用于将从所述直接波转换的电信号和从所述反射波转换的电信号转换成体素数据或像素数据,由此生成图像数据;
折叠数据生成单元,用于在所述图像数据上的与所述声波反射表面对应的位置处折叠所述图像数据,由此生成折叠的数据;和
叠加单元,用于将所述折叠的数据与被折叠之前的所述图像数据进行叠加。
4.根据权利要求1至3中任一项所述的测量系统,还包括电信号处理电路,所述电信号处理电路用于将通过声波接收器转换的电信号转换成数字信号,
其中,所述数据生成单元将已被转换成数字信号的电信号转换成体素数据或像素数据。
5.根据权利要求1至3中任一项所述的测量系统,其中,当声波反射表面是弯曲表面时,所述数据生成单元生成在执行从反射波转换的体素数据或像素数据的位置到当声波反射表面是平坦表面时获得的体素数据或像素数据的位置的坐标变换之后的图像数据。
6.根据权利要求4所述的测量系统,其中,当声波反射表面是弯曲表面时,所述数据生成单元生成在执行从反射波转换的体素数据或像素数据的位置到当声波反射表面是平坦表面时获得的体素数据或像素数据的位置的坐标变换之后的图像数据。
7.一种图像形成方法,用于通过声波接收器接收通过向样本照射光而产生的声波并将所接收的声波转换成电信号,以及用于通过使用所述电信号生成图像数据,所述图像形成方法包括:
数据生成步骤,将由从样本内的检测目标产生的声波中的直接到达所述声波接收器的直接波转换得到的电信号和由从检测目标产生的声波中的通过设置在样本上的声波反射表面反射的反射波转换得到的电信号转换成体素数据或像素数据,由此生成图像数据;
确定步骤,确定当在所述图像数据上观察时从所述直接波转换的体素数据或像素数据的位置与从所述反射波转换的体素数据或像素数据的位置关于所述图像数据上的与所述声波反射表面对应的位置是否是对称的;和
处理步骤,执行处理以区分在所述确定步骤中确定为非对称的体素数据或像素数据的图像。
8.一种图像形成方法,用于通过声波接收器接收通过向样本照射光而产生的声波并将所接收的声波转换成电信号,以及用于通过使用所述电信号生成图像数据,所述图像形成方法包括:
数据生成步骤,将由从样本内的检测目标产生的声波中的直接到达所述声波接收器的直接波转换得到的电信号和由从检测目标产生的声波中的通过设置在样本上的声波反射表面反射的反射波转换得到的电信号转换成体素数据或像素数据,由此生成图像数据;
确定步骤,确定当在所述图像数据上观察时从所述直接波转换的体素数据或像素数据的位置与从所述反射波转换的体素数据或像素数据的位置关于所述图像数据上的与所述声波反射表面对应的位置是否是对称的;和
处理步骤,执行处理以减少在所述确定步骤中确定为非对称的体素数据或像素数据的图像。
9.一种图像形成方法,用于通过声波接收器接收通过向样本照射光而产生的声波并将所接收的声波转换成电信号,以及用于通过使用所述电信号生成图像数据,所述图像形成方法包括:
数据生成步骤,将由从样本内的检测目标产生的声波中的直接到达所述声波接收器的直接波转换得到的电信号和由从检测目标产生的声波中的通过设置在样本上的声波反射表面反射的反射波转换得到的电信号转换成体素数据或像素数据,由此生成图像数据;
折叠步骤,在所述图像数据上的与声波反射表面对应的位置处折叠所述图像数据,由此生成折叠的数据;和
叠加步骤,将所述折叠的数据与被折叠之前的所述图像数据进行叠加。
10.根据权利要求7至9中任一项所述的图像形成方法,
其中,通过所述声波接收器转换的电信号被转换成数字信号,并且,
已被转换成数字信号的电信号在数据生成步骤中被转换成体素数据或像素数据。
11.根据权利要求7至9中任一项所述的图像形成方法,其中,当声波反射表面是弯曲表面时,在执行从反射波转换的体素数据或像素数据的位置到当声波反射表面是平坦表面时获得的体素数据或像素数据的位置的坐标变换之后在数据生成步骤中生成图像数据。
12.根据权利要求10所述的图像形成方法,其中,当声波反射表面是弯曲表面时,在执行从反射波转换的体素数据或像素数据的位置到当声波反射表面是平坦表面时获得的体素数据或像素数据的位置的坐标变换之后在数据生成步骤中生成图像数据。
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/JP2009/071062 WO2011074102A1 (ja) | 2009-12-17 | 2009-12-17 | 測定システム、画像形成方法及びプログラム |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102655815A CN102655815A (zh) | 2012-09-05 |
CN102655815B true CN102655815B (zh) | 2015-07-01 |
Family
ID=44150854
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN200980162927.2A Expired - Fee Related CN102655815B (zh) | 2009-12-17 | 2009-12-17 | 测量系统和图像形成方法 |
Country Status (5)
Country | Link |
---|---|
US (1) | US8814794B2 (zh) |
EP (1) | EP2514364B1 (zh) |
JP (1) | JP5518096B2 (zh) |
CN (1) | CN102655815B (zh) |
WO (1) | WO2011074102A1 (zh) |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7750536B2 (en) | 2006-03-02 | 2010-07-06 | Visualsonics Inc. | High frequency ultrasonic transducer and matching layer comprising cyanoacrylate |
JP5441781B2 (ja) * | 2010-03-25 | 2014-03-12 | キヤノン株式会社 | 光音響イメージング装置、光音響イメージング方法及びプログラム |
JP2014046069A (ja) * | 2012-09-03 | 2014-03-17 | Canon Inc | 被検体情報取得装置 |
JP6362301B2 (ja) * | 2013-04-30 | 2018-07-25 | キヤノン株式会社 | 被検体情報取得装置、および、被検体情報取得装置の作動方法 |
EP3116405A4 (en) | 2014-03-12 | 2018-01-03 | Fujifilm Sonosite, Inc. | High frequency ultrasound transducer having an ultrasonic lens with integral central matching layer |
US9937422B2 (en) | 2015-12-09 | 2018-04-10 | Microsoft Technology Licensing, Llc | Voxel-based, real-time acoustic adjustment |
US10293259B2 (en) | 2015-12-09 | 2019-05-21 | Microsoft Technology Licensing, Llc | Control of audio effects using volumetric data |
US10045144B2 (en) | 2015-12-09 | 2018-08-07 | Microsoft Technology Licensing, Llc | Redirecting audio output |
US10345267B2 (en) * | 2015-12-21 | 2019-07-09 | The Boeing Company | Composite inspection |
US11644442B2 (en) | 2016-04-01 | 2023-05-09 | The Board Of Regents Of The University Of Oklahoma | System and method for nanoscale photoacoustic tomography |
US10438382B2 (en) * | 2017-03-27 | 2019-10-08 | Canon Kabushiki Kaisha | Image processing apparatus and image processing method |
JP6882108B2 (ja) * | 2017-07-13 | 2021-06-02 | キヤノン株式会社 | 画像生成装置、画像生成方法、及びプログラム |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2001507952A (ja) * | 1996-10-04 | 2001-06-19 | オプトソニクス,インコーポレイテッド | 光音響胸部スキャナ |
CN1575770A (zh) * | 2003-07-02 | 2005-02-09 | 株式会社东芝 | 显示患者有关信息之影像的形成方法和仪器 |
CN1650794A (zh) * | 2004-02-06 | 2005-08-10 | 株式会社东芝 | 非探入式受检体信息成象方法和装置 |
CN101563035A (zh) * | 2006-12-19 | 2009-10-21 | 皇家飞利浦电子股份有限公司 | 组合光声和超声成像系统 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5840023A (en) | 1996-01-31 | 1998-11-24 | Oraevsky; Alexander A. | Optoacoustic imaging for medical diagnosis |
US20060116579A1 (en) * | 2004-11-29 | 2006-06-01 | Pai-Chi Li | Ultrasound imaging apparatus and method thereof |
JP4422626B2 (ja) * | 2005-01-25 | 2010-02-24 | 日本電信電話株式会社 | 生体画像化装置 |
GB2444078A (en) * | 2006-11-24 | 2008-05-28 | Ucl Business Plc | Ultrasonic sensor which detects direct and reflected signals emanating from sample being imaged |
JP5424846B2 (ja) * | 2009-12-11 | 2014-02-26 | キヤノン株式会社 | 光音響イメージング装置 |
-
2009
- 2009-12-17 CN CN200980162927.2A patent/CN102655815B/zh not_active Expired - Fee Related
- 2009-12-17 EP EP09852294.9A patent/EP2514364B1/en active Active
- 2009-12-17 WO PCT/JP2009/071062 patent/WO2011074102A1/ja active Application Filing
- 2009-12-17 JP JP2011545904A patent/JP5518096B2/ja not_active Expired - Fee Related
-
2010
- 2010-12-09 US US12/964,654 patent/US8814794B2/en not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2001507952A (ja) * | 1996-10-04 | 2001-06-19 | オプトソニクス,インコーポレイテッド | 光音響胸部スキャナ |
CN1575770A (zh) * | 2003-07-02 | 2005-02-09 | 株式会社东芝 | 显示患者有关信息之影像的形成方法和仪器 |
CN1650794A (zh) * | 2004-02-06 | 2005-08-10 | 株式会社东芝 | 非探入式受检体信息成象方法和装置 |
CN101563035A (zh) * | 2006-12-19 | 2009-10-21 | 皇家飞利浦电子股份有限公司 | 组合光声和超声成像系统 |
Non-Patent Citations (1)
Title |
---|
超声换能器带宽对光声成像的影响;谭毅等;《光学学报》;20050131;第25卷(第1期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
EP2514364A4 (en) | 2017-06-14 |
JP5518096B2 (ja) | 2014-06-11 |
US8814794B2 (en) | 2014-08-26 |
EP2514364A1 (en) | 2012-10-24 |
CN102655815A (zh) | 2012-09-05 |
JPWO2011074102A1 (ja) | 2013-04-25 |
EP2514364B1 (en) | 2020-02-12 |
WO2011074102A1 (ja) | 2011-06-23 |
US20110149680A1 (en) | 2011-06-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102655815B (zh) | 测量系统和图像形成方法 | |
CN102843960B (zh) | 被检体信息获得设备和方法 | |
CN103260502B (zh) | 被检体信息获取设备和被检体信息获取方法 | |
CN102525550B (zh) | 测量装置 | |
CN103025248B (zh) | 图像信息获取装置和图像信息获取方法 | |
US20200056986A1 (en) | Multifocal photoacoustic microscopy through an ergodic relay | |
US20160228090A1 (en) | Photoacoustic tracking and registration in interventional ultrasound | |
US20170172419A1 (en) | Photoacoustic apparatus and signal processing method | |
JP6362301B2 (ja) | 被検体情報取得装置、および、被検体情報取得装置の作動方法 | |
CN115177217B (zh) | 基于球形粒子光脉冲激发效应的光声信号仿真方法、装置 | |
CN104706323A (zh) | 一种高速大视场多光谱光声成像方法及其装置 | |
CN103389273A (zh) | 一种光声和光学融合的多模态成像系统 | |
WO2012035727A1 (en) | Photoacoustic imaging apparatus and control method thereof | |
CN108606777A (zh) | 基于可调聚焦型光纤传感器的光声计算层析成像系统 | |
CN103845078A (zh) | 对象信息获取设备和控制对象信息获取设备的方法 | |
CN201207035Y (zh) | 一种基于多元相控聚焦环阵列的光声成像装置 | |
Hakakzadeh et al. | Unipolar back-projection algorithm for photoacoustic tomography | |
US8611179B2 (en) | Acoustic wave measuring system and method of generating image data thereof | |
JP6184146B2 (ja) | 被検体情報取得装置およびその制御方法 | |
CN110907451A (zh) | 一种基于艾里光束和组织透明化技术的超深oct系统装置 | |
CN103654725A (zh) | 被检体信息获取设备 | |
JP2014188066A (ja) | 被検体情報取得装置およびその制御方法 | |
Gratt et al. | Photoacoustic imaging using a conical axicon detector | |
Webb et al. | Optical diffusion imaging using nonlinear multigrid optimization | |
CN105246396A (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150701 Termination date: 20161217 |