JPH10146325A

JPH10146325A - 画像補正処理方法

Info

Publication number: JPH10146325A
Application number: JP8310689A
Authority: JP
Inventors: Tetsuhiko Takahashi; 哲彦高橋
Original assignee: Hitachi Medical Corp
Current assignee: Hitachi Healthcare Manufacturing Ltd
Priority date: 1996-11-21
Filing date: 1996-11-21
Publication date: 1998-06-02
Anticipated expiration: 2016-11-21
Also published as: JP3689509B2

Abstract

(57)【要約】【課題】任意の形状、角度の画像に適用でき、エッジの
過補正のないシェーディング補正アルゴリズムを提供す
る。【解決手段】被写体と背景を有する２次元または、３次
元の濃淡画像である原画像のシェーディングを補正する
画像補正処理方法であって、原画像を被写体領域と背景
領域とに領域分割する手段（101）と、被写体と背景の
境界を抽出する手段（102）と、原画像から境界の近傍
の絵素を等方的に抽出する手段（103）と、抽出された
絵素のうち、背景の絵素の絵素値を、被写体の絵素値か
ら求めた信号値で置換し第２の画像を形成する手段（10
4）と、第２の画像に、低周波通過フィルタを作用させ
る手段（105）と、低周波通過フィルタ作用後の第３の
画像を使って原画像の輝度を補正する手段（106）とを
含む。

Description

【発明の詳細な説明】

【０００１】

【発明の属する利用分野】本発明は、濃淡画像のシェー
ディングを補正する方法に関し、特に被検体中の水素や
燐等からの核磁気共鳴（以下、「NMR」という）信号を
測定し、核の密度分布や緩和時間分布等を映像化する核
磁気共鳴撮影（ＭＲＩ）方法等により得られた人体断面
像等に適用される画像補正処理方法及びこのような画像
補正処理手段を備えたＭＲＩ装置に関する。

【０００２】

【従来の技術】ＭＲＩでは、主として、被検体の主たる
構成物質であるプロトンを撮影対象として、プロトン密
度の空間分布や、励起状態の緩和現象の空間分布を画像
化することで、人体頭部、腹部、四肢等の形態または、
機能を２次元もしくは３次元的に撮影する。

【０００３】このＭＲＩ画像は、臨床で広く活用されて
いるが、装置に起因するシェーディングが強く出る場合
がある。特に、局所ＲＦ受信コイルを使った場合、受信
コイルの感度分布に起因するシェーディングが顕著で、
診断がしにくい等の問題が生じる場合があった。例え
ば、図９に示すファントムの原画像のプロファイルは、
本来フラットなるべきところ凹となっており、受信コイ
ルの感度分布を反映している。これを回避する手段とし
て、シェーディング補正が提案されている。一例とし
て、画像の低周波成分を抽出し、これを装置起因のシェ
ーディングとみなし、原画像を補正する方法が提案され
ている（アクセル他、「表面コイルＭＲイメージングに
おける強度補正」アメリカン・ジャーナル・オブ・レン
トゲノロジー、148巻418〜420頁、1987年）。

【０００４】上記方法は簡便であるが、原画像の低周波
成分を抽出したものは、被写体と背景との境界近傍にお
いて実際のシェーディングとずれが生じ、これで原画像
を補正した場合、被写体の周辺が過補正になり、図９に
示すように補正後の画像で、被写体の周辺部が高輝度に
なる欠点があった。

【０００５】これを解決する方法として、被写体と背景
の境界から外側の領域を最も高い絵素値で置換した２次
画像を求め、この画像の低周波成分を抽出し、これによ
りシェーディング補正する方法が提案されている（ウォ
ルド他、「ヒト脳の高解像度ＭＲイメージングのため
の、フェイズドアレイ検出器及び自動強度補正アルゴリ
ズム」マグネティック・レゾナンス・イン・メディス
ン、34巻433〜439頁、1995年）。

【０００６】

【発明が解決しようとする課題】この方法では頭部横断
像で効果が示されているものの、被写体と背景の境界が
入組んでいるような複雑な形状の被写体には適用でき
ず、その適用は頭部のような単純な形状の被写体（円
形）に適用が止まっている。一方、ＭＲＩではあらゆる
部位をあらゆる角度から撮影するため、補正アルゴリズ
ムは、任意画像に対する適用が求められるが、これを可
能にするアルゴリズムはなかった。

