CN104013438A - 图像生成装置、图像生成方法和程序 - Google Patents

图像生成装置、图像生成方法和程序 Download PDF

Info

Publication number
CN104013438A
CN104013438A CN201410266998.XA CN201410266998A CN104013438A CN 104013438 A CN104013438 A CN 104013438A CN 201410266998 A CN201410266998 A CN 201410266998A CN 104013438 A CN104013438 A CN 104013438A
Authority
CN
China
Prior art keywords
supposition
target area
spread speed
view data
determining unit
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
Application number
CN201410266998.XA
Other languages
English (en)
Other versions
CN104013438B (zh
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.)
Canon Inc
Original Assignee
Canon Inc
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 Canon Inc filed Critical Canon Inc
Publication of CN104013438A publication Critical patent/CN104013438A/zh
Application granted granted Critical
Publication of CN104013438B publication Critical patent/CN104013438B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • 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
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/0059—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence
    • A61B5/0073—Measuring for diagnostic purposes; Identification of persons using light, e.g. diagnosis by transillumination, diascopy, fluorescence by tomography, i.e. reconstruction of 3D images from 2D projections
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/43—Detecting, measuring or recording for evaluating the reproductive systems
    • A61B5/4306—Detecting, measuring or recording for evaluating the reproductive systems for evaluating the female reproductive systems, e.g. gynaecological evaluations
    • A61B5/4312—Breast evaluation or disorder diagnosis
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7225—Details of analog processing, e.g. isolation amplifier, gain or sensitivity adjustment, filtering, baseline or drift compensation
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/46—Ultrasonic, sonic or infrasonic diagnostic devices with special arrangements for interfacing with the operator or the patient
    • A61B8/467—Ultrasonic, sonic or infrasonic diagnostic devices with special arrangements for interfacing with the operator or the patient characterised by special input means
    • A61B8/469—Ultrasonic, sonic or infrasonic diagnostic devices with special arrangements for interfacing with the operator or the patient characterised by special input means for selection of a region of interest
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5207—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of raw data to produce diagnostic data, e.g. for generating an image
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5215—Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/88—Sonar systems specially adapted for specific applications
    • G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
    • G01S15/8906—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/88—Sonar systems specially adapted for specific applications
    • G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
    • G01S15/8906—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
    • G01S15/8965—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using acousto-optical or acousto-electronic conversion techniques
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
    • G01S7/52046—Techniques for image enhancement involving transmitter or receiver
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/523—Details of pulse systems
    • G01S7/526—Receivers
    • G01S7/527—Extracting wanted echo signals
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/20—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for computer-aided diagnosis, e.g. based on medical expert systems

Landscapes

  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Public Health (AREA)
  • Medical Informatics (AREA)
  • Biomedical Technology (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Pathology (AREA)
  • General Health & Medical Sciences (AREA)
  • Animal Behavior & Ethology (AREA)
  • Veterinary Medicine (AREA)
  • Surgery (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Biophysics (AREA)
  • Molecular Biology (AREA)
  • Acoustics & Sound (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Signal Processing (AREA)
  • Gynecology & Obstetrics (AREA)
  • Reproductive Health (AREA)
  • Power Engineering (AREA)
  • Physiology (AREA)
  • Psychiatry (AREA)
  • Artificial Intelligence (AREA)
  • Data Mining & Analysis (AREA)
  • Databases & Information Systems (AREA)
  • Epidemiology (AREA)
  • Primary Health Care (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)
  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)

Abstract

本发明涉及图像生成装置、图像生成方法和程序。本发明的图像生成装置具有确定单元，该确定单元在被检体内部的一部分区域中设定目标区域，对于多个假定的速度，执行基于从所述检测元件中的每一个到目标区域的距离和假定的速度调整每一个检测信号的相位、以及计算经相位调整的多个检测信号的强度的离散度的处理，并且，确定所述多个假定的速度之中的强度的离散度最小化的速度作为传播速度。

Description

图像生成装置、图像生成方法和程序
本申请是申请日为2010年12月8日、申请号为201080055219.1、发明名称为“图像生成装置、图像生成方法和程序”的发明专利申请的分案申请。
技术领域
本发明涉及用于生成代表被检体内部的信息的图像的图像生成装置、图像生成方法和程序。
背景技术
在医学领域，对于将来自光源(例如，激光器)的光照射到诸如活体的被检体上并且将被检体内部的信息可视化的成像装置的研究一直在积极地进展。这种利用光的可视化技术的例子是光声断层成像(PAT)。光声断层成像装置在被检体周围的多个位置处检测从吸收了在被检体中传播并从中扩散的光的能量的活体的组织生成的声波(典型地为超声波)。然后，对获得的信号进行数学分析，并且，与被检体内部的光学性能有关的信息、特别是吸收系数分布被可视化。近来，正积极进行关于通过使用光声断层成像装置将小动物的血管成像的临床前研究和关于应用光声断层成像装置的原理以诊断乳癌等的临床研究。
在光声断层成像装置和常规上在医学领域中使用的超声诊断装置(用于检测在活体中反射的声波并生成图像的装置)的情况下，通常通过使用被检体的平均的声速(被检体内部的声波的声速、被检体内部的声波的传播速度或在被检体内部声波的传播速度)生成图像。一般地，基于实验值或基于文献的值确定声速。但是，传播速度具有个体差异，并且，例如，声速也根据保持被检体的方法改变。因此，如果用于生成图像的声速与实际的声速不同，那么图像的分辨率大大下降。
例如专利文献(PTL1)公开了一种用于解决该问题的方法。根据在专利文献(PTL1)中公开的技术，声速被确定为使得各像素(或体素(voxel))的亮度或对比度被最大化。由此，由于用于生成图像的声速与实际的声速的不匹配导致的图像质量的下降被抑制。
但是，在专利文献(PTL1)中的技术的情况下，由于各像素的亮度或对比度被最大化，因此，背景噪声的亮度或对比度也增大。并且，如果在检测信号中包含噪声，那么声速被确定为使得噪声分量和正常信号分量的总值被最大化，因此，不能获得精确的声速，并且，图像模糊。
(PTL1)日本专利申请公开No.2000-166925
发明内容
本发明提供可在抑制由于噪声和用于生成图像的声速与实际传播速度(声速)的不匹配导致的图像质量的下降的同时生成代表被检体内部的信息的图像的图像生成装置和图像生成方法。
本发明在其第一方面中提供一种图像生成装置，该图像生成装置包括：
具有检测从被检体内部传播的声波并输出检测信号的多个检测元件的探测器；
确定在被检体内部声波的传播速度的确定单元；
通过使用由确定单元确定的传播速度和从多个检测元件获得的多个检测信号生成代表被检体内部的信息的图像的图像生成单元，其中，
确定单元在被检体内部的一部分区域中设定目标区域，
对于多个假定的(tentative)速度，确定单元执行基于从所述检测元件中的每一个到目标区域的距离和假定的速度调整所述检测信号中的每一个的相位、以及计算经相位调整的多个检测信号的强度的离散度的处理，并且，
确定单元确定所述多个假定的速度之中的强度的离散度最小化的速度作为传播速度。
本发明在其第二方面中提供一种图像生成方法，该图像生成方法包括：
通过使用多个检测元件检测从被检体内部传播的声波并生成检测信号的步骤；
确定在被检体内部声波的传播速度的确定步骤；和
通过使用在所述确定步骤中确定的传播速度和从多个检测元件获得的多个检测信号生成代表被检体内部的信息的图像的图像生成步骤，其中，在所述确定步骤中，
在被检体内部的一部分区域中设定目标区域，
对于多个假定的速度，执行基于从所述检测元件中的每一个到目标区域的距离和假定的速度调整所述检测信号中的每一个的相位、以及计算经相位调整的多个检测信号的强度的离散度的处理，并且，
确定所述多个假定的速度之中的强度的离散度最小化的速度作为传播速度。
本发明在其第三方面中提供一种记录计算机程序的非暂态(non-transitory)计算机可读介质，该计算机程序用于使计算机执行包括以下的步骤的方法：
确定在被检体内部声波的传播速度的确定步骤；和
通过使用借助于用多个检测元件检测从被检体内部传播的声波生成的多个检测信号和在所述确定步骤中确定的传播速度来生成代表被检体内部的信息的图像的图像生成步骤，其中，
在所述确定步骤中，
在被检体内部的一部分区域中设定目标区域，
对于多个假定的速度，执行基于从所述检测元件中的每一个到目标区域的距离和假定的速度调整所述检测信号中的每一个的相位、以及计算经相位调整的多个检测信号的强度的离散度的处理，并且，
确定所述多个假定的速度之中的强度的离散度最小化的速度作为传播速度。
根据本发明，可在抑制由于噪声和用于生成图像的速度(声速)与实际传播速度(声速)的不匹配导致的图像质量的下降的同时，生成代表被检体内部的信息的图像。
从参照附图对示例性实施例的以下描述，本发明的其它特征将变得清晰。
附图说明
【图1】是示出本发明的原理的示图；
【图2】是示出根据本实施例的图像生成装置的配置的例子的示图；
【图3】是示出用于生成所生成的图像的方法的流程图；
【图4A】是示出检测元件与处理目标位置的位置关系的例子的示图；
【图4B】是示出检测元件的方向性的例子的示图；
【图5A】是示出根据例子1的模拟条件的例子的曲线图；
【图5B】是示出根据例子1的模拟的检测信号的例子的曲线图；
【图6】表示在根据例子1的模拟中生成的图像的例子；
【图7A】表示通过根据例子2的图像生成装置所生成的图像的例子；以及
【图7B】表示通过根据例子2的图像生成装置所生成的图像的例子。
具体实施方式
<原理>
现在将描述本发明的原理。
图1是示出本发明的原理的示图。在图1中，附图标记1表示声源(被检体内部)，附图标记2表示体素(或像素)，并且，附图标记3表示检测元件。各检测元件通过时间分割测量(time sharingmeasurement)检测从被检体内部传播的声波，并且输出检测信号。检测元件3的数量为N(N是2或更大的整数)。
一般地，使用声波(超声波)的图像生成装置通过使用从多个检测元件i获得的多个检测信号S(i,t)生成代表被检体内部的信息的图像。这种类型的图像被称为“生成图像”，并且，如果生成图像是三维图像(体素的集合)，则这种类型的图像也被称为“体数据(volumedata)”。如果生成图像是二维图像，则它也被称为“像素数据”。一般通过使用基于从各检测元件到与像素对应的位置的距离和传播速度(声速)(或在被检体内部声波的传播速度(声速))调整相位的检测信号来计算生成的图像的各像素(或体素)。但是，在傅立叶域方法的情况下，通过在频率空间中的操作生成图像。在图1中，i表示检测元件的数量(0～N-1范围中的整数)，并且，t表示时间。
以下，将具体描述作为一般图像生成方法的时域方法。首先，对于各检测元件，将从第i个检测元件到与生成的图像中的像素对应的位置的距离除以声速。通过使用该结果，计算当在该像素位置中生成声波时所述声波到达检测元件i的时间Ti(延迟时间)(生成时间被假定为t＝0)。对于各检测元件计算时间Ti的检测信号的强度S(i,Ti)，并且将结果相加，由此生成所述生成的图像的像素(通过以相同的方式对于多个位置生成像素来生成所述生成的图像)。根据在专利文献(PTL)1中公开的技术，声速即时间Ti在该生成图像的生成方法中被确定，使得，各体素(或像素)的数据被最大化。强度S(i,Ti)表示时间Ti的检测信号S(i,t)的强度(经相位调整的检测信号的强度)。
这里，如果作为声音速度估计的声速与实际的声速显著不同，那么强度S(i,Ti)大大分散。为了防止这一点，根据本发明，在被检体内部的一部分区域中设定目标区域5(要被关注的区域)。然后，对于多个假定的声速，执行基于从各检测元件i到目标区域5的距离和假定声速调整各检测信号S(i,t)的相位、以及计算经相位调整的多个检测信号的强度S(i,Ti)的离散度的处理。然后，多个假定的声速之中的将上述的强度S(i,Ti)的离散度最小化的声速被确定为声速，并且，基于该声速生成所述生成的图像。
例如，可使用由下式(101)给出的相干因子(CF)作为延迟信号值的离散度的指数。
[数学式1]
CF = | &Sigma; i = 0 N - 1 S ( i ) | 2 N &Sigma; i = 0 N - 1 | S ( i ) | 2 - - - ( 101 )
这里，S(i)是从检测元件i输出并且经相位调整的检测信号的强度，即，基于声速计算的时间Ti的检测信号的强度S(i,Ti)。如果S(i)与检测元件无关地变为相同的值，那么指数变为CF＝1，并且，随着值S(i)的离散度增大，指数CF更接近0。换句话说，当强度S(i,Ti)的离散度越大时，指数CF的值变得越接近0，并且，当离散度越小时，所述值变得越接近1。在本发明中，强度S(i,Ti)的离散度被最小化时的声速即指数CF被最大化的声速被确定为声速。指数不限于CF。可通过在统计领域中使用的诸如分散值和标准偏差值的尺度评价离散度。
因此，根据本发明，强度S(i,Ti)的离散度被最小化的声速被设为声速，由此可基于最接近实际声速的声速生成所述生成的图像。并且，在本实施例中，可通过限制强度S(i,Ti)的离散度最小化的区域来抑制背景噪声的增大。并且，根据本发明，强度S(i,Ti)的离散度被最小化的声速被视为声速，因此，即使在检测信号中随机包含噪声，噪声也可被均匀化(换句话说，噪声的影响可被最小化)。结果，即使在检测信号中包含噪声，也可精确地估计实际的声速。
因此，根据本发明，可以在抑制由于噪声和用于生成图像的声速与实际声速的不匹配导致的图像质量的下降的同时生成代表被检体内部的信息的图像。
<实施例>
现在将描述根据本实施例的图像生成装置和图像生成方法。图2表示根据本实施例的图像生成装置的配置的例子。这里，作为例子，将描述检测从被检体内部传播的声波并将生物信息可视化的使用光声断层成像技术的图像生成装置(光声断层成像装置)的情况。本发明也可被应用于通过从声波探测器传送声波(超声波)并检测在被检体内反射的反射声波(超声波)来生成代表被检体内部的信息的图像的使用超声诊断技术的图像生成装置(超声诊断装置)。
根据本实施例的图像生成装置由光源11、光学装置13、声波探测器17、信号收集器18、图像生成处理单元19、信号处理器20和显示装置21构成。从光源11发射的光12(经强度调制的光)通过光学装置13被照射到诸如活体的被检体15上。当在被检体15内传播的光的能量的一部分被诸如血管的光吸收体(声源)14吸收时，通过光吸收体14的热膨胀生成声波(一般为超声波)16。通过声波探测器17检测所生成的声波16，并且，通过随后的处理生成被检体的生物信息的图像。
光源11发射具有被构成活体的成分之中的特定的成分吸收的特定波长的光。光源11可与根据本实施例的图像生成装置一体化，或者可以是单独的单元。光源11由能够生成几纳秒到几百纳秒量级的脉冲光的一个或更多个脉冲光源构成。对于光源，由于激光器表现高的输出，因此，激光器是优选的，但是，作为激光器的替代，可以使用发光二极管等。对于激光器，可以使用各种激光器，包括固态激光器、气体激光器、染料激光器和半导体激光器。通过未示出的控制单元控制照射定时、光的波形、强度和其它的因素。
从光源11发射的光12通过光学装置13被引导到被检体(它也可被光学波导等引导)。例如光学装置13是反射光的反射镜以及收集或展开光的透镜。对于该光学装置，只要从光源发射的光12可在被检体15上被照射为希望的形状，就可以使用任何装置。一般地，考虑到活体的安全性并具有宽的诊断区域，优选将光12展开到一定的区域而不是通过透镜收集光12。优选被检体上的照射光的区域(照射区域)是可移动的(可变化的)。换句话说，优选本实施例的图像生成装置被构建为使得从光源生成的光可在被检体上移动。然后，可在更宽的范围上生成所述生成的图像。优选照射区域(在被检体上照射的光)可与声波探测器17同步地移动。可用于移动照射区域的方法为：使用可移动反射镜和在机械上移动光源自身等。
例如，对于人和动物的恶性肿瘤和血管疾病的诊断使用根据本实施例的图像生成装置。因此，作为被检体15，可以假定诸如人和动物的乳房、手指和肢体(limb)的诊断目标区域。在光声断层成像装置中，光吸收体(或声源)14是被检体的吸收系数高的部分。例如，在人体为测量目标的情况下，例如，氧合血红蛋白(oxygenatedhemoglobin)、还原血红蛋白(reduced hemoglobin)、包含大量的这些血红蛋白的血管、以及包含许多新血管的恶性肿瘤可以是光吸收体。在超声诊断装置的情况下，光源14是声学阻抗与周围区域不同的组织界面。
声波探测器17具有多个检测元件。检测元件由利用压电现象的换能器、利用光的共振的换能器、或利用电容的变化的换能器构成。配置不限于此，而可以是允许检测声波的任何配置。在声波探测器17中，一般在一维或二维配置中设置多个检测元件。通过使用这样设置的多个检测元件，可依次同时在多个位置中检测声波。因此，可减少检测时间并可使诸如被检体的振动之类的影响最小化。如果通过移动一个检测元件的位置在多个位置中检测声波，那么可以获得与使用多个检测元件的情况类似的信号(检测信号)。
信号收集器18放大通过声波探测器17(各检测元件)获得的电信号(模拟信号)，并且将其转换成数字信号(检测信号)。信号收集器18一般由放大器、A/D转换器和FPGA(场可编程门阵列)芯片和其它的部件构成。优选可同时处理多个检测信号。然后，可减少直到生成图像的时间。
图像生成处理单元19通过使用从多个检测元件(图像生成单元)获得的多个检测信号(从信号收集器18输出的多个数字信号)生成代表被检体内部的信息的图像(生成图像)。在一般的时域方法的各情况下，通过使用基于从各检测元件到与该像素对应的位置的距离调整相位的检测信号和通过后面描述的信号处理器20确定的声速，计算生成图像的各像素。
与图像生成处理单元19连接的信号处理器20确定被检体内部的声波的声速(确定单元)。根据本实施例，通过使用从信号收集器18输出的多个数字信号和通过图像生成处理单元19获得的图像来确定声速。
显示装置21是用于显示通过图像生成处理单元19生成的生成图像的装置。对于显示装置21，例如，可以使用液晶显示器、等离子体显示器、有机EL显示器和具有电子发射器件的显示器。
现在将参照图3中的流程图描述通过图像生成处理单元19和信号处理器20的处理。
首先，图像生成处理单元19基于被检体的估计的声速(预先确定的声速)生成假定的生成图像(体数据)(步骤301)。对于用于生成假定的生成图像的方法，可以使用常规上用于光声断层成像装置和超声诊断装置的方法。例如，可以使用时域或傅立叶域中的背投影方法。
在步骤301中，优选使用诸如超声诊断装置和光声断层成像装置的通过检测声波生成所述生成的图像的图像生成装置。但是，本发明不限于此，可通过基于诸如X射线CT和MRI的不同原理使用图像生成装置生成关于生物信息的体数据(生成的图像)。
然后，信号处理器20选择在步骤301中获得的体数据(生成的图像)中的强烈地反映生物信息的像素(目标像素)，并且确定包含与目标像素对应的位置的区域作为目标区域(步骤302)。在光声断层成像装置的情况下，目标像素是构成吸收体的像素，并且，在超声生成的图像的情况下，目标像素是超声波的反射高的区域中的像素。
可通过用户检查生成的图像来手动选择目标像素，或者可自动选择目标像素。在自动选择目标像素的情况下，选择例如图像中的亮度或对比度最高的像素作为目标像素。
基于假定的声速的后面描述的选择范围确定目标区域的尺寸。例如，考虑通过使用1500m/秒(m/sec)声速(估计的声速c0：预先确定的声速)获得生成图像、然后将声速切换到1400m/秒～1600m/秒范围中的另一声速(假定声速)并重新生成图像的情况。这里，像素节距(pitch)d为0.25mm，并且，被检体的厚度r(从被检体中的与目标像素对应的位置到声波探测器的距离)为40mm。在这种情况下，可能与目标像素对应的位置沿x方向、y方向和z方向最多偏移约10个像素。在这种情况下，选择与沿x方向、y方向和z方向从目标像素起包含10个像素的范围(例如，包含目标像素周围的21*21*21个像素的范围)对应的区域作为目标区域。
具体而言，与使用估计的声速c0(预先确定的声速)的情况下的目标像素对应的位置在改变使用的声速的情况下从该位置沿宽度改变约M像素。通过下式(102)获得值M。在式(102)中，cmin是假定声速的最小值，并且，cmax是假定声速的最大值。
[数学式2]
M = c 0 d ( r c min - r c max ) - - - ( 102 )
因此，与目标像素周围的通过式(102)获得的宽度M的范围对应的区域可被设为目标区域。换句话说，如果生成三维图像，那么选择与包含目标像素周围的M*M*M个体素的范围对应的区域作为目标区域。如果生成二维图像，那么选择与包含目标像素周围的M*M个像素的范围对应的区域作为目标区域。由此，包含诸如吸收体的声源的区域可以是目标区域。
特别地，式(102)中的被检体的厚度r是从与目标像素对应的位置到最远离该位置的检测元件的距离。但是，被检体的厚度不限于此，而可以是从与目标像素对应的位置到最接近该位置的检测元件的距离。或者，被检体的厚度可以是与目标像素对应的位置与各检测元件之间的距离的平均值。
然后，信号处理器20基于从各检测元件到目标区域(例如，与生成图像的像素对应的目标区域中的位置)的距离和假定的声速调整各检测信号的相位。然后，信号处理器20计算经相位调整的多个检测信号的强度的离散度(步骤303)。对于多个假定的声速执行步骤303中的处理。
然后，信号处理器20确定多个假定的声速之中的强度的离散度最小化的声速作为声速(步骤304)。图像生成处理单元19基于所确定的声速生成所述生成的图像。
通过以上的处理，可以在抑制由于噪声和用于生成图像的声速与实际声速的不匹配导致的图像质量的下降的同时生成代表被检体内部的信息的图像。
如果检测元件的面积大，那么检测信号强烈地受方向性影响。将参照图4A和图4B描述方向性。方向性是检测元件的检测敏感度根据检测元件与被检体内部的位置(与像素对应的位置：处理目标位置)之间的相对位置关系改变的特性。例如，如果检测元件如图4A所示的那样面向处理目标位置，那么检测来自该位置的声波的敏感度最大(检测信号的强度最大)。随着检测元件的方向与处理目标位置之间的偏差变大，检测敏感度降低(检测信号的强度降低)。当元件的宽度大时，或者当要被检测的声波的频率高时，方向性的该影响变得明显。方向性R(检测敏感度)由式(103)给出。
[数学式3]
R = | sin ( k &CenterDot; d sin &theta; ) k &CenterDot; d sin &theta; | - - - ( 103 )
d：检测元件的宽度
θ：由连接处理目标位置和检测元件的线与检测元件的方向形成的角度
k：2π/λ
λ：声波的波长
检测元件的形状被假定为正方形。
图4B表示元件宽度为2mm的检测元件的方向性的例子。这里，通过使用从检测元件的脉冲响应获得的频率特性计算各频率分量的大小，并且，通过使用式(103)计算方向性(通过最大值1标准化的检测敏感度)。如图4B所示，如果由连接处理目标位置和检测元件的线与检测元件的方向形成的角度增大，那么敏感度(检测信号的强度)突然下降。这种方向性增大检测信号间的离散度。
因此，根据与像素对应的位置与检测元件的相对位置关系，优选本实施例的图像生成装置还具有校正由于检测元件的方向性导致的检测信号的强度的下降的功能。并且，优选信号处理器20通过使用经强度校正的检测信号确定声速，并且，图像生成处理单元19通过使用经强度校正的检测信号生成图像。由此，可以抑制检测元件的方向性的影响。例如，可减小由于方向性导致的检测信号的离散度(由于声速导致的离散度变为支配性的)，因此可以更精确地确定声速。并且，由于方向性导致的检测信号的强度的离散度减小，因此，可以生成图像质量不下降的图像。
对于校正，例如，基于处理目标位置与检测元件的相对位置关系计算由连接处理目标位置和检测元件的线与检测元件的方向形成的角度，并且，从检测元件输出的检测信号的强度基于根据所述角度的方向性的值(例如，图4B中的Y轴的值)增大。具体而言，检测信号被乘以方向性的值的倒数。然后，可以抑制由于方向性导致的检测信号的强度的离散度。
在许多的情况下，通过执行软件(程序)的计算机(例如，工作站)的处理单元实现图像生成处理单元19和信号处理器20的处理。因此，可以不区分图像生成处理单元19和信号处理器20。可分别通过专用的芯片构建这些功能元件。
可通过硬件代替软件实现图像生成处理单元19和信号处理器20的处理。
根据本实施例，通过使用目标像素的位置和式(102)确定目标区域，但是，本发明不限于此。如果可估计目标区域的一般位置，诸如在测量预先确定的被检体的情况下，可事先在装置中存储关于目标区域的信息。可对于被检体的各种类型(测量区域)事先存储关于预先确定的目标区域的信息。在这种情况下，不需要生成本实施例的步骤301中的假定的生成图像。目标区域可以是任何尺寸(甚至一个像素尺寸)。例如，目标区域可被设定，使得包含与在步骤302中选择的目标像素对应的位置。由于声速被假定为具有一定的尺寸，因此，即使与目标像素对应的位置由于声速的变化而偏移到目标区域之外，这样设定的目标区域也可从声源接收声音。
<例子1>
模拟本发明，并且，验证其效果。将参照图5A、图5B和图6描述本例子。
作为被检体，如图5A所示，使用X方向为4cm、Y方向为4cm且Z方向为4cm的立方体假想体(phantom)43。在立方体假想体43的中心，设置直径为0.05cm并且高度为2cm的圆柱声源41，使得其轴向与X方向平行。这里，声源41在光声断层成像的情况下是光吸收体，并且，在超声诊断装置的情况下是超声反射器。假定被检体内部的声波的声速为1500m/秒，并且，在声源41内生成的初始生成的声压为3000Pa(帕斯卡)。
声波探测器42是二维阵列类型，并且，由20*20个正方形检测元件构成。各元件的一个边的长度为2mm，并且，节距距离为2mm。假定在1280个测量点处以20MHz的采样频率检测声波。
在以上的条件下，通过物理模拟生成通过声波探测器42的各元件检测的声波的检测信号。图5B表示通过探测器的中心的检测元件获得的检测信号的例子。图5B中的强的信号表示来自声源41的信号。假定检测信号包含三倍标准偏差为20Pa的白噪声。
通过使用以这种方式获得的检测信号和1580m/秒的声速(预先确定的声速)，生成所述生成的图像。对于生成方法，使用作为公知技术的根据时域方法的通用背投影。成像范围为X方向3.8cm、Y方向3.8cm和Z方向4cm，并且，像素节距(体素节距)为0.025cm。通过将成像范围除以体素节距，获取体素的总数即152*152*160。
图6表示生成的图像。图6所示的所有图像是沿声波探测器42的方向投影体素值的最大强度的MIP(最大强度投影)图像，并且，白区域表示高声压。在以上的各图像上示出Y＝2.0cm处的X方向的图像的声压变化。
在图6中，图像A表示被检体内部的声源的MIP图像。根据本例子，假定分辨率在图像更接近图像A时变得更高。
图像B是当使用1580m/秒声速作为声速时获得的生成图像。如该图像所示，如果使用与1500m/秒的实际声速不同的声速，那么图像模糊并且分辨率下降。
然后，在该图像(体数据)中，搜索亮度最大的体素。结果，检测到中心处的体素(所有152*152*160个体素中的第76*76*80个体素)。与以该体素为中心的周围21*21*21个体素对应的区域被设为目标区域。
然后，为了与现有技术相比较，在以2m/秒的间隔将使用的声速从1400m/秒变为1600m/秒的同时，对于目标区域中的各位置(体素)计算亮度值。换句话说，计算21*21*21*101个亮度值。结果，亮度值最大化的声速为1456m/秒。然后，通过使用该声速1456m/秒，重新生成图像。图像C是在此时获得的图像。该图像有些模糊，但是，与图像B相比，分辨率稍微改善了。最大亮度值也增大了。
为了实证本发明的效果，在以2m/秒的间隔将使用的声速从1400m/秒变为1600m/秒的同时，对于目标区域中的各位置(体素)计算检测信号的强度的离散度(指数CF)。换句话说，计算21*21*21*101个指数CF。结果，指数CF最大化的声速为1492m/秒。然后，通过使用该声速1492m/秒，重新生成图像。图像C′是在此时获得的图像。与作为现有技术的图像C相比，图像模糊更少并且分辨率提高。
[表1]
速度 图像质量(模糊量)
图像A(实际图像) 1500m/秒 无
图像B(估计) 1580m/秒 高
图像C(现有技术) 1456m/秒 中
图像C′(本例子) 1492m/秒 低
以这种方式，根据本发明，可以估计接近实际声速的声速(1492m/秒)。特别是在包含噪声的检测信号中，与通过现有技术获得的声速(1456m/秒)相比，可以估计更接近实际声速的声速。结果，可以获得没有分辨率的下降的图像。
虽然实际声速为1500m/秒，但是，由于噪声和检测元件的方向性的影响，声速被估计为1492m/秒。
在本例子中，对于目标区域中的各位置计算检测信号的强度的离散度，但是，可对于目标区域中的一部分位置计算所述离散度。例如，可仅对于目标区域的中心位置计算所述离散度。
<例子2>
现在将参照图2描述将本发明应用于光声断层成像装置的例子。在本例子中，使用对于约10纳秒生成具有1064nm波长的脉冲光的Q开关YAG激光器作为光源11。脉冲激光束(光12)的能量为0.6J，并且，该脉冲光通过使用诸如反射镜和扩束器的光学装置13展开为约2cm半径。对于被检体15，使用模拟活体的假想体。使用的假想体为通过琼脂(agar)凝固的1％脂肪乳(intralipid)。假定体的尺寸为宽度15cm、高度8cm和深度4cm。在该假想体中，在中心区域中，嵌入具有0.03cm的直径的黑色橡胶线作为光吸收体14。展开为2cm半径的光12被照射到该假定体上。对于声波探测器17，使用由PZT(锆钛酸铅)制成的超声换能器。该换能器是二维阵列类型，这里，元件的数量为18*18个，各元件的形状为正方形，并且，元件节距为2mm。各元件的宽度为约2mm。该声波探测器可与光照射区域同步地沿X方向和Y向移动以对大的区域成像。
如果脉冲光被照射到假想体的一个表面上，那么在假想体内扩散的光被橡胶线吸收，并且，生成光声波。该光声波被具有324个通道(channel)的各超声换能器检测，并且，通过使用由放大器、A/D转换器和FPGA构成的信号收集器18，对于各通道获得光声波的数字数据(检测信号)。为了提高信号的S/N比，激光被照射36次，并且，所有获得的检测信号被平均化。然后，获得的数字数据被传送到作为图像生成处理器19和信号处理器20的工作站(WS)，并且，WS存储该数据。在通过离散小波变换对于该存储的数据执行降噪处理之后，通过使用作为活体内部的声波的平均声速的1540m/秒声速生成图像。这里，通过使用作为时域方法的通常背投影方法生成体数据。在这种情况下使用的体素节距为0.05cm。成像范围被设为11.8cm*11.8cm*4.0cm。图7A表示获得的图像的例子。
然后，在以2m/秒间隔将使用的声速从1400m/秒变为1600m/秒的同时，对于目标区域中的各位置(体素)计算检测信号的强度的离散度(指数CF)。根据本例子，在获得的体数据中，包含在所述线中生成的声波(初始声压)的与体素组(21*21*21个体素)对应的区域被设为目标区域。并且，根据本例子，考虑基于检测元件的脉冲响应计算的方向性校正检测信号，并且，通过使用校正的检测信号计算指数CF。通过考虑方向性，最大指数CF从0.056提高到0.078。在本例子中，指数CF变得最大的声速为1454m/秒。然后，使用该1454m/秒声速和存储于PC中的检测信号，重新生成图像。在作为显示装置21的液晶显示器上显示获得的图像。图7B表示获得的图像的例子。
如以上的结果所示，在通过使用1454m/秒声速生成的图像(图7B)中，与通过使用1540m/秒声速生成的图像(图7A)相比，从橡胶线生成的声波(初始声压)的宽度明显更窄。显然几乎不存在任何图像模糊，换句话说，分辨率提高了。
以这种方式，根据本例子，可精确地估计被检体内部的声波的声速，并可因此提高生成图像的分辨率。
也可通过读出并执行记录在存储设备上的程序以执行上述的实施例的功能的系统或装置的计算机(或诸如CPU或MPU的设备)以及通过由系统或装置的计算机通过例如读出并执行记录在存储设备上的程序以执行上述的实施例的功能执行其各个步骤的方法，实现本发明的各方面。出于这种目的，例如通过网络或从用作存储设备的各种类型的记录介质(例如，非暂态计算机可读介质)向计算机提供程序。
虽然已参照示例性实施例说明了本发明，但应理解，本发明不限于公开的示例性实施例。所附权利要求的范围应被赋予最宽泛的解释以包含所有这样的修改以及等同的结构和功能。
本申请要求在2009年12月11日提交的日本专利申请No.2009-281452的权益，在此通过引用并入其全部内容。

Claims (25)

1.一种图像生成装置，包括：
多个检测元件，所述多个检测元件被配置为检测从被检体传播的声波并输出检测信号；
确定单元，所述确定单元被配置为确定声波在被检体内部的传播速度，以用于图像生成；和
生成单元，所述生成单元被配置为通过使用由所述确定单元确定的传播速度和从所述多个检测元件输出的多个检测信号来生成代表图像中的每个像素或体素的信息的图像数据，
其中，所述生成单元被配置为基于预先确定的传播速度和所述多个检测信号生成假定的图像数据，并且
其中，所述确定单元被配置为：
选择所述假定的图像数据之中的反映被检体信息的目标区域，
设定假定的传播速度，
基于从所述检测元件中的每一个到目标区域的距离和所述假定的传播速度，确定所述多个检测信号之中的对应于目标区域的多个信号，以及
基于对应于目标区域的所述多个信号，确定传播速度。
2.根据权利要求1的图像生成装置，其中所述确定单元被配置为选择所述假定的图像数据之中的具有最高亮度值的目标区域。
3.根据权利要求1的图像生成装置，其中，所述确定单元被配置为选择所述假定的图像数据之中的具有最高对比度的目标区域。
4.根据权利要求1的图像生成装置，其中，所述确定单元被配置为选择所述假定的图像数据之中的包含声源的区域作为目标区域。
5.根据权利要求1的图像生成装置，其中，所述确定单元被配置为：
对于多个假定的传播速度，执行确定对应于目标区域的多个信号的处理，并且
基于通过使用所述多个假定的传播速度确定的对应于目标区域的所述多个信号，确定所述多个假定的传播速度之中的一假定的传播速度作为传播速度。
6.根据权利要求1的图像生成装置，其中，所述确定单元被配置为：
基于从所述检测元件中的每一个到目标区域的距离和所述假定的传播速度，确定从目标区域生成的声波到达各检测元件的多个定时，并且
基于所述多个定时确定对应于目标区域的多个信号。
7.根据权利要求1的图像生成装置，还包括：
校正单元，所述校正单元被配置为校正由于检测元件的方向性导致的检测信号的强度的降低，
其中，所述确定单元被配置为通过使用经强度校正的检测信号来确定传播速度。
8.根据权利要求1的图像生成装置，其中，所述确定单元被配置为：
对于多个假定的传播，基于从所述多个检测元件到目标区域的距离和假定的传播速度，调整所述多个检测信号的相位，
获得对应于多个假定的传播速度之中的各个假定的传播速度的经相位调整的多个检测信号的强度的离散度，并且
基于对应于所述多个假定的传播速度的离散度，确定所述多个假定的传播速度之中的一假定的传播速度作为传播速度。
9.根据权利要求1的图像生成装置，其中，所述确定单元被配置为基于由下式表达的指数CF确定传播速度：
[数学式1]
CF = | &Sigma; i = 0 N - 1 S ( i ) | 2 N &Sigma; i = 0 N - 1 | S ( i ) | 2
这里，N是检测元件的数量，N是2或更大的整数，并且，S(i)是从所述检测元件之中的第i个检测元件输出的对应于目标区域的信号的强度，i是0到N-1的范围中的整数。
10.根据权利要求1的图像生成装置，还包括光源，
其中，从被检体传播的声波是通过用从光源发射的光照射被检体而生成的光声波。
11.根据权利要求10的图像生成装置，其中，所述确定单元被配置为选择所述假定的图像数据之中的包含光吸收体的区域作为目标区域。
12.根据权利要求10的图像生成装置，其中，所述确定单元被配置为选择所述假定的图像数据之中的包含血管的区域作为目标区域。
13.根据权利要求1的图像生成装置，其中，所述生成单元被配置为通过使用传播速度和所述多个检测信号来生成代表被检体内部的信息的三维图像数据。
14.一种图像生成装置，包括：
光源；
多个检测元件，所述多个检测元件被配置为检测通过用从光源发射的光照射被检体而生成的光声波并输出检测信号；
确定单元，所述确定单元被配置为确定光声波在被检体内部的传播速度；和
生成单元，所述生成单元被配置为通过使用由所述确定单元确定的传播速度和从所述多个检测元件获得的多个检测信号来生成代表被检体内的信息的图像数据，
其中，所述生成单元被配置为基于预先确定的传播速度和所述多个检测信号生成假定的图像数据，并且
其中，所述确定单元被配置为：
选择所述假定的图像数据之中的反映被检体信息的目标区域，
确定对应于目标区域的传播速度。
15.根据权利要求14的图像生成装置，其中，所述确定单元被配置为选择所述假定的图像数据之中的具有最高亮度值的区域作为目标区域。
16.根据权利要求14的图像生成装置，其中，所述确定单元被配置为选择所述假定的图像数据之中的具有最高对比度的区域作为目标区域。
17.根据权利要求14的图像生成装置，其中，所述确定单元被配置为选择所述假定的图像数据之中的包含声源的区域作为目标区域。
18.根据权利要求14的图像生成装置，其中，所述确定单元被配置为选择所述假定的图像数据之中的包含光吸收体的区域作为目标区域。
19.根据权利要求14的图像生成装置，其中，所述确定单元被配置为选择所述假定的图像数据之中的包含血管的区域作为目标区域。
20.一种图像生成方法，包括：
发射步骤，发射光；
检测步骤，检测通过用在发射步骤中发射的光照射被检体而生成的光声波；
确定步骤，确定光声波在被检体内部的传播速度；以及
生成步骤，通过使用由在所述确定步骤中确定的传播速度和多个检测信号来生成代表被检体内的信息的图像数据，
其中，在所述确定步骤中，
基于预先确定的传播速度和所述多个检测信号生成假定的图像数据，
从所述假定的图像数据之中选择反映被检体信息的目标区域，并且
基于对应于目标区域的多个信号，确定传播速度。
21.根据权利要求20所述的图像生成方法，其中，在所述确定步骤中，所述假定的图像数据之中的具有最高亮度值的区域被选择作为目标区域。
22.根据权利要求20所述的图像生成方法，其中，在所述确定步骤中，所述假定的图像数据之中的具有最高对比度的区域被选择作为目标区域。
23.根据权利要求20所述的图像生成方法，其中，在所述确定步骤中，所述假定的图像数据之中的包含声源的区域被选择作为目标区域。
24.根据权利要求20所述的图像生成方法，其中，在所述确定步骤中，所述假定的图像数据之中的包含光吸收体的区域被选择作为目标区域。
25.根据权利要求20所述的图像生成方法，其中，在所述确定步骤中，所述假定的图像数据之中的包含血管的区域被选择作为目标区域。
CN201410266998.XA 2009-12-11 2010-12-08 图像生成装置和图像生成方法 Active CN104013438B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
JP2009281452A JP5528083B2 (ja) 2009-12-11 2009-12-11 画像生成装置、画像生成方法、及び、プログラム
JP2009-281452 2009-12-11
CN201080055219.1A CN102640014B (zh) 2009-12-11 2010-12-08 图像生成装置和图像生成方法

Related Parent Applications (1)

Application Number Title Priority Date Filing Date
CN201080055219.1A Division CN102640014B (zh) 2009-12-11 2010-12-08 图像生成装置和图像生成方法

Publications (2)

Publication Number Publication Date
CN104013438A true CN104013438A (zh) 2014-09-03
CN104013438B CN104013438B (zh) 2016-09-07

Family

ID=43718662

Family Applications (2)

Application Number Title Priority Date Filing Date
CN201080055219.1A Active CN102640014B (zh) 2009-12-11 2010-12-08 图像生成装置和图像生成方法
CN201410266998.XA Active CN104013438B (zh) 2009-12-11 2010-12-08 图像生成装置和图像生成方法

Family Applications Before (1)

Application Number Title Priority Date Filing Date
CN201080055219.1A Active CN102640014B (zh) 2009-12-11 2010-12-08 图像生成装置和图像生成方法

Country Status (5)

Country Link
US (4) US8687868B2 (zh)
EP (1) EP2510382B1 (zh)
JP (1) JP5528083B2 (zh)
CN (2) CN102640014B (zh)
WO (1) WO2011070778A1 (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110403681A (zh) * 2018-04-27 2019-11-05 佳能医疗系统株式会社 超声波诊断装置以及穿刺辅助方法

Families Citing this family (28)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9052268B2 (en) * 2011-03-22 2015-06-09 Fujifilm Corporation Ultrasound diagnostic apparatus and method of producing ultrasound image
JP5281107B2 (ja) * 2011-03-25 2013-09-04 富士フイルム株式会社 超音波診断装置および超音波画像生成方法
JP5850633B2 (ja) 2011-04-12 2016-02-03 キヤノン株式会社 被検体情報取得装置
JP5863345B2 (ja) 2011-09-08 2016-02-16 キヤノン株式会社 被検体情報取得装置および被検体情報取得方法
CN104053402B (zh) * 2012-01-18 2017-04-19 佳能株式会社 被检体信息获得设备和关于被检体的信息的获得方法
JP6071260B2 (ja) 2012-06-13 2017-02-01 キヤノン株式会社 被検体情報取得装置および情報処理方法
US10422900B2 (en) * 2012-11-02 2019-09-24 Exxonmobil Upstream Research Company Analyzing seismic data
US20140182383A1 (en) 2012-12-28 2014-07-03 Canon Kabushiki Kaisha Object information obtaining device, display method, and non-transitory computer-readable storage medium
JP6222936B2 (ja) 2013-02-19 2017-11-01 キヤノン株式会社 装置および画像生成方法
JP2015167789A (ja) * 2014-03-10 2015-09-28 キヤノン株式会社 被検体情報取得装置および信号処理方法
EP3126969A4 (en) 2014-04-04 2017-04-12 Microsoft Technology Licensing, LLC Expandable application representation
CN105359055A (zh) 2014-04-10 2016-02-24 微软技术许可有限责任公司 计算设备的滑盖
JP6352050B2 (ja) 2014-05-19 2018-07-04 キヤノンメディカルシステムズ株式会社 超音波診断装置
JP6362420B2 (ja) * 2014-05-26 2018-07-25 キヤノン株式会社 被検体情報取得装置、被検体情報取得方法、および、プログラム
JP6371847B2 (ja) * 2014-06-30 2018-08-08 富士フイルム株式会社 光音響画像生成装置、信号処理装置、及び光音響画像生成方法
EP3174490B1 (en) 2014-07-28 2020-04-15 Intuitive Surgical Operations, Inc. Systems and methods for planning multiple interventional procedures
CN106662891B (zh) 2014-10-30 2019-10-11 微软技术许可有限责任公司 多配置输入设备
WO2016072080A1 (en) 2014-11-07 2016-05-12 Canon Kabushiki Kaisha Object information acquiring apparatus and method using a photoacoustic effect
JP6544910B2 (ja) 2014-11-07 2019-07-17 キヤノン株式会社 情報処理装置、被検体情報取得装置及び音速決定方法
JP6456129B2 (ja) * 2014-12-15 2019-01-23 キヤノン株式会社 被検体情報取得装置およびその制御方法ならびに光量制御方法
JP6521761B2 (ja) 2015-06-23 2019-05-29 キヤノン株式会社 情報処理装置および表示制御方法
JP6512969B2 (ja) * 2015-07-06 2019-05-15 キヤノン株式会社 処理装置、光音響装置、処理方法、及びプログラム
JP6632368B2 (ja) * 2015-12-21 2020-01-22 キヤノン株式会社 情報処理装置、光音響装置、情報処理方法、及びプログラム
JP6084313B2 (ja) * 2016-01-14 2017-02-22 キヤノン株式会社 被検体情報取得装置及び被検体情報取得方法
JP6400068B2 (ja) * 2016-12-20 2018-10-03 キヤノン株式会社 情報処理装置および情報処理方法
JP6537540B2 (ja) * 2017-01-25 2019-07-03 キヤノン株式会社 処理装置
CN110840484B (zh) * 2019-11-27 2022-11-11 深圳开立生物医疗科技股份有限公司 自适应匹配最优声速的超声成像方法、装置及超声设备
JP7479684B2 (ja) 2020-08-19 2024-05-09 株式会社Luxonus 画像生成方法、画像生成プログラムおよび画像生成装置

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030167002A1 (en) * 2000-08-24 2003-09-04 Ron Nagar Photoacoustic assay and imaging system
JP2008049063A (ja) * 2006-08-28 2008-03-06 Osaka Prefecture Univ 光トモグラフィ装置用プローブ
EP1935346A1 (en) * 2006-12-21 2008-06-25 Stichting voor de Technische Wetenschappen Imaging apparatus and method
CN101273903A (zh) * 2007-03-28 2008-10-01 株式会社东芝 超声波成像装置以及超声波速度优化方法
CN101406401A (zh) * 2007-10-10 2009-04-15 株式会社东芝 超声波图像取得装置

Family Cites Families (35)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS59171542A (ja) * 1983-03-18 1984-09-28 三菱電機株式会社 超音波診断装置
JPS61154651A (ja) 1984-12-28 1986-07-14 株式会社東芝 超音波診断装置
JPH03146039A (ja) * 1989-11-02 1991-06-21 Toshiba Corp 超音波診断装置
JP3642607B2 (ja) * 1995-05-26 2005-04-27 株式会社日立メディコ 超音波断層装置
JPH09133655A (ja) 1995-11-10 1997-05-20 Hitachi Ltd 光音響分光分析装置及び光音響分光分析方法
JP4334032B2 (ja) * 1998-02-23 2009-09-16 株式会社東芝 超音波診断装置
JP2000139914A (ja) 1998-11-04 2000-05-23 Aloka Co Ltd 超音波診断装置
KR20000038847A (ko) 1998-12-09 2000-07-05 이민화 초음파영상화시스템에서의 집속방법
JP4416256B2 (ja) * 2000-03-10 2010-02-17 株式会社日立メディコ 超音波撮像装置
JP4958348B2 (ja) * 2001-09-06 2012-06-20 株式会社日立メディコ 超音波撮像装置
US6705994B2 (en) * 2002-07-08 2004-03-16 Insightec - Image Guided Treatment Ltd Tissue inhomogeneity correction in ultrasound imaging
JP4406226B2 (ja) 2003-07-02 2010-01-27 株式会社東芝 生体情報映像装置
JP4643153B2 (ja) 2004-02-06 2011-03-02 株式会社東芝 非侵襲生体情報映像装置
US7862508B2 (en) * 2004-09-20 2011-01-04 Innervision Medical Technologies Inc. Systems and methods for ultrasound imaging
US20070083110A1 (en) 2005-10-09 2007-04-12 Sonowise, Inc. Programmable phase velocity in an ultrasonic imaging system
JP4755016B2 (ja) * 2006-05-11 2011-08-24 日本電信電話株式会社 成分濃度測定装置
CN101453939A (zh) 2006-05-25 2009-06-10 皇家飞利浦电子股份有限公司 光声成像方法
US20080173093A1 (en) * 2007-01-18 2008-07-24 The Regents Of The University Of Michigan System and method for photoacoustic tomography of joints
JP4739363B2 (ja) 2007-05-15 2011-08-03 キヤノン株式会社 生体情報イメージング装置、生体情報の解析方法、及び生体情報のイメージング方法
EP2002784B1 (en) 2007-06-11 2018-07-11 Canon Kabushiki Kaisha Intravital-information imaging apparatus
JP2009066110A (ja) * 2007-09-12 2009-04-02 Canon Inc 測定装置
CN101214156B (zh) 2008-01-10 2010-12-15 复旦大学 声速不均匀介质热声成像的重建算法
JP5284129B2 (ja) 2008-02-06 2013-09-11 キヤノン株式会社 イメージング装置、及び解析方法
JP5065153B2 (ja) 2008-05-21 2012-10-31 日本電産トーソク株式会社 電磁弁
CN101336832A (zh) 2008-08-12 2009-01-07 福建师范大学 脉冲式光声扫描软组织成像方法与装置
JP5460000B2 (ja) 2008-08-20 2014-04-02 キヤノン株式会社 イメージング装置およびイメージング方法
JP5496098B2 (ja) 2008-08-27 2014-05-21 キヤノン株式会社 被検体情報取得装置およびその制御方法
JP4900979B2 (ja) 2008-08-27 2012-03-21 キヤノン株式会社 光音響装置および光音響波を受信するための探触子
JP5451014B2 (ja) * 2008-09-10 2014-03-26 キヤノン株式会社 光音響装置
US20100081912A1 (en) * 2008-09-30 2010-04-01 Nellcor Puritan Bennett Llc Ultrasound-Optical Doppler Hemometer and Technique for Using the Same
JP2010088627A (ja) 2008-10-07 2010-04-22 Canon Inc 生体情報処理装置および生体情報処理方法
US9528966B2 (en) * 2008-10-23 2016-12-27 Washington University Reflection-mode photoacoustic tomography using a flexibly-supported cantilever beam
JP2011005042A (ja) 2009-06-26 2011-01-13 Canon Inc 光音響イメージング装置及び光音響イメージング方法
JP5675390B2 (ja) 2010-02-09 2015-02-25 キヤノン株式会社 測定装置
JP5574927B2 (ja) 2010-11-19 2014-08-20 キヤノン株式会社 測定装置

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030167002A1 (en) * 2000-08-24 2003-09-04 Ron Nagar Photoacoustic assay and imaging system
JP2008049063A (ja) * 2006-08-28 2008-03-06 Osaka Prefecture Univ 光トモグラフィ装置用プローブ
EP1935346A1 (en) * 2006-12-21 2008-06-25 Stichting voor de Technische Wetenschappen Imaging apparatus and method
CN101273903A (zh) * 2007-03-28 2008-10-01 株式会社东芝 超声波成像装置以及超声波速度优化方法
CN101406401A (zh) * 2007-10-10 2009-04-15 株式会社东芝 超声波图像取得装置

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110403681A (zh) * 2018-04-27 2019-11-05 佳能医疗系统株式会社 超声波诊断装置以及穿刺辅助方法
CN110403681B (zh) * 2018-04-27 2023-12-19 佳能医疗系统株式会社 超声波诊断装置以及图像显示方法

Also Published As

Publication number Publication date
JP2011120765A (ja) 2011-06-23
US20170071476A1 (en) 2017-03-16
EP2510382A1 (en) 2012-10-17
CN104013438B (zh) 2016-09-07
US20140148683A1 (en) 2014-05-29
CN102640014A (zh) 2012-08-15
US20190059740A1 (en) 2019-02-28
US9579085B2 (en) 2017-02-28
EP2510382B1 (en) 2018-04-11
US20120183190A1 (en) 2012-07-19
US8687868B2 (en) 2014-04-01
WO2011070778A1 (en) 2011-06-16
CN102640014B (zh) 2014-07-02
US10136821B2 (en) 2018-11-27
JP5528083B2 (ja) 2014-06-25

Similar Documents

Publication Publication Date Title
CN102640014B (zh) 图像生成装置和图像生成方法
US9560973B2 (en) Subject information acquiring device and subject information acquiring method
CN103260502B (zh) 被检体信息获取设备和被检体信息获取方法
JP6761755B2 (ja) 合成アパーチャ画像形成におけるコヒーレントなスペクトル拡散符号化波形
JP6504826B2 (ja) 情報処理装置および情報処理方法
JP5441795B2 (ja) イメージング装置及びイメージング方法
CN103025248B (zh) 图像信息获取装置和图像信息获取方法
US10064558B2 (en) Subject information acquisition device, method for controlling subject information acquisition device, and storage medium storing program therefor
JP2014054517A (ja) 光音響画像生成装置及び方法
CN105249993A (zh) 一种通过光声成像选取最佳声速组优化超声成像的方法
CN102821677A (zh) 光声成像装置、光声成像方法和程序
US20170176399A1 (en) Object information acquiring apparatus and control method thereof
CN107205720B (zh) 一种超声适应性波束形成方法及其对经颅成像的应用
Noda et al. Ultrasound imaging with a flexible probe based on element array geometry estimation using deep neural network
JP6296759B2 (ja) 被検体情報取得装置
US20120278010A1 (en) Object information acquiring apparatus
JP5940109B2 (ja) 画像生成装置、伝播速度決定方法、及び、プログラム
KR20230109615A (ko) 매질의 초음파 특성화 방법 및 시스템
JP6113330B2 (ja) 装置および画像生成方法
JP6381718B2 (ja) 装置および画像生成方法
Kuzmin et al. Fast low-cost single element ultrasound reflectivity tomography using angular distribution analysis

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