JP2001212119A - Image gaining apparatus and method - Google Patents
Image gaining apparatus and methodInfo
- Publication number
- JP2001212119A JP2001212119A JP2000028207A JP2000028207A JP2001212119A JP 2001212119 A JP2001212119 A JP 2001212119A JP 2000028207 A JP2000028207 A JP 2000028207A JP 2000028207 A JP2000028207 A JP 2000028207A JP 2001212119 A JP2001212119 A JP 2001212119A
- Authority
- JP
- Japan
- Prior art keywords
- image
- grid
- image acquisition
- frequency
- stripe
- 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
- 238000000034 method Methods 0.000 title claims description 18
- 238000005070 sampling Methods 0.000 claims abstract description 53
- 238000012545 processing Methods 0.000 claims abstract description 23
- 230000005855 radiation Effects 0.000 claims description 31
- 238000001914 filtration Methods 0.000 claims description 23
- 238000003384 imaging method Methods 0.000 claims description 9
- 230000005540 biological transmission Effects 0.000 claims description 6
- 230000001678 irradiating effect Effects 0.000 claims description 4
- 230000008569 process Effects 0.000 claims description 3
- 210000001015 abdomen Anatomy 0.000 claims description 2
- 210000000988 bone and bone Anatomy 0.000 claims description 2
- 210000004197 pelvis Anatomy 0.000 claims description 2
- 230000009467 reduction Effects 0.000 claims description 2
- 238000001228 spectrum Methods 0.000 description 19
- 238000010586 diagram Methods 0.000 description 14
- 238000006243 chemical reaction Methods 0.000 description 6
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 description 4
- 230000014509 gene expression Effects 0.000 description 4
- 230000007246 mechanism Effects 0.000 description 4
- 238000003745 diagnosis Methods 0.000 description 3
- 230000000903 blocking effect Effects 0.000 description 2
- 238000012952 Resampling Methods 0.000 description 1
- 230000035559 beat frequency Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 230000035807 sensation Effects 0.000 description 1
- 229910052709 silver Inorganic materials 0.000 description 1
- 239000004332 silver Substances 0.000 description 1
- -1 silver halide Chemical class 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5258—Devices using data or image processing specially adapted for radiation diagnosis involving detection or reduction of artifacts or noise
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Pathology (AREA)
- Heart & Thoracic Surgery (AREA)
- High Energy & Nuclear Physics (AREA)
- Physics & Mathematics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
Description
【０００１】[0001]
【発明の属する技術分野】本発明は、画像取得技術に関
し、特に被写体のＸ線透過分布に応じた画像を取得する
画像取得装置に関する。[0001] 1. Field of the Invention [0002] The present invention relates to an image acquiring technique, and more particularly to an image acquiring apparatus for acquiring an image corresponding to an X-ray transmission distribution of a subject.
【０００２】[0002]
【従来の技術】近年、医療用Ｘ線画像はデジタル化が進
み、Ｘ線強度の空間分布をデジタル画像として取得可能
になっている。たとえば、輝尽性蛍光体にＸ線エネルギ
ーにより潜像を作り、レーザによる励起光分布により画
像を取得する方式や、Ｘ線強度分布を光強度分布（蛍
光）に変換し、直接複数の画素を持つ面センサで電気信
号に変換した後デジタル画像に変換する方式や、Ｘ線強
度分布を直接電荷の分布に変換する方式などがある。2. Description of the Related Art In recent years, digitization of medical X-ray images has progressed, and the spatial distribution of X-ray intensity can be acquired as a digital image. For example, a method in which a latent image is formed on a stimulable phosphor by X-ray energy and an image is acquired by an excitation light distribution by a laser, or an X-ray intensity distribution is converted into a light intensity distribution (fluorescence), and a plurality of pixels are directly formed There are a method of converting the X-ray intensity distribution into a digital image after converting the electric signal into an electric signal by a surface sensor, and a method of directly converting an X-ray intensity distribution into a distribution of electric charges.
【０００３】Ｘ線画像をデジタル化することの利点は、 ・保存、転送の効率化。 ・デジタル画像処理により最適な画像が簡単に作り出せ
る（撮影時失敗を回復できる）。 ・診断の高効率化。 ・診断の低コスト化。 など、数々あげられる。[0003] The advantages of digitizing X-ray images are: • Efficient storage and transfer.・ Optimum images can be easily created by digital image processing (failure during shooting can be recovered).・ High efficiency of diagnosis.・ Low cost of diagnosis. And many more.
【０００４】しかし、Ｘ線が被写体を通過する時に発生
する散乱Ｘ線の問題はいまだに充分解決されておらず、
散乱Ｘ線の影響をなくし、コントラストの高い画像を得
るためには従来銀塩フィルムでおこなわれて来たと同様
の多数の鉛板を等方向にならべた散乱線除去グリッドを
用いるのが最良の手段となっている。However, the problem of scattered X-rays generated when X-rays pass through a subject has not yet been sufficiently solved.
In order to eliminate the influence of scattered X-rays and obtain a high-contrast image, it is best to use a scattered radiation removal grid in which a number of lead plates similar to those conventionally used for a silver halide film are arranged in the same direction. It has become.
【０００５】図１０はグリッドを使用する場合の断面図
を模式的に示したものであり、８１がＸ線管球のＸ線が
発生する点（Ｘ線焦点）であり、８２が被写体、８３が
グリッドであり、８４がＸ線強度分布を光強度もしくは
電荷量に変換するエネルギー変換部分を示し、８５がそ
の分布を空間的にサンプリングするセンサ部分である。
８３にはＸ線管球から直接到来するＸ線の他に被写体か
ら発せられる散乱Ｘ線も到達するが、８１のＸ線焦点に
向った鉛によるグリッド８３のためにその大半はカット
される。FIG. 10 schematically shows a cross-sectional view when a grid is used. Reference numeral 81 denotes a point (X-ray focal point) where an X-ray of an X-ray tube is generated; Denotes a grid, 84 denotes an energy conversion portion for converting the X-ray intensity distribution into light intensity or electric charge, and 85 denotes a sensor portion for spatially sampling the distribution.
The scattered X-rays emitted from the subject as well as the X-rays directly arriving from the X-ray tube reach the 83, but most of them are cut off by the grid 83 made of lead toward the X-ray focal point of the 81.
【０００６】グリッド８３の欠点は、散乱Ｘ線をカット
する一方、直接Ｘ線の一部も遮断する。その遮断の形態
はグリッドの鉛の配置に従い、通常縞状の損傷を画像に
残すことになる。Ｘ線画像は、（１）フィルムスクリー
ン系（アナログ画像）→（２）輝尽性蛍光体による潜像
をレーザスキャンで読み取りデジタル化→（３）Ｘ線量
の２次元空間分布を２次元空間で直接サンプリング（フ
ラットパネルセンサ）、と発展し、グリッド８３に使用
する鉛に起因する画像上の縞模様（グリッド像）に対す
る対策も変って来ている。A drawback of the grid 83 is that while cutting scattered X-rays, it also blocks some of the X-rays directly. The form of the block depends on the placement of the lead in the grid, and will usually leave banded damage in the image. The X-ray image is (1) film screen system (analog image) → (2) digitizing by reading the latent image by the stimulable phosphor by laser scanning → (3) two-dimensional spatial distribution of X-ray dose in two-dimensional space Direct sampling (flat panel sensor) has been developed, and measures against stripes (grid images) on an image caused by lead used for the grid 83 have been changed.
【０００７】上述の（１）のフィルムスクリーン系では
グリッド像を除去もしくは観察の邪魔にならないように
する方法に、以下の２つがある。 （ａ）Ｘ線曝射中にグリッドそのものを移動させ、散乱
線を除去しながらグリッド像を形成させないようにす
る。 （ｂ）グリッド縞の空間周波数を高め、グリッド像が画
像上に形成されたとしても、人眼には感知困難な状態に
するか、画像情報の周波数成分と重ならないようにす
る。In the film screen system of the above (1), there are the following two methods for removing the grid image or preventing it from hindering the observation. (A) The grid itself is moved during X-ray irradiation so that a grid image is not formed while removing scattered radiation. (B) The spatial frequency of grid stripes is increased so that even if a grid image is formed on an image, it is difficult for the human eye to perceive or a frequency component of image information does not overlap.
【０００８】Ｘ線画像取得のあらゆる場合に、上述の
（ａ）のグリッドそのものを移動させる手段は有効であ
る。しかし、移動にかかる駆動系などのコスト上昇・装
置の大型化、駆動タイミングとＸ線曝射タイミングの関
係、駆動速度の関係などの調整制御の手間などから、な
かなか採用できない。In all cases of acquiring an X-ray image, the means for moving the grid itself in (a) is effective. However, it cannot be easily adopted due to the cost increase of the drive system and the like, the size of the apparatus involved in the movement, the trouble of adjusting and controlling the relationship between the drive timing and the X-ray irradiation timing, the relationship between the drive speed, and the like.
【０００９】また、（ｂ）のグリッド縞の空間周波数を
高める手段にも限界がある。その理由はグリッド像を形
成されないような高空間周波数にグリッド縞の周波数を
設定すると、散乱線を遮断するための鉛板厚さはほぼ固
定されているため、直接線の通過する領域が狭まり、Ｘ
線量の利用効率が極端に低くなり、良好に撮影が行えな
いことにある。Also, there is a limit to the means for increasing the spatial frequency of the grid stripe in (b). The reason is that if the frequency of grid stripes is set to a high spatial frequency that does not form a grid image, the thickness of the lead plate for blocking scattered radiation is almost fixed, so the area where direct rays pass is narrowed, X
The use efficiency of the dose becomes extremely low, and good imaging cannot be performed.
【００１０】上述の（２）の輝尽性蛍光体による潜像を
レーザスキャンで読み取りデジタル化する時代になると
グリッド像を除去する方法にサンプリング前のアンチエ
リアシングフィルタと言う考え方を取り入れるようにな
る。輝尽性蛍光体による潜像をレーザスキャンで読み取
りデジタル化する場合は、レーザで１次元方向にスキャ
ンし、一旦ビデオ信号のような形状にして、時間軸上で
サンプリングを行う。グリッド縞の周波数をある程度高
くし、レーザスキャンをグリッド縞に直交する方向に行
い、グリッド縞をビデオ信号上での周期信号にする。こ
のビデオ信号の状態の時のアナログ信号の状態で低域通
過フィルタリングを行った後、時間軸でサンプリングを
行う通常のアンチエリアシングフィルタの考え方で、グ
リッドを除去できる。これに類するものとして特許第２
５０７６５９号では、グリッド像の存在および周波数を
予備的にサンプリングした画像をフーリエ変換すること
でもとめ、その結果に応じた低域通過フィルタを選択す
ることで、グリッド像を除去するものが開示されてい
る。In the age of digitizing the latent image formed by the stimulable phosphor by laser scanning in the above (2), the concept of an anti-aliasing filter before sampling is adopted in a method of removing a grid image. . When a latent image formed by a stimulable phosphor is read by a laser scan and digitized, the image is scanned in a one-dimensional direction with a laser, temporarily shaped like a video signal, and sampled on a time axis. The frequency of the grid stripes is increased to some extent, and laser scanning is performed in a direction orthogonal to the grid stripes, so that the grid stripes become periodic signals on the video signal. After performing low-pass filtering in the state of the analog signal in the state of the video signal, the grid can be removed by the concept of a normal anti-aliasing filter that performs sampling on the time axis. Patent No. 2 as similar
Japanese Patent No. 507659 discloses a technique in which an image obtained by preliminary sampling the existence and frequency of a grid image is determined by Fourier transform, and a low-pass filter according to the result is selected to remove the grid image. I have.
【００１１】さらに、アナログによる低域通過フィルタ
ではなく、時間軸のサンプリングを所望の間隔より短く
サンプリングし、グリッド縞情報のエリアシングを無く
し画像情報と分離した後に、デジタル的な低域通過フィ
ルタリングを行い、その後、デジタル的に間引いて（再
サンプリング）所望のサンプリング間隔による画像を得
るものも考えられた。これに類するものとして、特許第
２７５４０６８号、特開平８−０８８７６５が公開され
ている。Furthermore, instead of analog low-pass filtering, sampling on the time axis is performed at a shorter interval than a desired interval to eliminate aliasing of grid fringe information and separate it from image information. After that, it has been conceived to obtain an image at a desired sampling interval by digitally thinning out (resampling). Japanese Patent No. 2754068 and Japanese Patent Application Laid-Open No. 08-088765 have been disclosed as similar ones.
【００１２】さらに時代が進み、上述の（３）のＸ線量
の２次元空間分布を２次元空間で直接サンプリング（フ
ラットパネルセンサ）でデジタルＸ線画像を得るように
なると、上述のアンチエリアシングフィルタリングが使
えなくなって来る。すなわち、フラットパネルセンサの
２次元空間サンプリングピッチは、半導体による複数の
画素で構成され、そのピッチは技術的にもコスト的にも
必要以上に細かくできない。したがって、上述のアンチ
エリアシングフィルタリングの考え方が適応できない。
この２次元空間で直接サンプリングでデジタルＸ線画像
を得るものに対して、グリッド縞をとることを目的とす
るものに、特開平９−７５３３２があり、グリッド鉛の
間隔とサンプリングピッチを完全に一致させ、グリッド
縞による直接Ｘ線を遮断する領域と画素の隙間を一致さ
せグリッド縞を画像に出さないようにするものが開示さ
れている。As the era further advances and digital X-ray images are obtained by directly sampling (flat panel sensor) the two-dimensional spatial distribution of the X-ray amount of the above (3) in the two-dimensional space, the above-mentioned anti-aliasing filtering is performed. Comes to be unusable. That is, the two-dimensional spatial sampling pitch of the flat panel sensor is composed of a plurality of pixels made of semiconductor, and the pitch cannot be made finer than necessary in terms of technical and cost. Therefore, the above-described concept of anti-aliasing filtering cannot be applied.
Japanese Unexamined Patent Application Publication No. 9-75332 discloses a technique for obtaining a digital X-ray image by sampling directly in a two-dimensional space, which is intended to obtain grid stripes. There is disclosed an apparatus in which a gap between pixels and a region where direct X-rays are directly blocked by grid stripes is made coincident with each other to prevent grid stripes from appearing in an image.
【００１３】また、特開平９−７８９７０およびＵ．
Ｓ．Ｐａｔｅｎｔ５８０１３８５においては、グリッド
鉛間隔をサンプリングピッチより小さくし、１画素の持
つ受光部の開口の幅と同じくするか近づけグリッド縞の
コントラストを低減する方法が開示されている。また、
Ｕ．Ｓ．Ｐａｔｅｎｔ５，０５０，１９８においては、
複数の条件でグリッド像を撮影・記憶しておき、実際に
グリッドを用いて撮影した場合に、記憶してあるグリッ
ド像中でその条件にあったグリッド像にて除算を行うこ
とでグリッド像を除去する方法が開示されている。Further, Japanese Patent Application Laid-Open No. 9-78970 and U.S. Pat.
S. Patent 5801385 discloses a method in which the grid lead interval is made smaller than the sampling pitch, and the width of the aperture of the light receiving portion of one pixel is made equal to or close to the width, and the contrast of grid stripes is reduced. Also,
U. S. In Patent 5,050,198,
A grid image is photographed and stored under a plurality of conditions, and when a photograph is actually taken using a grid, the grid image is divided by a grid image that satisfies the conditions in the stored grid image to obtain a grid image. A method of removing is disclosed.
【００１４】[0014]
【発明が解決しようとする課題】上述の２次元空間でフ
ラットパネルセンサにて直接２次元空間におけるサンプ
リングでデジタルＸ線画像を得るものに対応するものと
して、特開平９−７５３３２のグリッド鉛間隔とサンプ
リングピッチを完全に一致させることは非常に困難であ
るという問題がある。また、特開平９−７８９７０およ
びＵ．Ｓ．Ｐａｔｅｎｔ５，８０１，３８５におけるグ
リッド鉛の間隔をサンプリングピッチより細かくし、１
画素の持つ受光部の開口の幅と同じくするか近づけるの
は有効であるが、フラットパネルセンサが高精細化し、
サンプリングピッチが０．１ｍｍ以下になると、グリッ
ド鉛の間隔も１ｍｍあたり１０本以上という非常に細か
なものが要求されるようになる。このような細かなもの
になると、散乱線を遮断するための鉛板厚さはほぼ固定
されているため、直接線の通過する領域が狭まり、Ｘ線
量の利用効率が極端に低くなり、良好に撮影が行えなく
なるという欠点がある。As a countermeasure for obtaining a digital X-ray image by sampling directly in a two-dimensional space with a flat panel sensor in a two-dimensional space as described above, Japanese Patent Application Laid-Open No. 9-75332 discloses a grid lead spacing. There is a problem that it is very difficult to completely match the sampling pitch. Also, JP-A-9-78970 and U.S. Pat. S. In the case of Patent 5,801,385, the interval between grid leads is made finer than the sampling pitch, and 1
It is effective to make the width equal to or close to the width of the aperture of the light receiving section of the pixel.
When the sampling pitch is 0.1 mm or less, very fine grid lead intervals of 10 or more per mm are required. With such a fine structure, since the thickness of the lead plate for blocking scattered radiation is almost fixed, the area through which direct rays pass is narrowed, and the use efficiency of X-ray dose becomes extremely low, and the There is a drawback that shooting cannot be performed.
【００１５】本発明の目的は、散乱線除去グリッドを使
用してＸ線を曝射した際に散乱線除去グリッドに起因す
る画像中の縞を必要に応じて除去することができる画像
取得装置及び画像取得方法を提供することである。SUMMARY OF THE INVENTION An object of the present invention is to provide an image acquisition apparatus capable of removing stripes in an image caused by a scattered radiation removal grid as required when X-rays are irradiated using the scattered radiation removal grid. An image acquisition method is provided.
【００１６】[0016]
【課題を解決するための手段】本発明の一観点によれ
ば、散乱線除去グリッドを使用してＸ線を曝射し、被写
体のＸ線透過分布を２次元的に所望の間隔でサンプリン
グして画像を取得する画像取得手段と、前記取得した画
像中から前記散乱線除去グリッドに起因する縞を画像処
理により除去する画像処理手段と、前記画像処理手段に
より前記散乱線除去グリッドに起因する縞を除去するか
否かの選択を可能にする選択手段とを有することを特徴
とする画像取得装置が提供される。According to one aspect of the present invention, X-rays are emitted using a scattered radiation removal grid, and an X-ray transmission distribution of a subject is sampled two-dimensionally at a desired interval. Image acquisition means for acquiring an image by performing image processing, image processing means for removing stripes caused by the scattered radiation removal grid from the acquired image by image processing, and stripes caused by the scattered radiation removal grid by the image processing means. And a selection unit that enables selection of whether or not to remove the image.
【００１７】本発明の他の観点によれば、（ａ）散乱線
除去グリッドを使用してＸ線を曝射し、被写体のＸ線透
過分布を２次元的に所望の間隔でサンプリングして画像
を取得するステップと、（ｂ）前記取得した画像中から
前記散乱線除去グリッドに起因する縞を除去するか否か
を選択するステップと、（ｃ）前記選択に応じて前記散
乱線除去グリッドに起因する縞を画像処理により除去す
るステップとを有することを特徴とする画像取得方法が
提供される。According to another aspect of the present invention, (a) X-rays are radiated using a scattered radiation removal grid, and the X-ray transmission distribution of a subject is sampled two-dimensionally at a desired interval. (B) selecting whether or not to remove fringes caused by the scattered radiation removal grid from the acquired image; and (c) selecting the scattered radiation removal grid in accordance with the selection. Removing the resulting fringe by image processing.
【００１８】本発明によれば、画像中から散乱線除去グ
リッドに起因する縞の除去を行うか否かの選択が可能で
あるので、必要な場合にのみ散乱線除去グリッドに起因
する縞の除去を行い、適切な画像を得ることができる。According to the present invention, it is possible to select whether or not to remove the fringes caused by the scattered radiation removal grid from the image. Therefore, it is possible to remove the fringes caused by the scattered radiation removal grid only when necessary. And an appropriate image can be obtained.
【００１９】[0019]
【発明の実施の形態】以下、本発明の実施形態を、実施
例に沿って図面に基づいて説明する。 （第１の実施例）本発明の第１の実施例による画像取得
装置では、取得した画像を再生する場合、グリッド像が
存在し、しかも画像成分にある程度干渉していても観察
者に対しては違和感を与え難い周波数の縞になるようグ
リッド鉛の間隔を設定することで、従来技術による問題
点を解決した。DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS The embodiments of the present invention will be described below with reference to the accompanying drawings. (First Embodiment) In an image acquisition apparatus according to a first embodiment of the present invention, when reproducing an acquired image, a grid image exists and even if it interferes to some extent with an image component, it is possible to observe an observer. Solved the problems of the prior art by setting the interval between the grid leads so that the stripes had a frequency that would not give a strange feeling.
【００２０】以下解説を行う。サンプリングされたデジ
タル画像において、グリッドによる画像のアーチファク
トを以下の２つに分類する。 （Ａ）低空間周波数のグリッド像 （Ｂ）低空間周波数のグリッド像をサンプリングしたこ
とによるビート像The following is a description. In a sampled digital image, image artifacts due to the grid are classified into the following two types. (A) Low spatial frequency grid image (B) Beat image obtained by sampling low spatial frequency grid image
【００２１】図７は、上述の（Ａ）、（Ｂ）のアーチフ
ァクトを説明するためにアーチファクトの様子を模式的
に示したものである。以下説明上、グリッド縞の２倍高
調波以上は、解像しないものとして、基本波長のみで説
明する。FIG. 7 schematically shows the state of the artifact to explain the above-mentioned artifacts (A) and (B). In the following description, only the fundamental wavelength will be described on the assumption that the second or higher harmonic of the grid stripe is not resolved.
【００２２】図７（１−ａ）は、上述の（Ａ）のアーチ
ファクトの一例を周波数領域（正領域のみ）での様子を
示したものであり、６１で示した領域が画像の成分の領
域であり、実質的な最高周波数をｆｉで示している。６
２がサンプリングにより数学的に生じる画像の高調波成
分である。ここでサンプリング周波数はｆｓであり、ナ
イキスト周波数ｆｎｑ＝ｆｓ／２である。６３の線スペ
クトルが周波数ｆｇをもつグリッド縞の成分であり、図
のように画像の成分６１と重なり、画像に損傷を与えて
いるとともに低周波であるため、観察者にとって縞模様
が違和感を与える。FIG. 7 (1-a) shows an example of the above-mentioned artifact (A) in the frequency domain (only the positive domain). The area 61 is the area of the image component. And the substantially highest frequency is indicated by fi. 6
2 is a harmonic component of the image that is mathematically generated by sampling. Here, the sampling frequency is fs, and the Nyquist frequency fnq = fs / 2. The line spectrum of 63 is a grid fringe component having a frequency fg, which overlaps with the image component 61 as shown in the figure, which damages the image and has a low frequency. .
【００２３】図７（１−ｂ）はサンプリングの様子を１
次元的に示したものであり、●（黒丸）がサンプリング
された点を示す。すなわち、この事からグリッド縞の周
波数ｆｇと画像の持つ周波数ｆｉ（ｆｉ＜ｆｓ／２）の
関係は、 ｆｇ ＞ ｆｉ ・・・（１） でなければならないことが分かる。FIG. 7 (1-b) shows the sampling state as 1.
This is shown in a dimensional manner, where ● (black circle) indicates a sampled point. That is, it can be understood from this that the relationship between the frequency fg of the grid stripe and the frequency fi (fi <fs / 2) of the image must be fg> fi (1).
【００２４】図７（２−ａ）は、上述の（Ａ）のアーチ
ファクトの他の例を周波数軸上（正領域のみ）での様子
を示したものであり、６１で示した領域が画像の成分の
領域であり、最高周波数をｆｉで示している。サンプリ
ング周波数ｆｓの１／２であるナイキスト周波数ｆｎｑ
よりグリッドの周波数ｆｇが大きいため、エリアシング
であるｆｓ−ｆｇの周波数成分６４が現れている。さら
に、図７（２−ａ）の場合、このエリアシングの周波数
成分６４が画像成分６１に重なり、画像に損傷を与えて
いるとともに観察者にとって縞模様が違和感を与える。
図７（２−ｂ）はサンプリングの様子を１次元的に示し
たものであり、●（黒丸）がサンプリングされた点を示
し、破線は、サンプリングされる前の信号形状を示す。
すなわち画像成分に損傷を与えないためには、グリッド
縞の周波数ｆｇと画像の持つ周波数ｆｉ（ｆｉ＜ｆｓ／
２）の関係を、 ｆｓ−ｆｇ ＞ ｆｉ ・・・（２） としなければならないことが分かる。FIG. 7 (2-a) shows another example of the above-described artifact (A) on the frequency axis (only the positive region), and the region indicated by 61 is the image. This is a component area, and the highest frequency is indicated by fi. Nyquist frequency fnq which is 1/2 of sampling frequency fs
Since the frequency fg of the grid is larger, a frequency component 64 of fs-fg, which is aliasing, appears. Further, in the case of FIG. 7 (2-a), the aliasing frequency component 64 overlaps the image component 61, causing damage to the image and giving the viewer a sense of discomfort due to the stripe pattern.
FIG. 7 (2-b) shows the state of sampling in a one-dimensional manner. A black dot indicates a sampled point, and a broken line indicates a signal shape before being sampled.
That is, in order not to damage the image components, the frequency fg of the grid stripe and the frequency fi (fi <fs /
It can be seen that the relationship of 2) must be fs-fg> fi (2).
【００２５】（１）式と（２）式によりグリッド縞の周
波数ｆｇは画像の持つ最高周波数成分ｆｉとサンプリン
グ周波数ｆｓとの関係で、 ｆｉ ＜ ｆｇ ＜ ｆｓ−ｆｉ （ｆｉ＜ｆｓ／２）・・・（３） という関係式の中にいれなければならない。しかし、こ
の場合でも上述の（Ｂ）のアーチファクトが発生してし
まう場合がある。According to the equations (1) and (2), the frequency fg of the grid stripe is determined by the relationship between the highest frequency component fi of the image and the sampling frequency fs: fi <fg <fs-fi (fi <fs / 2).・ ・ (3) must be in the relational expression. However, even in this case, the above-described artifact (B) may occur.
【００２６】図７（３−ａ）は、（３）式の条件を満た
すような空間周波数領域にグリッド縞の周波数を設定し
た場合の空間周波数上での様子を示したものである。図
７（３−ｂ）は、サンプリングの様子を１次元的に示し
たものであり、●（黒丸）がサンプリングされた点を示
し、破線は、サンプリングされる前の信号形状を示す。FIG. 7 (3-a) shows a state on the spatial frequency when the grid stripe frequency is set in the spatial frequency region satisfying the condition of the equation (3). FIG. 7 (3-b) shows the state of sampling in a one-dimensional manner. A black dot indicates a sampled point, and a broken line indicates a signal shape before sampling.
【００２７】図７（３−ａ）では、画像周波数成分６１
とグリッド像成分６３が重ならないため、画像へ与える
損傷はなく、違和感もないように見受けられる。しか
し、図７（３−ｂ）で示すような、ビート状の振幅の変
動成分が画像情報の周波数成分に重なる場合には、画像
の観察者には画像に損傷を与えるアーチファクトとして
認識され観察者にとって縞模様が違和感を与える。この
振幅変動の周波数は｜ｆｓ／２−ｆｇ｜である。画像へ
の何らかの非線形変換が行われると、実際に画像成分へ
損傷を与えるアーチファクトになる可能性もある。In FIG. 7 (3-a), the image frequency component 61
And the grid image component 63 do not overlap with each other, so that there is no damage to the image and there is no uncomfortable feeling. However, when the beat-like amplitude fluctuation component overlaps the frequency component of the image information as shown in FIG. 7 (3-b), the image is recognized as an artifact that damages the image to the observer of the image. The stripes give a sense of incongruity. The frequency of this amplitude variation is | fs / 2−fg |. Any non-linear transformation to the image can actually result in artifacts that damage image components.
【００２８】理想的にはすべてのアーチファクトの周波
数が画像の周波数ｆｉ外にあれば問題はない。この場
合、ｆｉを画像の最高周波数であるとしたが、画像を再
生した場合（ハードコピー、モニター表示など）に観察
者に対し違和感を与えなくなる周波数であるととらえて
も構わない。そこで、（３）式の関係と｜ｆｓ／２−ｆ
ｇ｜＞ｆｉの両方の関係を同時に満たさなければならな
い。これは、たとえばｆｓ＞ｆｇ＞ｆｓ／２の関係の位
置にｆｇを設定すると仮定すると、 ｆｉ＋ｆｓ／２ ＜ ｆｇ ＜ ｆｓ−ｆｉ ・・・（４） の条件の関係式が出る。この条件を満たすには、さらに
必然的にｆｉ＜ｆｓ／４という条件も満たさなければな
らない厳しいものである。（４）式においてｆｉ＜ｆｓ
／４の関係が満たされないと、左側不等号と右側不等号
の領域の合い重なる部分がない。したがって、理想的に
グリッドの影響が無い条件は、画像の周波数帯域が、ナ
イキスト周波数の半分以下であり、グリッド縞の周波数
が、ナイキスト周波数の半分の近辺にしなければならな
いという結果を得る。この条件は、取得画像が決まって
いれば、サンプリング周波数を取得画像の最大空間周波
数のすくなくとも４倍以上に設定しなければならないと
いうことを意味している。Ideally, there is no problem if the frequencies of all artifacts are outside the frequency fi of the image. In this case, fi is assumed to be the highest frequency of the image, but may be considered to be a frequency at which the observer does not feel uncomfortable when the image is reproduced (hard copy, monitor display, etc.). Therefore, the relation of equation (3) and | fs / 2−f
Both relations g |> fi must be satisfied simultaneously. For example, assuming that fg is set at the position of the relationship of fs>fg> fs / 2, the relational expression of the condition of fi + fs / 2 <fg <fs-fi (4) comes out. In order to satisfy this condition, the condition fi <fs / 4 must be satisfied. In equation (4), fi <fs
If the relationship of / 4 is not satisfied, there is no overlap between the left and right inequality regions. Therefore, under the condition that there is no influence of the grid ideally, the result is that the frequency band of the image is equal to or less than half of the Nyquist frequency, and the frequency of the grid stripe must be around half of the Nyquist frequency. This condition means that if the acquired image is determined, the sampling frequency must be set to at least four times the maximum spatial frequency of the acquired image.
【００２９】しかし、上述の論点で考慮していないの
は、ビート状の変動成分のパワーである。ビート状の変
動成分のパワーとグリッド縞のパワーを比較した結果、
および実際にグリッドを用いた画像を観察した結果か
ら、（４）式の関係をより緩やかなものにできることを
本発明者は考案し、実証した。However, what is not considered in the above-mentioned point is the power of the beat-like fluctuation component. As a result of comparing the power of the beat-like fluctuation component and the power of the grid stripe,
The present inventor has devised and verified that the relationship of the expression (4) can be made more gradual from the result of observation of an image actually using a grid.
【００３０】ビートである振幅変動は縞の周波数がナイ
キスト周波数から離れれば離れるほど高周波になり、観
察できなくなる。以下では、この縞の周波数がナイキス
トからどの程度離れれば、観察者にとって観察しにくく
なるかを考察する。The amplitude fluctuation, which is the beat, becomes higher as the frequency of the fringes becomes farther from the Nyquist frequency, and cannot be observed. In the following, the extent to which the frequency of the fringes is away from Nyquist will make it difficult for an observer to observe.
【００３１】サンプリングを行えば、必ずナイキスト周
波数を中心として鏡像関係の位置に線スペクトルが存在
し、それとの間に必ずビートが発生している。しかし、
図７（１−ｂ）及び図７（２−ｂ）ではビートが発生し
ているはずにもかかわらず、一見したところそのような
ものは存在しない。この考えられる一つの理由は周波数
であり、図７（２−ａ）もしくは図７（３−ａ）ではナ
イキスト周波数を中心として鏡像関係にある２つのスペ
クトルの距離（周波数）が長く、ビート周波数が高周波
になる。この時には、基本の正弦波のスペクトルはナイ
キスト周波数より充分低いところに存在し、この両者の
周波数およびパワーの差が大きいため、よりパワーの強
い基本正弦波のみが観察者に強く認識されていると考え
られる。When sampling is performed, a line spectrum always exists at a position related to a mirror image with respect to the Nyquist frequency, and a beat is always generated between the line spectrum and the line spectrum. But,
In FIG. 7 (1-b) and FIG. 7 (2-b), such a beat does not exist at first glance, although a beat should have occurred. One possible reason for this is frequency. In FIG. 7 (2-a) or FIG. 7 (3-a), the distance (frequency) between two spectra that are in a mirror image relationship around the Nyquist frequency is long, and the beat frequency is high. High frequency. At this time, the spectrum of the basic sine wave exists sufficiently lower than the Nyquist frequency, and the difference between the frequency and the power of the two is large, so that only the fundamental sine wave with higher power is strongly recognized by the observer. Conceivable.
【００３２】ビートの起こる様子を数式で示す。空間周
波数ｆｇの正弦波をサンプリング周波数ｆｓでサンプリ
ングした場合を考える。この時実際に使用するグリッド
を考えてｆｇ＞ｆｓ／２と仮定する。ただし、これは必
要な条件ではない。この時、図８（ａ）で示すように、
正弦波の周波数に相当する線スペクトルペアが発生す
る。線スペクトルのピークをａ／２と仮定し、さらにも
とのグリッド像をａ×ｃｏｓ（２πｆｇｘ）と仮定す
る。周波数ｆｇの余弦波を周波数ｆｓでサンプリングし
た場合には、ｆｓ以下には２つの余弦波が現れる。The manner in which a beat occurs is shown by mathematical expressions. Consider a case where a sine wave having a spatial frequency fg is sampled at a sampling frequency fs. At this time, it is assumed that fg> fs / 2 in consideration of a grid actually used. However, this is not a necessary condition. At this time, as shown in FIG.
A line spectrum pair corresponding to the frequency of the sine wave is generated. Assume that the peak of the line spectrum is a / 2 and that the original grid image is a × cos (2πfgx). When a cosine wave of frequency fg is sampled at frequency fs, two cosine waves appear below fs.
【００３３】 ｇ（ｘ）＝ａ｛ｃｏｓ（２πｆｇｘ） ＋ｃｏｓ（２π（ｆｓ−ｆｇ）ｘ）｝ ＝２ａ×ｃｏｓ（２π（ｆｇ−ｆｓ／２）ｘ） ×ｃｏｓ（２πｆｓｘ／２） ・・・（５）G (x) = a {cos (2πfgx) + cos (2π (fs−fg) x)} = 2a × cos (2π (fg−fs / 2) x) × cos (2πfsx / 2) (5)
【００３４】（５）式はビートをあらわし、２つの正弦
波の差の周波数の正弦波によって変調された振幅変調に
相当する。（ｆｇ／２−ｆｓ）が０ではないが小さな値
になると低周波の不安定な振幅変動（ビート）を起こ
す。Equation (5) represents a beat and corresponds to amplitude modulation modulated by a sine wave having a frequency of a difference between two sine waves. When (fg / 2-fs) is not 0 but becomes a small value, unstable low-frequency amplitude fluctuation (beat) occurs.
【００３５】シャノンのサンプリング定理によれば、ナ
イキスト周波数以下でサンプリングされたものは理想的
なフィルタリング手段（ナイキスト周波数以下を通過さ
せるフィルタ）を用いれば完全に復元でき、そこには線
スペクトルペアによるビートは存在しないはずである。
すなわち、図７（３−ｂ）で観察者がビートを強く認識
した理由は、そのフィルタリング手段によるものである
可能性が高い。According to Shannon's sampling theorem, a signal sampled below the Nyquist frequency can be completely restored by using ideal filtering means (a filter that passes the signal below the Nyquist frequency). Should not exist.
That is, the reason why the observer strongly recognized the beat in FIG. 7 (3-b) is likely to be due to the filtering means.
【００３６】通常ではサンプルした点を直線などでつな
ぎ、サンプリングした信号を復元しようとしている。こ
れはサンプリング定理の理想的なフィルタ手段とは異な
るものである。Usually, the sampled points are connected by a straight line or the like to restore the sampled signal. This is different from the ideal filter means of the sampling theorem.
【００３７】しかし、一般的な感覚をもつ人もしくは表
示装置は、サンプリング定理で言うような理想的なフィ
ルタ（ｓｉｎｃ関数をカーネルとするようなコンボリュ
ーション）はかけず、近傍の点を直接直線でつなぐよう
な補間をおこなうのであり、画像として観察する場合も
同様である。つまり、このような一般的な感覚とサンプ
リング定理とのずれがビートとなってあらわれているこ
とに気付く。However, a person or a display device having a general sensation does not apply an ideal filter (convolution using a sinc function as a kernel) as in the sampling theorem, but directly points nearby points by a straight line. Interpolation is performed as if it were connected, and the same applies to observation as an image. In other words, we notice that such a gap between the general sense and the sampling theorem appears as a beat.
【００３８】直線で補間することは図８（ｂ）の７１で
示すような特性を持つフィルタリングを行うことに相当
する。Interpolation by a straight line is equivalent to performing filtering having a characteristic indicated by reference numeral 71 in FIG. 8B.
【００３９】図８（ｂ）の７１でフィルタの形状ｓ
（ｆ）は、 ｓ（ｆ）＝ｓｉｎ2（πｆ／ｆｓ）／（πｆ／ｆｓ）2 ・・・（６） であり、このフィルタをかけられた後の鏡像にある線ス
ペクトルペアの高さをｃ／２およびｄ／２とする。The filter shape s at 71 in FIG.
(F) is s (f) = sin 2 (πf / fs) / (πf / fs) 2 ... (6), and the height of the line spectrum pair in the mirror image after being filtered. Are c / 2 and d / 2.
【００４０】この時の正弦波の和を（５）式と同様にあ
らわす。 ｇ（ｘ）＝ｄ×ｃｏｓ（２πｆｇｘ） ＋ｃ×ｃｏｓ（２π（ｆｓ−ｆｇ）ｘ） ＝２ｄ×ｃｏｓ（２π（ｆｇ／２−ｆｓ）ｘ） ×ｃｏｓ（２πｆｓｘ／２） ＋（ｃ−ｄ）×ｃｏｓ（２π（ｆｓ−ｆｇ）ｘ）・・・（７）The sum of the sine waves at this time is expressed in the same manner as in equation (5). g (x) = d × cos (2πfgx) + c × cos (2π (fs−fg) x) = 2d × cos (2π (fg / 2-fs) x) × cos (2πfsx / 2) + (cd) ) × cos (2π (fs−fg) x) (7)
【００４１】（７）式の第１項がビートの成分をあらわ
し、第２項が通常の正弦波の成分をあらわしている。ま
た、（５）式より以下の関係がある。The first term of the equation (7) represents a beat component, and the second term represents a normal sine wave component. Further, the following relationship is obtained from equation (5).
【００４２】[0042]
【数１】 (Equation 1)
【００４３】ここで、（６）式においてビート成分であ
る第１項のパワーより通常の正弦波成分である第２項の
パワーが上回れば、観察者としてはビート成分が認識し
にくくなると思われる。Here, if the power of the second term, which is a normal sine wave component, exceeds the power of the first term, which is a beat component, in equation (6), it will be difficult for an observer to recognize the beat component. .
【００４４】以上の仮説のもとに（６）式の第１項と第
２項のパワーの比をとると、以下のようになる。 （ｃ−ｄ）2／２ｄ2 ＞ １ （第２項のパワーが第１
項を上回る条件）When the power ratio of the first and second terms in the equation (6) is calculated based on the above hypothesis, the following is obtained. (Cd) 2 / 2d 2 > 1 (the power of the second term is the first
Condition that exceeds the term)
【００４５】これを変形すると、次のようになる。 ｃ／ｄ ＞ ２1/2＋１ これに、（７）式を代入すると、次式になる。When this is modified, the following is obtained. c / d> 21/2 + 1 By substituting equation (7) into this, the following equation is obtained.
【００４６】[0046]
【数２】 (Equation 2)
【００４７】（９）式の結果より、サンプリングされる
正弦波の周波数がサンプリングする周波数の６０．８％
より大きければ（エリアシングは起こす）、ビートが観
測されにくいサンプリングが行われることがわかる。こ
の場合、グリッドの縞情報はナイキスト周波数（ｆｓ／
２）の８０％以下の成分として現れる。これは、サンプ
リング周波数の４０％以下の成分と等価である。したが
って、グリッド縞がナイキスト周波数の８０％以下であ
れば、サンプリングした場合のビートが目立たず安定し
た、縞模様が観察できることになる。From the result of equation (9), the frequency of the sine wave to be sampled is 60.8% of the frequency to be sampled.
If it is larger (aliasing occurs), it is understood that sampling in which a beat is hardly observed is performed. In this case, the grid stripe information is the Nyquist frequency (fs /
It appears as a component of 80% or less of 2). This is equivalent to a component of 40% or less of the sampling frequency. Therefore, if the grid stripes are 80% or less of the Nyquist frequency, the beat when sampled is inconspicuous and a stable stripe pattern can be observed.
【００４８】以上の考察は、グリッドに起因するナイキ
スト周波数以下に現れる縞（グリッド縞）の上限の空間
周波数を規定するものであり、具体的にはナイキスト周
波数の８０％以下（サンプリング周波数の４０％以下）
にグリッド縞の周波数を設定するというものであるが、
実際には、低ければ低いほどよいというものではなく、
前述の縞自体によるアーチファクトすなわち、画像成分
に縞情報が重なってしまうことは許されないため、画像
成分の最高周波数が、グリッド縞の下限の周波数より小
さい必要がある。The above consideration defines the upper limit spatial frequency of fringes (grid fringes) appearing below the Nyquist frequency due to the grid. Specifically, the spatial frequency is 80% or less of the Nyquist frequency (40% of the sampling frequency). Less than)
Is to set the frequency of the grid stripes,
In fact, the lower the better, the better.
Since the above-described artifact due to the fringes itself, that is, the fringe information is not allowed to overlap with the image component, the highest frequency of the image component needs to be lower than the lower limit frequency of the grid fringe.
【００４９】一般的に、画像を代表する信号の最高周波
数成分は正確に規定できるものではない。画像の場合の
評価基準の例を以下に列挙する。 ・取得した画像を表示装置もしくは記録装置で再現した
場合に、観察者に充分であると認識させるような周波
数。 ・最高周波数を想定したとき、その１．５〜２倍を最高
周波数と考えて、サンプリングピッチを決定する（中溝
他著、「計数・測定」（培風館）参照）。In general, the highest frequency component of a signal representing an image cannot be accurately defined. Examples of evaluation criteria for images are listed below. A frequency that, when reproduced on a display or recording device, causes an observer to perceive that it is sufficient. Assuming the highest frequency, 1.5 to 2 times the highest frequency is considered the highest frequency, and the sampling pitch is determined (see Nakamizo et al., "Counting and Measurement" (Baifukan)).
【００５０】特に、後者の条件はよく用いられるもので
あり、一般に必要とする空間周波数がナイキスト周波数
の６０％以下（サンプリング周波数の３０％以下）にな
るように、サンプリングピッチを設定する場合が多い。
すなわち、サンプリング周波数の３０％以上にグリッド
縞情報が現れるように設定する。In particular, the latter condition is often used, and the sampling pitch is often set so that the generally required spatial frequency is 60% or less of the Nyquist frequency (30% or less of the sampling frequency). .
That is, it is set so that grid stripe information appears at 30% or more of the sampling frequency.
【００５１】たとえば、サンプリングピッチが０．１ｍ
ｍ（ｆｓ＝１０ｃｙｃ／ｍｍ）である場合を想定する。
ナイキスト周波数は５ｃｙｃ／ｍｍであり、通常用いる
画像成分の周波数は上述の一般的な条件からサンプリン
グ周波数の３０％以下すなわち３ｃｙｃ／ｍｍであるた
め、３ｃｙｃ／ｍｍが縞情報の下限であり、上限はサン
プリング周波数の４０％、すなわち４ｃｙｃ／ｍｍであ
る。したがって、グリッド縞の周波数は３〜４ｃｙｃ／
ｍｍの範囲に来るように設定すれば、グリッド縞模様に
はビート成分が目立たず、観察にも邪魔にならないよう
になる。For example, if the sampling pitch is 0.1 m
Assume that m (fs = 10 cyc / mm).
The Nyquist frequency is 5 cyc / mm, and the frequency of a normally used image component is 30% or less of the sampling frequency, that is, 3 cyc / mm from the above general conditions, so that 3 cyc / mm is the lower limit of the stripe information, and the upper limit is the upper limit. It is 40% of the sampling frequency, that is, 4 cyc / mm. Therefore, the frequency of the grid stripe is 3 to 4 cyc /
If it is set so as to fall within the range of mm, the beat component is not conspicuous in the grid stripe pattern, and does not disturb the observation.
【００５２】以上の条件は、グリッド縞模様の周波数の
決定であり、この数値に縞模様が来るグリッド本体の鉛
の周波数はサンプリングされることにより異なったもの
になる。このとき、グリッド本体の周波数範囲ＲＧ［ｃ
ｙｃ／ｍｍ］は、図９に示すように、次の範囲になる。
ただし、ｎは０以上の整数である。５（２ｎ＋１）−２
〜５（２ｎ＋１）−１もしくは５（２ｎ＋１）＋１〜
（２ｎ＋１）＋２ ；ｎ＝０，１，２，・・・となる。
（図９参照）The above condition is the determination of the frequency of the grid stripe pattern, and the lead frequency of the grid body where the stripe pattern comes to this numerical value becomes different due to sampling. At this time, the frequency range RG [c
yc / mm] is in the following range as shown in FIG.
Here, n is an integer of 0 or more. 5 (2n + 1) -2
~ 5 (2n + 1) -1 or 5 (2n + 1) + 1 ~
(2n + 1) +2; n = 0, 1, 2,...
(See Fig. 9)
【００５３】さらに、この式は、サンプリングピッチが
０．１ｍｍ（ｆｓ＝１０ｃｙｃ／ｍｍ）の時の計算値で
あるが、一般的にサンプリング周波数ｆｓ［ｃｙｃ／ｍ
ｍ］（サンプリングピッチ１／ｆｓ［ｍｍ］）の場合に
は、グリッド本体の周波数範囲は、 ｆｓ（ｎ＋０．３）〜ｆｓ（ｎ＋０．４） もしくはｆｓ（ｎ＋０．６）〜ｆｓ（ｎ＋０．８）［ｃ
ｙｃ／ｍｍ］ と表される。Further, although this equation is a calculated value when the sampling pitch is 0.1 mm (fs = 10 cyc / mm), generally, the sampling frequency fs [cyc / m
m] (sampling pitch 1 / fs [mm]), the frequency range of the grid body is fs (n + 0.3) to fs (n + 0.4) or fs (n + 0.6) to fs (n + 0.8). ) [C
yc / mm].
【００５４】ここでさらに、グリッドの本来の目的であ
る散乱線除去の性能、および平面センサの解像力などを
考慮して、上記の範囲の中からグリッド本体の周波数
（グリッド本数）が決定される。通常では、センサの解
像力を考え、第２高調波を解像しにくく、散乱線除去率
の高い６〜７ｃｙｃ／ｍｍのグリッド縞周波数が選定で
きる。Here, the frequency (the number of grids) of the grid body is determined from the above range in consideration of the performance of scattered radiation removal, which is the original purpose of the grid, and the resolution of the flat sensor. Normally, in consideration of the resolution of the sensor, it is difficult to resolve the second harmonic, and a grid fringe frequency of 6 to 7 cyc / mm with a high scattered radiation removal rate can be selected.
【００５５】本実施例では、使用するグリッドの周波数
を上記の計算から選定し、グリッド縞を除去しなくても
観察者の慣れもしくは使用目的によりある程度満足でき
るように選定することで対応し、問題点を克服した。In the present embodiment, the frequency of the grid to be used is selected from the above calculation, and is selected so as to be somewhat satisfied depending on the familiarity of the observer or the purpose of use without removing grid stripes. Overcame the point.
【００５６】グリッド縞は空間周波数が固定しているた
め、フィルタリング作業によりある程度除去可能であ
る。そこで、上述のようにグリッド縞周波数を設定する
ことにより、グリッド縞を除去する操作において、完全
な除去が不可能であっても、グリッド縞の強度が低減さ
れると、さらに観察者への影響が最小限にとどめられ
る。Since the spatial frequency of the grid stripe is fixed, it can be removed to some extent by a filtering operation. Therefore, by setting the grid fringe frequency as described above, even if it is impossible to completely remove the grid fringes, if the intensity of the grid fringes is reduced, the influence on the observer is further reduced Is minimized.
【００５７】第１の実施例の模式図を図１に示す。同図
はテーブル上に横たわる人体を撮影するためのシステム
の模式図であり、１はＸ線管球、２は被写体である人
体、１１は散乱Ｘ線除去用のグリッドを示し、被写体２
から発生する散乱Ｘ線を除去し、着脱可能な構成になっ
ている。３はＸ線強度分布（Ｘ線透過分布）を電荷分布
に変換し、２次元的に所望の間隔でサンプリングして順
次出力するＸ線センサパネル、５はアナログ／デジタル
変換器、４はＸ線の曝射と画像取得のタイミングを制御
するコントローラ、６は画像を一旦記憶するメモリをあ
らわす。Ｘ線センサパネル３には画素ごとにオフセット
およびゲインのばらつきがある。このばらつきを補正す
るために８で示すメモリにはＸ線を曝射しないで取得し
た画像であるオフセット値、９には被写体２およびグリ
ッド１１がない状態で取得されたゲイン値を対数変換し
たものを記憶しておく。７は対数変換を行う変換装置で
具体的には参照テーブル（ルックアップテーブル）であ
る。取得された人体の画像は８のオフセット値を減じた
（除去した）後、対数変換され、９のゲインとの差分
（割り算）を行い、ゲイン値のばらつきを補正したＸ線
の強度分布画像を得て、１０で示すメモリに一旦記憶さ
れる。その後、この記憶された画像を取り出し、画像の
保存、画像処理、画像表示、ハードコピーなどの行為が
行われ、診断などに用いられる。FIG. 1 is a schematic diagram of the first embodiment. FIG. 1 is a schematic diagram of a system for photographing a human body lying on a table, 1 is an X-ray tube, 2 is a human body as a subject, 11 is a grid for removing scattered X-rays,
It is configured to remove scattered X-rays generated from and remove it. An X-ray sensor panel 3 converts an X-ray intensity distribution (X-ray transmission distribution) into a charge distribution, samples two-dimensionally at a desired interval, and sequentially outputs the sampled data. The controller 6, which controls the timing of the exposure of the image and the acquisition of the image, represents a memory for temporarily storing the image. The X-ray sensor panel 3 has variations in offset and gain for each pixel. In order to correct this variation, an offset value which is an image obtained without exposing X-rays to the memory indicated by 8 and a gain value obtained in the absence of the subject 2 and the grid 11 are logarithmically converted to 9 are obtained. Is stored. Reference numeral 7 denotes a conversion device that performs logarithmic conversion, and is specifically a reference table (lookup table). The acquired image of the human body is subjected to logarithmic conversion after subtracting (removing) the offset value of 8 and performing a difference (division) with the gain of 9 to obtain an X-ray intensity distribution image in which the variation of the gain value is corrected. Then, it is temporarily stored in the memory indicated by 10. After that, the stored image is taken out, and actions such as image storage, image processing, image display, and hard copy are performed, and used for diagnosis and the like.
【００５８】１２で示されるブロックは、画像処理（フ
ィルタリング）によりグリッド縞を除去する画像処理手
段（フィルタリング手段）を示し、１０に記憶された画
像を用いて空間フィルタリングを施し、グリッド縞成分
を除去する。１３で示される機構は、外部からのグリッ
ド除去操作実効設定手段（機構）１６又は操作パネル２
０の出力に応じて、操作者が１２のフィルタリングを行
うかどうかを選択する選択手段（スイッチ）である。ス
イッチ１３の選択に応じて、１２のブロックの動作をス
キップできるように信号の流れを変えることができる。Block 12 indicates image processing means (filtering means) for removing grid fringes by image processing (filtering), and performs spatial filtering using the image stored in 10 to remove grid fringe components. I do. The mechanism indicated by 13 is an external grid removal operation effective setting means (mechanism) 16 or the operation panel 2
Selection means (switch) for the operator to select whether to perform twelve filterings according to the output of zero. According to the selection of the switch 13, the signal flow can be changed so that the operation of the 12 blocks can be skipped.
【００５９】ここで、３で示すＸ線センサパネルは２次
元空間的に複数の画素をもち、０．１ｍｍピッチで縦横
に分布しており、この機構で２次元的かつ離散的なサン
プリングが行われる。前述のグリッドの縞の周波数の設
定で述べたように、１１のグリッド縞の周波数（グリッ
ド本数）は６〜７ｃｙｃ／ｍｍにしてあるため、観察者
はグリッド縞をフィルタリングで除去しなくても違和感
なく画像を観察可能である。Here, the X-ray sensor panel indicated by 3 has a plurality of pixels in a two-dimensional space and is distributed vertically and horizontally at a pitch of 0.1 mm, and two-dimensional and discrete sampling is performed by this mechanism. Will be As described in the setting of the frequency of the grid stripes described above, since the frequency (the number of grids) of the 11 grid stripes is set to 6 to 7 cyc / mm, the observer can feel uncomfortable without removing the grid stripes by filtering. Without observing the image.
【００６０】１０で示すメモリに記憶された画像のグリ
ッド縞に直交する方向の１次元の振幅スペクトルの様子
を模式的に示したのが図２（ａ）である。図２（ａ）
で、３２が画像成分のスペクトル、３１がグリッド縞結
成分のスペクトルであり、ノイズを無視した実質的なス
ペクトル形状を示している。３１のグリッド縞成分が存
在しても、観察者には安定したこの周波数成分しか観察
されず、これ以外のビートに関する成分はなくなるた
め、慣れもしくは縞の認識により観察者にはグリッド縞
の存在が比較的邪魔にならない。FIG. 2A schematically shows the state of a one-dimensional amplitude spectrum in a direction orthogonal to the grid stripes of the image stored in the memory indicated by reference numeral 10. FIG. 2 (a)
, 32 is the spectrum of the image component, and 31 is the spectrum of the grid stripe component, and shows a substantial spectrum shape ignoring noise. Even if there are 31 grid fringe components, only the stable frequency components are observed by the observer, and there are no other components related to beats. Relatively out of the way.
【００６１】しかし、後段の画像処理の関係もしくは画
像再生機構の関係で操作者もしくは観察者がこの画像か
らグリッド縞を除去したいと望めば、操作パネル２０か
らグリッド除去選択を行う。However, if the operator or the observer desires to remove grid stripes from this image due to the relationship between the subsequent image processing or the image reproducing mechanism, the operator selects a grid removal from the operation panel 20.
【００６２】図２（ｂ）はグリッド除去フィルタリング
を行った後のスペクトルの様子を模式的に示したもので
あり、３３がフィルタ特性の例、３４がフィルタ通過後
の画像スペクトルである。フィルタは空間的な特性を安
定にするため、急峻なものはかけられない。したがっ
て、３１のグリッド縞成分を除去しようとすると必然的
に画像成分の一部も応答が低下してしまう。操作者はこ
の事を把握した上で、グリッド除去を行うかどうかの選
択を行う。本実施例で、９に入るゲイン値をしめす画像
取得時には、グリッド１１を設置したままで取得しても
グリッドの周波数は不変であるので構わない。FIG. 2B schematically shows the state of the spectrum after grid removal filtering is performed. Reference numeral 33 denotes an example of a filter characteristic, and reference numeral 34 denotes an image spectrum after passing through the filter. Since the filter stabilizes the spatial characteristics, a steep filter cannot be applied. Therefore, when trying to remove the 31 grid stripe components, the response of some of the image components is inevitably reduced. After grasping this, the operator selects whether or not to remove the grid. In this embodiment, at the time of acquiring an image showing a gain value falling within 9, even if the image is acquired with the grid 11 installed, the frequency of the grid does not change.
【００６３】図３は、本実施例をソフトウエアで行う場
合のフローチャートである。同図で処理ブロック（ステ
ップ）をＣ１〜Ｃ１１までの動作に分割する。ブロック
Ｃ１の動作はゲイン値の画像を取得するもので、被写体
のない状態でＸ線を曝射し取得した画像を画像Ａとす
る。それを、ブロックＣ２では対数変換して、画像Ｂと
する。ブロックＣ３ではオフセット値を取得するもの
で、Ｘ線を曝射しない状態で画像を取得し、画像Ｃとす
る。次のブロックＣ４が実際の被写体の画像を取得する
ブロックであり、グリッドを設置し、被写体がある状態
でＸ線を曝射し取得した画像を画像Ｄとする。ブロック
Ｃ５では画像Ｄから画像Ｃを減じ、オフセットの補正が
なされた画像Ｅを得る。ブロックＣ６では画像Ｅを対数
変換し、画像Ｆとし、ブロックＣ７では画像Ｆから画像
Ｂを減じ、ゲインの補正がなされた被写体画像Ｇを得
る。次のブロックＣ８では、操作パネル２０による操作
者の操作（指示）に応じて、グリッド除去を行うかどう
かによって分岐が行われる。ブロックＣ９ではグリッド
除去をおこなう指示であるので、画像Ｇに対しグリッド
除去のフィルタリング処理を行い、グリッドが除去され
た被写体画像Ｈを得る。次のブロックＣ１０では、被写
体画像Ｈを出力する。グリッド除去を指示されなけれ
ば、ブロックＣ１１で被写体画像Ｇをそのまま出力す
る。FIG. 3 is a flowchart when the present embodiment is performed by software. In the figure, a processing block (step) is divided into operations C1 to C11. The operation of block C1 is to acquire an image of a gain value, and assume that an image acquired by irradiating X-rays without a subject is an image A. It is logarithmically converted into an image B in block C2. In block C3, an offset value is acquired, and an image is acquired without X-ray irradiation, which is referred to as an image C. The next block C4 is a block for acquiring an image of an actual subject, and an image obtained by arranging a grid and irradiating X-rays with the subject present is referred to as an image D. In block C5, the image C is subtracted from the image D to obtain the image E with the offset corrected. In a block C6, the image E is logarithmically converted into an image F, and in a block C7, the image B is subtracted from the image F to obtain a subject image G whose gain has been corrected. In the next block C8, branching is performed depending on whether or not grid removal is to be performed in accordance with the operation (instruction) of the operator on the operation panel 20. In block C9, since the instruction is to perform grid removal, filtering processing for grid removal is performed on the image G to obtain the subject image H from which the grid has been removed. In the next block C10, the subject image H is output. If grid removal is not instructed, the subject image G is output as it is in block C11.
【００６４】また、本実施例では図１の１０のメモリ手
段を磁気ディスクなどの不揮発性の記憶媒体にし、常に
グリッドを含む画像を記憶しておき、操作者が必要に応
じて、グリッド除去画像もしくはグリッドがある画像両
者を選択もしくは同時に出力することも可能である。ま
た、画像取得系がなく、画像保存系のみの構成でも本実
施例は成立する。In this embodiment, the memory means 10 shown in FIG. 1 is replaced by a non-volatile storage medium such as a magnetic disk, and an image including a grid is always stored. Alternatively, both images having a grid can be selected or output simultaneously. Further, the present embodiment can be realized even in a configuration having no image acquisition system and only an image storage system.
【００６５】（第２の実施例）図４は第２の実施例を示
すブロック図であり、グリッド縞の周波数・Ｘ線センサ
パネルのサンプリングピッチなどの設定は図１と同様で
ある。図４では１４で示される使用目的のテーブルが用
意してある。操作者は、画像を取得したのち、１０に示
されるメモリもしくは磁気ディスクに貯えられた画像に
対し、１７で示す使用目的設定手段によって使用目的を
選択すると、自動的にグリッド除去を行うかどうかの選
択がなされる。スイッチ１３がＡ側を選択すると、記憶
装置１０の画像がフィルタリング装置１２によりフィル
タリングされて出力され、スイッチ１３がＢ側を選択す
ると、記憶装置１０の画像がフィルタリングされずにそ
のまま出力される。(Second Embodiment) FIG. 4 is a block diagram showing a second embodiment. The setting of the frequency of grid stripes, the sampling pitch of the X-ray sensor panel, and the like are the same as those in FIG. In FIG. 4, a table for use indicated by 14 is prepared. After acquiring the image, the operator selects the purpose of use of the image stored in the memory or the magnetic disk indicated by 10 by the use purpose setting means indicated by 17, and determines whether to automatically remove the grid. A choice is made. When the switch 13 selects the A side, the image of the storage device 10 is filtered and output by the filtering device 12, and when the switch 13 selects the B side, the image of the storage device 10 is output as it is without filtering.
【００６６】前述のグリッドの縞の周波数の設定で述べ
たように、１１のグリッド縞の周波数（グリッド本数）
は６〜７ｃｙｃ／ｍｍにしてあるため、観察者はグリッ
ド縞をフィルタリングで除去しなくても違和感なく画像
を観察可能である。As described in the setting of the frequency of grid stripes, the frequency of 11 grid stripes (the number of grids)
Is set to 6 to 7 cyc / mm, so that an observer can observe an image without a sense of incongruity without removing grid stripes by filtering.
【００６７】このテーブル１４は、たとえば、空間周波
数強調処理などの高い空間周波数を強調するような処理
ではグリッド像が邪魔になるので、除去する（１３のス
イッチをＡ側）。等倍以上の拡大表示もしくはハードコ
ピーであれば、画像のぼけをなるべく起こさないよう
に、スイッチ１３をＢ側に設定し、グリッド像の除去を
行わない。縮小を伴う表示又はハードコピーであれば、
スイッチ１３をＡ側に設定し、グリッド像の除去を行
う。別の記憶手段に画像を保存する目的であれば保存さ
れた画像に対して、除去処理が可能になるため、情報量
を多くし、グリッド像の除去を行わないようにするた
め、スイッチ１３をＢ側に設定する。The table 14 is removed (for example, the switch 13 is set to the A side) because the grid image becomes a hindrance in a process that emphasizes a high spatial frequency, such as a spatial frequency emphasis process. In the case of a magnified display or hard copy of the same size or larger, the switch 13 is set to the B side so that the image is not blurred as much as possible, and the grid image is not removed. If it is a display or hard copy with reduction,
The switch 13 is set to the A side to remove the grid image. If the purpose is to store the image in another storage means, the stored image can be removed, so that the switch 13 is set in order to increase the amount of information and not to remove the grid image. Set to B side.
【００６８】（第３の実施例）図５は第３の実施例を示
すブロック図であり、グリッド縞の周波数・Ｘ線センサ
パネルのサンプリングピッチなどの設定は図１と同様で
ある。図５では１８で示される撮影部位のテーブルが用
意してある。操作者は、画像撮影時に１０に示されるメ
モリもしくは磁気ディスクに貯えられた画像に対し、１
９で示す撮影部位設定手段によって使用目的を選択する
と、自動的にグリッド除去を行うかどうかの選択がなさ
れる。(Third Embodiment) FIG. 5 is a block diagram showing a third embodiment. The setting of the frequency of grid stripes, the sampling pitch of the X-ray sensor panel, and the like are the same as in FIG. In FIG. 5, a table of imaging parts indicated by 18 is prepared. At the time of capturing the image, the operator adds 1 to the image stored in the memory or the magnetic disk indicated by 10.
When the purpose of use is selected by the imaging region setting means 9, whether or not grid removal is to be performed is automatically selected.
【００６９】前述のグリッドの縞の周波数の設定で述べ
たように、１１のグリッド縞の周波数（グリッド本数）
は６〜７ｃｙｃ／ｍｍにしてあるため、観察者はグリッ
ド縞をフィルタリングで除去しなくても違和感なく画像
を観察可能である。As described above in the setting of the frequency of grid stripes, the frequency of 11 grid stripes (the number of grid lines)
Is set to 6 to 7 cyc / mm, so that an observer can observe an image without a sense of incongruity without removing grid stripes by filtering.
【００７０】このテーブル１８は、画像の観察者が、画
像に対して高い空間周波数まで必要とする骨盤や関節等
の骨部の画像であれば、スイッチ１３をＢ側に設定し、
グリッド像を除去せず、ぼけのない状態で画像を観察す
る。胸部（胸部正面）や腹部などの高い空間周波数まで
は必要とせず画像の淡い濃淡を観察しやすい状態が必要
な場合は、スイッチ１３をＡ側に設定し、グリッド像を
除去して観察することができる。If the observer of the image is an image of a bone such as a pelvis or a joint that requires a high spatial frequency for the image, the switch 13 is set to the B side.
Observe the image without blurring without removing the grid image. If you do not need high spatial frequencies such as the chest (front of the chest) or the abdomen and want to easily observe the light and shade of the image, set the switch 13 to the A side and remove the grid image for observation. Can be.
【００７１】（第４の実施例）グリッド周波数を比較的
邪魔にならないように設定したが、グリッド像のコント
ラストが強ければ、やはり観察者の邪魔になる。グリッ
ド像のコントラストは使用するＸ線の条件（エネルギー
など）によって異なるものである。さらに、グリッドが
ない状態で被写体を撮影する場合もある。(Fourth Embodiment) The grid frequency is set so as to be relatively unobtrusive. However, if the contrast of the grid image is strong, it will also obstruct the observer. The contrast of the grid image differs depending on the X-ray conditions (energy and the like) used. Further, the subject may be photographed without the grid.
【００７２】この事を判断するため、取得した画像の適
当な部分のグリッドに直交する方向のスペクトルを計算
し、グリッド成分のピーク値（対数画像であればそのま
まコントラスト値に対応する）により、グリッド除去を
行うかどうかの選択を行う。To judge this, the spectrum in the direction orthogonal to the grid of an appropriate part of the acquired image is calculated, and the peak value of the grid component (corresponding to the contrast value as it is in the case of a logarithmic image) is used. Select whether or not to remove.
【００７３】図６は、第４の実施例のブロック図であ
り、グリッド縞の周波数・Ｘ線センサパネルのサンプリ
ングピッチなどの設定は図１と同様である。前述のグリ
ッドの縞の周波数の設定で述べたように、１１のグリッ
ド縞の周波数（グリッド本数）は６〜７ｃｙｃ／ｍｍに
してあるため、観察者はグリッド縞をフィルタリングで
除去しなくても違和感なく画像を観察可能である。FIG. 6 is a block diagram of the fourth embodiment. The setting of the frequency of grid stripes, the sampling pitch of the X-ray sensor panel, and the like are the same as those in FIG. As described in the setting of the frequency of the grid stripes described above, since the frequency (the number of grids) of the 11 grid stripes is set to 6 to 7 cyc / mm, the observer can feel uncomfortable without removing the grid stripes by filtering. Without observing the image.
【００７４】図６の１５で示されるブロックはソフトウ
エアを含むものであり、内部にフローチャートを示して
いる。フローチャートのブロックＣ２１では、メモリ１
０から画像の任意のラインを取得し、ブロックＣ２２で
はこれを１次元フーリエ変換し、振幅スペクトルを計算
する。次に、ブロックＣ２３では振幅スペクトル値から
グリッド周波数に相当する振幅値（スペクトル値）Ｖｐ
を測定する。ブロックＣ２４では、振幅値Ｖｐを、あら
かじめ設定してあるしきい値ＴＨと比較する。ＶｐがＴ
Ｈより大きければ、除去しなければならないグリッド像
であり、グリッド像の除去を行うため、ブロックＣ２５
でスイッチ１３をＡ側に設定し、ＶｐがＴＨより大きく
なければ、ブロックＣ２６でスイッチ１３をＢ側に設定
する。The block indicated by 15 in FIG. 6 includes software, and internally shows a flowchart. In block C21 of the flowchart, the memory 1
An arbitrary line of the image is obtained from 0, and is subjected to a one-dimensional Fourier transform at block C22 to calculate an amplitude spectrum. Next, in block C23, the amplitude value (spectral value) Vp corresponding to the grid frequency is calculated from the amplitude spectrum value.
Is measured. In block C24, the amplitude value Vp is compared with a preset threshold value TH. Vp is T
If it is larger than H, it is a grid image that must be removed, and the block C25 is used to remove the grid image.
, The switch 13 is set to the A side, and if Vp is not greater than TH, the switch 13 is set to the B side in block C26.
【００７５】なお、グリッド周波数に存在するグリッド
縞のコントラストの大きさに応じて、スイッチ１３の選
択を決めてもよい。コントラストが所定のしきい値より
も大きいときには、スイッチ１３をＡ側に設定し、グリ
ッド像の除去を行う。一方、コントラストが所定のしき
い値よりも小さいときには、スイッチ１３をＢ側に設定
し、グリッド像の除去を行わない。The selection of the switch 13 may be determined according to the magnitude of the contrast of the grid stripe existing at the grid frequency. When the contrast is larger than the predetermined threshold, the switch 13 is set to the A side to remove the grid image. On the other hand, when the contrast is smaller than the predetermined threshold, the switch 13 is set to the B side, and the grid image is not removed.
【００７６】以上のように、本実施例によれば、Ｘ線画
像を２次元的にサンプリングし、デジタル画像を構成す
る系において、画像上の散乱線除去のためのグリッドに
起因する縞画像の空間周波数が安定した状態になるサン
プリング周波数の４０％以下でありかつ観察者の画像観
察に邪魔になり難く、画像成分に重ならない周波数（通
常ナイキスト周波数の６０％）以上になるようグリッド
本体の鉛の周波数を設定することにより、画像上からグ
リッド縞を除去しなくても観察者にとっての違和感の少
ない画像を生成できる。As described above, according to the present embodiment, in a system which samples an X-ray image two-dimensionally and forms a digital image, a stripe image caused by a grid for removing scattered radiation on the image is obtained. The lead of the grid body should be 40% or less of the sampling frequency at which the spatial frequency becomes stable, hardly obstruct the observer's image observation, and be higher than a frequency (usually 60% of the Nyquist frequency) which does not overlap with image components. By setting the frequency of, it is possible to generate an image with less discomfort for the observer without removing grid stripes from the image.
【００７７】また、フィルタリング手段１２は、スイッ
チ１３の選択に応じて、グリッド縞の除去を行うか否か
の選択が可能である。グリッド縞の除去は、操作者（観
察者）の操作、画像の使用目的、撮影部位、グリッド縞
の振幅（強度）に応じて、自動的に行われるので、必要
な場合にのみグリッド縞を除去する処理を行い、適切な
画像を得ることができる。The filtering means 12 can select whether or not to remove grid stripes according to the selection of the switch 13. The grid stripes are automatically removed according to the operation of the operator (observer), the purpose of use of the image, the imaging site, and the amplitude (intensity) of the grid stripes. Therefore, the grid stripes are removed only when necessary. And an appropriate image can be obtained.
【００７８】なお、上記実施例は、何れも本発明を実施
するにあたっての具体化のほんの一例を示したものに過
ぎず、これらによって本発明の技術的範囲が限定的に解
釈されてはならないものである。すなわち、本発明はそ
の精神、またはその主要な特徴から逸脱することなく、
様々な形で実施することができる。The above embodiments are merely examples of embodiments for carrying out the present invention, and the technical scope of the present invention should not be interpreted in a limited manner. It is. That is, the present invention does not depart from the spirit or the main features thereof,
It can be implemented in various forms.
【００７９】[0079]
【発明の効果】以上説明したように本発明によれば、画
像中から散乱線除去グリッドに起因する縞の除去を行う
か否かの選択が可能であるので、必要な場合にのみ散乱
線除去グリッドに起因する縞の除去を行い、適切な画像
を得ることができる。As described above, according to the present invention, it is possible to select whether or not to remove the fringes caused by the scattered radiation removal grid from the image, so that the scattered radiation can be removed only when necessary. Striping caused by the grid is removed, and an appropriate image can be obtained.
【図１】本発明の第１の実施例を説明するためのブロッ
ク図である。FIG. 1 is a block diagram for explaining a first embodiment of the present invention.
【図２】画像中の空間スペクトルを説明するための図で
ある。FIG. 2 is a diagram for explaining a spatial spectrum in an image.
【図３】第１の実施例をソフトウエアで実行するための
フローチャートである。FIG. 3 is a flowchart for executing the first embodiment by software.
【図４】本発明の第２の実施例を説明するためのブロッ
ク図である。FIG. 4 is a block diagram for explaining a second embodiment of the present invention.
【図５】本発明の第３の実施例を説明するためのブロッ
ク図である。FIG. 5 is a block diagram for explaining a third embodiment of the present invention.
【図６】本発明の第４の実施例を説明するためのブロッ
ク図である。FIG. 6 is a block diagram for explaining a fourth embodiment of the present invention.
【図７】グリッド縞に起因するアーチファクトを説明す
るための図である。FIG. 7 is a diagram for explaining an artifact caused by grid fringes.
【図８】グリッド縞をサンプリングした時のビートにつ
いて説明するための図である。FIG. 8 is a diagram for explaining beats when grid stripes are sampled.
【図９】グリッド本体の空間周波数の範囲の例を示す図
である。FIG. 9 is a diagram showing an example of a spatial frequency range of a grid body.
【図１０】散乱線除去グリッドを説明するための図であ
る。FIG. 10 is a diagram for explaining a scattered radiation removal grid.
１ Ｘ線管球 ２ 被写体 ３ Ｘ線画像取得用のセンサ ４ コントローラ ５ アナログ／デジタル変換器 ６ メモリ（記憶装置） ７ ルックアップテーブル（対数変換装置） ８ 減算器 ９ 記憶装置 １０ 記憶装置 １１ 散乱線除去グリッド １２ フィルタリング手段（グリッド縞除去手段） １３ スイッチ １４ テーブル １５ 処理ブロック １６ グリッド除去操作実効設定手段 １７ 使用目的設定手段 １８ テーブル １９ 撮影部位設定手段 ２０ 操作パネル ８１ Ｘ線焦点（Ｘ線管球） ８２ 被写体 ８３ グリッド ８４ エネルギー変換部 ８５ センサ DESCRIPTION OF SYMBOLS 1 X-ray tube 2 Subject 3 X-ray image acquisition sensor 4 Controller 5 Analog / digital converter 6 Memory (storage device) 7 Look-up table (logarithmic conversion device) 8 Subtractor 9 Storage device 10 Storage device 11 Scattered radiation Removal grid 12 Filtering means (grid stripe removal means) 13 Switch 14 Table 15 Processing block 16 Grid removal operation effective setting means 17 Purpose of use setting means 18 Table 19 Imaging part setting means 20 Operation panel 81 X-ray focal point (X-ray tube) 82 subject 83 grid 84 energy conversion unit 85 sensor
Claims (21)
射し、被写体のＸ線透過分布を２次元的に所望の間隔で
サンプリングして画像を取得する画像取得手段と、 前記取得した画像中から前記散乱線除去グリッドに起因
する縞を画像処理により除去する画像処理手段と、 前記画像処理手段により前記散乱線除去グリッドに起因
する縞を除去するか否かの選択を可能にする選択手段と
を有することを特徴とする画像取得装置。1. An image acquisition means for irradiating X-rays using a scattered radiation removal grid and sampling an X-ray transmission distribution of a subject two-dimensionally at a desired interval to acquire an image. An image processing unit that removes stripes caused by the scattered radiation removal grid from the image by image processing; and a selection that enables the image processing unit to select whether to remove the stripes caused by the scattered radiation removal grid. And an image acquisition device.
理により前記散乱線除去グリッドに起因する縞を除去す
る請求項１記載の画像取得装置。2. The image acquisition apparatus according to claim 1, wherein the image processing unit removes a stripe caused by the scattered radiation removal grid by a filtering process.
前記縞を除去するか否かの選択を行う請求項１又は２記
載の画像取得装置。3. The image acquisition apparatus according to claim 1, wherein the selection unit selects whether or not to remove the stripe according to an operation of an operator.
て前記縞を除去するか否かの選択を行う請求項１又は２
記載の画像取得装置。4. The image processing apparatus according to claim 1, wherein the selection unit selects whether or not to remove the stripe according to a purpose of use of the image.
The image acquisition device according to the above.
処理を目的とするときには前記縞を除去する選択を行う
請求項４記載の画像取得装置。5. The image acquisition apparatus according to claim 4, wherein the selection unit performs a selection for removing the stripe when performing spatial frequency enhancement processing on the image.
表示又はハードコピーを目的とするときには前記縞を除
去しない選択を行う請求項４記載の画像取得装置。6. The image acquisition apparatus according to claim 4, wherein the selection unit performs a selection not to remove the stripes when the purpose is an enlarged display or a hard copy of an image that is equal to or larger than the same size.
又はハードコピーを目的とするときには前記縞を除去す
る選択を行う請求項４記載の画像取得装置。7. The image acquisition apparatus according to claim 4, wherein the selection unit performs a selection for removing the stripes when a display or a hard copy accompanied by reduction of an image is intended.
るときには前記縞を除去しない選択を行う請求項４記載
の画像取得装置。8. The image acquisition apparatus according to claim 4, wherein the selection unit performs selection of not removing the stripe when saving an image.
縞を除去するか否かの選択を行う請求項１又は２記載の
画像取得装置。9. The image acquisition apparatus according to claim 1, wherein the selection unit selects whether or not to remove the stripe according to an imaging region.
るときには前記縞を除去しない選択を行う請求項９記載
の画像取得装置。10. The image acquisition apparatus according to claim 9, wherein the selection unit performs a selection not to remove the stripe when the imaging part is a bone part.
関節であるときには前記縞を除去しない選択を行う請求
項９記載の画像取得装置。11. The image acquisition apparatus according to claim 9, wherein the selection unit performs a selection not to remove the stripe when the imaging part is a pelvis or a joint.
腹部であるときには前記縞を除去する選択を行う請求項
９記載の画像取得装置。12. The image acquisition apparatus according to claim 9, wherein the selection unit performs selection for removing the stripe when the imaging part is a chest or abdomen.
ッドに起因する縞の周波数の振幅の大きさに応じて前記
縞を除去するか否かの選択を行う請求項１又は２記載の
画像取得装置。13. The image acquisition apparatus according to claim 1, wherein the selection unit selects whether or not to remove the fringes according to the magnitude of the frequency amplitude of the fringes caused by the scattered radiation removal grid. apparatus.
ッドに起因する縞の周波数の振幅が所定のしきい値より
も大きいときには前記縞を除去する選択を行う請求項１
３記載の画像取得装置。14. The apparatus according to claim 1, wherein said selection means selects said stripe to be removed when the amplitude of the frequency of the stripe caused by said scattered radiation removal grid is larger than a predetermined threshold value.
3. The image acquisition device according to 3.
ッドに起因する縞のコントラストの大きさに応じて前記
縞を除去するか否かの選択を行う請求項１又は２記載の
画像取得装置。15. The image acquisition apparatus according to claim 1, wherein the selection unit selects whether or not to remove the stripe according to the magnitude of the contrast of the stripe caused by the scattered radiation removal grid.
画像からオフセット値を除去するオフセット値除去手段
を有する請求項１〜１５のいずれかに記載の画像取得装
置。16. The image acquisition apparatus according to claim 1, further comprising an offset value removing unit that removes an offset value from an image acquired by exposing the X-ray.
線を曝射して取得した画像からＸ線を曝射せずに取得し
た画像を減算する第１の減算手段を含む請求項１６記載
の画像取得装置。17. The method according to claim 17, wherein the offset value removing unit is configured to:
17. The image acquisition device according to claim 16, further comprising a first subtraction unit that subtracts an image acquired without exposing X-rays from an image acquired by exposing the lines.
イン値のばらつきを補正するゲイン値補正手段を有する
請求項１〜１５のいずれかに記載の画像取得装置。18. The image acquisition apparatus according to claim 1, further comprising a gain value correction unit configured to correct a variation in a gain value with respect to the acquired image.
して、前記ゲイン値補正手段は、前記取得した画像に対
して被写体がない状態で取得した画像の割り算を行う割
り算手段を有する請求項１８記載の画像取得装置。19. A method according to claim 18, wherein said gain value correcting means has a dividing means for dividing the obtained image by an image obtained without a subject. The image acquisition device according to the above.
た画像を対数変換した画像から、被写体がない状態で取
得した画像を対数変換した画像を減算する第２の減算手
段を含む請求項１８記載の画像取得装置。20. The gain value correcting means includes a second subtracting means for subtracting an image obtained by logarithmically converting an image obtained without a subject from an image obtained by logarithmically converting the obtained image. Image acquisition device.
Ｘ線を曝射し、被写体のＸ線透過分布を２次元的に所望
の間隔でサンプリングして画像を取得するステップと、 （ｂ）前記取得した画像中から前記散乱線除去グリッド
に起因する縞を除去するか否かを選択するステップと、 （ｃ）前記選択に応じて前記散乱線除去グリッドに起因
する縞を画像処理により除去するステップとを有するこ
とを特徴とする画像取得方法。21. (a) irradiating X-rays using a scattered radiation removal grid and sampling an X-ray transmission distribution of a subject two-dimensionally at a desired interval to obtain an image; And (c) selecting whether or not to remove stripes caused by the scattered radiation removal grid from the acquired image; and (c) removing the stripes caused by the scattered radiation removal grid by image processing according to the selection. Performing an image acquisition method.
Priority Applications (5)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2000028207A JP4612754B2 (en) | 2000-02-04 | 2000-02-04 | Image acquisition apparatus and image acquisition method |
| US09/773,095 US20010033638A1 (en) | 2000-02-04 | 2001-01-31 | Image acquisition method and apparatus |
| EP20010102319 EP1126289B1 (en) | 2000-02-04 | 2001-02-01 | Image acquisition method and apparatus |
| US10/210,739 US6826256B2 (en) | 2000-02-04 | 2002-08-01 | Apparatus and method for a radiation image through a grid |
| US10/719,331 US6868137B2 (en) | 2000-02-04 | 2003-11-20 | Apparatus and method for a radiation image through a grid |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2000028207A JP4612754B2 (en) | 2000-02-04 | 2000-02-04 | Image acquisition apparatus and image acquisition method |
Publications (3)
| Publication Number | Publication Date |
|---|---|
| JP2001212119A true JP2001212119A (en) | 2001-08-07 |
| JP2001212119A5 JP2001212119A5 (en) | 2007-03-22 |
| JP4612754B2 JP4612754B2 (en) | 2011-01-12 |
Family
ID=18553658
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2000028207A Expired - Fee Related JP4612754B2 (en) | 2000-02-04 | 2000-02-04 | Image acquisition apparatus and image acquisition method |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP4612754B2 (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2004329783A (en) * | 2003-05-12 | 2004-11-25 | Canon Inc | Radiation imaging apparatus and method |
| JP2005274576A (en) * | 2004-03-25 | 2005-10-06 | Mitsutoyo Corp | Image inspection method and image inspection system |
| US7433537B2 (en) | 2001-06-19 | 2008-10-07 | Fujifilm Corporation | Image reproduction system for removing noise of a periodic pattern |
| JP2012200455A (en) * | 2011-03-25 | 2012-10-22 | Fujifilm Corp | Radiographic imaging system and program |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2507659B2 (en) * | 1990-04-04 | 1996-06-12 | 株式会社東芝 | Radiation image reader |
| JP2754068B2 (en) * | 1989-06-09 | 1998-05-20 | 富士写真フイルム株式会社 | Radiation image signal generation method and radiation image reading device |
| JP2981701B2 (en) * | 1992-09-08 | 1999-11-22 | 富士写真フイルム株式会社 | Radiation imaging method and bone mineral quantitative analysis method using the same |
| JP2000023954A (en) * | 1998-06-02 | 2000-01-25 | General Electric Co <Ge> | Method for removing grid line artifact from x-ray image |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS61133847A (en) * | 1984-12-04 | 1986-06-21 | Shoichi Yasui | Promation of image clear of scattered line in x ray image |
| JPH02250186A (en) * | 1988-08-02 | 1990-10-05 | Toshiba Corp | Method and device for reconstituting picture |
| JPH0486157A (en) * | 1990-07-30 | 1992-03-18 | Hitachi Medical Corp | Picture digitizing and compressing device |
| JPH09160149A (en) * | 1995-12-12 | 1997-06-20 | Fuji Photo Film Co Ltd | Image processing method and device |
| JPH1098587A (en) * | 1996-09-20 | 1998-04-14 | Fuji Photo Film Co Ltd | Radiation picture reading device |
| JPH11113889A (en) * | 1997-10-16 | 1999-04-27 | Canon Inc | Radiography unit, radiography device, and radiography system |
-
2000
- 2000-02-04 JP JP2000028207A patent/JP4612754B2/en not_active Expired - Fee Related
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2754068B2 (en) * | 1989-06-09 | 1998-05-20 | 富士写真フイルム株式会社 | Radiation image signal generation method and radiation image reading device |
| JP2507659B2 (en) * | 1990-04-04 | 1996-06-12 | 株式会社東芝 | Radiation image reader |
| JP2981701B2 (en) * | 1992-09-08 | 1999-11-22 | 富士写真フイルム株式会社 | Radiation imaging method and bone mineral quantitative analysis method using the same |
| JP2000023954A (en) * | 1998-06-02 | 2000-01-25 | General Electric Co <Ge> | Method for removing grid line artifact from x-ray image |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7433537B2 (en) | 2001-06-19 | 2008-10-07 | Fujifilm Corporation | Image reproduction system for removing noise of a periodic pattern |
| JP2004329783A (en) * | 2003-05-12 | 2004-11-25 | Canon Inc | Radiation imaging apparatus and method |
| JP4497837B2 (en) * | 2003-05-12 | 2010-07-07 | キヤノン株式会社 | Radiation imaging equipment |
| JP2005274576A (en) * | 2004-03-25 | 2005-10-06 | Mitsutoyo Corp | Image inspection method and image inspection system |
| JP4641441B2 (en) * | 2004-03-25 | 2011-03-02 | 株式会社ミツトヨ | Image inspection method and image inspection system |
| JP2012200455A (en) * | 2011-03-25 | 2012-10-22 | Fujifilm Corp | Radiographic imaging system and program |
| US8735829B2 (en) | 2011-03-25 | 2014-05-27 | Fujifilm Corporation | Radiographic image capturing system, program storage medium, and method |
Also Published As
| Publication number | Publication date |
|---|---|
| JP4612754B2 (en) | 2011-01-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP3903027B2 (en) | Radiation image processing method and apparatus, and grid selection method and apparatus | |
| JP3193806B2 (en) | Noise reduction method and apparatus | |
| US5878108A (en) | Method for generating X-ray image and apparatus therefor | |
| JP4393483B2 (en) | Radiation image processing apparatus, image processing system, radiation image processing method, storage medium, and program | |
| US6868137B2 (en) | Apparatus and method for a radiation image through a grid | |
| KR20130128124A (en) | Panorama image data providing method and apparatus | |
| JP2000023954A (en) | Method for removing grid line artifact from x-ray image | |
| JP3557567B2 (en) | X-ray CT system | |
| EP1126289B1 (en) | Image acquisition method and apparatus | |
| JP4500400B2 (en) | Image acquisition apparatus and image acquisition method | |
| Wang et al. | Film digitization aliasing artifacts caused by grid line patterns | |
| JP3793053B2 (en) | Radiation image processing apparatus, image processing system, radiation image processing method, storage medium, and program | |
| JP4137499B2 (en) | Phase information restoration method, phase information restoration device, and phase information restoration program | |
| JP4612754B2 (en) | Image acquisition apparatus and image acquisition method | |
| JP3223195B2 (en) | X-ray CT system | |
| JP2005040506A (en) | Method and apparatus for processing x-ray image | |
| JP4746761B2 (en) | Radiation image processing apparatus, radiation image processing method, storage medium, and program | |
| JP3755144B2 (en) | Radiography equipment | |
| JP3783116B2 (en) | Radiation image enhancement processing method and apparatus | |
| JP3825989B2 (en) | Radiation image processing apparatus, image processing system, radiation image processing method, computer-readable storage medium, and computer program | |
| JP3445258B2 (en) | Radiation image processing apparatus, image processing system, radiation image processing method, storage medium, program, radiation imaging apparatus, and radiation imaging system | |
| JP2002325755A (en) | Device, system and method for processing radiograph, computer readable recording medium recording program and program | |
| JP4677339B2 (en) | Radiographic image acquisition apparatus, radiographic image acquisition method and design method | |
| JP3884929B2 (en) | Radiation image acquisition apparatus and design method | |
| JP4393436B2 (en) | Radiation image processing apparatus, image processing system, radiation image processing method, storage medium, and program |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20070201 |
|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20070201 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20091208 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20091215 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20100209 |
|
| A02 | Decision of refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A02 Effective date: 20100413 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20100430 |
|
| A911 | Transfer to examiner for re-examination before appeal (zenchi) |
Free format text: JAPANESE INTERMEDIATE CODE: A911 Effective date: 20100713 |
|
| 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: 20101005 |
|
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20101016 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20131022 Year of fee payment: 3 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 4612754 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| LAPS | Cancellation because of no payment of annual fees |