【０００７】そこで本発明は、任意の画像に適用できる
シェーディング補正アルゴリズムを提供することを目的
とする。

【０００８】

【課題を解決するための手段】上記課題は、被写体と背
景を有する２次元または、３次元の濃淡画像である原画
像のシェーディングを補正する画像補正処理方法におい
て、原画像を被写体領域と背景領域とに領域分割する手
段（101）と、被写体と背景の境界を抽出する手段（10
2）と、原画像から境界の近傍の絵素を等方的に抽出す
る手段（103）と、抽出された絵素のうち、背景の絵素
の絵素値を、被写体の信号値で置換し第２の画像を形成
する手段（104）と、第２の画像に、低周波通過フィル
タを作用させる手段（105）と、低周波通過フィルタ作
用後の第３の画像を使って原画像の輝度を補正する手段
（106）とを含む本発明の画像補正処理方法によって解
決される。

【０００９】境界近傍の背景を等方的に抽出することに
より、あらゆる角度の画像に適用できる。ここで等方的
に抽出するとは、２次元画像であれば、境界点を中心に
固有の大きさの正方領域或いは円形領域以内の絵素を抽
出する。また３次元画像であれば、境界点を中心に固有
の大きさの立方領域或いは球形領域以内の絵素を抽出す
る。

【００１０】また背景の絵素の絵素値を、被写体の信号
値で置換する手段は、被写体の信号値を被写体である絵
素の絵素値を元に計算して求め、この信号値を背景の絵
素の絵素値とする。この場合、被写体の信号値として
は、被写体全体の絵素値を元に求めておいてもよいが、
等方的に抽出した絵素のうち被写体である絵素の絵素値
を元に平均値或いはメディアン値等を計算することによ
り求めることができる。

【００１１】このように境界近傍の背景の絵素を被写体
の信号値で置換した第２の画像に低周波通過フィルタを
作用させた画像は、境界近傍において実際のシェーディ
ングとほぼ一致したものとなり、且つこのような処理が
等方的になされているので、原画像がどの方向にシェー
ディングを持っていても、この画像を用いて原画像を補
正することにより、過補正の防止された理想的な補正を
行うことができる。

【００１２】また本発明のＭＲＩ装置は、静磁場、傾斜
磁場、高周波磁場の各磁場を発生する磁場発生手段と、
被検体を構成する組織の原子核スピンに核磁気共鳴を生
じさせる高周波磁場を照射するために磁気発生手段を制
御するシーケンサと、被検体から生じた核磁気共鳴信号
を検出する手段と、核磁気共鳴信号を信号処理して画像
再構成する画像処理手段と、画像を表示する手段とを備
え、画像処理手段は、上記画像補正処理方法を実行する
シェーディング補正手段を備えている。

【００１３】本発明のＭＲＩ装置は、このようなシェー
ディング補正手段を備えていることにより、検出する手
段の感度分布に起因するＭＲ画像のシェーディングを過
補正なく補正することができ、診断に有効な画像を提供
できる。

【００１４】

【発明の実施の形態】本発明の画像補正処理方法の実施
例を図１に示すフロー図を参照して以下説明する。図２
（ａ）は、ＭＲＩ画像（原画像）200を模式的に示した
図で、被写体201部分は高輝度で背景部分202は低輝度の
濃淡画像となっている。

【００１５】本発明の画像補正処理方法では、まずこの
ような濃淡画像200の高輝度領域を被写体として抽出
し、低輝度領域を背景として抽出して領域分割する（手
段101）。この領域分割は、例えば絵素値の最大値の５
０％を閾値として、画像を２値化し、ゼロを背景、1を
被写体とする。このように得られた２値化画像を図２
（ｂ）に示す。尚、領域分割する際の閾値は、最大値の
５０％に限るものではなく、閾値を決める手法として、
モード法、百分率濃度分離法と公知の手法を採用するこ
とができる。

【００１６】次に２値化画像300を元に被写体201と背景
202の境界を抽出する（手段102）。境界401を抽出する
手段102としては、２値化画像300に微分フィルタを作用
させることができ、好適にはラプラスフィルタを用い
る。これにより図３に示すような境界401の絵素値が１
で他（背景402）の絵素値はゼロである画像（境界画
像）400が得られる。この場合、必要に応じて公知の細
線化処理を行ってもよい。

