JPS63174642A - Mri apparatus - Google Patents

Mri apparatus

Info

Publication number
JPS63174642A
JPS63174642A JP62005580A JP558087A JPS63174642A JP S63174642 A JPS63174642 A JP S63174642A JP 62005580 A JP62005580 A JP 62005580A JP 558087 A JP558087 A JP 558087A JP S63174642 A JPS63174642 A JP S63174642A
Authority
JP
Japan
Prior art keywords
data
dimensional
distance
filter
reference point
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
JP62005580A
Other languages
Japanese (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.)
Shimadzu Corp
Original Assignee
Shimadzu Corp
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 Shimadzu Corp filed Critical Shimadzu Corp
Priority to JP62005580A priority Critical patent/JPS63174642A/en
Publication of JPS63174642A publication Critical patent/JPS63174642A/en
Pending legal-status Critical Current

Links

Abstract

(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。
(57) [Summary] This bulletin contains application data before electronic filing, so abstract data is not recorded.

Description

【発明の詳細な説明】[Detailed description of the invention] 【産業上の利用分野】[Industrial application field]

この発明は、核磁気共鳴現象(NMR現象)を利用して
イメージングを行うMHI装置に関し、特にフーリエ変
換を行って画像を再構成するMRI装置に関する。
The present invention relates to an MHI apparatus that performs imaging using a nuclear magnetic resonance phenomenon (NMR phenomenon), and more particularly to an MRI apparatus that performs Fourier transform to reconstruct an image.

【従来の技術】[Conventional technology]

収集されたデータに対し2次元(あるいは3次元)のフ
ーリエ変換を行って画像(2次元または3次元画像)を
再構成する場合、フーリエ変換はその性質からサンプリ
ングしたデータ列を周期関数として取り扱う。ところが
MRI装置によって実際に収集したデータは元来周期関
数でない上、比較的少ないポイント数でサンプリングが
行われることが多い。このような場合には、打ち切り幅
の両端で不連続な点が現れることになり、フーリエ変換
の結果つまり周波数領域に不要な周波数成分がつけ加わ
る。これはリーケイジ(Lea、kage)と呼ばれ、
これを抑えるために、一般に標本点列に対するフィルタ
リングが行われている。このフィルタとして代表的なも
のにHannigのフィルタが知られている。
When reconstructing an image (two-dimensional or three-dimensional image) by performing two-dimensional (or three-dimensional) Fourier transform on collected data, the Fourier transform treats the sampled data string as a periodic function due to its nature. However, data actually collected by an MRI apparatus is not originally a periodic function, and sampling is often performed using a relatively small number of points. In such a case, discontinuous points will appear at both ends of the truncation width, and unnecessary frequency components will be added to the Fourier transform result, that is, the frequency domain. This is called Lea, kage.
In order to suppress this, filtering is generally performed on the sample point sequence. A Hannig filter is known as a typical example of this filter.

【発明が解決しようとする問題点】[Problems to be solved by the invention]

しかし、MR画偶は2次元あるいは3次元のフーリエ変
換によって再構成されるので、フィルタリングも2次元
的にあるいは3次元的に行う必要があって、そのための
フィルタとしての情報を2次元(あるいは3次元)の形
に作成し、記憶しておかなければならない。2次元のフ
ィルタ情報の作成には多くの時間を要し、またこれを記
憶させるにも大容量のメモリを必要とする。 そのため、従来では、フィルタ処理を省略したり、ある
いはフィルタ処理を行う場合でも1次元のフィルタにと
どめる方法を採用することが多い。 このように省略したりするのでは、前述のリーケイジを
十分に抑えることができず、このリーケイジのために再
構成画像にアーティファクトを生じる原因ともなってい
る。 この発明は、収集データに対するフィルタリングを2次
元(あるいは3次元)のフィルタ情報を必要とせずに行
え、且つ高速に理想的なフィルタリングを行うことがで
きるよう改善したMRI装置を提供することを目的とす
る。
However, since MR figurines are reconstructed by two-dimensional or three-dimensional Fourier transformation, it is necessary to perform filtering two-dimensionally or three-dimensionally. dimension) and memorize it. It takes a lot of time to create two-dimensional filter information, and a large capacity memory is also required to store it. Therefore, conventionally, a method is often adopted in which filter processing is omitted, or even when filter processing is performed, it is limited to a one-dimensional filter. Omitting it in this way does not allow the aforementioned leakage to be sufficiently suppressed, and this leakage also causes artifacts to be produced in the reconstructed image. An object of the present invention is to provide an improved MRI apparatus that can perform filtering on collected data without requiring two-dimensional (or three-dimensional) filter information and can perform ideal filtering at high speed. do.

