JPH10262959A - Method for correcting and preserving digital data of x-ray radiographic images - Google Patents

Method for correcting and preserving digital data of x-ray radiographic images

Info

Publication number
JPH10262959A
JPH10262959A JP9075745A JP7574597A JPH10262959A JP H10262959 A JPH10262959 A JP H10262959A JP 9075745 A JP9075745 A JP 9075745A JP 7574597 A JP7574597 A JP 7574597A JP H10262959 A JPH10262959 A JP H10262959A
Authority
JP
Japan
Prior art keywords
digital data
ray
image
correction
test chart
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
Application number
JP9075745A
Other languages
Japanese (ja)
Other versions
JP3863963B2 (en
Inventor
Masanao Watanabe
正直 渡辺
Kiyoshi Matsuo
清 松尾
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Dai Nippon Printing Co Ltd
Original Assignee
Dai Nippon Printing Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Dai Nippon Printing Co Ltd filed Critical Dai Nippon Printing Co Ltd
Priority to JP07574597A priority Critical patent/JP3863963B2/en
Publication of JPH10262959A publication Critical patent/JPH10262959A/en
Application granted granted Critical
Publication of JP3863963B2 publication Critical patent/JP3863963B2/en
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Landscapes

  • Measuring And Recording Apparatus For Diagnosis (AREA)
  • Radiography Using Non-Light Waves (AREA)
  • Apparatus For Radiation Diagnosis (AREA)

Abstract

PROBLEM TO BE SOLVED: To unify the contrast and size of digital radiographic image among plural images. SOLUTION: Inside the area to be irradiated with X-ray beams 26 generated from an X-ray generator 24, a test chart 10B showing different X-ray absorption coefficients in plural divided areas is arranged at a certain position so as not to be overlapped with a patient 22, and radiographed together with the patient 22 and the density of chart image in the divided area corresponding to the test chart 10B in the digital image provided by an X-ray film 30 is read. Then, the image is corrected so that this density can be coincident with a reference standard value, and the same diseased part can be displayed with the same density in the plural images.

Description

【発明の詳細な説明】DETAILED DESCRIPTION OF THE INVENTION

【0001】[0001]

【発明の属する技術分野】本発明は、例えば医療の分野
において、X線撮影された画像中に、被写体あるいはそ
の近傍領域でテストチャートにより指標を残すと共に、
該指標に基づいて、撮影画像のデジタルデータを濃度、
寸法補正して撮影条件の異なる画像を統一して保存でき
るようにしたX線撮影画像の補正、保存方法に関するも
のである。
BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to, for example, in the field of medical treatment, in an X-ray photographed image, an index is left by a test chart in a subject or a region near the subject,
Based on the index, the digital data of the captured image is
The present invention relates to a method for correcting and storing an X-ray photographed image in which images having different photographing conditions can be unified and stored by performing dimension correction.

【0002】[0002]

【従来の技術】従来の、この種のX線撮影評価・補正用
テストチャートとしては、例えば特公昭62−5377
0号公報に開示されるように、X線透過性の良好な材質
からなる基板に、形状がC字型又はO字型をなす金属箔
性の複数個のテストピースを配設し、このテストピース
の外形寸法を所定方向に順次段階的に異ならせ、更に、
この異なったそれぞれの外形寸法のものに対応して、そ
の厚さを順次段階的に異ならせたものを前記所定方向と
は別方向に、前記基板上に配設してなるX線写真評価用
テストチャートがある。前記テストピースの厚さの調整
は、重ね合わせやエッチングで行い、これによりX線吸
収係数を変化させている。
2. Description of the Related Art A conventional test chart for this type of X-ray imaging evaluation / correction is disclosed, for example, in Japanese Patent Publication No. Sho 62-5377.
As disclosed in Japanese Unexamined Patent Publication No. 0, a plurality of test pieces made of a metal foil having a C-shape or an O-shape are arranged on a substrate made of a material having good X-ray transparency. The outer dimensions of the pieces are made to differ gradually in a predetermined direction, and further,
An X-ray photographic evaluation device is provided in which the thicknesses of the different external dimensions are sequentially and stepwise varied on the substrate in a direction different from the predetermined direction. There is a test chart. The thickness of the test piece is adjusted by superposition or etching, thereby changing the X-ray absorption coefficient.

【0003】又、例えば特開昭61−248665号公
報や、特開平5−329141号公報に開示されるもの
がある。
Further, there are, for example, those disclosed in JP-A-61-248665 and JP-A-5-329141.

【0004】前者は、被写体をX線撮影する際に、写し
込み装置によって、撮影年月日、被写体の氏名等に加え
て、標準濃度チャートを写し込み、画像読取り後、その
チャートの情報に基づいてフィルムの特性曲線を算出す
るようにしたものである。
[0004] In the former, when a subject is X-ray photographed, a standard density chart is photographed by an imprinting device in addition to a photographing date, a subject name, and the like. Thus, the characteristic curve of the film is calculated.

【0005】又後者は、被写体とX線検出器との間に標
準ファントムを配置し、X線検出器を走査して測定した
X線透過画像を用いて計算した標準ファントムの測定値
を用い、正しい計算ができるように補正を行うX線診断
装置及び補正方法である。
In the latter, a standard phantom is arranged between a subject and an X-ray detector, and a measured value of the standard phantom calculated using an X-ray transmission image measured by scanning the X-ray detector is used. An X-ray diagnostic apparatus and a correction method for performing correction so that correct calculation can be performed.

【0006】[0006]

【発明が解決しようとする課題】ここで、医療分野にお
いては、X線撮影は、同一の患者について複数個所角度
を変えたりしながら撮影することが多い。更には、同一
の患者について、時間の経過に伴って複数回撮影するこ
とも多いが、その都度、X線管の印加電圧、電流、フィ
ルムもメーカーによる特性の相違、更には同一メーカー
のX線フィルムでも、製造ロッドによって特性が異なる
ことが多く、撮影された患者が同一人であっても、X線
撮影画像は異なる濃度となってしまうことがある。
Here, in the medical field, X-ray photography is often performed while changing the angle of a same patient at a plurality of positions. Furthermore, the same patient is often photographed a plurality of times with the lapse of time. Each time, the applied voltage, current, and film of the X-ray tube have different characteristics depending on the manufacturer, and further, X-rays of the same manufacturer. Films often have different characteristics depending on the manufacturing rod, and even if the imaged patient is the same person, the X-ray images may have different densities.

【0007】更には、X線管、患者、フィルムの三者の
相対的距離を常に一定とすることは困難であるので、同
一人の同一患部を撮影しても、その大きさが異なること
もある。
Further, since it is difficult to keep the relative distance between the X-ray tube, the patient, and the film constant at all times, even if the same affected part is photographed by the same person, the size may differ. is there.

【0008】これに対して、得られたX線撮影画像をコ
ンピュータ処理して、その濃度、大きさを補正すること
が考えられるが、基準が不明確であるために補正が恣意
的に行われ易く、濃度の評価や相対的な比較に正確さを
欠くという問題点がある。
On the other hand, it is conceivable to correct the density and size of the obtained radiographic image by computer processing, but the correction is arbitrarily performed because the standard is not clear. There is a problem that the method is easy and the evaluation of the concentration and the relative comparison lack accuracy.

【0009】又、、撮影したX線画像の分解能(解像
力)についても、同一患者の同一患部についての複数の
X線撮影画像において、統一しない場合は、画像から患
部組織の状態を正確に把握することができないという問
題点がある。
In addition, if the resolution (resolution) of the taken X-ray image is not unified in a plurality of X-ray images of the same affected part of the same patient, the state of the affected part tissue is accurately grasped from the image. There is a problem that it is not possible.

【0010】更に、上記従来のようなテストチャート、
フィルム画像読取装置あるいはX線診断装置及び補正方
法では、X線撮影下に、被写体の特定した位置を、患部
近傍に指標として残すことが困難であり、又残したとし
ても、患部組織の異なる多岐のX線吸収特性に対応し
て、X線撮影画像を評価・補正し、これによってX線撮
影画像から患部の位置と大きさを特定するということが
困難であるという問題点があった。
[0010] Further, the conventional test chart as described above,
With a film image reading apparatus, an X-ray diagnostic apparatus, and a correction method, it is difficult to leave a specified position of a subject as an index near an affected part under X-ray imaging. There is a problem that it is difficult to evaluate and correct an X-ray radiographic image corresponding to the X-ray absorption characteristics described above, and thereby to specify the position and size of the affected part from the X-ray radiographic image.

【0011】この発明は、上記従来の問題点を解消すべ
くなされたものであって、X線撮影下に、被写体の特定
した位置を、大きさ、及びX線吸収係数が規定されたテ
ストチャートにより患部近傍に指標として残し、且つこ
の指標が、患部組織の異なる多岐のX線吸収特性に対応
することができるようにできるようにし、これにより、
複数のX線撮影画像における濃淡の程度(コントラス
ト)、表面の粗密状態(粒状性)、撮影方向、撮影個所
の大きさ、厚さ、検出エレメントの感度特性、分解能
(解像力)の補正・評価及び保存をする方法を提供する
ことを目的とする。
SUMMARY OF THE INVENTION The present invention has been made to solve the above-mentioned conventional problems, and is a test chart in which a specified position of a subject, a size, and an X-ray absorption coefficient are defined under X-ray photography. To leave an index near the affected area, and this index can correspond to a variety of different X-ray absorption characteristics of the affected tissue, so that
Correction and evaluation of the degree of contrast (contrast), the density of the surface (granularity), the imaging direction, the size and thickness of the imaging location, the sensitivity characteristics of the detection element, and the resolution (resolution) of a plurality of X-ray images, and The purpose is to provide a way to save.

【0012】[0012]

【課題を解決するための手段】本発明は、請求項1記載
のように、X線を吸収しない、又は、X線吸収係数が低
いシート上に、X線吸収係数が高い金属を含み、そのX
線吸収係数が所定規格となる膜厚で、且つ、大きさ及び
撮影方向を特定できるパターンのX線吸収膜を形成し、
このX線吸収膜を被って、被写体へのシート貼付用の粘
着層を形成してなるX線撮影評価・補正用テストチャー
トを、前記粘着層により被写体及びX線照射領域内であ
って、被写体にかぶらない位置に貼着して、これらのテ
ストチャートを被写体と共にX線撮影する過程と、撮影
した画像をデジタルデータとして記憶装置に取り込む過
程と、前記記憶装置に取り込まれたデジタルデータを、
ディスプレイ上に表示する過程と、ディスプレイ上にお
ける前記被写体にかぶっていないテストチャートのチャ
ート画像の濃度を読み取る過程と、読み取った濃度のデ
ジタルデータを標準値に合わせるように補正する過程
と、この補正値に基づき画像全体のデジタルデータを補
正する過程と、補正したデジタルデータを前記記憶装置
に保存する過程と、を有してなるX線撮影画像のデジタ
ルデータ補正、保存方法により、上記目的を達成するも
のである。
According to the present invention, a sheet which does not absorb X-rays or has a high X-ray absorption coefficient is provided on a sheet which does not absorb X-rays. X
Forming an X-ray absorption film having a film thickness whose line absorption coefficient becomes a predetermined standard, and a pattern capable of specifying the size and the imaging direction,
An X-ray imaging evaluation / correction test chart formed by forming an adhesive layer for sticking a sheet to a subject by covering the X-ray absorbing film is used to obtain a test chart for the subject and the X-ray irradiation area by the adhesive layer. Affixing the test chart together with the subject at a position where it does not cover the subject, X-ray imaging of these test charts, capturing the captured image as digital data in a storage device, and digital data captured in the storage device,
Displaying on the display, reading the density of the chart image of the test chart not covering the subject on the display, correcting the digital data of the read density to a standard value, and correcting the correction value. The above object is achieved by a method of correcting and storing digital data of an X-ray radiographic image, comprising the steps of: correcting digital data of the entire image based on the image data; and storing the corrected digital data in the storage device. Things.

【0013】請求項2の発明は、請求項1の発明におい
て、前記チャート画像の濃度の読み取りを白、黒、グレ
ーの3色の色値を読み取ることにより行い、これらのう
ち少なくとも白及び黒の色値の目標値を予め設定してお
き、前記ディスプレイ上の画像から読み取った色値のデ
ジタルデータと前記目標値とから補正値を算出するもの
である。
According to a second aspect of the present invention, in the first aspect of the present invention, the reading of the density of the chart image is performed by reading three color values of white, black, and gray. A target value of a color value is set in advance, and a correction value is calculated from digital data of a color value read from an image on the display and the target value.

【0014】請求項3の発明は、請求項1の発明におい
て、ディスプレイ上における前記被写体に貼着したテス
トチャートのチャート画像の大きさのデジタルデータに
基づくアウトプットの解像度と、該テストチャートの、
前記と同一個所の寸法を前記記憶装置に登録したときの
デジタルデータに基づくアウトプットの解像度とを一致
させて、該テストチャート近傍における被写体の寸法の
デジタルデータを補正する過程と、補正したデジタルデ
ータを前記記憶装置に保存する過程と、を含んで構成し
たものである。
According to a third aspect of the present invention, in the first aspect of the present invention, an output resolution based on digital data of a size of a chart image of a test chart attached to the subject on a display, and an output resolution of the test chart,
Correcting the digital data of the dimensions of the subject in the vicinity of the test chart by matching the resolution of the output based on the digital data when the dimensions of the same place are registered in the storage device; and In the storage device.

【0015】請求項4の発明は、請求項1、2又は3の
発明において、前記被写体を前記テストチャートと共に
異なる角度で複数回X線撮影する過程と、各X線撮影に
よる画像をデジタルデータとして記憶装置に取り込む過
程と、これら異なる撮影角度の複数の画像のデジタルデ
ータを、3次元グラフィックソフトで数式処理して、前
記被写体の立体画像を得て、そのデジタルデータを前記
記憶装置に記憶させる過程と、を含んで構成したもので
ある。
According to a fourth aspect of the present invention, in the first, second or third aspect of the present invention, the step of X-ray photographing the subject together with the test chart at different angles a plurality of times and an image obtained by each X-ray photographing as digital data. A step of capturing the digital data of a plurality of images at different photographing angles with a three-dimensional graphic software to obtain a three-dimensional image of the subject, and storing the digital data in the storage device. And is included.

【0016】請求項5の発明は、請求項1乃至4のいず
れかの発明において、前記テストチャートにおけるX線
吸収膜をリング状に形成し、該リングは円周方向に複数
の領域に分割され、且つ、各領域毎に異なる膜厚とした
ものである。
According to a fifth aspect of the present invention, in any one of the first to fourth aspects, the X-ray absorbing film in the test chart is formed in a ring shape, and the ring is divided into a plurality of regions in a circumferential direction. In addition, the thickness is different for each region.

【0017】請求項6の発明は、請求項1乃至5のいず
れかの発明において、前記補正値を、前記補正したデジ
タルデータと独立して保存する過程を設けてなり、次回
以降の同一撮影条件の撮像時に、前記テストチャートの
被写体への貼付過程、及び、テストチャートのチャート
画像の濃度読み取り過程の省略可能としたものである。
According to a sixth aspect of the present invention, in any one of the first to fifth aspects of the present invention, a step of storing the correction value independently of the corrected digital data is provided. At the time of imaging, the process of attaching the test chart to the subject and the process of reading the density of the chart image of the test chart can be omitted.

【0018】[0018]

【発明の実施の形態】以下本発明の実施の形態の例を図
面を参照して詳細に説明する。
Embodiments of the present invention will be described below in detail with reference to the drawings.

【0019】本発明方法に用いるX線撮影評価・補正用
のテストチャート10は、図1及び図2に示されるよう
に、X線を吸収しない、又はX線吸収係数が低い、例え
ば透明易接着PETフィルム12を基板用シートとし
て、その処理面にX線吸収係数が高い金属を含み、その
X線吸収係数が所定規格となる膜厚で、且つ大きさ及び
撮影方向が特定されるパターンのX線吸収膜14を形成
し、更にその外側を被写体(例えば人体)の表皮に粘着
できるようにするための粘着層16を形成し、更にこの
粘着層16を離型紙18によって被ったものである。
As shown in FIGS. 1 and 2, a test chart 10 for evaluating and correcting X-ray photography used in the method of the present invention does not absorb X-rays or has a low X-ray absorption coefficient. The PET film 12 is used as a substrate sheet. The processed surface contains a metal having a high X-ray absorption coefficient. A line absorbing film 14 is formed, and an adhesive layer 16 for allowing the outside to adhere to the skin of a subject (for example, a human body) is formed. The adhesive layer 16 is further covered with a release paper 18.

【0020】前記X線吸収膜14は、その平面形状がリ
ング状とされ、中心角が90°の4つの分割領域20A
〜20Dに区画されている。
The X-ray absorbing film 14 has a ring shape in plan view and four divided regions 20A having a central angle of 90 °.
2020D.

【0021】これら分割領域20A〜20Dは、この順
で膜厚が段階的に厚くなるようにされて、これにより、
X線吸収係数が段階的に4階調となるようにされてい
る。
The thickness of each of the divided regions 20A to 20D is increased stepwise in this order.
The X-ray absorption coefficient is set to be four gradations step by step.

【0022】又、リング状のX線吸収膜14の内側に
は、粘着層16及びPETフィルム12からなる小円形
の透過領域20Eが形成され、これによって、テストチ
ャート10は、全体として5階調のX線吸収係数の領域
から構成されている。
Further, inside the ring-shaped X-ray absorbing film 14, a small circular transmission region 20E composed of the adhesive layer 16 and the PET film 12 is formed, whereby the test chart 10 has five gradations as a whole. Of the X-ray absorption coefficient.

【0023】前記X線吸収膜14に含まれるX線吸収係
数が高い金属は、鉛、ビスマス、バリウム、これらの化
合物あるいはこれらの混合物であって、いずれも粉末状
のものを、バインダーに混ぜてX線遮断インキとして、
PETフィルム12上にパターン印刷する。
The metal having a high X-ray absorption coefficient contained in the X-ray absorption film 14 is lead, bismuth, barium, a compound thereof, or a mixture thereof. As X-ray blocking ink,
A pattern is printed on the PET film 12.

【0024】あるいは、前記金属粉末を、光硬化性樹脂
からなるバインダーに混合してX線遮断インキとし、パ
ターン印刷後に硬化用光源から光を照射して、該パター
ンを硬化させることにより形成する。
Alternatively, the metal powder is mixed with a binder made of a photocurable resin to form an X-ray blocking ink, and after the pattern is printed, light is irradiated from a curing light source to cure the pattern.

【0025】前記バインダーとしては、アクリル系、ウ
レタン系、塩化酢酸ビニール系、ポリエステル系のう
ち、少なくとも1つの樹脂から構成する。又、バインダ
ーとしての光硬化性樹脂は、例えば紫外線硬化されるウ
レタンアクリレート、ポリエステルアクリレートあるい
はこれらの混合物とする。
The binder is made of at least one resin selected from the group consisting of acrylic, urethane, vinyl chloride, and polyester. The photo-curable resin as the binder is, for example, urethane acrylate, polyester acrylate or a mixture thereof which is cured by ultraviolet rays.

【0026】上記のような金属粉末、バインダー、溶
剤、P/V比は、例えば次の表1のようにする。
The above-mentioned metal powder, binder, solvent, and P / V ratio are, for example, as shown in Table 1 below.

【0027】[0027]

【表1】 [Table 1]

【0028】次に、上記テストチャート10の製造方法
について説明する。
Next, a method of manufacturing the test chart 10 will be described.

【0029】まず、スクリーン印刷用インキSSから表
1に示されるような、X線遮断インキ〜を、3本ロ
ールで混練形成する。
First, X-ray blocking inks as shown in Table 1 are kneaded with three rolls from the screen printing ink SS.

【0030】次に、図3に示されるパターンP1で、レ
ジスト厚100μmに構成したシルクスクリーン版(テ
トロン150メッシュ、株式会社ミノグループ製)で、
厚さが100μmの透明易接着PETフィルム(東レ株
式会社製)の処理面に、110μmの塗膜厚に印刷す
る。図3の符号h1は、空白(孔)となる部分を示す。
Next, a silk screen plate (Tetron 150 mesh, manufactured by Mino Group Co., Ltd.) having a pattern P1 shown in FIG.
On a treated surface of a 100 μm-thick transparent easy-adhesion PET film (manufactured by Toray Industries, Inc.), a coating film thickness of 110 μm is printed. The reference numeral h1 in FIG. 3 indicates a blank portion (hole).

【0031】膜厚及び塗膜状態は、レーザフォーカス変
位計(株式会社キーエンス製)と、デジマチックインジ
ケータ(ミツトヨ株式会社製)で確認の上、オーブンで
100℃、1時間乾燥し、90〜100μの膜厚のパタ
ーンを得る。
The film thickness and the state of the coating film were confirmed with a laser focus displacement meter (manufactured by Keyence Corporation) and a digimatic indicator (manufactured by Mitutoyo Corporation), dried in an oven at 100 ° C. for 1 hour, and dried at 90 to 100 μm. Is obtained.

【0032】次に、このパターンP1上に図4に示され
るパターンP2を前記と同様の手順で重ねて印刷・乾燥
する。更に同様に、図5及び図6に示されるパターンP
3、P4を、前記パターンP1、P2に順次重ねて印刷
・乾燥する。なお、図4〜6において符号h2〜h4は
いずれも空白(孔)を示す。
Next, the pattern P2 shown in FIG. 4 is overlaid on the pattern P1 in the same procedure as described above, and then printed and dried. Similarly, the pattern P shown in FIGS.
3 and P4 are sequentially printed on the patterns P1 and P2 and dried. 4 to 6, reference numerals h2 to h4 indicate blanks (holes).

【0033】次に、これらパターンP1〜P4によって
形成されたX線吸収膜14に、両面テープをラミネーシ
ョンパッカーで貼り合わせ、これによって粘着層16及
び離型紙18の層を形成し、最後に、図7に示されるパ
ターンP5で抜き加工して、図8(PETフィルム12
を除いた状態)、図9に示されるパターンのチャートシ
ールSを得る。
Next, a double-sided tape is adhered to the X-ray absorbing film 14 formed by the patterns P1 to P4 with a lamination packer, thereby forming a layer of an adhesive layer 16 and a release paper 18, and finally, as shown in FIG. 8 (PET film 12).
), And the chart seal S having the pattern shown in FIG. 9 is obtained.