【００１７】次いで原画像200の境界近傍の絵素を等方
的に抽出する（手段103）。このために、まず図４
（ａ）に示すように手段102で得られた境界画像400を順
次走査し、絵素値が１である番地（即ち、境界の番地）
では、その周囲の特定の範囲にある絵素の番地を抽出す
る。ここで抽出される特定の範囲としては、例えば５×
５或いは７×７の正方領域とすることができる。図５
（ａ）に５×５の抽出された領域400aを示す。ここでは
境界301（斜線）である番地［３．３］の周囲に２４の
番地が抽出されている。

【００１８】図５（ｂ）はこの抽出された領域（番地）
の原画像200aを示すもので、次の手段（104）では、こ
の領域内の原画像の絵素のうち、背景の絵素値を、境界
近傍の原画像の被写体の信号値で置き換える。この処理
を以下、エッジフィルと呼ぶ。このエッジフィル処理の
ために、まず抽出された領域の絵素を被写体と背景とに
分ける。これは各番地と対応する２値化画像300の絵素
値が、１ならば被写体、０なら背景とする。図５（ｂ）
では背景の絵素を斜線で示している。次に領域内の原画
像の被写体の絵素値から信号値を計算する。被写体の信
号値は、例えば絵素値の単純平均値或いはメディアン値
を計算して求め、この信号値を原画像の背景の絵素値の
値とする。図５（ｂ）に示す実施例では、被写体とされ
た１０の絵素の絵素値ａ25、ａ34、ａ35、ａ43、ａ44、
ａ45、ａ52、ａ53、ａ54、ａ55から被写体の信号値が計
算され、背景である１５の絵素の絵素値と置換される。

【００１９】上述したエッジフィル処理は、図４
（ｂ）、（ｃ）に示すように境界全体にわたって行わ
れ、最終的に、図４（ｃ）に示すような境界全体にわた
って背景の絵素値を置換した画像（第２の画像）500が
得られる。尚、エッジフィル処理は境界となる絵素全部
について行ってもよいが、離散的に処理してもよい。

【００２０】このようにして得られた第２の画像500を
低周波通過フィルタ（以下、ＬＰＦと略す）を作用させ
る（手段105）。このため第２の画像を２次元フーリエ
変換し、周波数空間画像とし、これに２次元バターワー
スフィルタ、２次元ガウスフィルタ、２次元ハニングフ
ィルタ等のＬＰＦを作用し、再度２次元逆フーリエ変換
する。このように第２の画像500にＬＰＦを作用させた
第３の画像は、原画像の外側にエッジフィルを行ってい
るので、被写体の周辺において本来のシェーディングに
近いシェーディングとなる。

【００２１】最後に手段105で得られた第３の画像（シ
ェーディング画像）を用いて原画像200の輝度を補正す
る（手段106）。この補正は、例えば、原画像200をシェ
ーディング画像で除することにより行われる。これによ
りシェーディングを補正した補正画像を得ることができ
る。

【００２２】尚、以上の説明では２次元画像に対し、本
発明の補正処理方法を適用した例について説明したが、
本発明は３次元画像に対しても同様に適用できる。この
場合、輪郭（境界）を抽出後、境界となる絵素の周囲に
等方的に３次元領域を抽出し、この領域についてエッジ
フィル処理を行う。この場合にも、境界となる絵素全部
ではなく離散的にエッジフィルを施してもよい。

【００２３】また等方的に抽出する領域は、２次元の場
合、正方領域ではなく円形領域としてもよい。３次元の
場合には立方領域や球形領域とすることができる。

【００２４】更に本発明の補正処理方法は、装置に起因
してシェーディングを有する画像の補正に適用できる
が、ＭＲＩ装置における画像処理に適用する場合につい
て以下説明する。

【００２５】図６は、本発明が適用されるＭＲＩ装置の
全体概要を示す図で、このＭＲＩ装置は、大別すると、
中央処理装置（ＣＰＵ）１と、シーケンサ２と、送信系
３と、静磁場発生磁石４と、傾斜磁場発生系２１と、受
信系５と、信号処理系６とを備えている。

【００２６】ＣＰＵ１は、予めキーボード２２等の入力
装置から撮影パラメータを入力し、プログラムに従って
シーケンサ２、送信系３、受信系５、信号処理系６の各
々を制御するものである。シーケンサ２は、ＣＰＵ１か
らの制御指令に基づいて動作し、被検体７の断層画像の
データ収集に必要な種々の命令を送信系３、傾斜磁場発
生系２１及び受信系５に送るようにしている。

