(57)【要約】【課題】スペックルの除去と組織構造の明確化という二つの要請を同時に満たすことが可能な技術を提供する。【解決手段】超音波画像のノイズを除去し、ノイズ除去後の画像に対し、モルフォロジー処理を施す。モルフォロジー処理は、ダイレーション／エロージョンを行う第一演算と、オープニング／クロージングを行う第二演算と、ノイズ除去後の画像に対し第一演算を行った結果を用いて、モルフォロジーの第二演算に用いる構造要素の値を決定する。【選択図】図１
2020110296.tif 000002
【特許請求の範囲】
【請求項１】
検査対象に超音波を送信するとともに反射波である超音波を受信し、超音波信号を取得する超音波計測部と、
前記超音波計測部が取得した超音波信号を用いて、検査対象の超音波画像を生成する画像生成部と、
前記画像生成部が作成した超音波画像を処理する画像処理部と、を備え、
前記画像処理部は、前記超音波画像のノイズを除去するノイズ除去部と、前記ノイズ除去部が生成するノイズ除去後画像についてモルフォロジー処理を行うモルフォロジー処理部と、前記ノイズ除去後画像を用いて前記モルフォロジー処理に用いる構造要素の値を決定する構造要素決定部と、を備えることを特徴とする超音波撮像装置。
【請求項２】
請求項１に記載の超音波撮像装置であって、
前記モルフォロジー処理部は、第一の構造要素を用いたダイレーション及びエロージョンの少なくとも一方を行う第一演算部と、第二の構造要素を用いたオープニング及びクロージングの少なくとも一方を行う第二演算部とを有し、
前記構造要素決定部は、前記ノイズ除去後画像を用いて、前記第二の構造要素の値を決定することを特徴とする超音波撮像装置。
【請求項３】
請求項２に記載の超音波撮像装置であって、
前記構造要素決定部は、前記第一演算部でダイレーション処理された画像とエロージョン処理された画像との差分により、前記第二の構造要素の値を決定することを特徴とする超音波撮像装置。
【請求項４】
請求項２に記載の超音波撮像装置であって、
前記第一演算部は、第一の構造要素を用いたダイレーション及びエロージョンを行い、
前記第二演算部は、前記第一の構造要素を用いたダイレーション結果に対し、前記第二の構造要素を用いてクロージングを行い、前記第一の構造要素を用いたエロージョン結果に対しオープニングを行い、
前記モルフォロジー処理部は、前記第二の構造要素を用いたクロージング結果と前記第一の構造要素を用いたクロージング結果との差分、及び、前記第二の構造要素を用いたオープニング結果と前記第一の構造要素を用いたオープニング結果との差分を、前記ノイズ除去後画像に加算し、モルフォロジー処理画像を生成することを特徴とする超音波撮像装置。
【請求項５】
請求項１に記載の超音波撮像装置であって、
前記ノイズ除去部は、エッジ保存型平滑化フィルタを用いてノイズ除去を行うことを特徴とする超音波撮像装置。
【請求項６】
請求項１に記載の超音波撮像装置であって、
前記超音波画像の座標を極座標から直交座標に変換する座標変換部をさらに備え、
前記座標変換部は、前記モルフォロジー処理部による処理前の超音波画像に対し座標変換を行うことを特徴とする超音波撮像装置。
【請求項７】
請求項１に記載の超音波撮像装置であって、
前記超音波画像の座標を極座標から直交座標に変換する座標変換部をさらに備え、
前記座標変換部は、前記モルフォロジー処理部による処理後の超音波画像に対し座標変換を行うことを特徴とする超音波撮像装置。
【請求項８】
請求項１に記載の超音波撮像装置であって、
前記画像処理部は、前記構造要素のサイズを決定するサイズ決定部をさらに備えることを徴とする超音波撮像装置。
【請求項９】
請求項８に記載の超音波撮像装置であって、
前記サイズ決定部は、超音波画像の分解能及びサンプリング間隔をもとに前記構造要素のサイズを決定することを特徴とする超音波撮像装置。
【請求項１０】
請求項９に記載の超音波撮像装置であって、
前記画像処理部は、超音波計測の撮像条件と超音波画像の分解能との関係を示すテーブルを格納する記憶部をさらに備え、
前記サイズ決定部は、前記記憶部に格納された前記テーブルを参照して前記構造要素のサイズを決定することを特徴とする超音波撮像装置。
【請求項１１】
請求項８に記載の超音波撮像装置であって、
前記サイズ決定部は、超音波画像の画素位置ごとに前記構造要素のサイズを決定することを特徴とする超音波撮像装置。
【請求項１２】
請求項８に記載の超音波撮像装置であって、
前記サイズ決定部は、超音波画像の画素位置ごとに前記構造要素のサイズを算出し、その最大値、最小値、及び平均値のいずれかを、前記構造要素のサイズとして決定することを特徴とする超音波撮像装置。
【請求項１３】
請求項１に記載の超音波撮像装置であって、
前記モルフォロジー処理部による処理の選択を受付けるユーザーインターフェイス部をさらに備えたことを特徴とする超音波撮像装置。
【請求項１４】
検査対象の内部に穿針を挿入し、穿針の進路をガイドする支援画像を撮像する超音波撮像装置であって、検査対象に超音波を送信するとともに反射波である超音波を受信し、超音波信号を取得する超音波計測部と、前記超音波計測部が取得した超音波信号を用いて、検査対象の超音波画像を生成する画像生成部と、前記画像生成部が作成した超音波画像を処理する画像処理部と、を備え、
前記画像処理部は、前記画像生成部が作成した超音波画像のノイズを除去するノイズ除去部と、前記ノイズ除去部が生成するノイズ除去後画像についてモルフォロジー処理を行うモルフォロジー処理部、前記ノイズ除去後画像を用いて前記モルフォロジー処理部で用いる構造要素の値を決定する構造要素決定部と、を備え、
前記穿針とそれ以外の組織との境界を強調した前記支援画像を生成することを特徴とする超音波撮像装置。
【請求項１５】
超音波撮像装置で受信した超音波信号を用いて生成した超音波画像を処理する画像処理方法であって、
前記超音波画像のノイズを除去するステップと、
ノイズ除去後の画像に対し、モルフォロジー処理を施すステップと、を含み
前記ノイズ除去後の画像を用いて、前記モルフォロジー処理に用いる構造要素の値を決定するステップを含むことを特徴とする画像処理方法。
【請求項１６】
請求項１５に記載の超音波画像を処理する画像処理方法であって、
前記モルフォロジー処理は、第一の構造要素を用いたダイレーション及びエロージョンと、第二の構造要素を用いたオープニング及びクロージングの少なくとも一方とを含み、
前記第二の構造要素の値を、前記ダイレーション後の画像と前記エロージョン後の画像との差分により決定することを特徴とする画像処理方法。
【請求項１７】
請求項１５に記載の画像処理方法であって、
ノイズ除去前又はモルフォロジー処理前に、前記超音波画像の座標を極座標から直交座標に変換する座標変換ステップをさらに含むことを特徴とする画像処理方法。
【請求項１８】
請求項１５に記載の画像処理方法であって、
モルフォロジー処理後に、前記超音波画像の座標を極座標から直交座標に変換する座標変換ステップをさらに含むことを特徴とする画像処理方法。
【請求項１９】
請求項１５に記載の画像処理方法であって、
前記モルフォロジー処理で用いる構造要素のサイズを決定するステップをさらに含むことを特徴とする画像処理方法。
【請求項２０】
請求項１９に記載の画像処理方法であって、
前記構造要素のサイズを決定するステップは、超音波画像の分解能とサンプリング間隔とを用いて前記サイズを決定することを特徴とする画像処理方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は超音波撮像装置に係り、特に超音波撮像装置で得た超音波画像を対象とする画像処理に関する。
【背景技術】
【０００２】
超音波撮像装置は簡便に且つリアルタイムで検査対象の画像を取得し表示できる装置であり、医療分野等に広く用いられている。しかし、超音波撮像によって得られる画像は、ＣＴ装置等の医用撮像装置に比べ解像度が低く検査者の熟練度により検査対象の構造の特定や診断に差を生じやすい。また超音波画像の画質を劣化させる原因として、超音波の受信信号間の干渉により生じるスペックルと呼ばれるノイズがある。このスペックルは深度が深いほど細長くなるなど形状やサイズが画像内で異なることに加え、スペックルのサイズが生体の構造と同等となる場合があるため、通常のフィルタ処理等のノイズ除去技術では除去することが困難である。このためスペックルが画像内に残り、組織構造と見分けにくいという問題がある。
【０００３】
このような超音波画像の画質を向上する技術の一つとして、モルフォロジー技術を利用した技術が提案されている（非特許文献１、特許文献１）。モルフォロジーでは、ダイレーション（拡張処理）、エロージョン（減少処理）、オープニング処理及びクロージング処理の４つの処理を組み合わせて、原画像を処理することで、構造の境界であるエッジを強調することができる。オープニングは、エロージョン後の画像に対してダイレーションを行う処理であり、クロージングは、ダイレーション後の画像に対してエロージョンを行う処理である。これらの処理には、処理対象とする画像の画素値をその周囲の画素の画素値を用いて決定するために構造要素と呼ばれる関数が用いられる。構造要素のサイズは処理対象画素を中心として周囲の画素を含む領域のサイズであり、適切なサイズ及び値のものを原画像の各画素に適用して、ダイレーションやエロージョンを行うことで、所望のエッジ強調効果を得ることができる。
【０００４】
非特許文献１には、上述したモルフォロジー技術を超音波画像に適用する際に、ダイレーション／エロージョン後のオープニング／クロージングで用いる構造要素を超音波画像の位置に応じて異ならせる技術が開示されている。具体的には、処理対象である画像の位置が、エッジ領域か一様領域かに応じて構造要素の値を変化させ、これにより高いエッジ強調効果を図っている。構造要素の値の決定には、組織境界（エッジ領域）ではその他の領域に比べて輝度分散が大きくなる傾向にあることを利用し、構造要素内の輝度分散情報を活用している。
【０００５】
また特許文献１には、スペックル低減処理を行う際にアーチファクトを生じ難い画像処理を行うために、マスク領域内の位置に応じて設定される第１群の重み付け係数を乗じてダイレーションを行い、その結果に対して第２群の重み付け係数を乗じてエロージョンを行う方法が開示されている。この技術は入力画像に対して重み付け係数を乗じる手法であり、構造要素の値は画像内の位置によらず同一である。
【先行技術文献】
【非特許文献】
【０００６】
【非特許文献１】ＩＥＥＪ  Ｔｒａｎｓ．ＥＩＳ，椿井  正義ら「２重構造要素による多値モルフォロジーと超音波画像のコントラスト強調への応用」、Ｖｏｌ．１２４、Ｎｏ．９、２００４
【特許文献】
【０００７】
【特許文献１】特開２００５−２０５１９９号公報
【発明の開示】
【発明が解決しようとする課題】
【０００８】
非特許文献１に記載された技術では、オープニングやクロージングに用いる構造要素の値を画素の位置によって変化させることで、高いエッジ強調を図るとともに、平滑部においてはスペックルノイズを効果的に除去するというものであるが、構造とスペックルの区別が不十分な場合、スペックルのエッジが強調されてしまい、平滑化効果が限定的となってしまう場合がある。
【０００９】
特許文献１に記載された技術では、エッジ領域か一様領域かによらず同一の処理を行っているため、得られるエッジ強調効果は限定的である。
【００１０】
本発明は、ノイズ除去後画像に基づいて、構造要素の値を適応的に決定することにより、スペックルの除去と組織構造の明確化という二つの要請を同時に満たすことが可能な技術を提供することを課題とする。
【課題を解決するための手段】
【００１１】
上記課題を解決するため、本発明は、超音波画像に対し、ノイズ処理を行い、ノイズ処理後の画像に対し、その画像のダイレーション結果やエロージョン結果をもとに値を決定した構造要素を用いてモルフォロジー処理を行う。
【００１２】
即ち、本発明の超音波撮像装置は、検査対象に超音波を送信するとともに反射波である超音波を受信し、超音波信号を取得する超音波計測部と、前記超音波計測部が取得した超音波信号を用いて、検査対象の超音波画像を生成する画像生成部と、前記画像生成部が作成した超音波画像を処理する画像処理部と、を備え、前記画像処理部は、前記超音波画像のノイズを除去するノイズ除去部と、前記ノイズ除去部が生成するノイズ除去後画像についてモルフォロジー処理を行うモルフォロジー処理部と、前記ノイズ除去後画像を用いて前記モルフォロジー処理に用いる構造要素の値を決定する構造要素決定部と、を備える。
【００１３】
また本発明の画像処理方法は、超音波撮像装置で受信した超音波信号を用いて生成した超音波画像を処理する画像処理方法であって、前記超音波画像のノイズを除去するステップと、ノイズ除去後の画像に対し、モルフォロジー処理を施すステップと、前記ノイズ除去後の画像を用いて、前記モルフォロジー処理に用いる構造要素の値を決定するステップを含む。
【発明の効果】
【００１４】
本発明は、超音波画像における構造の特徴と超音波画像特有のノイズとを弁別した処理を行うことができ、ノイズが低減し且つ検査対象の構造を明確にした超音波画像を得ることができる。
【図面の簡単な説明】
【００１５】
【図１】実施形態の超音波撮像装置の全体構成を示す図。
【図２】画像処理の概要を示すフロー図。
【図３】第一実施形態の超音波撮像装置の画像処理部の構成を示すブロック図。
【図４】第一実施形態の画像処理部の処理を示すフロー図。
【図５】モルフォロジー処理を説明する図。
【図６】（ａ）〜（ｅ）は、画像処理結果を説明する図。
【図７】（ａ）〜（ｄ）は、第一実施形態の効果を説明する図。
【図８】（ａ）は、第一実施形態における処理結果、（ｂ）は従来手法による処理結果を示す図。
【図９】第二実施形態の超音波撮像装置の画像処理部の構成を示すブロック図。
【図１０】第二実施形態の画像処理の概要を示すフロー図。
【図１１】第二実施形態の画像処理部の処理を示すフロー図。
【図１２】第三実施形態の超音波撮像装置のＵＩの一例を示す図。
【発明を実施するための形態】
【００１６】
以下、本発明の超音波撮像装置及び画像処理方法の実施の形態を、図面を参照して説明する。
【００１７】
まず本発明が適用される超音波撮像装置について簡単に説明する。超音波撮像装置１００は、図１に示すように、主たる構成として、検査対象に超音波を送信するとともに反射波である超音波を受信し、超音波信号を取得する超音波計測部１０と、超音波計測部１０が取得した超音波信号を用いて、検査対象の超音波画像を生成する画像生成部２０と、画像生成部２０が作成した超音波画像を処理する画像処理部３０と、を備える。また付属装置として、検査対象に接触させて超音波の送受信を行う超音波探触子（超音波プローブ）８０が備えられる。超音波撮像装置１００は、さらに、装置の各部の動作を制御する制御部４０を備え、ユーザーが計測条件や計測の開始／停止などの指示を入力するための入力部７０、画像生成部２０や画像処理部３０が作成した画像や画像処理部３０の演算に必要なデータなどを格納する記憶部５０、画像やＧＵＩ（グラフィックユーザーインターフェイス）などを表示する表示部６０、スキャンコンバータ（不図示）などが備えられる。
【００１８】
超音波計測部１０の構成は、公知の超音波撮像装置と同様であり、検査対象の目的の内部領域に対し、超音波プローブ８０から所定周波数の超音波を送信するために、送信信号のビームフォーミングを行う送信ＢＦ部１１、検査対象の内部領域からの反射波である超音波信号を受信し、増幅・整相等の処理を行う受信ＢＦ部１２、送受信の切替を行う送受信切替部１３、などを備えている。これらの要素の一部は、超音波プローブ８０（探触子）内に設けられる場合もある。
【００１９】
画像生成部２０は、受信ＢＦ部１２から超音波信号を受け取り、輝度変換や周波数解析などを行い、表示部６０に表示される画像を生成するもので、信号処理部２１、Ｂモード画像などの画像データを生成する断層像演算部２３などを備えている。
【００２０】
画像処理部３０は、画像生成部２０が生成した画像（断層像）に対し、ノイズ除去やエッジ強調などの処理を施し、画像処理された画像を生成するもので、超音波画像のノイズを除去するノイズ除去部３１と、ノイズ除去部３１が生成するノイズ除去後画像についてモルフォロジー処理を行うモルフォロジー処理部３３と、画像データの座標を変換する座標変換部３２備え、さらに本実施形態の画像処理部３０は、ノイズ除去後画像を用いてモルフォロジー処理に用いる構造要素の値を決定する構造要素決定部３５を備えている。画像変換部３２は、画像の座標を、超音波ビームの深度方向ｒを縦軸、角度方向θを横軸とする極座標（ｒ−θ座標）から、実空間の座標（ｘ−ｙ直交座標）に変換する。座標変換部３２はノイズ除去部３１の前段であってもよいし、モルフォロジー処理部３３の後段であってもよい。
【００２１】
上述した超音波撮像装置における制御部４０及び画像処理部３０の主な機能は、ＣＰＵやＧＰＵとメモリを備えた計算機で実現することが可能であるが、画像処理部３０の機能の一部あるいは全部を超音波撮像装置とは別の計算機で実現することも可能であり、そのような実施形態も本発明に包含される。さらに計算機のソフトウェアで実現する機能の一部をＡＳＩＣやＦＰＧＡ等のハードウェアで実現してもよい。
【００２２】
上記構成の超音波撮像装置における処理の概要を図２に示す。図示するように、受信ＢＦ部１２が超音波プローブ８０を介して超音波信号を受信すると、画像生成部２０は１フレーム分の信号を画像データとしてメモリに格納する（Ｓ２１）。メモリには１ないし複数フレームの画像データが蓄積され、フレームが更新されるたびにメモリも更新される。
【００２３】
各フレームの画像データが、画像処理部３０に入力されると、画像処理部３０は、入力された画像データに対し、ノイズ除去処理を行った後（Ｓ２２）、ノイズ除去後の画像に基づき構造要素値を決定された適応的構造要素を用いてモルフォロジー処理を行う（Ｓ２３、Ｓ２４）。適応的構造要素とは、処理対象とする画像の画素位置に応じて構造要素の値及びサイズの少なくとも一方を変化させながら適用する構造要素である。このようなモルフォロジー処理を行うことにより、スペックルの偽強調がなく組織構造が明確な画像が得られる。モルフォロジー処理後の画像データは、座標変換部３２により実空間座標（ｘ−ｙ座標）に変換された後（Ｓ２５）、スキャンコンバータ７０により表示画像に変換され、撮像条件や被検体情報などの付帯情報とともに表示部６０に表示される（Ｓ２６）。
【００２４】
なお図２では、構造要素決定のステップＳ２４を、モルフォロジー処理Ｓ２５の前段に示したが、ステップＳ２４は、モルフォロジー処理Ｓ２５の間に行われる場合もある。
また、図２では、画像座標を変換する処理Ｓ２５をモルフォロジー処理後に行う場合を例示したが、画像変換処理Ｓ２５はノイズ処理前に行ってもよいし、ノイズ処理前及びモルフォロジー処理後に行ってもよい。
【００２５】
以下、上述した超音波撮像装置の構成と動作をもとにして、画像処理部３０の構成と画像処理方法の各実施形態を説明する。
【００２６】
＜第一実施形態＞
本実施形態は、構造要素決定部３５が、モルフォロジーの第一演算の結果を用いて構造要素の値を決定することが特徴である。
本実施形態の画像処理部３０は、図３に示すように、モルフォロジー処理部３３が、第一の構造要素を用いたダイレーション及びエロージョンの少なくとも一方を行う第一演算部３３Ａと、第二の構造要素を用いたオープニング及びクロージングの少なくとも一方を行う第二演算部３３Ｂとを有し、構造要素決定部３４は、ノイズ除去後画像を用いて、第二演算部３３Ｂで用いる第二の構造要素の値を決定する。本実施形態では、構造要素決定部３５を構造要素値決定部という。また第一演算部３３Ａ及び第二演算部３３Ｂの処理結果を用いてモルフォロジー処理の出力画像を生成する合成部３７をさらに備えている。なお図３では座標変換部３２は図示を省略している。
【００２７】
本実施形態の画像処理部３０は、図４に示すように、画像生成部２０が生成した画像データを受け取り、ノイズ処理（Ｓ２２）、モルフォロジー処理の第一演算（Ｓ２４１）、構造要素決定（Ｓ２３）、モルフォロジー処理の第二演算（Ｓ２４２）、合成処理（Ｓ２４３）を順次行う。モルフォロジー処理では、構造要素を画像の適応的に設定し、第一演算において、ダイレーション、エロージョンを行い、第二演算においてオープニング、クロージングの４つの処理を行う。構造要素決定ステップでは、第二演算に用いる構造要素を決定する。
【００２８】
以下、画像処理部３０における各処理の詳細を説明する。
［ノイズ除去処理：Ｓ２２］
画像データ（原画像）に対し、所定のフィルタを用いてノイズ除去を行う。ノイズ除去フィルタとしては、エッジ領域のボケを抑制しつつノイズ除去を行うエッジ保存型平滑化フィルタが好ましい。具体的には、目的画素を中心とする局所領域の画素値を、目的画素の重み係数を大きくして加重平均する加重平均フィルタ、エッジなどの構造物の方向成分を検出し, その方向に1次元平滑化処理をする方向依存型フィルタ、入力画像データをスケーリング関数で展開し、その結果を異なる解像度のウェブレット関数に展開する多重解像度解析などを採用することができ、特に多重解像度解析が好ましい。
【００２９】
［モルフォロジー処理（第一演算）：Ｓ２４１］
ノイズ除去後の画像データ（ノイズ除去後画像）に対し、ダイレーション及びエロージョンの各処理を行う。この処理では、一定のサイズで固定値を持つ構造要素を用いる。処理内容は公知のダイレーション及びエロージョンと同じであり、それぞれ、式（１）、（２）で表すことができる。
【００３０】
【数１】
2020110296.tif 000003
式中、ｆは入力データである画像データ、ｘは画像データにおける位置、ｇは構造要素、ｕは構造要素の位置を表す。「＋」及び「−」を○で囲った符号は、それぞれ、構造要素ｇを用いた画像データｆのダイレーション、エロージョンの処理を表す演算子である。
【００３１】
即ち、図５に示すように、例えば、入力画像ｆの各画素に対し、固定値を持つ構造要素ｇ（図には、サイズ：３×３、値：１を示している）を作用させて、処理対象の画素の画素値を式（１）、式（２）により決定する。これによりダイレーションでは対応する画素の構造要素の値を足して、最大値を取り処置対象の画素の画素値とし、エロージョンでは対応する画素の構造要素の値を減じて、最小値を取り、処理対象の画素の画素値とする。例えば入力画像を図６（ａ）に示すように一次元のグラフで示した場合、ダイレーションの結果は図６（ｂ）に点線で示すようになり、エロージョンの結果は一点鎖線で示すようになる。図６（ａ）の下に構造要素ｇのサイズを表している。
【００３２】
［構造要素決定：Ｓ２３］
構造要素値決定部３５が上記第一演算の結果を用いて、次の第二演算で用いる構造要素を決定する。第二演算の構造要素は、処理対象である画素の位置に応じて値を異ならせた適応的構造要素とし、エッジ領域で値が大きく、一様領域で値が小さくなるよう構造要素の値を決定する。このため、上記ステップＳ２４１のダイレーションの結果（出力）Ｙｄとエロージョンの結果（出力）Ｙｅを差分する。
【数２】
2020110296.tif 000004
【００３３】
差分した結果は、図６（ｂ）に太線で示すように、画素値が急激に変化する領域すなわちエッジ領域において値が大きく、一様な領域では値が小さくなる。本実施形態では、この差分である関数の値を構造要素値とする構造要素ｇ’を決定する。構造要素のサイズは、第一演算の構造要素と同一でもよいし、異ならせてもよい。
【００３４】
［モルフォロジー処理（第二演算）：Ｓ２４２］
ステップＳ２３で決定した構造要素ｇ’を用いて、ステップＳ２４１のダイレーション結果及びエロージョン結果に対し、それぞれ、クロージング処理とオープニング処理を行う。オープニングはエロージョン後にダイレーションする処理であり、クロージングはダイレーション後にエロージョンする処理であり、ダイレーションとエロージョンを表す演算子を用いて、式（４）、（５）で表すことができる。
【数３】
2020110296.tif 000005
【００３５】
図６（ｂ）に示す例について、オープニングとクロージングを行った結果を、図６（ｃ）、（ｄ）にそれぞれ示す。図示するように、これらの処理により、境界領域の両側で、それぞれ原画像よりも急峻になりエッジが強調されていることがわかる。
本実施形態では、この適応的構造要素ｇ’を用いたオープニング及びクロージングとは別に、第一演算で用いた構造要素ｇを用いたオープニングとクロージングを行う。その結果は次の合成処理で用いることができる。
【００３６】
［合成処理：Ｓ２４３］
最後に、合成部３７がステップＳ２４２の処理結果と入力画像とを式（６）により合成する。
【数４】
2020110296.tif 000006
合成した結果を図６（ｅ）に示す。
【００３７】
この例では、適用的構造要素ｇ’を用いたオープニング結果と構造要素ｇを用いたオープニング結果の差分、及び適用的構造要素ｇ’を用いたクロージング結果と構造要素ｇを用いたクロージング結果の差分を入力画像に加算し、出力画像としている。これによりエッジ強調効果を高めながら、ノイズが平滑化された画像を得ることができる。
【００３８】
但し、合成の仕方は式（６）に限定されず、例えば、適用的構造要素ｇ’を用いたオープニング結果と構造要素ｇを用いたオープニング結果の差分のみ、或いは、適用的構造要素ｇ’を用いたクロージング結果と構造要素ｇを用いたクロージング結果の差分のみを、入力画像に加算してもよいし、また差分をとるのではなく、適用的構造要素ｇ’を用いたオープニング結果とクロージング結果を加算してもよい。
【００３９】
以上説明したステップＳ２４１〜Ｓ２４３により画像処理部３０における処理（図２：Ｓ２３、Ｓ２４）が完了する。合成後の画像は出力画像として、表示部６０に表示される。或いは記憶部５０に格納される。
【００４０】
本実施形態の効果を、図７を用いて説明する。図７（ａ）は入力画像を模式的に示す図である。この入力画像に対し、そのままモルフォロジー処理を適用してエッジ強調を行った場合には、図７（ｂ）に示すように、エッジ部分も強調されるが、構造要素に適さないスペックルが存在した場合、一様部分（平坦な部分）のスペックルも同時に強調される場合がある。従ってその後にノイズ除去等の平滑化を行ったとしても、偽強調されたスペックルが残り、組織の境界などと区別できないことになる。これに対し、本実施形態では、図７（ｃ）、（ｄ）に示すように、ノイズ除去処理後の画像（図７（ｃ））をもとにモルフォロジー処理を行うので、本来強調すべきエッジのみ強調された画像（図７（ｄ））を得ることができる。
【００４１】
図８に、実際に超音波撮像装置で取得した画像について、本実施形態の処理を行った画像（ａ）と、ノイズ除去処理を行うことなくモルフォロジー処理を行った画像（ｂ）を示す。この図からも本実施形態の処理を行うことで、スペックルが除去され画像中央の組織が周辺組織と明確に区別できていることがわかる。
【００４２】
以上説明したように、本実施形態によれば、原画像に対しノイズ除去処理を行った後の画像に対し、適応構造要素を用いたモルフォロジー処理を行うことで、一様な領域について十分なノイズ除去を実現しつつ、エッジ領域について高いエッジ強調効果を得ることができる。これによりスペックルノイズに起因する画質の低下を防止できる。
【００４３】
＜第二実施形態＞
第一実施形態では、モルフォロジーに用いる適応的構造要素の値をモルフォロジーの第一演算結果を用いて決定したが、本実施形態は、位置に応じてサイズを異ならせた構造要素を決定し、用いる。
【００４４】
一般にモルフォロジー処理は構造要素内の成分を平滑化する作用があるため、構造要素のサイズをスペックルのサイズに相当するサイズにすることが望ましい。スペックルのサイズは、撮像に用いる超音波プローブの種類、超音波信号の周波数、撮像条件などによって画像毎に異なる。また深度やサンプリング間隔などによってもスペックルのサイズは変化するため、同一画像内でも位置によって異なる。本実施形態では、撮像条件に応じて構造要素のサイズを決定することにより、効果的にスペックルノイズの除去を行う。
【００４５】
本実施形態の画像処理部の構成は、第一実施形態の画像処理部の構成と基本的には同じであるが、図９に示すように、サイズ決定部３９が追加される。図９において図３と同じ要素は同じ符号で示し、重複する説明は省略する。また図９でも座標変換部３２は図示を省略している。
【００４６】
本実施形態の画像処理部３０は、図１０に示すように、超音波画像データを取得した後（Ｓ４１）、後のモルフォロジー処理に用いる構造要素のサイズを決定するために必要な情報（例えば撮像パラメータ）を取得し（Ｓ４２）、サイズを決定する（Ｓ４３）。その後、ノイズ除去処理、モルフォロジー処理（第一演算、第二演算）（Ｓ４４）、画像変換処理を行う。
【００４７】
以下、本実施形態の特徴であるサイズ決定部３９の処理を説明する。
サイズ決定部３９は、画像データを受け取ると、その画像データを取得したときの撮像条件等から、画像の分解能に関する情報とサンプリング間隔を取得する。画像の分解能は深度及び方位により異なり、それらは、超音波プローブ、周波数、焦点位置、撮像モードなどの条件によって決まる。超音波プローブの種類や撮像条件ごとに、予め深度と分解能との関係を音場数値計算や画像シミュレーションや実測などにより得ておくことで、画像のパラメータで決まる深度から分解能を得ることができる。深度と分解能との関係は例えテーブルＬＵＴとして記憶部５０に格納しておくことが好ましく、ＬＵＴを参照して取得した深度に対応する分解能を得る。
【００４８】
サンプリング間隔は（ｒ−θ座標系またはｘ−ｙ座標系における）入力画像上で隣接するデータ間隔の実空間における距離であり、その情報は入力画像から直接得ることができる。例えば、入力画像がｒ−θ座標系（座標変換前）であり、入力画像上でのθ方向のデータ間隔が１degreeであった場合、ｒ＝３０ｍｍにおけるθ方向のサンプリング間隔Δｓは
Δｓ＝３０[mm]×(１／１８０×π)＝０．５[mm]
となる。
【００４９】
サイズ決定部３９は、分解能をＲ、サンプリング間隔をΔｓとすると、例えば次式（７）により構造要素のサイズＮを決定する。
Ｎ＝Ｒ／Δｓ        （７）
【００５０】
分解能Ｒは画素の位置によって異なるので、画素位置に応じてサイズを変化させることとしてもよいが、全画素について求めた「Ｒ／Δｓ」の最大値、最小値、或いは平均値などの統計値を算出し、いずれかを構造要素のサイズと決定してもよい。或いは全画素を「Ｒ／Δｓ」の値によって複数のグループに分けて、各グループの「Ｒ／Δｓ」（＝Ｎ）の例えば平均値を算出し、平均値を各ブループに適用する構造要素のサイズＮとしてもよい。平均値ではなく、最大値や最小値を適用してもよい。
【００５１】
決定した構造要素のサイズは、モルフォロジー処理の第一演算及び第二演算の両方で用いる構造要素に適用することができる。モルフォロジーの演算は第一実施形態と同様である。モルフォロジー処理の第一演算部３３Ａは、サイズ決定部３９が決定したサイズの構造要素を用いて、ノイズ除去後画像に対し、ダイレーション及びエロージョンを行う。ついで第二演算部３３Ｂは、サイズ決定部３９が決定したサイズの構造要素を用いて、オープニング及びクロージングを行う。最後に、オープニング／クロージング結果と入力画像とを前掲の式（６）により合成する。
【００５２】
あるいは、第一実施形態の構造要素値決定処理（Ｓ２３）を行わない場合、第一演算のダイレーション結果Ｙｄ及びエロージョン結果Ｙｅに対して、Ｙｄ・Ｙｅのうち原画像に近い方を画素ごとに重み付け選択して合成してもよい。
【００５３】
各演算に用いる構造要素値は、固定値でもよいし、位置に応じて固定値を異ならせてもよい。本実施形態の画像処理部３０が、第一実施形態の画像処理部３０と同様の構造要素決定部３５を備える場合には、図１１に示すように、モルフォロジー処理の第一演算（Ｓ２４１）におけるダイレーション結果Ｙｄとエロージョン結果Ｙｅとの差分を用いて、画像の位置によって値が変わる構造要素ｇ’を決定し（Ｓ２３）、このような適応的構造要素ｇ’を用いて第二演算を行ってもよい（Ｓ２４２）。その場合には、式（６）に示したように、第一演算と同じ構造要素ｇを用いたオープニング結果と適応的構造要素ｇ’を用いたオープニング結果との差分、及び、第一演算と同じ構造要素ｇを用いたクロージング結果と適応的構造要素ｇ’を用いたクロージング結果との差分、をそれぞれ、ノイズ除去後画像に加算して最終的な画像を作成する。
【００５４】
なお、本実施形態においても、画像座標を変換する処理はモルフォロジー処理の前、後のいずれでもよい。ただし、モルフォロジー処理後に画像変換を行う場合には、構造要素のサイズを決定する分解能は、ｒ−θ座標における分解能の値を用いる。
【００５５】
本実施形態によれば、第一実施形態と同様の、ノイズ除去後画像に対しモルフォロジー処理を行うことにより、スペックルの偽強調を防止したエッジ強調と平滑化をともに実現できるという効果に加え、スペックルサイズに相当するサイズの構造要素を用いてモルフォロジー処理に用いることにより、スペックルノイズの影響を最小限にしてエッジ領域のエッジ強調と一様領域の平滑化を適切に行うことができ、画質の良好な画像を得ることができる。
【００５６】
また本実施形態によれば、さらに第一実施形態の適応的構造要素（位置に応じて値を変化させる構造要素）を用いた場合には、エッジ領域でのエッジ強調効果と一様領域での平滑化をさらに向上することができる。但し、このような適応的構造要素を用いない場合も本実施形態に包含される。
【００５７】
なお第一実施形態及び第二実施形態において、画像処理の対象とする超音波画像は、Ｂモード画像のみならず、ＴＨＩ（組織ハーモニック像）や血流情報を含む断層像などでもよい。
【００５８】
＜第三実施形態＞
本実施形態は、構造要素の値やサイズの決定を、ユーザーが選択可能にする構成を追加したことが特徴である。このため本実施形態の超音波撮像装置は、図１に示したように、例えば制御部４０にＵＩ（Ｕｓｅｒ  Ｉｎｔｅｒｆａｃｅ）部が含まれる。
【００５９】
ＵＩ部は、ユーザー選択を可能にするグラフィックＵＩ（ＧＵＩ）を表示部６０に表示するとともに、入力部７０に備えられたマウスやキーボード、タッチパネル等の入力デバイスの操作によってＧＵＩ上のボタン等の操作を受付け、受付けたユーザー指示に従い、画像処理部３０の処理を決定する。
【００６０】
ＧＵＩの一例を図１２に示す。この例では、ノイズ除去後画像について第一実施形態或いは第二実施形態のモルフォロジー処理を行うか否かの選択をユーザーから受付ける。即ち、表示部６０の画面６００には、画像を表示する画像表示ブロック６１０と、撮像条件やパラメータを表示するパラメータ表示ブロック６２０と、ユーザーが操作可能なボタンを表示する受付部６３０とが含まれる。例えば、画像表示ブロック６１０にノイズ除去後のＢモード画像が表示されると、ユーザーはその画像から診断に十分な情報が得られていると判断した場合やそれ以上の画像処理を行いたくないと判断した場合には、受付け部６３０に表示されたエッジ強調処理のオンオフボタンを操作し、「ＯＦＦ」にする。エッジ強調処理が必要と判断した場合には、エッジ強調処理を「ＯＮ」にする。これにより、ノイズ除去後画像に対し、モルフォロジー処理部３３は例えば第一実施形態に基づくモルフォロジー処理を行い、その結果画像を画像表示ブロック６１０に表示する。この際、処理前のノイズ除去後画像と処理後の画像とを並列に表示してもよい。
【００６１】
さらに「平滑」の程度（Ｌｅｖｅｌ）を調整するための調整用ダイヤルや調整バーを設けてもよく、受付けた平滑の程度に応じて、構造要素サイズを調整してもよい。この場合にも、調整した構造要素のサイズを用いて行ったモルフォロジー処理の結果を、単独で或いは処理前のノイズ除去後画像とともに画像表示ブロック６１０に表示してもよい。
【００６２】
本実施形態によれば、ユーザーにより処理の自由度を高めることができ、またエッジ強調処理の前後や調整度合いを比較することで、より適切な画像処理を行うことができる。
【００６３】
＜第四実施形態＞
本実施形態は、上述した各実施形態の画像処理部の構成は、特に、検査対象の内部に穿針を挿入し、穿針の進路をガイドする支援画像を撮像する超音波撮像装置に適用した実施形態である。
【００６４】
生体組織とともに穿刺器具を画面上に表示する超音波撮像装置では、穿刺針は超音波画像上で生体組織よりも相対的に高輝度な点として描出されることで視認されるが、生体内に刺入された穿刺針は、針からの反射波が入射波の方向と異なる方向に戻ってしまうため視認性が十分でない場合がある。そこで超音波ガイド下穿刺を行う場合に、入力画像に対して、ノイズ除去とその結果を利用した適応構造要素によるエッジ強調とを適用することで、スペックルノイズを低減したうえで穿刺針のエッジ成分がより鮮鋭化された超音波画像を得る。
【００６５】
即ち、本実施形態の超音波撮像装置は、穿針の進路をガイドする支援画像を撮像する超音波撮像装置であり、検査対象に超音波を送信するとともに反射波である超音波を受信し、超音波信号を取得する超音波計測部と、超音波計測部が取得した超音波信号を用いて、検査対象の超音波画像を生成する画像生成部と、画像生成部が作成した超音波画像を処理する画像処理部と、を備える。画像処理部は、画像生成部が作成した超音波画像のノイズを除去するノイズ除去部と、ノイズ除去部が生成するノイズ除去後画像についてモルフォロジー処理を行うモルフォロジー処理部と、ノイズ除去後画像を用いてモルフォロジー処理部で用いる構造要素の値を決定する構造要素決定部とを備え、穿針とそれ以外の組織との境界を強調した支援画像を生成する。
【００６６】
超音波計測部、画像生成部及び画像処理部（ノイズ除去部、モルフォロジー処理部、構造要素決定部）の構成は、第一実施形態の超音波撮像装置における対応する要素と同じであり、画像処理のフローは図４で説明したものと同様である。また画像処理部は、撮像条件（画像の分解能やサンプリング間隔）に応じて構造要素のサイズを決定するサイズ決定部３９（図９）を備えていてもよく、その場合の画像処理のフローは図１０又は図１１で説明したものと同様である。さらに必要に応じて、方向依存性平滑処理などを用いた公知の穿刺針強調処理等を加えてもよい。
【００６７】
本実施形態によれば、穿刺ガイド支援用の超音波撮像装置に本発明の画像処理を適用することで、穿刺針のエッジ成分がより鮮鋭化された超音波画像を提供することができ、穿刺時の針の視認性を向上させることができる。
【００６８】
以上、本発明の超音波撮像装置と画像処理方法の実施形態を説明したが、各実施形態で説明した要素のうち、いくつかの要素を省略したり、公知の構成を追加したりすることも本発明に包含される。
【符号の説明】
【００６９】
１００：超音波撮像装置、１０：超音波計測部、２０：画像生成部、３０：画像処理部、３１：ノイズ除去部、３２：座標変換部、３３：モルフォロジー処理部、３３Ａ：第一演算部、３３Ｂ：第二演算部、３５：構造要素決定部、３７：合成部、３９：サイズ決定部、４０：制御部、５０：記憶部、６０：表示部、７０：入力部、８０：超音波プローブ。
【図１】
2020110296.tif 000007
【図２】
2020110296.tif 000008
【図３】
2020110296.tif 000009
【図４】
2020110296.tif 000010
【図５】
2020110296.tif 000011
【図６】
2020110296.tif 000012
【図７】
2020110296.tif 000013
【図８】
2020110296.tif 000014
【図９】
2020110296.tif 000015
【図１０】
2020110296.tif 000016
【図１１】
2020110296.tif 000017
【図１２】
2020110296.tif 000018