【問題点を解決するための手段】[Means to solve the problem]

この発明によるMHI装置は、収集されたデータマトリ
クスにおいて各データ点と所定の基準点との距離を求め
る手段と、基準点からの距離の関数として与えられるフ
ィルタ値がその距離毎に予め記憶させられている記憶手
段と、上記の距離により該当するフィルタ値を上記記憶
手段から読み出し、各データ点のデータに作用させる手
段とを有する。
The MHI device according to the present invention includes means for determining the distance between each data point and a predetermined reference point in the collected data matrix, and a filter value given as a function of the distance from the reference point is stored in advance for each distance. and means for reading a filter value corresponding to the distance from the storage means and applying it to the data of each data point.

【作  用】[For production]

リーケイジを抑えるための2次元(または3次元)のフ
ィルタリングを行うことを検討してみると、データマト
リクス上で基準点から等距離の点にあるデータに対して
は同一のフィルタ値(重み)を作用していることが分か
る。 そこで、このことに着目して、収集したデータマI・リ
クス上での各データと基進点との距離毎にフィルタ値を
記憶しておき、データマトリクス上の各データ点と基準
点との距離を求めて、この距離で上記の記憶していたフ
ィルタ値を読み出し、その点のデータに作用させれば、
必要なフィルタリングが行えたことになる。 こうすると、データマトリクス上での距離を求めること
と、記憶手段からの読み出しと、読み出したフィルタ値
をデータに作用させることだけを行えばよいので、2次
元または3次元のフィルタ情報を作成するのに多くの時
間をかけたり、記憶するのに大容量のメモリを必要とす
ることを避けることができる。そのため、理想的なフィ
ルタリングを高速に行うことができる。 [実 施 例] 第1図において、レジスタ1.2にはホストコンピュー
タあるいはマイクロプロセッサなどから送られた、デー
タマトリクス上のデータ位置を表す列番号iと行番号j
とがセットされる。そして、このi、jによりメモリ3
.4がアクセスされ、メモリ3.4に記憶させられてい
たLx2テーブル、Ly2テーブル上の対応するLx2
、Ly2の各位が読み出される。これらの各位は加算器
5で加算され、その結果得られるL2(=Lx2+Ly
2)がレジスタ6に格納される。つぎにこのL2により
メモリ7の平方根テーブル上の対応するしの値が読み出
され、レジスタ8に送られる。このしによりフィルタ値
く重み)Wが、メモリ9のフィルタテーブルより読み出
され、このWがレジスタ10に格納される。 ここで、MRI装置によって収集されたデータが第2図
に示すように2次元のマトリクス状に配列されるものと
する。このマトリクスでX方向は周波数方向であり、X
方向は位相方向である(実際の空間における方向とは異
なる)。説明を簡単にするためマトリクスは8×8とし
である。理想的なフィルタリングは基準点0から等距離
に位置しているデータに対して同一のフィルタ値(重み
)を作用させる。すなわち、基準点Oから距ILに位置
しているデータに対し、Hanningフィルタの場合
、第3図で示すような、 W = 1/2+ (1/2)CO8((Z L )の
フィルタ値Wを作用させる。 そのため、まず、基準点Oと各データの位置との間の距
離りを求める必要がある。そこで、各データ位置を表す
列番号iと行番号jとから、基準点0からのX方向距離
Lxとy方向距’4 L yとを求めて、それらを2乗
し、 L=’rロア萼工]τ の計算を行う。メモリ3.4には、予めi、jに対応し
て計算されたLi2、Li2の各位が記憶させられてい
るので、これらを読み出し、加算器5で加算した上、平
方根の各位がテーブル化して記憶させられているメモリ
7を読み出すことによってこの計算を行うことができる
。 こうしてデータ位置(i、j)と基準点Oとの距離りが
求められるので、このLの値をメモリ9のフィルタテー
ブルをアクセスするためのアドレス情報として用いる。 すなわち、第3図で示される関数W (L)にLの値を
与えて、そのLに対応するフィルタ値Wを求めることに
相当する操作が行われる。こうしてメモリ9から読み出
されレジスタ10によって保持される値Wはデータ位置
くi、j)に対して重み付けするためのフィルタ値とな
る。したがって、この値をそのデータに対して作用させ
ればフィルタリングを行うことができる。そこで、この
操作を上記のマトリクス状に配列された全てのデータに
対して行うことによってリーケイジを抑えることができ
る。 なお、メモリ3.4.9に記憶させられたLz2、Li
2、Wの各テーブルは任意に書き換えることもできるの
で、この部分を制御するマイクロプロセッサあるいはホ
ストコンピュータ等によって書き換えるようにすれば、
基準点0の位置を変更できるほか、作用させるフィルタ
の形状も上記とは異なる任意のものを自由に選択できる
。 さらに、収集されたデータが3次元的なマトリクスに配
列される場合にも拡張できる。すなわち、第4図に示す
ように、第1図の構成にレジスタ11とメモリ12と加
算器13とを追加するのである。ここで、レジスタ11
にはデータ面の番号kが格納されるものとする。つまり
kは第5図のように2次元にデータが配列された2次元
マトリクス(データ面)が何枚も3次元的に並べられて
いる場合のデータ面の番号である。このように3次元的
なデータマトリクスの場合も理想的なフィルタリングは
、基準点0から等距離に位置しているデータに対して同
一のフィルタ値を作用させることで実現される。すなわ
ち、基準点0を中心とする同一半径の球の表面に位置す
るデータについては同一のフィルタ値を作用させる。デ
ータ面にのある位置Pのデータについては、基準点0が
ら点Pまでの距NLを求め、それに対応するフィルタ値
Wを探し、この位置PのデータにWを作用させればよい
。上記の実施例の場合とはこのLを求める構成だけが異
なり、Lが求められた後は同じである。正確に言えばL
zを求めれば後は同じである。 L2=Lx2+Ly2+Lz2 であって、この計算を行えばよい。Lzは基準点0から
点Pまでの2方向の距離であるから、データ面の番号k
によってメモリ12のアドレスを指定する。このメモリ
12には予め基準点0がら各データ面までの2方向距離
の2乗に対応する値のテーブルが記憶させられている。 そこで、このメモリ12からLz2の値を読み出し、こ
れを加算器13に送り、加算器5からのLx2+Ly2
と加えて、その結果をレジスタ6に保持させれば、上記
の計算が行える。 このように各メモリに平方値、平方根値、フィルタ値を
記憶させ、これを読み出すだけであり、計算としては単
なる加算だけであるから、フィルタ値を2次元、3次元
情報として計算、保持する場合に比較して、その計算量
及び必要とするメモリ量を大幅に削減できる。たとえば
、データマトリクスを256X256とする場合、メモ
リ容量は1750〜1/100程度に少なくできる。3
次元のデータマトリクスの場合は、この改善率は一層大
きく、上記のさらに1710以下となる。また、作用さ
せるフィルタの特性を変更することも容易である。
When considering performing two-dimensional (or three-dimensional) filtering to suppress leakage, it is possible to apply the same filter value (weight) to data at points equidistant from the reference point on the data matrix. I can see that it is working. Therefore, focusing on this, filter values are stored for each distance between each data on the collected data matrix and the reference point, and the distance between each data point on the data matrix and the reference point is , read out the filter value stored above at this distance, and apply it to the data at that point.
This means that the necessary filtering has been performed. In this way, you only need to calculate the distance on the data matrix, read it from the storage means, and apply the read filter value to the data, making it easy to create two-dimensional or three-dimensional filter information. You can avoid spending a lot of time on it or requiring a large amount of memory to store it. Therefore, ideal filtering can be performed at high speed. [Embodiment] In FIG. 1, register 1.2 contains a column number i and a row number j indicating a data position on a data matrix, which are sent from a host computer or microprocessor.
is set. Then, with these i and j, memory 3
.. 4 was accessed and stored in memory 3.4, the corresponding Lx2 on the Ly2 table
, Ly2 are read out. Each of these positions is added by an adder 5, and the resultant L2 (=Lx2+Ly
2) is stored in register 6. Next, the corresponding value on the square root table in the memory 7 is read out by this L2 and sent to the register 8. Accordingly, the filter value (weight) W is read out from the filter table in the memory 9, and this W is stored in the register 10. Here, it is assumed that the data collected by the MRI apparatus are arranged in a two-dimensional matrix as shown in FIG. In this matrix, the X direction is the frequency direction, and
The direction is the phase direction (different from the direction in real space). To simplify the explanation, the matrix is assumed to be 8×8. Ideal filtering applies the same filter value (weight) to data located equidistant from reference point 0. That is, for data located at a distance IL from the reference point O, in the case of a Hanning filter, the filter value W of W = 1/2 + (1/2) CO8 ((Z L ) as shown in FIG. Therefore, it is first necessary to find the distance between the reference point O and the position of each data.Therefore, from the column number i and row number j representing each data position, calculate the distance from the reference point 0. Find the distance Lx in the X direction and the distance Lx in the y direction, square them, and calculate L='rLower Calyx]τ.Memory 3.4 contains data corresponding to i and j in advance. Since each part of Li2 and Li2 calculated by Since the distance between the data position (i, j) and the reference point O is determined in this way, this value of L is used as address information for accessing the filter table in the memory 9. An operation equivalent to giving the value of L to the function W (L) shown in FIG. The value W becomes a filter value for weighting the data position i, j). Therefore, filtering can be performed by applying this value to the data. Therefore, leakage can be suppressed by performing this operation on all the data arranged in the matrix. In addition, Lz2, Li stored in memory 3.4.9
2. Each table of W can be rewritten arbitrarily, so if this part is rewritten by the microprocessor or host computer that controls it,
In addition to being able to change the position of the reference point 0, the shape of the filter to be applied can also be freely selected from the above. Furthermore, it can be extended to cases where the collected data is arranged in a three-dimensional matrix. That is, as shown in FIG. 4, a register 11, a memory 12, and an adder 13 are added to the configuration of FIG. Here, register 11
It is assumed that the data plane number k is stored in . In other words, k is the number of the data plane when a number of two-dimensional matrices (data planes) in which data is arranged two-dimensionally are arranged three-dimensionally as shown in FIG. In this way, even in the case of a three-dimensional data matrix, ideal filtering is achieved by applying the same filter value to data located equidistantly from the reference point 0. That is, the same filter value is applied to data located on the surface of a sphere having the same radius centered on the reference point 0. For data at a position P on the data plane, it is sufficient to find the distance NL from the reference point 0 to the point P, find the filter value W corresponding to it, and apply W to the data at this position P. This embodiment differs from the above embodiment only in the configuration for determining L, but after L is determined, it is the same. To be exact, L
Once you find z, the rest is the same. L2=Lx2+Ly2+Lz2, and this calculation can be performed. Since Lz is the distance in two directions from reference point 0 to point P, the data plane number k
The address of the memory 12 is specified by . This memory 12 stores in advance a table of values corresponding to the square of the distance in two directions from the reference point 0 to each data surface. Therefore, the value of Lz2 is read from this memory 12 and sent to the adder 13, and the value of Lx2+Ly2 from the adder 5 is
In addition, the above calculation can be performed by holding the result in register 6. In this way, the square value, square root value, and filter value are stored in each memory and only read out, and the calculation is just a simple addition, so when the filter value is calculated and stored as two-dimensional or three-dimensional information. The amount of calculation and amount of memory required can be significantly reduced compared to . For example, when the data matrix is 256×256, the memory capacity can be reduced to about 1750 to 1/100. 3
In the case of a dimensional data matrix, this improvement rate is even greater, even below the above-mentioned value of 1710. Further, it is also easy to change the characteristics of the filter to be applied.

【発明の効果】【Effect of the invention】

この発明のMRI装置によれば、データマトリクス上で
基準点から等距離にあるデータには同一のフィルタ値を
作用させればよいことに着目して、2次元あるいは3次
元の理想的なフィルタ処理を高速に行うことができる。 そのため、再構成画像のアーティファクトを軽減するの
に効果的である。
According to the MRI apparatus of the present invention, two-dimensional or three-dimensional ideal filter processing is performed by focusing on the fact that the same filter value should be applied to data equidistant from a reference point on a data matrix. can be done quickly. Therefore, it is effective in reducing artifacts in reconstructed images.

【図面の簡単な説明】[Brief explanation of the drawing]

第1図はこの発明の一実施例のブロック図、第2図は距
離を変数とするフィルタ関数を表すグラフ、第、3図は
2次元データマトリクスを表す図、第4図は他の実施例
のブロック図、第5図は3次元データマトリクスを表す
図である。 1.2.6.8.10.11・・・レジスタ、3.4.
7.9.12・・・メモリ、5.13・・・加算器。
Figure 1 is a block diagram of one embodiment of this invention, Figure 2 is a graph representing a filter function with distance as a variable, Figures 3 and 3 are diagrams representing a two-dimensional data matrix, and Figure 4 is another embodiment. The block diagram of FIG. 5 is a diagram representing a three-dimensional data matrix. 1.2.6.8.10.11...Register, 3.4.
7.9.12...Memory, 5.13...Adder.

Claims (1)

【特許請求の範囲】[Claims] (1)収集されたデータマトリクスにおいて各データ点
と所定の基準点との距離を求める手段と、基準点からの
距離の関数として与えられるフィルタ値がその距離毎に
予め記憶させられている記憶手段と、上記の距離により
該当するフィルタ値を上記記憶手段から読み出し、各デ
ータ点のデータに作用させる手段とを有するMRI装置
(1) A means for determining the distance between each data point and a predetermined reference point in the collected data matrix, and a storage means in which a filter value given as a function of the distance from the reference point is stored in advance for each distance. and means for reading a filter value corresponding to the distance from the storage means and applying it to data of each data point.
JP62005580A 1987-01-13 1987-01-13 Mri apparatus Pending JPS63174642A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP62005580A JPS63174642A (en) 1987-01-13 1987-01-13 Mri apparatus

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP62005580A JPS63174642A (en) 1987-01-13 1987-01-13 Mri apparatus

Publications (1)

Publication Number Publication Date
JPS63174642A true JPS63174642A (en) 1988-07-19

Family

ID=11615172

Family Applications (1)

Application Number Title Priority Date Filing Date
JP62005580A Pending JPS63174642A (en) 1987-01-13 1987-01-13 Mri apparatus

Country Status (1)

Country Link
JP (1) JPS63174642A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH08506755A (en) * 1993-10-06 1996-07-23 ピ,ダニエル Cartridge for applying drug to eyes

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH08506755A (en) * 1993-10-06 1996-07-23 ピ,ダニエル Cartridge for applying drug to eyes

Similar Documents

Publication Publication Date Title
EP0105700B1 (en) Nuclear magnetic resonance methods
WO1995024682A1 (en) Integrating texture memory and interpolation logic
CN112014782A (en) Magnetic resonance imaging method, magnetic resonance imaging apparatus, and computer storage medium
JPS63174642A (en) Mri apparatus
US4841247A (en) Method for interpolative extension of nuclear magnetic resonance image
JPH08251400A (en) Digital image interpolation circuit provided with plurality of interpolation kernels
EP0786738B1 (en) An image processing system and method
JPH0750124B2 (en) Method and apparatus for compressing digital time series data in a digital oscilloscope
JPS6012043A (en) Nmr imaging method
JPS63111848A (en) Magnetic resonance imaging apparatus
JPH0410095B2 (en)
JPH0475637A (en) Magnetic resonance video device
Feng et al. Fast direct fourier reconstruction of radial and PROPELLER MRI data using the chirp transform algorithm on graphics hardware
JPS6149246A (en) Image rotation processing system
JPS6055253A (en) Image formation method by nmr technique
JPS6126128A (en) Displacing method of two-dimensional data
JPH0748204B2 (en) Matrix calculation processing method
JP3850150B2 (en) Method and digital signal processor for performing digital signal processing operations
Scollar Interactive processing of geophysical data from archaeological sites
CN116047388A (en) Magnetic resonance imaging method, apparatus, computer device and storage medium
JPH02126364A (en) Processor for multidimensional data
JP2714926B2 (en) Image processing method and image processing apparatus
JPS62231648A (en) Nmr imaging method
JP2852050B2 (en) Image processing device
JPH05101153A (en) Three-dimensional data processing system