【００２７】送信系３は、高周波発信器８と変調器９と
高周波コイルとしての照射コイル１１を有し、シーケン
サ２の指令により高周波発信器８からの高周波パルスを
変調器９で振幅変調し、この振幅変調された高周波パル
スを高周波増幅器１０を介し増幅して照射コイル１１に
供給することにより、所定のパルス状の電磁波を被検体
７に照射するようにしている。

【００２８】静磁場発生磁石４は、被検体７の回りに任
意の方向に均一な静磁場を発生させるためのものであ
る。この静磁場発生磁石の内部には、照射コイル１１の
他、傾斜磁場を発生させる傾斜磁場コイル１３と、受信
系５の受信コイル１４が設置されている。傾斜磁場発生
系２１は互いに直交するデカルト座標軸方向にそれぞれ
独立に傾斜磁場を印加できる構成を有する３方向の傾斜
磁場コイル１３と傾斜磁場コイルに電流を供給する傾斜
磁場電源１２と、傾斜磁場電源１２を制御するシーケン
サ２により構成される。

【００２９】受信系５は、高周波コイルとしての受信コ
イル１４と該受信コイル１４に接続された増幅器１５と
直交位相検波器１６とＡ／Ｄ変換器１７とを有し、被検
体７からのＮＭＲ信号を受信コイル１４が検出すると、
その信号を増幅器１５、直交位相検波器１６、Ａ／Ｄ変
換器１７を介しデジタル量に変換するとともに、シーケ
ンサ２からの指令によるタイミングで直交位相検波器１
６によってサンプリングされた二系列の収集データに変
換してＣＰＵ１に送るようにしている。

【００３０】信号処理系６は、光ディスク１９、磁気デ
ィスク２０等の外部記憶装置と、ＣＲＴ等からなるディ
スプレイ１８とを有し、受信系５からのデータがＣＰＵ
１に入力されると、ＣＰＵ１が信号処理、画像再構成等
の処理を実行し、その結果の被検体７の所望の断面像を
ディスプレイ１８に表示するとともに、外部記憶装置の
磁気ディスク２０等に記録する。

【００３１】ＣＰＵ１は、また図１に示すような画像補
正処理アルゴリズムを有しており、受信コイル１４の感
度分布に起因する画像のシェーディングを補正する。

【００３２】このようなＭＲＩ装置においてソレノイド
型受信コイルを用いて撮影したファントム画像を例とし
て、本発明による補正処理方法と、原画像にＬＰＦを作
用させたものをシェーディングとする従来の補正処理方
法とを比較して説明する。図７に示すグラフは、ファン
トム600のＡ−Ａ’断面のｘ方向位置を横軸をとし、信
号値を縦軸としたグラフであり、図において当然ながら
信号はファントム部分だけであり、原画像のプロファイ
ルは受信コイル700の感度分布（シェーディング）に起
因して凹となっている。従来法ではこの原画像に直接Ｌ
ＰＦを作用させてコイルの感度分布を推定しているの
で、この場合には、図中点線で示すように被写体のエッ
ジ部分（斜線で示す部分）で感度分布が低くなり、本来
の感度分布との差が生じる。これに対し、本発明の補正
処理方法では、原画像にエッジフィル処理を行い、被写
体の外側に被写体の信号値を等方的に拡張しているた
め、ＬＰＦを作用しても、被写体の周辺では推定感度分
布が低下することなく実際の感度分布にほぼ一致してい
る。

【００３３】図８は、本発明の補正処理方法により推定
された感度分布を用いて、ファントムのＭＲ画像のシェ
ーディングを補正した結果を示すものである。原画像の
プロフィアルは不均一性｛（最大信号値−最小信号値）
／最小信号値＝ｂ／ａ｝が０．３３であったが、シェー
ディング補正により、０．１１（ｂ’／ａ’）に改善で
き、図９に示すような従来の補正画像におけるエッジの
強調は見られない。

【００３４】本発明の補正処理方法は、頸部画像、肩部
画像、顎間接画像、頭部画像、脊椎画像等あらゆる部位
のＭＲ画像に有効であり、また各種局所コイル、マルチ
プルコイルのシェーディング補正の有効である。

【００３５】