【0034】このチャートシールSには、共通の離型紙
18上に粘着層16によって貼り付けられている複数の
テストチャート10が形成されていることになる。
On the chart seal S, a plurality of test charts 10 adhered to the common release paper 18 by the adhesive layer 16 are formed.

【0035】表1のに示される、UV硬化型樹脂から
なるバインダーを用いたインキによりテストチャートを
印刷する場合は、表1の組成でUVスクリーンX線遮
断インキを混練し、前記第1の実施の形態の例と同条件
でパターンP1〜パターンP4を重ねて印刷し、且つ各
パターン印刷毎に120Wメタルハライドランプによっ
て紫外線を照射して硬化させ、最終的に図8及び図9に
示されると同様のパターンのチャートシールSを得る。
When printing a test chart with an ink using a binder made of a UV-curable resin as shown in Table 1, a UV screen X-ray blocking ink having the composition shown in Table 1 was kneaded, and the first embodiment was performed. The pattern P1 to the pattern P4 are superimposed and printed under the same conditions as in the example of the embodiment, and each pattern is printed and cured by irradiating ultraviolet rays with a 120 W metal halide lamp, and finally the same as shown in FIGS. 8 and 9 To obtain a chart seal S having the following pattern.

【0036】上記テストチャート10を、図10に示さ
れるように、鼻骨骨折の患者(被写体)22の鼻頭頂部
にX線発生器24と平行に添着し、又、X線発生器24
からのX線ビーム26の照射領域内で患者22にかぶら
ない位置のベッド28上にもテストチャート10Bを添
着した。図10の符号30はX線フィルムを示す。
As shown in FIG. 10, the test chart 10 is attached to the top of the nose of a patient (subject) 22 having a nose bone fracture in parallel with an X-ray generator 24.
The test chart 10B was also attached to a bed 28 at a position where the X-ray beam 26 did not cover the patient 22 within the irradiation area. Reference numeral 30 in FIG. 10 indicates an X-ray film.

【0037】上記のような条件でX線撮影し、現像した
結果、図11に示されるようになり、X線フィルム30
を2波長スキャニングデンシトメータ32及びレーザス
キャナ34(図12参照)で測定し、コンピュータ50
における記憶装置51の、元画像データ領域51Aにデ
ジタル画像データとして保管した。ここで、現像したX
線フィルム30からは、前記テストチャート10におけ
る分割領域20A〜20Dの各濃度段階を読み取ること
ができた。
As a result of X-ray photography and development under the above conditions, the result is as shown in FIG.
Is measured with a two-wavelength scanning densitometer 32 and a laser scanner 34 (see FIG. 12), and a computer 50 is used.
Was stored as digital image data in the original image data area 51A of the storage device 51. Here, the developed X
From the line film 30, each density step of the divided areas 20A to 20D in the test chart 10 could be read.

【0038】又、図13に示されるように、ベット28
を水平面から角度θだけ傾けて前記と同一条件で、テス
トチャート10、10Bと共に患者22をX線撮影し
て、X線フィルム30を現像した結果、図14に示され
るようになった。
Further, as shown in FIG.
Was tilted from the horizontal plane by an angle θ, and under the same conditions as above, the patient 22 was X-ray photographed together with the test charts 10 and 10B, and the X-ray film 30 was developed. As a result, the result was as shown in FIG.

【0039】この図14の現像されたX線フィルム30
の画像も、前述と同様に、2波長スキャニングデンシト
メータ32及びレーザスキャナ34によって測定し、そ
のデジタル画像データを、前記記憶装置51における元
画像データ領域51Aに保存する。
The developed X-ray film 30 shown in FIG.
Is measured by the two-wavelength scanning densitometer 32 and the laser scanner 34 in the same manner as described above, and the digital image data is stored in the original image data area 51A in the storage device 51.

【0040】なお、図11及び図14において、患者2
2及びテストチャート10、10Bの画像の周囲の白抜
き部分は、実際は最も黒化している部分である。
Note that in FIG. 11 and FIG.
2 and the blank portions around the images of the test charts 10 and 10B are actually the blackest portions.

【0041】次に、前記元画像データ領域51Aに保存
されたX線撮影画像のデジタルデータに基づいて、画像
の濃淡の程度(コントラスト)、分解能(解像力)とを
補正し、同一基準にてデジタル保存する過程について説
明する。
Next, based on the digital data of the X-ray photographed image stored in the original image data area 51A, the degree of image density (contrast) and the resolution (resolution) are corrected, and the The process of saving will be described.

【0042】まず、記憶装置51に取り込まれているデ
ジタルデータをディスプレイ52上に表示し、該ディス
プレイ52上の、前記テストチャート10Bのチャート
画像11Bにおける、最もX線吸収係数の大きい前記分
割領域20Dに対応する部分を「白」として、その色値
を例えばポインティングデバイス54で指示して、コン
ピュータ50におけるチャート値読取部55に取り込
む。
First, the digital data taken in the storage device 51 is displayed on the display 52, and the divided area 20D having the largest X-ray absorption coefficient in the chart image 11B of the test chart 10B is displayed on the display 52. Is set to “white”, the color value thereof is indicated by, for example, the pointing device 54 and is taken into the chart value reading unit 55 of the computer 50.

【0043】同様に、チャート画像11Bにおける中心
の透過領域20E(X線を吸収しないか、又はほとんど
吸収しない領域)に対応する部分を「黒」としてポイン
ティングデバイス54により色値を読み取り、チャート
値読取部55に取り込み、更に、チャート画像11Bに
おける、分割領域20Bに対応する部分を「グレー」と
して、ポインティングデバイス54によりその色値を読
み取り、チャート値読取部55に取り込む。
Similarly, the color value is read by the pointing device 54 with the portion corresponding to the center transmission region 20E (the region that does not absorb or hardly absorbs X-rays) in the chart image 11B read by the pointing device 54, and the chart value is read. The color value is read by the pointing device 54, and the color value is read by the pointing device 54, and is captured by the chart value reading unit 55.

【0044】これら読み取ったテストチャート10Bの
色値のデジタルデータと、目標値設定部62によって、
予め設定され、且つ補正カーブ作成部56に入力されて
いる各色毎の目標値とにより、補正カーブ作成部56に
おいて、これら目標値又は設定値と読み取った色値とを
比較することにより、各色毎の補正値が算出される。
The digital data of the color values of the read test chart 10 B and the target value setting unit 62
The correction curve creation unit 56 compares these target values or set values with the read color values in accordance with the target values for each color which are set in advance and input to the correction curve creation unit 56. Is calculated.

【0045】この算出された補正値により、前記記憶装
置51の元画像データ領域51Aに取り込まれているデ
ジタルデータを、例えば画像データ変換部58により補
正して、記憶装置51の補正画像データ領域51Bに取
り込み、保存する。
With the calculated correction value, the digital data taken into the original image data area 51A of the storage device 51 is corrected by, for example, the image data conversion section 58, and the corrected image data area 51B of the storage device 51 is corrected. And save it.

【0046】図12の符号60は画像補正部全体、64
は3次元画像作成部、68は補正パラメーター保存部、
51Cは記憶装置51における3D画像データ領域をそ
れぞれ示す。
The reference numeral 60 in FIG.
Is a three-dimensional image creation unit, 68 is a correction parameter storage unit,
51C indicates a 3D image data area in the storage device 51, respectively.

【0047】前記目標値又は設定値は、テストチャート
10における各分割領域20A〜20DにおけるX線吸
収膜20及び透過領域20Eによる、既知のX線吸収係
数に基づいて設定する。
The target value or set value is set based on a known X-ray absorption coefficient by the X-ray absorbing film 20 and the transmission region 20E in each of the divided regions 20A to 20D in the test chart 10.

【0048】前記補正カーブ作成部56においては、図
15に示されるように、補正値を縦軸、読取値を横軸と
する2次元のグラフにおける補正カーブを作成し、この
補正カーブに基づいて前記取り込まれたデジタルデータ
を補正するものである。
As shown in FIG. 15, the correction curve creating section 56 creates a correction curve in a two-dimensional graph with the correction value on the vertical axis and the read value on the horizontal axis, and based on this correction curve. The digital data is corrected.

【0049】具体的には、図15に示されるように、ス
テップ101で被写体である患者22及びベッド28に
テストチャート10、10Bを貼り、ステップ102に
おいて被写体をテストチャート10、10Bと共に撮影
し、ステップ103において撮影した画像をデジタル化
して記憶し、ステップ104において記憶した画像をデ
ィスプレイ52上に表示し、ステップ105においてポ
インティングデバイス54により、該ディスプレイ52
上のチャート画像11Bにおける白、黒及びグレーの各
区画の色値を読み取る。
More specifically, as shown in FIG. 15, the test charts 10 and 10B are attached to the patient 22 and the bed 28 as the subject in step 101, and the subject is photographed together with the test charts 10 and 10B in step 102. In step 103, the captured image is digitized and stored. In step 104, the stored image is displayed on the display 52. In step 105, the display 52
The color values of the white, black, and gray sections in the upper chart image 11B are read.

【0050】次に、この読み取ったデジタルデータと、
前記テストチャート10Bにおける白、黒、グレーの予
め設定されている目標値のデジタルデータとにより、補
正カーブ作成部56において、図16に示されるよう
に、補正値を縦軸、読み取り値を横軸とする2次元のグ
ラフにおける補正カーブを形成し(ステップ106参
照)、この補正カーブによる補正パラメーターを、補正
パラメーター保存部68に保存する(ステップ106A
参照)と共に、補正カーブに基づいて前記取り込まれた
デジタルデータを補正し(ステップ107参照)、ステ
ップ108において補正されたデジタルデータを保存す
る。前記パラメーター保存部68に保存された補正デー
タは次回以降の同一条件での撮影の際に利用される。こ
の場合、テストチャートの被写体への貼付けが不要とな
る。
Next, the read digital data and
Based on the digital data of the preset target values of white, black, and gray in the test chart 10B, in the correction curve creation unit 56, as shown in FIG. A correction curve in a two-dimensional graph is formed (see step 106), and correction parameters based on this correction curve are stored in the correction parameter storage unit 68 (step 106A).
), The digital data taken in is corrected based on the correction curve (see step 107), and the digital data corrected in step 108 is stored. The correction data stored in the parameter storage unit 68 is used at the time of the next and subsequent imaging under the same conditions. In this case, it is not necessary to attach the test chart to the subject.

【0051】上記の過程を更に詳細に説明すると、図1
2に示される、目標値設定部62により補正カーブ作成
部56に予め白、黒、グレーの目標値OH 、OL 、OM
を設定しておき、ディスプレイ52上に画像を開いた
後、前述のように、チャート画像11B中の白、黒及び
グレーの領域でのRGB値IH (R) 、IH (G) 、I
H (B) 、IL (R) 、IL (G) 、IL (B) 、IM (R) 、I
M (G) 、IM (B) をそれぞれ読み取り、前記2次元グラ
フの点(0,0)−(IL ,OL )−(IM ,OM )−
(IH ,OH )−(255,255)を通る補正曲線を
形成し、この補正曲線に基づいて、既に取り込まれてい
るRGB各デジタルデータを変換し、このデジタルデー
タを記憶装置51に保存する。
The above process will be described in more detail.
The white, black, and gray target values O H , OL , and O M are previously stored in the correction curve creation unit 56 by the target value setting unit 62 shown in FIG.
After the image is opened on the display 52, the RGB values I H (R), I H (G), and I H in the white, black, and gray areas in the chart image 11B are set as described above.
H (B), IL (R), IL (G), IL (B), IM (R),
M (G), read I M a (B) respectively, the point of the two-dimensional graph (0,0) - (I L, O L) - (I M, O M) -
A correction curve passing (I H , O H )-(255, 255) is formed, and based on the correction curve, the RGB digital data that has already been captured is converted, and the digital data is stored in the storage device 51. I do.

【0052】例えば、透過領域、分割領域(以下領域)
のRGB値(入力値)が、 白:IH (R) =208 ,IH (G) =188 ,IH (B) =217 黒:IH (R) = 29 ,IH (G) = 16 ,IH (B) = 19 グレー:IH (R) =126 ,IH (G) = 90 ,IH (B) =
138 目標とするRGB値(目標値)が、 白:OH (R) =OH (G) =OH (B) =228 黒:OL (R) =OL (G) =OL (B) = 35 グレー:OM (R) =OM (G) =OM (B) =114 の場合、補正カーブは、 Rデータ:(0,0) −(29,35) −(126,114) −(208,228)
−(255,255) Gデータ:(0,0) −(16,35) −( 90,114)−(188,228)
−(255,255) Bデータ:(0,0) −(19,35) −(138,114) −(217,228)
−(255,255) となる(図16参照)。
For example, a transmission area, a division area (hereinafter, area)
RGB values (input values) are: white: IH (R) = 208, IH (G) = 188, IH (B) = 217 black: IH (R) = 29, IH (G) = 16, IH (B) = 19 Gray: IH (R) = 126, IH (G) = 90, IH (B) =
138 The target RGB values (target values) are: White: OH (R) = OH (G) = OH (B) = 228 Black: OL (R) = OL (G) = OL ( B) = 35 gray: for O M (R) = O M (G) = O M (B) = 114, the correction curve, R data: (0,0) - (29, 35) - (126,114) − (208,228)
− (255,255) G data: (0,0) − (16,35) − (90,114) − (188,228)
− (255,255) B data: (0,0) − (19,35) − (138,114) − (217,228)
− (255,255) (see FIG. 16).

【0053】前記図15に示されるフローチャート1の
補正方法の例を補正アルゴリズム1とすると、この補正
アルゴリズム1では、白、黒及びグーの各領域の読取り
デジタルデータを同時に補正変換しているので、白、
黒、グレー毎に順次補正した場合のように、先に行った
補正の結果が、後に行う補正によって影響を受けるとい
う欠点がないので理想的である。
Assuming that the example of the correction method of the flowchart 1 shown in FIG. 15 is a correction algorithm 1, in this correction algorithm 1, the read digital data of the white, black and goo areas are simultaneously corrected and converted. White,
This is ideal because there is no disadvantage that the result of the correction performed earlier is affected by the correction performed later, as in the case where the correction is performed sequentially for each of black and gray.

【0054】但し、コンピュータの構成あるいはソフト
ウェアの仕様等によって、アルゴリズム1の補正方法を
とることが困難な場合、次の補正アルゴリズム2〜4の
ような補正方法をとってもよい。
However, if it is difficult to use the correction method of Algorithm 1 due to the configuration of the computer or the specifications of software, correction methods such as the following correction algorithms 2 to 4 may be used.

【0055】これらアルゴリズム2、3及び4は、図1
7に示されるフローチャートの手順で取込まれたデジタ
ルデータを補正するものである。
These algorithms 2, 3 and 4 are shown in FIG.
This is to correct the digital data taken in the procedure of the flowchart shown in FIG.

【0056】まず、補正アルゴリズム2を図17を参照
して説明する。
First, the correction algorithm 2 will be described with reference to FIG.

【0057】予め、白、黒、グレーの目標値OH
L 、OM を設定しておき、ステップ201で画像を表
示し、ステップ2でチャート画像11B中の白領域のR
GB値I H (R) 、IH (G) 、IH (B) を読取り、ステッ
プ203で補正アルゴリズム(0,0)−(IH
H )−(255,255)を通るハイライト補正カー
ブを作成し(図18参照)、次にこのハイライト補正カ
ーブにより、ステップ204においてハイライト補正を
するようにRGB各データを変換し、次にステップ20
5において画像中の黒領域のRGB値IL (R) 、I
L (G) 、IL (B) を読取り、ステップ206において
(0,0)−(IL ,OL )−(255,255)を通
るシャドウ補正カーブを作成し、ステップ207におい
てシャドウ補正するようにRGB各データを変換する。
The target values O for white, black and gray are set in advance.H,
OL, OMIs set, and the image is displayed in step 201.
In step 2, R of the white area in the chart image 11B
GB value I H(R), IH(G), IHRead (B) and
In step 203, the correction algorithm (0,0)-(IH,
OH)-Highlight correction car passing through (255, 255)
(See FIG. 18) and then highlight highlight correction
The highlight correction in step 204.
RGB data is converted so that
5, the RGB value I of the black area in the imageL(R), I
L(G), IL(B) is read and in step 206
(0,0)-(IL, OL)-(255, 255)
Create a shadow correction curve
RGB data is converted to perform shadow correction.

【0058】次に画像中のグレー領域のRGB値I
M (R) 、IM (G) 、IM (B) を読取り、ステップ209
において(0,0)−(IM ,OM )−(255,25
5)を通るグレー補正カーブを作成しこれに基づき、グ
レー補正するようにRGB各データを変換して、ステッ
プ211で画像(デジタルデータ)を保存する。
Next, the RGB value I of the gray area in the image
M (R), IM (G) and IM (B) are read, and step 209 is executed.
In (0,0) - (I M, O M) - (255,25
A gray correction curve passing through 5) is created, and based on the gray correction curve, RGB data are converted so as to perform gray correction, and an image (digital data) is stored in step 211.

【0059】例えば、チャートのRGB値(入力値)
が、 白:IH (R) =208 ,IH (G) =188 ,IH (B) =217 黒:IH (R) = 29 ,IH (G) = 16 ,IH (B) = 19 グレー:IH (R) =126 ,IH (G) = 90 ,IH (B) =
138 目標とするRGB値(目標値)が、 白:OH (R) =OH (G) =OH (B) =228 黒:OL (R) =OL (G) =OL (B) = 35 グレー:OM (R) =OM (G) =OM (B) =114 の場合、ハイライト補正カーブは、 Rデータ:(0,0) −(208,228) −(255,255) Gデータ:(0,0) −(188,228) −(255,255) Bデータ:(0,0) −(217,228) −(255,255) シャドウ補正カーブは、 Rデータ:(0,0) −(29,35) −(255,255) Gデータ:(0,0) −(16,35) −(255,255) Bデータ:(0,0) −(19,35) −(255,255) グレー補正カーブは、 Rデータ:(0,0) −(126,114) −(255,255) Gデータ:(0,0) −( 90,114)−(255,255) Bデータ:(0,0) −(138,114) −(255,255) となる(図18〜20参照)。
For example, RGB values (input values) of a chart
White: IH (R) = 208, IH (G) = 188, IH (B) = 217 Black: IH (R) = 29, IH (G) = 16, IH (B) = 19 gray: I H (R) = 126, I H (G) = 90, I H (B) =
138 The target RGB values (target values) are: White: OH (R) = OH (G) = OH (B) = 228 Black: OL (R) = OL (G) = OL ( B) = 35 gray: for O M (R) = O M (G) = O M (B) = 114, highlight correction curve, R data: (0,0) - (208 and 228) - (255, 255) G data: (0,0)-(188,228)-(255,255) B data: (0,0)-(217,228)-(255,255) Shadow correction curve: R data: (0,0)-(29,35) )-(255,255) G data: (0,0)-(16,35)-(255,255) B data: (0,0)-(19,35)-(255,255) The gray correction curve is the R data: ( 0,0)-(126,114)-(255,255) G data: (0,0)-(90,114)-(255,255) B data: (0,0)-(138,114)-(255,255) 20).

【0060】補正アルゴリズム3は、前記ステップ20
9及び210において補正アルゴリズム2と異なる。補
正アルゴリズム3におけるステップ209は、読取った
グレー領域のRGB値IM (R) 、IM (G) 、IM (B) の
平均値Ave (IM )を計算して、(0,0)−IM 、A
ve (IM )−(255,255)を通るグレー補正カー
ブを作成し、ステップ210でこの補正カーブに基づい
てRGB各デジタルデータを変換し、ステップ211で
変換された画像のデジタルデータを保存するものであ
る。
The correction algorithm 3 is based on the step 20
9 and 210 differ from the correction algorithm 2. Step 209 in the correction algorithm 3 calculates an average value Ave (I M ) of the read RGB values I M (R), I M (G), and I M (B) of the gray area to obtain (0, 0). −I M , A
ve (I M ) —Create a gray correction curve passing through (255, 255), convert the RGB digital data based on the correction curve in step 210, and save the digital data of the converted image in step 211. Things.

【0061】チャートのRGB値(入力値)が、 白:IH (R) =208 ,IH (G) =188 ,IH (B) =217 黒:IH (R) = 29 ,IH (G) = 16 ,IH (B) = 19 グレー:IH (R) =126 ,IH (G) = 90 ,IH (B) =
138 目標とするRGB値(目標値)が、 白:OH (R) =OH (G) =OH (B) =228 黒:OL (R) =OL (G) =OL (B) = 35 の場合、ハイライト補正カーブは、 Rデータ:(0,0) −(208,228) −(255,255) Gデータ:(0,0) −(188,228) −(255,255) Bデータ:(0,0) −(217,228) −(255,255) シャドウ補正カーブは、 Rデータ:(0,0) −(29,35) −(255,255) Gデータ:(0,0) −(16,35) −(255,255) Bデータ:(0,0) −(19,35) −(255,255) グレー補正カーブは、 Rデータ:(0,0) −(126,118) −(255,255) Gデータ:(0,0) −( 90,118)−(255,255) Bデータ:(0,0) −(138,118) −(255,255) となる(図18〜20参照)。
The RGB values (input values) of the chart are: white: I H (R) = 208, I H (G) = 188, I H (B) = 217 black: I H (R) = 29, I H (G) = 16, I H (B) = 19 Gray: I H (R) = 126, I H (G) = 90, I H (B) =
138 The target RGB values (target values) are: White: OH (R) = OH (G) = OH (B) = 228 Black: OL (R) = OL (G) = OL ( When B) = 35, the highlight correction curve is as follows: R data: (0,0)-(208,228)-(255,255) G data: (0,0)-(188,228)-(255,255) B data: (0,0) , 0) − (217,228) − (255,255) The shadow correction curve is as follows: R data: (0,0) − (29,35) − (255,255) G data: (0,0) − (16,35) − ( 255,255) B data: (0,0)-(19,35)-(255,255) Gray correction curve is as follows: R data: (0,0)-(126,118)-(255,255) G data: (0,0)- (90,118) − (255,255) B data: (0,0) − (138,118) − (255,255) (see FIGS. 18 to 20).

【0062】次に、補正アルゴリズム4について説明す
る。
Next, the correction algorithm 4 will be described.

【0063】この補正アルゴリズムは、前記補正アルゴ
リズム2と、ステップ201〜208が同一であり、ス
テップ209では、次式によってO′M を算出する。
[0063] The correction algorithm, and the correction algorithm 2, step 201 to 208 are identical, in step 209, calculates the O 'M by the following equation.

【0064】 O′M (R) =ΣIM *(OM (R) /ΣOM ) O′M (G) =ΣIM *(OM (G) /ΣOM ) O′M (B) =ΣIM *(OM (B) /ΣOM ) 更に算出したO′M に基づいて、(0,0)−(IM
O′M )−(255,255)を通るグレー補正カーブ
を作成し、ステップ210において、前記グレー補正カ
ーブに基づいてRGB各デジタルデータを変換して、ス
テップ211において変換された画像のデジタルデータ
を保存するものである。
[0064] O 'M (R) = ΣI M * (O M (R) / ΣO M) O' M (G) = ΣI M * (O M (G) / ΣO M) O 'M (B) = ΣI M * (O M (B) / ΣO M ) Further, based on the calculated O ′ M , (0,0) − (I M ,
O ′ M )-(255, 255), a gray correction curve is created, and in step 210, each digital data of RGB is converted based on the gray correction curve, and the digital data of the image converted in step 211 is converted. What you save.

【0065】更に同様に、チャートのRGB値(入力
値)が、 白:IH (R) =208 ,IH (G) =188 ,IH (B) =217 黒:IH (R) = 29 ,IH (G) = 16 ,IH (B) = 19 グレー:IH (R) =126 ,IH (G) = 90 ,IH (B) =
138 目標とするRGB値(目標値)が、 白:OH (R) =OH (G) =OH (B) =228 黒:OL (R) =OL (G) =OL (B) = 35 グレー:OH (R) =OH (G) =OH (B) =114 の場合、ハイライト補正カーブは、 Rデータ:(0,0) −(208,228) −(255,255) Gデータ:(0,0) −(188,228) −(255,255) Bデータ:(0,0) −(217,228) −(255,255) シャドウ補正カーブは、 Rデータ:(0,0) −(29,35) −(255,255) Gデータ:(0,0) −(16,35) −(255,255) Bデータ:(0,0) −(19,35) −(255,255) グレー補正カーブは、 Rデータ:(0,0) −(126,118) −(255,255) Gデータ:(0,0) −( 90,118)−(255,255) Bデータ:(0,0) −(138,118) −(255,255) と同じ結果となる(図18〜20参照)。
Similarly, the RGB values (input values) of the chart are as follows: white: I H (R) = 208, I H (G) = 188, I H (B) = 217 black: I H (R) = 29, IH (G) = 16, IH (B) = 19 Gray: IH (R) = 126, IH (G) = 90, IH (B) =
138 The target RGB values (target values) are: White: OH (R) = OH (G) = OH (B) = 228 Black: OL (R) = OL (G) = OL ( B) = 35 Gray: When O H (R) = O H (G) = O H (B) = 114, the highlight correction curve is R data: (0,0)-(208,228)-(255,255) G data: (0,0)-(188,228)-(255,255) B data: (0,0)-(217,228)-(255,255) Shadow correction curve: R data: (0,0)-(29,35) )-(255,255) G data: (0,0)-(16,35)-(255,255) B data: (0,0)-(19,35)-(255,255) The gray correction curve is the R data: ( 0,0)-(126,118)-(255,255) G data: (0,0)-(90,118)-(255,255) B data: (0,0)-(138,118)-Same result as (255,255) 18 to 20).

【0066】前記補正アルゴリズム2、3及び4は、
白、黒、グレーの順で、即ちハイライト補正、シャドウ
補正、グレー補正を独立で行うために、先の補正の結果
が後の補正によって影響を受ける場合があるが、実際上
大きな不都合はない。
The correction algorithms 2, 3, and 4 are as follows:
In order to perform white correction, black correction, and gray correction independently, that is, highlight correction, shadow correction, and gray correction, the result of the previous correction may be affected by the subsequent correction, but there is no significant inconvenience in practice. .

【0067】又、上記の如く、補正アルゴリズム2、
3、4において、補正をハイライト(白)、シャドウ
(黒)、グレーの順にすると、白は後の黒、グレーの補
正によって変化する可能性があるが、見た目の変化が少
ないので、最も良好な補正結果が得られた。
As described above, the correction algorithm 2
When the correction is performed in the order of highlight (white), shadow (black), and gray in 3 and 4, white may change due to the subsequent correction of black and gray. Results were obtained.

【0068】上記補正アルゴリズム1〜4では、画像デ
ータとしてRGBデータを例としてあげているが、これ
はCMYデータでもよい。
In the above correction algorithms 1 to 4, RGB data is taken as an example of image data, but this may be CMY data.

【0069】次に、L* * * やL* * * 、YI
Q等、明度(輝度)と色差が分離した表色系の画像デー
タに対する補正アルゴリズム5を、L* * * データ
を例に、図21のフローチャートを参照して説明する。
Next, L * a * b * , L * u * v * , YI
A correction algorithm 5 for color system image data in which lightness (luminance) and color difference are separated such as Q will be described with reference to a flowchart of FIG. 21 using L * a * b * data as an example.

【0070】予め、白、黒、グレーの目標OH 、OL
M を設定しておき、ステップ301で画像を表示し、
ステップ302でチャート画像11B中の白、黒領域の
*値IH (L) 、IL (L) 及びグレー領域のL* *
* 値IM (L) 、IM (a) 、I M (b) をそれぞれ読み取
り、ステップ303で(0,0)−(IL (L) ,O
L (L) )−(IM (L) ,OM (L) )−(IH (L) ,OH
(L) )−(100,100)を通るL* 補正カーブを作
成し、ステップ304でこの補正カーブに基づいて、既
に取り込まれているL* データを変換し、ステップ30
5で補正値ΔM(a) 、ΔM(b) を次式によって算出す
る。
The white, black and gray targets OH, OL,
OMIs set, and an image is displayed in step 301,
In step 302, the white and black areas in the chart image 11B are
L*Value IH(L), IL(L) and L in the gray area*a*b
*Value IM(L), IM(a), I M(b)
In step 303, (0,0)-(IL(L), O
L(L))-(IM(L), OM(L))-(IH(L), OH
(L))-L passing through (100,100)*Create a correction curve
At step 304 based on the correction curve.
L taken in*Transform the data, step 30
In step 5, the correction values ΔM (a) and ΔM (b) are calculated by the following equations.
You.

【0071】ΔM(a) =OM (a)−IM (a) ΔM(b) =OM (b)−IM (b) そして、ステップ306で既に取り込まれているa*
びb* データに前記補正値ΔM(a) 、ΔM(B) をそれぞ
れ加算し、ステップ307で変換された画像のデジタル
データを保存するものである。
[0071] ΔM (a) = O M ( a) -I M (a) ΔM (b) = O M (b) -I M (b) Then, already captured in step 306 a * and b * The correction values ΔM (a) and ΔM (B) are added to the data, and the digital data of the image converted in step 307 is stored.

【0072】例えば、チャートのL* * * 値(入力
値)が、 白:IH (L) =84 黒:IL (L) =12 グレー:IM (L) =55,IM (a) =24,IM (b) =
−19 目標とするL* * * 値(目標値)が、 白:OH (L) =93 黒:OL (L) =22 グレー:OM (L) =58,OM (a) =0,OM (b) =0 の場合、L* 補正カーブは、 (0,0)−(12,22)−(55,58)−(8
4,93)−(100,100) となり(図22参照)、補正値ΔM(a) 、ΔM(b) は、 ΔM(a) =0−24=−24 ΔM(b) =0−(−19)=19 となる。
For example, the L * a * b * values (input values) of the chart are as follows: white: I H (L) = 84 black: I L (L) = 12 gray: I M (L) = 55, I M (a) = 24, I M (b) =
L * a * b * values to -19 target (target value), white: O H (L) = 93 Black: O L (L) = 22 gray: O M (L) = 58 , O M (a ) = 0, if the O M (b) = 0, L * correction curve, (0,0) - (12,22) - (55, 58) - (8
4,93) − (100,100) (see FIG. 22), and the correction values ΔM (a) and ΔM (b) are ΔM (a) = 0−24 = −24 ΔM (b) = 0 − (− 19) = 19.

【0073】上記補正アルゴリズム5では、画像データ
としてL* * * を例としてあげているが、これはL
* * * やYIQ等、明度(輝度)と色差が分離した
表色系の画像データでもよい。
In the correction algorithm 5, L * a * b * is taken as an example of image data.
* u * v * , YIQ, or other color system image data in which lightness (luminance) and color difference are separated may be used.

【0074】又、前記補正方法では、ディスプレイ上に
表示された画像毎に補正カーブを作成して色値を補正す
るようにしているが、撮影条件が同一の複数の画像の場
合は、同一の補正カーブを利用することができる。
In the above-described correction method, a correction curve is created for each image displayed on the display to correct the color value. However, when a plurality of images have the same shooting conditions, the same value is used. A correction curve can be used.

【0075】従って、この場合、初回の補正操作のとき
作成し、保存した補正カーブを利用することにより、図
15、図17及び図21のフローチャートにおける補正
カーブ作成のステップを省略して、補正作業を短縮化す
ることができる。
Therefore, in this case, by using the correction curve created and stored at the time of the first correction operation, the correction curve creation steps in the flowcharts of FIGS. 15, 17 and 21 can be omitted, and the correction work can be performed. Can be shortened.

【0076】次に、画像サイズを補正し、且つ解像度を
調整する過程について説明する。
Next, the process of correcting the image size and adjusting the resolution will be described.

【0077】図23に示されるように、ステップ401
で、前記元画像データ51A又は補正画像データ領域5
1Bの画像データにより、画像をディスプレイ52上に
表示して、ステップ402において、ディスプレイ52
上の画像におけるチャート画像11Aの最も大きい直径
部分に直線を引く。
As shown in FIG.
The original image data 51A or the corrected image data area 5
1B, the image is displayed on the display 52, and in step 402, the display 52
A straight line is drawn on the largest diameter portion of the chart image 11A in the above image.

【0078】次に、ステップ403において、その直線
上でのチャート画像11Aの距離Dを測定、記憶し、ス
テップ404において、この記憶した距離Dと実際のテ
ストチャート10の直径dとを比較して、d/Dを、演
算部66において演算し、ステップ405において、記
憶する。
Next, in step 403, the distance D of the chart image 11A on the straight line is measured and stored, and in step 404, the stored distance D is compared with the actual test chart 10 diameter d. , D / D are calculated in the calculation unit 66 and stored in step 405.

【0079】このD/dは、患者22の患部近くにおけ
る水平方向の画像上の寸法と実際の寸法との比であるの
で、テストチャート10近傍の患部の実寸は、画像上の
寸法にD/dの逆数をかけることによって求めることが
できる。
Since D / d is the ratio of the size on the horizontal image near the affected part of the patient 22 to the actual size, the actual size of the affected part near the test chart 10 is represented by D / D. It can be obtained by multiplying the reciprocal of d.

【0080】又、複数の画像において、このD/dの数
値を統一すれば、異なる条件で撮影された画像データも
同一の縮尺で比較することができる。
If the numerical value of D / d is unified in a plurality of images, image data shot under different conditions can be compared at the same scale.

【0081】上記演算の結果のD/dに基づいて、ステ
ップ406において、画像データ変換部58にて画像サ
イズ補正を行い、補正された画像データは、前記計算値
D/dと共に、記憶装置51の補正画像データ領域51
Bに取り込まれる。
In step 406, based on the result D / d of the above operation, the image size is corrected by the image data converter 58. The corrected image data is stored in the storage device 51 together with the calculated value D / d. Corrected image data area 51
B is taken in.

【0082】次に、ステップ408において、予め演算
部66に入力されている解像度データを、前記画像サイ
ズ補正に伴って、演算部66により演算して、解像度を
複数の画像データにおいて同一となるように調整する。
Next, in step 408, the resolution data previously input to the calculation unit 66 is calculated by the calculation unit 66 in accordance with the image size correction so that the resolution becomes the same for a plurality of image data. Adjust to

【0083】即ち、(補正前の解像度)×(チャート画
像11A上で測定した直径の距離D)/(サイズ補正後
の前記直径部分の長さD1 )=補正後の解像度とする。
That is, (resolution before correction) × (distance D of diameter measured on chart image 11A) / (length D 1 of diameter portion after size correction) = resolution after correction.

【0084】具体的には、補正前の画像の解像度=72
pixels/inchとすると、チャート画像上の直径D=7.
48cmのとき、サイズ補正後のチャート画像の直径D
1 =4cmにしたい場合は、72×7.48/4=13
4.64に変更する。
Specifically, the resolution of the image before correction = 72
Assuming pixels / inch, the diameter D on the chart image = 7.
At 48 cm, the diameter D of the chart image after size correction
If you want 1 = 4 cm, 72 × 7.48 / 4 = 13
Change to 4.64.

【0085】この解像度を、前記画像データ変換部58
を介して、補正画像データ領域51Bにある前記サイズ
補正された画像データと共に取り込んでおけば、所定の
寸法に画像を表示する際に、その補正値に応じて、解像
度を自動的に変更することができる(ステップ409参
照)。
The resolution is determined by the image data converter 58.
If the image data is captured together with the size-corrected image data in the corrected image data area 51B via the, when the image is displayed in a predetermined size, the resolution is automatically changed according to the correction value. (See step 409).

【0086】又、ステップ410に示されるように、前
記解像度の補正の際における補正パラメーターを前記補
正パラメーター保存部68に独立して保存しておけば、
次回以降の同一条件での撮影の際に、テストチャートを
被写体に貼る過程、及びその画像データの読み取りが不
要となる。
Further, as shown in step 410, if the correction parameters for the resolution correction are stored independently in the correction parameter storage section 68,
In the next and subsequent shootings under the same conditions, the process of attaching the test chart to the subject and the reading of the image data are not required.

【0087】次に、複数のX線撮影画像から、3次元画
像を作成する場合について説明する。
Next, the case where a three-dimensional image is created from a plurality of X-ray images will be described.

【0088】図24に示されるように、ステップ501
において、補正画像データ領域51Bからサイズ補正さ
れた画像データを読み出し、ディスプレイ52に表示す
る。
As shown in FIG. 24, step 501
, The size-corrected image data is read out from the corrected image data area 51B and displayed on the display 52.

【0089】ステップ502においてディスプレイ52
での画像中におけるチャート画像11の短径方向に直線
を引き、ステップ503で、その直線上の短径の大きさ
Yを計測して記憶する、実物のテストチャート10の大
きさをd、画像上でのチャート画像11の中心から、患
部までの距離をZとすると、実際の患者22の患部とテ
ストチャート10の中央部からの深さ方向の距離は、Z
×d/Y tanθとなるので、演算部66でこれを計算す
る(ステップ504)。
In step 502, the display 52
In step 503, a straight line is drawn in the minor axis direction of the chart image 11 in the image at step 503, and the size Y of the minor axis on the straight line is measured and stored. Assuming that the distance from the center of the chart image 11 to the affected part is Z, the actual distance in the depth direction from the affected part of the patient 22 to the center of the test chart 10 is Z
Xd / Y tan θ, which is calculated by the calculation unit 66 (step 504).

【0090】従って、図10における、ベッド28を傾
けない状態での患者22の患部とテストチャート10の
上下方向の距離が判明することになる。
Therefore, the vertical distance between the affected part of the patient 22 and the test chart 10 in a state where the bed 28 is not tilted in FIG. 10 is determined.

【0091】このようにして得られた患部の高さ方向の
距離(深さ方向の距離)と、水平方向の距離とが判明す
ることによって、コンピュータ50における3次元画像
作成部64において、3Dグラフィックソフト上で数式
処理し(ステップ505)、オブジェクト表現すると、
簡易に3D画像としてディスプレイ52やプリンタ(図
示省略)に表示することができた(ステップ506)。
この3D画像は、3Dグラフィックソフト上で簡単に回
転等ができるので、患部を特定することも可能である。
The distance in the height direction (distance in the depth direction) and the distance in the horizontal direction of the affected part obtained in this manner are determined, and the three-dimensional image forming unit 64 in the computer 50 determines the 3D graphic. Formula processing is performed on the software (step 505).
It could be simply displayed on the display 52 or a printer (not shown) as a 3D image (step 506).
Since the 3D image can be easily rotated on the 3D graphic software, the affected part can be specified.

【0092】なお、上記実施の形態の例は、いずれも金
属微粉末を含むインキにより、パターンを印刷してテス
トチャートを形成したものであるが、本発明はこれに限
定されるものでなく、X線吸収係数が段階的に異なるよ
うに膜厚を調整したテストチャートであればよい。
In each of the above embodiments, the test chart is formed by printing a pattern with an ink containing fine metal powder, but the present invention is not limited to this. Any test chart may be used as long as the film thickness is adjusted so that the X-ray absorption coefficient varies stepwise.

【0093】従って、例えば、X線吸収係数が高い金属
板上にフォトレジストを、撮影方向及び大きさの特定が
できるパターン状に微細印刷し、該パターンを、エッチ
ング液の腐食強度、エッチング時間、及び回数で制御し
て、パターンにおける分割領域のX線吸収係数が段階的
に異なるようにしてもよい。
Therefore, for example, a photoresist is finely printed on a metal plate having a high X-ray absorption coefficient in a pattern shape capable of specifying a photographing direction and a size, and the pattern is subjected to etching strength, etching time, And the number of times, the X-ray absorption coefficients of the divided regions in the pattern may be changed stepwise.

【0094】[0094]

【発明の効果】本発明は上記のように構成したので、X
線撮影時に、被写体と同一画面に、撮影中・撮影後の大
きさ、方向及びX線吸収段階の特定ができるので、正確
なX線撮影画像の補正をすることができ、又、複数のX
線撮影画像データにおいて、濃淡の程度や、サイズを統
一して保存、表示することができるという優れた効果を
有する。
The present invention is constructed as described above.
At the time of X-ray imaging, the size, direction, and X-ray absorption stage during and after imaging can be specified on the same screen as the subject, so that accurate X-ray imaging images can be corrected.
In the line image data, there is an excellent effect that the degree of shading and the size can be unified and stored and displayed.

【図面の簡単な説明】[Brief description of the drawings]

【図1】本発明の実施に用いるX線撮影画像の補正用テ
ストチャートの例を示す平面図
FIG. 1 is a plan view showing an example of a test chart for correcting an X-ray radiographic image used in an embodiment of the present invention.

【図2】図1のII−II線に沿う拡大断面図FIG. 2 is an enlarged sectional view taken along the line II-II of FIG.

【図3】テストチャートを重ね印刷する場合のパターン
を示す平面図
FIG. 3 is a plan view showing a pattern when a test chart is overprinted.

【図4】テストチャートを重ね印刷する場合のパターン
を示す平面図
FIG. 4 is a plan view showing a pattern when a test chart is overprinted.

【図5】テストチャートを重ね印刷する場合のパターン
を示す平面図
FIG. 5 is a plan view showing a pattern when a test chart is overprinted.

【図6】テストチャートを重ね印刷する場合のパターン
を示す平面図
FIG. 6 is a plan view showing a pattern when a test chart is overprinted.

【図7】テストチャート形成過程で用いる抜き型を示す
平面図
FIG. 7 is a plan view showing a punch used in a test chart forming process.

【図8】完成したチャートシールを示す平面図FIG. 8 is a plan view showing a completed chart seal.

【図9】同側面図FIG. 9 is a side view of the same.

【図10】上記テストチャートを用いて患者をX線撮影
する形態を示す略示正面図
FIG. 10 is a schematic front view showing a form in which a patient is X-rayed using the test chart.

【図11】同患者をX線撮影した現像フィルムを示す平
面図
FIG. 11 is a plan view showing a developed film obtained by radiographing the patient.

【図12】本発明方法を実施するためのシステムを示す
ブロック図
FIG. 12 is a block diagram showing a system for implementing the method of the present invention.

【図13】患者のベッドを傾けてX線撮影する形態を示
す略示正面図
FIG. 13 is a schematic front view showing a form in which a patient's bed is tilted to perform X-ray imaging.

【図14】同患者をX線撮影した現像フィルムを示す平
面図
FIG. 14 is a plan view showing a developed film obtained by radiographing the patient.

【図15】本発明の画像を補正する手順を示すフローチ
ャート
FIG. 15 is a flowchart showing a procedure for correcting an image according to the present invention.

【図16】図15のフローチャートで示された補正アル
ゴリズムによって作成されたRGB補正カーブを示す線
FIG. 16 is a diagram showing an RGB correction curve created by the correction algorithm shown in the flowchart of FIG.

【図17】本発明に係る補正方法の他の補正アルゴリズ
ムを示すフローチャート
FIG. 17 is a flowchart showing another correction algorithm of the correction method according to the present invention.

【図18】図17に示されるフローチャートのアルゴリ
ズムによって作成されたハイライト補正カーブを示す線
FIG. 18 is a diagram showing a highlight correction curve created by the algorithm of the flowchart shown in FIG. 17;

【図19】同シャドウ補正カーブを示す線図FIG. 19 is a diagram showing the shadow correction curve;

【図20】同グレー補正カーブを示す線図FIG. 20 is a diagram showing the gray correction curve;

【図21】本発明に係る補正方法の更に他の補正アルゴ
リズムを示すフローチャート
FIG. 21 is a flowchart showing still another correction algorithm of the correction method according to the present invention.

【図22】図21に示されるフローチャートのアルゴリ
ズムによって作成された補正カーブを示す線図
FIG. 22 is a diagram showing a correction curve created by the algorithm of the flowchart shown in FIG. 21;

【図23】本発明方法により、X線撮影デジタル画像の
サイズを補正する手順を示すフローチャート
FIG. 23 is a flowchart showing a procedure for correcting the size of a radiographic digital image according to the method of the present invention.

【図24】本発明方法により、3次元立体画像を作成す
る手順を示すフローチャート
FIG. 24 is a flowchart showing a procedure for creating a three-dimensional stereoscopic image according to the method of the present invention.

【符号の説明】[Explanation of symbols]

10、10B…テストチャート 11、11B…チャートパターン 12…PETフィルム 14…X線吸収膜 16…粘着層 18…離型紙 20A〜20D…分割領域 20E…透過領域 22…患者 22A…頭部 24…X線発生器 26…X線ビーム 28…ベッド 30…X線フィルム P1〜P4…パターン 50…コンピュータ 51…記憶装置 52…ディスプレイ 54…ポインティングデバイス 56…補正カーブ作成部 58…画像データ変換部 64…三次元画像作成部 66…演算部 10, 10B ... test chart 11, 11B ... chart pattern 12 ... PET film 14 ... X-ray absorbing film 16 ... adhesive layer 18 ... release paper 20A to 20D ... divided area 20E ... transmission area 22 ... patient 22A ... head 24 ... X Ray generator 26 X-ray beam 28 Bed 30 X-ray film P1 to P4 Pattern 50 Computer 51 Storage device 52 Display device 56 Pointing device 56 Correction curve generator 58 Image data converter 64 Tertiary Original image creation unit 66 ... Calculation unit

Claims (6)

【特許請求の範囲】[Claims] 【請求項1】X線を吸収しない、又は、X線吸収係数が
低いシート上に、X線吸収係数が高い金属を含み、その
X線吸収係数が所定規格となる膜厚で、且つ、大きさ及
び撮影方向を特定できるパターンのX線吸収膜を形成
し、このX線吸収膜を被って、被写体へのシート貼付用
の粘着層を形成してなるX線撮影評価・補正用テストチ
ャートを、前記粘着層により被写体及びX線照射領域内
であって、被写体にかぶらない位置に貼着して、これら
のテストチャートを被写体と共にX線撮影する過程と、 撮影した画像をデジタルデータとして記憶装置に取り込
む過程と、 前記記憶装置に取り込まれたデジタルデータを、ディス
プレイ上に表示する過程と、 ディスプレイ上における前記被写体にかぶっていないテ
ストチャートのチャート画像の濃度を読み取る過程と、 読み取った濃度のデジタルデータを標準値に合わせるよ
うに補正する過程と、 この補正値に基づき画像全体のデジタルデータを補正す
る過程と、 補正したデジタルデータを前記記憶装置に保存する過程
と、を有してなるX線撮影画像のデジタルデータ補正、
保存方法。
1. A sheet which does not absorb X-rays or contains a metal having a high X-ray absorption coefficient on a sheet having a low X-ray absorption coefficient, and has a film thickness and a large size such that the X-ray absorption coefficient becomes a predetermined standard. An X-ray imaging evaluation / correction test chart is formed by forming an X-ray absorption film having a pattern capable of specifying the direction and the imaging direction, and forming an adhesive layer for attaching a sheet to a subject by covering the X-ray absorption film. Attaching the test chart to the subject and the X-ray irradiation area in a position that does not cover the subject by the adhesive layer, and X-ray photographing the test chart together with the subject; and storing the photographed image as digital data. Capturing the digital data captured by the storage device on a display; and displaying the digital image of the test chart not covering the subject on the display. Reading the digital data, correcting the digital data of the read density to the standard value, correcting the digital data of the entire image based on the corrected value, and storing the corrected digital data in the storage device. Digital data correction of the radiographic image, comprising:
Preservation method.
【請求項2】請求項1において、前記チャート画像の濃
度の読み取りを白、黒、グレーの3色の色値を読み取る
ことにより行い、これらのうち少なくとも白及び黒の色
値の目標値を予め設定しておき、前記ディスプレイ上の
画像から読み取った色値のデジタルデータと前記目標値
とから補正値を算出することを特徴とするX線撮影画像
のデジタルデータ補正、保存方法。
2. The method according to claim 1, wherein the reading of the density of the chart image is performed by reading three color values of white, black, and gray, and at least the target values of the color values of white and black are determined in advance. A digital data correction and storage method for an X-ray photographed image, wherein a correction value is calculated from digital data of color values read from an image on the display and the target value.
【請求項3】請求項1において、ディスプレイ上におけ
る前記被写体に貼着したテストチャートのチャート画像
の大きさのデジタルデータに基づくアウトプットの解像
度と、該テストチャートの、前記と同一個所の寸法を前
記記憶装置に登録したときのデジタルデータに基づくア
ウトプットの解像度とを一致させて、該テストチャート
近傍における被写体の寸法のデジタルデータを補正する
過程と、 補正したデジタルデータを前記記憶装置に保存する過程
と、を有してなるX線撮影画像のデジタルデータ補正、
保存方法。
3. The method according to claim 1, wherein the output resolution based on the digital data of the size of the chart image of the test chart attached to the subject on the display and the dimensions of the same portion of the test chart as described above are determined. Correcting the digital data of the dimensions of the subject in the vicinity of the test chart by matching the output resolution based on the digital data registered in the storage device, and storing the corrected digital data in the storage device Digital data correction of the radiographic image, comprising:
Preservation method.
【請求項4】請求項1、2又は3において、前記被写体
を前記テストチャートと共に異なる角度で複数回X線撮
影する過程と、 各X線撮影による画像をデジタルデータとして記憶装置
に取り込む過程と、 これら異なる撮影角度の複数の画像のデジタルデータ
を、3次元グラフィックソフトで数式処理して、前記被
写体の立体画像を得て、そのデジタルデータを前記記憶
装置に記憶させる過程と、を有してなることを特徴とす
るX線撮影画像のデジタルデータ補正、保存方法。
4. The method according to claim 1, 2 or 3, wherein the subject is X-ray-photographed a plurality of times together with the test chart at different angles; Mathematically processing digital data of the plurality of images at these different photographing angles with three-dimensional graphic software to obtain a three-dimensional image of the subject, and storing the digital data in the storage device. A method for correcting and storing digital data of an X-ray photographed image.
【請求項5】請求項1乃至4のいずれかにおいて、前記
テストチャートにおけるX線吸収膜は、複数の領域に分
割され、且つ、各分割領域における膜厚が異なるX線吸
収係数となるようにされたことを特徴とするX線撮影画
像のデジタルデータ補正、保存方法。
5. The X-ray absorbing film according to claim 1, wherein the X-ray absorbing film in the test chart is divided into a plurality of regions, and the film thickness in each divided region has a different X-ray absorption coefficient. A method for correcting and storing digital data of an X-ray photographed image.
【請求項6】請求項1乃至5のいずれかにおいて、前記
補正値を、前記補正したデジタルデータと独立して保存
する過程を有してなり、次回以降の同一撮影条件の撮像
時に、前記テストチャートの被写体への貼付過程、及
び、テストチャートのチャート画像の濃度読み取り過程
の省略可能としたことを特徴とするX線撮影画像のデジ
タルデータ補正、保存方法。
6. The method according to claim 1, further comprising the step of storing the correction value independently of the corrected digital data, and performing the test at the next and subsequent imaging under the same imaging condition. A digital data correction and storage method for an X-ray photographed image, wherein a process of attaching a chart to a subject and a process of reading the density of a chart image of a test chart can be omitted.
JP07574597A 1997-03-27 1997-03-27 Digital data correction and storage method and apparatus for X-ray image Expired - Lifetime JP3863963B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP07574597A JP3863963B2 (en) 1997-03-27 1997-03-27 Digital data correction and storage method and apparatus for X-ray image

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP07574597A JP3863963B2 (en) 1997-03-27 1997-03-27 Digital data correction and storage method and apparatus for X-ray image

Publications (2)

Publication Number Publication Date
JPH10262959A true JPH10262959A (en) 1998-10-06
JP3863963B2 JP3863963B2 (en) 2006-12-27

Family

ID=13585128

Family Applications (1)

Application Number Title Priority Date Filing Date
JP07574597A Expired - Lifetime JP3863963B2 (en) 1997-03-27 1997-03-27 Digital data correction and storage method and apparatus for X-ray image

Country Status (1)

Country Link
JP (1) JP3863963B2 (en)

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7379529B2 (en) 2000-08-29 2008-05-27 Imaging Therapeutics, Inc. Methods and devices for quantitative analysis of x-ray images
US7660453B2 (en) 2000-10-11 2010-02-09 Imaging Therapeutics, Inc. Methods and devices for analysis of x-ray images
US7676023B2 (en) 2000-08-29 2010-03-09 Imaging Therapeutics, Inc. Methods and devices for quantitative analysis of x-ray images
JP2012120906A (en) * 2012-03-28 2012-06-28 Konica Minolta Holdings Inc Information processing apparatus and program
JP2014000390A (en) * 2012-06-11 2014-01-09 Planmeca Oy Dental surface model
JP2014054301A (en) * 2012-09-11 2014-03-27 Hitachi Aloka Medical Ltd Medical x-ray measuring apparatus
US8913818B2 (en) 2000-10-11 2014-12-16 Imatx, Inc. Methods and devices for evaluating and treating a bone condition based on X-ray image analysis
US8939917B2 (en) 2009-02-13 2015-01-27 Imatx, Inc. Methods and devices for quantitative analysis of bone and cartilage
US8965087B2 (en) 2004-09-16 2015-02-24 Imatx, Inc. System and method of predicting future fractures
US8965075B2 (en) 2002-09-16 2015-02-24 Imatx, Inc. System and method for predicting future fractures
JP2015142719A (en) * 2013-12-25 2015-08-06 東芝メディカルシステムズ株式会社 Medical image processor, x-ray diagnostic apparatus, phantom and medical image processing program
US9155501B2 (en) 2003-03-25 2015-10-13 Imatx, Inc. Methods for the compensation of imaging technique in the processing of radiographic images
US9267955B2 (en) 2001-05-25 2016-02-23 Imatx, Inc. Methods to diagnose treat and prevent bone loss

Cited By (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7676023B2 (en) 2000-08-29 2010-03-09 Imaging Therapeutics, Inc. Methods and devices for quantitative analysis of x-ray images
US7379529B2 (en) 2000-08-29 2008-05-27 Imaging Therapeutics, Inc. Methods and devices for quantitative analysis of x-ray images
US7660453B2 (en) 2000-10-11 2010-02-09 Imaging Therapeutics, Inc. Methods and devices for analysis of x-ray images
US9767551B2 (en) 2000-10-11 2017-09-19 Imatx, Inc. Methods and devices for analysis of x-ray images
US9275469B2 (en) 2000-10-11 2016-03-01 Imatx, Inc. Methods and devices for evaluating and treating a bone condition on x-ray image analysis
US8913818B2 (en) 2000-10-11 2014-12-16 Imatx, Inc. Methods and devices for evaluating and treating a bone condition based on X-ray image analysis
US9267955B2 (en) 2001-05-25 2016-02-23 Imatx, Inc. Methods to diagnose treat and prevent bone loss
US9460506B2 (en) 2002-09-16 2016-10-04 Imatx, Inc. System and method for predicting future fractures
US8965075B2 (en) 2002-09-16 2015-02-24 Imatx, Inc. System and method for predicting future fractures
US9155501B2 (en) 2003-03-25 2015-10-13 Imatx, Inc. Methods for the compensation of imaging technique in the processing of radiographic images
US8965087B2 (en) 2004-09-16 2015-02-24 Imatx, Inc. System and method of predicting future fractures
US8939917B2 (en) 2009-02-13 2015-01-27 Imatx, Inc. Methods and devices for quantitative analysis of bone and cartilage
JP2012120906A (en) * 2012-03-28 2012-06-28 Konica Minolta Holdings Inc Information processing apparatus and program
JP2014000390A (en) * 2012-06-11 2014-01-09 Planmeca Oy Dental surface model
JP2014054301A (en) * 2012-09-11 2014-03-27 Hitachi Aloka Medical Ltd Medical x-ray measuring apparatus
JP2015142719A (en) * 2013-12-25 2015-08-06 東芝メディカルシステムズ株式会社 Medical image processor, x-ray diagnostic apparatus, phantom and medical image processing program
US10169845B2 (en) 2013-12-25 2019-01-01 Toshiba Medical Systems Corporation Medical image processing apparatus, x-ray diagnostic apparatus, phantom, and medical image processing method

Also Published As

Publication number Publication date
JP3863963B2 (en) 2006-12-27

Similar Documents

Publication Publication Date Title
US5008947A (en) Method and apparatus for correcting extension rates of images
US5357549A (en) Method of dynamic range compression of an X-ray image and apparatus effectuating the method
US7092581B2 (en) Balancing areas of varying density in a digital image
JPH10262959A (en) Method for correcting and preserving digital data of x-ray radiographic images
EP0170270A2 (en) Method of adjusting radiation image read-out conditions
CN1057114A (en) The X ray image pick-up method
JP6853376B2 (en) How to reconstruct a 2D image from multiple X-ray images
CN115018720A (en) Geometric correction method and device for multilayer flat panel detector
JPH1085206A (en) Method for operating medical x-ray device
US7394924B2 (en) Scatter correction in scanning imaging systems
US4149250A (en) Method and apparatus for producing a cross sectional image of a body
Reinstein et al. Contrast enhancement of high-energy radiotherapy films
JPH11253429A (en) Method and device for taking and displaying radiograph
Wong et al. Calibration procedure in dual‐energy scanning using the basis function technique
JP4318763B2 (en) Test chart for ionizing radiation imaging evaluation / correction and manufacturing method thereof
CN113587810A (en) Method and device for generating light source position
JPS6139573Y2 (en)
JP3356356B2 (en) X-ray equipment
CN108742664B (en) Method, system, device and storage medium for calculating tissue density of mammary gland
JP2004093368A (en) X-ray mask manufacturing method for x-ray fluoroscopic inspection device
JP3131503B2 (en) Sensitivity compensation intensifying screen and method of manufacturing the same
JP2614219B2 (en) Manufacturing method of sensitivity compensation intensifying screen
CN117115069A (en) Hardening correction method, system and storage medium
DE102022100509A1 (en) X-RAY SYSTEM, X-RAY DEVICE, INFORMATION PROCESSING DEVICE AND INFORMATION PROCESSING METHOD
US4900644A (en) Unsharp masking of diagnostic radiation intensifying screens

Legal Events

Date Code Title Description
A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20060124

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20060327

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20061002

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20091006

Year of fee payment: 3

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20101006

Year of fee payment: 4

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20111006

Year of fee payment: 5

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20121006

Year of fee payment: 6