CN105832299A - 被检体信息获取装置和被检体信息获取方法 - Google Patents
被检体信息获取装置和被检体信息获取方法 Download PDFInfo
- Publication number
- CN105832299A CN105832299A CN201610156516.4A CN201610156516A CN105832299A CN 105832299 A CN105832299 A CN 105832299A CN 201610156516 A CN201610156516 A CN 201610156516A CN 105832299 A CN105832299 A CN 105832299A
- Authority
- CN
- China
- Prior art keywords
- image
- roi
- unit
- interest
- area
- 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
Links
- 238000000034 method Methods 0.000 title abstract description 36
- 238000012545 processing Methods 0.000 claims abstract description 40
- 230000003287 optical effect Effects 0.000 claims abstract description 6
- 238000002604 ultrasonography Methods 0.000 claims description 19
- 238000012800 visualization Methods 0.000 claims description 8
- 238000005259 measurement Methods 0.000 claims description 6
- 230000015572 biosynthetic process Effects 0.000 claims description 4
- 238000003786 synthesis reaction Methods 0.000 claims description 4
- 230000002194 synthesizing effect Effects 0.000 abstract 1
- 238000003384 imaging method Methods 0.000 description 49
- 230000008569 process Effects 0.000 description 20
- 230000006870 function Effects 0.000 description 15
- 230000008859 change Effects 0.000 description 14
- 230000001186 cumulative effect Effects 0.000 description 12
- 238000010586 diagram Methods 0.000 description 12
- 230000005540 biological transmission Effects 0.000 description 9
- 108010054147 Hemoglobins Proteins 0.000 description 5
- 102000001554 Hemoglobins Human genes 0.000 description 5
- 241000222065 Lycoperdon Species 0.000 description 4
- 241000768494 Polymorphum Species 0.000 description 4
- 238000002835 absorbance Methods 0.000 description 4
- 238000009825 accumulation Methods 0.000 description 4
- 238000010521 absorption reaction Methods 0.000 description 3
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 3
- 230000006399 behavior Effects 0.000 description 3
- 239000003086 colorant Substances 0.000 description 3
- 239000004615 ingredient Substances 0.000 description 3
- 239000001301 oxygen Substances 0.000 description 3
- 229910052760 oxygen Inorganic materials 0.000 description 3
- 230000002463 transducing effect Effects 0.000 description 3
- 241000208340 Araliaceae Species 0.000 description 2
- 108010064719 Oxyhemoglobins Proteins 0.000 description 2
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 2
- 235000003140 Panax quinquefolius Nutrition 0.000 description 2
- 239000008280 blood Substances 0.000 description 2
- 210000004369 blood Anatomy 0.000 description 2
- 238000003745 diagnosis Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 235000008434 ginseng Nutrition 0.000 description 2
- 238000001727 in vivo Methods 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- VYMDGNCVAMGZFE-UHFFFAOYSA-N phenylbutazonum Chemical compound O=C1C(CCCC)C(=O)N(C=2C=CC=CC=2)N1C1=CC=CC=C1 VYMDGNCVAMGZFE-UHFFFAOYSA-N 0.000 description 2
- 108010022579 ATP dependent 26S protease Proteins 0.000 description 1
- 206010059245 Angiopathy Diseases 0.000 description 1
- 206010006187 Breast cancer Diseases 0.000 description 1
- 208000026310 Breast neoplasm Diseases 0.000 description 1
- 244000297179 Syringa vulgaris Species 0.000 description 1
- 235000004338 Syringa vulgaris Nutrition 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 201000008275 breast carcinoma Diseases 0.000 description 1
- 201000011510 cancer Diseases 0.000 description 1
- 239000003990 capacitor Substances 0.000 description 1
- 210000004027 cell Anatomy 0.000 description 1
- 238000013098 chemical test method Methods 0.000 description 1
- 238000002512 chemotherapy Methods 0.000 description 1
- 238000004040 coloring Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 239000007789 gas Substances 0.000 description 1
- 239000012528 membrane Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 229910052573 porcelain Inorganic materials 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000010079 rubber tapping Methods 0.000 description 1
- 229920006395 saturated elastomer Polymers 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 201000000849 skin cancer Diseases 0.000 description 1
- 201000008261 skin carcinoma Diseases 0.000 description 1
- 238000003325 tomography Methods 0.000 description 1
- 210000004881 tumor cell Anatomy 0.000 description 1
- 238000012285 ultrasound imaging Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
- 239000011800 void material Substances 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/13—Tomography
- A61B8/14—Echo-tomography
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/0033—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room
- A61B5/0035—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room adapted for acquisition of images from more than one imaging mode, e.g. combining MRI and optical tomography
-
- 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/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/74—Details of notification to user or communication with user or patient ; user input means
- A61B5/742—Details of notification to user or communication with user or patient ; user input means using visual displays
- A61B5/7425—Displaying combinations of multiple images regardless of image source, e.g. displaying a reference anatomical image with a live image
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/74—Details of notification to user or communication with user or patient ; user input means
- A61B5/742—Details of notification to user or communication with user or patient ; user input means using visual displays
- A61B5/7435—Displaying user selection data, e.g. icons in a graphical user interface
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/42—Details of probe positioning or probe attachment to the patient
- A61B8/4272—Details of probe positioning or probe attachment to the patient involving the acoustic interface between the transducer and the tissue
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/50—Image enhancement or restoration using two or more images, e.g. averaging or subtraction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/90—Dynamic range modification of images or parts thereof
- G06T5/94—Dynamic range modification of images or parts thereof based on local image properties, e.g. for local contrast enhancement
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/74—Details of notification to user or communication with user or patient ; user input means
- A61B5/7475—User input or interface means, e.g. keyboard, pointing device, joystick
- A61B5/748—Selection of a region of interest, e.g. using a graphics tablet
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10132—Ultrasound image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20092—Interactive image processing based on input by user
- G06T2207/20104—Interactive definition of region of interest [ROI]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20212—Image combination
- G06T2207/20221—Image fusion; Image merging
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30101—Blood vessel; Artery; Vein; Vascular
- G06T2207/30104—Vascular flow; Blood flow; Perfusion
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Surgery (AREA)
- Public Health (AREA)
- Pathology (AREA)
- Veterinary Medicine (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Acoustics & Sound (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Human Computer Interaction (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physiology (AREA)
- Psychiatry (AREA)
- Signal Processing (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
本公开涉及被检体信息获取装置和被检体信息获取方法。一种被检体信息获取装置包括:光声图像获取单元,其被配置为产生与被检体内的光学特性相关的第一图像;超声图像获取单元,其被配置为产生与被检体内的声学特性相关的第二图像;感兴趣区域指定单元,其被配置为接收关于第一图像对感兴趣区域的指定;图像处理单元,其被配置为对第一图像在感兴趣区域内部和感兴趣区域外部分别使用不同的图像处理参数执行图像处理;以及图像合成单元,其被配置为叠加并合成已经经过了图像处理的第一图像和第二图像。
Description
本申请是基于申请号为201310719112.8、申请日为2013年12月24日、发明名称为“被检体信息获取装置和被检体信息获取方法”的专利申请的分案申请。
技术领域
本发明涉及在被检体信息获取装置中显示图像数据的技术。
背景技术
以前关于使用光对被检体的断层图像进行成像的技术提出了各种提议。在这样的提议之中,存在使用光声层析成像(PAT)技术的光声断层图像成像装置(以下,称之为“光声成像装置”)。
光声成像装置将诸如脉冲激光束的测量光发射到被检体,接收当测量光被被检体中的活体组织吸收时所产生的声波,并对这些声波执行分析处理以便使与活体内部的光学特性相关的信息(机能信息)可视化。
动脉血中所含有的大量氧合血红蛋白和静脉血中所含有的大量还原血红蛋白吸收激光束并产生声波,但是激光束的吸收率根据波长而不同。例如,氧合血红蛋白具有高的吸收805nm或更小的光的速率,还原血红蛋白具有高的吸收805nm或更大的光的速率。
因此,通过发射不同波长的激光束并测量相应声波,可以使氧合血红蛋白和还原血红蛋白的分布状态可视化,并且通过分析所获得的信息来计算血红蛋白的数量或氧饱和度。因为这种类型的机能信息可以用作与肿瘤细胞附近的血管化相关的信息,所以已知光声成像装置对于皮肤癌和乳腺癌的诊断特别有效。
同时,超声成像装置又被已知为可以与光声成像装置一样无暴露地、无创伤性地执行成像的图像诊断装置。超声成像装置将超声波发射到活体,并接收作为在被检体内传播的超声波反射出活体组织中的具有不同声学特性(声学阻抗)的组织界面的结果而产生的声波。另外,通过对所接收的声波执行分析处理,使与作为被检体的活体内部的声学特性相关的信息(形状信息)可视化。可视化的形状信息的独特之处在于,它可以提供活体组织的形状的指示。
尽管光声成像装置可以获取机能信息,但是仅通过机能信息,难以确定这样的机能信息是从活体组织的哪个部分产生的。因此,提出了将超声成像单元合并在光声成像装置内、同时获取形状信息的技术。例如,日本专利申请公开No.2005-21580公开了一种活体信息成像装置,该活体信息成像装置获取光声图像和超声图像这两者,并通过叠加这两个图像数据或者彼此相邻地显示这两个图像数据来便利于被检体内的位置的理解。
当对机能信息进行成像和显示时,存在这样的问题,即,感兴趣区域(ROI)内部的对比度由于ROI外部的无用图像分量(从与皮肤的边界产生的强噪声和伪像)而变得不足。
例如,机能信息之中的作为无用图像分量的来自皮肤表面的强反射波和基于多次反射的伪像有时变为等于或大于ROI内部的强信号。当对机能信息进行成像时,因为根据输入信号分配像素值,所以当基于整个图像的信号水平决定像素值时,存在ROI内部的对比度变得不足的情况。另外,当叠加并显示具有两种不同类型的特性的图像信息(诸如机能信息和形状信息)时,如果在ROI内部没有获得足够的对比度,则变得难以区分这两个图像。
发明内容
鉴于前述问题,本发明的一个目的是提供一种能够产生在感兴趣区域中保证有足够对比度的光声图像的被检体信息获取装置。
本发明在其一方面提供了一种被检体信息获取装置,该被检体信息获取装置包括:光声图像获取单元,其被配置为将测量光发射到被检体,接收在被检体中产生的光声波,并基于光声波产生使与被检体内的光学特性相关的信息可视化的第一图像;超声图像获取单元,其被配置为将超声波发送到被检体,接收在被检体中反射的超声回声,并基于超声回声产生使与被检体内的声学特性相关的信息可视化的第二图像;感兴趣区域指定单元,其被配置为接收关于第一图像对感兴趣区域的指定;图像处理单元,其被配置为对第一图像在所指定的感兴趣区域内部和所指定的感兴趣区域外部分别使用不同的图像处理参数执行图像处理;以及图像合成单元,其被配置为叠加并合成已经经过了图像处理的第一图像和第二图像。
本发明在其另一方面提供了一种被检体信息获取装置,该被检体信息获取装置包括:光声图像获取单元,其被配置为将不同波长的测量光发射到被检体,对于每个波长,接收在被检体中产生的光声波,并且对于每个波长,基于光声波产生使与被检体内的光学特性相关的信息可视化的图像;感兴趣区域指定单元,其被配置为接收感兴趣区域的指定;图像处理单元,其被配置为对多个图像中的每个在感兴趣区域内部和外部分别使用不同的图像处理参数执行图像处理;以及图像合成单元,被配置为叠加并合成已经经过了图像处理的所述多个图像。
根据本发明,可以提供一种能够产生在感兴趣区域内保证有足够对比度的光声图像的被检体信息获取装置。
从以下参照附图对示例性实施例的描述,本发明的进一步的特征将变得清楚。
附图说明
图1是示出根据第一实施例的光声成像装置的总体配置的示图;
图2是示出根据第一实施例的光声成像装置的修改例子的示图;
图3是示出根据第一实施例的ROI指定模式的GUI显示例子的示图;
图4是示出根据第一实施例的叠加图像显示模式的GUI显示例子的示图;
图5是示出ROI内部部分的光声图像的例子的示图;
图6是示出ROI外部部分的光声图像的例子的示图;
图7是示出超声图像的例子的示图;
图8是示出叠加图像的例子的示图;
图9A和图9B是示出第一实施例中的控制流程图的示图;
图10是示出根据第二实施例的光声成像装置的总体配置的示图;和
图11是示出根据第二实施例的GUI显示例子的示图。
具体实施方式
现在参照附图详细解释本发明的实施例。指出,作为一般规则,相同的组成元件被给予相同的标号,并且省去它们的多余的解释。
(第一实施例)
<系统配置>
首先,参照图1解释根据第一实施例的光声成像装置的配置。根据本发明的第一实施例的光声成像装置是用于对作为被检体的活体的信息进行成像以用于恶性肿瘤和血管疾病的诊断或化疗的随访的装置。活体的信息例如是基于光的照射而产生的声波(以下称之为“光声波”)的产生源分布、活体中的初始声压分布、或从其推导的光能吸收密度分布。换句话讲,根据第一实施例的光声成像装置也可被称为被检体信息获取装置。
根据第一实施例的光声成像装置具有将测量光发射到被检体并分析光声波以使与光学特性相关的机能信息可视化或成像的光声成像功能。而且,光声成像装置还具有将超声波发射到被检体并分析在被检体内部反射的超声波(以下称之为“超声回声”)以使与声学特性相关的形状信息成像的超声成像功能。而且,光声成像装置还具有叠加并合成(以下简单地称之为“叠加”)所获得的图像并显示叠加图像的功能。在随后的解释中,通过光声成像而获得的图像被称为光声图像,通过超声成像而获得的图像被称为超声图像。
根据本发明的第一实施例的光声成像装置1由光声图像获取单元10、超声图像获取单元20、图像产生单元30、图像显示单元40、操作输入单元50和控制器单元60构成。指出,标号2表示作为被检体的活体的一部分。现在解释显示图像的方法的概要,同时解释构成根据第一实施例的光声成像装置的各个单元。
<<光声图像获取单元10>>
光声图像获取单元10是用于通过光声成像产生光声图像的单元。例如,可以获取表示氧饱和度的图像,氧饱和度是活体的机能信息。光声图像获取单元10由光照射控制单元11、光照射单元12、光声信号测量单元13、光声信号处理单元14、光声图像累积单元15和超声探测器16构成。
光照射单元12是用于产生将发射到作为被检体的活体的近红外测量光的单元,光照射控制单元11是用于控制光照射单元12的单元。
优选地,从光照射单元12产生被构成活体的组成部分之中的特定组成部分吸收的特定波长的光。具体地讲,优选使用的是能够产生几纳秒至几百纳秒的数量级的脉冲光的脉冲光源。尽管光源优选为用于产生激光束的光源,但是还可以使用发光二极管来代替激光束源。当使用激光器时,可以使用诸如固态激光器、气体激光器、染料激光器或半导体激光器的各种激光器。
而且,激光束的波长优选地在活体内吸收率低的700nm至1100nm的区域中。然而,当获得相对接近活体表面的活体组织的光学特性值分布时,还可以使用比前述波长区域的范围宽的波长区域;例如,还可以使用400nm至1600nm的波长区域。在前述范围内的光中,特定波长可以基于将被测量的组成部分来选择。
超声探测器16是用于检测在作为被检体的活体内产生的光声波、并将所检测的光声波换能为模拟电信号的单元。因为从活体产生的光声波是100KHz至100MHz的超声波,所以用作超声探测器16的是能够接收前述频段的超声换能器。具体地讲,所使用的是利用压电陶瓷(PZT)的感测元件或麦克风类型的电容感测元件。
而且,还可以使用电容型电容微加工超声换能器(CMUT)、使用磁性膜的磁性MUT(MMUT)、或使用压电薄膜的压电MUT(PMUT)。
指出,任何类型的感测元件可以用作超声探测器16,只要它可以将声波信号换能为电信号即可。
被超声探测器16换能的模拟电信号被光声信号测量单元13放大,被转换为数字信号、然后被光声信号处理单元14转换为图像数据。在本发明中,该图像数据是第一图像。所产生的图像数据被存储在光声图像累积单元15中。
<<超声图像获取单元20>>
超声图像获取单元20是用于通过超声成像获取活体的形状信息并产生超声图像的单元。超声图像可以是B模式图像、或者基于多普勒方法或弹性成像而产生的图像。超声图像获取单元20由超声发送控制单元21、超声探测器22、超声信号测量单元23、信号处理单元24、超声图像累积单元25和超声发送/接收开关26构成。
超声探测器22与超声探测器16一样是包括感测元件的探测器,并且可以将超声波束发送到被检体。
超声发送控制单元21是用于产生将施加于构建到超声探测器22中的各个声学元件的信号、并控制将发送的超声波的频率和声压的单元。
因为超声信号测量单元23、信号处理单元24和超声图像累积单元25分别是执行与光声信号测量单元13、光声信号处理单元14和光声图像累积单元15类似处理的单元,所以省去它们的详细解释。唯一的不同之处是将被处理的信号是在被检体内部产生的光声波、还是其中超声波在被检体内部反射的超声回声。而且,在本发明中,超声图像获取单元20产生的图像数据是第二图像。
超声发送/接收开关26是由超声发送控制单元21控制的开关,并且是用于切换与超声探测器22的超声波发送和接收的单元。超声发送控制单元21在将超声发送/接收开关26切换到“发送”的状态下发送超声波,并且通过在给定时间过去之后切换到“接收”,接收从被检体内部返回的超声回声。
<<图像产生单元30>>
图像产生单元30是用于对在光声图像累积单元15中累积的光声图像执行图像处理的单元。而且,图像产生单元30还是用于执行叠加处理后的光声图像和在超声图像累积单元25中累积的超声图像的处理、并产生将呈现给用户的图像的单元。
图像产生单元30由光声图像处理单元31和图像合成单元32构成。
光声图像处理单元31是用于对在光声图像累积单元15中累积的光声图像执行图像处理的单元。稍后将解释处理内容的细节。
图像合成单元32是用于叠加已经经受了光声图像处理单元31的图像处理的光声图像和在超声图像累积单元25中累积的超声图像并产生单张图像的单元。在随后的解释中,图像合成单元32产生的图像被称为叠加图像。
指出,图像产生单元30还具有产生用于执行图像处理的操作用GUI并将所产生的操作用GUI输出到图像显示单元40的功能。
<<图像显示单元40>>
图像显示单元40是用于向用户呈现图像产生单元30产生的操作用GUI的单元。图像产生单元30产生的叠加图像与操作用GUI一起呈现给用户。
<<操作输入单元50>>
操作输入单元50是用于从用户接收操作输入的单元。用于操作输入的单元可以是诸如鼠标或手写平板的定点设备或键盘等。而且,操作输入单元50还可以是与图像显示单元40整体地形成的、诸如触摸板或触摸屏的设备。
<<控制器单元60>>
控制器单元60是由CPU、DRAM、非易失性存储器、控制部分等(所有这些都未示出)构成的计算机。作为存储在非易失性存储器中的程序被CPU执行的结果,光声成像装置1的各个模块被控制。尽管在本实施例中控制器单元是计算机,但是控制器单元还可以是专门设计的硬件。
<<超声探测器的布置例子>>
图1示出光声图像获取单元10所使用的超声探测器16和超声图像获取单元20所使用的超声探测器22互相独立的例子。然而,因为用于光声成像的超声探测器和用于超声成像的超声探测器互相接收同一频段的超声波,所以它们可以被共享。
因此,还可以省去超声探测器16,并且基于分时控制使光声图像获取单元10和超声图像获取单元20共享超声探测器22。图2是示出共享超声探测器22的例子的系统配置图。指出,因为光声信号测量单元13也可以与超声信号测量单元23被共享,所以在图2中省去它。
<操作用GUI>
现在解释用于将指令给予光声成像装置并显示图像的操作用GUI。图3示出图像产生单元30产生的、在图像显示单元40上显示的操作用GUI的例子。这里,解释配置操作用GUI的各个界面。
<<用于显示图像的界面>>
图像显示区域41是显示光声图像或叠加图像的区域。在本实施例中,假设超声图像是具有40mm的宽度×30mm的高度(深度)以及每像素12位灰度分级(gradation)(4096灰度分级)的B模式图像。而且,还假设光声图像类似地是具有40mm的宽度×30mm的高度(深度)以及每像素12位灰度分级(4096灰度分级)的图像。
指出,尽管超声图像和光声图像都是灰度级图像,但是为了提高可见性,对光声图像进行将不同颜色分配给各个像素的每个像素值(即,亮度值)的彩色显示。例如,通过将红色分配给高亮度侧、将黄绿色分配给中间值并将蓝色分配给低亮度侧来显示光声图像。稍后将解释分配颜色的方法。指出,在随后的解释中,通过使用术语“亮度值”将光声图像解释为上色之前的灰度图像。
<<用于调整光声图像的对比度的界面>>
亮度值指定界面42是用于执行所获取的光声图像的对比度调整的界面。具体地讲,亮度值指定界面42是用于指定当调整对比度时的亮度值的上限/下限的界面。低端表示最低亮度,高端表示最高亮度。
这里,基于先前已经对于光声图像指定了ROI的假设,解释这个界面。稍后将解释指定ROI的方法。
在亮度值指定界面42上重叠并显示了两种类型的滑杆。一个是亮度值上限滑杆421,另一个是亮度值下限滑杆422。在操作用GUI的初始屏幕上,各个滑杆分别布置在表示存在于光声图像的ROI内部的像素(以下称之为“ROI内部的像素”)之中的最高亮度的位置处和表示这些像素之中的最低亮度的位置处。
使用通过各个滑杆指定的亮度值的范围来重新分配光声图像的所有像素的亮度值。例如,考虑ROI中所包含的像素的最低亮度值为n并且最高亮度值为m的情况。亮度值下限滑杆422设置在表示亮度值n的位置处,亮度值上限滑杆421设置在表示亮度值m的位置处。另外,亮度值n至m被重新分配给最小亮度值至最大亮度值。最小亮度值或者最大亮度值被分配给亮度值为n或更小或者m或更大的像素。换句话讲,对整个光声图像执行最大程度地强调ROI内部的对比度的图像处理。
指出,可手动地将各个滑杆的位置变到任意位置。当滑杆的位置改变时,再一次执行对比度调整;也就是说,基于新的位置分配亮度值。
这里,现在更详细地解释当滑杆的位置改变时对比度的变化。
例如,假设用户基于使用鼠标的拖动操作将滑杆421从初始值向上移动。如上所述,通过执行将用滑杆指定的亮度值的范围重新分配给最小亮度值至最大亮度值的处理来执行对比度调整。因此,执行使整个图像的对比度弱化的处理。
相反,假设将滑杆421从初始值向下移动。因为在前述情况下也执行类似的处理,所以执行强调整个图像的对比度的图像处理。因为最大亮度值被分配给亮度值大于用滑杆421指定的亮度值的所有像素,所以显示器将变得饱和。
接着,考虑将滑杆422从初始值向下移动的情况。在前述情况下,与向上移动滑杆421的情况一样,执行使整个图像的对比度弱化的图像处理。
相反,假设从初始值向上移动滑杆422。在前述情况下,与向下移动滑杆421的情况一样,执行强调整个图像的对比度的图像处理。最小亮度值被分配给亮度值小于用滑杆422指定的亮度值的所有像素。
如上所述,可通过设置在亮度值指定界面42上的两个滑杆来调整整个图像的对比度,以使得ROI内部的可见性变为最高。图像产生单元30产生的、通过操作输入单元50操作的亮度值指定界面42构成本发明中的像素值范围指定单元。
<<用于调整光声图像的ROI外部的不透明度的界面>>
ROI外部透明度指定界面43是用于调整所获取的光声图像的ROI外部的像素的不透明度的界面。关于ROI外部透明度指定界面43,下侧表示低不透明度(也就是说,更透明),上侧表示高不透明度(也就是说,更不透明)。
在ROI外部透明度指定界面43上叠加并显示了一种类型的滑杆(ROI外部不透明度指定滑杆431)。滑杆431是用于指定被指定为ROI的区域外部的区域的像素(以下称之为“ROI外部的像素”)的不透明度的滑杆。在初始画面上,滑杆431设置在预先设置的值(例如,50%的不透明度)处。
ROI外部的像素的不透明度被设置为使得它变为通过滑杆指示的值。例如,当滑杆在指示50%的位置处时,对光声图像的ROI外部的像素执行将不透明度设置为50%的图像处理。
指出,滑杆431可以用于通过使用鼠标的拖动操作任意地改变值。
这里,考虑从初始值向下拖动滑杆431的情况。在前述情况下,执行降低ROI外部的不透明度的图像处理。换句话讲,当叠加图像时,提高ROI外部的像素的透射率,并且背景图像(在本实施例中,超声图像)变得可见性更高。
而且,当从初始值向上拖动滑杆431时,执行提高ROI外部的不透明度的图像处理。换句话讲,当叠加图像时,降低ROI外部的像素的透射率,并且背景图像变得可见性更低。
<<用于指定光声图像的ROI的界面>>
现在解释用于指定光声图像的ROI的用户界面。
ROI指定单元45是用于指定光声图像的ROI的界面。ROI指定单元45由ROI指定按钮451和ROI半径显示单元452构成。通过用鼠标点击ROI指定按钮451,模式变为ROI指定模式。而且,通过再次点击ROI指定按钮451,模式变为叠加图像显示模式。
首先解释ROI指定模式。ROI指定模式是使得能够实现指定ROI的操作的模式。图3是ROI指定模式的画面显示例子。
在ROI指定模式下,显示在图像显示区域41上的是光声图像和作为用于显示ROI范围的图形的ROI显示46。ROI显示46被显示为使用与其他UI中所使用的颜色不同的颜色(例如,浅紫色)的虚线圆圈。可通过用鼠标拖动ROI显示46来移动它。
而且,在ROI指定模式下,总共在八个位置(即,表示ROI的圆圈的顶部、底部、左侧、右侧、左上、左下、右上和右下)处显示ROI半径指定手柄461。用户可通过使用鼠标拖动ROI半径指定手柄461中的一个来改变ROI半径。
这里,基于拖动操作而变化的ROI半径同时还显示在ROI半径显示单元452上。而且,相反,还可以通过直接将ROI半径的数值输入到ROI半径显示单元452中来指定ROI半径。在前述情况下,输入的ROI半径被反映,并且ROI显示46被更新。图像产生单元30产生的、通过操作输入单元50操作的ROI指定单元45和ROI显示46构成本发明中的感兴趣区域指定单元。
现在解释叠加图像显示模式。叠加图像显示模式是在图像显示区域41上显示图像处理之后的光声图像(即,在对比度和不透明度已经被调整之后的光声图像)和超声图像的叠加图像的模式。图4是叠加图像显示模式下的画面显示例子。指出,为了使可见性更好,图4仅示出了光声图像。在叠加图像显示模式下显示表示ROI的圆圈时,不显示ROI半径指定手柄461,并且不可以移动ROI或改变半径。
<<其他UI>>
现在参照图4解释其他UI的例子。
标号44示出显示表示超声图像的亮度值的标度的区域。通过将最大亮度值分配给白色来显示最大亮度值,通过将中间值分配给灰色来显示该中间值,通过将最小亮度值分配给黑色来显示最小亮度值。
标号47示出用于命令光声图像获取单元10和超声图像获取单元20分别获取图像的图像获取按钮。
标号48示出用于命令光声成像装置1结束它的操作的按钮。
标号49示出用于显示关于光声图像的ROI内部和外部的像素的亮度值直方图的直方图显示区域。这里,ROI内部的像素的亮度值直方图用黑色显示,ROI外部的像素的亮度值直方图用灰色显示。
<<图像处理操作>>
现在参照图4解释图像产生单元30对光声图像执行的图像处理的细节。
图像产生单元30首先获取关于所指定的ROI的信息,然后产生作为ROI内部的像素的亮度值直方图(频率分布)的ROI内部直方图491、以及作为ROI外部的像素的亮度值直方图的ROI外部直方图493。
图像产生单元30从ROI内部直方图491提取ROI内部的像素的最大亮度值和最小亮度值,将最大亮度值设置为滑杆421的值,并将最小亮度值设置为滑杆422的值。在随后的解释中,滑杆421所指示的亮度值被表示为ROImax,滑杆422所指示的亮度值被表示为ROImin。
指出,在亮度值指定界面42所表示的区域之中,在滑杆421上方的区域中以及在滑杆422下方的区域中显示大意是具有亮度值的像素不存在于ROI内部的消息。相应区域例如被用灰色填充。
随后,对于光声图像中的所有像素,使用ROImax和ROImin重新分配亮度值。具体地讲,具有ROImin或更小的值的像素的亮度值被分配最低亮度值,具有ROImax或更大的值的像素的亮度值被分配最高亮度值,并且中间值通过线性插值被分配。指出,还可以通过诸如直方图平坦化或伽马校正的方法来分配亮度值。
随后,执行用于改进图像的可见性的颜色分配。
当亮度值被重新分配时,光声图像处理单元31对于光声图像用暗红色取代具有最大亮度值的像素,并用深蓝色取代具有最低亮度值的像素。关于中间亮度值,可以分配任意颜色显示。
示出颜色分配方法的例子。考虑这样的情况,即,RGB的各个颜色和用8位显示不透明度α的颜色坐标被定义为(R,G,B,α),并且从最低亮度值开始按顺序分配深蓝色、蓝色、浅蓝色、绿色、黄色、橙色、红色和暗红色。可以如下表示各个颜色的颜色坐标:
深蓝色:(0,0,128,255),蓝色:(0,0,255,255)
浅蓝色:(0,255,255,255),绿色(0,255,0,255)
黄色:(255,255,0,255),橙色:(255,128,0,255)
红色:(255,0,0,255),暗红色:(128,0,0,255)。
换句话讲,只有B坐标在深蓝色与蓝色之间在128至255的范围内变化,只有G坐标在蓝色与浅蓝色之间在0至255的范围内变化,并且只有B坐标在浅蓝色与绿色之间在255至0的范围内变化。而且,只有R坐标在绿色与黄色之间在0至255的范围内变化,并且只有G坐标在黄色、橙色与红色之间在255至0的范围内变化。只有R坐标在红色与暗红色之间在255至128的范围内变化。换句话讲,存在1280种模式的颜色坐标。
在本实施例中,尽管光声图像是12位灰度分级(4096灰度分级)的,但是因为存在1280种模式的取代颜色坐标,所以基于对比度调整,原始亮度值被1280灰度分级取代。通过对原始亮度值Vpix进行对比度调整并且用1280灰度分级取代而获得的值Vroi将如公式1中所示。
(1)当Vpix≥ROImax时,Vroi=1280
(2)当ROImin<Vpix<ROImax时,Vroi=1280×(Vpix-ROImin)/(4096×(ROImax-ROImin))
(3)当Vpix≤ROImin时,Vroi=0
...公式1(0≤Vroi≤1280)
首先解释通过使用所确定的Vroi确定ROI内部的像素的像素值的方法。当所确定的Vroi被颜色坐标取代时,实现以下。
(1)当0≤Vroi<127时,(R,G,B,α)=(0,0,Vroi+128,255)
(2)当127≤Vroi<382时,(R,G,B,α)=(0,Vroi-127,255,255)
(3)当382≤Vroi<637时,(R,G,B,α)=(0,255,637-Vroi,255)
(4)当637≤Vroi<892时,(R,G,B,α)=(Vroi-637,255,0,255)
(5)当892≤Vroi<1147时,(R,G,B,α)=(0,1147-Vroi,255,255)
(6)当1147≤Vroi≤1280时,(R,G,B,α)=(1402-Vroi,0,0,255)
...公式2
因此,ROI内部的所有像素都可以在调整对比度之后被转换为彩色显示。指出,可以在亮度值指定界面42上将原始亮度值和所分配的颜色的对应关系显示为颜色标度。
光声图像在ROI外部的各个像素的像素值也基于与ROI内部的像素相同的方法来确定。然而,因为无用噪声分量和伪像常常存在于ROI外部,所以希望的是另外对ROI外部的像素执行降低可见性的处理。
因此,除了对ROI内部的像素执行的对比度调整之外,还通过对ROI外部的像素降低不透明度来降低可见性。这里,不透明度α被设置,并且不透明度α被设置给ROI外部的所有像素。不透明度α是由滑杆431指定的值。初始值为50%(也就是说,α=128)。
这里,当所指定的不透明度为αext时,ROI外部的像素的颜色坐标将如公式3中所示。公式3与公式2相比的不同之处仅在于不透明度的指定。
(1)当0≤Vroi<127时,(R,G,B,α)=(0,0,Vroi+128,αext)
(2)当127≤Vroi<382时,(R,G,B,α)=(0,Vroi-127,255,αext)
(3)当382≤Vroi<637时,(R,G,B,α)=(0,255,637-Vroi,αext)
(4)当637≤Vroi<892时,(R,G,B,α)=(Vroi-637,255,0,αext)
(5)当892≤Vroi<1147时,(R,G,B,α)=(0,1147-Vroi,255,αext)
(6)当1147≤Vroi≤1280时,(R,G,B,α)=(1402-Vroi,0,0,αext)
...公式3
图5示出应用公式2并提高ROI内部的像素的可见性的光声图像的例子。而且,图6示出应用公式3并降低ROI外部的像素的可见性的光声图像的例子。在这个例子中,尽管为了便利于解释而分别在图5和图6中示出了图像,但是作为图像处理的结果而产生的光声图像是单个光声图像。
而且,图7示出了超声图像的例子,图8示出了叠加并显示已经经过图像处理的光声图像和超声图像的例子。
如上所述,根据第一实施例的光声成像装置可以执行这样的图像处理,该图像处理基于对比度调整来提高ROI内部的像素的可见性,并通过另外执行不透明度调整来降低ROI外部的像素的可见性。
<处理流程图>
现在参照图9A和图9B解释根据第一实施例的光声成像装置产生叠加图像的处理,图9A和图9B是处理流程图。
在步骤S1中,在光声成像装置1的电源开启并且各种初始化被执行之后,图像产生单元30在图像显示单元40上显示图3中所示的操作用GUI。
在步骤S2中,确定图像获取按钮47是否被点击。当点击事件发生时,例行程序前进到步骤S3,当点击事件还没有发生时,所述处理等待事件发生。
在步骤S3中,光声图像获取单元10获取光声图像,超声图像获取单元20获取超声图像。光声图像被存储在光声图像累积单元15中,超声图像被存储在超声图像累积单元25中。
在步骤S4中,光声图像处理单元31设置操作参数中的初始值。操作参数是由当前模式(叠加图像显示模式或ROI指定模式)、ROI的中心点坐标和ROI半径构成的信息。例如,模式被设置为叠加图像显示模式,并且ROI的中心点坐标被设置为图像显示区域的中心。而且,ROI半径被设置为5mm。
在步骤S5中,光声图像处理单元31获取操作参数。模式、ROI的中心点坐标和ROI半径从而被阐述,并且ROI被识别。
在步骤S6中,光声图像处理单元31使用在步骤S5中识别的ROI信息,并产生ROI内部的像素的直方图和ROI外部的像素的直方图。所产生的直方图显示在用标号49所示的区域中。
而且,滑杆421、422的位置分别被设置为ROI内部的像素的最大亮度值和最小亮度值。而且,当滑杆421、422已经在所设置的ROI中被手动移动时,省去这个处理。
随后,用滑杆421、422所指定的亮度值替代ROImax和ROImin。而且,用滑杆431所指定的不透明度替代αext。如果滑杆431从未被操作,则αext为128。
在步骤S7中,对在步骤S3中获取的光声图像执行图像处理。具体地讲,ROI的中心点坐标和ROI半径用于确定构成在步骤S3中获取的光声图像的像素是在ROI内部、还是在ROI外部,并且公式1用于调整像素的亮度值,公式2和3用于分配颜色。因此,获得经过图像处理之后的光声图像。所获得的图像被临时存储。
而且,在步骤S7中,在亮度值指定界面42上将基于公式1和2分配给各个亮度值的颜色显示为颜色标度。不存在于ROI内部的亮度值用灰色显示。
在步骤S8中,图像合成单元32叠加已经在步骤S7中经过图像处理的光声图像与在步骤S3中获取的超声图像,并在图像显示区域41上将叠加图像与ROI显示46一起显示。这里,当模式是ROI指定模式时,显示ROI半径指定手柄461。当模式是叠加图像显示模式时,不显示ROI半径指定手柄。
步骤S9是等待对于构成操作用GUI的各个部分的诸如点击或拖动的事件的发生的步骤。一旦事件发生,例行程序就前进到图9B的步骤S10。
步骤S10是确定发生的事件的类型的步骤。现在解释各个事件。
当结束按钮48被点击(S11)时,例行程序前进到步骤S12,并且关闭光声成像装置1以结束该处理。
当ROI指定按钮451被点击(S20)时,例行程序前进到步骤S21,并且通过更新指示模式的操作参数来切换模式。在当前模式是叠加图像显示模式时,模式切换到ROI指定模式,在当前模式是ROI指定模式时,模式切换到叠加图像显示模式。指出,仅在当前模式是ROI指定模式时,才使得能够拖动ROI显示46和ROI半径指定手柄461并且将数值输入到ROI半径显示单元452中。当这个处理结束时,例行程序前进到步骤S5。
当ROI半径指定手柄461被拖动(S30)时,例行程序前进到步骤S32,并且改变ROI半径。具体地讲,从当拖动完成时的手柄坐标和ROI的中心点坐标计算ROI半径,并更新指示ROI半径的操作参数。
而且,在ROI半径显示单元452中反映所计算的ROI半径,并更新ROI显示46。当这个处理结束时,例行程序前进到步骤S5。
当数值被输入到ROI半径显示单元452中(S31)时,所述处理还前进到步骤S32,并且改变ROI半径。具体地讲,用作为ROI半径的值的输入数值更新指示ROI半径的操作参数。而且,根据新的ROI半径更新ROI显示46。当这个处理结束时,例行程序前进到步骤S5。
当ROI显示46被拖动(S40)时,例行程序前进到步骤S41,并且移动ROI。具体地讲,获取当拖动完成时ROI显示46的中心点坐标,并使用所获取的中心点坐标来更新指示ROI的中心点的操作参数。而且,根据中心点坐标更新ROI显示46。当这个处理结束时,例行程序前进到步骤S5。
当亮度值上限滑杆421被拖动(S50)时,或者当亮度值下限滑杆422被拖动(S51)时,例行程序前进到步骤S52,并且更新各个滑杆的位置。当这个处理结束时,例行程序前进到步骤S5。
而且,当ROI外部不透明度指定滑杆431被拖动(S53)时,例行程序前进到步骤S54,并且更新滑杆431的位置。当这个处理结束时,例行程序前进到步骤S5。
当各个滑杆被拖动时,在步骤S6中重新设置ROImax、ROImin和αext,并且在步骤S7中使用所设置的值来执行图像处理。
在步骤S8中,当事件没有发生或者除了上述那些事件之外的事件发生时,不执行所述处理,并且例行程序等待。
如以上所解释的,在第一实施例中,在叠加并显示光声图像和超声图像的光声成像装置中,在感兴趣区域内部和感兴趣区域外部分别通过使用不同的图像处理参数执行图像处理。因此,可以改进ROI内部的信号的可见性并且使ROI外部的信号(噪声、伪像)变得不明显。
指出,当然,将被分配给光声图像的各个像素的颜色可以不同于所示的颜色。例如,最大值侧可以被分配白色,最小值侧可以被分配黑色,以实现黑白显示,或者其他颜色显示可以被分配。
(第二实施例)
第二实施例是将多个波长的测量光发射到被检体、获取多个光声图像并对各个光声图像执行图像处理的实施例。
例如,分别对通过发射作为第一波长的750nm附近的激光束而获取的第一光声图像和通过发射作为第二波长的830nm附近的激光束而获取的第二光声图像执行图像处理,并且叠加并显示所获得的两个图像。对各个图像执行的图像处理的内容与第一实施例相同。
图10是示出根据第二实施例的光声成像装置的总体配置的示图。
尽管光照射单元18类似于根据第一实施例的光照射单元12,但是它的不同之处在于它可以发射两个不同波长的激光束。而且,尽管光照射控制单元17类似于根据第一实施例的光照射控制单元11,但是它的不同之处在于它可以向光照射单元18发出波长切换命令。
而且,光声信号处理单元14与第一实施例的不同之处在于,在第一光声图像累积单元15中累积通过发射第一波长而获得的第一光声图像,并且在第二光声图像累积单元19中累积通过发射第二波长而获得的第二光声图像。
而且,根据第二实施例的光声成像装置1不包括超声图像获取单元20。因为其他单元与第一实施例相同,所以省去它们的解释。
图11示出根据第二实施例的光声成像装置中的操作用GUI显示的例子。这里,解释与第一实施例中的操作用GUI显示的不同之处。第二实施例中的操作用GUI显示与第一实施例的不同之处在于,分别包括两个直方图显示区域、两个亮度值指定界面和两个ROI外部透明度指定界面。各个区域和界面对应于第一光声图像和第二光声图像。
直方图显示区域49是用于显示第一光声图像的ROI内部和ROI外部的亮度值直方图的直方图显示区域。而且,直方图显示区域4a是用于显示第二光声图像的ROI内部和ROI外部的亮度值直方图的直方图显示区域。
而且,亮度值指定界面42是用于调整第一光声图像的亮度值的界面,亮度值指定界面4b是用于调整第二光声图像的亮度值的界面。
而且,ROI外部透明度指定界面43是用于调整第一光声图像的ROI外部的像素的不透明度的界面,ROI外部透明度指定界面4c是用于调整第二光声图像的ROI外部的像素的不透明度的界面。因为各个操作与第一实施例相同,所以省去它们的解释。
在第一实施例中,通过基于像素的亮度值分配不同颜色来执行颜色显示,但是在第二实施例中,因为光声图像被叠加,所以如果采用相同方法,则相同颜色将被分配给不同图像,并且图像的区分将变得困难。
因此,在第二实施例中,在第一光声图像和第二光声图像中使用不同色调进行上色。具体地讲,第一光声图像基于红色,并且通过提高高亮度侧的明亮度并降低低亮度侧的明亮度来分配颜色。此外,第二光声图像基于蓝色,并且通过提高高亮度侧的明亮度并降低低亮度侧的明亮度来分配颜色。从而可以区分这两个图像。
现在解释将颜色分配给像素的方法。
首先,从第一光声图像的ROI内部的直方图提取最大值ROI1max和最小值ROI1min,并将淡红色(255,191,191,255)分配给ROI1max,将暗红色(128,0,0,255)分配给ROI1min。
在暗红色与淡红色之前,R坐标首先在128至255的范围内变化,随后G坐标和B坐标同时在0至191的范围内变化。换句话讲,存在分配给第一光声图像的320种模式的颜色坐标。
类似地,从第二光声图像的ROI内部的直方图提取最大值ROI2max和最小值ROI2min,并将淡蓝色(191,191,255,255)分配给ROI2max,将深蓝色(0,0,128,255)分配给ROI2min。
在深蓝色与浅蓝色之间,B坐标首先在128至255的范围内变化,随后R坐标和G坐标同时在0至191的范围内变化。换句话讲,类似地,存在分配给第二光声图像的320种模式的颜色坐标。
在第二实施例中,因为存在320种模式的取代颜色坐标,所以基于对比度调整,原始亮度值被320灰度分级替代。通过对第一光声图像的亮度值V1pix进行对比度调整并且用320灰度分级取代而获得的值V1roi将如公式4中所示。
(1)当V1pix≥ROI1max时,V1roi=319
(2)当ROI1min<V1pix<ROI1max时,V1roi=319×(V1pix-ROI1min)/(4096×(ROI1max-ROI1min))
(3)当V1pix≤ROI1min时,V1roi=0
...公式4(0≤V1roi≤319)
当V1roi被颜色坐标取代时,实现以下。
(1)当0≤V1roi<128时,(R,G,B,α)=(V1roi+128,0,0,αext)
(2)当128≤V1roi≤319时,(R,G,B,α)=(255,V1roi-128,V1roi-128,αext)
...公式5
然而,当目标像素是ROI内部的像素时,αext=255,并且当目标像素是ROI外部的像素时,αext被设置为显示在ROI外部透明度指定界面43上的ROI外部不透明度指定滑杆所指定的值。
类似地,通过对作为每一像素12位灰度分级(4096灰度分级)的第二光声图像的各个亮度值V2pix进行对比度调整而获得的值V2roi将如公式6中所示。
(1)当V2pix≥ROI2max时,V2roi=319
(2)当ROI2min<V2pix<R012max时,V2roi=319×(V2pix-ROI2min)/(4096×(ROI2max-ROI2min))
(3)当V2pix≤ROI2min时,V2roi=0
...公式6(0≤V2roi≤319)
当V2roi被颜色坐标取代时,实现以下。
(1)当0≤V2roi<128时,(R,G,B,α)=(0,0,V2roi+128,αext)
(2)当128≤V2roi≤319时,(R,G,B,α)=(V2roi-128,V2roi-128,255,αext)
...公式7
然而,当目标像素是ROI内部的像素时,αext=255,并且,当目标像素是ROI外部的像素时,αext被设置为显示在ROI外部透明度指定界面4c上的ROI外部不透明度指定滑杆所指定的值。
如上所述,通过使用公式4至7将颜色分配给两种类型的光声图像的ROI内部的各个像素,可以执行对比度调整和不透明度调整。指出,当然,分配颜色的方法可以不同于所示的颜色显示分配。
根据第二实施例的光声成像装置叠加已经经过如上所述的对比度调整和不透明度调整的第一光声图像和第二光声图像,并在图像显示区域41上显示叠加图像。
如以上所解释的,本发明不限于叠加光声图像和超声图像,还可以应用于叠加并显示不同的光声图像的情况。
通过第二实施例,可以在对多个光声图像单个地执行对比度调整和不透明度调整的情况下叠加并显示这些光声图像,从而改进ROI内部的信号的可见性并且使ROI外部的信号(噪声、伪像)变得不明显。
指出,尽管第二实施例例示说明了提供两个UI(每个用于执行对比度调整和不透明度调整)并对两个图像执行处理的情况,但是还可以对三个或更多个图像中的每个执行对比度调整和不透明度调整并且随后叠加这些图像。
指出,各个实施例的解释是用于解释本发明的举例说明,本发明可以通过适当地将实施例改变或组合到这样的改变或组合将不偏离本发明的主旨的程度来实现。
例如,尽管实施例解释了通过指定灰度级图像中的亮度值范围来执行对比度调整的情况,但是输入图像还可以不同于灰度级图像。在前述情况下,还可以基于像素值(即,各个颜色的亮度值)执行对比度调整。
本发明可以被实现为控制被检体信息获取装置的、包括前述处理中的至少一部分的方法。前述处理和手段可以通过自由组合实现,只要没有技术矛盾发生即可。
尽管已经参照示例性实施例描述了本发明,但是要理解本发明不限于所公开的示例性实施例。权利要求的范围要遵循最广泛的解释,以便包含所有这样的修改以及等同的结构和功能。
Claims (10)
1.一种被检体信息获取装置,包括:
光声图像获取单元,所述光声图像获取单元被配置为将测量光发射到被检体,接收在所述被检体中产生的光声波,并基于光声波产生使与所述被检体内的光学特性相关的信息可视化的第一图像;
超声图像获取单元,所述超声图像获取单元被配置为将超声波发送到所述被检体,接收在所述被检体中反射的超声回声,并基于超声回声产生使与所述被检体内的声学特性相关的信息可视化的第二图像;
感兴趣区域指定单元,所述感兴趣区域指定单元被配置为接收关于第一图像对感兴趣区域的指定;
图像处理单元,所述图像处理单元被配置为对第一图像在所指定的感兴趣区域内部和所指定的感兴趣区域外部分别使用不同的图像处理参数执行图像处理;以及
图像合成单元,所述图像合成单元被配置为叠加并合成已经经过了所述图像处理的第一图像和第二图像。
2.根据权利要求1所述的被检体信息获取装置,
其中,所述图像处理单元被配置为获取第一图像的感兴趣区域内部的像素值的频率分布,并基于所述频率分布对第一图像执行对比度调整。
3.根据权利要求2所述的被检体信息获取装置,
其中,所述图像处理单元被配置为使用第一图像的感兴趣区域内部的像素值的最大值和最小值对第一图像执行对比度调整。
4.根据权利要求1至3中的任何一个所述的被检体信息获取装置,还包括:
像素值范围指定单元,所述像素值范围指定单元被配置为接收在第一图像中将被强调的像素值的范围的指定,
其中,所述图像处理单元被配置为使用所指定的像素值范围对第一图像执行对比度调整。
5.根据权利要求1所述的被检体信息获取装置,还包括:
透明度指定单元,所述透明度指定单元被配置为接收第一图像的感兴趣区域外部的透明度的指定,
其中,所述图像处理单元被配置为对第一图像的感兴趣区域外部的像素设置所指定的透明度。
6.一种被检体信息获取装置,包括:
光声图像获取单元,所述光声图像获取单元被配置为将不同波长的测量光发射到被检体,对于每个波长,接收在所述被检体中产生的光声波,并且对于每个波长,基于光声波产生使与所述被检体内的光学特性相关的信息可视化的图像;
感兴趣区域指定单元,所述感兴趣区域指定单元被配置为接收对感兴趣区域的指定;
图像处理单元,所述图像处理单元被配置为对多个图像中的每个在感兴趣区域内部和外部分别使用不同的图像处理参数执行图像处理;以及
图像合成单元,所述图像合成单元被配置为叠加并合成已经经过了所述图像处理的所述多个图像。
7.根据权利要求6所述的被检体信息获取装置,
其中,所述图像处理单元被配置为对于所述多个图像中的每个,获取感兴趣区域内部的像素值的频率分布,并基于所述频率分布对各个图像执行对比度调整。
8.根据权利要求7所述的被检体信息获取装置,
其中,所述图像处理单元被配置为对所述各个图像使用对所述多个图像中的每个获取的感兴趣区域内部的像素值的最大值和最小值来执行对比度调整。
9.根据权利要求6至8中的任何一个所述的被检体信息获取装置,
像素值范围指定单元,所述像素值范围指定单元被配置为接收对于所述多个图像中的每个将被强调的像素值的范围的指定,
其中,所述图像处理单元被配置为分别使用所指定的像素值范围对各个图像执行对比度调整。
10.根据权利要求6所述的被检体信息获取装置,还包括:
透明度指定单元,所述透明度指定单元被配置为对所述多个图像中的每个,接收感兴趣区域外部的透明度的指定,
其中,所述图像处理单元被配置为对各个图像的感兴趣区域外部的像素设置所指定的透明度。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2012-286546 | 2012-12-28 | ||
JP2012286546A JP6103931B2 (ja) | 2012-12-28 | 2012-12-28 | 被検体情報取得装置、被検体情報取得方法 |
CN201310719112.8A CN103908292B (zh) | 2012-12-28 | 2013-12-24 | 被检体信息获取装置和被检体信息获取方法 |
Related Parent Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310719112.8A Division CN103908292B (zh) | 2012-12-28 | 2013-12-24 | 被检体信息获取装置和被检体信息获取方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN105832299A true CN105832299A (zh) | 2016-08-10 |
Family
ID=49726531
Family Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310719112.8A Active CN103908292B (zh) | 2012-12-28 | 2013-12-24 | 被检体信息获取装置和被检体信息获取方法 |
CN201610156516.4A Pending CN105832299A (zh) | 2012-12-28 | 2013-12-24 | 被检体信息获取装置和被检体信息获取方法 |
Family Applications Before (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310719112.8A Active CN103908292B (zh) | 2012-12-28 | 2013-12-24 | 被检体信息获取装置和被检体信息获取方法 |
Country Status (4)
Country | Link |
---|---|
US (2) | US20140187936A1 (zh) |
EP (1) | EP2749208B2 (zh) |
JP (1) | JP6103931B2 (zh) |
CN (2) | CN103908292B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108095689A (zh) * | 2016-11-25 | 2018-06-01 | 佳能株式会社 | 光声装置、信息处理方法及存储程序的非暂时性存储介质 |
WO2020082270A1 (zh) * | 2018-10-24 | 2020-04-30 | 中国医学科学院北京协和医院 | 一种成像方法以及成像系统 |
CN116942200A (zh) * | 2023-09-20 | 2023-10-27 | 杭州励影光电成像有限责任公司 | 一种非复用式超声多模态成像系统及方法 |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10143381B2 (en) | 2013-04-19 | 2018-12-04 | Canon Kabushiki Kaisha | Object information acquiring apparatus and control method therefor |
JP6215057B2 (ja) * | 2014-01-09 | 2017-10-18 | 富士通株式会社 | 可視化装置、可視化プログラムおよび可視化方法 |
EP3220831B1 (en) * | 2014-11-18 | 2019-06-26 | Koninklijke Philips N.V. | Visualization apparatus for property change of a tissue |
USD860253S1 (en) * | 2016-09-08 | 2019-09-17 | Canon Kabushiki Kaisha | Display screen with icon |
JP6708529B2 (ja) * | 2016-10-07 | 2020-06-10 | キヤノン株式会社 | 制御装置、制御方法、制御システム及びプログラム。 |
WO2018116891A1 (ja) * | 2016-12-19 | 2018-06-28 | オリンパス株式会社 | 画像処理装置、超音波診断システム、画像処理装置の作動方法、及び画像処理装置の作動プログラム |
WO2019138758A1 (ja) * | 2018-01-11 | 2019-07-18 | 古野電気株式会社 | レーダー装置 |
FR3081699B1 (fr) * | 2018-06-01 | 2020-10-02 | Cartimage Medical | Sonde photoacoustique endoscopique |
CN110832840A (zh) * | 2018-08-29 | 2020-02-21 | 深圳市大疆创新科技有限公司 | 一种图像的画质调节方法及系统、自主移动平台 |
CN111292676B (zh) * | 2018-11-20 | 2021-09-07 | 群创光电股份有限公司 | 电子装置 |
WO2021092250A1 (en) * | 2019-11-05 | 2021-05-14 | California Institute Of Technology | Spatiotemporal antialiasing in photoacoustic computed tomography |
US11830176B2 (en) * | 2021-09-12 | 2023-11-28 | Nanya Technology Corporation | Method of measuring a semiconductor device |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5042077A (en) * | 1987-10-02 | 1991-08-20 | General Electric Company | Method of highlighting subtle contrast in graphical images |
US20030068074A1 (en) * | 2001-10-05 | 2003-04-10 | Horst Hahn | Computer system and a method for segmentation of a digital image |
CN1575770A (zh) * | 2003-07-02 | 2005-02-09 | 株式会社东芝 | 显示患者有关信息之影像的形成方法和仪器 |
US7062714B1 (en) * | 2000-07-28 | 2006-06-13 | Ge Medical Systems Global Technology Company, Llc | Imaging system having preset processing parameters adapted to user preferences |
Family Cites Families (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0646263B1 (en) | 1993-04-20 | 2000-05-31 | General Electric Company | Computer graphic and live video system for enhancing visualisation of body structures during surgery |
US5797397A (en) * | 1996-11-25 | 1998-08-25 | Hewlett-Packard Company | Ultrasound imaging system and method using intensity highlighting to facilitate tissue differentiation |
DE10141186A1 (de) | 2001-08-22 | 2003-03-20 | Siemens Ag | Einrichtung zum Bearbeiten von Bildern, insbesondere medizinischen Bildern |
US6973219B2 (en) * | 2001-11-21 | 2005-12-06 | Ge Medical Systems Global Technology Company, Llc | Segmentation driven image noise reduction filter |
AU2003300052A1 (en) | 2002-12-31 | 2004-08-10 | John Herbert Cafarella | Multi-sensor breast tumor detection |
JP4426225B2 (ja) | 2003-07-02 | 2010-03-03 | Hoya株式会社 | 蛍光観察内視鏡システム及び蛍光観察内視鏡用光源装置 |
KR20090088909A (ko) | 2006-12-19 | 2009-08-20 | 코닌클리케 필립스 일렉트로닉스 엔.브이. | 결합된 광음향 및 초음파 이미징 시스템 |
GB0712432D0 (en) | 2007-06-26 | 2007-08-08 | Isis Innovation | Improvements in or relating to determination and display of material properties |
US8023760B1 (en) * | 2007-12-06 | 2011-09-20 | The United States Of America As Represented By The Secretary Of The Navy | System and method for enhancing low-visibility imagery |
WO2010009412A2 (en) * | 2008-07-18 | 2010-01-21 | University Of Rochester Medical Center | Low-cost device for c-scan photoacoustic imaging |
JP5257375B2 (ja) * | 2009-03-13 | 2013-08-07 | オムロン株式会社 | 画像処理装置および画像処理方法 |
JP5627328B2 (ja) | 2010-07-28 | 2014-11-19 | キヤノン株式会社 | 光音響診断装置 |
US20130338501A1 (en) * | 2012-06-13 | 2013-12-19 | Seno Medical Instruments, Inc. | System and method for storing data associated with the operation of a dual modality optoacoustic/ultrasound system |
JP5502686B2 (ja) | 2010-09-30 | 2014-05-28 | 富士フイルム株式会社 | 光音響画像診断装置、画像生成方法、及びプログラム |
JP5777394B2 (ja) * | 2011-04-28 | 2015-09-09 | 富士フイルム株式会社 | 光音響画像化方法および装置 |
JP5683383B2 (ja) * | 2011-05-24 | 2015-03-11 | 富士フイルム株式会社 | 光音響撮像装置およびその作動方法 |
KR102090840B1 (ko) | 2011-10-12 | 2020-03-18 | 세노 메디컬 인스투르먼츠 인코포레이티드 | 광음향 데이터를 획득하며 그것의 파라메트릭 맵들을 생성하기 위한 시스템 및 방법 |
-
2012
- 2012-12-28 JP JP2012286546A patent/JP6103931B2/ja active Active
-
2013
- 2013-12-03 EP EP13195412.5A patent/EP2749208B2/en active Active
- 2013-12-20 US US14/135,705 patent/US20140187936A1/en not_active Abandoned
- 2013-12-24 CN CN201310719112.8A patent/CN103908292B/zh active Active
- 2013-12-24 CN CN201610156516.4A patent/CN105832299A/zh active Pending
-
2019
- 2019-03-15 US US16/355,409 patent/US20190209126A1/en not_active Abandoned
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5042077A (en) * | 1987-10-02 | 1991-08-20 | General Electric Company | Method of highlighting subtle contrast in graphical images |
US7062714B1 (en) * | 2000-07-28 | 2006-06-13 | Ge Medical Systems Global Technology Company, Llc | Imaging system having preset processing parameters adapted to user preferences |
US20030068074A1 (en) * | 2001-10-05 | 2003-04-10 | Horst Hahn | Computer system and a method for segmentation of a digital image |
CN1575770A (zh) * | 2003-07-02 | 2005-02-09 | 株式会社东芝 | 显示患者有关信息之影像的形成方法和仪器 |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108095689A (zh) * | 2016-11-25 | 2018-06-01 | 佳能株式会社 | 光声装置、信息处理方法及存储程序的非暂时性存储介质 |
WO2020082270A1 (zh) * | 2018-10-24 | 2020-04-30 | 中国医学科学院北京协和医院 | 一种成像方法以及成像系统 |
CN116942200A (zh) * | 2023-09-20 | 2023-10-27 | 杭州励影光电成像有限责任公司 | 一种非复用式超声多模态成像系统及方法 |
CN116942200B (zh) * | 2023-09-20 | 2024-02-06 | 杭州励影光电成像有限责任公司 | 一种非复用式超声多模态成像系统及方法 |
Also Published As
Publication number | Publication date |
---|---|
EP2749208B1 (en) | 2019-04-24 |
CN103908292B (zh) | 2016-04-06 |
EP2749208B2 (en) | 2022-08-03 |
JP6103931B2 (ja) | 2017-03-29 |
EP2749208A1 (en) | 2014-07-02 |
CN103908292A (zh) | 2014-07-09 |
JP2014128318A (ja) | 2014-07-10 |
US20140187936A1 (en) | 2014-07-03 |
US20190209126A1 (en) | 2019-07-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105832299A (zh) | 被检体信息获取装置和被检体信息获取方法 | |
US11715202B2 (en) | Analyzing apparatus and analyzing method | |
US10772606B2 (en) | Method and apparatus for displaying ultrasound images | |
US6641538B2 (en) | Ultrasonic diagnostic apparatus and method of controlling a ultrasonic diagnostic apparatus | |
JP4657106B2 (ja) | 超音波診断装置 | |
US20090015587A1 (en) | Ultrasonic imaging apparatus | |
US10799215B2 (en) | Ultrasound systems, methods and apparatus for associating detection information of the same | |
JP2018529453A (ja) | 心拍出量を測定するための方法およびシステム | |
JP6608126B2 (ja) | 超音波診断装置 | |
WO2007043310A1 (ja) | 画像表示方法及び医用画像診断システム | |
CN108283505A (zh) | 超声波诊断装置、图像处理装置及图像处理方法 | |
JP7490839B2 (ja) | 解析装置及び解析プログラム | |
WO2021042298A1 (zh) | 一种vti测量装置及方法 | |
JP5113548B2 (ja) | 超音波画像処理装置 | |
WO2012039192A1 (ja) | 超音波診断装置および超音波画像の表示方法 | |
US20180146954A1 (en) | Method of ultrasound apparatus parameters configuration and an ultrasound apparatus of using the same | |
JP5865251B2 (ja) | 超音波診断装置および画像処理装置 | |
KR101014559B1 (ko) | 3차원 초음파 영상 제공 초음파 시스템 및 방법 | |
CN106999154B (zh) | 医学介入设备与组织的基于单模态的视觉区分 | |
JP7286025B2 (ja) | 胎盤を評価するためのシステム及び方法 | |
JP6786260B2 (ja) | 超音波診断装置及び画像生成方法 | |
KR102578754B1 (ko) | 초음파 이미지 디스플레이 방법 및 이를 위한 장치 | |
JP7304947B2 (ja) | 超音波制御ユニット | |
CN116407149A (zh) | 超声成像方法、超声成像系统及非暂态计算机可读介质 | |
JP2019155004A (ja) | 光音響装置および被検体情報取得方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20160810 |
|
RJ01 | Rejection of invention patent application after publication |