JP2018161467A - 画像処理装置および画像処理方法 - Google Patents
画像処理装置および画像処理方法 Download PDFInfo
- Publication number
- JP2018161467A JP2018161467A JP2018045346A JP2018045346A JP2018161467A JP 2018161467 A JP2018161467 A JP 2018161467A JP 2018045346 A JP2018045346 A JP 2018045346A JP 2018045346 A JP2018045346 A JP 2018045346A JP 2018161467 A JP2018161467 A JP 2018161467A
- Authority
- JP
- Japan
- Prior art keywords
- subject
- dimensional data
- image
- image processing
- voxels
- 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
- 238000012545 processing Methods 0.000 title claims abstract description 74
- 238000003672 processing method Methods 0.000 title claims description 18
- 238000004364 calculation method Methods 0.000 claims abstract description 7
- 238000000034 method Methods 0.000 claims description 81
- 230000008569 process Effects 0.000 claims description 42
- 238000001914 filtration Methods 0.000 claims description 39
- 238000003860 storage Methods 0.000 claims description 11
- 230000007423 decrease Effects 0.000 claims description 4
- 238000001514 detection method Methods 0.000 abstract description 4
- 238000009826 distribution Methods 0.000 description 35
- 239000000523 sample Substances 0.000 description 17
- 239000006096 absorbing agent Substances 0.000 description 16
- 230000003287 optical effect Effects 0.000 description 15
- 238000010521 absorption reaction Methods 0.000 description 14
- 210000004204 blood vessel Anatomy 0.000 description 12
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 11
- 239000001301 oxygen Substances 0.000 description 11
- 229910052760 oxygen Inorganic materials 0.000 description 11
- 210000000481 breast Anatomy 0.000 description 8
- 238000009877 rendering Methods 0.000 description 8
- 102000001554 Hemoglobins Human genes 0.000 description 7
- 108010054147 Hemoglobins Proteins 0.000 description 7
- 210000004369 blood Anatomy 0.000 description 6
- 239000008280 blood Substances 0.000 description 6
- 230000006870 function Effects 0.000 description 6
- 239000000463 material Substances 0.000 description 6
- 238000005259 measurement Methods 0.000 description 6
- 238000006243 chemical reaction Methods 0.000 description 5
- 210000003743 erythrocyte Anatomy 0.000 description 5
- 239000000126 substance Substances 0.000 description 5
- XUMBMVFBXHLACL-UHFFFAOYSA-N Melanin Chemical compound O=C1C(=O)C(C2=CNC3=C(C(C(=O)C4=C32)=O)C)=C2C4=CNC2=C1C XUMBMVFBXHLACL-UHFFFAOYSA-N 0.000 description 4
- 230000031700 light absorption Effects 0.000 description 4
- 238000010895 photoacoustic effect Methods 0.000 description 4
- 210000003462 vein Anatomy 0.000 description 4
- 210000001367 artery Anatomy 0.000 description 3
- 201000011510 cancer Diseases 0.000 description 3
- 238000003384 imaging method Methods 0.000 description 3
- 238000003325 tomography Methods 0.000 description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 3
- WQZGKKKJIJFFOK-GASJEMHNSA-N Glucose Natural products OC[C@H]1OC(O)[C@H](O)[C@@H](O)[C@@H]1O WQZGKKKJIJFFOK-GASJEMHNSA-N 0.000 description 2
- 241001465754 Metazoa Species 0.000 description 2
- 206010028980 Neoplasm Diseases 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000003745 diagnosis Methods 0.000 description 2
- 238000009792 diffusion process Methods 0.000 description 2
- 239000008103 glucose Substances 0.000 description 2
- 230000010365 information processing Effects 0.000 description 2
- 230000001678 irradiating effect Effects 0.000 description 2
- HFGPZNIAWCZYJU-UHFFFAOYSA-N lead zirconate titanate Chemical compound [O-2].[O-2].[O-2].[O-2].[O-2].[Ti+4].[Zr+4].[Pb+2] HFGPZNIAWCZYJU-UHFFFAOYSA-N 0.000 description 2
- 229910052451 lead zirconate titanate Inorganic materials 0.000 description 2
- 239000013307 optical fiber Substances 0.000 description 2
- 230000001902 propagating effect Effects 0.000 description 2
- PXFBZOLANLWPMH-UHFFFAOYSA-N 16-Epiaffinine Natural products C1C(C2=CC=CC=C2N2)=C2C(=O)CC2C(=CC)CN(C)C1C2CO PXFBZOLANLWPMH-UHFFFAOYSA-N 0.000 description 1
- 102000008186 Collagen Human genes 0.000 description 1
- 108010035532 Collagen Proteins 0.000 description 1
- 238000005481 NMR spectroscopy Methods 0.000 description 1
- 239000002033 PVDF binder Substances 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 229920001436 collagen Polymers 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000007789 gas Substances 0.000 description 1
- 238000001727 in vivo Methods 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 229920000642 polymer Polymers 0.000 description 1
- 238000003825 pressing Methods 0.000 description 1
- 102000004169 proteins and genes Human genes 0.000 description 1
- 108090000623 proteins and genes Proteins 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 208000019553 vascular disease Diseases 0.000 description 1
- 238000012800 visualization Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/008—Specific post-processing after tomographic reconstruction, e.g. voxelisation, metal artifact correction
-
- 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
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T15/00—3D [Three Dimensional] image rendering
- G06T15/08—Volume rendering
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/70—Denoising; Smoothing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/73—Deblurring; Sharpening
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B2576/00—Medical imaging apparatus involving image processing or analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10084—Hybrid tomography; Concurrent acquisition with multiple different tomographic modalities
-
- 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/20004—Adaptive image processing
- G06T2207/20012—Locally adaptive
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20024—Filtering details
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30024—Cell structures in vitro; Tissue sections in vitro
-
- 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/30088—Skin; Dermal
-
- 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
-
- 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
-
- 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
- G16H30/00—ICT specially adapted for the handling or processing of medical images
- G16H30/40—ICT specially adapted for the handling or processing of medical images for processing medical images, e.g. editing
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Medical Informatics (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Acoustics & Sound (AREA)
- Heart & Thoracic Surgery (AREA)
- Pathology (AREA)
- Biophysics (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Biomedical Technology (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Quality & Reliability (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Computer Graphics (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
【課題】三次元データを二次元画像に変換する画像処理装置において、奥行き方向の情報を区別しやすい画像を提供する。【解決手段】被検体の特性情報を複数のボクセルによって表した三次元データを取得する情報取得手段と、前記被検体の表面形状に関する情報を取得する形状情報取得手段と、前記表面形状に関する情報に基づいて、前記ボクセルのそれぞれについて、前記被検体の表面から当該ボクセルに対応する前記被検体内の位置までの距離を算出する距離算出手段と、前記ボクセルのそれぞれに対して、前記算出した距離に応じたぼかし処理を含むフィルタ処理を行うフィルタ手段と、前記フィルタ処理後の前記三次元データに基づいて二次元画像を生成する画像生成手段と、を有する。【選択図】図１
Description
本発明は、被検体を撮像して得られた画像を処理する装置に関する。
近年、医療分野において、被検体内の構造情報や、生理的情報、すなわち機能情報をイメージングするための研究が進められている。このような技術の一つとして、近年、光音響トモグラフィ（ＰＡＴ：PhotoAcoustic Tomography）が提案されている。
レーザ光などの光を被検体である生体に照射すると、光が被検体内の生体組織で吸収される際に音響波（典型的には超音波）が発生する。この現象を光音響効果と呼び、光音響効果により発生した音響波を光音響波と呼ぶ。被検体を構成する組織は、光エネルギーの吸収率がそれぞれ異なるため、発生する光音響波の音圧も異なったものとなる。ＰＡＴでは、発生した光音響波を探触子で受信し、受信信号を数学的に解析することにより、被検体内の特性情報を取得することができる。
レーザ光などの光を被検体である生体に照射すると、光が被検体内の生体組織で吸収される際に音響波（典型的には超音波）が発生する。この現象を光音響効果と呼び、光音響効果により発生した音響波を光音響波と呼ぶ。被検体を構成する組織は、光エネルギーの吸収率がそれぞれ異なるため、発生する光音響波の音圧も異なったものとなる。ＰＡＴでは、発生した光音響波を探触子で受信し、受信信号を数学的に解析することにより、被検体内の特性情報を取得することができる。
特性情報とは、初期音圧分布、光吸収エネルギー密度分布、光吸収係数分布などの光学特性値の分布である。また、これらの情報を複数波長の光を用いて取得することにより、被検体内の特定物質の濃度（血液中のヘモグロビン濃度や、血液の酸素飽和度など）を定量的に計測することができる。
光音響トモグラフィでは、被検体情報を三次元情報として取得することができる。
一般的に、画像の表示は二次元平面上にて行うため、三次元情報を二次元平面に投影する必要がある。被検体を測定して得られた三次元データを二次元画像に変換する方法として、奥行き方向において最大値を持つボクセルを投影する最大値投影法、視点から見て最も手前に位置する信号のみを可視化するサーフェイスレンダリング法などが知られている。
一般的に、画像の表示は二次元平面上にて行うため、三次元情報を二次元平面に投影する必要がある。被検体を測定して得られた三次元データを二次元画像に変換する方法として、奥行き方向において最大値を持つボクセルを投影する最大値投影法、視点から見て最も手前に位置する信号のみを可視化するサーフェイスレンダリング法などが知られている。
特許文献１には、光音響トモグラフィにより取得された三次元画像データを、最大値投影法にて二次元平面に投影する方法が記載されている。
ところが、特許文献１に記載の方法では、三次元データを二次元画像に変換する場合、奥行き方向の情報を区別しにくいという課題がある。
本発明はこのような従来技術の課題に鑑みてなされたものであり、三次元データを二次元画像に変換する画像処理装置において、奥行き方向の情報を区別しやすい画像を提供することを目的とする。
本発明に係る画像処理装置は、
被検体の特性情報を複数のボクセルによって表した三次元データを取得する情報取得手段と、前記被検体の表面形状に関する情報を取得する形状情報取得手段と、前記表面形状に関する情報に基づいて、前記ボクセルのそれぞれについて、前記被検体の表面から当該
ボクセルに対応する前記被検体内の位置までの距離を算出する距離算出手段と、前記ボクセルのそれぞれに対して、前記算出した距離に応じたぼかし処理を含むフィルタ処理を行うフィルタ手段と、前記フィルタ処理後の前記三次元データに基づいて二次元画像を生成する画像生成手段と、を有することを特徴とする。
被検体の特性情報を複数のボクセルによって表した三次元データを取得する情報取得手段と、前記被検体の表面形状に関する情報を取得する形状情報取得手段と、前記表面形状に関する情報に基づいて、前記ボクセルのそれぞれについて、前記被検体の表面から当該
ボクセルに対応する前記被検体内の位置までの距離を算出する距離算出手段と、前記ボクセルのそれぞれに対して、前記算出した距離に応じたぼかし処理を含むフィルタ処理を行うフィルタ手段と、前記フィルタ処理後の前記三次元データに基づいて二次元画像を生成する画像生成手段と、を有することを特徴とする。
また、本発明に係る画像処理方法は、
被検体の特性情報を複数のボクセルによって表した三次元データを取得する情報取得ステップと、前記被検体の表面形状に関する情報を取得する形状情報取得ステップと、前記表面形状に関する情報に基づいて、前記ボクセルのそれぞれについて、前記被検体の表面から当該ボクセルに対応する前記被検体内の位置までの距離を算出する距離算出ステップと、前記ボクセルのそれぞれに対して、前記算出した距離に応じたぼかし処理を含むフィルタ処理を行うフィルタステップと、前記フィルタ処理後の前記三次元データに基づいて二次元画像を生成する画像生成ステップと、を含むことを特徴とする。
被検体の特性情報を複数のボクセルによって表した三次元データを取得する情報取得ステップと、前記被検体の表面形状に関する情報を取得する形状情報取得ステップと、前記表面形状に関する情報に基づいて、前記ボクセルのそれぞれについて、前記被検体の表面から当該ボクセルに対応する前記被検体内の位置までの距離を算出する距離算出ステップと、前記ボクセルのそれぞれに対して、前記算出した距離に応じたぼかし処理を含むフィルタ処理を行うフィルタステップと、前記フィルタ処理後の前記三次元データに基づいて二次元画像を生成する画像生成ステップと、を含むことを特徴とする。
また、本発明に係る非一時的記憶媒体は、
被検体の特性情報を複数のボクセルによって表した三次元データを取得する情報取得ステップと、前記被検体の表面形状に関する情報を取得する形状情報取得ステップと、前記表面形状に関する情報に基づいて、前記ボクセルのそれぞれについて、前記被検体の表面から当該ボクセルに対応する前記被検体内の位置までの距離を算出する距離算出ステップと、前記ボクセルのそれぞれに対して、前記算出した距離に応じたぼかし処理を含むフィルタ処理を行うフィルタステップと、前記フィルタ処理後の前記三次元データに基づいて二次元画像を生成する画像生成ステップと、を含む画像処理方法をコンピュータに実行させるためのプログラムが記録されたことを特徴とする。
被検体の特性情報を複数のボクセルによって表した三次元データを取得する情報取得ステップと、前記被検体の表面形状に関する情報を取得する形状情報取得ステップと、前記表面形状に関する情報に基づいて、前記ボクセルのそれぞれについて、前記被検体の表面から当該ボクセルに対応する前記被検体内の位置までの距離を算出する距離算出ステップと、前記ボクセルのそれぞれに対して、前記算出した距離に応じたぼかし処理を含むフィルタ処理を行うフィルタステップと、前記フィルタ処理後の前記三次元データに基づいて二次元画像を生成する画像生成ステップと、を含む画像処理方法をコンピュータに実行させるためのプログラムが記録されたことを特徴とする。
本発明によれば、三次元データを二次元画像に変換する画像処理装置において、奥行き方向の情報を区別しやすい画像を提供することができる。
以下に図面を参照しつつ、本発明の好適な実施の形態について説明する。ただし、以下に記載されている構成部品の寸法、材質、形状およびそれらの相対配置などは、発明が適用される装置の構成や各種条件により適宜変更されるべきものである。よって、この発明の範囲を以下の記載に限定する趣旨のものではない。
本発明は、被検体から伝播する音響波に基づいて生成した被検体内部の特性情報を、二次元画像に変換する技術に関する。よって本発明は、画像処理装置または画像処理方法として捉えることができる。また、当該画像処理装置を含む被検体情報取得装置、画像処理方法を含む被検体情報取得方法として捉えることもできる。
本発明はまた、これらの方法をＣＰＵやメモリ等のハードウェア資源を備える情報処理装置に実行させるプログラムや、そのプログラムを格納した記憶媒体や、情報処理装置として捉えることもできる。
本発明はまた、これらの方法をＣＰＵやメモリ等のハードウェア資源を備える情報処理装置に実行させるプログラムや、そのプログラムを格納した記憶媒体や、情報処理装置として捉えることもできる。
本発明の画像処理装置は、被検体に光（電磁波）を照射することにより被検体内で発生した音響波を受信して、被検体の特性情報を画像データとして取得する光音響効果を利用した装置に適用することができる。この場合、特性情報とは、光音響波を受信することにより得られる受信信号を用いて生成される、被検体内の複数位置のそれぞれに対応する特性値の情報である。
光音響測定により取得される特性情報は、光エネルギーの吸収率を反映した値である。例えば、光照射によって生じた音響波の発生源、被検体内の初期音圧、あるいは初期音圧から導かれる光エネルギー吸収密度や吸収係数、組織を構成する物質の濃度を含む。また、物質濃度として酸化ヘモグロビン濃度と還元ヘモグロビン濃度を求めることにより、酸素飽和度分布を算出できる。また、グルコース濃度、コラーゲン濃度、メラニン濃度、脂肪や水の体積分率なども求められる。
被検体内の各位置の特性情報に基づいて、二次元または三次元の特性情報分布が得られる。分布データは画像データとして生成され得る。特性情報は、数値データとしてではなく、被検体内の各位置の分布情報として求めてもよい。すなわち、初期音圧分布、エネルギー吸収密度分布、吸収係数分布や酸素飽和度分布などの分布情報である。
本明細書における音響波とは、典型的には超音波であり、音波、音響波と呼ばれる弾性波を含む。探触子等により音響波から変換された電気信号を音響信号とも呼ぶ。ただし、本明細書における超音波または音響波という記載には、それらの弾性波の波長を限定する意図はない。光音響効果により発生した音響波は、光音響波または光超音波と呼ばれる。光音響波に由来する電気信号を光音響信号とも呼ぶ。
（第一の実施形態）
<システム構成>
図１は、第一の実施形態に係る被検体情報取得装置の構成を示す機能ブロック図である。
第一の実施形態に係る被検体情報取得装置は光音響装置であり、光源１０１、光学系１０２、音響整合材１０５、音響波探触子１０６、信号処理部１０７、データ処理部１０８、表示部１０９を有している。なお、符号１０３は被検体であり、符号１０４は光吸収体である。
<システム構成>
図１は、第一の実施形態に係る被検体情報取得装置の構成を示す機能ブロック図である。
第一の実施形態に係る被検体情報取得装置は光音響装置であり、光源１０１、光学系１０２、音響整合材１０５、音響波探触子１０６、信号処理部１０７、データ処理部１０８、表示部１０９を有している。なお、符号１０３は被検体であり、符号１０４は光吸収体である。
<<被検体１０３>>
被検体１０３は、本実施形態に係る被検体情報取得装置の一部を構成するものではないが、以下に説明する。本実施形態に係る被検体情報取得装置は、人や動物の悪性腫瘍や血管疾患、血糖値などの診断や化学治療の経過観察を行う装置である。よって、被検体としては、生体、具体的には人体や動物の乳房や指、手足などが想定される。
被検体の内部には、生体であれば水や脂肪、タンパク質、酸化ヘモグロビン、還元ヘモグロビンなど、光吸収係数の大きい光吸収体が存在するため、光照射に起因して光音響波が発生する。被検体としてファントムを用いる場合、光学特性を模擬した物質を光吸収体として内部に封入することで、光音響波の発生および測定が可能となる。
被検体１０３は、本実施形態に係る被検体情報取得装置の一部を構成するものではないが、以下に説明する。本実施形態に係る被検体情報取得装置は、人や動物の悪性腫瘍や血管疾患、血糖値などの診断や化学治療の経過観察を行う装置である。よって、被検体としては、生体、具体的には人体や動物の乳房や指、手足などが想定される。
被検体の内部には、生体であれば水や脂肪、タンパク質、酸化ヘモグロビン、還元ヘモグロビンなど、光吸収係数の大きい光吸収体が存在するため、光照射に起因して光音響波が発生する。被検体としてファントムを用いる場合、光学特性を模擬した物質を光吸収体として内部に封入することで、光音響波の発生および測定が可能となる。
<<光源１０１>>
光源１０１は、被検体に照射するパルス光を発生させる装置である。光源は、大出力を得るためレーザ光源であることが望ましいが、レーザの代わりに発光ダイオードやフラッシュランプ等を用いることもできる。光源としてレーザを用いる場合、固体レーザ、ガスレーザ、色素レーザ、半導体レーザなど様々なものが使用できる。
また、パルス光の波長は、被検体を構成する成分のうち特定の成分に吸収される特定の
波長であって、被検体内部まで光が伝搬する波長であることが望ましい。具体的には、被検体が生体である場合、７００ｎｍ以上１２００ｎｍ以下であることが望ましい。この領域の光は比較的生体深部まで到達することができ、深部の情報を得ることができる。
また、光音響波を効果的に発生させるためには、被検体の熱特性に応じて十分短い時間に光を照射させなければならない。被検体が生体である場合、光源から発生するパルス光のパルス幅は数１０ナノ秒以下が好適である。
なお、光照射のタイミング、波形、強度等は、不図示の制御手段によって制御される。
光源１０１は、被検体に照射するパルス光を発生させる装置である。光源は、大出力を得るためレーザ光源であることが望ましいが、レーザの代わりに発光ダイオードやフラッシュランプ等を用いることもできる。光源としてレーザを用いる場合、固体レーザ、ガスレーザ、色素レーザ、半導体レーザなど様々なものが使用できる。
また、パルス光の波長は、被検体を構成する成分のうち特定の成分に吸収される特定の
波長であって、被検体内部まで光が伝搬する波長であることが望ましい。具体的には、被検体が生体である場合、７００ｎｍ以上１２００ｎｍ以下であることが望ましい。この領域の光は比較的生体深部まで到達することができ、深部の情報を得ることができる。
また、光音響波を効果的に発生させるためには、被検体の熱特性に応じて十分短い時間に光を照射させなければならない。被検体が生体である場合、光源から発生するパルス光のパルス幅は数１０ナノ秒以下が好適である。
なお、光照射のタイミング、波形、強度等は、不図示の制御手段によって制御される。
<<光学系１０２>>
光学系１０２は、光源１０１から発せられたパルス光を伝送する部材である。光源から出射された光は、レンズやミラーなどの光学部品により、所定の光分布形状に加工されながら被検体に導かれ、照射される。なお、光ファイバなどの光導波路などを用いて光を伝搬させることも可能である。
光学系１０２は、例えば、レンズ、ミラー、プリズム、光ファイバ、拡散板、シャッターなどの光学機器などを含んでいてもよい。このような光学部品は、光源から発せられた光が被検体に所望の形状で照射されれば、どのようなものを用いてもかまわない。なお、光はレンズで集光させるより、ある程度の面積に広げる方が、生体への安全性ならびに診断領域を広げられるという観点で好ましい。
光学系１０２は、光源１０１から発せられたパルス光を伝送する部材である。光源から出射された光は、レンズやミラーなどの光学部品により、所定の光分布形状に加工されながら被検体に導かれ、照射される。なお、光ファイバなどの光導波路などを用いて光を伝搬させることも可能である。
光学系１０２は、例えば、レンズ、ミラー、プリズム、光ファイバ、拡散板、シャッターなどの光学機器などを含んでいてもよい。このような光学部品は、光源から発せられた光が被検体に所望の形状で照射されれば、どのようなものを用いてもかまわない。なお、光はレンズで集光させるより、ある程度の面積に広げる方が、生体への安全性ならびに診断領域を広げられるという観点で好ましい。
<<音響波探触子１０６>>
音響波探触子１０６は、被検体１０３の内部から到来する音響波を受信して、電気信号に変換する手段である。音響波検出素子は、探触子、音響波探触子、音響波検出器、音響波受信器、トランスデューサとも呼ばれる。
生体から発生する音響波は、１００ＫＨｚから１００ＭＨｚの超音波であるため、音響波検出素子には、上記の周波数帯を受信できる素子を用いる。具体的には、圧電現象を用いたトランスデューサ、光の共振を用いたトランスデューサ、容量の変化を用いたトランスデューサなどを用いることができる。
音響波探触子１０６は、被検体１０３の内部から到来する音響波を受信して、電気信号に変換する手段である。音響波検出素子は、探触子、音響波探触子、音響波検出器、音響波受信器、トランスデューサとも呼ばれる。
生体から発生する音響波は、１００ＫＨｚから１００ＭＨｚの超音波であるため、音響波検出素子には、上記の周波数帯を受信できる素子を用いる。具体的には、圧電現象を用いたトランスデューサ、光の共振を用いたトランスデューサ、容量の変化を用いたトランスデューサなどを用いることができる。
また、音響素子は、感度が高く、周波数帯域が広いものを用いることが望ましい。具体的にはＰＺＴ（チタン酸ジルコン酸鉛）などを用いた圧電素子、ＰＶＤＦ（ポリフッ化ビニリデン）などの高分子圧電膜材料、ＣＭＵＴ（容量性マイクロマシン超音波トランスデューサ）、ファブリペロー干渉計を用いたものなどが挙げられる。ただし、ここに挙げたものだけに限定されず、探触子としての機能を満たすものであれば、どのようなものであってもよい。
音響波探触子１０６と被検体１０３との間には、音響インピーダンスを整合させるための部材である音響整合材１０５が配置される。音響整合材１０５には、ジェルや水、オイルなどを利用することができる。
信号処理部１０７は、取得した電気信号を増幅してデジタル信号に変換する手段である。
信号処理部１０７は、受信信号を増幅する増幅器、アナログの受信信号をデジタル変換するＡ／Ｄ変換器、受信信号を記憶するＦＩＦＯ等のメモリと、ＦＰＧＡチップ等の演算回路を用いて構成されてもよい。また、信号処理部１０７は、複数のプロセッサや演算回路から構成されていてもよい。
信号処理部１０７は、受信信号を増幅する増幅器、アナログの受信信号をデジタル変換するＡ／Ｄ変換器、受信信号を記憶するＦＩＦＯ等のメモリと、ＦＰＧＡチップ等の演算回路を用いて構成されてもよい。また、信号処理部１０７は、複数のプロセッサや演算回路から構成されていてもよい。
データ処理部１０８は、デジタル変換された信号（以下、光音響信号）に基づいて、被検体の内部の光吸収係数や酸素飽和度等といった被検体情報を取得する手段である。具体的には、収集された電気信号から三次元の被検体内の初期音圧分布を生成する。初期音圧分布の生成には、例えば、ユニバーサルバックプロジェクション（Ｕｎｉｖｅｒｓａｌ
Ｂａｃｋ−ｐｒｏｊｅｃｔｉｏｎ、以下ＵＢＰ）アルゴリズムやディレイアンドサム（Ｄｅｌａｙ ａｎｄ Ｓｕｍ）アルゴリズムを用いることができる。
Ｂａｃｋ−ｐｒｏｊｅｃｔｉｏｎ、以下ＵＢＰ）アルゴリズムやディレイアンドサム（Ｄｅｌａｙ ａｎｄ Ｓｕｍ）アルゴリズムを用いることができる。
また、データ処理部１０８は、被検体に照射される光量に関する情報に基づいて、被検体内の三次元の光強度分布を生成する。三次元の光強度分布は、二次元の光強度分布に関する情報から光拡散方程式を解くことで取得できる。また、光音響信号から生成された被検体内の初期音圧分布と、三次元の光強度分布とを用いて、被検体内の吸収係数分布を得ることができる。また、複数の波長における吸収係数分布を演算することで、被検体内の酸素飽和度分布を得ることができる。
また、データ処理部１０８は、後述する方法で、三次元の吸収係数分布を二次元画像に変換する。
データ処理部１０８が、本発明における形状情報取得手段、距離算出手段、画像生成手段、表示制御手段である。
また、データ処理部１０８は、後述する方法で、三次元の吸収係数分布を二次元画像に変換する。
データ処理部１０８が、本発明における形状情報取得手段、距離算出手段、画像生成手段、表示制御手段である。
データ処理部１０８は、ＣＰＵとＲＡＭ、および不揮発メモリ、制御ポートを有するコンピュータで構成してもよい。不揮発メモリに格納されたプログラムがＣＰＵで実行されることにより、各モジュールの制御が行われる。データ処理部１０８は、汎用コンピュータや専用に設計されたワークステーションであってもよい。
また、データ処理部１０８の演算機能を担うユニットは、ＣＰＵやＧＰＵ（Ｇｒａｐｈｉｃｓ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）等のプロセッサ、ＦＰＧＡ（Ｆｉｅｌｄ Ｐｒｏｇｒａｍｍａｂｌｅ Ｇａｔｅ Ａｒｒａｙ）チップ等の演算回路で構成されていてもよい。これらのユニットは、単一のプロセッサや演算回路から構成されるだけでなく、複数のプロセッサや演算回路から構成されていてもよい。
また、データ処理部１０８の記憶機能を担うユニットは、ＲＯＭ（Ｒｅａｄ ｏｎｌｙ
ｍｅｍｏｒｙ）、磁気ディスクやフラッシュメモリなどの非一時記憶媒体や、ＲＡＭ（Ｒａｎｄｏｍ Ａｃｃｅｓｓ Ｍｅｍｏｒｙ）などの揮発性の媒体であってもよい。なお、プログラムが格納される記憶媒体は、非一時記憶媒体である。なお、これらのユニットは、１つの記憶媒体から構成されるだけでなく、複数の記憶媒体から構成されていてもよい。データ処理部１０８の制御機能を担うユニットは、ＣＰＵなどの演算素子で構成される。
また、データ処理部１０８の演算機能を担うユニットは、ＣＰＵやＧＰＵ（Ｇｒａｐｈｉｃｓ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）等のプロセッサ、ＦＰＧＡ（Ｆｉｅｌｄ Ｐｒｏｇｒａｍｍａｂｌｅ Ｇａｔｅ Ａｒｒａｙ）チップ等の演算回路で構成されていてもよい。これらのユニットは、単一のプロセッサや演算回路から構成されるだけでなく、複数のプロセッサや演算回路から構成されていてもよい。
また、データ処理部１０８の記憶機能を担うユニットは、ＲＯＭ（Ｒｅａｄ ｏｎｌｙ
ｍｅｍｏｒｙ）、磁気ディスクやフラッシュメモリなどの非一時記憶媒体や、ＲＡＭ（Ｒａｎｄｏｍ Ａｃｃｅｓｓ Ｍｅｍｏｒｙ）などの揮発性の媒体であってもよい。なお、プログラムが格納される記憶媒体は、非一時記憶媒体である。なお、これらのユニットは、１つの記憶媒体から構成されるだけでなく、複数の記憶媒体から構成されていてもよい。データ処理部１０８の制御機能を担うユニットは、ＣＰＵなどの演算素子で構成される。
表示部１０９は、データ処理部で取得した情報及びその加工情報を表示する手段であり、典型的にはディスプレイ装置である。表示部１０９は、複数の表示部を備え、並列表示が可能な構成であってもよい。
<被検体の測定方法>
次に、本実施形態に係る被検体情報取得装置によって、被検体である生体を測定する方法について説明する。
まず、光源１０１から発せられたパルス光を被検体１０３に照射する。被検体の内部を伝搬した光のエネルギーの一部が血液などの光吸収体に吸収されると、熱膨張により当該光吸収体から音響波が発生する。生体内にがんが存在する場合は、がんの新生血管において他の正常部の血液と同様に光が特異的に吸収され、音響波が発生する。生体内で発生した光音響波は、音響波探触子１０６で受信される。
音響波探触子１０６が受信した信号は、信号処理部１０７で変換されたのち、データ処理部１０８で解析される。解析結果は、生体内の特性情報（例えば、初期音圧分布や吸収係数分布）を表す画像データに変換され、表示部１０９を介して出力される。本実施形態では、データ処理部１０８が吸収係数分布を表す三次元データをまず生成し、当該三次元データを二次元画像に変換して表示部１０９に出力する。
次に、本実施形態に係る被検体情報取得装置によって、被検体である生体を測定する方法について説明する。
まず、光源１０１から発せられたパルス光を被検体１０３に照射する。被検体の内部を伝搬した光のエネルギーの一部が血液などの光吸収体に吸収されると、熱膨張により当該光吸収体から音響波が発生する。生体内にがんが存在する場合は、がんの新生血管において他の正常部の血液と同様に光が特異的に吸収され、音響波が発生する。生体内で発生した光音響波は、音響波探触子１０６で受信される。
音響波探触子１０６が受信した信号は、信号処理部１０７で変換されたのち、データ処理部１０８で解析される。解析結果は、生体内の特性情報（例えば、初期音圧分布や吸収係数分布）を表す画像データに変換され、表示部１０９を介して出力される。本実施形態では、データ処理部１０８が吸収係数分布を表す三次元データをまず生成し、当該三次元データを二次元画像に変換して表示部１０９に出力する。
次に、データ処理部１０８が、取得した三次元の吸収係数分布（以下、三次元データ）
を、二次元画像に変換する方法について説明する。図２は、データ処理部１０８が行う処理のフローチャート図である。
を、二次元画像に変換する方法について説明する。図２は、データ処理部１０８が行う処理のフローチャート図である。
まず、ステップＳ１１で、被検体情報である三次元データを取得する。当該三次元データは、三次元空間上のボクセルによって表現されたデータである。本実施形態では、被検体情報に対応する三次元空間内の座標を（Ｘ，Ｙ，Ｚ）と表す。
次に、ステップＳ１２で、被検体１０３の正面となる方向（正面方向）を決定する。正面方向とは、三次元空間内において、被検体画像を提供するユーザの視点が位置する方向である。当該方向は固定されていてもよいし、ユーザの入力（例えば、画面上において被検体を移動／回転させる操作）に応じて変更されてもよい。
次に、ステップＳ１３で、ユーザの視点を基準とする座標系を設定し、三次元データに対して座標変換を行う。ここでは、設定された正面方向と平行な軸であるＺ軸、当該Ｚ軸と直交するＸ軸およびＹ軸からなる座標系を設定し、三次元データの座標変換を行う。ここでは、Ｚ軸の正方向をユーザの視点側であるものとする。また、変換後の座標系における座標を（Ｘ’，Ｙ’，Ｚ’）と表す。
次に、ステップＳ１４で、変換後の座標系におけるＸ’およびＹ’の組み合わせのそれぞれについて、Ｚ軸における位置を正方向から負方向に変化させながら走査し、最初に値が得られた（すなわち信号が存在する）ボクセルの座標を記録する。これにより、被検体の表面形状を表すデータ（以下、形状データ）を得ることができる。
ここで、被検体表面の座標を得る方法の具体例を、図３を参照しながら説明する。
符号３０１〜３０４は、Ｘ’−Ｙ’平面をＺ’軸方向に動かす様子を示したものである。図３の例では、符号３０１から符号３０４へ進むにつれてＺ軸の座標が小さくなる。
符号３０１〜３０４は、Ｘ’−Ｙ’平面をＺ’軸方向に動かす様子を示したものである。図３の例では、符号３０１から符号３０４へ進むにつれてＺ軸の座標が小さくなる。
ここで、符号３０２の位置において、ボクセル３０５の位置に信号が観測されたものとする。この場合、ボクセル３０５の座標にＺ’座標を加えたものを、三次元空間中における被検体表面の座標であると判定する。この処理を、Ｘ’−Ｙ’平面上のボクセル全てについて行うことで、被検体表面の立体的な形状データを取得することができる。
ただし、あるボクセルに信号が観測されたとしても、隣接するボクセルに信号が観測されない場合（例えばボクセル３０６）、当該ボクセルについては信号を観測していないものとして扱う。これは、ノイズを考慮した処理である。すなわち、隣接ボクセルと連続的に信号が観測された位置を、被検体の表面とみなす。例えば、データ処理部１０８は、隣接するボクセル間のボクセル値の差が所定の値よりも小さい場合に、対応する位置を被検体の表面とみなす処理を行うことができる。
なお、ここでは光音響信号を用いて被検体の形状データを取得する例を挙げたが、他の手段によって形状データを取得してもよい。例えば、超音波、核磁気共鳴、Ｘ線などを用いて測定されたデータを、被検体の形状データとして利用してもよい。
なお、データ処理部１０８が、取得した形状データが示す形状が光音響信号を取得したときの被検体の形状と異なる場合、形状データを補正してもよい。すなわち、データ処理部１０８が、取得した形状データの示す形状と光音響信号を取得したときの被検体の形状とが対応するように、形状データを補正してもよい。例えば、データ処理部１０８が、形状データの示す形状をアフィン変換等にしたがって変形させ、変形後の形状を示す形状データを取得してもよい。
なお、データ処理部１０８が、取得した形状データが示す形状が光音響信号を取得したときの被検体の形状と異なる場合、形状データを補正してもよい。すなわち、データ処理部１０８が、取得した形状データの示す形状と光音響信号を取得したときの被検体の形状とが対応するように、形状データを補正してもよい。例えば、データ処理部１０８が、形状データの示す形状をアフィン変換等にしたがって変形させ、変形後の形状を示す形状データを取得してもよい。
図２に戻って説明を続ける。
ステップＳ１５で、被検体の内部（表面より内側）に存在する各ボクセルについて、被検体表面までの距離を求める。以下、ここで求めた距離を、ボクセル深度と称する。なお、被検体表面とボクセルとの距離とは、被検体表面における法線方向の距離のことを指す。
ステップＳ１５で、被検体の内部（表面より内側）に存在する各ボクセルについて、被検体表面までの距離を求める。以下、ここで求めた距離を、ボクセル深度と称する。なお、被検体表面とボクセルとの距離とは、被検体表面における法線方向の距離のことを指す。
次に、ステップＳ１６で、各ボクセルが保持している値に対して、対応するボクセル深度を用いてフィルタ処理を行う。ここでは、フィルタの方法として以下の２つの方法を例示する。もちろん、例示したフィルタ以外を適用してもよい。
（１）ボクセル深度が大きいほど、ボクセル値が小さくなるフィルタ処理を行う
当該方法によると、被検体の表面（例えば人の皮膚）から奥行き方向に遠ざかる（すなわち、ユーザの視点から遠ざかる）ほど、ボクセルに対応づいた値が小さくなる（二次元画像に変換した際の輝度が小さくなる）ようなフィルタ処理が行われる。例えば、ボクセル深度が大きいほど、ボクセル値に対する重みを小さくするフィルタ処理を適用してもよい。
（２）ボクセル深度が大きいほど、ボケの強度が強くなるフィルタ処理を行う
当該方法によると、被検体の表面から奥行き方向に遠ざかるほど、ボクセルに対応づいた値がより強く散乱する（二次元画像に変換した際のボケの強度が大きくなる）ようなフィルタ処理が行われる。例えば、ガウスぼかしを適用し、被検体表面からの距離に基づいて、ぼかし半径を変化させることで、被検体表面からの距離が遠いボクセルの画像を強くぼかすように処理する。すなわち、被検体表面からの距離が遠いほど、ガウシアン分布の半径を大きくして、三次元データにガウシアン分布をコンボリューションする。
（１）ボクセル深度が大きいほど、ボクセル値が小さくなるフィルタ処理を行う
当該方法によると、被検体の表面（例えば人の皮膚）から奥行き方向に遠ざかる（すなわち、ユーザの視点から遠ざかる）ほど、ボクセルに対応づいた値が小さくなる（二次元画像に変換した際の輝度が小さくなる）ようなフィルタ処理が行われる。例えば、ボクセル深度が大きいほど、ボクセル値に対する重みを小さくするフィルタ処理を適用してもよい。
（２）ボクセル深度が大きいほど、ボケの強度が強くなるフィルタ処理を行う
当該方法によると、被検体の表面から奥行き方向に遠ざかるほど、ボクセルに対応づいた値がより強く散乱する（二次元画像に変換した際のボケの強度が大きくなる）ようなフィルタ処理が行われる。例えば、ガウスぼかしを適用し、被検体表面からの距離に基づいて、ぼかし半径を変化させることで、被検体表面からの距離が遠いボクセルの画像を強くぼかすように処理する。すなわち、被検体表面からの距離が遠いほど、ガウシアン分布の半径を大きくして、三次元データにガウシアン分布をコンボリューションする。
次に、ステップＳ１７で、フィルタ処理後の三次元データを用いて、所望のレンダリング方法により二次元画像を生成し、表示部１０９に出力する。例えば、レンダリング方法としては、最大値投影法、最小値投映法、ボリュームレンダリング、サーフェイスレンダリングなどの公知のレンダリング方法を採用することができる。
例えば、上記（１）のフィルタ処理後の三次元データを最大値投影法により投影すると、ボクセル深度の大きい位置の透明度が高い二次元画像を表示させることができる。また、上記（２）のフィルタ処理後の三次元データを最大値投影法により投影すると、ボクセル深度の大きい位置でのボケ具合が高い二次元画像を表示させることができる。このようにして、奥行き感を表現した二次元画像を表示させることができる。その結果、医師等の使用者は、奥行き方向の情報を区別しやすい二次元画像を確認して診断を行うことができる。
本実施形態では（１）または（２）のフィルタ処理のように三次元データのボクセル値を変更するフィルタ処理を適用した。ただし、これらのフィルタ処理以外にも、ボクセル深度に応じて透明度、または、ボケ具合を変更することのできるフィルタ処理であれば、いかなる方法を採用してもよい。
なお、フィルタ処理後の三次元データをメモリに保存した後に、メモリに保存されたフィルタ処理後の三次元データを投影することにより二次元画像を生成してもよい。また、三次元データを投影する際にフィルタ処理を適用してもよい。このような処理も、フィルタ処理後の三次元データに基づいて二次元画像を生成するものと等価である。
例えば、上記（１）のフィルタ処理後の三次元データを最大値投影法により投影すると、ボクセル深度の大きい位置の透明度が高い二次元画像を表示させることができる。また、上記（２）のフィルタ処理後の三次元データを最大値投影法により投影すると、ボクセル深度の大きい位置でのボケ具合が高い二次元画像を表示させることができる。このようにして、奥行き感を表現した二次元画像を表示させることができる。その結果、医師等の使用者は、奥行き方向の情報を区別しやすい二次元画像を確認して診断を行うことができる。
本実施形態では（１）または（２）のフィルタ処理のように三次元データのボクセル値を変更するフィルタ処理を適用した。ただし、これらのフィルタ処理以外にも、ボクセル深度に応じて透明度、または、ボケ具合を変更することのできるフィルタ処理であれば、いかなる方法を採用してもよい。
なお、フィルタ処理後の三次元データをメモリに保存した後に、メモリに保存されたフィルタ処理後の三次元データを投影することにより二次元画像を生成してもよい。また、三次元データを投影する際にフィルタ処理を適用してもよい。このような処理も、フィルタ処理後の三次元データに基づいて二次元画像を生成するものと等価である。
ここで、図４に示したような円錐形である被検体４０１に対して光音響測定を行い、生成された三次元データに対してフィルタ処理を行う例を挙げる。被検体４０１には、底面と平行な直線状の光吸収体４本と、表面に沿った円状の光吸収体４つが配置されている。本例では、被検体の頂点側を手前とし、底面側と奥側として光音響測定を行う。
図５に示した画像５０２は、取得した三次元データを、従来の最大値投影法を用いて二次元化した画像である。画像５０２では、すべての線が同じ輝度で描写されていることがわかる。また、画像５０３は、三次元データを、従来のボリュームレンダリングの手法を
用いて二次元化した画像である。画像５０３では、視点からの距離に応じて線の輝度が変化していることがわかる。
用いて二次元化した画像である。画像５０３では、視点からの距離に応じて線の輝度が変化していることがわかる。
画像５０４は、取得した三次元データを、本発明に係る方法によって二次元化した画像である。画像５０４では、最も近い被検体表面からの距離に基づいて、ボクセルに対してぼかしフィルタを適用する。このようしてぼかしフィルタを適用した三次元ボクセルを最大値投映法により二次元画像として表示する。
画像５０５，５０６，５０７は、それぞれ画像５０２，５０３，５０４の中心付近を拡大したものである。図からもわかるように、画像５０５からは、光吸収体の位置関係を把握することはできない。これに対して、画像５０６および５０７では、光吸収体の位置関係を把握することが容易になっている。
同様に、画像５０８，５０９，５１０は、それぞれ画像５０２，５０３，５０４の右下付近を拡大したものである。図からもわかるように、画像５０８および５０９と比較して、画像５１０では、光吸収体の深度（表面からの距離）が把握しやすくなっている。
なお、図５の例では、被検体表面と光吸収体の位置関係を示すため、被検体の表面に近い光吸収体を強調し、被検体の深部にある光吸収体をぼかす処理を行ったが、他の方法も採用可能である。例えば、被検体の深部にある光吸収体に対してより小さい輝度を与え、表面付近にある光吸収体に対してより大きい輝度を与えるようなフィルタ処理を行ってもよい。
図６は、表示部１０９に表示されるユーザインタフェースの例である。当該ユーザインタフェースは、マウス等のポインティングシステムを用いて操作するウィンドウシステムによって表示される。
ウィンドウ６０１には、取得した光音響信号を被検体ごとに一覧表示するリスト表示部６０２と、光音響信号を画像化して表示する画像表示部６０３が配置されている。また、画像表示部における画像の表示方法を変更するための切り替えボタン６０４が配置されている。
画像表示部６０３には、リスト表示部６０２にて選択した被検体に対応する光音響画像が表示される。画像表示部６０３には、切り替えボタン６０４にて選択された画像処理方法によって二次元化された被検体画像が表示される。
本実施形態では、切り替えボタン６０４を押下することで、任意の画像処理方法を選択することができる。選択可能な画像処理方法として、例えば、断層画像の三断面表示、ボリュームレンダリングデータの二次元射影、最大値投影法、本発明に係る表示方法などが挙げられる。
なお、リスト表示部６０２と画像表示部６０３は、同じウィンドウ上に配置してあっても、別のウィンドウ上に配置してあってもよい。また、画像処理方法を選択するためのボタンは、切り替えボタン６０４のようにウィンドウ上に配置されなくてもよい。例えば、ウィンドウに付属するメニュー６０５に配置してもよい。
また、本実施形態では、画像表示部６０３にて表示した画像を、ポインティングデバイスなどを用いて回転させることができる。ここで、任意の向きを正面として決定すると、切り替えボタン６０４にて選択した方法を用いて、三次元データが二次元画像に変換される。また、切り替えボタン６０４によって、本発明に係る方法を選択していた場合、決定
した姿勢に応じた座標変換がステップＳ１３で行われる。
した姿勢に応じた座標変換がステップＳ１３で行われる。
（第二の実施形態）
第二の実施形態は、半球状の音響波探触子を用いて人の乳房に対する測定を行い、乳房内の血管を可視化する実施形態である。また、第二の実施形態に係る被検体情報取得装置は、複数の波長を用いて測定を行うことで、乳房内の酸素飽和度分布を取得する。これにより、乳房内における動脈と静脈の位置を判断することができる。
第二の実施形態は、半球状の音響波探触子を用いて人の乳房に対する測定を行い、乳房内の血管を可視化する実施形態である。また、第二の実施形態に係る被検体情報取得装置は、複数の波長を用いて測定を行うことで、乳房内の酸素飽和度分布を取得する。これにより、乳房内における動脈と静脈の位置を判断することができる。
第二の実施形態では、音響波探触子１０６が、半球状の支持部材に取り付けられている。なお、図７では、音響波探触子が支持部材の一部に取り付けられている例を示したが、音響波探触子１０６は、支持部材上を移動可能な構成であってもよい。また、複数の音響素子を支持部材の内側に分散して配置することで音響波探触子１０６を構成してもよい。
その他の構成は第一の実施形態と同様であるため、説明は省略する。
その他の構成は第一の実施形態と同様であるため、説明は省略する。
光源１０１から、光学系１０２を通じて被検体１０３に相当する乳房７０３にパルス光を照射すると、吸収体１０４に相当する血管内の赤血球７０４で光の吸収が起きる。
赤血球７０４は、酸素化ヘモグロビンと脱酸素化ヘモグロビンを含み、照射される波長ごとに吸収量が異なる。このことを利用して、光源１０１から複数の波長の光を照射することで、被検体内の酸素飽和度の分布を得ることができる。複数の波長とは、例えば、７５６ナノメートルと、７９７ナノメートルの二種類の波長である。
赤血球７０４は、酸素化ヘモグロビンと脱酸素化ヘモグロビンを含み、照射される波長ごとに吸収量が異なる。このことを利用して、光源１０１から複数の波長の光を照射することで、被検体内の酸素飽和度の分布を得ることができる。複数の波長とは、例えば、７５６ナノメートルと、７９７ナノメートルの二種類の波長である。
赤血球７０４から得られた光音響波は、音響整合材１０５を通じて音響波探触子１０６に到達し、電気信号（光音響信号）に変換される。得られた光音響信号は、信号処理部１０７を介してデータ処理部１０８に送信される。
データ処理部１０８では、複数の波長ごとに得られた光音響信号を用いて、ボクセルごとの酸素飽和度を計算する。
また、乳房７０３の形状データを取得した後、光音響信号が検出されたボクセルと被検体表面との距離を計算し、得られた距離に基づいて、各ボクセルに対するフィルタ処理を行なう。具体的には、皮膚表面からの距離が大きいボクセルは透明化処理を行う。フィルタ処理後の画像は、奥行き方向に重なるボクセルでそれぞれ信号が有る場合には、より深い位置のボクセルが透明化されているために、被検体表面に近いボクセルのみを投影されることで、立体感のある二次元画像として表示部１０９に出力される。
データ処理部１０８では、複数の波長ごとに得られた光音響信号を用いて、ボクセルごとの酸素飽和度を計算する。
また、乳房７０３の形状データを取得した後、光音響信号が検出されたボクセルと被検体表面との距離を計算し、得られた距離に基づいて、各ボクセルに対するフィルタ処理を行なう。具体的には、皮膚表面からの距離が大きいボクセルは透明化処理を行う。フィルタ処理後の画像は、奥行き方向に重なるボクセルでそれぞれ信号が有る場合には、より深い位置のボクセルが透明化されているために、被検体表面に近いボクセルのみを投影されることで、立体感のある二次元画像として表示部１０９に出力される。
赤血球７０４は、血管を流れる血液中に存在しているため、乳房内の赤血球を画像化することで血管像を描出することができる。また、複数の波長を用いて酸素飽和度を計算すると、血管内の酸素飽和度分布を得ることができるため、関心領域内の血管が動脈であるか、静脈であるかを把握できるようになる。
光音響測定によって得られた画像上において血管の動静脈の区別をする場合、血管の連続性が視認しやすいことが重要となる。本発明に係る方法によって二次元画像を生成することで、被検体組織の深さと血管の連続性、血管同士の重なりを判別することが可能となる。すなわち、画像化された血管の動静脈の区別が容易になる。
なお、本実施形態において、複数の波長に対応する複数の三次元データのそれぞれにフィルタ処理を実行した後に、フィルタ処理後の複数の三次元データに基づいて酸素飽和度分布を算出してもよい。このとき、複数の三次元データの夫々に異なるフィルタ処理を実行してもよい。例えば、短波長側の波長ではメラニンの吸収が強くなり、深部で発生する音響波の振幅が小さくなるため、深部の画質が低下する場合がある。そこで、長波長側の波長に比べて、短波長側の波長に対応する三次元データに対して、表面からの距離が大きくなるにつれてボケの強度を強くしたり、輝度を弱くしたりしてもよい。これにより、波長間における三次元データの画質の違いを小さくすることができる。
光音響測定によって得られた画像上において血管の動静脈の区別をする場合、血管の連続性が視認しやすいことが重要となる。本発明に係る方法によって二次元画像を生成することで、被検体組織の深さと血管の連続性、血管同士の重なりを判別することが可能となる。すなわち、画像化された血管の動静脈の区別が容易になる。
なお、本実施形態において、複数の波長に対応する複数の三次元データのそれぞれにフィルタ処理を実行した後に、フィルタ処理後の複数の三次元データに基づいて酸素飽和度分布を算出してもよい。このとき、複数の三次元データの夫々に異なるフィルタ処理を実行してもよい。例えば、短波長側の波長ではメラニンの吸収が強くなり、深部で発生する音響波の振幅が小さくなるため、深部の画質が低下する場合がある。そこで、長波長側の波長に比べて、短波長側の波長に対応する三次元データに対して、表面からの距離が大きくなるにつれてボケの強度を強くしたり、輝度を弱くしたりしてもよい。これにより、波長間における三次元データの画質の違いを小さくすることができる。
なお、第二の実施形態に、被検体の形状を測定する手段を追加してもよい。例えば、超音波を被検体に送信し、反射波に基づいて当該被検体の形状を取得する検出装置７１３を追加してもよい。
（その他の実施形態）
なお、各実施形態の説明は本発明を説明する上での例示であり、本発明は、発明の趣旨を逸脱しない範囲で適宜変更または組み合わせて実施することができる。
例えば、本発明は、上記処理の少なくとも一部を実施する被検体情報取得装置として実施することもできる。また、本発明は、上記処理の少なくとも一部を含む被検体情報取得方法として実施することもできる。
また、本発明は、被検体情報取得装置が取得した三次元データに対して画像処理を行う画像処理装置として実施することもできる。また、本発明は、当該画像処理装置が行う画像処理方法として実施することもできる。
また、上記方法をコンピュータに実行させるためのプログラムや、当該プログラムが記憶された非一時的記憶媒体として実施することもできる。
上記処理や手段は、技術的な矛盾が生じない限りにおいて、自由に組み合わせて実施することができる。
なお、各実施形態の説明は本発明を説明する上での例示であり、本発明は、発明の趣旨を逸脱しない範囲で適宜変更または組み合わせて実施することができる。
例えば、本発明は、上記処理の少なくとも一部を実施する被検体情報取得装置として実施することもできる。また、本発明は、上記処理の少なくとも一部を含む被検体情報取得方法として実施することもできる。
また、本発明は、被検体情報取得装置が取得した三次元データに対して画像処理を行う画像処理装置として実施することもできる。また、本発明は、当該画像処理装置が行う画像処理方法として実施することもできる。
また、上記方法をコンピュータに実行させるためのプログラムや、当該プログラムが記憶された非一時的記憶媒体として実施することもできる。
上記処理や手段は、技術的な矛盾が生じない限りにおいて、自由に組み合わせて実施することができる。
また、実施形態の説明では、対象のボクセルが被検体の表面に近いほど視認性を高める（値を大きく補正する）フィルタ処理を例示したが、被検体の表面との距離に基づいてフィルタの強度を決定するものであれば、他のフィルタ処理を行ってもよい。例えば、被検体の深部であるほど値を大きく補正してもよい。
また、本発明における「複数のボクセル」は、取得した三次元データの全領域である必要はない。例えば、取得した三次元データを構成する一部のボクセルに対してのみ、実施形態で説明したようなフィルタ処理を行うようにしてもよい。また、残りのボクセルについては、フィルタ処理の結果を反映させる（例えば補間などを行う）ようにしてもよい。
本発明は、以下の処理を実行することによっても実現される。即ち、上述した各実施形態の１以上の機能を実現するプログラムを、ネットワーク又は各種記憶媒体を介してシステム或いは装置に供給し、そのシステム或いは装置のコンピュータにおける一つ以上のプロセッサがプログラムを読み出して実行する処理でも実現可能である。また、１以上の機能を実現する回路（例えば、ＦＰＧＡやＡＳＩＣ）によっても実現可能である。
１０７：信号処理部、１０８：データ処理部
Claims (15)
- 被検体の特性情報を複数のボクセルによって表した三次元データを取得する情報取得手段と、
前記被検体の表面形状に関する情報を取得する形状情報取得手段と、
前記表面形状に関する情報に基づいて、前記ボクセルのそれぞれについて、前記被検体の表面から当該ボクセルに対応する前記被検体内の位置までの距離を算出する距離算出手段と、
前記ボクセルのそれぞれに対して、前記算出した距離に応じたぼかし処理を含むフィルタ処理を行うフィルタ手段と、
前記フィルタ処理後の前記三次元データに基づいて二次元画像を生成する画像生成手段と、
を有することを特徴とする、画像処理装置。 - 前記フィルタ手段は、前記ボクセルのそれぞれに対して、前記被検体の表面からの距離が大きいほどボケの強度を強くする前記フィルタ処理を行う
ことを特徴とする、請求項１に記載の画像処理装置。 - 前記画像生成手段は、前記フィルタ処理後の前記三次元データを用いて、投影により前記二次元画像を生成する
ことを特徴とする、請求項１に記載の画像処理装置。 - 前記形状情報取得手段は、隣接ボクセルと連続的に信号が観測された位置を、被検体の表面とみなす処理によって、前記被検体の表面形状に関する情報を取得する
ことを特徴とする、請求項１に記載の画像処理装置。 - 前記三次元データは、前記被検体への光照射により発生した音響波に基づいて生成されたデータである
ことを特徴とする、請求項１に記載の画像処理装置。 - 前記三次元データは、複数の波長の光照射に対応する複数の三次元データを含み、
前記フィルタ手段は、前記複数の三次元データのそれぞれに互いに異なるフィルタ処理を行う
ことを特徴とする、請求項５に記載の画像処理装置。 - 前記フィルタ手段は、前記複数の三次元データに対して、光照射の波長が小さいほどボケの強度を強くするフィルタ処理を行う
ことを特徴とする、請求項６に記載の画像処理装置。 - 被検体の特性情報を複数のボクセルによって表した三次元データを取得する情報取得ステップと、
前記被検体の表面形状に関する情報を取得する形状情報取得ステップと、
前記表面形状に関する情報に基づいて、前記ボクセルのそれぞれについて、前記被検体の表面から当該ボクセルに対応する前記被検体内の位置までの距離を算出する距離算出ステップと、
前記ボクセルのそれぞれに対して、前記算出した距離に応じたぼかし処理を含むフィルタ処理を行うフィルタステップと、
前記フィルタ処理後の前記三次元データに基づいて二次元画像を生成する画像生成ステップと、
を含むことを特徴とする、画像処理方法。 - 前記フィルタステップでは、前記ボクセルのそれぞれに対して、前記被検体の表面からの距離が大きいほどボケの強度を強くする前記フィルタ処理を行う
ことを特徴とする、請求項８に記載の画像処理方法。 - 前記画像生成ステップでは、前記フィルタ処理後の前記三次元データを用いて、
投影により前記二次元画像を生成する
ことを特徴とする、請求項８に記載の画像処理方法。 - 前記形状情報取得ステップでは、隣接ボクセルと連続的に信号が観測された位置を、被検体の表面とみなす処理によって、前記被検体の表面形状に関する情報を取得する
ことを特徴とする、請求項８に記載の画像処理方法。 - 前記三次元データは、前記被検体への光照射により発生した音響波に基づいて生成されたデータである
ことを特徴とする、請求項８に記載の画像処理方法。 - 前記三次元データは、複数の波長の光照射に対応する複数の三次元データを含み、
前記フィルタステップでは、前記複数の三次元データのそれぞれに互いに異なるフィルタ処理を行う
ことを特徴とする、請求項１２に記載の画像処理方法。 - 前記フィルタステップでは、前記複数の三次元データに対して、光照射の波長が小さいほどボケの強度を強くするフィルタ処理を行う
ことを特徴とする、請求項１３に記載の画像処理方法。 - 被検体の特性情報を複数のボクセルによって表した三次元データを取得する情報取得ステップと、
前記被検体の表面形状に関する情報を取得する形状情報取得ステップと、
前記表面形状に関する情報に基づいて、前記ボクセルのそれぞれについて、前記被検体の表面から当該ボクセルに対応する前記被検体内の位置までの距離を算出する距離算出ステップと、
前記ボクセルのそれぞれに対して、前記算出した距離に応じたぼかし処理を含むフィルタ処理を行うフィルタステップと、
前記フィルタ処理後の前記三次元データに基づいて二次元画像を生成する画像生成ステップと、
を含む画像処理方法をコンピュータに実行させるためのプログラムが記録された非一時的記憶媒体。
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2017061164 | 2017-03-27 | ||
| JP2017061164 | 2017-03-27 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JP2018161467A true JP2018161467A (ja) | 2018-10-18 |
Family
ID=63581146
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2018045346A Pending JP2018161467A (ja) | 2017-03-27 | 2018-03-13 | 画像処理装置および画像処理方法 |
Country Status (2)
| Country | Link |
|---|---|
| US (1) | US10438382B2 (ja) |
| JP (1) | JP2018161467A (ja) |
Family Cites Families (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7447344B2 (en) * | 2004-04-16 | 2008-11-04 | Siemens Medical Solutions Usa, Inc. | System and method for visualization of pulmonary emboli from high-resolution computed tomography images |
| JP4832927B2 (ja) * | 2006-03-14 | 2011-12-07 | オリンパスメディカルシステムズ株式会社 | 医療用画像処理装置及び医療用画像処理方法 |
| JP5599709B2 (ja) * | 2007-09-03 | 2014-10-01 | コーニンクレッカ フィリップス エヌ ヴェ | ボクセルデータの可視化 |
| EP2312534A3 (en) * | 2009-10-15 | 2011-07-06 | Hitachi Aloka Medical, Ltd. | Ultrasonic volume data processing device |
| JP5518096B2 (ja) * | 2009-12-17 | 2014-06-11 | キヤノン株式会社 | 測定システム、画像形成方法及びプログラム |
| GB2478329B (en) * | 2010-03-03 | 2015-03-04 | Samsung Electronics Co Ltd | Medical image processing |
| US20130109963A1 (en) * | 2011-10-31 | 2013-05-02 | The University Of Connecticut | Method and apparatus for medical imaging using combined near-infrared optical tomography, fluorescent tomography and ultrasound |
| JP5762995B2 (ja) | 2012-02-28 | 2015-08-12 | 富士フイルム株式会社 | 光音響画像生成装置及び方法 |
-
2018
- 2018-03-12 US US15/917,951 patent/US10438382B2/en not_active Expired - Fee Related
- 2018-03-13 JP JP2018045346A patent/JP2018161467A/ja active Pending
Also Published As
| Publication number | Publication date |
|---|---|
| US20180276855A1 (en) | 2018-09-27 |
| US10438382B2 (en) | 2019-10-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP6504826B2 (ja) | 情報処理装置および情報処理方法 | |
| US10531798B2 (en) | Photoacoustic information acquiring apparatus and processing method | |
| JP2010088627A (ja) | 生体情報処理装置および生体情報処理方法 | |
| JP6335612B2 (ja) | 光音響装置、処理装置、処理方法、及びプログラム | |
| JP2013255697A (ja) | 被検体情報取得装置およびその制御方法 | |
| US20170325693A1 (en) | Photoacoustic apparatus and control method of photoacoustic apparatus | |
| JP6732476B2 (ja) | 被検体情報取得装置 | |
| US20170265750A1 (en) | Information processing system and display control method | |
| EP3329843B1 (en) | Display control apparatus, display control method, and program | |
| KR101899838B1 (ko) | 광음향 장치 및 정보 취득장치 | |
| US20180228377A1 (en) | Object information acquiring apparatus and display method | |
| JP6486056B2 (ja) | 光音響装置および光音響装置の処理方法 | |
| JP6469133B2 (ja) | 処理装置、光音響装置、処理方法、およびプログラム | |
| JP2018126389A (ja) | 情報処理装置、情報処理方法、およびプログラム | |
| WO2018207692A1 (en) | Display control apparatus, image display method, and program | |
| JP6776115B2 (ja) | 処理装置および処理方法 | |
| JP2019088346A (ja) | 光音響装置および被検体情報取得方法 | |
| US20200275840A1 (en) | Information-processing apparatus, method of processing information, and medium | |
| US20190099083A1 (en) | Information processing apparatus and information processing method | |
| US20200085345A1 (en) | Object information acquisition apparatus and method of controlling the same | |
| JP2018161467A (ja) | 画像処理装置および画像処理方法 | |
| WO2018230409A1 (ja) | 情報処理装置、情報処理方法、及びプログラム | |
| JP2019083887A (ja) | 情報処理装置および情報処理方法 | |
| US11599992B2 (en) | Display control apparatus, display method, and non-transitory storage medium | |
| WO2017222033A1 (en) | Apparatus and information processing method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| RD02 | Notification of acceptance of power of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7422 Effective date: 20181116 |