JP6744764B2 - 画像診断装置、及び画像取得方法 - Google Patents

画像診断装置、及び画像取得方法 Download PDF

Info

Publication number
JP6744764B2
JP6744764B2 JP2016105023A JP2016105023A JP6744764B2 JP 6744764 B2 JP6744764 B2 JP 6744764B2 JP 2016105023 A JP2016105023 A JP 2016105023A JP 2016105023 A JP2016105023 A JP 2016105023A JP 6744764 B2 JP6744764 B2 JP 6744764B2
Authority
JP
Japan
Prior art keywords
correction
image
unit
diagnostic apparatus
noise reduction
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
Application number
JP2016105023A
Other languages
English (en)
Other versions
JP2017209329A5 (ja
JP2017209329A (ja
Inventor
喜実 野口
喜実 野口
昌宏 荻野
昌宏 荻野
毅倫 村瀬
毅倫 村瀬
高橋 哲彦
哲彦 高橋
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hitachi Ltd
Original Assignee
Hitachi Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Hitachi Ltd filed Critical Hitachi Ltd
Priority to JP2016105023A priority Critical patent/JP6744764B2/ja
Priority to PCT/JP2017/014966 priority patent/WO2017203875A1/ja
Priority to US16/093,680 priority patent/US10732245B2/en
Priority to CN201780031293.1A priority patent/CN109152549B/zh
Publication of JP2017209329A publication Critical patent/JP2017209329A/ja
Publication of JP2017209329A5 publication Critical patent/JP2017209329A5/ja
Application granted granted Critical
Publication of JP6744764B2 publication Critical patent/JP6744764B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/54Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/5608Data processing and visualization specially adapted for MR, e.g. for feature analysis and pattern recognition on the basis of measured MR data, segmentation of measured MR data, edge contour detection on the basis of measured MR data, for enhancing measured MR data in terms of signal-to-noise ratio by means of noise filtering or apodization, for enhancing measured MR data in terms of resolution by means for deblurring, windowing, zero filling, or generation of gray-scaled images, colour-coded images or images displaying vectors instead of pixels
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/05Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves 
    • A61B5/055Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves  involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/02Devices for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03Computerised tomographs
    • A61B6/032Transmission computed tomography [CT]
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/52Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5205Devices using data or image processing specially adapted for radiation diagnosis involving processing of raw data to produce diagnostic data
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus for radiation diagnosis, e.g. combined with radiation therapy equipment
    • A61B6/52Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5258Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5207Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of raw data to produce diagnostic data, e.g. for generating an image
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5269Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving detection or reduction of artifacts
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/54Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/54Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/565Correction of image distortions, e.g. due to magnetic field inhomogeneities
    • G01R33/56509Correction of image distortions, e.g. due to magnetic field inhomogeneities due to motion, displacement or flow, e.g. gradient moment nulling
    • G06T5/70
    • G06T5/73
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48NMR imaging systems
    • G01R33/54Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/546Interface between the MR system and the user, e.g. for controlling the operation of the MR system or for the design of pulse sequences
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10072Tomographic images
    • G06T2207/10088Magnetic resonance imaging [MRI]

Description

本発明は画像診断装置に係り、画像処理における高画質化技術に関する。
医用画像診断装置の一つである磁気共鳴イメージング(Magnetic Resonance Imaging: MRI)は、核磁気共鳴(Nuclear Magnetic Resonance: NMR)現象を利用して生体などの被検体の情報を画像にする方法である。臨床において対象となる部位を正確に読影するためには十分な分解能が必要となるが。MRI装置はその原理上、分解能を上げると信号対雑音比(Signal to Noise Ratio: SNR)が低下する課題がある。そこで一般的なMRI装置においては同一部位、すなわち同一位置を複数回撮像し加算する信号加算回数(Number of EXcitations: NEX)を増加させることにより、SNRを向上させる。このような画像処理の先行技術としては、例えば特許文献1がある。特許文献1では、視覚特性を利用したRetinex理論を用いた画像の補正処理を開示している。
特開2005-4506公報
上述したNEXと撮像時間は比例関係にあるため、NEXを増加させると撮像時間が長大化する課題がある。撮像時間を短縮するには、撮像で得られた観察データに対するポスト処理においてノイズを低減する必要があるが、ノイズを低減すると本来の信号であるエッジ(輪郭)や輝度差(コントラスト)も低減してしまう。特許文献1では、Retinex理論を用いた補正処理として、複数のボケフィルタを利用することにより、適応的な補正を実現しているが、処理量が多く、処理に時間がかかる課題がある。
本発明の目的は、上記の課題を解決し、加算回数が少ないMRI画像で画質の維持と処理時間の低減を可能とする画像診断装置、及び画像取得方法を提供することにある。
上記目的を達成するため、本発明においては、画像診断装置であって、画像に変換された観測データのノイズを低減するノイズ低減部と、ノイズ低減部で得られたノイズ低減データに対し、人間の視覚特性を利用した補正処理を行う画像補正部とを備えた構成の画像診断装置を提供する。
また、上記の目的を達成するため、本発明においては、画像診断装置であって、画像に変換された観測データのノイズを低減するノイズ低減部と、ノイズ低減部で得られたノイズ低減データから大局的輝度成分と局所的変動成分を分離する分離部と、観測データ、ノイズ低減データ、大局的輝度成分を用いて補正度を算出する補正度算出部と、観測データ、局所的変動成分、並びに補正度を用いて補正を行う補正処理部と、観測データ取得時のパラメータを用いて、ノイズ低減部、補正度算出部、及び補正処理部を制御する補正制御部とを備える構成の画像診断装置を提供する。
更に、上記の目的を達成するため、本発明においては、画像診断装置の画像取得手法であって、取得された観測データからノイズを低減するノイズ低減処理と、ノイズ低減処理で得られたノイズ低減データから、大局的輝度成分と局所的変動成分を分離する分離処理と、観測データ、ノイズ低減データ、および大局的輝度成分を用いて補正度を算出する補正度算出処理と、観測データ、局所的変動成分、および補正度を用いて補正処理を行う画像取得手法を提供する。
本発明によれば、画像診断装置における撮影時間を短縮しつつ、高画質な画像を取得することが可能となる。
各実施例に係るMRI装置の全体構成の一例を示す図である。 実施例1における画像処理部の構成の一例を示す図である。 実施例1におけるノイズ低減処理の処理フローの一例を示す図である。 実施例1における画像補正部の構成の一例を示す図である。 実施例1における補正度マップの作成手法を説明するための図である。 実施例1における入力部の提示例を示す図である。 実施例1における表示部の表示例を示す図である。 実施例2における画像処理部の構成の一例を示す図である。 実施例3における超音波診断装置の全体構成の一例を示すブロック図である。 実施例4におけるCT装置の全体構成の一例を示すブロック図である。
以下、本発明の種々の実施形態を図面に従い説明する。
実施例1は、画像に変換された観測データのノイズを低減するノイズ低減部と、ノイズ低減部で得られたノイズ低減データに対し、人間の視覚特性を利用した補正処理を行う画像補正部とを備えた構成の画像診断装置、及び画像診断装置の画像取得手法であって、取得された観測データからノイズを低減するノイズ低減処理と、ノイズ低減処理で得られたノイズ低減データから、大局的輝度成分と局所的変動成分を分離する分離処理と、観測データ、ノイズ低減データ、および大局的輝度成分を用いて補正度を算出する補正度算出処理と、観測データ、局所的変動成分、および補正度を用いて補正処理を行う画像取得手法の実施例である。
図1は、実施例1に係るMRI装置の全体構成の一例を示す図である。図1に明らかなように、医用画像診断装置であるMRI装置は大きく分けて、被検体を観測し、観測データを出力する観測部100、及び観測データに基づき被検体の画像を再構成する再構成部106から構成される。
観測部100は、静磁場発生系102、傾斜磁場発生系103、送信系104、受信系105、シーケンサ107、中央処理装置(CPU)108から構成される。静磁場発生系102は、被検体101の周囲の空間に、均一な磁場を発生させる。均一な磁場発生のため、静磁場発生系102には、永久磁石または常電導方式もしくは超電導方式の磁場発生手段が配置される。傾斜磁場発生系103は、傾斜磁場コイル109、および傾斜磁場コイル109を駆動する傾斜磁場電源110から構成され、傾斜磁場を被検体101に印加する。
シーケンサ107は、高周波磁場パルス(RFパルス)と傾斜磁場パルスをある所定のパルスシーケンスで繰り返し印加する制御手段であり、CPU108の制御で動作し、被検体101の断層画像のデータ収集に必要な種々の命令を送信系104、傾斜磁場発生系103、受信系105へ送信する。送信系104は、高周波発生器111と変調器112、増幅器113、高周波コイル114aから構成され、被検体101を構成する原子の原子核スピンに核磁気共鳴を起こさせるRFパルスを照射する。受信系105は、高周波コイル114b、増幅器115、直交位相検波器116、アナログデジタル(A/D)変換器117から構成され、原子核スピンの核磁気共鳴により放出されるエコー信号を受信し、観測データを再構成部106に送信する。
再構成部106は、画像処理部118と、キーボード、マウス、タッチパネル、ボタン等を含む入力部119と、ディスプレイ、プリンタ等を含む出力部120と、磁気ディスク、光ディスク等を含み、データや所望のプログラムを記憶する記憶装置121で構成される。画像処理部118は受信系105から観測データが入力されると、画像を再構成し、出力部120により表示させるとともに、記憶装置121に記録する。図1に示すように、この画像処理部118はCPU108によるプログラム処理で実現することが可能であるが、再構成部106に、CPU108とは別の中央処理装置(CPU)を設置して利用することや、画像処理用の専用ハードウェアで構成して処理することも可能である。
次に図1に示したMRI装置の画像処理部118の処理内容について説明する。図2は、本実施例のMRI装置における画像処理部118で行われる処理の一例を示す機能ブロック図である。本実施例の画像処理部118は、観測部100によって観測された観測データをフーリエ変換により画像に変換する変換部200と、画像に変換された観測データのノイズを低減するノイズ低減部201と、人間の視覚特性を利用した補正処理を行う画像補正部202とから構成される。人間の視覚特性を利用した補正処理は、例えば、Retinex理論などを利用した補正処理であり、エッジ強調、コントラスト変更等が含まれる。人間の視覚特性を利用した補正処理については後述する。
ノイズ低減部201では、観測データを変換して得られた画像のノイズを低減する。ノイズ低減処理には重みつき平滑化フィルタ、バイラテラルフィルタ、ノンローカルミーンフィルタ、ガイデッドフィルタ等の手法がある。本実施例では画像のスパース性(sparsity)を利用したノイズ低減処理について説明する。画像のスパース性、すなわち疎性とは、画像もしくは任意の基底関数により変換された係数にゼロ成分が多く含まれる性質を示す。
図3は本実施例のノイズ低減部201の処理の一例を示す図である。図3に示すようにノイズ低減部201は、ステップST300において初期化を行った後、ステップST301〜ST305までを繰返し実行することにより、画像に変換された観測データに含まれるノイズを低減する。
以下、k+1回目の繰り返しについて説明する。ステップST301において式(1)が計算され、推定結果uk+1が算出される。
Figure 0006744764
ここで、fkは直前(k回目)の繰り返しにより更新された画像の周波数成分、Φはフーリエ変換、ΦTはΦの逆変換を示す。INはすべての要素が1であり、fkと同サイズの配列である。またuC k、u k 、bC k、b kは直前(k回目)の繰り返しにおいて算出された変更成分である。またμはパラメータとしての正の定数である。
次に、ステップST402において、uC k+1、u k+1を式(2)および(3)により算出する。
Figure 0006744764
Figure 0006744764
ここで、Ψ T、Ψ TはそれぞれCurvelet逆変換およびWavelet逆変換である。ここではCurvelet変換およびWavelet変換を利用しているが、これ以外にも、TV(Total Variation)、Ridgelet変換等を用いてもかまわない。また、これらを組み合わせて用いてもよい。
SCおよびSはSoft Shrinkageと呼ばれる処理を示す。SCおよびSはすべての要素についてそれぞれ式(4)、(5)で示す処理を行う。ここでΨ、ΨはそれぞれCurvelet変換およびWavelet変換である。またλはパラメータとしての定数である。
Figure 0006744764
Figure 0006744764
次にステップST303において、bC k+1、b k+1が式(6)、(7)を用いて算出される。
Figure 0006744764
Figure 0006744764
次にステップST304において終了判定がなされる。終了判定としては変更成分u、bの最大値、合計値等が規定値を下回った場合、もしくは繰返し回数が一定回数に達した場合等により終了を判定する。以上のステップを繰返し実行することにより、変換部200から入力された観測データである画像のノイズが低減される。このように、本実施例のノイズ低減部201において、観測データのスパース性を利用してノイズを低減することが可能となる。
次に本実施例の画像補正部202について説明する。画像補正部202は人間の視覚特性を利用した補正処理、すなわち人間の視覚特性に基づいてMRI画像の補正を行う。本明細書において、人間の視覚特性を利用した補正処理とは、例えば、人間の目が照明光に関係なく色や明るさを感じられるとする色恒常性や明るさ恒常性を持つとするRetinex理論などを利用した補正処理を言う。
図4は本実施例の画像補正部202の処理の一例を示す図である。図4に示すように画像補正部202は、大局的輝度成分と局所的変動成分を分離する分離部400と、入力画像の各位置における補正度を算出する補正度算出部401と、算出された各成分および補正度を用いて画像を補正する画像補正部402から構成される。
分離部400、補正度算出部401、画像補正部402各々の処理について説明する。分離部400において、式(8)により大局的輝度成分Lを算出する。
Figure 0006744764
式(8)において、I(x,y)はノイズ低減部201から入力されるノイズ低減画像、G(x,y)はボケ関数、*は畳み込み演算である。
ここでボケ関数Gには移動平均フィルタ、ガウシアンフィルタ、バイラテラルフィルタ等が利用できる。式(8)においては畳み込み演算として記述しているが、フーリエ変換を利用し、積和演算として畳み込み処理も用いても構わない。
次に、式(9)をもちいて、局所的変動成分Rを算出する。
Figure 0006744764
次に補正度算出部401について説明する。補正度算出部401では、入力画像の各位置における補正度を算出し、補正度マップCを作成する。補正度の算出には、輝度、勾配、入力画像とノイズ低減画像の差分値、MRI撮像パラメータ等のうちいずれかひとつまたはそれらの組み合わせを評価指標として用いることができる。補正度算出部401は、補正度の算出に利用する評価指標について、それぞれ補正度マップを作成する。
図5は、補正度算出部401が作成した補正度マップCの概要を説明するための模式図である。ここでは、評価指標として輝度を用いる例を示す。図5において、画像処理部118の変換部200で観測データを変換した観測画像500に基づき、ノイズ低減部201で算出されたノイズ低減画像501が画像補正部202の分離部400に入力される。分離部400により大局的輝度成分と局所的変動成分に分離され、大局的輝度成分502が算出される。この大局的輝度成分502が入力される補正度算出部401は、たとえば低輝度領域におけるノイズの抑制のために、大局的輝度成分502において既定された閾値を下回る領域については補正度を0、上回る領域については補正度を1とするような輝度を評価指標とする補正度マップ503を作成する。また、th1〜th2といった輝度の範囲を設定することで、輝度範囲補正度マップ504が算出される。なお、図5においては、補正度を0(黒)もしくは1(白)の2値としているが、どのような値であってもかまわない。
評価指標として勾配を用いる場合は、たとえば差分法、ソーベルフィルタ等により勾配を算出し、勾配の大きな領域においては大きな補正度を設定することにより、勾配補正度マップ505が算出される。また、評価指標として、入力される観測画像500とノイズ低減画像501の差分値を算出し、ノイズとして低減された領域を推定し補正度マップ506を算出することもできる。このように本実施例の補正度算出部401は1枚もしくはn枚の補正度マップ503〜506を算出する。複数の補正度マップを算出することにより、コントラスト補正やエッジ強調を行う領域および各補正効果の強弱を制御可能となる。すなわち、本実施例の装置においては、補正度算出部401において、複数の評価指標に基づき複数の補正度を算出することにより、より柔軟な補正処理を行うことができる。
次に補正処理部402について説明する。補正処理部402では、分離部400で分離された局所的変動成分Rや補正度算出部401により算出された補正度マップCを用いて、式(10)によりノイズ低減画像を補正する。
Figure 0006744764
ここで、O(x,y)は補正度算出部401で補正された補正画像、Ci(x,y)は算出された補正度マップ、αは補正の割合を制御する正のパラメータである。
式(10)により各画素位置における最終的な補正度を算出し、入力画像を補正する。ここでは補正度マップを乗算しているが、加算することで代用してもかまわない。補正度が0〜1の範囲で規格化されている場合においては乗算を用いることによって最終的な補正度が0〜1の範囲内に収まるため、より簡易に補正が可能である。また、加算、もしくは加算と乗算を組み合わせた多項式とする構成とすることにより、より柔軟な補正度の算出が可能である。すなわち、本実施例の補正処理部402において、補正度算出部401で算出した複数の補正度を乗算あるいは加算もしくは加算と乗算と組み合わせて補正処理を行うことにより、柔軟な対応を可能とすることができる。
このように、本実施例の画像補正部202は、変換部で画像に変換された観測データ、ノイズ低減部から出力されるノイズ低減データ、並びにノイズ低減データ中の大局的輝度成分を用いて補正度を算出する補正度算出部と、観測データ、ノイズ低減データ中の局所的変動成分、並びに補正度を用いて補正を行う補正処理部を含む構成により、より柔軟な補正度の算出が可能である。
本実施例の構成においては、再構成部106の入力部119を用いて、ノイズ低減処理および補正処理におけるパラメータをユーザに設定させる。図6は入力部119におけるユーザへの提示画面の一例である。図6に示すように、数値による入力600、ラジオボタンによる入力601、スライダーによる入力602等が考えられる。このときパラメータのセットをいくつか用意しておくことが望ましい。また、被検体の撮像部位や撮像モードごとにユーザが選択できるようにしても良い。
出力部120はディスプレイ等の出力装置を備え、図2に示した各機能ブロックで得られる観測画像、ノイズ低減画像、補正画像を表示する。図7は出力部120のディスプレイ等の出力装置における表示例である。図7に示すように、出力装置700上に、観測画像701、ノイズ低減画像702、補正画像703のほかに、補正度マップの組み合わせや、補正処理の異なる複数の補正結果704〜706を同時に表示できるようにし、ユーザが選択できるようにしても良い。また、処理に用いた補正度マップの種類、パラメータ等を表示させる領域707を表示しても良い。
以上、本実施例の構成によれば、人間の視覚特性を利用した補正処理により、臨床に耐えうる高画質なMRI画像を高速に撮像可能であり、加算回数の低減に基づく撮像時間の短縮による患者負担の軽減と、診断効率の向上が可能である。
実施例2は、MRI装置の各種パラメータを利用して、ノイズ除去および補正を制御することで、ユーザが簡便に好ましい画像を取得することが可能な実施例である。すなわち、本実施例は、画像診断装置であって、画像に変換された観測データのノイズを低減するノイズ低減部と、ノイズ低減部で得られたノイズ低減データから大局的輝度成分と局所的変動成分を分離する分離部と、観測データ、ノイズ低減データ、大局的輝度成分を用いて補正度を算出する補正度算出部と、記観測データ、局所的変動成分、並びに補正度を用いて補正を行う補正処理部と、観測データ取得時のパラメータを用いて、ノイズ低減部、補正度算出部、及び補正処理部を制御する補正制御部とを備える構成の画像診断装置及び画像取得方法の実施例である。
本実施例においてもMRI装置の全体構成は実施例1同様、図1で表される。本実施例が実施例1と異なる点は、人間の視覚特性を利用した補正処理を行うためのノイズ低減部201ならびに画像補正部202において、観測部100の観測パラメータを利用することにある。
図8は実施例2における要部構成の一例を示すブロック図である。実施例2においては、画像処理部118にノイズ低減部201および画像補正部202を制御する補正制御部700を追加する。
補正制御部800について説明する。補正制御部800は、観測部100から観測時に用いたパラメータ情報が入力される。観測時に用いたパラメータ情報とは、たとえば、MRI装置のTR(Repetition Time)、フリップ角、スライス厚、FOV(Field of View)、マトリクスサイズ、複数回撮像し加算する信号の加算回数、バンド幅、部位情報等である。MRI装置においてはこれらのパラメータによりSNRが変動するため、これらのパラメータを利用して適切なノイズ低減および補正を行うことで、ユーザはより好適な画像を簡便に取得可能となる。例えば、加算回数が少ない場合、SNRが低下することが予想されるため、ノイズ低減効果および補正効果を大きくすることが必要となる。補正制御部800は、上記のパラメータから適切なノイズ低減および補正パラメータを算出し、ノイズ低減部201および画像補正部202にそれぞれ送信する。
まずノイズ低減部201に送信されるパラメータは、式(1)で利用されるパラメータμ、および式(4)、(5)で使用されるパラメータλである。一般的に、SNRが低下するのは、TRが短い場合、フリップ角が小さい場合、スライス圧が薄い場合、FOVが小さい場合、マトリクスサイズが大きい場合、加算回数が少ない場合、バンド幅が広い場合である。実際にはこれらの組み合わせを用いて最終的な補正パラメータを決定することが望ましい。また、これらのパラメータに対する補正パラメータのテーブルを用意し、CPU108の記憶部や、再構成部106の記憶装置121に記憶しておき利用しても良い。SNRの低下が予想される場合、補正制御部800はパラメータλおよびパラメータμを大きくすることでノイズ低減効果を向上させることができる。
一方、補正制御部800から画像補正部202に送信されるパラメータは、補正度マップ算出用パラメータと、補正処理パラメータの2つである。補正度マップ算出用パラメータとしては、たとえば閾値、輝度の幅および、それぞれについての補正度情報等が送信される。すなわち、MRI撮影におけるモード情報、部位情報を用いて、臨床上必要な輝度範囲を事前に定義し、定義された補正パラメータを算出することによってより好適な補正を行う。たとえばMRI画像のうち、観察対象の輝度が概ねth1〜th2の範囲にある場合、th1およびth2が補正度算出部401に送信され、補正度マップが作成される。作成される補正度マップとしては、th1およびth2を用いた輝度範囲補正度マップや、th1およびth2付近の勾配を用いた勾配補正度マップが挙げられる。
また、補正処理パラメータは、式(10)で用いられるパラメータαである。補正処理パラメータαは画像補正部202の補正処理部402に送信される。前述したようにSNRが低下する場合、ノイズ低減効果が大きくなる。これにより必要なエッジ情報がボケてしまうことがある。これに対応し、補正制御部800は補正処理パラメータαを大きくすることにより、補正がより強く実行され、好適な画像を取得することを可能とする。
以上、本実施例によれば、MRI装置の各種パラメータを利用してユーザはより簡便に好適なMRI画像を取得可能である。
実施例3は、画像診断装置が、超音波診断装置の場合の実施例である。
図9は、医用画像診断装置である超音波診断装置の一構成例を示すブロック図である。図9において、図1に示した構成要素と同一の要素に関しては、同一の符号を付してその説明を省略する。同図において、超音波診断装置を構成する超音波観測部900は、超音波探触子901と、送信部902と、受信部903と、超音波送受信制御部904と、整相加算部905から構成される。
送信部902は超音波探触子901を介して被検体101に時間間隔をおいて超音波を繰り返し送信する。受信部903は被検体101から発生する時系列の反射エコー信号を受信する。超音波送受信制御部904は送信部902と受信部903を制御する。整相加算部905は受信された反射エコー信号を整相加算して、RF信号のフレームデータを時系列に生成する。整相加算部905はアナログデジタル(A/D)変換器を内蔵し、RF信号フレームデータを観測データとして、再構成部106の画像処理部118に出力され、画像処理部118はRFフレームデータからなる観測データを用いて超音波エコー画像を生成する。
本実施例の超音波観測部900は、観測データを画像処理部118に送信する。本実施例における画像処理部118も第1および第2の実施例同様、図2で表される。ただし、本実施例において、画像処理部118の変換部200はRFフレームデータを画像に変換する。変換された画像に対し、実施例1または実施例2に示した人間の視覚特性を利用した補正処理補正処理を行い、高画質な超音波画像を算出することができる。
本実施例の超音波診断装置によれば、高画質な超音波画像を取得することが可能となる。
実施例4は、画像診断装置が、高画質なCT(Computed Tomography)画像を取得することが可能なCT装置の実施例である。
図10は、医用画像診断装置としてのCT(Computed Tomography)装置の一構成例を示すブロック図である。同図においても、図1に示した構成要素と同一の要素に関しては、同一の符号を付してその説明を省略する。
図10において、CT観測部1000は、X線管装置1001と、回転円盤1002と、コリメータ1003と、X線検出器1006と、データ収集装置1007と、寝台1005と、システム制御部1008から構成される。X線管装置1001は、寝台1005上に載置された被検体にX線を照射する装置である。コリメータ1003はX線管装置1001から照射されるX線の放射範囲を制限する装置である。回転円盤1002は、寝台1005上に載置された被検体が入る開口部1004を備えるとともに、X線管装置1001とX線検出器1006を搭載し、被検体の周囲を回転するものである。
X線検出器1006は、X線管装置1001と対向配置され被検体を透過したX線を検出することにより透過X線の空間的な分布を計測する装置であり、多数のX線検出素子を回転円盤1002の回転方向に配列したもの、若しくは回転円盤1002の回転方向と回転軸方向との2次元に配列したものである。データ収集装置1007は、X線検出器1006で検出されたX線量をデジタルデータとして収集する装置である。また、システム制御部1008は、回転円盤1002の回転、寝台1005の上下前後左右動、X線管装置1001に入力される電力等を制御するものである。
本実施例のCT装置のCT観測部1000は、観測データを再構成部106の画像処理部118に送信する。本実施例における画像処理部118も第1および第2の実施例同様、図2で表される。ただし、本実施例において、変換部200は、観測データをフィルタ補正逆投影法、逐次近似画像再構成法等の処理によって画像に変換する。変換された画像に対し、実施例1または実施例2に示した人間の視覚特性を利用した補正処理を行い、高画質なCT画像を算出する。
以上のように、本実施例によれば、高画質なCT画像を取得することが可能となる。更に、従来よりもX線量を低減することが可能となり、被ばく量の低減が期待できる。
なお、本発明は上記した実施例に限定されるものではなく、様々な変形例が含まれる。例えば、上記した実施例は本発明のより良い理解のために詳細に説明したのであり、必ずしも説明の全ての構成を備えるものに限定されものではない。また、ある実施例の構成の一部を他の実施例の構成に置き換えることが可能であり、また、ある実施例の構成に他の実施例の構成を加えることが可能である。また、各実施例の構成の一部について、他の構成の追加・削除・置換をすることが可能である。
更に、上述した各構成、機能、処理部等は、それらの一部又は全部を実現するCPUのプログラムを作成する例を説明したが、それらの一部又は全部を例えば集積回路で設計する等によりハードウェアで実現しても良いことは言うまでもない。
100 観測部
101 被検体
102 静磁場発生系
103 傾斜磁場発生系
104 送信系
105 受信系
106 再構成部
107 シーケンサ
108 中央処理装置(CPU)
109 傾斜磁場コイル
110 傾斜磁場電源
111 高周波発生器
112 変調器
113、115 増幅器
114 高周波コイル
116 直交位相検波器
117 アナログデジタル(A/D)変換器
118 画像処理部
119 入力部
120 出力部
121 記憶装置
200 変換部
201 ノイズ低減部
202 画像補正部
400 分離部
401 補正度算出部
402 補正処理部
500 観測画像
501 ノイズ低減画像
502 大局的輝度成分
503〜505 補正度マップ
600 数値による入力
601 ラジオボタンによる入力
602 スライダーによる入力
700 表示装置
701 観測画像
702 ノイズ低減画像
703 補正画像
704〜706 ノイズ低減画像または補正画像
707 表示領域
800 補正制御部
900 超音波観測部
901 探触子
902 送信部
903 受信部
904 超音波送受信制御部
905 整相加算部
1000 CT観測部
1001 X線管装置
1002 回転円盤
1003 コリメータ
1004 開口部
1005 寝台
1006 X線検出器
1007 データ収集装置
1008 システム制御部

Claims (14)

  1. 画像診断装置であって、
    画像に変換された観測データのノイズを低減するノイズ低減部と、
    前記ノイズ低減部で得られたノイズ低減データに対し、人間の視覚特性を利用した補正処理を行う画像補正部とを備え
    前記画像補正部は、
    前記観測データ、前記ノイズ低減データ、並びに前記ノイズ低減データ中の大局的輝度成分を用いて補正度を算出する補正度算出部を含む
    ことを特徴とする画像診断装置。
  2. 請求項に記載の画像診断装置であって、
    前記画像補正部は、
    前記観測データ、前記ノイズ低減データ中の局所的変動成分、並びに前記補正度を用いて補正を行う補正処理部を更に含む
    ことを特徴とする画像診断装置。
  3. 請求項に記載の画像診断装置であって、
    前記補正度算出部は、
    複数の評価指標に基づき複数の前記補正度を算出する
    ことを特徴とする画像診断装置。
  4. 請求項に記載の画像診断装置であって、
    前記補正処理部は、
    複数の前記補正度を積算し補正処理を行う
    ことを特徴とする画像診断装置。
  5. 請求項に記載の画像診断装置であって、
    前記ノイズ低減部において、前記観測データのスパース性を利用してノイズを低減することを特徴とする画像診断装置。
  6. 請求項1記載の画像診断装置であって、
    前記観測データは核磁気共鳴イメージング(以下、MRIと称す)装置で得られたデータである
    ことを特徴とする画像診断装置。
  7. 画像診断装置であって、
    画像に変換された観測データのノイズを低減するノイズ低減部と、
    前記ノイズ低減部で得られたノイズ低減データから大局的輝度成分と局所的変動成分を分離する分離部と、
    前記観測データ、前記ノイズ低減データ、前記大局的輝度成分を用いて補正度を算出する補正度算出部と、
    前記観測データ、前記局所的変動成分、並びに前記補正度を用いて補正を行う補正処理部と、
    前記観測データ取得時のパラメータを用いて、前記ノイズ低減部、前記補正度算出部、及び前記補正処理部を制御する補正制御部と、を備える
    ことを特徴とする画像診断装置。
  8. 請求項に記載の画像診断装置であって、
    前記補正度算出部は、
    複数の評価指標に基づき複数の前記補正度を算出する
    ことを特徴とする画像診断装置。
  9. 請求項に記載の画像診断装置であって、
    前記補正処理部は、
    複数の前記補正度を積算し補正処理を行う
    ことを特徴とする画像診断装置。
  10. 請求項に記載の画像診断装置であって、
    前記ノイズ低減部において、前記観測データのスパース性を利用してノイズを低減することを特徴とする画像診断装置。
  11. 請求項記載の画像診断装置であって、
    前記観測データはMRI装置で得られたデータである
    ことを特徴とする画像診断装置。
  12. 請求項11に記載の画像診断装置であって、
    前記補正制御部で用いるパラメータとして、前記MRI装置で複数回撮像して加算する加算回数を利用する
    ことを特徴とする画像診断装置。
  13. 画像診断装置の画像取得手法であって、
    取得された観測データからノイズを低減するノイズ低減処理と、
    前記ノイズ低減処理で得られたノイズ低減データから、大局的輝度成分と局所的変動成分を分離する分離処理と、
    前記観測データ、前記ノイズ低減データ、および前記大局的輝度成分を用いて補正度を算出する補正度算出処理と、
    前記観測データ、前記局所的変動成分、および前記補正度を用いて補正処理を行う
    ことを特徴とする画像取得手法。
  14. 請求項13に記載の画像取得手法であって、
    前記観測データ取得時のパラメータを用いて、前記ノイズ低減処理、前記補正度算出処理、前記補正処理を制御する制御処理を更に含む
    ことを特徴とする、画像取得手法。
JP2016105023A 2016-05-26 2016-05-26 画像診断装置、及び画像取得方法 Active JP6744764B2 (ja)

Priority Applications (4)

Application Number Priority Date Filing Date Title
JP2016105023A JP6744764B2 (ja) 2016-05-26 2016-05-26 画像診断装置、及び画像取得方法
PCT/JP2017/014966 WO2017203875A1 (ja) 2016-05-26 2017-04-12 画像診断装置、及び画像取得方法
US16/093,680 US10732245B2 (en) 2016-05-26 2017-04-12 Diagnostic imaging device and image acquisition method
CN201780031293.1A CN109152549B (zh) 2016-05-26 2017-04-12 图像诊断装置以及图像取得方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2016105023A JP6744764B2 (ja) 2016-05-26 2016-05-26 画像診断装置、及び画像取得方法

Publications (3)

Publication Number Publication Date
JP2017209329A JP2017209329A (ja) 2017-11-30
JP2017209329A5 JP2017209329A5 (ja) 2018-12-20
JP6744764B2 true JP6744764B2 (ja) 2020-08-19

Family

ID=60412333

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2016105023A Active JP6744764B2 (ja) 2016-05-26 2016-05-26 画像診断装置、及び画像取得方法

Country Status (4)

Country Link
US (1) US10732245B2 (ja)
JP (1) JP6744764B2 (ja)
CN (1) CN109152549B (ja)
WO (1) WO2017203875A1 (ja)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6896177B2 (ja) * 2018-08-10 2021-06-30 三菱電機株式会社 画像処理装置及び画像処理方法
JP7186604B2 (ja) * 2018-12-25 2022-12-09 キヤノンメディカルシステムズ株式会社 医用画像処理装置、医用画像処理方法、およびプログラム
JP7455508B2 (ja) 2018-12-26 2024-03-26 キヤノンメディカルシステムズ株式会社 磁気共鳴イメージング装置および医用複素数画像処理装置
JP7297485B2 (ja) * 2019-03-25 2023-06-26 キヤノンメディカルシステムズ株式会社 超音波診断装置、医用画像処理装置および医用画像処理プログラム
DE102019217220A1 (de) * 2019-11-07 2021-05-12 Siemens Healthcare Gmbh Computerimplementiertes Verfahren zur Bereitstellung eines Ausgangsdatensatzes
CN113469919A (zh) * 2021-07-27 2021-10-01 深圳市赛禾医疗技术有限公司 超声图像的处理方法、装置及电子设备

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4345366B2 (ja) * 2003-06-12 2009-10-14 コニカミノルタホールディングス株式会社 画像処理プログラムおよび画像処理装置
DE10351928A1 (de) * 2003-11-07 2005-06-16 MAX-PLANCK-Gesellschaft zur Förderung der Wissenschaften e.V. Messung am visuellen System eines Probanden
EP2533688B1 (en) * 2010-02-08 2013-11-20 Koninklijke Philips N.V. Apparatus and method for influencing and/or detecting magnetic particles in a field of view having an array of single-sided transmit coil sets
US8787638B2 (en) * 2011-04-07 2014-07-22 The Chinese University Of Hong Kong Method and device for retinal image analysis
JP6009180B2 (ja) * 2012-03-07 2016-10-19 日立マクセル株式会社 医療用画像表示システム
JP6071444B2 (ja) * 2012-11-07 2017-02-01 キヤノン株式会社 画像処理装置及びその作動方法、プログラム
CN105125176B (zh) * 2015-09-23 2018-07-10 博联众科(武汉)科技有限公司 一种静脉血管导航装置和方法

Also Published As

Publication number Publication date
US10732245B2 (en) 2020-08-04
US20190137589A1 (en) 2019-05-09
CN109152549B (zh) 2022-08-02
CN109152549A (zh) 2019-01-04
WO2017203875A1 (ja) 2017-11-30
JP2017209329A (ja) 2017-11-30

Similar Documents

Publication Publication Date Title
JP6744764B2 (ja) 画像診断装置、及び画像取得方法
KR101604812B1 (ko) 의료 영상 처리 장치 및 그에 따른 의료 영상 처리 방법
JP4662852B2 (ja) 画像処理方法、画像処理装置、医用画像診断支援システム、及び時間軸方向フィルタリング方法
US11419555B2 (en) Image capturing apparatus and method
JP4820582B2 (ja) ヘリカルマルチスライスctのための回復ノイズを伴うヘリカルウィンドミルアーチファクトを低減する方法
JP5248010B2 (ja) データ補正装置、データ補正方法、磁気共鳴イメージング装置およびx線ct装置
JP5349673B2 (ja) データ補正装置
JP6753798B2 (ja) 医用撮像装置、画像処理方法及びプログラム
US10429466B2 (en) MRI apparatus, image processing apparatus, and image processing method
US20160012569A1 (en) Image processing device, magnetic resonance imaging apparatus and image processing method
JP7122301B2 (ja) 画像データ内のノイズ低減
EP3349655B1 (en) Tomography apparatus and controlling method for the same
KR101525014B1 (ko) 자기 공명 영상 장치 및 그에 따른 자기 공명 영상의 이미징 방법
KR20160120963A (ko) 단층 촬영 장치 및 그에 따른 단층 영상 복원 방법
CN112120722A (zh) X射线断层合成装置、图像处理装置以及计算机可读记录介质
KR101681313B1 (ko) 의료 영상 제공 장치 및 그에 따른 의료 영상 제공 방법
JP6377712B2 (ja) 超音波診断装置、及び画像撮像方法
KR101601011B1 (ko) 의료 영상 장치에서의 혈류 속도 산출 방법 및 그 의료 영상 장치
US20220019850A1 (en) Medical data processing apparatus and method
CN116630213A (zh) 4d-mri图像的获取方法、装置及存储介质

Legal Events

Date Code Title Description
A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20181106

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20181106

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20191210

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20200203

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: 20200707

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20200731

R150 Certificate of patent or registration of utility model

Ref document number: 6744764

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

S111 Request for change of ownership or part of ownership

Free format text: JAPANESE INTERMEDIATE CODE: R313111

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250