JPH07385A - Medical image processing system - Google Patents
Medical image processing systemInfo
- Publication number
- JPH07385A JPH07385A JP5165876A JP16587693A JPH07385A JP H07385 A JPH07385 A JP H07385A JP 5165876 A JP5165876 A JP 5165876A JP 16587693 A JP16587693 A JP 16587693A JP H07385 A JPH07385 A JP H07385A
- Authority
- JP
- Japan
- Prior art keywords
- image
- data
- filter
- post
- table memory
- 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.)
- Granted
Links
- 239000011159 matrix material Substances 0.000 claims abstract description 43
- 230000006870 function Effects 0.000 claims abstract description 27
- 230000001186 cumulative effect Effects 0.000 claims abstract description 20
- 238000001914 filtration Methods 0.000 claims description 12
- 238000005259 measurement Methods 0.000 claims description 12
- 238000007781 pre-processing Methods 0.000 claims description 6
- 230000005855 radiation Effects 0.000 claims 1
- 238000000034 method Methods 0.000 abstract description 14
- 238000010586 diagram Methods 0.000 description 4
- 230000000052 comparative effect Effects 0.000 description 2
- 238000004846 x-ray emission Methods 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
Abstract
Description
【０００１】[0001]
【産業上の利用分野】本発明は、例えばＸ線ＣＴ装置に
おいて被検体のある診断部位について収集した計測デー
タを取り込んで画像再構成すると共に画像表示する医用
画像処理装置に関し、特に画像再構成されたデータに対
して後補正処理としてマトリックスフィルタをかける際
にその画像に適するマトリックスフィルタを自動的に選
択して後補正処理ができる医用画像処理装置に関する。BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to a medical image processing apparatus that captures measurement data collected for a certain diagnostic region of a subject in an X-ray CT apparatus, reconstructs an image and displays the image, and particularly reconstructs the image. The present invention relates to a medical image processing apparatus capable of automatically selecting a matrix filter suitable for an image when applying a matrix filter as post-correction processing to the obtained data and performing the post-correction processing.
【０００２】[0002]
【従来の技術】従来のこの種の医用画像処理装置は、図
５に示すように、被検体の周りに互いに対向配置された
Ｘ線源及び多チャンネルの検出器を回転して検出すると
共にディジタル化された計測データを入力して前処理を
行う補正手段１と、この補正処理されたデータについて
画像のボケ補正のためのフィルタリング処理を行う手段
２と、このフィルタリング処理をされたデータについて
逆投影処理を行って２次元画像にする手段３と、この２
次元画像にされたデータについて後補正処理を行う手段
４と、この後補正をされた画像データを表示する手段５
とを有して成っていた。ここで、上記の後補正処理手段
４における補正処理としては、ＣＴ値補正などを行う
が、最近では画質の向上のために画像フィルタをかける
ことが行われている。この画像フィルタとしては、２次
元画像にかけるマトリックスフィルタが用いられてい
る。2. Description of the Related Art A conventional medical image processing apparatus of this type, as shown in FIG. 5, rotates and detects an X-ray source and a multi-channel detector, which are arranged to face each other around a subject, and digitally detects them. Correction means 1 for inputting the converted measurement data and performing preprocessing, means 2 for performing filtering processing for image blur correction on the corrected data, and backprojection for the filtered data A means 3 for processing to form a two-dimensional image;
Means 4 for performing post-correction processing on the data formed into the three-dimensional image, and means 5 for displaying the post-correction image data.
It was made up of and. Here, as the correction processing in the post-correction processing means 4, CT value correction or the like is performed, but recently image filtering is performed to improve image quality. As this image filter, a matrix filter applied to a two-dimensional image is used.
【０００３】上記のマトリックスフィルタの一例を図６
に示す。図６において、それぞれの四角形は画像の画素
を示しており、例えば９点あるとする（実際は25点ぐら
いあることが多い）。そして、注目する画素を真中のＰ
0とすると、この画素Ｐ0に対してその周囲の８点の画素
Ｐ1〜Ｐ8を用いてフィルタをかけることにより、新たな
画素Ｐ0′を求めるものである。これを数式を用いて示
すと、次の式（１）のようになる。 Ｐ0′＝Ｗ0Ｐ0＋Ｗ1Ｐ1＋Ｗ2Ｐ2＋…＋Ｗ8Ｐ8 …（１） ここで、Ｗ0，Ｗ1，Ｗ2，…，Ｗ8は、注目画素及びその
周囲の画素に対する重み付けの値であり、この値により
フィルタの特性が定まる。この重み付けの値（フィルタ
関数とも言う）は、予め数種類〜数十種類用意してお
き、色々な画質の画像に対処できるようになっている。An example of the above matrix filter is shown in FIG.
Shown in. In FIG. 6, each quadrangle represents a pixel of the image, and is assumed to have, for example, 9 points (actually, there are often 25 points). Then, the pixel of interest is set to P in the middle.
If it is set to 0 , a new pixel P 0 ′ is obtained by applying a filter to this pixel P 0 using the pixels P 1 to P 8 at eight points around it. If this is shown using a mathematical expression, the following Expression (1) is obtained. P 0 ′ = W 0 P 0 + W 1 P 1 + W 2 P 2 + ... + W 8 P 8 (1) Here, W 0 , W 1 , W 2 , ..., W 8 are the target pixel and its surroundings. It is a weighting value for a pixel, and this value determines the characteristics of the filter. Several to several tens of weighting values (also referred to as filter functions) are prepared in advance so that images of various image qualities can be dealt with.
【０００４】[0004]
【発明が解決しようとする課題】しかし、このような従
来の医用画像処理装置においては、図５に示す後補正処
理手段４で行う補正処理に用いる画像フィルタのフィル
タ関数は、操作者が人手によって画像に適すると思われ
るものを選択して適用していた。すなわち、まず、何も
画像フィルタをかけない状態で画像を再構成すると共に
表示し、この表示画像を観察して予め用意された複数種
類のフィルタ関数のうちいずれを使用するか操作者が判
断して選択し、それにより図６に示すようにマトリック
スフィルタをかけた画像を作り直して表示していた。次
に、この作り直した画像を観察してまだ十分でない場合
は、さらに別のフィルタ関数を選択して、再度図６に示
すようにマトリックスフィルタをかけた画像を作り直し
て表示していた。その後、さらにこの作り直した画像を
観察してまだ十分でない場合は、さらに別のフィルタ関
数を選択して画像を作り直して表示していた。このよう
に、従来装置においては、後補正処理に用いる画像フィ
ルタのフィルタ関数は、操作者が人手によって試行錯誤
的に選択して適用していたので、画像データに後補正処
理を施すのに多くの時間と手間がかかるものであった。
また、操作者の個人的な判断によるフィルタ関数の選択
であるので、後補正処理後の画像の画質として個人差が
出ることがあり、複数人で読影する場合に見づらいこと
があった。However, in such a conventional medical image processing apparatus, the filter function of the image filter used for the correction processing performed by the post-correction processing means 4 shown in FIG. I selected and applied what seemed to be suitable for the image. That is, first, the image is reconstructed and displayed without applying any image filter, and the displayed image is observed, and the operator judges which of a plurality of types of filter functions prepared in advance is to be used. Then, the image on which the matrix filter has been applied is recreated and displayed as shown in FIG. Next, if the regenerated image is observed and it is still insufficient, another filter function is selected and the image on which the matrix filter is applied as shown in FIG. 6 is regenerated and displayed. After that, if the reconstructed image is observed further and it is still not sufficient, another filter function is selected and the image is reconstructed and displayed. As described above, in the conventional apparatus, the filter function of the image filter used for the post-correction process is manually selected by the operator on a trial and error basis and applied, so that the post-correction process is often performed on the image data. It took a lot of time and effort.
Moreover, since the filter function is selected by the operator's personal judgment, the image quality of the image after the post-correction process may vary from person to person, and it may be difficult for a plurality of people to interpret the image.
【０００５】そこで、本発明は、このような問題点に対
処し、画像再構成されたデータに対して後補正処理とし
てマトリックスフィルタをかける際にその画像に適する
マトリックスフィルタを自動的に選択して後補正処理が
できる医用画像処理装置を提供することを目的とする。Therefore, the present invention addresses such a problem, and when a matrix filter is applied as post-correction processing to image reconstructed data, a matrix filter suitable for the image is automatically selected. An object of the present invention is to provide a medical image processing device that can perform post-correction processing.
【０００６】[0006]
【課題を解決するための手段】上記目的を達成するため
に、本発明による医用画像処理装置は、被検体の周りに
互いに対向配置されたＸ線源及び多チャンネルの検出器
を回転して検出すると共にディジタル化された計測デー
タを入力して前処理を行う補正手段と、この補正処理さ
れたデータについて画像のボケ補正のためのフィルタリ
ング処理を行う手段と、このフィルタリング処理をされ
たデータについて逆投影処理を行って２次元画像にする
手段と、この２次元画像にされたデータについて後補正
処理を行う手段と、この後補正をされた画像データを表
示する手段とを有する医用画像処理装置において、上記
フィルタリング処理手段におけるフーリエ変換後のデー
タを取り込んでそのデータを周波数成分ごとに区切って
上記の検出器によるビュー方向に加算する累積加算器
と、上記の周波数成分ごとにある域値幅を有する比較値
を並べた表を複数種類備えた比較表メモリと、上記累積
加算器からの加算結果を入力し比較表メモリからの表デ
ータと比較していずれの表と一致するかを調べる比較器
と、上記比較表メモリ内の複数種類の表にそれぞれ対応
する複数種類のマトリックスフィルタ関数を備え上記比
較器からの一致信号によりそれに対応する一つのマトリ
ックスフィルタ関数を選択するフィルタテーブルメモリ
とから成る後補正処理の自動制御系を設け、このフィル
タテーブルメモリで選択されたマトリックスフィルタ関
数を上記後補正処理手段に送出して画像データにマトリ
ックスフィルタをかけるようにしたものである。In order to achieve the above object, a medical image processing apparatus according to the present invention detects an object by rotating an X-ray source and a multi-channel detector which are arranged to face each other around an object. In addition, correction means for inputting digitized measurement data and performing preprocessing, means for performing filtering processing for image blur correction on the corrected data, and inverse processing for the filtered data In a medical image processing apparatus having means for performing projection processing to form a two-dimensional image, means for performing post-correction processing on the data converted into the two-dimensional image, and means for displaying the post-corrected image data , The data after the Fourier transform in the filtering processing means is taken in, and the data is sectioned for each frequency component, A cumulative adder for adding in the view direction, a comparative table memory having a plurality of types of tables in which comparative values having a certain threshold width for each frequency component are arranged, and a comparison table for inputting the addition results from the cumulative adder. Comparing the table data from the memory with a comparator that checks which table matches and a matrix filter function that corresponds to each of the tables in the comparison table memory. An automatic control system for post-correction processing, which comprises a filter table memory that selects one matrix filter function corresponding to the signal, is provided, and the matrix filter function selected in this filter table memory is sent to the post-correction processing means. The matrix filter is applied to the image data.
【０００７】[0007]
【作用】このように構成された医用画像処理装置は、累
積加算器と比較表メモリと比較器とフィルタテーブルメ
モリとから成る後補正処理の自動制御系を設けたことに
より、フィルタリング処理手段におけるフーリエ変換後
のデータを上記累積加算器に取り込み、上記フィルタテ
ーブルメモリで選択されたマトリックスフィルタ関数を
後補正処理手段に送出して、画像データにマトリックス
フィルタをかけるように動作する。これにより、画像再
構成されたデータに対して後補正処理としてマトリック
スフィルタをかける際にその画像に適するマトリックス
フィルタを自動的に選択して後補正処理ができる。The medical image processing apparatus thus constructed is provided with an automatic control system for post-correction processing which comprises a cumulative adder, a comparison table memory, a comparator and a filter table memory, so that the Fourier processing in the filtering processing means is performed. The converted data is taken into the cumulative adder, the matrix filter function selected in the filter table memory is sent to the post-correction processing means, and the image data is subjected to a matrix filter operation. As a result, when the matrix filter is applied to the image-reconstructed data as the post-correction process, the matrix filter suitable for the image is automatically selected and the post-correction process can be performed.
【０００８】[0008]
【実施例】以下、本発明の実施例を添付図面に基づいて
詳細に説明する。図１は本発明による医用画像処理装置
の実施例を示すブロック図である。この医用画像処理装
置は、例えばＸ線ＣＴ装置において被検体のある診断部
位について収集した計測データを取り込んで画像再構成
すると共に画像表示するもので、図１に示すように、補
正演算器１１と、フーリエ変換器１２と、フィルタ乗算
器１３と、逆フーリエ変換器１４と、逆投影演算器１５
と、画像フィルタ演算器１６と、表示装置１７とを有
し、さらに累積加算器１８と、比較表メモリ１９と、比
較器２０と、フィルタテーブルメモリ２１とを備えて成
る。Embodiments of the present invention will now be described in detail with reference to the accompanying drawings. FIG. 1 is a block diagram showing an embodiment of a medical image processing apparatus according to the present invention. This medical image processing apparatus, for example, captures measurement data collected for a certain diagnosis site of a subject in an X-ray CT apparatus, reconstructs an image, and displays an image. As shown in FIG. , Fourier transformer 12, filter multiplier 13, inverse Fourier transformer 14, and backprojection calculator 15
It has an image filter calculator 16 and a display device 17, and further comprises a cumulative adder 18, a comparison table memory 19, a comparator 20, and a filter table memory 21.
【０００９】上記補正演算器１１は、図示外の被検体を
中心として互いに対向配置されたＸ線源及び多チャンネ
ルの検出器を上記被検体の周りに回転して検出すると共
にディジタル化された計測データを入力して前処理を行
う補正手段となるもので、入力系の特性のばらつきなど
を補正するようになっている。The correction calculator 11 detects an X-ray source and a multi-channel detector, which are arranged opposite to each other around a subject (not shown) as a center, by rotating the detector around the subject and performing digitized measurement. It serves as a correction unit that inputs data and performs preprocessing, and corrects variations in the characteristics of the input system.
【００１０】フーリエ変換器１２は、上記補正演算器１
１で前処理された計測データを入力してフーリエ変換し
実空間から周波数空間へ変換するもので、図２（ａ）に
示すように縦軸に前記多チャンネルの検出器によるＸ線
放射に対するビュー方向をとり、横軸に上記検出器の各
チャンネルをとって表せる実空間の入力データを、同図
（ｂ）に示すように同じく縦軸に上記検出器によるＸ線
放射に対するビュー方向をとり、横軸には該検出器のチ
ャンネル方向を周波数として表せる周波数空間のデータ
に変換するようになっている。また、フィルタ乗算器１
３は、上記フーリエ変換器１２で周波数空間に変換され
た周波数データを入力してボケ補正のためのフィルタデ
ータを乗算するものである。さらに、逆フーリエ変換器
１４は、上記フィルタ乗算器１３でボケ補正がされた周
波数データを入力して逆フーリエ変換し周波数空間から
実空間へ変換するものである。そして、これらフーリエ
変換器１２とフィルタ乗算器１３と逆フーリエ変換器１
４とで、上記補正演算器１１で前処理された計測データ
について画像のボケ補正のためのフィルタリング処理を
行う手段２２を構成している。The Fourier transformer 12 is the correction calculator 1 described above.
The measurement data pre-processed in 1 is input and Fourier-transformed to convert from a real space to a frequency space. As shown in FIG. 2 (a), the vertical axis is a view for X-ray emission by the multi-channel detector. Direction, and the horizontal axis represents input data in real space that can be represented by taking each channel of the detector, and the vertical axis also represents the view direction for X-ray emission by the detector, as shown in FIG. On the horizontal axis, the channel direction of the detector is converted into frequency space data that can be expressed as frequency. In addition, the filter multiplier 1
3 is for inputting the frequency data converted into the frequency space by the Fourier transformer 12 and multiplying it by the filter data for blur correction. Further, the inverse Fourier transformer 14 inputs the frequency data subjected to the blur correction by the filter multiplier 13 and inverse Fourier transforms it to transform it from the frequency space to the real space. Then, these Fourier transformer 12, filter multiplier 13 and inverse Fourier transformer 1
4 constitutes a means 22 for performing filtering processing for image blur correction on the measurement data preprocessed by the correction calculator 11.
【００１１】逆投影演算器１５は、上記フィルタリング
処理手段２２でフィルタリング処理をされたデータを入
力して逆投影処理を行って２次元画像に再構成する手段
となるものである。また、画像フィルタ演算器１６は、
上記逆投影演算器１５で２次元画像に再構成された画像
データを入力して後補正処理を行う手段となるもので、
上記再構成された２次元画像に対して画質の向上のため
にマトリックスフィルタから成る画像フィルタをかける
ようになっている。さらに、表示装置１７は、上記画像
フィルタ演算器１６で後補正処理をされた画像データを
入力してアナログ信号に変換し画像として表示する手段
となるもので、例えばテレビモニタから成る。The back-projection calculator 15 serves as means for inputting the data filtered by the filtering processing means 22 and performing back-projection processing to reconstruct a two-dimensional image. In addition, the image filter calculator 16
It serves as a means for performing post-correction processing by inputting image data reconstructed into a two-dimensional image by the back projection calculator 15.
An image filter composed of a matrix filter is applied to the reconstructed two-dimensional image in order to improve the image quality. Further, the display device 17 is a means for inputting the image data which has been subjected to the post-correction processing by the image filter arithmetic unit 16, converting it into an analog signal and displaying it as an image, and is composed of, for example, a television monitor.
【００１２】ここで、本発明においては、上記フィルタ
リング処理手段２２のフーリエ変換器１２で変換された
後のデータを取り込み、所要の動作により選択されたマ
トリックスフィルタ関数を上記画像フィルタ演算器１６
へ送出する後補正処理の自動制御系２３が設けられてい
る。そして、この自動制御系２３は、累積加算器１８
と、比較表メモリ１９と、比較器２０と、フィルタテー
ブルメモリ２１とから成る。Here, in the present invention, the data after being transformed by the Fourier transformer 12 of the filtering processing means 22 is fetched, and the matrix filter function selected by the required operation is applied to the image filter calculator 16
An automatic control system 23 for post-correction processing is provided for sending to the. Then, this automatic control system 23 uses the cumulative adder 18
And a comparison table memory 19, a comparator 20, and a filter table memory 21.
【００１３】上記累積加算器１８は、前記フーリエ変換
器１２でフーリエ変換した後のデータを取り込んでその
データを周波数成分ごとに区切って前記多チャンネルの
検出器によるビュー方向に加算するもので、前述の図２
（ｂ）に示すフーリエ変換後のデータをビュー方向に加
算を行い、図３（ａ）に示すように縦軸に加算値をと
り、横軸に周波数をとって表される加算結果を、同図上
である周波数範囲ＡＢごとにまとめて加算（斜線部参
照）するようになっている。これを周波数の全範囲にわ
たってある周波数成分ごとに区切って加算した様子を示
したのが図３（ｂ）である。The accumulator 18 takes in the data after the Fourier transform by the Fourier transformer 12, divides the data into frequency components, and adds them in the view direction by the multi-channel detector. Figure 2
The Fourier-transformed data shown in (b) is added in the view direction, and as shown in FIG. 3 (a), the addition value is plotted on the vertical axis and the frequency is plotted on the horizontal axis. The frequency ranges AB in the figure are collectively added (see shaded portions). FIG. 3B shows a state in which this is divided into certain frequency components over the entire frequency range and added.
【００１４】比較表メモリ１９は、上記累積加算器１８
の図３（ｂ）に示す周波数成分ごとにある域値幅を有す
る比較値を並べた表（テーブル）を複数種類格納したも
ので、例えば図４（ａ）,（ｂ）に示すように、縦軸に
比較値をとり、横軸にある範囲の周波数成分ごとに区切
った周波数をとり、各周波数成分の幅に対して上下方向
の矢印で示すある域値幅を有する比較値を並べたもので
ある。そして、上記比較値の並び方のパターンは、周波
数に対して図４（ａ）,（ｂ）に示すようにそれぞれ異
なっており、色々な特性の表が用意されている。The comparison table memory 19 includes the cumulative adder 18 described above.
3 (b), a plurality of types of tables in which comparison values having a certain range of threshold values are arranged for each frequency component are stored. For example, as shown in FIGS. 4 (a) and 4 (b), The comparison value is plotted on the axis, the frequency is divided for each frequency component in the range on the horizontal axis, and the comparison values having a certain range width indicated by an up and down arrow with respect to the width of each frequency component are arranged. . The patterns of arrangement of the comparison values are different with respect to the frequencies as shown in FIGS. 4A and 4B, and various characteristic tables are prepared.
【００１５】また、比較器２０は、上記累積加算器１８
からの加算結果を入力し比較表メモリ１９からの表デー
タと比較していずれの表と一致するかを調べるもので、
累積加算器１８から入力する図３（ｂ）に示すようなパ
ターンの加算結果と、比較表メモリ１９から入力する図
４（ａ）,（ｂ）に示すようなパターンの色々な特性の
表とを比較し、図３（ｂ）に示すパターンの加算結果と
略一致するパターンの比較表を調べ、その一致信号Ｓを
後述のフィルタテーブルメモリ２１へ送出するようにな
っている。Further, the comparator 20 is the cumulative adder 18 described above.
The addition result from is input and compared with the table data from the comparison table memory 19 to check which table matches,
A result of addition of the pattern as shown in FIG. 3B input from the cumulative adder 18 and a table of various characteristics of the pattern as shown in FIGS. 4A and 4B input from the comparison table memory 19. Are compared with each other, a comparison table of patterns that substantially match the addition result of the patterns shown in FIG. 3B is checked, and the matching signal S is sent to the filter table memory 21 described later.
【００１６】さらに、フィルタテーブルメモリ２１は、
上記比較表メモリ１９内の複数種類の表にそれぞれ対応
する複数種類のマトリックスフィルタ関数を備え上記比
較器２０からの一致信号Ｓによりそれに対応する一つの
マトリックスフィルタ関数を選択するもので、前述の式
（１）に示されるような画素に対する重み付けの値
Ｗ0，Ｗ1，Ｗ2，…，Ｗ8から成るマトリックスフィルタ
関数が、図４（ａ）,（ｂ）に示す複数の表に対応させ
てそれぞれ格納されている。このマトリックスフィルタ
関数の対応関係は、予め多数の画像を調査して、どの比
較表にどのフィルタ関数を対応させるかを決めておけば
よい。そして、このフィルタテーブルメモリ２１で選択
されたマトリックスフィルタ関数は、前記画像フィルタ
演算器１６へ送出されるようになっている。Further, the filter table memory 21 is
A plurality of types of matrix filter functions respectively corresponding to a plurality of types of tables in the comparison table memory 19 are provided, and one matrix filter function corresponding thereto is selected by the coincidence signal S from the comparator 20. A matrix filter function consisting of weighting values W 0 , W 1 , W 2 , ..., W 8 for pixels as shown in (1) corresponds to a plurality of tables shown in FIGS. 4 (a) and 4 (b). Are stored separately. As for the correspondence relationship of this matrix filter function, a large number of images may be investigated in advance to determine which filter function is associated with which comparison table. Then, the matrix filter function selected in the filter table memory 21 is sent to the image filter calculator 16.
【００１７】次に、このように構成された医用画像処理
装置の動作について説明する。まず、図１において、デ
ィジタル化されて入力した計測データは、補正演算器１
１で前処理としての補正を行う。この前処理後の計測デ
ータは、次のフーリエ変換器１２へ入力して、図２
（ａ）に示すチャンネル−ビュー空間（実空間）から周
波数−ビュー空間（周波数空間）に変換される。その
後、フィルタ乗算器１３でボケ補正が行われ、逆フーリ
エ変換器１４で周波数空間から実空間へ変換される。そ
して、逆投影演算器１５で２次元画像に再構成されると
共に、次の画像フィルタ演算器１６で後補正処理として
画質向上のための画像フィルタがかけられる。その後、
この後補正処理された画像データが表示装置１７に画像
として表示される。Next, the operation of the medical image processing apparatus configured as described above will be described. First, in FIG. 1, the digitized measurement data is input to the correction calculator 1
At 1, correction is performed as preprocessing. The measurement data after this pre-processing is input to the next Fourier transformer 12,
The channel-view space (real space) shown in (a) is converted to the frequency-view space (frequency space). Thereafter, the filter multiplier 13 performs blur correction, and the inverse Fourier transformer 14 transforms the frequency space into a real space. Then, the back-projection computing unit 15 reconstructs a two-dimensional image, and the next image-filter computing unit 16 applies an image filter for improving the image quality as a post-correction process. afterwards,
After this, the corrected image data is displayed as an image on the display device 17.
【００１８】このとき、上記後補正処理のために、フー
リエ変換器１２でフーリエ変換されたデータは累積加算
器１８にも取り込まれる。この累積加算器１８では、上
記フーリエ変換後のデータを図３（ａ）に示すようにビ
ュー方向に加算を行い、さらに同図（ｂ）に示すように
周波数成分ごとに区切ってビュー方向に累積加算を行
う。そして、この加算結果を次の比較器２０へ送る。比
較器２０では、上記の加算結果と、比較表メモリ１９か
ら読み出した図４（ａ）,（ｂ）に示すような色々な特
性の表とを比較し、例えば図３（ｂ）に示すパターンの
加算結果と略一致するパターンの比較表を調べる。この
とき、周波数を基準にして総ての周波数成分において、
図３（ｂ）に示すまとめ加算値が図４（ａ）,（ｂ）に
示す比較値の域値幅内にあるかどうかを比較し、当ては
まるものを選択する。ここでは、図４（ｂ）の比較表が
略一致するのでこれを選択し、この一致信号Ｓを出力す
る。At this time, the data Fourier-transformed by the Fourier transformer 12 is also fetched by the cumulative adder 18 for the post-correction processing. In the cumulative adder 18, the data after the Fourier transform is added in the view direction as shown in FIG. 3A, and further divided into frequency components as shown in FIG. 3B and accumulated in the view direction. Add. Then, this addition result is sent to the next comparator 20. The comparator 20 compares the above addition result with the tables of various characteristics read from the comparison table memory 19 as shown in FIGS. 4 (a) and 4 (b), for example, the pattern shown in FIG. 3 (b). Examine the comparison table for patterns that approximately match the result of addition. At this time, in all frequency components based on the frequency,
It is compared whether the summary addition value shown in FIG. 3B is within the threshold range of the comparison values shown in FIGS. 4A and 4B, and the one that applies is selected. Here, since the comparison table of FIG. 4 (b) substantially matches, it is selected and the matching signal S is output.
【００１９】すると、この一致信号Ｓは、次のフィルタ
テーブルメモリ２１へ入力する。このフィルタテーブル
メモリ２１は、上記一致信号Ｓの入力によりこの一致信
号Ｓが表す比較表に対応する一つのマトリックスフィル
タ関数を選択し、この選択されたマトリックスフィルタ
関数を前記画像フィルタ演算器１６に送出する。これに
より、上記画像フィルタ演算器１６は、逆投影演算器１
５で２次元画像に再構成された画像データに対して上記
選択して入力されたマトリックスフィルタ関数を適用
し、前述の式（１）に示される演算を実行してマトリッ
クスフィルタをかける。この結果、当該画像に適したマ
トリックスフィルタを自動的に選択して後補正処理が実
行される。Then, the coincidence signal S is input to the next filter table memory 21. The filter table memory 21 selects one matrix filter function corresponding to the comparison table represented by the match signal S by the input of the match signal S, and sends the selected matrix filter function to the image filter calculator 16. To do. As a result, the image filter arithmetic unit 16 operates in the back projection arithmetic unit 1.
The selected and input matrix filter function is applied to the image data reconstructed into the two-dimensional image in 5, and the matrix filter is applied by executing the operation represented by the above-mentioned equation (1). As a result, the post-correction process is executed by automatically selecting the matrix filter suitable for the image.
【００２０】[0020]
【発明の効果】本発明は以上のように構成されたので、
累積加算器と比較表メモリと比較器とフィルタテーブル
メモリとから成る後補正処理の自動制御系を設けたこと
により、フィルタリング処理手段におけるフーリエ変換
後のデータを上記累積加算器に取り込み、上記フィルタ
テーブルメモリで選択されたマトリックスフィルタ関数
を後補正処理手段に送出して、画像データにマトリック
スフィルタをかけることができる。これにより、画像再
構成されたデータに対して後補正処理としてマトリック
スフィルタをかける際にその画像に適するマトリックス
フィルタを自動的に選択して後補正処理ができる。従っ
て、従来のように後補正処理に用いる画像フィルタのフ
ィルタ関数を操作者が人手により試行錯誤的に選択する
のをなくし、画像データの後補正処理を短時間で容易に
実行することができる。また、後補正処理後の画質を均
一化することができ、個人差をなくして見易い画像とす
ることができる。Since the present invention is constructed as described above,
By providing an automatic control system for post-correction processing consisting of a cumulative adder, a comparison table memory, a comparator and a filter table memory, the data after the Fourier transform in the filtering processing means is taken into the cumulative adder and the filter table is used. The matrix filter function selected in the memory can be sent to the post-correction processing means to apply the matrix filter to the image data. As a result, when the matrix filter is applied to the image-reconstructed data as the post-correction process, the matrix filter suitable for the image is automatically selected and the post-correction process can be performed. Therefore, it is possible to prevent the operator from manually selecting the filter function of the image filter used for the post-correction processing by trial and error as in the related art, and to easily execute the post-correction processing of the image data in a short time. In addition, the image quality after the post-correction processing can be made uniform, and an image that is easy to see can be obtained by eliminating individual differences.
【図１】本発明による医用画像処理装置の実施例を示す
ブロック図である。FIG. 1 is a block diagram showing an embodiment of a medical image processing apparatus according to the present invention.
【図２】フーリエ変換器による計測データの周波数空間
への変換を説明するためのグラフである。FIG. 2 is a graph for explaining conversion of measurement data into a frequency space by a Fourier transformer.
【図３】累積加算器による計測データの周波数成分ごと
に区切ってビュー方向に累積加算する状態を説明するた
めのグラフである。FIG. 3 is a graph for explaining a state in which a cumulative addition is performed for each frequency component of measurement data and cumulative addition is performed in a view direction.
【図４】比較表メモリに格納された色々な特性を有する
比較表の例を示すグラフである。FIG. 4 is a graph showing an example of a comparison table having various characteristics stored in a comparison table memory.
【図５】従来の医用画像処理装置を示すブロック図であ
る。FIG. 5 is a block diagram showing a conventional medical image processing apparatus.
【図６】マトリックスフィルタの一例を示す説明図であ
る。FIG. 6 is an explanatory diagram showing an example of a matrix filter.
１１…補正演算器 １２…フーリエ変換器 １３…フィルタ乗算器 １４…逆フーリエ変換器 １５…逆投影演算器 １６…画像フィルタ演算器 １７…表示装置 １８…累積加算器 １９…比較表メモリ ２０…比較器 ２１…フィルタテーブルメモリ ２２…フィルタリング処理手段 ２３…後補正処理の自動制御系 Ｓ…一致信号 11 ... Correction calculator 12 ... Fourier transformer 13 ... Filter multiplier 14 ... Inverse Fourier transformer 15 ... Back projection calculator 16 ... Image filter calculator 17 ... Display device 18 ... Cumulative adder 19 ... Comparison table memory 20 ... Comparison Device 21 ... Filter table memory 22 ... Filtering processing means 23 ... Automatic control system of post-correction processing S ... Matching signal
Claims (1)
線源及び多チャンネルの検出器を回転して検出すると共
にディジタル化された計測データを入力して前処理を行
う補正手段と、この補正処理されたデータについて画像
のボケ補正のためのフィルタリング処理を行う手段と、
このフィルタリング処理をされたデータについて逆投影
処理を行って２次元画像にする手段と、この２次元画像
にされたデータについて後補正処理を行う手段と、この
後補正をされた画像データを表示する手段とを有する医
用画像処理装置において、上記フィルタリング処理手段
におけるフーリエ変換後のデータを取り込んでそのデー
タを周波数成分ごとに区切って上記の検出器によるビュ
ー方向に加算する累積加算器と、上記の周波数成分ごと
にある域値幅を有する比較値を並べた表を複数種類備え
た比較表メモリと、上記累積加算器からの加算結果を入
力し比較表メモリからの表データと比較していずれの表
と一致するかを調べる比較器と、上記比較表メモリ内の
複数種類の表にそれぞれ対応する複数種類のマトリック
スフィルタ関数を備え上記比較器からの一致信号により
それに対応する一つのマトリックスフィルタ関数を選択
するフィルタテーブルメモリとから成る後補正処理の自
動制御系を設け、このフィルタテーブルメモリで選択さ
れたマトリックスフィルタ関数を上記後補正処理手段に
送出して画像データにマトリックスフィルタをかけるよ
うにしたことを特徴とする医用画像処理装置。1. Xs arranged to face each other around a subject.
Correcting means for performing preprocessing by rotating and detecting the radiation source and the multi-channel detector and inputting digitized measurement data, and filtering processing for image blur correction of the corrected data. Means to do
Back-projection processing is performed on the filtered data to form a two-dimensional image, post-correction processing is performed on the two-dimensional image data, and post-corrected image data is displayed. In the medical image processing device having means, a cumulative adder for taking in the data after the Fourier transform in the filtering processing means, dividing the data for each frequency component and adding in the view direction by the detector, and the frequency A comparison table memory having a plurality of types of tables in which comparison values having a certain threshold value width are arranged for each component, and the addition result from the cumulative adder is input and compared with the table data from the comparison table memory to determine which table Comparator to check if they match, and multiple types of matrix filter functions corresponding to the multiple types of tables in the comparison table memory An automatic control system for post-correction processing, which comprises a filter table memory for selecting one matrix filter function corresponding to the match signal from the comparator, is provided, and the matrix filter function selected in this filter table memory is A medical image processing apparatus characterized in that a matrix filter is applied to image data sent to a correction processing means.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP16587693A JP3383010B2 (en) | 1993-06-14 | 1993-06-14 | Medical image processing equipment |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP16587693A JP3383010B2 (en) | 1993-06-14 | 1993-06-14 | Medical image processing equipment |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH07385A true JPH07385A (en) | 1995-01-06 |
| JP3383010B2 JP3383010B2 (en) | 2003-03-04 |
Family
ID=15820662
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP16587693A Expired - Fee Related JP3383010B2 (en) | 1993-06-14 | 1993-06-14 | Medical image processing equipment |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3383010B2 (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7460700B2 (en) | 2003-07-24 | 2008-12-02 | Nihon University | Image processing method and computer-readable recording medium containing image processing program |
| DE102020212382B3 (en) | 2020-09-30 | 2022-01-20 | Siemens Healthcare Gmbh | Reconstruction of images with filter matrices |
-
1993
- 1993-06-14 JP JP16587693A patent/JP3383010B2/en not_active Expired - Fee Related
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7460700B2 (en) | 2003-07-24 | 2008-12-02 | Nihon University | Image processing method and computer-readable recording medium containing image processing program |
| DE102020212382B3 (en) | 2020-09-30 | 2022-01-20 | Siemens Healthcare Gmbh | Reconstruction of images with filter matrices |
Also Published As
| Publication number | Publication date |
|---|---|
| JP3383010B2 (en) | 2003-03-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US4189775A (en) | Method and apparatus for subsection analysis using computed tomography | |
| JP6214128B2 (en) | Image processing apparatus, image processing method, and storage medium | |
| US6704440B1 (en) | Method and apparatus for processing a medical image containing clinical and non-clinical regions | |
| US5606585A (en) | Methods and apparatus for multislice helical image reconstruction in a computer tomography system | |
| DE3918354A1 (en) | METHOD FOR SHORTERING THE IMAGE RECONSTRUCTION TIME IN COMPUTER TOMOGRAPHY WITH A LIMITED ANGLE | |
| JP5828649B2 (en) | Image processing apparatus, image processing method, and computer program | |
| DE102007028828A1 (en) | Image data processing method for e.g. computer tomography system, involves providing image with subsets of pixels, and modifying subsets between iterations such that pixel composition and size of subsets differs with each iteration | |
| DE102010029281A1 (en) | Method and image reconstruction device for the reconstruction of image data | |
| JP6987352B2 (en) | Medical image processing equipment and medical image processing method | |
| US11291422B2 (en) | Reconstructing cardiac frequency phenomena in angiographic data | |
| CN105916445B (en) | X ray CT device and X ray CT device image operation device | |
| DE69933338T2 (en) | Computer-controlled multi-frame tomographic frame reconstruction method and apparatus for spiral scanning | |
| KR20050042727A (en) | Ct image producing method and x-ray ct apparatus | |
| JPH07385A (en) | Medical image processing system | |
| US6647084B1 (en) | Method and apparatus for filtering projection data of a helical scan | |
| US5761267A (en) | Methods and apparatus for simplified filtering of scan data in an imaging system | |
| CN114943780A (en) | Image reconstruction processing method and system and medical imaging equipment | |
| JP3018194B2 (en) | X-ray CT scanner | |
| DE102016206297A1 (en) | Reconstruction of image data in sequence imaging | |
| AU2008313191A1 (en) | Image processing program, image processing method and image processor | |
| JPS63183045A (en) | Method and apparatus for re-constituting tomographic image | |
| JP3754807B2 (en) | Image processing apparatus, image processing method, and computer-readable storage medium | |
| JPH08272945A (en) | Image processor | |
| US20190383755A1 (en) | X-ray ct device and image display method | |
| JP3733263B2 (en) | Spatial filter creation method, radiation image processing apparatus, and recording medium |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20081220 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20091220 Year of fee payment: 7 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20101220 Year of fee payment: 8 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20101220 Year of fee payment: 8 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20111220 Year of fee payment: 9 |
|
| LAPS | Cancellation because of no payment of annual fees |