【発明の効果】以上の説明からわかるように、本発明の
補正処理方法によれば、原画像の境界を抽出するととも
に、この境界近傍の背景の絵素を等方的に抽出して、信
号値を置換する処理（エッジフィル処理）をすることに
より、エッジ部分の過補正がなくなり、理想的なシェー
ディング補正が可能となる。また被写体領域が、等法的
に拡張されるので、どの向きに被写体がシェーディング
を持っていても、理想的な補正ができ、あらゆる部位、
あらゆる角度の被写体にも適用できる。

【図面の簡単な説明】

【図１】本発明の一実施例を示すフロー図。

【図２】図１のフロー図における手段101の処理を説明
する図で、（ａ）は原画像を示す図、（ｂ）は２値化画
像を示す図。

【図３】図１のフロー図における手段102の処理を説明
する図。

【図４】図１のフロー図における手段103、104の処理を
説明する図で、（ａ）は境界画像からの領域抽出を示す
図、（ｂ）及び（ｃ）はそれぞれ境界近傍の背景の絵素
値の置換を示す図。

【図５】本発明によるエッジフィルの一実施例を説明す
る図。

【図６】本発明が適用されるＭＲＩ装置の全体概要を示
すブロック図。

【図７】本発明の原理を説明する図。

【図８】本発明の補正処理方法による処理の一実施例を
示すグラフ。

【図９】従来の補正処理方法による処理を示すグラフ。

【符号の説明】

１・・・・・・ＣＰＵ（画像処理手段）２・・・・・・シーケンサ２１・・・・傾斜磁場発生系（磁場発生手段）３・・・・・・送信系４・・・・・・静磁場磁石（磁場発生手段）５・・・・・・受信系６・・・・・・信号処理系（画像処理手段） 200・・・・・原画像 201・・・・・被写体 202・・・・・背景 300・・・・・２値化画像 400・・・・・境界画像 401・・・・・境界 500・・・・・第２の画像

Claims

【特許請求の範囲】

【請求項１】被写体と背景を有する２次元または、３次
元の濃淡画像である原画像のシェーディングを補正する
画像補正処理方法において、前記原画像を被写体領域と
背景領域とに領域分割する手段（101）と、前記被写体
と前記背景の境界を抽出する手段（102）と、前記原画
像から前記境界の近傍の絵素を等方的に抽出する手段
（103）と、前記抽出された絵素のうち、背景の絵素の
絵素値を、被写体の信号値で置換し第２の画像を形成す
る手段（104）と、前記第２の画像に、低周波通過フィ
ルタを作用させる手段（105）と、前記低周波通過フィ
ルタ作用後の第３の画像を使って前記原画像の輝度を補
正する手段（106）とを含む画像補正処理方法。
【請求項２】前記境界の近傍の絵素を等方的に抽出する
手段（103）は、境界点を中心に固有の大きさの正方領
域以内の絵素を抽出する手段である請求項１記載の画像
補正処理方法。
【請求項３】前記被写体の信号値は、等方的に抽出され
た絵素のうち被写体である絵素の絵素値の平均値或いは
メディアン値であることを特徴とする請求項１又は２記
載の画像補正処理方法。
【請求項４】静磁場、傾斜磁場、高周波磁場の各磁場を
発生する磁場発生手段と、被検体を構成する組織の原子
核スピンに核磁気共鳴を生じさせる高周波磁場を照射す
るために前記磁場発生手段を制御するシーケンサと、前
記被検体から生じた核磁気共鳴信号を検出する手段と、
前記核磁気共鳴信号を信号処理して画像再構成する画像
処理手段と、画像を表示する手段とを備えた磁気共鳴イ
メージング装置において、前記画像処理手段は、前記原
画像を被検体領域と背景領域とに領域分割する手段（10
1）と、前記被検体と前記背景の境界を抽出する手段（1
02）と、前記原画像から前記境界の近傍の絵素を等方的
に抽出する手段（103）と、前記抽出された絵素のう
ち、背景の絵素の絵素値を、被検体の信号値で置換し第
２の画像を形成する手段（104）と、前記第２の画像
に、低周波通過フィルタを作用させる手段（105）と、
前記低周波通過フィルタ作用後の第３の画像を使って前
記原画像の輝度を補正する手段（106）とを含むシェー
ディング補正手段を備えたことを特徴とする磁気共鳴イ
メージング装置。