JP7479684B2 - 画像生成方法、画像生成プログラムおよび画像生成装置 - Google Patents
画像生成方法、画像生成プログラムおよび画像生成装置 Download PDFInfo
- Publication number
- JP7479684B2 JP7479684B2 JP2020138622A JP2020138622A JP7479684B2 JP 7479684 B2 JP7479684 B2 JP 7479684B2 JP 2020138622 A JP2020138622 A JP 2020138622A JP 2020138622 A JP2020138622 A JP 2020138622A JP 7479684 B2 JP7479684 B2 JP 7479684B2
- Authority
- JP
- Japan
- Prior art keywords
- image
- interest
- region
- sound speed
- reconstructing
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims description 163
- PWPJGUXAGUPAHP-UHFFFAOYSA-N lufenuron Chemical compound C1=C(Cl)C(OC(F)(F)C(C(F)(F)F)F)=CC(Cl)=C1NC(=O)NC(=O)C1=C(F)C=CC=C1F PWPJGUXAGUPAHP-UHFFFAOYSA-N 0.000 title 1
- 238000005259 measurement Methods 0.000 claims description 118
- 238000012360 testing method Methods 0.000 claims description 80
- 230000008569 process Effects 0.000 claims description 65
- 230000008859 change Effects 0.000 claims description 53
- 238000004364 calculation method Methods 0.000 claims description 51
- 238000012545 processing Methods 0.000 claims description 33
- 238000003384 imaging method Methods 0.000 description 52
- 239000000463 material Substances 0.000 description 29
- 238000003860 storage Methods 0.000 description 18
- 238000012986 modification Methods 0.000 description 16
- 230000004048 modification Effects 0.000 description 16
- 210000004204 blood vessel Anatomy 0.000 description 14
- 238000009826 distribution Methods 0.000 description 13
- 238000000926 separation method Methods 0.000 description 13
- 230000007246 mechanism Effects 0.000 description 12
- 230000006870 function Effects 0.000 description 11
- 238000002604 ultrasonography Methods 0.000 description 10
- 238000001914 filtration Methods 0.000 description 9
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 9
- 238000001514 detection method Methods 0.000 description 8
- 230000000694 effects Effects 0.000 description 8
- 230000003287 optical effect Effects 0.000 description 8
- 238000005516 engineering process Methods 0.000 description 7
- 238000010521 absorption reaction Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 6
- 238000011160 research Methods 0.000 description 6
- 238000003325 tomography Methods 0.000 description 6
- 238000004458 analytical method Methods 0.000 description 5
- 230000003247 decreasing effect Effects 0.000 description 5
- 210000001519 tissue Anatomy 0.000 description 5
- 230000007423 decrease Effects 0.000 description 4
- 238000010895 photoacoustic effect Methods 0.000 description 4
- 238000005070 sampling Methods 0.000 description 4
- 206010028980 Neoplasm Diseases 0.000 description 3
- 239000006096 absorbing agent Substances 0.000 description 3
- 201000011510 cancer Diseases 0.000 description 3
- 239000012530 fluid Substances 0.000 description 3
- 230000004044 response Effects 0.000 description 3
- 210000004872 soft tissue Anatomy 0.000 description 3
- 239000007787 solid Substances 0.000 description 3
- 239000002033 PVDF binder Substances 0.000 description 2
- 239000004698 Polyethylene Substances 0.000 description 2
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 2
- 238000003745 diagnosis Methods 0.000 description 2
- 230000001678 irradiating effect Effects 0.000 description 2
- 229910052451 lead zirconate titanate Inorganic materials 0.000 description 2
- 239000007788 liquid Substances 0.000 description 2
- 239000001301 oxygen Substances 0.000 description 2
- 229910052760 oxygen Inorganic materials 0.000 description 2
- 229920000573 polyethylene Polymers 0.000 description 2
- -1 polyethylene terephthalate Polymers 0.000 description 2
- 229920000139 polyethylene terephthalate Polymers 0.000 description 2
- 239000005020 polyethylene terephthalate Substances 0.000 description 2
- 229920002981 polyvinylidene fluoride Polymers 0.000 description 2
- 238000002360 preparation method Methods 0.000 description 2
- 238000003825 pressing Methods 0.000 description 2
- 238000003672 processing method Methods 0.000 description 2
- 230000000644 propagated effect Effects 0.000 description 2
- 239000000523 sample Substances 0.000 description 2
- 238000012800 visualization Methods 0.000 description 2
- 206010006187 Breast cancer Diseases 0.000 description 1
- 208000026310 Breast neoplasm Diseases 0.000 description 1
- 241001465754 Metazoa Species 0.000 description 1
- 239000000654 additive Substances 0.000 description 1
- 230000000996 additive effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000036760 body temperature Effects 0.000 description 1
- 210000000481 breast Anatomy 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 239000003086 colorant Substances 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000002316 cosmetic surgery Methods 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 239000007789 gas Substances 0.000 description 1
- 238000010191 image analysis Methods 0.000 description 1
- 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 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
- 238000012827 research and development Methods 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 230000001960 triggered effect Effects 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
Landscapes
- Ultra Sonic Daignosis Equipment (AREA)
Description
本発明は、画像生成方法、画像生成プログラムおよび画像生成装置に関する。
被検体の所定の被検部位を測定した音響波の信号データに基づいて画像を再構成する画像生成装置が開発されている。
例えば、光源（レーザーなど）からの光を生体などの被検体に照射し、該被検体の内部の情報を可視化する画像生成装置の研究が医療分野で積極的に進められている。このような光による可視化技術の一つとして、Ｐｈｏｔｏａｃｏｕｓｔｉｃ Ｔｏｍｏｇｒａｐｈｙ（ＰＡＴ：光音響トモグラフィー、光超音波トモグラフィーと称することもある）がある。光音響トモグラフィーを用いたイメージング装置では、照射した光が被検体内を伝播し、拡散した光のエネルギーを吸収した光吸収性を有する生体組織から発生される音響波（典型的に超音波である）を、被検体を取り囲む複数の箇所で検出する。そして、得られた信号を数学的に解析処理し、被検体内部の光学特性値、特に吸収係数分布に関連した情報を可視化する。近年、この光音響トモグラフィー装置を用いて、小動物の血管像をイメージングする非臨床研究や、この原理を乳がんなどの画像診断あるいは形成外科領域の術前計画に応用するための臨床研究が積極的に進められている。臨床研究用光音響トモグラフィー装置として、例えば内閣府の革新的研究開発推進プログラム（ＩｍＰＡＣＴプログラム）で研究開発された非特許文献１の装置が、良好な画質の三次元画像を提供できる装置として知られている。
このような光音響イメージング装置や、従来から医療分野で使われている超音波診断装置（生体内で反射した音響波を検出し、画像を再構成する装置）では、通常、被検体の平均的な音速（被検体内部における音響波の伝播速度）を用いて、画像が再構成される。一般的に、音響波の伝播速度は経験的な値や文献値などにより決定される。しかしながら、伝播速度には個体差があり、また、伝播速度は同一被検体であったとしても被検体の体温など撮影時の状況により変化する。そのため、画像再構成時に用いた速度と実際の伝播速度が異なる場合、装置固有の分解能が得られないだけでなく、画像コントラストも低下する、形状の再現性が低下するなどの現象が起き、画質が大きく劣化するという課題があった。
このような課題の解決方法は、例えば特許文献１に開示されている。特許文献１に開示の技術では、超音波伝搬速度を変更し得る速度値変更手段と、この手段による速度値変更に従って送受信回路における遅延時間を修正する遅延時間修正手段とを具備することで、最適音速に設定することが可能である。
しかしながら、特許文献１に開示の技術では、超音波Ｂモード画像を撮影しながら速度値変更手段であるボタンを押して最適値を探索する手法となる。一般にプローブを把持して画像を取得する超音波診断装置では超音波診断技師の技量によって画像が左右されることが課題として挙げられる。それに加えて音速も撮影者の技量に依存する構成とすると、さらなる画質の技師による差異が増加することが課題となる。
別の解決方法は、例えば特許文献２に開示されている。特許文献２に開示の技術では、整相された複数の受信信号に基づいて設定音速と実際の媒質音速との誤差を推定し、その誤差に基づいて媒質音速を推定する処理を被検体の複数の計測領域について行なう。この計測により得られた複数の媒質音速推定値について確からしさを評価し、最も評価の高い媒質音速推定値を複数の計測領域に対する共通の設定音速として設定する。これによって、整相処理に用いる設定音速を実際の媒質音速と同一又は近い値に自動的に調整することが可能である。
自動調整であれば技師依存の課題は解決されるが、特許文献２に開示の技術は、いわゆるイテレーション処理と称される反復計算技術である。このような技術は、一般に計算リソースが大きいため、リアルタイム処理には不向きである。そのため超音波検査技師が患者を超音波検査しながら、その場で最適音速の画像を得るのは難しい。
また本発明者らの検討により、以下の知見が得られている。
非特許文献１に記載した光音響イメージング装置では、撮影を行いながらリアルタイムで画像再構成が行われ、画像が操作パネル上に表示される。この画像再構成の際に、計算に用いる音速が実際の値とずれていると、フォーカスが合わずに、ぼやけた画像になってしまう。この装置では、被検体からセンサに至るまでの音響伝搬媒体（音響整合材）として水を使用している。このため、撮影時の温度を計測して、その温度に対応した水の音速を計算に用いることによって、おおよそ適切な音速が設定される。
しかしながら、この装置の持つ分解能はサブミリメートルオーダーである。より正確に音速を調整しないと、装置が持つ解像度の画像を提供することができない。すなわち、皮下直下の血管の画像を再構成する場合であっても、単に水の音速を適用するだけでは、ぼやけた画像が再構成されてしまう。この理由として、水の層の中で温度分布が存在する影響で、温度センサの値が全体の水の平均水温を反映しない可能性がある。または、皮膚の音速が水の音速と異なるために、音速の違いが皮下の血管のイメージング性能に影響する可能性がある。
そのため、この装置を用いた臨床研究を行う場合には、いわゆる「オフライン画像再構成」を行っている。すなわち、撮影しながらリアルタイムで表示される画像を参照するとともに、取得された光音響信号の生データを保存する。一通りの撮影が終了してから、音速を変えて画像再構成する。このオフライン画像再構成では、正確な音速は事前にわからない。このため、複数の音速を設定し、バッチ処理による画像再構成を行い、事後的に最適な音速を知り画像を得る。
しかしながら、例えば画像のサイズを非特許文献１に記載の最大サイズである横１８０ｍｍ、縦２７０ｍｍの範囲とし、画像再構成の最小単位である立方体（ボクセル）の一辺を０．１２５ｍｍとした場合、一つの音速条件で、画像再構成の時間は約１０分以上を要する。したがって、複数の音速条件をバッチ処理する場合は、画像再構成の時間は数時間を要することもある。
そして、被検体の位置によって最適な音速が異なるので、診断や治療に特に影響する関心領域で適した音速が用いられることが期待される。また、皮下直下で最適な音速が設定されたとしても、例えば皮下１ｃｍ以上の深部では別の音速が最適となることがある。一つの音速条件では、被検体全体として音速が最適だと言えない場合がある。このような場合では、一般的なカメラにおける被写界深度が狭い場合と同様の現象が生じる可能性がある。例えば、生体の深部方向に対して限定された深さのみでフォーカスが合い、それ以外の深さではぼやけてしまうという現象が発生する。
したがって、最適な音速をいち早く求めることが望まれる。また、被検体の関心領域の画像再構成が適切な音速を用いて行われ、高画質の画像を再現できることが望まれる。
Ｎａｇａｅ Ｋ， Ａｓａｏ Ｙ， Ｓｕｄｏ Ｙ ｅｔ ａｌ． Ｒｅａｌ－ｔｉｍｅ ３Ｄ Ｐｈｏｔｏａｃｏｕｓｔｉｃ Ｖｉｓｕａｌｉｚａｔｉｏｎ Ｓｙｓｔｅｍ ｗｈｉｔｈ ａ Ｗｉｄｅ Ｆｉｅｌｄ ｏｆ Ｖｉｅｗ ｆｗｏｒ Ｉｍａｇｉｎｇ Ｈｕｍａｎ Ｌｉｍｂｓ． Ｆ１０００Ｒｅｓｅａｒｃｈ ２０１９， ７：１８１３
本発明の目的は、関心領域について画質が特に良好な音響画像を早く再構成することである。
本発明の一態様によれば、
所定の被検部位を測定した音響波の信号データに基づいて、音響画像を生成する画像生成方法であって、
前記信号データに基づいて、音響画像からなる仮の測定画像を再構成する工程と、
前記仮の測定画像内で部分的な関心領域の指定を受け付ける工程と、
前記関心領域において、音速パラメータを変更しながら、当該変更された音速パラメータに基づいて、音響画像からなる関心領域画像を順次再構成する工程と、
前記関心領域画像に基づいて前記音速パラメータの適切な値を決定する工程と、
前記音速パラメータの前記適切な値に基づいて、前記関心領域よりも広い領域の音響画像を再構成する工程と、
を有し、
前記仮の測定画像を再構成する工程では、
前記仮の測定画像を再構成するための単位面積当たりの計算量を、前記関心領域画像を再構成するための単位面積当たりの計算量よりも少なくする
画像生成方法が提供される。
所定の被検部位を測定した音響波の信号データに基づいて、音響画像を生成する画像生成方法であって、
前記信号データに基づいて、音響画像からなる仮の測定画像を再構成する工程と、
前記仮の測定画像内で部分的な関心領域の指定を受け付ける工程と、
前記関心領域において、音速パラメータを変更しながら、当該変更された音速パラメータに基づいて、音響画像からなる関心領域画像を順次再構成する工程と、
前記関心領域画像に基づいて前記音速パラメータの適切な値を決定する工程と、
前記音速パラメータの前記適切な値に基づいて、前記関心領域よりも広い領域の音響画像を再構成する工程と、
を有し、
前記仮の測定画像を再構成する工程では、
前記仮の測定画像を再構成するための単位面積当たりの計算量を、前記関心領域画像を再構成するための単位面積当たりの計算量よりも少なくする
画像生成方法が提供される。
本発明の他の態様によれば、
所定の被検部位を測定した音響波の信号データに基づいて、音響画像を生成する画像生成方法であって、
前記被検部位の撮影領域内で部分的な関心領域の指定を受け付ける工程と、
前記関心領域において、音速パラメータを変更しながら、当該変更された音速パラメータに基づいて、音響画像からなる関心領域画像を順次再構成する工程と、
前記関心領域画像に基づいて前記音速パラメータの適切な値を決定する工程と、
前記音速パラメータの前記適切な値に基づいて、前記関心領域よりも広い領域の音響画像を再構成する工程と、
を有し、
前記関心領域の指定を受け付ける工程では、
複数の関心領域の指定を受け付け、
前記関心領域画像を再構成する工程では、
前記複数の関心領域において、共通に設定された音速パラメータに基づいて、複数の関心領域画像を再構成し、
前記音速パラメータの前記適切な値を決定する工程では、
前記複数の関心領域において、共通する音速パラメータの適切な値を決定する
画像生成方法が提供される。
所定の被検部位を測定した音響波の信号データに基づいて、音響画像を生成する画像生成方法であって、
前記被検部位の撮影領域内で部分的な関心領域の指定を受け付ける工程と、
前記関心領域において、音速パラメータを変更しながら、当該変更された音速パラメータに基づいて、音響画像からなる関心領域画像を順次再構成する工程と、
前記関心領域画像に基づいて前記音速パラメータの適切な値を決定する工程と、
前記音速パラメータの前記適切な値に基づいて、前記関心領域よりも広い領域の音響画像を再構成する工程と、
を有し、
前記関心領域の指定を受け付ける工程では、
複数の関心領域の指定を受け付け、
前記関心領域画像を再構成する工程では、
前記複数の関心領域において、共通に設定された音速パラメータに基づいて、複数の関心領域画像を再構成し、
前記音速パラメータの前記適切な値を決定する工程では、
前記複数の関心領域において、共通する音速パラメータの適切な値を決定する
画像生成方法が提供される。
本発明のさらに他の態様によれば、
所定の被検部位を測定した音響波の信号データに基づいて、音響画像を生成する画像生成方法であって、
前記被検部位の撮影領域内で部分的な関心領域の指定を受け付ける工程と、
前記関心領域において、音速パラメータを変更しながら、当該変更された音速パラメータに基づいて、音響画像からなる関心領域画像を順次再構成する工程と、
前記関心領域画像に基づいて前記音速パラメータの適切な値を決定する工程と、
前記音速パラメータの前記適切な値に基づいて、前記関心領域よりも広い領域の音響画像を再構成する工程と、
を有し、
前記関心領域の指定を受け付ける工程では、
複数の関心領域の指定を受け付け、
前記関心領域画像を再構成する工程では、
前記複数の関心領域において、それぞれ個別に設定された音速パラメータに基づいて、複数の関心領域画像を再構成する
画像生成方法が提供される。
所定の被検部位を測定した音響波の信号データに基づいて、音響画像を生成する画像生成方法であって、
前記被検部位の撮影領域内で部分的な関心領域の指定を受け付ける工程と、
前記関心領域において、音速パラメータを変更しながら、当該変更された音速パラメータに基づいて、音響画像からなる関心領域画像を順次再構成する工程と、
前記関心領域画像に基づいて前記音速パラメータの適切な値を決定する工程と、
前記音速パラメータの前記適切な値に基づいて、前記関心領域よりも広い領域の音響画像を再構成する工程と、
を有し、
前記関心領域の指定を受け付ける工程では、
複数の関心領域の指定を受け付け、
前記関心領域画像を再構成する工程では、
前記複数の関心領域において、それぞれ個別に設定された音速パラメータに基づいて、複数の関心領域画像を再構成する
画像生成方法が提供される。
本発明のさらに他の態様によれば、
所定の被検部位を測定した音響波の信号データに基づいて、音響画像を生成する画像生成プログラムであって、
前記信号データに基づいて、音響画像からなる仮の測定画像を再構成する手順と、
前記仮の測定画像内で部分的な関心領域の指定を受け付ける手順と、
前記関心領域において、音速パラメータを変更しながら、当該変更された音速パラメータに基づいて、音響画像からなる関心領域画像を順次再構成する手順と、
前記関心領域画像に基づいて前記音速パラメータの適切な値を決定する手順と、
前記音速パラメータの前記適切な値に基づいて、前記関心領域よりも広い領域の音響画像を再構成する手順と、
をコンピュータに実行させ、
前記仮の測定画像を再構成する手順では、
前記仮の測定画像を再構成するための単位面積当たりの計算量を、前記関心領域画像を再構成するための単位面積当たりの計算量よりも少なくする
画像生成プログラムが提供される。
所定の被検部位を測定した音響波の信号データに基づいて、音響画像を生成する画像生成プログラムであって、
前記信号データに基づいて、音響画像からなる仮の測定画像を再構成する手順と、
前記仮の測定画像内で部分的な関心領域の指定を受け付ける手順と、
前記関心領域において、音速パラメータを変更しながら、当該変更された音速パラメータに基づいて、音響画像からなる関心領域画像を順次再構成する手順と、
前記関心領域画像に基づいて前記音速パラメータの適切な値を決定する手順と、
前記音速パラメータの前記適切な値に基づいて、前記関心領域よりも広い領域の音響画像を再構成する手順と、
をコンピュータに実行させ、
前記仮の測定画像を再構成する手順では、
前記仮の測定画像を再構成するための単位面積当たりの計算量を、前記関心領域画像を再構成するための単位面積当たりの計算量よりも少なくする
画像生成プログラムが提供される。
本発明のさらに他の態様によれば、
所定の被検部位を測定した音響波の信号データに基づいて、音響画像を生成する画像生成プログラムであって、
前記被検部位の撮影領域内で部分的な関心領域の指定を受け付ける手順と、
前記関心領域において、音速パラメータを変更しながら、当該変更された音速パラメータに基づいて、音響画像からなる関心領域画像を順次再構成する手順と、
前記関心領域画像に基づいて前記音速パラメータの適切な値を決定する手順と、
前記音速パラメータの前記適切な値に基づいて、前記関心領域よりも広い領域の音響画像を再構成する手順と、
をコンピュータに実行させ、
前記関心領域の指定を受け付ける手順では、
複数の関心領域の指定を受け付け、
前記関心領域画像を再構成する手順では、
前記複数の関心領域において、共通に設定された音速パラメータに基づいて、複数の関心領域画像を再構成し、
前記音速パラメータの前記適切な値を決定する手順では、
前記複数の関心領域において、共通する音速パラメータの適切な値を決定する
画像生成プログラムが提供される。
所定の被検部位を測定した音響波の信号データに基づいて、音響画像を生成する画像生成プログラムであって、
前記被検部位の撮影領域内で部分的な関心領域の指定を受け付ける手順と、
前記関心領域において、音速パラメータを変更しながら、当該変更された音速パラメータに基づいて、音響画像からなる関心領域画像を順次再構成する手順と、
前記関心領域画像に基づいて前記音速パラメータの適切な値を決定する手順と、
前記音速パラメータの前記適切な値に基づいて、前記関心領域よりも広い領域の音響画像を再構成する手順と、
をコンピュータに実行させ、
前記関心領域の指定を受け付ける手順では、
複数の関心領域の指定を受け付け、
前記関心領域画像を再構成する手順では、
前記複数の関心領域において、共通に設定された音速パラメータに基づいて、複数の関心領域画像を再構成し、
前記音速パラメータの前記適切な値を決定する手順では、
前記複数の関心領域において、共通する音速パラメータの適切な値を決定する
画像生成プログラムが提供される。
本発明のさらに他の態様によれば、
所定の被検部位を測定した音響波の信号データに基づいて、音響画像を生成する画像生成プログラムであって、
前記被検部位の撮影領域内で部分的な関心領域の指定を受け付ける手順と、
前記関心領域において、音速パラメータを変更しながら、当該変更された音速パラメータに基づいて、音響画像からなる関心領域画像を順次再構成する手順と、
前記関心領域画像に基づいて前記音速パラメータの適切な値を決定する手順と、
前記音速パラメータの前記適切な値に基づいて、前記関心領域よりも広い領域の音響画像を再構成する手順と、
をコンピュータに実行させ、
前記関心領域の指定を受け付ける手順では、
複数の関心領域の指定を受け付け、
前記関心領域画像を再構成する手順では、
前記複数の関心領域において、それぞれ個別に設定された音速パラメータに基づいて、複数の関心領域画像を再構成する
画像生成プログラムが提供される。
所定の被検部位を測定した音響波の信号データに基づいて、音響画像を生成する画像生成プログラムであって、
前記被検部位の撮影領域内で部分的な関心領域の指定を受け付ける手順と、
前記関心領域において、音速パラメータを変更しながら、当該変更された音速パラメータに基づいて、音響画像からなる関心領域画像を順次再構成する手順と、
前記関心領域画像に基づいて前記音速パラメータの適切な値を決定する手順と、
前記音速パラメータの前記適切な値に基づいて、前記関心領域よりも広い領域の音響画像を再構成する手順と、
をコンピュータに実行させ、
前記関心領域の指定を受け付ける手順では、
複数の関心領域の指定を受け付け、
前記関心領域画像を再構成する手順では、
前記複数の関心領域において、それぞれ個別に設定された音速パラメータに基づいて、複数の関心領域画像を再構成する
画像生成プログラムが提供される。
本発明のさらに他の態様によれば、
所定の被検部位を測定した音響波の信号データに基づいて、音響画像を生成する処理部を備え、
前記処理部は、
前記信号データに基づいて、音響画像からなる仮の測定画像を再構成する処理と、
前記仮の測定画像内で部分的な関心領域の指定を受け付ける処理と、
前記関心領域において、音速パラメータを変更しながら、当該変更された音速パラメータに基づいて、音響画像からなる関心領域画像を順次再構成する処理と、
前記関心領域画像に基づいて前記音速パラメータの適切な値を決定する処理と、
前記音速パラメータの前記適切な値に基づいて、前記関心領域よりも広い領域の音響画像を再構成する処理と、
を実行し、
前記仮の測定画像を再構成する処理では、
前記仮の測定画像を再構成するための単位面積当たりの計算量を、前記関心領域画像を再構成するための単位面積当たりの計算量よりも少なくする
画像生成装置が提供される。
所定の被検部位を測定した音響波の信号データに基づいて、音響画像を生成する処理部を備え、
前記処理部は、
前記信号データに基づいて、音響画像からなる仮の測定画像を再構成する処理と、
前記仮の測定画像内で部分的な関心領域の指定を受け付ける処理と、
前記関心領域において、音速パラメータを変更しながら、当該変更された音速パラメータに基づいて、音響画像からなる関心領域画像を順次再構成する処理と、
前記関心領域画像に基づいて前記音速パラメータの適切な値を決定する処理と、
前記音速パラメータの前記適切な値に基づいて、前記関心領域よりも広い領域の音響画像を再構成する処理と、
を実行し、
前記仮の測定画像を再構成する処理では、
前記仮の測定画像を再構成するための単位面積当たりの計算量を、前記関心領域画像を再構成するための単位面積当たりの計算量よりも少なくする
画像生成装置が提供される。
本発明のさらに他の態様によれば、
所定の被検部位を測定した音響波の信号データに基づいて、音響画像を生成する処理部を備え、
前記処理部は、
前記被検部位の撮影領域内で部分的な関心領域の指定を受け付ける処理と、
前記関心領域において、音速パラメータを変更しながら、当該変更された音速パラメータに基づいて、音響画像からなる関心領域画像を順次再構成する処理と、
前記関心領域画像に基づいて前記音速パラメータの適切な値を決定する処理と、
前記音速パラメータの前記適切な値に基づいて、前記関心領域よりも広い領域の音響画像を再構成する処理と、
を実行し、
前記関心領域の指定を受け付ける処理では、
複数の関心領域の指定を受け付け、
前記関心領域画像を再構成する処理では、
前記複数の関心領域において、共通に設定された音速パラメータに基づいて、複数の関心領域画像を再構成し、
前記音速パラメータの前記適切な値を決定する処理では、
前記複数の関心領域において、共通する音速パラメータの適切な値を決定する
画像生成装置が提供される。
所定の被検部位を測定した音響波の信号データに基づいて、音響画像を生成する処理部を備え、
前記処理部は、
前記被検部位の撮影領域内で部分的な関心領域の指定を受け付ける処理と、
前記関心領域において、音速パラメータを変更しながら、当該変更された音速パラメータに基づいて、音響画像からなる関心領域画像を順次再構成する処理と、
前記関心領域画像に基づいて前記音速パラメータの適切な値を決定する処理と、
前記音速パラメータの前記適切な値に基づいて、前記関心領域よりも広い領域の音響画像を再構成する処理と、
を実行し、
前記関心領域の指定を受け付ける処理では、
複数の関心領域の指定を受け付け、
前記関心領域画像を再構成する処理では、
前記複数の関心領域において、共通に設定された音速パラメータに基づいて、複数の関心領域画像を再構成し、
前記音速パラメータの前記適切な値を決定する処理では、
前記複数の関心領域において、共通する音速パラメータの適切な値を決定する
画像生成装置が提供される。
本発明のさらに他の態様によれば、
所定の被検部位を測定した音響波の信号データに基づいて、音響画像を生成する処理部を備え、
前記処理部は、
前記被検部位の撮影領域内で部分的な関心領域の指定を受け付ける処理と、
前記関心領域において、音速パラメータを変更しながら、当該変更された音速パラメータに基づいて、音響画像からなる関心領域画像を順次再構成する処理と、
前記関心領域画像に基づいて前記音速パラメータの適切な値を決定する手順と、
前記音速パラメータの前記適切な値に基づいて、前記関心領域よりも広い領域の音響画像を再構成する処理と、
を実行し、
前記関心領域の指定を受け付ける処理では、
複数の関心領域の指定を受け付け、
前記関心領域画像を再構成する処理では、
前記複数の関心領域において、それぞれ個別に設定された音速パラメータに基づいて、複数の関心領域画像を再構成する
画像生成装置が提供される。
所定の被検部位を測定した音響波の信号データに基づいて、音響画像を生成する処理部を備え、
前記処理部は、
前記被検部位の撮影領域内で部分的な関心領域の指定を受け付ける処理と、
前記関心領域において、音速パラメータを変更しながら、当該変更された音速パラメータに基づいて、音響画像からなる関心領域画像を順次再構成する処理と、
前記関心領域画像に基づいて前記音速パラメータの適切な値を決定する手順と、
前記音速パラメータの前記適切な値に基づいて、前記関心領域よりも広い領域の音響画像を再構成する処理と、
を実行し、
前記関心領域の指定を受け付ける処理では、
複数の関心領域の指定を受け付け、
前記関心領域画像を再構成する処理では、
前記複数の関心領域において、それぞれ個別に設定された音速パラメータに基づいて、複数の関心領域画像を再構成する
画像生成装置が提供される。
本発明によれば、関心領域について画質が特に良好な音響画像を早く再構成することができる。
＜用語＞
本明細書における用語は、例えば、以下のように定義される。
本明細書における用語は、例えば、以下のように定義される。
「被検体（被験体）１００」とは、例えば、検査対象の生体（人体）である。被検体１００となる人のことを「被検者」という。
「被検部位（被験部位）１１０」とは、被検体１００を構成する所定部位であって、音響波を測定（検出）する部位のことを意味する。被検部位１１０は、例えば、手、足、顔、体幹部、乳房などである。
「撮影領域」とは、被検部位１１０内で音響波を測定（撮影、撮像）した実空間としての領域のことを意味する。「撮影領域」は、被検部位１１０の一部であってもよいし、被検部位１１０の全体であってもよい。
「ユーザ」とは、画像生成装置を操作する人のことを意味し、装置の操作のみを行う人に必ずしも限られない。被検者自身が装置を操作する場合は、当該被検者をユーザとして考えてもよい。
「被検部位（被験部位）１１０」とは、被検体１００を構成する所定部位であって、音響波を測定（検出）する部位のことを意味する。被検部位１１０は、例えば、手、足、顔、体幹部、乳房などである。
「撮影領域」とは、被検部位１１０内で音響波を測定（撮影、撮像）した実空間としての領域のことを意味する。「撮影領域」は、被検部位１１０の一部であってもよいし、被検部位１１０の全体であってもよい。
「ユーザ」とは、画像生成装置を操作する人のことを意味し、装置の操作のみを行う人に必ずしも限られない。被検者自身が装置を操作する場合は、当該被検者をユーザとして考えてもよい。
「音響波」とは、所定の媒体を伝播する弾性波（粗密波）のことを意味する。
「光音響効果」とは、以下の現象のことを意味する。すなわち、被検体１００の所定の被検部位１１０に対して光を照射する。被検部位１１０の組織で光が吸収されると、光を吸収した部分が熱を放出し、体積膨張によって音響波を発生させる。このようにして音響波が生じる現象のことを「光音響効果」という。
「光音響波」とは、光音響効果により発生する音響波のことを意味し、「光超音波」と称されることもある。
音響波の「信号」または「検出信号Ｓ（ｉ，ｔ）」とは、検出（受信）した音響波を変換した電気信号（例えばデジタル信号）のことを意味し、「音響信号」ということもある。
「信号データ」とは、音響信号の電子データのことを意味し、「信号強度データ」ともいう。光音響波の信号データを「光音響信号データ」ともいう。
「信号データ量」とは、信号データが含む情報量のことを意味し、「データ容量」ということもある。「信号データ量」は、例えば、画像再構成するときの画素サイズ、同一画素２２の再構成に用いる信号の重畳数、被検部位１１０の所定方向における信号データの範囲、信号を受信したセンサ数（チャンネル数）などに依存する。
「計算量」とは、例えば、画像再構成するためにコンピュータ７０９が計算する処理量のことを意味する。「計算量」は、例えば、信号データ量、処理方法および処理数などに依存する。
「光音響効果」とは、以下の現象のことを意味する。すなわち、被検体１００の所定の被検部位１１０に対して光を照射する。被検部位１１０の組織で光が吸収されると、光を吸収した部分が熱を放出し、体積膨張によって音響波を発生させる。このようにして音響波が生じる現象のことを「光音響効果」という。
「光音響波」とは、光音響効果により発生する音響波のことを意味し、「光超音波」と称されることもある。
音響波の「信号」または「検出信号Ｓ（ｉ，ｔ）」とは、検出（受信）した音響波を変換した電気信号（例えばデジタル信号）のことを意味し、「音響信号」ということもある。
「信号データ」とは、音響信号の電子データのことを意味し、「信号強度データ」ともいう。光音響波の信号データを「光音響信号データ」ともいう。
「信号データ量」とは、信号データが含む情報量のことを意味し、「データ容量」ということもある。「信号データ量」は、例えば、画像再構成するときの画素サイズ、同一画素２２の再構成に用いる信号の重畳数、被検部位１１０の所定方向における信号データの範囲、信号を受信したセンサ数（チャンネル数）などに依存する。
「計算量」とは、例えば、画像再構成するためにコンピュータ７０９が計算する処理量のことを意味する。「計算量」は、例えば、信号データ量、処理方法および処理数などに依存する。
「音響画像」とは、音響波の信号データに基づいて再構成した２次元または３次元の画像のことを意味する。「音響画像」は、光音響波の信号データに基づいて再構成した「光音響画像」、および超音波の信号データに基づいて再構成した「超音波画像」を含んでいる。このような画像は「再構成画像」とも呼ばれる。なお、再構成画像が３次元画像（ボクセルの集合体）の場合には、再構成画像はボリュームデータとも呼ばれる。
また、「音響画像」は、例えば、被検部位１１０内における２次元または３次元の特性情報分布を示している。具体的な「特性情報」としては、例えば、音響波の発生源の位置、被検部位１１０内の初期音圧、初期音圧に基づいて求められるエネルギー吸収密度ならびに吸収係数、被検部位１１０の組織を構成する物質の濃度などが挙げられる。具体的な「特性情報分布」としては、例えば、初期音圧分布、エネルギー吸収密度分布、吸収係数分布、酸素飽和度分布などが挙げられる。
また、「音響画像」は、例えば、被検部位１１０内における２次元または３次元の特性情報分布を示している。具体的な「特性情報」としては、例えば、音響波の発生源の位置、被検部位１１０内の初期音圧、初期音圧に基づいて求められるエネルギー吸収密度ならびに吸収係数、被検部位１１０の組織を構成する物質の濃度などが挙げられる。具体的な「特性情報分布」としては、例えば、初期音圧分布、エネルギー吸収密度分布、吸収係数分布、酸素飽和度分布などが挙げられる。
「仮の測定画像２０」とは、関心領域２４の指定を受け付けるために用いられる簡易的な音響画像のことを意味する。
「関心領域２４」とは、仮の測定画像２０内でユーザが関心を持っている領域のことを意味する。
「関心領域画像４０」とは、関心領域２４において再構成した音響画像のことを意味する。
「広域画像５０」とは、関心領域２４よりも広い領域を有する音響画像のことを意味する。
これらについては、具体的な実施形態において詳細を説明する。
「関心領域２４」とは、仮の測定画像２０内でユーザが関心を持っている領域のことを意味する。
「関心領域画像４０」とは、関心領域２４において再構成した音響画像のことを意味する。
「広域画像５０」とは、関心領域２４よりも広い領域を有する音響画像のことを意味する。
これらについては、具体的な実施形態において詳細を説明する。
「画素２２」（ボクセルまたはピクセル）とは、被検部位１１０内の単位空間領域の特性情報を、音響画像内に所定の画素値（輝度値）で表した単位領域のことを意味する。つまり、音響画像内の画素２２と、被検部位１１０内の単位空間領域とが相互に対応している。
「解像度」とは、再構成した音響画像における画素２２の密度（単位体積当たりの画素数）のことを意味する。例えば、「低解像度」または「高解像度」とは、それぞれ、画素２２の密度が低いこと、または画素２２の密度が高いことを意味する。
「画質」とは、コントラスト、解像度、アーチファクト（虚像）の程度、所定の注目部位の見え方などのことを意味する。
「解像度」とは、再構成した音響画像における画素２２の密度（単位体積当たりの画素数）のことを意味する。例えば、「低解像度」または「高解像度」とは、それぞれ、画素２２の密度が低いこと、または画素２２の密度が高いことを意味する。
「画質」とは、コントラスト、解像度、アーチファクト（虚像）の程度、所定の注目部位の見え方などのことを意味する。
「音速」とは、所定の媒体内を伝播する音響波の伝播速度のことを意味する。
「音速パラメータ」とは、音速に関連したパラメータのことを意味する。当該パラメータが音速に関連していれば、音速パラメータは「音速」に限られるものではない。音速以外の音速パラメータとしては、例えば、単位時間当たりに伝播する音響波の距離である「波長」（音速を所定のサンプリング周波数（固定値）で除したパラメータ）、単位長当たりの波の個数である「波数」（上記パラメータの逆数）などが挙げられる。
「音速パラメータの適切な値」とは、所定の条件下で適切であると判断された音速パラメータの特定値のことを意味する。
「音速パラメータ」とは、音速に関連したパラメータのことを意味する。当該パラメータが音速に関連していれば、音速パラメータは「音速」に限られるものではない。音速以外の音速パラメータとしては、例えば、単位時間当たりに伝播する音響波の距離である「波長」（音速を所定のサンプリング周波数（固定値）で除したパラメータ）、単位長当たりの波の個数である「波数」（上記パラメータの逆数）などが挙げられる。
「音速パラメータの適切な値」とは、所定の条件下で適切であると判断された音速パラメータの特定値のことを意味する。
「ＸＹ方向」とは、センサ３４０側から被検部位１１０を見たときの縦横に直交する水平方向のことを意味し、音響画像では、沿面方向（画面の画素配列方向）に相当する。
「Ｚ方向」とは、センサ３４０側から被検部位１１０を見たときの深さ方向のことを意味し、音響画像では、奥行方向（画面の法線方向）に相当する。
「Ｚ方向」とは、センサ３４０側から被検部位１１０を見たときの深さ方向のことを意味し、音響画像では、奥行方向（画面の法線方向）に相当する。
＜本発明の第１実施形態＞
以下、本発明の第１実施形態について図面を参照しながら説明する。
以下、本発明の第１実施形態について図面を参照しながら説明する。
（１）光音響イメージング装置
図１に示すように、本実施形態に係る光音響イメージング装置１０は、例えば、画像生成装置の一態様として構成され、所定の被検部位１１０に光を照射して測定した光音響波の信号データに基づいて、音響画像を生成するよう構成されている。具体的には、光音響イメージング装置１０は、例えば、支持台（基台）２００と、センサユニット（検出ユニット、受信ユニット）３００と、光源６２０と、光学系６４０と、走査機構（移動機構）３８０と、分離部４００と、処理部７００と、表示部７２０と、入力部７４０と、を備えている。
図１に示すように、本実施形態に係る光音響イメージング装置１０は、例えば、画像生成装置の一態様として構成され、所定の被検部位１１０に光を照射して測定した光音響波の信号データに基づいて、音響画像を生成するよう構成されている。具体的には、光音響イメージング装置１０は、例えば、支持台（基台）２００と、センサユニット（検出ユニット、受信ユニット）３００と、光源６２０と、光学系６４０と、走査機構（移動機構）３８０と、分離部４００と、処理部７００と、表示部７２０と、入力部７４０と、を備えている。
（支持台）
支持台２００は、例えば、被検体１００が載置される基台として構成されている。具体的には、支持台２００は、例えば、支持面２１０と、開口２２０と、を有している。
支持台２００は、例えば、被検体１００が載置される基台として構成されている。具体的には、支持台２００は、例えば、支持面２１０と、開口２２０と、を有している。
支持面２１０は、例えば、被検体１００のうち被検部位１１０以外の部分を支持している。支持台２００の支持面２１０の下には、空間があけられており、後述のセンサユニット３００などが設けられている。
開口２２０は、例えば、被検体１００の所定の被検部位１１０を測定するために、支持面２１０に開設されている。開口２２０は、所定の被検部位１１０からの音響波を測定するため、被検部位１１０よりも広く設けられている。開口２２０の平面形状は、例えば、四角形である。
（分離部）
分離部４００は、例えば、被検体１００側とセンサユニット３００側とを分離するよう構成されている。本実施形態では、分離部４００は、例えば、分離フィルム４２０を有している。
分離部４００は、例えば、被検体１００側とセンサユニット３００側とを分離するよう構成されている。本実施形態では、分離部４００は、例えば、分離フィルム４２０を有している。
分離フィルム４２０は、例えば、音響整合材３１０を浸透させないよう構成されている。また、分離フィルム４２０は、例えば、光源６２０からの光を透過するよう構成されている。さらに、分離フィルム４２０は、例えば、被検部位１１０からの音響波を伝播させることができるように、被検体１００と整合する音響インピーダンスを有している。具体的には、上述の要件を満たす分離フィルム４２０の材料としては、例えば、ポリエチレンテレフタラート（ＰＥＴ）、ポリエチレン（ＰＥ）などが挙げられる。
分離フィルム４２０の厚さは、被検部位１１０からの音響波の周波数帯域と分離フィルム４２０内の縦波音速とに基づいて決定される。
分離フィルム４２０は、例えば、支持台２００の開口２２０を塞ぐ（覆う）ように設けられ、支持台２００に固定されている。分離フィルム４２０は、例えば、凹部（符号不図示）を有し、音響整合材４１０を凹部内に収容可能に構成されている。
音響整合材４１０は、例えば、液状またはゲル状であり、被検体１００と整合する音響インピーダンスを有している。本実施形態でいう音響インピーダンスが「被検体１００と整合する」とは、被検体１００の音響インピーダンスと完全に一致する場合だけでなく、被検体１００の音響インピーダンスに所定の誤差で近似される場合も含んでいる。具体的には、「被検体１００と整合する音響インピーダンス」は、例えば、被検体１００の音響インピーダンスの０．５倍以上２倍以下の範囲内である。具体的な音響整合材３１０としては、例えば、水、油などである。
分離フィルム４２０は、例えば、後述の容器３２０内に収容された音響整合材３１０に接触している。
（センサユニット）
センサユニット３００は、例えば、被検体１００の所定の被検部位１１０からの音響波を受信するよう構成されている。本実施形態のセンサユニット３００は、例えば、容器３２０と、センサ（探触子、変換素子）３４０と、素子保持部３６０と、を有している。
センサユニット３００は、例えば、被検体１００の所定の被検部位１１０からの音響波を受信するよう構成されている。本実施形態のセンサユニット３００は、例えば、容器３２０と、センサ（探触子、変換素子）３４０と、素子保持部３６０と、を有している。
［容器］
容器３２０は、例えば、支持面２１０よりも鉛直下側に設けられている。容器３２０は、例えば、音響整合材３１０を収容（貯留）可能に構成されている。
容器３２０は、例えば、支持面２１０よりも鉛直下側に設けられている。容器３２０は、例えば、音響整合材３１０を収容（貯留）可能に構成されている。
音響整合材３１０は、例えば、音響整合材４１０と同様に、液状であり、被検体１００と整合する音響インピーダンスを有している。具体的な音響整合材３１０としては、例えば、水、油などである。
本実施形態では、容器３２０は、例えば、音響整合材３１０を固定せず不定形に変化させることが可能な状態で音響整合材３１０を収容しており、すなわち、流動性を有する状態で音響整合材３１０を収容している。
また、本実施形態では、容器３２０内において、分離フィルム４２０に接する位置まで、音響整合材３１０が充填される。これにより、被検部位１１０からセンサ３４０までの音響波の伝播経路に空気が介在することを抑制することができる。
［センサ］
センサ３４０は、例えば、支持面２１０よりも鉛直下側に設けられている。センサ３４０は、例えば、被検部位１１０から発生する音響波を受信するよう構成されている。
センサ３４０は、例えば、支持面２１０よりも鉛直下側に設けられている。センサ３４０は、例えば、被検部位１１０から発生する音響波を受信するよう構成されている。
また、センサ３４０は、例えば、受信した音響波を電気信号に変換するよう構成されている。センサ３４０は、例えば、１００ｋＨｚ以上１０００ＭＨｚ以下の周波数を有する音響波を受信可能に構成されている。より好ましくは、センサ３４０は、例えば、１００ｋＨｚ以上５０ＭＨｚ以下の周波数を有する音響波を受信可能に構成されている。具体的なセンサ３４０としては、例えば、チタン酸ジルコン酸鉛（ＰＺＴ）などからなる圧電素子、ポリフッ化ビニリデン（ＰＶＤＦ）などの高分子圧電膜材料、容量性マイクロマシン超音波トランスデューサ（ＣＭＵＴ）、ファブリペロー干渉計などが挙げられる。
本実施形態では、センサ３４０は、例えば、複数設けられている。複数のセンサ３４０により音響波を受信することで、測定精度を向上させることができる。例えば、被検部位１１０内の特性情報の測定位置精度を向上させることができる。
［素子保持部］
素子保持部３６０は、例えば、複数のセンサ３４０を保持している。素子保持部３６０は、例えば、鉛直下側に向けて凹んだ半球状（お椀状）に構成されている。ここでいう「半球状」とは、平坦な断面で分割された真円球の形状、平坦な断面で分割された楕円球の形状、またはそれらに所定の誤差で近似される形状のことを意味している。素子保持部３６０が構成する球面の中心角は、例えば、１４０°以上１８０°以下である。
素子保持部３６０は、例えば、複数のセンサ３４０を保持している。素子保持部３６０は、例えば、鉛直下側に向けて凹んだ半球状（お椀状）に構成されている。ここでいう「半球状」とは、平坦な断面で分割された真円球の形状、平坦な断面で分割された楕円球の形状、またはそれらに所定の誤差で近似される形状のことを意味している。素子保持部３６０が構成する球面の中心角は、例えば、１４０°以上１８０°以下である。
素子保持部３６０は、例えば、複数のセンサ３４０のそれぞれの指向軸が半球面の曲率中心付近に集中するように、半球面に沿ってアレイ状に複数のセンサ３４０を保持している。これにより、半球面の曲率中心付近において、高分解能を得ることができる。
本実施形態では、例えば、後述の音響整合材４１０の中に被検部位１１０を浸漬させたとき、または音響整合材４１０の中に被検部位１１０を載置させたとき、素子保持部３６０の半球面の曲率中心が被検部位１１０内に位置するように設定されている。これにより、所定の被検部位１１０内で高分解能の測定を行うことができる。
本実施形態では、素子保持部３６０は、例えば、容器３２０の底部に設けられ、容器３２０に一体として固定されている。素子保持部３６０内には、上述の音響整合材３１０が収容される。これにより、センサ３４０は、音響整合材３１０を介して音響波を受信するようになっている。
また、本実施形態では、上述のように容器３２０が、流動性を有する状態で音響整合材３１０を収容していることで、素子保持部３６０が複雑な形状を有していても、素子保持部３６０内に、空気を介在させることなく音響整合材３１０を密に充填することができる。
（光源）
光源６２０は、例えば、所定の被検部位１１０に対して光を照射するよう構成されている。光源６２０は、例えば、パルス光を出射可能に構成されている。具体的には、光源６２０は、例えば、レーザ、発光ダイオード、またはフラッシュランプである。レーザとしては、例えば、ガスレーザ、固体レーザ、色素レーザ、半導体レーザなどが挙げられる。
光源６２０は、例えば、所定の被検部位１１０に対して光を照射するよう構成されている。光源６２０は、例えば、パルス光を出射可能に構成されている。具体的には、光源６２０は、例えば、レーザ、発光ダイオード、またはフラッシュランプである。レーザとしては、例えば、ガスレーザ、固体レーザ、色素レーザ、半導体レーザなどが挙げられる。
光源６２０は、例えば、光音響効果が得られる条件下で、光を出射するよう構成されている。
光源６２０から出射される光の波長は、例えば、被検部位１１０の組織を構成する所定の吸収体に吸収される波長であり、かつ、被検部位１１０の内部まで伝播可能な波長である。具体的には、光の波長は、例えば、５００ｎｍ以上１２００ｎｍ以下である。
なお、光源６２０は、例えば、異なる波長の光を出射可能に構成されていてもよい。異なる波長の光を被検部位１１０に照射することで、異なる波長における吸収係数の違いに基づいて特性情報の分布を得ることができる。例えば、酸素飽和度分布などを得ることができる。
光源６２０から出射される光のパルス幅は、いわゆる熱閉じ込め条件およびストレス閉じ込め条件の両方を満たしている。すなわち、パルス幅は、被検部位１１０内の所定の吸収体から熱が伝播して逃げる前に光照射が終了する時間幅であり、かつ、音響波が吸収体内を通過する前に光照射が終了する時間幅である。具体的には、パルス幅は、例えば、１ｎｓ以上１００ｎｓ以下である。
（光学系）
光学系６４０は、例えば、光源６２０からの光を伝送するよう構成されている。光学系６４０は、例えば、レンズならびにミラーなどの光学部品、光ファイバなどにより構成されている。
光学系６４０は、例えば、光源６２０からの光を伝送するよう構成されている。光学系６４０は、例えば、レンズならびにミラーなどの光学部品、光ファイバなどにより構成されている。
光学系６４０の終端における光出射口６６０は、光源６２０から伝送された光を被検部位１１０に向けて出射するよう構成されている。光出射口６６０は、例えば、素子保持部３６０の底部に設けられている。光出射口６６０がセンサ３４０とともに素子保持部３６０に設けられていることで、被検部位１１０内の広い範囲における光音響波を測定することができる。
光出射口６６０は、例えば、出射される光量が、生体としての被検体１００に照射可能な光エネルギーの最大許容露光量（Ｍａｘｉｍｕｍ Ｐｅｒｍｉｓｓｉｖｅ Ｅｘｐｏｓｕｒｅ、ＭＰＥ）を超えないように、光学的に設計されていることが好ましい。具体的には、光出射口６６０は、例えば、凹レンズなどを有することで、照射範囲が広げられていることが好ましい。
（走査機構）
走査機構３８０は、例えば、支持台２００上に載置された被検体１００に対して、センサ３４０を相対的に走査（移動）させるよう構成されている。本実施形態では、走査機構３８０は、例えば、容器３２０とセンサ３４０とを有するセンサユニット３００を一体として走査させるよう構成されている。
走査機構３８０は、例えば、支持台２００上に載置された被検体１００に対して、センサ３４０を相対的に走査（移動）させるよう構成されている。本実施形態では、走査機構３８０は、例えば、容器３２０とセンサ３４０とを有するセンサユニット３００を一体として走査させるよう構成されている。
走査機構３８０は、少なくとも所定の１方向にセンサ３４０を走査させるよう構成されている。走査機構３８０がセンサ３４０を走査させる方向は、例えば、２次元方向（ＸＹ方向）であっても、または３次元方向（ＸＹＺ方向）であってもよい。走査の方向は一方向の直線的な動きに限らず、回転させる動きを採用してもよい。本実施形態では、走査機構３８０は、例えば、支持面２１０に平行な水平面上をＸＹ方向にセンサ３４０を走査させるよう構成されている。
本実施形態では、上述のように容器３２０が、流動性を有する状態で音響整合材３１０を収容していることで、走査機構３８０によりセンサユニット３００を走査させても、センサ３４０が音響整合材３１０と接した状態を維持することができる。
（供給部）
供給部（図示せず）は、例えば、容器３２０内に供給管を介して音響整合材３１０を供給するよう構成されている。供給管は、例えば、素子保持部３６０の一部に接続されている。供給部からの音響整合材３１０の供給により、容器３２０内の音響整合材３１０の上面が所定位置に維持される。
供給部（図示せず）は、例えば、容器３２０内に供給管を介して音響整合材３１０を供給するよう構成されている。供給管は、例えば、素子保持部３６０の一部に接続されている。供給部からの音響整合材３１０の供給により、容器３２０内の音響整合材３１０の上面が所定位置に維持される。
（処理部）
処理部７００は、例えば、光音響イメージング装置１０の各部を制御し、取得した信号データに基づいて画像を再構成するよう構成され、すなわち、被検部位１１０内の特性情報を処理するよう構成されている。
処理部７００は、例えば、光音響イメージング装置１０の各部を制御し、取得した信号データに基づいて画像を再構成するよう構成され、すなわち、被検部位１１０内の特性情報を処理するよう構成されている。
図２に示すように、処理部７００は、例えば、コンピュータ７０９および信号処理部７１０を有している。具体的には、コンピュータ７０９は、例えば、ＣＰＵ７０１（Ｃｅｎｔｒａｌ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）、ＧＰＵ７０２（Ｇｒａｐｈｉｃｓ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）、ＲＡＭ７０３（Ｒａｎｄｏｍ Ａｃｃｅｓｓ Ｍｅｍｏｒｙ）、記憶装置７０４、およびＩ／Ｏポート７０５を有している。ＲＡＭ７０３、記憶装置７０４、およびＩ／Ｏポート７０５は、ＣＰＵ７０１とデータ交換可能に構成されている。Ｉ／Ｏポート７０５は、例えば、所定の増幅器、ＡＤ変換器および演算回路などの信号処理部７１０を介してセンサ３４０のそれぞれに接続され、さらに、光源６２０、走査機構３８０、表示部７２０および入力部７４０に接続されている。
記憶装置７０４は、例えば、光音響波測定に係るプログラム、画像再構成に係るプログラム（画像生成プログラム）、信号データ、および被検部位１１０内の特性情報などを記憶するよう構成されている。記憶装置７０４は、例えば、ＨＤＤ（Ｈａｒｄ Ｄｉｓｋ Ｄｒｉｖｅ）、フラッシュメモリなどである。ＲＡＭ７０３は、ＣＰＵ７０１またはＧＰＵ７０２によって記憶装置７０４から読み出される情報やプログラム等が一時的に保持されるよう構成されている。
ＣＰＵ７０１は、記憶装置７０４に格納された所定のプログラムを実行することで、光音響イメージング装置１０の各部を制御し、光音響波測定に係る処理および信号処理などを実行するように構成されている。ＧＰＵ７０２は、記憶装置７０４に格納された所定のプログラムを実行することで、画像再構成に係る処理などを実行するように構成されている。なお、ＧＰＵ７０２は、画像再構成に係る処理を単独で実行してもよいし、ＣＰＵ７０１と連携して実行してもよい。
表示部７２０は、例えば、所定のプログラムの実行によって再構成された音響画像などを表示するよう構成されている。表示部７２０は、例えば、液晶ディスプレイ、有機ＥＬ（ＯＬＥＤ）ディスプレイ、ヘッドマウントディスプレイ、直視型立体ディスプレイなどである。
入力部７４０は、例えば、ユーザが所定の操作を行う情報をコンピュータ７０９に入力可能に構成されている。入力部７４０は、例えば、マウス、キーボード、トラックボール、ジョイスティック、タッチパネル、マイク、視線入力デバイスなどである。
なお、本実施形態における画像生成プログラムは、音響画像の再構成を高速に行うために高性能なコンピュータ７０９を用いることが好ましい。特に近年では、ＧＰＵの発展が目覚ましく、ＧＰＵを汎用演算に用いるＧＰＧＰＵ（Ｇｅｎｅｒａｌ－Ｐｕｒｐｏｓｅ ｃｏｍｐｕｔｉｎｇ ｏｎ ＧＰＵ）が広く用いられる。本実施形態においてＧＰＵ７０２としてＧＰＧＰＵを用いることによって、データ量の多いボリュームデータを再構成する場合であっても、１ショット当たり数十ミリ秒以内で計算することが可能となる。
（２）画像再構成の原理
次いで、図３を用い、画像再構成の原理について説明する。図３において、音響画像を再構成したときの画素２２を被検部位１１０と重ね合わせて示している。また、図３において、所定の音響波を発生させる音源１２０の１つを黒点として示している。各センサ３４０は、被検体１００内部から伝播してくる音響波を検出して、検出信号を出力する。なお、センサ３４０の数はＮ（Ｎは２以上の整数）個とする。
次いで、図３を用い、画像再構成の原理について説明する。図３において、音響画像を再構成したときの画素２２を被検部位１１０と重ね合わせて示している。また、図３において、所定の音響波を発生させる音源１２０の１つを黒点として示している。各センサ３４０は、被検体１００内部から伝播してくる音響波を検出して、検出信号を出力する。なお、センサ３４０の数はＮ（Ｎは２以上の整数）個とする。
音響波（超音波）を利用した画像生成装置では、複数のセンサ３４０ｉから得られる複数の検出信号Ｓ（ｉ，ｔ）に基づいて被検体１００内部の特性情報を表す画像を再構成する。再構成画像の各画素２２（ボクセル）の輝度値（または画素値）は、各センサ３４０からその画素２２に対応する位置（単位空間領域）までの距離と被検体１００内の音速とにより位相が調整された検出信号に基づいて算出される。なお、検出信号Ｓ（ｉ，ｔ）において、ｉはセンサ３４０の番号（０からＮ－１の範囲の整数）、ｔは時間を示す。
以下、画像再構成法の一例であるタイムドメイン法について具体的に説明する。まず、センサ３４０毎に、ｉ番目のセンサ３４０から再構成画像内の或る画素２２に対応する位置までの距離を伝播速度（音速）で除算する。それにより、その画素２２に対応する位置で音響波が発生した場合において、該音響波がセンサ３４０ｉに到達するまでの時間τｉ（遅延時間）を算出する（発生時刻をｔ＝０とする）。そして、センサ３４０毎に、時間τｉのときの検出信号（の強度）Ｓ（ｉ，τｉ）を算出し、それらを足し合わせることにより、再構成画像の画素値を生成する。同様に複数の位置について画素２２を生成することにより、再構成画像が生成される。本実施形態では、例えば、光音響イメージング装置１０のタイムドメインでの画像再構成法の一つであり、以下の式（１）で表されるＵｎｉｖｅｒｓａｌ ｂａｃｋ－ｐｒｏｊｅｃｔｉｏｎ（ＵＢＰ）法を使用できる。
ここで、ｒ０は再構成する位置（再構成位置、注目位置とも呼ぶ）を示す位置ベクトル、ｐ０（ｒ０，ｔ）は再構成する位置の初期音圧、ｃは伝搬経路の音速を示す。また、ΔΩｉは再構成する位置からｉ番目のセンサ３４０を見込む立体角、Ｎは再構成に用いるセンサ３４０の個数を示す。式（１）は、受信信号ｐ（ｒｉ，ｔ）に微分等の処理を行い、それらに立体角の加重をかけて整相加算すること（逆投影）を示している。式（１）のｔは、注目位置とセンサ３４０とを結ぶ音線を光音響波が伝搬する時間（伝搬時間）である。なお、ｂ（ｒｉ、ｔ）の計算においては、他にも演算処理を施してもよい。他の演算処理としては、例えば、周波数フィルタリング（ローパス、ハイパス、バンドパス等）、デコンボリューション、包絡線検波、ウェーブレットフィルタリング、等が挙げられる。また、本発明においては、センサ３４０と注目位置とを結ぶ音線の伝搬時間を求めて再構成する方法であれば、どのような再構成アルゴリズムを用いても良い。例えば、タイムドメインでの逆投影法として、Ｆｉｌｔｅｒｅｄ ｂａｃｋ－ｐｒｏｊｅｃｔｉｏｎなどを採用してもよい。タイムドメインではなく、フーリエドメインで再構成してもよい。
（３）画像生成方法
次いで、本実施形態に係る画像生成方法について説明する。
次いで、本実施形態に係る画像生成方法について説明する。
［光音響イメージング装置１０による基本的な画像生成方法］
光音響イメージング装置１０によって得られた光音響信号データを処理部７００内にある記憶装置７０４に保存する。光音響信号データは、サンプリング時間とサンプリング周波数の積によってデータ数が決まる１チャンネル（１センサ）当たりの信号強度データ、及びこれとセンサ数との積のデータ容量を有する１ショットデータ、及びこれとトータルショット数との積のデータ容量を有する全体のスキャンデータからなる。すなわち、光音響イメージング装置１０では、光源６２０としてのレーザを１回照射すると、被検体１００から光音響波が発生し、当該光音響波を複数チャンネルのセンサ３４０で受信することによって、１ショットのレーザ照射によって得られるＮ個分の信号強度データを得ることができる。１ショットデータを取得したのち、走査機構３８０によりセンサユニット３００の位置を変化させて、別の位置での１ショットデータを取得する。これらを繰り返すことによって、撮影領域全体の信号データ（スキャンデータ）を完成させる。このようにして得られた信号データは、例えば、センサユニット３００の位置を示す座標情報と関連づけられている。
光音響イメージング装置１０によって得られた光音響信号データを処理部７００内にある記憶装置７０４に保存する。光音響信号データは、サンプリング時間とサンプリング周波数の積によってデータ数が決まる１チャンネル（１センサ）当たりの信号強度データ、及びこれとセンサ数との積のデータ容量を有する１ショットデータ、及びこれとトータルショット数との積のデータ容量を有する全体のスキャンデータからなる。すなわち、光音響イメージング装置１０では、光源６２０としてのレーザを１回照射すると、被検体１００から光音響波が発生し、当該光音響波を複数チャンネルのセンサ３４０で受信することによって、１ショットのレーザ照射によって得られるＮ個分の信号強度データを得ることができる。１ショットデータを取得したのち、走査機構３８０によりセンサユニット３００の位置を変化させて、別の位置での１ショットデータを取得する。これらを繰り返すことによって、撮影領域全体の信号データ（スキャンデータ）を完成させる。このようにして得られた信号データは、例えば、センサユニット３００の位置を示す座標情報と関連づけられている。
１ショットデータに基づいて画像再構成することによって、所定の範囲（Ｆｉｅｌｄ ｏｆ ｖｉｅｗ，ＦＯＶ）のボリュームデータを生成することができる。センサユニット３００の位置を変化させたのちに、別の１ショットデータに基づいて画像再構成することによって、別の場所のボリュームデータを生成することができる。これら二つのボリュームデータが部分的にでもオーバーラップしていれば、センサユニット３００の位置情報を参照しながら、ボリュームデータ同士を加算させて音響画像を構成する。この加算効果によってアーチファクトを減少させ、高画質な音響画像を得ることができる。このように、センサユニット３００の位置情報を参照しながら、すべての信号データを重ね合わせることによって、広い範囲で高画質なボリュームデータを作成することができる。
［本実施形態の概略］
本実施形態では、音速パラメータの適切な値を決定するにあたり、まず関心領域２４を指定し、その関心領域２４の画像を見ながらユーザが即時的に音速パラメータの適切な値を探索する方法を採用する。すなわち、音速パラメータを変更するのに応じて画像を再構成し、再構成画像を表示し、更なる変更に応じて同様に再構成および表示を行うことにより、即時的に音速パラメータの適切な値を探索することができる。
本実施形態では、音速パラメータの適切な値を決定するにあたり、まず関心領域２４を指定し、その関心領域２４の画像を見ながらユーザが即時的に音速パラメータの適切な値を探索する方法を採用する。すなわち、音速パラメータを変更するのに応じて画像を再構成し、再構成画像を表示し、更なる変更に応じて同様に再構成および表示を行うことにより、即時的に音速パラメータの適切な値を探索することができる。
［本実施形態の具体的な画像生成方法］
以下、図４～図８を用い、本実施形態の具体的な画像生成方法について説明する。
以下、図４～図８を用い、本実施形態の具体的な画像生成方法について説明する。
図４に示すように、本実施形態の画像生成方法は、例えば、信号データ準備工程Ｓ１２０～Ｓ１４０、仮の測定画像再構成工程Ｓ２２０～Ｓ２４０、関心領域指定受付工程Ｓ３２０、関心領域画像再構成工程Ｓ４２０～Ｓ４６０、適切値決定工程Ｓ４８０、および広域画像再構成工程Ｓ５４０を有している。各工程における処理は、処理部７００により実行される。
（Ｓ１２０～１４０：信号データ準備工程）
まず、所定の被検部位１１０を測定した音響波の信号データを準備する。
まず、所定の被検部位１１０を測定した音響波の信号データを準備する。
本実施形態では、例えば、信号データとして、同一の画素２２に再構成される位置に複数の信号が重畳されたデータを準備する。また、例えば、信号データとして、複数のセンサ３４０により受信した複数の信号を含むデータを準備する。
具体的には、例えば、上述の光音響イメージング装置１０を用い、被検体１００の所定の被検部位１１０に対して光を照射し、被検部位１１０から発生する音響波をセンサ３４０により受信することで、光音響信号を取得する（Ｓ１２０）。光音響信号を取得したら、信号データを記憶装置７０４に保存する（Ｓ１４０）。
（Ｓ２２０～Ｓ２４０：仮の測定画像再構成工程）
次に、図５に示すように、上述の信号データに基づいて、音響画像からなる仮の測定画像２０を再構成する。
次に、図５に示すように、上述の信号データに基づいて、音響画像からなる仮の測定画像２０を再構成する。
具体的には、記憶装置７０４から信号データを読み出し（Ｓ２２０）、信号データに基づいて仮の測定画像２０を再構成し、表示部７２０の所定の画面上に表示させる（Ｓ２４０）。
このとき、本実施形態の仮の測定画像２０は、後述の関心領域２４の指定を受け付けるために用いられる。仮の測定画像２０において、関心領域２４の位置（例えば、所定の血管などの注目位置）を把握できればよいため、画質を低くして（粗くして）仮の測定画像２０を再構成する。
すなわち、本実施形態では、例えば、仮の測定画像２０を再構成するための単位面積当たりの信号データ量を、関心領域画像４０を再構成するための単位面積当たりの信号データ量よりも少なくする。ここでいう「単位面積」とは、音響画像を平面視したとき（すなわち実空間でＸＹ面を見たときに相当）の単位面積のことを意味する。このように仮の測定画像２０を再構成するときの信号データ量を少なくすることで、単位面積当たりの計算量を少なくすることができる。その結果、仮の測定画像２０を早く再構成することができる。
具体的には、例えば、仮の測定画像２０の同一画素２２の再構成に用いる信号データの重畳数を少なくする。例えば、センサユニット３００の位置を移動させながら所定の測定間隔で取得した全信号データから、測定間隔よりも広い間隔で信号データを減少させる。なお、このとき、画素サイズは変化させない。上述のような方法により、同一画素２２の再構成に用いる信号データの重畳数を少なくすることができる。これにより、仮の測定画像２０を再構成するときの単位面積当たりの信号データ量を少なくすることができる。その結果、仮の測定画像２０を再構成するときの単位面積当たりの計算量を少なくすることができる。
このとき、仮の測定画像２０として再構成される領域を、例えば、関心領域２４の指定を受け付けることが可能な程度に広い領域とする。
具体的には、例えば、図５に示すように、撮影領域の全体に亘って、仮の測定画像２０を再構成する。ここでいう「撮影領域の全体」とは、再構成可能な信号データを得た全ての単位空間領域のうちで、最も外側の単位空間領域によって囲まれた範囲内のことを意味する。また、「撮影領域の全体に亘って画像を再構成する」とは、最も外側の単位空間領域によって囲まれた範囲内で画像を再構成することを意味し、全ての単位空間領域から得られた全信号データに基づいて画像を再構成することに限定されるものではない。
なお、仮の測定画像２０が関心領域２４の指定を受け付けることが可能な程度に広ければ、必ずしも撮影領域の全体に亘って仮の測定画像２０を再構成しなくてもよく、撮影領域の一部に基づいて仮の測定画像２０を再構成してもよい。
また、このとき、仮の測定画像２０を再構成するときの音速パラメータは、厳密性（正確性）を欠いても構わない。例えば、仮の測定画像２０を再構成するときの音速パラメータとして、既知の水の音速パラメータを用いてもよい。
また、このとき、後述の関心領域２４を指定しやすくするために、仮の測定画像２０として、例えば、最大強度投影法（Ｍａｘｉｍｕｍ Ｉｎｔｅｎｓｉｔｙ Ｐｒｏｊｅｃｔｉｏｎ，ＭＩＰ）による２次元画像を生成し、表示部７２０の画面上に表示させてもよい。
（Ｓ３２０：関心領域指定受付工程）
次に、図６に示すように、仮の測定画像２０内で、部分的な関心領域２４の指定を受け付ける。
次に、図６に示すように、仮の測定画像２０内で、部分的な関心領域２４の指定を受け付ける。
本実施形態の画像生成プログラムは、上述の方法で得た仮の測定画像２０をユーザが見ながら、関心領域２４の指定をユーザから受け付けるビューワ機能およびグラフィカルユーザーインターフェース（ＧＵＩ）機能を有する。関心領域２４指定の手段は、入力部７４０を用いる。入力部７４０による関心領域２４の具体的な指定方法としては、例えば、マウス、トラックボール、あるいはジョイスティックで指定する方法、タッチパネルに指でタッチして指定する方法、キーボードで入力して指定する方法、マイクにより音声入力して指定する方法、視線入力デバイスを利用して指定する方法などが挙げられる。
このとき、仮の測定画像２０内で、例えば、連続する複数の画素２２を含む関心領域２４の指定を受け付ける。関心領域２４が複数の画素２２を含むことで、所定の血管などの注目部位を確実に関心領域２４内に含めることができる。
また、このとき、仮の測定画像２０内で、例えば、仮の測定画像２０のサイズよりも小さい任意のサイズ（大きさ）を有する関心領域２４の指定を受け付ける。これにより、所定の血管などの注目部位を確実に関心領域２４内に含めること、および、後述の関心領域画像４０を早く再構成することなどの、ユーザが優先する条件に応じて、関心領域２４を指定することができる。
なお、このとき、関心領域２４は、固定されたサイズで位置だけを指定してもよい。
また、このとき、仮の測定画像２０内で、例えば、任意の形状を有する関心領域２４の指定を受け付けてもよい。関心領域２４の形状としては、例えば、円形、楕円形、多角形などが挙げられる。
また、このとき、仮の測定画像２０内で、例えば、複数の関心領域２４の指定を受け付ける。これにより、ユーザが注目する部位が離れた位置にあっても、画質が良好となるように音速パラメータを調整することができる。また、複数の関心領域２４の一部を重ね合わせることで、結合された複数の関心領域２４全体として、上述の例として挙げた形状に限られない任意の形状に関心領域２４を調整することができる。
（Ｓ４２０～Ｓ４６０：関心領域画像再構成工程）
次に、図７に示すように、関心領域２４において、音速パラメータを変更しながら、当該変更された音速パラメータに基づいて、音響画像からなる関心領域画像４０を順次再構成する。
次に、図７に示すように、関心領域２４において、音速パラメータを変更しながら、当該変更された音速パラメータに基づいて、音響画像からなる関心領域画像４０を順次再構成する。
具体的には、所定の音速パラメータに基づいて関心領域画像４０を再構成し、表示部７２０の所定の画面上に表示させる（Ｓ４２０）。関心領域画像４０以外の領域は再構成せず、表示部７２０には、例えば、黒色を表示させる。なお、最初のＳ４２０においては、後述するユーザによる音速パラメータの指定を受け付けてもよいし、予め設定された音速パラメータを用いてもよい。次に、再構成された関心領域画像４０の画質をユーザにより判断する（Ｓ４４０）。関心領域画像４０の画質がＮＧである場合には（Ｓ４４０でＮＧ）、音速パラメータの任意の変更をユーザから受け付ける（Ｓ４６０）。なお、関心領域画像４０の画質がＯＫである場合には（Ｓ４４０でＯＫ）、後述の適切値決定工程Ｓ４８０に進む。
上述のように、本実施形態では、音速パラメータの任意の変更を受け付けながら、当該変更された音速パラメータに基づいて、関心領域画像４０を順次再構成する。これにより、ユーザが関心領域画像４０を視認しながら、音速パラメータを調整することができる。
このとき、本実施形態では、関心領域画像４０を、例えば、仮の測定画像２０の画質よりも高画質に再構成する。すなわち、関心領域画像４０を再構成するための単位面積当たりの信号データ量を、仮の測定画像２０を再構成するための単位面積当たりの信号データ量よりも多くする。例えば、センサユニット３００の位置を移動させながら所定の測定間隔で取得した関心領域２４内のそのままの信号データに基づいて、すなわち、関心領域２４内の信号データを減少させることなく、関心領域画像４０を再構成する。したがって、関心領域画像４０を再構成するための単位面積当たりの計算量を、仮の測定画像２０を再構成するための単位面積当たりの計算量よりも多くすることとなる。
このとき、ユーザによって指定された関心領域画像４０は、センサ３４０の位置が異なる複数ショットでのオーバーラップ領域がある場合、ボリュームデータの重ね合わせで生成される。例えば、異なるセンサ３４０の位置で得られたボリュームデータを、指定された関心領域２４に限定して順次再構成して、各画像を重畳して表示する。このようにして得られた関心領域画像４０としてのボリュームデータは、例えば、ＭＩＰによる２次元画像に変換して、表示部７２０の画面上に表示される。あるいは、例えば、ヘッドマウントディスプレイまたは直視型立体ディスプレイを用いて、３次元ボリュームデータとして関心領域画像４０を表示してもよい。当該ボリュームデータは、Ｓ４２０～Ｓ４６０を含むサイクルを順次繰り返しながら再構成され、関心領域画像４０が音速パラメータの変更に応じて更新される。なお、音速パラメータなどの画像形成のためのパラメータに変更が無い場合には、関心領域画像４０の再構成処理を停止してもよい。
また、このとき用いられる画像生成プログラムは、例えば、図７で示すような音速パラメータを変更するための音速変更ＧＵＩ３０を有する。図７における音速変更ＧＵＩ３０では、例えば、数値ボタンの上に、矢印カーソルを乗せると、カーソルの形状が矢印から楕円に変化する。以下、当該変化後のカーソルを「楕円カーソル３１」と呼ぶ。通常の矢印カーソルを数値ボタンの上に乗せると、数値が見づらくなる。これに対し、矢印カーソルを楕円カーソル３１に変化させることで、数値の視認性を向上させることができる。
図７において、楕円カーソル３１は数値ボタンの下に表示されているが、楕円カーソル３１の指示位置（ｈｏｔｓｐｏｔ）は数値ボタンの中央となる。楕円カーソル３１は、マウスのホイールを模したものである。数値ボタン上で、ホイールを回転させることにより、音速パラメータを増減させることができる。
なお、音速パラメータは、例えば、マウスクリックにより変更可能であってもよい。例えば、マウスの左クリックで音速パラメータを－１減少させることができ、右クリックで音速パラメータを＋１増加させることができる。
さらに、音速パラメータは、例えば、マウスのドラッグにより変更可能であってもよい。例えば、楕円カーソル３１の左側のマイナス（－）は、音速パラメータの減少に対応するボタンおよびドラッグ方向を示している。一方で、楕円カーソル３１の右側のプラス（＋）は、音速パラメータの増加に対応するボタンおよびドラッグ方向を示している。音速パラメータの数値を大きく変えたいときには、マウスの左ボタンで楕円カーソル３１を左方向または右方向にドラッグすることで、音速パラメータを大幅に増減させることができる。
なお、音速変更ＧＵＩ３０は、図７の例に限られない。その他、音速変更ＧＵＩ３０による方法としては、例えば、テキストボックス内に音速パラメータの数値を直接入力する方法、上述の例以外のプラス（＋）キーまたはマイナス（－）キーをクリックすることで音速パラメータを増減させる方法、上述の例以外のスライダバーをクリックしたままドラッグすることで音速パラメータを変更する方法、上述の例以外の態様でマウスホイールを回転させて音速パラメータを変更する方法などが挙げられる。なお、ここで挙げた方法に限られるものではなく、音速パラメータを変化させることができるのであれば、いずれの方法を用いてもよい。
上述のように音速パラメータの変更を受け付けると、変更された音速パラメータに基づいて、即座に画像再構成処理が行われ、ＭＩＰ表示処理を施したのちに画像が表示される。ＧＰＧＰＵによる高速画像再構成により、ユーザは、音速パラメータを変えることによる画像の先鋭度（フォーカス）をリアルタイムで調整することが可能となる。
なお、このとき、異なるセンサ３４０の位置のボリュームデータを順次再構成しているため、音速パラメータを変えた直後の関心領域画像４０は、異なる音速パラメータに基づいて再構成された画像同士の加算画像となる。１ショット当たり数百ミリ秒以内で画像が再構成されている場合には、視認性にはほぼ問題がない。一方で、関心領域２４が大きいと、関心領域画像４０全体の更新に時間がかかる。このため、コンピュータ７０９のスペックが低い場合には、関心領域２４の境界をマウスクリックで変更可能状態とし、関心領域２４を任意のサイズまで狭くすることで、画像再構成におけるリアルタイム性および操作性を確保することが好ましい。関心領域２４のサイズは、コンピュータ７０９のスペックによるが、１秒以内で関心領域画像４０全体が更新される程度の大きさに設定することが好ましい。
上述のように、音速パラメータを変化させることで、被検体１００の画像再構成結果のフォーカス状態をリアルタイムに変えることができる。これにより、ユーザが慣れ親しんだカメラのフォーカス合わせのような使用感で所望のフォーカスとなるように設定することが可能となる。
さらに、このとき、本実施形態では、複数の関心領域２４において、共通に設定された音速パラメータに基づいて、複数の関心領域画像４０を再構成する。例えば、音速変更ＧＵＩ３０を用いて、音速パラメータの変更をユーザから受け付け、当該変更された音速パラメータを複数の関心領域２４において共通する音速パラメータとして設定する。その後、共通に設定された音速パラメータに基づいて、複数の関心領域画像４０を再構成する。このような方法により、複数の関心領域２４において共通する音速パラメータの適切な値を早く見出すことができる。
また、このとき、例えば、実空間での位置関係と等しい位置関係で画像再構成したときの複数の関心領域２４同士の間隔よりも近づけて、複数の関心領域画像４０を所定の画面上に表示させる。すなわち、本実施形態では、仮の測定画像２０における複数の関心領域２４同士の間隔よりも近づけて、複数の関心領域画像４０を表示部７２０の所定の画面上に表示させる。複数の関心領域画像４０同士の間隔は、等しい間隔であってもよいし、それぞれ異なる間隔であってもよい。また、少なくとも２つの関心領域画像４０の一部が重なり合っていてもよい。一方で、複数の関心領域画像４０のうちで、近づけられていない２つ以上の関心領域画像４０があってもよい。上述のように、複数の関心領域画像４０を近づけることで、複数の関心領域画像４０同士を容易に比較しながら、音速パラメータを容易に調整することができる。
（Ｓ４８０：適切値決定工程）
上述の関心領域画像４０をユーザが確認し、関心領域画像４０の画質がＯＫである場合には（Ｓ４４０でＯＫ）、当該関心領域画像４０に基づいて音速パラメータの適切な値を決定する。
上述の関心領域画像４０をユーザが確認し、関心領域画像４０の画質がＯＫである場合には（Ｓ４４０でＯＫ）、当該関心領域画像４０に基づいて音速パラメータの適切な値を決定する。
音速パラメータの適切な値を決定するときの、関心領域画像４０の画質がＯＫである場合としては、例えば、所定の血管などの注目部位のコントラストが高い場合、関心領域画像４０の全体として各部位の視認性が一様に向上している場合（複数の血管の視認性が同一画像内で向上している場合）などが挙げられる。
また、このとき、本実施形態では、共通に設定された音速パラメータに基づいて再構成された複数の関心領域画像４０の画質を判断することで、複数の関心領域２４において、共通する音速パラメータの適切な値を決定する。これにより、全ての関心領域２４で許容可能な音速パラメータの適切な値を早く決定することができる。
上述のように音速パラメータの適切な値を決定したら、当該音速パラメータの適切な値を記憶装置７０４に保存する（Ｓ４８０）。
（Ｓ５４０：広域画像再構成工程）
次に、図８に示すように、上述の音速パラメータの適切な値に基づいて、関心領域２４よりも広い領域の広域画像５０を再構成する。
次に、図８に示すように、上述の音速パラメータの適切な値に基づいて、関心領域２４よりも広い領域の広域画像５０を再構成する。
具体的には、記憶装置７０４から音速パラメータの適切な値および信号データを読み出し、音速パラメータの適切な値および信号データに基づいて広域画像５０を再構成し、表示部７２０の所定の画面上に表示させる。
このとき、広域画像５０として再構成される領域を、例えば、関心領域２４よりも広く、かつ、関心領域２４を含む領域とすることが好ましい。例えば、図８に示すように、撮影領域の全体に亘って、広域画像５０を再構成することが好ましい。
このとき、信号データを減少させることなく、高画質の広域画像５０を再構成する。例えば、センサユニット３００の位置を移動させながら撮影領域全体に亘って所定の測定間隔で取得したそのままの全信号データに基づいて、すなわち、信号データを減少させることなく、広域画像５０を再構成する。
上述のように、関心領域２４の音速パラメータの適切な値が決まれば、当該音速パラメータの適切な値に基づいて被検部位１１０全体を再構成することによって、高画質な全体の再構成画像を得ることが可能となる。
なお、広域画像５０として再構成される領域が関心領域２４よりも広ければ、必ずしも撮影領域の全体に亘って広域画像５０を再構成しなくてもよく、撮影領域の一部に基づいて広域画像５０を再構成してもよい。この場合、広域画像５０の境界をマウスクリックで変更可能状態とし、広域画像５０のサイズおよび形状を変更可能にすることが好ましい。
以上により、本実施形態の画像生成工程を終了する。
（４）本実施形態により得られる効果
本実施形態によれば、以下に示す１つまたは複数の効果が得られる。
本実施形態によれば、以下に示す１つまたは複数の効果が得られる。
（ａ）本実施形態では、関心領域２４の指定を受け付けるために、仮の測定画像２０を再構成する。仮の測定画像２０においては、関心領域２４の位置を把握できればよいため、画質を低くして（粗くして）仮の測定画像２０を再構成することができる。すなわち、仮の測定画像２０を再構成するための単位面積当たりの計算量を、関心領域画像４０を再構成するための単位面積当たりの計算量よりも少なくすることができる。これにより、仮の測定画像２０を早く再構成することができる。
その後、関心領域２４において、単位面積当たりの計算量を多くして、関心領域画像４０を再構成する。このとき、仮の測定画像２０の領域よりも小さい領域を有する関心領域画像４０を再構成することで、音速パラメータを変更しながらであっても、高画質の関心領域画像４０を順次再構成することができる。これにより、関心領域画像４０の画質に基づいて、音速パラメータの適切な値を効率よく決定することができる。
これらの結果、上述のようにして決定された音速パラメータの適切な値に基づいて、広い範囲に亘って画質が良好な広域画像５０を早く再構成することが可能となる。
（ｂ）仮の測定画像再構成工程Ｓ２２０～Ｓ２４０では、仮の測定画像２０の同一画素２２の再構成に用いる信号データの重畳数を少なくする。これにより、仮の測定画像２０を再構成するときの単位面積当たりの信号データ量を少なくすることができる。その結果、仮の測定画像２０を再構成するための単位面積当たりの計算量を、関心領域画像４０を再構成するための単位面積当たりの計算量よりも少なくすることができる。
また、信号データの重畳数を少なくすることで、アーチファクトの程度は悪化するものの、仮の測定画像２０の解像度（画素２２の密度）を高く維持することができる。これにより、仮の測定画像２０における関心領域２４の位置の視認性低下を抑制することができる。
（ｃ）関心領域指定受付工程Ｓ３２０では、仮の測定画像２０内で、複数の関心領域２４の指定を受け付ける。これにより、ユーザが注目する複数の部位が互いに離れた位置にあっても、これらの部位で画質が良好となるように音速パラメータを調整することができる。また、互いに離れた位置にある複数の部位を含む大きな関心領域２４を指定する必要が無いので、関心領域画像４０の再構成を早く行うことができる。
また、複数の関心領域２４の一部を重ね合わせることができる。これにより、結合された複数の関心領域２４全体として、任意の形状に関心領域２４を調整することができる。例えば、血管などの複雑な形状に合わせて関心領域２４の形状を調整することができる。
このように、複数の関心領域２４の指定を受け付けることで、ユーザが優先する条件（位置または形状など）を満たす複数の関心領域画像４０を容易に得ることが可能となる。さらに、上述のようにして得られた複数の関心領域画像４０に基づいて音速パラメータを調整することで、ユーザの関心を充分に満たし、広い範囲に亘って画質が良好な広域画像５０を、早く再構成することが可能となる。
（ｄ）本実施形態では、複数の関心領域２４において、共通に設定された音速パラメータに基づいて、複数の関心領域画像４０を再構成し、当該複数の関心領域２４において、共通する音速パラメータの適切な値を決定する。
ここで、後述の第２実施形態のように、各場所に適切な値を探して、様々な場所での音速パラメータを適切な値に調整することは、画質の観点では理想的ではあるものの、処理部７００の計算負荷が大きくなるのが課題である。
これに対し、本実施形態では、共通に設定された音速パラメータに基づいて、複数の関心領域画像４０を再構成することにより、複数の場所の全てで、ベストではないかもしれないが許容可能な画質が得られるように、音速パラメータの適切な値を決定することが可能である。すなわち、一つの音速パラメータの適切な値に基づいて、あらゆる場所で画質が許容可能な広域画像５０を容易に得ることができる。その結果、処理部７００の計算負荷を軽減することが可能となる。
（５）本実施形態の変形例
上述の実施形態の仮の測定画像再構成工程Ｓ２２０～Ｓ２４０では、仮の測定画像２０の同一画素の再構成に用いる信号データの重畳数を少なくする場合について説明したが、必要に応じて、以下に示す変形例のように変更することができる。
上述の実施形態の仮の測定画像再構成工程Ｓ２２０～Ｓ２４０では、仮の測定画像２０の同一画素の再構成に用いる信号データの重畳数を少なくする場合について説明したが、必要に応じて、以下に示す変形例のように変更することができる。
以下、上述の実施形態と異なる要素についてのみ説明し、上述の実施形態で説明した要素と実質的に同一の要素には、同一の符号を付してその説明を省略する。なお、以下の第２実施形態～第４実施形態などについても、本変形例と同様に説明を省略する。
＜変形例１＞
変形例１の仮の測定画像再構成工程Ｓ２２０～Ｓ２４０では、例えば、仮の測定画像２０の画素サイズを大きくする。具体的には、仮の測定画像２０の画素サイズを、例えば、被検部位１１０内の複数の単位空間領域に相当するサイズに大きくする。このように仮の測定画像２０の画素サイズを大きくすることで、仮の測定画像２０の解像度は低くなる。また、仮の測定画像２０では、大きくした１つの画素２２を、代表的な１つの単位領域に相当する信号のみに基づいて再構成する。
変形例１の仮の測定画像再構成工程Ｓ２２０～Ｓ２４０では、例えば、仮の測定画像２０の画素サイズを大きくする。具体的には、仮の測定画像２０の画素サイズを、例えば、被検部位１１０内の複数の単位空間領域に相当するサイズに大きくする。このように仮の測定画像２０の画素サイズを大きくすることで、仮の測定画像２０の解像度は低くなる。また、仮の測定画像２０では、大きくした１つの画素２２を、代表的な１つの単位領域に相当する信号のみに基づいて再構成する。
変形例１によれば、仮の測定画像２０の画素サイズを大きくすることで、仮の測定画像２０を再構成するときの単位面積当たりの信号データ量を少なくすることができる。例えば、ＸＹＺ方向にそれぞれに画素サイズ（ボクセルサイズ）を大きくすることで、単位面積当たりの信号データ量を容易に少なくすることができる。これにより、仮の測定画像２０を再構成するときの単位面積当たりの計算量を少なくすることができる。その結果、仮の測定画像２０の再構成速度を向上させることができる。
具体的には、例えば、画素２２としてのボクセルの一辺の長さを０．１ｍｍから０．４ｍｍに変更した場合であっても、表示部７２０の画面上の視認性をあまり変化させることなく、仮の測定画像２０の再構成速度を大幅に向上させることが可能となる。
＜変形例２＞
変形例２の仮の測定画像再構成工程Ｓ２２０～Ｓ２４０では、例えば、被検部位１１０の深さ方向（すなわちＺ方向）における信号データの範囲を狭くする。例えば、被検部位１１０において血管が集中する表面に近い部分に、信号データの範囲を限定する。なお、このとき、画素サイズは変化させない。
変形例２の仮の測定画像再構成工程Ｓ２２０～Ｓ２４０では、例えば、被検部位１１０の深さ方向（すなわちＺ方向）における信号データの範囲を狭くする。例えば、被検部位１１０において血管が集中する表面に近い部分に、信号データの範囲を限定する。なお、このとき、画素サイズは変化させない。
変形例２によれば、Ｚ方向における信号データの範囲を狭くすることで、ＸＹ面の単位面積当たりでＺ方向の範囲全体を積算した総信号データ量を少なくすることができる。すなわち、仮の測定画像２０を再構成するときの単位面積当たりの信号データ量を少なくすることができる。これにより、仮の測定画像２０を再構成するときの単位面積当たりの計算量を少なくすることができる。その結果、仮の測定画像２０の再構成速度を向上させることができる。
また、限られたＺ方向の範囲内では、解像度を高く維持することができる。これにより、高解像な断面像（に近い像）としての仮の測定画像２０を再構成することができる。その結果、仮の測定画像２０内で関心領域２４を指定する際の位置視認精度を向上させることができる。
＜変形例３＞
変形例３の仮の測定画像再構成工程Ｓ２２０～Ｓ２４０では、例えば、複数のセンサ３４０により受信した複数の信号のうちの一部に基づいて、仮の測定画像２０を再構成する。なお、このとき、画素サイズは変化させない。
変形例３の仮の測定画像再構成工程Ｓ２２０～Ｓ２４０では、例えば、複数のセンサ３４０により受信した複数の信号のうちの一部に基づいて、仮の測定画像２０を再構成する。なお、このとき、画素サイズは変化させない。
このとき、仮の測定画像２０を再構成するための信号が選択される一部のセンサ３４０は、例えば、半球状の素子保持部３６０の一か所に集中しているよりは、素子保持部３６０の全体に亘って分散していることが好ましく、素子保持部３６０の全体に亘って所定の等しい間隔で互いに離れて配置されていることがより好ましい。例えば、センサ３４０の数を１／２にする場合に、半球状の素子保持部３６０の底面付近のみを選択し、素子保持部３６０の側面のセンサ３４０を選択しない場合、素子保持部３６０を半球状にした効果が得られ難くなる。これに対し、選択される一部のセンサ３４０を、素子保持部３６０の全体に亘って分散させることで、素子保持部３６０を半球状にした効果を得ることができる。
変形例３によれば、複数のセンサ３４０により受信した複数の信号のうちの一部に信号データを限定することで、仮の測定画像２０を再構成するときの単位面積当たりの信号データ量を少なくすることができる。これにより、仮の測定画像２０を再構成するときの単位面積当たりの計算量を少なくすることができる。
また、変形例３によれば、仮の測定画像２０を再構成するための信号が選択される一部のセンサ３４０を素子保持部３６０の全体に亘って分散させることで、上述のように、素子保持部３６０を半球状にした効果を得ることができる。すなわち、広い立体角で受信した信号に基づいて、仮の測定画像２０を再構成することができる。その結果、アーチファクトの程度は悪化するものの、仮の測定画像２０の解像度（画素２２の密度）を高く維持することが可能となる。
＜変形例４＞
変形例４の仮の測定画像再構成工程Ｓ２２０～Ｓ２４０では、例えば、以下のように、仮の測定画像２０を再構成するための計算量が少ない方法、すなわち簡易的な処理方法を採用する。
変形例４の仮の測定画像再構成工程Ｓ２２０～Ｓ２４０では、例えば、以下のように、仮の測定画像２０を再構成するための計算量が少ない方法、すなわち簡易的な処理方法を採用する。
例えば、ノイズ除去のためのフィルタリングを実空間のフィルタで行う。ここで、通常の画像再構成処理、即ち、関心領域画像再構成工程Ｓ４２０および広域画像再構成工程Ｓ５４０では、ノイズ除去のためのフィルタリングを周波数空間で行っている。しかしながら、周波数空間でのフィルタリングでは、フーリエ変換およびフーリエ逆変換が必要となる。このため、画像再構成のための計算量が多くなる傾向がある。これに対し、本変形例では、仮の測定画像２０の再構成処理のノイズ除去のためのフィルタリングを実空間のフィルタで行うことで、フーリエ変換およびフーリエ逆変換を不要とすることができる。このような実空間のフィルタリングにより、信号データ量を減らさずとも、計算量の増加を抑制することができる。また、このような手法であっても、関心領域２４の指定を受け付けるための仮の測定画像２０の画質を、充分に確保することができる。これらの結果、所望の仮の測定画像２０を高速に得ることができる。
変形例４によれば、ノイズ除去のためのフィルタリングを実空間のフィルタで行うことで仮の測定画像２０を再構成する。これにより、仮の測定画像２０を再構成するときの単位面積当たりの計算量を少なくすることができる。その結果、高速で仮の測定画像２０を得ることが可能となる。
なお、仮の測定画像２０を高速に得る方法は、以上述べた実施形態、変形例１から４に限らず、例えば、データサンプリングの周波数を下げたり、上位ビットのみを用いたりしても良い。また、これらの実施形態、変形例１から４などを組合せても構わない。
＜本発明の第２実施形態＞
次に、本発明の第２実施形態について説明する。
次に、本発明の第２実施形態について説明する。
（１）画像生成方法
［本実施形態の概略］
本実施形態では、例えば、複数の関心領域２４において、それぞれ個別に設定された音速パラメータに基づいて、複数の関心領域画像４０を再構成する。
［本実施形態の概略］
本実施形態では、例えば、複数の関心領域２４において、それぞれ個別に設定された音速パラメータに基づいて、複数の関心領域画像４０を再構成する。
［本実施形態の具体的な画像生成方法］
図９を用い、本実施形態の具体的な画像生成方法について説明する。
図９を用い、本実施形態の具体的な画像生成方法について説明する。
（Ｓ３２０：関心領域指定受付工程）
仮の測定画像２０を再構成したら、仮の測定画像２０内で、部分的な関心領域２４の指定を受け付ける。本実施形態では、当該工程の１回において、複数の関心領域２４の指定を受け付けてもよいし、１つの関心領域２４のみの指定を受け付けてもよい。
仮の測定画像２０を再構成したら、仮の測定画像２０内で、部分的な関心領域２４の指定を受け付ける。本実施形態では、当該工程の１回において、複数の関心領域２４の指定を受け付けてもよいし、１つの関心領域２４のみの指定を受け付けてもよい。
（Ｓ４２０～Ｓ４６０：関心領域画像再構成工程）
次に、関心領域２４において、音速パラメータを変更しながら、当該変更された音速パラメータに基づいて、音響画像からなる関心領域画像４０を順次再構成する。
次に、関心領域２４において、音速パラメータを変更しながら、当該変更された音速パラメータに基づいて、音響画像からなる関心領域画像４０を順次再構成する。
このとき、本実施形態では、上述の関心領域２４において、個別に設定された音速パラメータに基づいて、関心領域画像４０を再構成する。
具体的には、例えば、個別に設定された音速パラメータに基づいて、１つの関心領域画像４０を再構成し、表示部７２０の所定の画面上に表示させる（Ｓ４２０）。次に、再構成された１つの関心領域画像４０の画質をユーザにより判断する（Ｓ４４０）。当該関心領域画像４０の画質がＮＧである場合には（Ｓ４４０でＮＧ）、１つの関心領域画像４０に個別に設定する音速パラメータの任意の変更をユーザから受け付ける（Ｓ４６０）。
（Ｓ４８０：適切値決定工程）
上述の関心領域画像４０をユーザが確認し、関心領域画像４０の画質がＯＫである場合には（Ｓ４４０でＯＫ）、１つの関心領域画像４０に対して個別に音速パラメータの適切な値を決定し、記憶装置７０４に記憶する。
上述の関心領域画像４０をユーザが確認し、関心領域画像４０の画質がＯＫである場合には（Ｓ４４０でＯＫ）、１つの関心領域画像４０に対して個別に音速パラメータの適切な値を決定し、記憶装置７０４に記憶する。
（Ｓ４９０：全関心領域終了判断工程）
次に、全ての関心領域２４において、それぞれ個別に音速パラメータの適切な値を決定したか否かを判断する。
次に、全ての関心領域２４において、それぞれ個別に音速パラメータの適切な値を決定したか否かを判断する。
（Ｓ３４０：別の関心領域指定受付工程）
全ての関心領域２４において、それぞれ個別に音速パラメータの適切な値を決定していない場合（すなわち、適切な値が決定していない関心領域２４が残っている場合）には（Ｓ４９０でＮｏ）、仮の測定画像２０内で、別の関心領域２４の指定を受け付ける。
全ての関心領域２４において、それぞれ個別に音速パラメータの適切な値を決定していない場合（すなわち、適切な値が決定していない関心領域２４が残っている場合）には（Ｓ４９０でＮｏ）、仮の測定画像２０内で、別の関心領域２４の指定を受け付ける。
別の関心領域２４の指定を受け付けたら、当該別の関心領域２４について、関心領域画像再構成工程Ｓ４２０～Ｓ４６０を行う。
（Ｓ５２０：データ読み出し工程）
上述のＳ４２０～Ｓ４９０を含むサイクルを繰り返した後に、全ての関心領域２４において、それぞれ個別に音速パラメータの適切な値を決定した場合には（Ｓ４９０でＹｅｓ）、広域画像５０の再構成のために、信号データと、全ての関心領域２４にそれぞれ個別に決定された音速パラメータの適切な値と、を記憶装置７０４から読み出す。
上述のＳ４２０～Ｓ４９０を含むサイクルを繰り返した後に、全ての関心領域２４において、それぞれ個別に音速パラメータの適切な値を決定した場合には（Ｓ４９０でＹｅｓ）、広域画像５０の再構成のために、信号データと、全ての関心領域２４にそれぞれ個別に決定された音速パラメータの適切な値と、を記憶装置７０４から読み出す。
（Ｓ５４０：広域画像再構成工程）
次に、上述の音速パラメータの適切な値に基づいて、広域画像５０を再構成する。本実施形態では、上述のように音速パラメータの適切な値が関心領域２４ごとに決定されているため、広域画像５０の再構成方法としては、例えば、以下の２つの方法が考えられる。
次に、上述の音速パラメータの適切な値に基づいて、広域画像５０を再構成する。本実施形態では、上述のように音速パラメータの適切な値が関心領域２４ごとに決定されているため、広域画像５０の再構成方法としては、例えば、以下の２つの方法が考えられる。
例えば、まず、複数の関心領域２４をそれぞれ含む複数の領域に広域画像５０を区画する。広域画像５０を区画したら、区画された複数の領域においてそれぞれ個別に決定された音速パラメータの適切な値を適用し、区画された各領域の画像を再構成する。その後、区画された各領域の画像を結合し、全体としての広域画像５０を生成する。
或いは、例えば、全ての関心領域２４にそれぞれ個別に決定された音速パラメータの適切な値に基づいて、音速パラメータの１つの選択値を決定してもよい。１つの選択値としては、例えば、複数の音速パラメータのうちの代表値、複数の音速パラメータの最小値以上最大値以下の範囲内の１つの値、複数の音速パラメータの平均値などが挙げられる。このように決定された音速パラメータの選択値に基づいて、広域画像５０を再構成してもよい。
（Ｓ５６０：広域画像再構成工程）
次に、必要に応じて、広域画像５０内に、音速パラメータの適切な値の分布を表示させる。音速パラメータの適切な値の分布を表示させる方法としては、例えば、音速パラメータの所定範囲ごとに色を異ならせる方法、音速パラメータの所定範囲ごとに色の濃さを異ならせる方法などが挙げられる。
次に、必要に応じて、広域画像５０内に、音速パラメータの適切な値の分布を表示させる。音速パラメータの適切な値の分布を表示させる方法としては、例えば、音速パラメータの所定範囲ごとに色を異ならせる方法、音速パラメータの所定範囲ごとに色の濃さを異ならせる方法などが挙げられる。
（２）本実施形態により得られる効果
本実施形態では、複数の関心領域２４において、それぞれ個別に設定された音速パラメータに基づいて、複数の関心領域画像４０を再構成する。これにより、複数の関心領域２４においてそれぞれ個別に音速パラメータの適切な値を決定することができる。その結果、被検部位１１０の場所によって音速パラメータのばらつきが生じていたとしても、それぞれの場所に応じた音速パラメータの適切な値に基づいて、全体に亘って画質が良好な広域画像５０を再構成することが可能となる。
本実施形態では、複数の関心領域２４において、それぞれ個別に設定された音速パラメータに基づいて、複数の関心領域画像４０を再構成する。これにより、複数の関心領域２４においてそれぞれ個別に音速パラメータの適切な値を決定することができる。その結果、被検部位１１０の場所によって音速パラメータのばらつきが生じていたとしても、それぞれの場所に応じた音速パラメータの適切な値に基づいて、全体に亘って画質が良好な広域画像５０を再構成することが可能となる。
更に、それぞれの場所ごとに個別に音速パラメータを変更しながら、関心領域画像４０を順次再構成することで、被検部位１１０の面内方向による違いだけでなく、被検部位１１０の奥行き方向の違いも、画像再構成中にリアルタイムに把握することが可能となる。
また、画像生成プログラム中に、それぞれの場所ごとに個別に決定した音速パラメータの適切な値を記録する機能を付与することによって、被検体１００の音速パラメータ分布を求めて、画像として提示することも可能となる。例えば、軟組織中に癌が発生した状態を想定すると、軟組織中の音速は遅いが、組織が硬い癌の中の音速は速い。こうした軟組織中の癌に起因する音速の違いを描出することも可能となる。
＜本発明の第３実施形態＞
次に、本発明の第３実施形態について説明する。
次に、本発明の第３実施形態について説明する。
（１）画像生成方法
［本実施形態の概略］
本実施形態では、例えば、音速パラメータの変更範囲および変更幅のうち少なくともいずれかの変更条件を受け付け、その後、当該変更条件下で音速パラメータを自動で変更しながら、当該変更された音速パラメータに基づいて、関心領域画像４０を順次再構成する。
［本実施形態の概略］
本実施形態では、例えば、音速パラメータの変更範囲および変更幅のうち少なくともいずれかの変更条件を受け付け、その後、当該変更条件下で音速パラメータを自動で変更しながら、当該変更された音速パラメータに基づいて、関心領域画像４０を順次再構成する。
［本実施形態の具体的な画像生成方法］
図１０を用い、本実施形態の具体的な画像生成方法について説明する。
図１０を用い、本実施形態の具体的な画像生成方法について説明する。
（Ｓ４１２～Ｓ４４０：関心領域画像再構成工程）
仮の測定画像２０内で関心領域２４の指定を受け付けたら、後述の関心領域画像４０の再構成の際に用いる音速パラメータの変更範囲および変更幅のうち少なくともいずれかの変更条件を受け付ける（Ｓ４１２）。
仮の測定画像２０内で関心領域２４の指定を受け付けたら、後述の関心領域画像４０の再構成の際に用いる音速パラメータの変更範囲および変更幅のうち少なくともいずれかの変更条件を受け付ける（Ｓ４１２）。
すなわち、本実施形態では、音速パラメータを変更するための音速変更ＧＵＩ３０を用いて、ユーザが音速パラメータを変えるのではなく、音速パラメータの変更範囲と変更幅をあらかじめ決めておいて、その後、自動的に音速パラメータを変化させる。ここでいう「変更範囲」とは、音速パラメータを変化させるときの最小値以上最大値以下の範囲のことを意味し、「変更幅」とは、音速パラメータを自動でインクリメントさせる際の刻み幅のことを意味する。
音速パラメータの変更条件を受け付ける方法としては、例えば、第１実施形態での音速変更ＧＵＩ３０による方法と同様に、テキストボックス内に数値を直接入力する方法、数値を増減するボタンをクリックする方法、スライダバーをクリックしたままドラッグすることで条件を変更する方法、マウスホイールを回転させて条件を変更する方法などが挙げられる。なお、上記方法に限られるものではなく、変更条件に係る数値を変化させることができるのであれば、いずれの方法を用いてもよい。あるいは、ＧＵＩによる方法ではなく、設定ファイルを別途用意してその設定情報ファイルを読み込む方法を用いてもよい。
音速パラメータの変更条件を受け付けたら、当該変更条件下で音速パラメータを自動で変更しながら、当該変更された音速パラメータに基づいて、関心領域画像４０を順次再構成する。
具体的には、音速パラメータの変更条件の受け付け後、ユーザが表示部７２０の画面上でスタートボタンを押すなどの操作をトリガーとして、自動的に音速パラメータを変更し、所定の音速パラメータを設定する（Ｓ４１６）。自動で変更された音速パラメータに基づいて関心領域画像４０を再構成し、表示部７２０の所定の画面上に表示させる（Ｓ４２０）。次に、再構成された関心領域画像４０の画質をユーザにより判断する（Ｓ４４０）。関心領域画像４０の画質がＮＧである場合には（Ｓ４４０でＮＧ）、音速パラメータの自動変更および設定を再度行い（Ｓ４１０）、それ以降のサイクルを繰り返す。
（Ｓ４７０：適切値決定工程）
ユーザは、連続的に変化する関心領域画像４０を見ながら、関心領域画像４０の画質がＯＫな状態となったタイミングで（Ｓ４４０でＯＫ）、表示部７２０の画面上でストップボタンを押す。ユーザによる当該操作をトリガーとして、音速パラメータの自動変更と、関心領域画像４０の順次再構成と、を停止させる。このように自動で変更される音速パラメータを停止させたときの当該音速パラメータを、音速パラメータの適切な値として決定する。
ユーザは、連続的に変化する関心領域画像４０を見ながら、関心領域画像４０の画質がＯＫな状態となったタイミングで（Ｓ４４０でＯＫ）、表示部７２０の画面上でストップボタンを押す。ユーザによる当該操作をトリガーとして、音速パラメータの自動変更と、関心領域画像４０の順次再構成と、を停止させる。このように自動で変更される音速パラメータを停止させたときの当該音速パラメータを、音速パラメータの適切な値として決定する。
（２）本実施形態により得られる効果
本実施形態では、ユーザによって指定された変更条件下で音速パラメータを自動で変更しながら、関心領域画像４０を順次再構成する。ユーザは、最初に変更条件を指定するだけでよく、関心領域画像４０の再構成の際に、第１実施形態のような音速変更ＧＵＩ３０を入力または操作する必要がない。つまり、音速パラメータの変更におけるユーザの操作負担を軽減することが可能となる。
本実施形態では、ユーザによって指定された変更条件下で音速パラメータを自動で変更しながら、関心領域画像４０を順次再構成する。ユーザは、最初に変更条件を指定するだけでよく、関心領域画像４０の再構成の際に、第１実施形態のような音速変更ＧＵＩ３０を入力または操作する必要がない。つまり、音速パラメータの変更におけるユーザの操作負担を軽減することが可能となる。
また、本実施形態では、ユーザが関心領域画像４０の画質に基づいて音速パラメータの自動変更を停止するだけで、音速パラメータの適切な値を決定することができる。つまり、音速パラメータの適切な値の決定においても、ユーザの負担を軽減することが可能となる。
（３）本実施形態の変形例
上述の実施形態の関心領域画像再構成工程Ｓ４１２～Ｓ４４０では、自動で変更された音速パラメータに基づいて関心領域画像４０を順次再構成して所定の画面上に表示させる場合について説明したが、この場合に限られない。
上述の実施形態の関心領域画像再構成工程Ｓ４１２～Ｓ４４０では、自動で変更された音速パラメータに基づいて関心領域画像４０を順次再構成して所定の画面上に表示させる場合について説明したが、この場合に限られない。
例えば、ユーザが表示部７２０の画面上でスタートボタンを押すと、それぞれの音速パラメータに基づいて得られた関心領域２４のボリュームデータのＭＩＰ画像が。静止画像として記憶装置７０４に保存される。この静止画は、汎用の画像ビューワを用いて閲覧することが可能である。この機能によって、第３実施形態で決定した音速パラメータを事後に検証することが可能となる。
なお、本変形例で述べた音速パラメータごとにＭＩＰ画像が記憶装置７０４に保存される機能は、他の実施形態においても追加機能として採用することが可能である。
＜本発明の第４実施形態＞
次に、本発明の第４実施形態について説明する。
次に、本発明の第４実施形態について説明する。
（１）画像生成方法
［本実施形態の概略］
本実施形態では、例えば、ユーザがスタートボタンを押すと、それぞれの音速パラメータで得られた関心領域２４について、自動的に音速パラメータの適切な値を決定できるオートフォーカス機能を搭載する。このオートフォーカスは、ＲＡＭ７０３上の画像を用いてリアルタイムに解析することが好ましい。ただし、ＭＩＰ画像を静止画像として記憶装置７０４に保存し、一連の静止画像の保存が完了したのち、画像解析によって自動的にフォーカスを決定する（すなわち音速パラメータの適切な値を決定する）方法を採用してもよい。このオートフォーカスを実現するために、ユーザは関心領域２４を指定するとともに、関心領域２４内のどの画像（部位）に注目すべきか（対象が血管の場合、どの血管に注目してオートフォーカス処理を行うか）を指定する。
［本実施形態の概略］
本実施形態では、例えば、ユーザがスタートボタンを押すと、それぞれの音速パラメータで得られた関心領域２４について、自動的に音速パラメータの適切な値を決定できるオートフォーカス機能を搭載する。このオートフォーカスは、ＲＡＭ７０３上の画像を用いてリアルタイムに解析することが好ましい。ただし、ＭＩＰ画像を静止画像として記憶装置７０４に保存し、一連の静止画像の保存が完了したのち、画像解析によって自動的にフォーカスを決定する（すなわち音速パラメータの適切な値を決定する）方法を採用してもよい。このオートフォーカスを実現するために、ユーザは関心領域２４を指定するとともに、関心領域２４内のどの画像（部位）に注目すべきか（対象が血管の場合、どの血管に注目してオートフォーカス処理を行うか）を指定する。
［本実施形態の具体的な画像生成方法］
図１１を用い、本実施形態の具体的な画像生成方法について説明する。
図１１を用い、本実施形態の具体的な画像生成方法について説明する。
（Ｓ４１２～Ｓ４５４：関心領域画像再構成工程）
仮の測定画像２０内で関心領域２４の指定を受け付けたら、関心領域画像４０の再構成の際に用いる音速パラメータの変更範囲および変更幅のうち少なくともいずれかの変更条件を受け付ける（Ｓ４１２）。
仮の測定画像２０内で関心領域２４の指定を受け付けたら、関心領域画像４０の再構成の際に用いる音速パラメータの変更範囲および変更幅のうち少なくともいずれかの変更条件を受け付ける（Ｓ４１２）。
音速パラメータの変更条件を受け付けたら、関心領域２４内で被検部位１１０の注目部位の指定を受け付ける（Ｓ４１４）。ここでいう「注目部位」とは、例えば、所定の血管およびその周辺部位などのことを意味する。注目部位の指定受付方法としては、特に限定されるものではないが、例えば、ユーザが注目部位をポインタでクリックする方法などが挙げられる。
注目部位の指定を受け付けたら、上述のオートフォーカスを行う。すなわち、音速パラメータの変更条件下で音速パラメータを自動で変更しながら、当該変更された音速パラメータに基づいて、関心領域画像４０を順次再構成する。このように関心領域画像４０を順次再構成しながら、関心領域画像４０の解像度を除く画質に基づいて、音速パラメータの適切な値を自動で決定する。例えば、関心領域画像４０における注目部位の解像度を除く画質に基づいて、音速パラメータの適切な値を自動で決定する。
具体的には、音速パラメータの変更条件の受け付け後、ユーザが表示部７２０の画面上でスタートボタンを押すなどの操作をトリガーとして、自動的に音速パラメータを変更し、所定の音速パラメータを設定する（Ｓ４１６）。自動で変更された音速パラメータに基づいて関心領域画像４０を再構成し、表示部７２０の所定の画面上に表示させる（Ｓ４２０）。次に、順次再構成される関心領域画像４０において、注目部位の解像度を除く画質を自動で解析し、当該自動解析結果が良好か否かを判断する（Ｓ４５２）。
オートフォーカスにおける自動解析処理の一例としては、以下の処理が考えられる。まず、指定された複数の血管とその周辺の部位に関して、初期値の音速パラメータで再構成された関心領域画像４０でのコントラストを求める。次に、コントラストの総和値を求めて、ＲＡＭ７０３に記憶させる。次に、音速パラメータを自動で変更し、関心領域画像４０の再構成を行い、再度コントラストの総和値を求める。最新のコントラストの総和値として、前回処理のコントラストの総和値よりも高い値が得られれば、そのときの音速パラメータをより好ましい音速パラメータとして記憶する。当該サイクルを繰り返すことにより、コントラストの総和値が最大となったときの音速パラメータを音速パラメータの適切な値として決定し、その値を記憶装置７０４に記憶させる。また、音速パラメータの適切な値に基づいて関心領域画像４０を再構成し、表示部７２０の所定の画面上に表示させる。以上のようにして、例えば、注目部位と背景とのコントラストが高く、エッジがシャープになるなどの関心領域画像４０の特徴を利用して、画像処理によって最適フォーカス位置（すなわち、音速パラメータの適切な値）を検出することができる。
しかしながら、オートフォーカスにおける自動解析処理では、画質の自動判定がしばしばローカルミニマムに落ち込んで真の最適条件に到達しないことがある。このため、オートフォーカス後に誤った数値が出力されたときに修正して再計算させる機能を付与しておくことが好ましい。
具体的には、注目部位の画質の自動解析結果がＮＧであった場合には（Ｓ４５２でＮＧ）、音速パラメータの自動変更および設定を再度行い（Ｓ４１６）、それ以降のサイクルを繰り返す。
一方で、注目部位の画質の自動解析結果がＯＫであった場合には（Ｓ４５２でＮＧ）、ユーザが関心領域画像４０を確認する。関心領域画像４０がＮＧであったとユーザが判断した場合には（Ｓ４５４でＮＧ）、音速パラメータの変更条件を再度受け付け（Ｓ４１２）、それ以降のサイクルを繰り返す。
（Ｓ５４０：広域画像再構成工程）
関心領域画像４０がＯＫであったとユーザが判断した場合には（Ｓ４５４でＯＫ）、上述のオートフォーカスによって決定された音速パラメータの適切な値に基づいて、広域画像５０を再構成する。
関心領域画像４０がＯＫであったとユーザが判断した場合には（Ｓ４５４でＯＫ）、上述のオートフォーカスによって決定された音速パラメータの適切な値に基づいて、広域画像５０を再構成する。
（２）本実施形態により得られる効果
本実施形態では、オートフォーカス機能により、関心領域画像４０において注目部位の解像度を除く画質に基づいて、音速パラメータの適切な値を自動で決定する。これにより、音速パラメータの適切な値の決定におけるユーザの負担を軽減することが可能となる。
本実施形態では、オートフォーカス機能により、関心領域画像４０において注目部位の解像度を除く画質に基づいて、音速パラメータの適切な値を自動で決定する。これにより、音速パラメータの適切な値の決定におけるユーザの負担を軽減することが可能となる。
また、注目部位の画質に基づいて音速パラメータの適切な値を自動で決定することで、ユーザが目視確認することなく、注目部位の視認性を向上させた広域画像５０を容易に再構成することができる。
また、注目部位の画質を自動で解析するときの基準を数値化しておけば、被検体１００が代わった場合であっても、複数の被検体１００を撮影した広域画像５０同士を、共通する画質条件下で比較することが可能となる。
＜他の実施形態＞
以上、本発明の実施形態を具体的に説明した。しかしながら、本発明は上述の実施形態に限定されるものではなく、その要旨を逸脱しない範囲で種々変更可能である。なお、以下において、特に限定せずにいう「上述の実施形態」とは、全ての実施形態を含んでいる。
以上、本発明の実施形態を具体的に説明した。しかしながら、本発明は上述の実施形態に限定されるものではなく、その要旨を逸脱しない範囲で種々変更可能である。なお、以下において、特に限定せずにいう「上述の実施形態」とは、全ての実施形態を含んでいる。
上述の実施形態では、一つの音速パラメータで媒体中を音波が伝搬することを仮定した再構成処理を述べてきたが、この場合に限られない。従来の超音波装置で報告がなされているように、多層構成の場合にそれぞれの音速パラメータを用いて画像再構成を行ってもよい。
上述の実施形態では、光音響イメージング装置１０がＰＡＴ装置として構成されている場合について説明したが、光音響イメージング装置１０は、音響波を測定可能であれば、ＰＡＴ装置以外の装置として構成されていてもよい。例えば、光音響イメージング装置１０は、被検体１００の所定の被検部位１１０に対して音響波（超音波）を照射し、照射された部分から反射または散乱された音響波（反射波）を受信する超音波エコーイメージング装置として構成されていてもよい。
上述の実施形態では、光音響イメージング装置１０が音響波を測定する機能と画像生成を行う機能との両方を有する場合について説明したが、上述した画像生成処理を実行する画像生成装置は、音響波測定装置と別に設けられていてもよい。言い換えれば、上述の画像生成処理は、音響波測定装置に接続されているコンピュータを用いてもよいし、別のコンピュータを用いてもよい。
上述の第３実施形態および第４実施形態では、ユーザから受け付けた音速パラメータの変更条件下で音速パラメータを自動で変更しながら、関心領域画像４０を順次再構成する場合について説明したが、この場合に限られない。例えば、予め設定された変更範囲および変更幅で音速パラメータを自動で変更しながら、当該変更された音速パラメータに基づいて、関心領域画像４０を順次再構成してもよい。これにより、ユーザの負担をさらに軽減することができる。
上述の実施形態では、仮の測定画像２０内で、複数の関心領域２４の指定を受け付ける場合について説明したが、この場合に限られない。例えば、仮の測定画像２０を再構成する際に、仮の測定画像２０を再構成するための単位面積当たりの計算量を、関心領域画像４０を再構成するための単位面積当たりの計算量よりも少なくするだけでも、画質が良好な広域画像５０を早く再構成することが可能となる。したがって、この場合では、仮の測定画像２０内で、１つの関心領域２４だけの指定を受け付けてもよい。
上述の実施形態では、仮の測定画像２０を再構成する際に、仮の測定画像２０を再構成するための単位面積当たりの計算量を、関心領域画像４０を再構成するための単位面積当たりの計算量よりも少なくする場合について説明したが、この場合に限られない。例えば、被検部位１１０の撮影領域内で、複数の関心領域２４の指定を受け付け、複数の関心領域画像４０を再構成するだけでも、画質が良好な広域画像５０を早く再構成することが可能となる。したがって、この場合では、仮の測定画像２０を再構成するための単位面積当たりの計算量を少なくせずに、仮の測定画像２０を再構成してもよい。また、この場合では、複数の関心領域２４の指定を受け付ける際には、領域指定のための信号データと最終的な画像再構成のための信号データとが一致する必要は無い。このため、必ずしも、上述の実施形態で説明した仮の測定画像２０内で、複数の関心領域２４の指定を受け付けなくてもよい。
＜本発明の好ましい態様＞
以下、本発明の好ましい態様について付記する。
以下、本発明の好ましい態様について付記する。
（付記１）
所定の被検部位を測定した音響波の信号データに基づいて、音響画像を生成する画像生成方法であって、
前記信号データに基づいて、音響画像からなる仮の測定画像を再構成する工程と、
前記仮の測定画像内で部分的な関心領域の指定を受け付ける工程と、
前記関心領域において、音速パラメータを変更しながら、当該変更された音速パラメータに基づいて、音響画像からなる関心領域画像を順次再構成する工程と、
前記関心領域画像に基づいて前記音速パラメータの適切な値を決定する工程と、
前記音速パラメータの前記適切な値に基づいて、前記関心領域よりも広い領域の音響画像を再構成する工程と、
を有し、
前記仮の測定画像を再構成する工程では、
前記仮の測定画像を再構成するための単位面積当たりの計算量を、前記関心領域画像を再構成するための単位面積当たりの計算量よりも少なくする
画像生成方法。
所定の被検部位を測定した音響波の信号データに基づいて、音響画像を生成する画像生成方法であって、
前記信号データに基づいて、音響画像からなる仮の測定画像を再構成する工程と、
前記仮の測定画像内で部分的な関心領域の指定を受け付ける工程と、
前記関心領域において、音速パラメータを変更しながら、当該変更された音速パラメータに基づいて、音響画像からなる関心領域画像を順次再構成する工程と、
前記関心領域画像に基づいて前記音速パラメータの適切な値を決定する工程と、
前記音速パラメータの前記適切な値に基づいて、前記関心領域よりも広い領域の音響画像を再構成する工程と、
を有し、
前記仮の測定画像を再構成する工程では、
前記仮の測定画像を再構成するための単位面積当たりの計算量を、前記関心領域画像を再構成するための単位面積当たりの計算量よりも少なくする
画像生成方法。
（付記２）
前記仮の測定画像を再構成する工程では、
前記仮の測定画像の画素サイズを大きくする
付記１に記載の画像生成方法。
前記仮の測定画像を再構成する工程では、
前記仮の測定画像の画素サイズを大きくする
付記１に記載の画像生成方法。
（付記３）
前記信号データとして、同一の画素に再構成される位置に複数の信号が重畳されたデータを準備する工程をさらに有し、
前記仮の測定画像を再構成する工程では、
前記仮の測定画像の同一画素の再構成に用いる前記信号データの重畳数を少なくする
付記１に記載の画像生成方法。
前記信号データとして、同一の画素に再構成される位置に複数の信号が重畳されたデータを準備する工程をさらに有し、
前記仮の測定画像を再構成する工程では、
前記仮の測定画像の同一画素の再構成に用いる前記信号データの重畳数を少なくする
付記１に記載の画像生成方法。
（付記４）
前記仮の測定画像を再構成する工程では、
前記被検部位の深さ方向における前記信号データの範囲を狭くする
付記１に記載の画像生成方法。
前記仮の測定画像を再構成する工程では、
前記被検部位の深さ方向における前記信号データの範囲を狭くする
付記１に記載の画像生成方法。
（付記５）
前記信号データとして、複数のセンサにより受信した複数の信号を含むデータを準備する工程をさらに有し、
前記仮の測定画像を再構成する工程では、
前記複数の信号のうちの一部に基づいて前記仮の測定画像を再構成する
付記１に記載の画像生成方法。
前記信号データとして、複数のセンサにより受信した複数の信号を含むデータを準備する工程をさらに有し、
前記仮の測定画像を再構成する工程では、
前記複数の信号のうちの一部に基づいて前記仮の測定画像を再構成する
付記１に記載の画像生成方法。
（付記６）
前記仮の測定画像を再構成する工程では、
ノイズ除去のためのフィルタリングを実空間のフィルタで行うことにより、前記仮の測定画像を再構成する
付記１に記載の画像生成方法。
前記仮の測定画像を再構成する工程では、
ノイズ除去のためのフィルタリングを実空間のフィルタで行うことにより、前記仮の測定画像を再構成する
付記１に記載の画像生成方法。
（付記７）
前記関心領域の指定を受け付ける工程では、
複数の関心領域の指定を受け付ける
付記１～６のいずれか１つに記載の画像生成方法。
前記関心領域の指定を受け付ける工程では、
複数の関心領域の指定を受け付ける
付記１～６のいずれか１つに記載の画像生成方法。
（付記８）
前記関心領域画像を再構成する工程では、
前記複数の関心領域において、共通に設定された音速パラメータに基づいて、複数の関心領域画像を再構成し、
前記音速パラメータの前記適切な値を決定する工程では、
前記複数の関心領域において、共通する音速パラメータの適切な値を決定する
付記７に記載の画像生成方法。
前記関心領域画像を再構成する工程では、
前記複数の関心領域において、共通に設定された音速パラメータに基づいて、複数の関心領域画像を再構成し、
前記音速パラメータの前記適切な値を決定する工程では、
前記複数の関心領域において、共通する音速パラメータの適切な値を決定する
付記７に記載の画像生成方法。
（付記９）
所定の被検部位を測定した音響波の信号データに基づいて、音響画像を生成する画像生成方法であって、
前記被検部位の撮影領域内で部分的な関心領域の指定を受け付ける工程と、
前記関心領域において、音速パラメータを変更しながら、当該変更された音速パラメータに基づいて、音響画像からなる関心領域画像を順次再構成する工程と、
前記関心領域画像に基づいて前記音速パラメータの適切な値を決定する工程と、
前記音速パラメータの前記適切な値に基づいて、前記関心領域よりも広い領域の音響画像を再構成する工程と、
を有し、
前記関心領域の指定を受け付ける工程では、
複数の関心領域の指定を受け付け、
前記関心領域画像を再構成する工程では、
前記複数の関心領域において、共通に設定された音速パラメータに基づいて、複数の関心領域画像を再構成し、
前記音速パラメータの前記適切な値を決定する工程では、
前記複数の関心領域において、共通する音速パラメータの適切な値を決定する
画像生成方法。
所定の被検部位を測定した音響波の信号データに基づいて、音響画像を生成する画像生成方法であって、
前記被検部位の撮影領域内で部分的な関心領域の指定を受け付ける工程と、
前記関心領域において、音速パラメータを変更しながら、当該変更された音速パラメータに基づいて、音響画像からなる関心領域画像を順次再構成する工程と、
前記関心領域画像に基づいて前記音速パラメータの適切な値を決定する工程と、
前記音速パラメータの前記適切な値に基づいて、前記関心領域よりも広い領域の音響画像を再構成する工程と、
を有し、
前記関心領域の指定を受け付ける工程では、
複数の関心領域の指定を受け付け、
前記関心領域画像を再構成する工程では、
前記複数の関心領域において、共通に設定された音速パラメータに基づいて、複数の関心領域画像を再構成し、
前記音速パラメータの前記適切な値を決定する工程では、
前記複数の関心領域において、共通する音速パラメータの適切な値を決定する
画像生成方法。
（付記１０）
前記関心領域画像を再構成する工程では、
前記複数の関心領域において、それぞれ個別に設定された音速パラメータに基づいて、複数の関心領域画像を再構成する
付記７に記載の画像生成方法。
前記関心領域画像を再構成する工程では、
前記複数の関心領域において、それぞれ個別に設定された音速パラメータに基づいて、複数の関心領域画像を再構成する
付記７に記載の画像生成方法。
（付記１１）
所定の被検部位を測定した音響波の信号データに基づいて、音響画像を生成する画像生成方法であって、
前記被検部位の撮影領域内で部分的な関心領域の指定を受け付ける工程と、
前記関心領域において、音速パラメータを変更しながら、当該変更された音速パラメータに基づいて、音響画像からなる関心領域画像を順次再構成する工程と、
前記関心領域画像に基づいて前記音速パラメータの適切な値を決定する工程と、
前記音速パラメータの前記適切な値に基づいて、前記関心領域よりも広い領域の音響画像を再構成する工程と、
を有し、
前記関心領域の指定を受け付ける工程では、
複数の関心領域の指定を受け付け、
前記関心領域画像を再構成する工程では、
前記複数の関心領域において、それぞれ個別に設定された音速パラメータに基づいて、複数の関心領域画像を再構成する
画像生成方法。
所定の被検部位を測定した音響波の信号データに基づいて、音響画像を生成する画像生成方法であって、
前記被検部位の撮影領域内で部分的な関心領域の指定を受け付ける工程と、
前記関心領域において、音速パラメータを変更しながら、当該変更された音速パラメータに基づいて、音響画像からなる関心領域画像を順次再構成する工程と、
前記関心領域画像に基づいて前記音速パラメータの適切な値を決定する工程と、
前記音速パラメータの前記適切な値に基づいて、前記関心領域よりも広い領域の音響画像を再構成する工程と、
を有し、
前記関心領域の指定を受け付ける工程では、
複数の関心領域の指定を受け付け、
前記関心領域画像を再構成する工程では、
前記複数の関心領域において、それぞれ個別に設定された音速パラメータに基づいて、複数の関心領域画像を再構成する
画像生成方法。
（付記１２）
前記信号データに基づいて、音響画像からなる仮の測定画像を再構成する工程をさらに有し、
前記関心領域の指定を受け付ける工程では、
前記仮の測定画像内で前記関心領域の指定を受け付ける
付記９又は１１に記載の画像生成方法。
前記信号データに基づいて、音響画像からなる仮の測定画像を再構成する工程をさらに有し、
前記関心領域の指定を受け付ける工程では、
前記仮の測定画像内で前記関心領域の指定を受け付ける
付記９又は１１に記載の画像生成方法。
（付記１３）
前記関心領域画像を再構成する工程では、
実空間での位置関係と等しい位置関係で画像再構成したときの前記複数の関心領域同士の間隔よりも近づけて、前記複数の関心領域画像を所定の画面上に表示させる
付記７～１２のいずれか１つに記載の画像生成方法。
前記関心領域画像を再構成する工程では、
実空間での位置関係と等しい位置関係で画像再構成したときの前記複数の関心領域同士の間隔よりも近づけて、前記複数の関心領域画像を所定の画面上に表示させる
付記７～１２のいずれか１つに記載の画像生成方法。
（付記１４）
前記関心領域画像を再構成する工程では、
前記音速パラメータの任意の変更を受け付けながら、当該変更された音速パラメータに基づいて、前記関心領域画像を順次再構成する
付記１～１３のいずれか１つに記載の画像生成方法。
前記関心領域画像を再構成する工程では、
前記音速パラメータの任意の変更を受け付けながら、当該変更された音速パラメータに基づいて、前記関心領域画像を順次再構成する
付記１～１３のいずれか１つに記載の画像生成方法。
（付記１５）
前記関心領域画像を再構成する工程は、
前記音速パラメータの変更範囲および変更幅のうち少なくともいずれかの変更条件を受け付ける工程と、
前記変更条件下で前記音速パラメータを自動で変更しながら、当該変更された音速パラメータに基づいて、前記関心領域画像を順次再構成する工程と、
を有する
付記１～１３のいずれか１つに記載の画像生成方法。
前記関心領域画像を再構成する工程は、
前記音速パラメータの変更範囲および変更幅のうち少なくともいずれかの変更条件を受け付ける工程と、
前記変更条件下で前記音速パラメータを自動で変更しながら、当該変更された音速パラメータに基づいて、前記関心領域画像を順次再構成する工程と、
を有する
付記１～１３のいずれか１つに記載の画像生成方法。
（付記１６）
前記関心領域画像を再構成する工程では、
予め設定された変更範囲および変更幅で前記音速パラメータを自動で変更しながら、当該変更された音速パラメータに基づいて、前記関心領域画像を順次再構成する
付記１～１３のいずれか１つに記載の画像生成方法。
前記関心領域画像を再構成する工程では、
予め設定された変更範囲および変更幅で前記音速パラメータを自動で変更しながら、当該変更された音速パラメータに基づいて、前記関心領域画像を順次再構成する
付記１～１３のいずれか１つに記載の画像生成方法。
（付記１７）
前記音速パラメータの前記適切な値を決定する工程では、
自動で変更される前記音速パラメータを停止させたときの当該音速パラメータを、前記適切な値として決定する
付記１５又は１６に記載の画像生成方法。
前記音速パラメータの前記適切な値を決定する工程では、
自動で変更される前記音速パラメータを停止させたときの当該音速パラメータを、前記適切な値として決定する
付記１５又は１６に記載の画像生成方法。
（付記１８）
前記音速パラメータの前記適切な値を決定する工程では、
前記関心領域画像の解像度を除く画質に基づいて、前記音速パラメータの前記適切な値を自動で決定する
付記１５又は１６に記載の画像生成方法。
前記音速パラメータの前記適切な値を決定する工程では、
前記関心領域画像の解像度を除く画質に基づいて、前記音速パラメータの前記適切な値を自動で決定する
付記１５又は１６に記載の画像生成方法。
（付記１９）
前記関心領域内で前記被検部位の注目部位の指定を受け付ける工程をさらに有し、
前記音速パラメータの前記適切な値を決定する工程では、
前記注目部位の解像度を除く画質に基づいて、前記音速パラメータの前記適切な値を自動で決定する
付記１８に記載の画像生成方法。
前記関心領域内で前記被検部位の注目部位の指定を受け付ける工程をさらに有し、
前記音速パラメータの前記適切な値を決定する工程では、
前記注目部位の解像度を除く画質に基づいて、前記音速パラメータの前記適切な値を自動で決定する
付記１８に記載の画像生成方法。
（付記２０）
所定の被検部位を測定した音響波の信号データに基づいて、音響画像を生成する画像生成プログラムであって、
前記信号データに基づいて、音響画像からなる仮の測定画像を再構成する手順と、
前記仮の測定画像内で部分的な関心領域の指定を受け付ける手順と、
前記関心領域において、音速パラメータを変更しながら、当該変更された音速パラメータに基づいて、音響画像からなる関心領域画像を順次再構成する手順と、
前記関心領域画像に基づいて前記音速パラメータの適切な値を決定する手順と、
前記音速パラメータの前記適切な値に基づいて、前記関心領域よりも広い領域の音響画像を再構成する手順と、
をコンピュータに実行させ、
前記仮の測定画像を再構成する手順では、
前記仮の測定画像を再構成するための単位面積当たりの計算量を、前記関心領域画像を再構成するための単位面積当たりの計算量よりも少なくする
画像生成プログラム、
または、該プログラムを記録したコンピュータ読み取り可能な記録媒体。
所定の被検部位を測定した音響波の信号データに基づいて、音響画像を生成する画像生成プログラムであって、
前記信号データに基づいて、音響画像からなる仮の測定画像を再構成する手順と、
前記仮の測定画像内で部分的な関心領域の指定を受け付ける手順と、
前記関心領域において、音速パラメータを変更しながら、当該変更された音速パラメータに基づいて、音響画像からなる関心領域画像を順次再構成する手順と、
前記関心領域画像に基づいて前記音速パラメータの適切な値を決定する手順と、
前記音速パラメータの前記適切な値に基づいて、前記関心領域よりも広い領域の音響画像を再構成する手順と、
をコンピュータに実行させ、
前記仮の測定画像を再構成する手順では、
前記仮の測定画像を再構成するための単位面積当たりの計算量を、前記関心領域画像を再構成するための単位面積当たりの計算量よりも少なくする
画像生成プログラム、
または、該プログラムを記録したコンピュータ読み取り可能な記録媒体。
（付記２１）
所定の被検部位を測定した音響波の信号データに基づいて、音響画像を生成する画像生成プログラムであって、
前記被検部位の撮影領域内で部分的な関心領域の指定を受け付ける手順と、
前記関心領域において、音速パラメータを変更しながら、当該変更された音速パラメータに基づいて、音響画像からなる関心領域画像を順次再構成する手順と、
前記関心領域画像に基づいて前記音速パラメータの適切な値を決定する手順と、
前記音速パラメータの前記適切な値に基づいて、前記関心領域よりも広い領域の音響画像を再構成する手順と、
をコンピュータに実行させ、
前記関心領域の指定を受け付ける手順では、
複数の関心領域の指定を受け付け、
前記関心領域画像を再構成する手順では、
前記複数の関心領域において、共通に設定された音速パラメータに基づいて、複数の関心領域画像を再構成し、
前記音速パラメータの前記適切な値を決定する手順では、
前記複数の関心領域において、共通する音速パラメータの適切な値を決定する
画像生成プログラム、
または、該プログラムを記録したコンピュータ読み取り可能な記録媒体。
所定の被検部位を測定した音響波の信号データに基づいて、音響画像を生成する画像生成プログラムであって、
前記被検部位の撮影領域内で部分的な関心領域の指定を受け付ける手順と、
前記関心領域において、音速パラメータを変更しながら、当該変更された音速パラメータに基づいて、音響画像からなる関心領域画像を順次再構成する手順と、
前記関心領域画像に基づいて前記音速パラメータの適切な値を決定する手順と、
前記音速パラメータの前記適切な値に基づいて、前記関心領域よりも広い領域の音響画像を再構成する手順と、
をコンピュータに実行させ、
前記関心領域の指定を受け付ける手順では、
複数の関心領域の指定を受け付け、
前記関心領域画像を再構成する手順では、
前記複数の関心領域において、共通に設定された音速パラメータに基づいて、複数の関心領域画像を再構成し、
前記音速パラメータの前記適切な値を決定する手順では、
前記複数の関心領域において、共通する音速パラメータの適切な値を決定する
画像生成プログラム、
または、該プログラムを記録したコンピュータ読み取り可能な記録媒体。
（付記２２）
所定の被検部位を測定した音響波の信号データに基づいて、音響画像を生成する画像生成プログラムであって、
前記被検部位の撮影領域内で部分的な関心領域の指定を受け付ける手順と、
前記関心領域において、音速パラメータを変更しながら、当該変更された音速パラメータに基づいて、音響画像からなる関心領域画像を順次再構成する手順と、
前記関心領域画像に基づいて前記音速パラメータの適切な値を決定する手順と、
前記音速パラメータの前記適切な値に基づいて、前記関心領域よりも広い領域の音響画像を再構成する手順と、
をコンピュータに実行させ、
前記関心領域の指定を受け付ける手順では、
複数の関心領域の指定を受け付け、
前記関心領域画像を再構成する手順では、
前記複数の関心領域において、それぞれ個別に設定された音速パラメータに基づいて、複数の関心領域画像を再構成する
画像生成プログラム、
または、該プログラムを記録したコンピュータ読み取り可能な記録媒体。
所定の被検部位を測定した音響波の信号データに基づいて、音響画像を生成する画像生成プログラムであって、
前記被検部位の撮影領域内で部分的な関心領域の指定を受け付ける手順と、
前記関心領域において、音速パラメータを変更しながら、当該変更された音速パラメータに基づいて、音響画像からなる関心領域画像を順次再構成する手順と、
前記関心領域画像に基づいて前記音速パラメータの適切な値を決定する手順と、
前記音速パラメータの前記適切な値に基づいて、前記関心領域よりも広い領域の音響画像を再構成する手順と、
をコンピュータに実行させ、
前記関心領域の指定を受け付ける手順では、
複数の関心領域の指定を受け付け、
前記関心領域画像を再構成する手順では、
前記複数の関心領域において、それぞれ個別に設定された音速パラメータに基づいて、複数の関心領域画像を再構成する
画像生成プログラム、
または、該プログラムを記録したコンピュータ読み取り可能な記録媒体。
（付記２３）
所定の被検部位を測定した音響波の信号データに基づいて、音響画像を生成する処理部を備え、
前記処理部は、
前記信号データに基づいて、音響画像からなる仮の測定画像を再構成する処理と、
前記仮の測定画像内で部分的な関心領域の指定を受け付ける処理と、
前記関心領域において、音速パラメータを変更しながら、当該変更された音速パラメータに基づいて、音響画像からなる関心領域画像を順次再構成する処理と、
前記関心領域画像に基づいて前記音速パラメータの適切な値を決定する処理と、
前記音速パラメータの前記適切な値に基づいて、前記関心領域よりも広い領域の音響画像を再構成する処理と、
を実行し、
前記仮の測定画像を再構成する処理では、
前記仮の測定画像を再構成するための単位面積当たりの計算量を、前記関心領域画像を再構成するための単位面積当たりの計算量よりも少なくする
画像生成装置。
所定の被検部位を測定した音響波の信号データに基づいて、音響画像を生成する処理部を備え、
前記処理部は、
前記信号データに基づいて、音響画像からなる仮の測定画像を再構成する処理と、
前記仮の測定画像内で部分的な関心領域の指定を受け付ける処理と、
前記関心領域において、音速パラメータを変更しながら、当該変更された音速パラメータに基づいて、音響画像からなる関心領域画像を順次再構成する処理と、
前記関心領域画像に基づいて前記音速パラメータの適切な値を決定する処理と、
前記音速パラメータの前記適切な値に基づいて、前記関心領域よりも広い領域の音響画像を再構成する処理と、
を実行し、
前記仮の測定画像を再構成する処理では、
前記仮の測定画像を再構成するための単位面積当たりの計算量を、前記関心領域画像を再構成するための単位面積当たりの計算量よりも少なくする
画像生成装置。
（付記２４）
所定の被検部位を測定した音響波の信号データに基づいて、音響画像を生成する処理部を備え、
前記処理部は、
前記被検部位の撮影領域内で部分的な関心領域の指定を受け付ける処理と、
前記関心領域において、音速パラメータを変更しながら、当該変更された音速パラメータに基づいて、音響画像からなる関心領域画像を順次再構成する処理と、
前記関心領域画像に基づいて前記音速パラメータの適切な値を決定する処理と、
前記音速パラメータの前記適切な値に基づいて、前記関心領域よりも広い領域の音響画像を再構成する処理と、
を実行し、
前記関心領域の指定を受け付ける処理では、
複数の関心領域の指定を受け付け、
前記関心領域画像を再構成する処理では、
前記複数の関心領域において、共通に設定された音速パラメータに基づいて、複数の関心領域画像を再構成し、
前記音速パラメータの前記適切な値を決定する処理では、
前記複数の関心領域において、共通する音速パラメータの適切な値を決定する
画像生成装置。
所定の被検部位を測定した音響波の信号データに基づいて、音響画像を生成する処理部を備え、
前記処理部は、
前記被検部位の撮影領域内で部分的な関心領域の指定を受け付ける処理と、
前記関心領域において、音速パラメータを変更しながら、当該変更された音速パラメータに基づいて、音響画像からなる関心領域画像を順次再構成する処理と、
前記関心領域画像に基づいて前記音速パラメータの適切な値を決定する処理と、
前記音速パラメータの前記適切な値に基づいて、前記関心領域よりも広い領域の音響画像を再構成する処理と、
を実行し、
前記関心領域の指定を受け付ける処理では、
複数の関心領域の指定を受け付け、
前記関心領域画像を再構成する処理では、
前記複数の関心領域において、共通に設定された音速パラメータに基づいて、複数の関心領域画像を再構成し、
前記音速パラメータの前記適切な値を決定する処理では、
前記複数の関心領域において、共通する音速パラメータの適切な値を決定する
画像生成装置。
（付記２５）
所定の被検部位を測定した音響波の信号データに基づいて、音響画像を生成する処理部を備え、
前記処理部は、
前記被検部位の撮影領域内で部分的な関心領域の指定を受け付ける処理と、
前記関心領域において、音速パラメータを変更しながら、当該変更された音速パラメータに基づいて、音響画像からなる関心領域画像を順次再構成する処理と、
前記関心領域画像に基づいて前記音速パラメータの適切な値を決定する手順と、
前記音速パラメータの前記適切な値に基づいて、前記関心領域よりも広い領域の音響画像を再構成する処理と、
を実行し、
前記関心領域の指定を受け付ける処理では、
複数の関心領域の指定を受け付け、
前記関心領域画像を再構成する処理では、
前記複数の関心領域において、それぞれ個別に設定された音速パラメータに基づいて、複数の関心領域画像を再構成する
画像生成装置。
所定の被検部位を測定した音響波の信号データに基づいて、音響画像を生成する処理部を備え、
前記処理部は、
前記被検部位の撮影領域内で部分的な関心領域の指定を受け付ける処理と、
前記関心領域において、音速パラメータを変更しながら、当該変更された音速パラメータに基づいて、音響画像からなる関心領域画像を順次再構成する処理と、
前記関心領域画像に基づいて前記音速パラメータの適切な値を決定する手順と、
前記音速パラメータの前記適切な値に基づいて、前記関心領域よりも広い領域の音響画像を再構成する処理と、
を実行し、
前記関心領域の指定を受け付ける処理では、
複数の関心領域の指定を受け付け、
前記関心領域画像を再構成する処理では、
前記複数の関心領域において、それぞれ個別に設定された音速パラメータに基づいて、複数の関心領域画像を再構成する
画像生成装置。
１０ 光音響イメージング装置
２０ 測定画像
２２ 画素
２４ 関心領域
４０ 関心領域画像
５０ 広域画像
１００ 被検体
１１０ 被検部位
１２０ 音源
２００ 支持台
２１０ 支持面
２２０ 開口
３００ センサユニット
３１０ 音響整合材
３２０ 容器
３４０（３４０ｉ） センサ
３６０ 素子保持部
３８０ 走査機構
４００ 分離部
４１０ 音響整合材
４２０ 分離フィルム
６２０ 光源
６４０ 光学系
６６０ 光出射口
７００ 処理部
７０１ ＣＰＵ
７０３ ＲＡＭ
７０４ 記憶装置
７０５ Ｉ／Ｏポート
７０９ コンピュータ
７１０ 信号処理部
７２０ 表示部
７４０ 入力部
２０ 測定画像
２２ 画素
２４ 関心領域
４０ 関心領域画像
５０ 広域画像
１００ 被検体
１１０ 被検部位
１２０ 音源
２００ 支持台
２１０ 支持面
２２０ 開口
３００ センサユニット
３１０ 音響整合材
３２０ 容器
３４０（３４０ｉ） センサ
３６０ 素子保持部
３８０ 走査機構
４００ 分離部
４１０ 音響整合材
４２０ 分離フィルム
６２０ 光源
６４０ 光学系
６６０ 光出射口
７００ 処理部
７０１ ＣＰＵ
７０３ ＲＡＭ
７０４ 記憶装置
７０５ Ｉ／Ｏポート
７０９ コンピュータ
７１０ 信号処理部
７２０ 表示部
７４０ 入力部
Claims (17)
- 所定の被検部位を測定した音響波の信号データに基づいて、音響画像を生成する画像生成方法であって、
前記信号データに基づいて、音響画像からなる仮の測定画像を再構成する工程と、
前記仮の測定画像内で部分的な関心領域の指定を受け付ける工程と、
前記関心領域において、音速パラメータを変更しながら、当該変更された音速パラメータに基づいて、音響画像からなる関心領域画像を順次再構成する工程と、
前記関心領域画像に基づいて前記音速パラメータの適切な値を決定する工程と、
前記音速パラメータの前記適切な値に基づいて、前記関心領域よりも広い領域の音響画像を再構成する工程と、
を有し、
前記仮の測定画像を再構成する工程では、
前記仮の測定画像の画素サイズを前記関心領域画像の画素サイズよりも大きくすることにより、
前記仮の測定画像を再構成するための単位面積当たりの計算量を、前記関心領域画像を再構成するための単位面積当たりの計算量よりも少なくする
画像生成方法。 - 前記関心領域の指定を受け付ける工程では、
複数の関心領域の指定を受け付ける
請求項１に記載の画像生成方法。 - 前記関心領域画像を再構成する工程では、
前記複数の関心領域において、共通に設定された音速パラメータに基づいて、複数の関心領域画像を再構成し、
前記音速パラメータの前記適切な値を決定する工程では、
前記複数の関心領域において、共通する音速パラメータの適切な値を決定する
請求項２に記載の画像生成方法。 - 前記関心領域画像を再構成する工程では、
前記複数の関心領域において、それぞれ個別に設定された音速パラメータに基づいて、複数の関心領域画像を再構成する
請求項２に記載の画像生成方法。 - 前記関心領域画像を再構成する工程では、
実空間での位置関係と等しい位置関係で画像再構成したときの前記複数の関心領域同士の間隔よりも近づけて、前記複数の関心領域画像を所定の画面上に表示させる
請求項２～４のいずれか１項に記載の画像生成方法。 - 前記関心領域画像を再構成する工程では、
前記音速パラメータの任意の変更を受け付けながら、当該変更された音速パラメータに基づいて、前記関心領域画像を順次再構成する
請求項１～５のいずれか１項に記載の画像生成方法。 - 前記関心領域画像を再構成する工程は、
前記音速パラメータの変更範囲および変更幅のうち少なくともいずれかの変更条件を受け付ける工程と、
前記変更条件下で前記音速パラメータを自動で変更しながら、当該変更された音速パラメータに基づいて、前記関心領域画像を順次再構成する工程と、
を有する
請求項１～５のいずれか１項に記載の画像生成方法。 - 前記関心領域画像を再構成する工程では、
予め設定された変更範囲および変更幅で前記音速パラメータを自動で変更しながら、当該変更された音速パラメータに基づいて、前記関心領域画像を順次再構成する
請求項１～５のいずれか１項に記載の画像生成方法。 - 前記関心領域画像を再構成する工程では、
前記仮の測定画像を再構成する工程で用いられた前記音響波の第１信号データの一部を共通データとして含む前記音響波の第２信号データを用いる
請求項１～８のいずれか１項に記載の画像生成方法。 - 所定の被検部位を測定した音響波の信号データに基づいて、音響画像を生成する画像生成プログラムであって、
前記信号データに基づいて、音響画像からなる仮の測定画像を再構成する手順と、
前記仮の測定画像内で部分的な関心領域の指定を受け付ける手順と、
前記関心領域において、音速パラメータを変更しながら、当該変更された音速パラメータに基づいて、音響画像からなる関心領域画像を順次再構成する手順と、
前記関心領域画像に基づいて前記音速パラメータの適切な値を決定する手順と、
前記音速パラメータの前記適切な値に基づいて、前記関心領域よりも広い領域の音響画像を再構成する手順と、
をコンピュータに実行させ、
前記仮の測定画像を再構成する手順では、
前記仮の測定画像の画素サイズを前記関心領域画像の画素サイズよりも大きくすることにより、
前記仮の測定画像を再構成するための単位面積当たりの計算量を、前記関心領域画像を再構成するための単位面積当たりの計算量よりも少なくする
画像生成プログラム。 - 前記関心領域画像を再構成する手順では、
前記仮の測定画像を再構成する工程で用いられた前記音響波の第１信号データの一部を共通データとして含む前記音響波の第２信号データを用いる
請求項１０に記載の画像生成プログラム。 - 前記関心領域の指定を受け付ける手順では、
複数の関心領域の指定を受け付ける
請求項１０または１１に記載の画像生成プログラム。 - 前記関心領域画像を再構成する手順では、
実空間での位置関係と等しい位置関係で画像再構成したときの前記複数の関心領域同士の間隔よりも近づけて、前記複数の関心領域画像を所定の画面上に表示させる
請求項１２に記載の画像生成プログラム。 - 所定の被検部位を測定した音響波の信号データに基づいて、音響画像を生成する処理部を備え、
前記処理部は、
前記信号データに基づいて、音響画像からなる仮の測定画像を再構成する処理と、
前記仮の測定画像内で部分的な関心領域の指定を受け付ける処理と、
前記関心領域において、音速パラメータを変更しながら、当該変更された音速パラメータに基づいて、音響画像からなる関心領域画像を順次再構成する処理と、
前記関心領域画像に基づいて前記音速パラメータの適切な値を決定する処理と、
前記音速パラメータの前記適切な値に基づいて、前記関心領域よりも広い領域の音響画像を再構成する処理と、
を実行し、
前記仮の測定画像を再構成する処理では、
前記仮の測定画像の画素サイズを前記関心領域画像の画素サイズよりも大きくすることにより、
前記仮の測定画像を再構成するための単位面積当たりの計算量を、前記関心領域画像を再構成するための単位面積当たりの計算量よりも少なくする
画像生成装置。 - 前記関心領域画像を再構成する処理では、
前記仮の測定画像を再構成する工程で用いられた前記音響波の第１信号データの一部を共通データとして含む前記音響波の第２信号データを用いる
請求項１４に記載の画像生成装置。 - 前記関心領域の指定を受け付ける処理では、
複数の関心領域の指定を受け付ける
請求項１４または１５に記載の画像生成装置。 - 前記関心領域画像を再構成する処理では、
実空間での位置関係と等しい位置関係で画像再構成したときの前記複数の関心領域同士の間隔よりも近づけて、前記複数の関心領域画像を所定の画面上に表示させる
請求項１６に記載の画像生成装置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2020138622A JP7479684B2 (ja) | 2020-08-19 | 2020-08-19 | 画像生成方法、画像生成プログラムおよび画像生成装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2020138622A JP7479684B2 (ja) | 2020-08-19 | 2020-08-19 | 画像生成方法、画像生成プログラムおよび画像生成装置 |
Publications (3)
| Publication Number | Publication Date |
|---|---|
| JP2022034766A JP2022034766A (ja) | 2022-03-04 |
| JP2022034766A5 JP2022034766A5 (ja) | 2023-08-01 |
| JP7479684B2 true JP7479684B2 (ja) | 2024-05-09 |
Family
ID=80443840
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2020138622A Active JP7479684B2 (ja) | 2020-08-19 | 2020-08-19 | 画像生成方法、画像生成プログラムおよび画像生成装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP7479684B2 (ja) |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2011120765A (ja) | 2009-12-11 | 2011-06-23 | Canon Inc | 画像生成装置、画像生成方法、及び、プログラム |
| JP2012157387A (ja) | 2011-01-28 | 2012-08-23 | Toshiba Corp | 超音波診断装置及び画像生成制御プログラム |
| JP2013244204A (ja) | 2012-05-25 | 2013-12-09 | Fujifilm Corp | 超音波画像診断装置 |
| JP2014068755A (ja) | 2012-09-28 | 2014-04-21 | Fujifilm Corp | 超音波検査装置、超音波検査装置の信号処理方法およびプログラム |
| JP2016101415A (ja) | 2014-11-28 | 2016-06-02 | キヤノン株式会社 | 被検体情報取得装置 |
| WO2019044075A1 (ja) | 2017-08-31 | 2019-03-07 | 富士フイルム株式会社 | 画像生成装置および作動方法 |
-
2020
- 2020-08-19 JP JP2020138622A patent/JP7479684B2/ja active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2011120765A (ja) | 2009-12-11 | 2011-06-23 | Canon Inc | 画像生成装置、画像生成方法、及び、プログラム |
| US20140148683A1 (en) | 2009-12-11 | 2014-05-29 | Canon Kabushiki Kaisha | Image generating apparatus, image generating method, and program |
| JP2012157387A (ja) | 2011-01-28 | 2012-08-23 | Toshiba Corp | 超音波診断装置及び画像生成制御プログラム |
| JP2013244204A (ja) | 2012-05-25 | 2013-12-09 | Fujifilm Corp | 超音波画像診断装置 |
| JP2014068755A (ja) | 2012-09-28 | 2014-04-21 | Fujifilm Corp | 超音波検査装置、超音波検査装置の信号処理方法およびプログラム |
| JP2016101415A (ja) | 2014-11-28 | 2016-06-02 | キヤノン株式会社 | 被検体情報取得装置 |
| WO2019044075A1 (ja) | 2017-08-31 | 2019-03-07 | 富士フイルム株式会社 | 画像生成装置および作動方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| JP2022034766A (ja) | 2022-03-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP5528083B2 (ja) | 画像生成装置、画像生成方法、及び、プログラム | |
| JP6174656B2 (ja) | 対象体の立体実時間光音響撮像用の手持ち式装置及び方法 | |
| US20140196544A1 (en) | Object information acquiring apparatus | |
| JP5917037B2 (ja) | 被検体情報取得装置および被検体情報取得方法 | |
| JP6544910B2 (ja) | 情報処理装置、被検体情報取得装置及び音速決定方法 | |
| WO2012011242A1 (en) | Image information acquiring apparatus, image information acquiring method and image information acquiring program | |
| JP5984547B2 (ja) | 被検体情報取得装置およびその制御方法 | |
| JP5847454B2 (ja) | 被検体情報取得装置、表示制御方法およびプログラム | |
| KR102218308B1 (ko) | 초음파 영상 처리 장치 및 방법 | |
| JP6742734B2 (ja) | 被検体情報取得装置および信号処理方法 | |
| JP6656229B2 (ja) | 光音響装置 | |
| CN108601536A (zh) | 信息获取装置和控制方法 | |
| CN107847141B (zh) | 用于获取与多个图像数据集的位置位移有关的信息的设备、方法和程序 | |
| JP7479684B2 (ja) | 画像生成方法、画像生成プログラムおよび画像生成装置 | |
| JP6234518B2 (ja) | 情報処理装置および情報処理方法 | |
| WO2023047601A1 (ja) | 画像生成方法、画像生成プログラムおよび画像生成装置 | |
| JP6072206B2 (ja) | 被検体情報取得装置および表示方法 | |
| KR20180062393A (ko) | 표시 제어장치, 표시 제어방법, 및 기억매체 | |
| JP6486056B2 (ja) | 光音響装置および光音響装置の処理方法 | |
| JP5940109B2 (ja) | 画像生成装置、伝播速度決定方法、及び、プログラム | |
| Mustafa et al. | In vivo three-dimensional raster scan optoacoustic mesoscopy using frequency domain inversion | |
| JP7077384B2 (ja) | 被検体情報取得装置 | |
| JP2018192309A (ja) | 情報処理装置、情報処理方法およびプログラム | |
| JP6942847B2 (ja) | 被検体情報取得装置および信号処理方法 | |
| JP2018012027A (ja) | 記録データのデータ構造 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20230724 |
|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20230724 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20240216 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20240319 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20240417 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 7479684 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |