JP3538244B2 - ラジオグラフィー像の表示において飽和画素を捜し出す方法 - Google Patents
ラジオグラフィー像の表示において飽和画素を捜し出す方法Info
- Publication number
- JP3538244B2 JP3538244B2 JP30999094A JP30999094A JP3538244B2 JP 3538244 B2 JP3538244 B2 JP 3538244B2 JP 30999094 A JP30999094 A JP 30999094A JP 30999094 A JP30999094 A JP 30999094A JP 3538244 B2 JP3538244 B2 JP 3538244B2
- Authority
- JP
- Japan
- Prior art keywords
- value
- image
- signal
- density
- conversion
- 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.)
- Expired - Lifetime
Links
- 229920006395 saturated elastomer Polymers 0.000 title claims description 42
- 238000006243 chemical reaction Methods 0.000 claims description 53
- 238000000034 method Methods 0.000 claims description 53
- 230000005855 radiation Effects 0.000 claims description 29
- 230000007704 transition Effects 0.000 claims description 27
- 238000012545 processing Methods 0.000 claims description 16
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 claims description 12
- 230000003247 decreasing effect Effects 0.000 claims 1
- 238000005286 illumination Methods 0.000 claims 1
- 238000013507 mapping Methods 0.000 description 11
- 230000009466 transformation Effects 0.000 description 10
- 230000008569 process Effects 0.000 description 8
- 230000006870 function Effects 0.000 description 6
- 238000012986 modification Methods 0.000 description 4
- 230000004048 modification Effects 0.000 description 4
- 230000004936 stimulating effect Effects 0.000 description 3
- 230000006399 behavior Effects 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 230000007423 decrease Effects 0.000 description 2
- 238000003745 diagnosis Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 230000004313 glare Effects 0.000 description 2
- 238000002601 radiography Methods 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 241001479434 Agfa Species 0.000 description 1
- 238000005481 NMR spectroscopy Methods 0.000 description 1
- 208000027418 Wounds and injury Diseases 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000001444 catalytic combustion detection Methods 0.000 description 1
- 230000006378 damage Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000007435 diagnostic evaluation Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000012854 evaluation process Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 230000004907 flux Effects 0.000 description 1
- 208000014674 injury Diseases 0.000 description 1
- 239000003550 marker Substances 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000000414 obstructive effect Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000000638 stimulation Effects 0.000 description 1
- 238000003325 tomography Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
- 238000002604 ultrasonography Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/90—Dynamic range modification of images or parts thereof
- G06T5/94—Dynamic range modification of images or parts thereof based on local image properties, e.g. for local contrast enhancement
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10116—X-ray image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20016—Hierarchical, coarse-to-fine, multiscale or multiresolution image processing; Pyramid transform
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
- Conversion Of X-Rays Into Visible Images (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Measuring And Recording Apparatus For Diagnosis (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
- Image Generation (AREA)
- Facsimile Scanning Arrangements (AREA)
- Facsimile Image Signal Circuits (AREA)
- Picture Signal Circuits (AREA)
- Closed-Circuit Television Systems (AREA)
Description
ラフィーの分野におけるものである。より詳しくは、こ
の発明は、ラジオグラフィー像を表示し、且つその表示
された像における飽和画素を探索する方法及び装置に関
する。
いては、広範囲にわたる多様な画像捕捉技術、例えば、
コンピュータ化断層写真法、核磁気共鳴、超音波、CCD
センサ又はビデオカメラによる放射線像の検出、ラジオ
グラフィーフイルムの走査等、が開発されている。これ
らの技術によりラジオグラフィー像のディジタル表示が
得られる。
ば物体より送出されるX 線像は光刺激性リン光体、例え
ば‘92年9月16日付欧州特許公開503702及び
米国Ser.07/842603に記述されているリン光体
の一つ、から成るスクリーンに記憶される。記憶された
放射線像を読み取る技術は、適当な波長を有するレーザ
光線のような刺激を与える放射線でスクリーンを走査
し、刺激によって放射された光を検出し、そしてその放
射光を例えば光電子増倍管により電気表示に変換し、最
後にその信号をディジタル化することから成る。
られたディジタル像は、広範囲にわたる多様な画像処理
技術で処理することができる。
ば、必要とされるだけ何回でも、種々の処理パラメータ
値を用いて行う処理は勿論、オフラインで別の種類の画
像処理技術により処理することもできる。
示装置に印加することができ及び/又はレーザ記録計又
はその類の画像記録計によって再生できる。
で、信号値は、最大及び最小の表示可能な又は再生可能
な濃度値の間で所定の形状を有する信号対濃度マッピン
グ曲線に合った表示又はハードコピーで考えられる濃度
値に変換される。
は、通常、ウインドウ/レベリング過程で処理される。
得られた信号値の全範囲以内で部分範囲が限定される。
この範囲以内では信号値は前記の最小及び最大濃度値の
間の濃度値にマッピングされ、またこの範囲外では、信
号値は、もしそれらが前記部分範囲の最小値より小さけ
れば最小濃度値に、及びもしそれらが前記部分範囲の最
大値より大きければ最大濃度値にマッピングされる。
値の全範囲に対応するその位置とその幅は、例えば特別
の試験方式に関する特定特性として事前に限定されてよ
い。
れる、その範囲の幅と位置は、得られた信号値を例えば
ワークステーションの一部であるモニタにかけ、次いで
表示装置によりオペレータに満足のいく像が得られるま
で前記モニタの輝度及びコントラストを可変して行われ
る視覚制御の下で決定することができる。
た信号範囲に比例して信号対濃度変換曲線がシフトす
る。
断上関連のある領域で像細部が識別される可視像が得ら
れるよう調節される。
対濃度マッピング曲線の形状がラジオグラフィーフイル
ムの特性曲線の形状と類似しているということは、放射
線技師がその特性、より詳しくは従来のラジオグラフィ
ーフイルムの感度測定を経験している故、好都合であ
る。
する補正を信号対濃度転換に組み入れることはさらに有
益である。
イナミックレンジを有する表示装置でのディジタル化処
理のため、1より大きいオリジナル画素値は表示装置に
対しては単一の入力値にマッピングされる。この表示装
置の入力値は、その後、単一濃度値に変換され、その結
果最終的に異なったオリジナル画素値が単一濃度値にマ
ッピングされることになる。異なったオリジナル画素値
から生ずるそのような画素は、以後「飽和画素」と呼
ぶ。
る画素から生ずる濃度値の間の差異は何等識別されず、
その結果貴重な情報が失われてしまうかも知れない。
イナミックレンジに限界があるため、飽和された像領域
をスクリーン上で識別することは、不飽和像又は画素に
比して困難なので、飽和画素は「飽和」されていると認
識されないかも知れない。
ー像を評価する時、放射線技師は、関連画素が飽和状態
になっており且つ事実上これらの画素の有効濃度は表示
された像に見られる濃度とは異なっているという事実に
気付かないかも知れない。
て)写真フィルムのようなハードコピー再生媒体上に飽
和領域が存在する像のハードコピーを作る時、表示装置
のダイナミックレンジに比してハードコピーのそれが大
きいため、その飽和領域は見えるようになる。飽和領域
は像解釈中はむしろその妨げをする。
ル信号表示の形で得られ且つその放射線像が表示装置上
に可視化される系において飽和されている画素の位置を
決定する方法を提供することにある。
射線像の診断に役立つ評価プロセスにとって重要な画素
が飽和されないような方法で像を表示することにある。
ンの表示モニタに応用できるよう前述の方法を提供する
ことにある。
ン光体スクリーンに記憶され且つ前記スクリーンが刺激
を加える照射で走査され、照射で放出された光が検出さ
れてディジタル表示に変換される系に応用するために前
述の方法を提供することにある。
ろう。
前記像のディジタル画素値が下記の処置から成る第1の
信号対濃度転換に従って濃度値に転換される時に飽和さ
れる画素の位置を決定することによって達成される: − 前記放射線像のディジタル画素値を、前記第1の転
換に従って変換する時飽和されるであろう前記像の画素
に対応する画素値に少なくとも一つの不連続を含んでい
ること以外は第1の転換と実質上同一の第2の信号対濃
度転換に従って濃度値に転換することと; − 前記第2の転換に従って変換された画素値を表示装
置に印加すること。
は、表示装置の入力時の個別化処理中表示装置に対する
単一の入力値に変換される異なった画素値から生ずる画
素のことである。これらの異なった画素値はその後さら
に単一の濃度値に変換される。
その転換に従って変換される時飽和画素に帰着するであ
ろう画素値を見つけ出すことができる。しかし、像上の
飽和画素の位置は簡単な転換の評価では見つけ出すこと
はできない。この発明の方法処置に従って像を表示する
ことにより飽和画素の位置が目に見えるようになる。
は、一般に、信号値を表示装置に対する別々の入力値に
マップする転換と、表示装置の入力値を濃度値に変換す
る転換との組み合わせで構成され、後者の転換はしばし
ば付加的に補正可能な非線形転換である。
の横座標値が転換されるであろう)縦座標値に第1の転
換によって転換されるであろう横座標値(複数の場合あ
り)に不連続(複数の場合あり)を有している前記第2
の転換に従って転換することにより、(表示に相応する
画素は「飽和画素」である)、及び前記第2の転換に従
って転換される画素を表示することにより、飽和画素の
位置は表示中で識別できるようになる。
い。従って、1より多い不連続が個別画素に存在してい
てよく、さらに従って画素値の範囲にわたって転換が不
連続であってよい。
の縦座標値は次のように設定できる:連続的に進展する
転換内で不連続は存在でき、そこでは横座標値は最大表
示可能濃度値(又は強度値)にマッピングされ、従って
飽和画素は表示装置上では黒い画素として現れる。また
その上、変換関数が不連続故、横座標値は最小表示可能
値と等しい縦座標値に転換され、従って飽和画素は表示
では白い画素として表されることになる。
何れにしろ不連続のある縦座標値がこの不連続に隣合う
画素の縦座標値から大きく異なることは最も好都合であ
り、その結果、第2の変換転換に従って像を変換する
時、飽和画素はそれら自体の濃度とその周辺濃度との間
に存在する大きな差異により非常に良好に目に見えるよ
うになる。
射線像のディジタル表示は、既に像が記憶されている光
刺激性リン光体スクリーンを刺激を加える照射で走査
し、刺激によって放出された光を検出し、そして検出さ
れた光をディジタル信号表示に変換することにより得ら
れるものである。
1の信号対濃度転換は、放射線技師が経験している写真
フィルムの特性曲線に類似している転換曲線と、表示装
置の非線形出力特性と、さらにこの非線形特性を補正す
る曲線との組合わせとして得られる形状を持っている。
により表すことができ、それにより所定の最小値より小
さい信号値は最小表示可能濃度値に等しい値に転換さ
れ、所定の最大値を越える信号値は最大表示可能濃度値
に等しい値に転換される。
ミュレーションである要素の構成法は、欧州特開549
009に記述されている。転換がS形を有する範囲の限
界である予定値は、表示に向いており且つさらに診断に
適切な信号範囲の限界でもある。
の不連続があり且つ前記最大値に第2の不連続があるこ
とを除けば、実質上第1の転換と同一の形状をむしろ有
している。
り小さい画素値の縦座標値は、好ましいことには表示可
能濃度値の最大値と等しく、前記第2の不連続より大き
い画素値の縦座標値は、好ましいことには表示可能濃度
値の最小値と等しく、従って縦座標値は周囲の縦座標値
とは大きく異なり、結果として第1の転換が適用される
とき飽和される画素の位置が最適条件で可視化される。
上で検出されれば、放射線技師又はオペレータは、配置
画素が飽和されるという事実によって診断上の評価プロ
セスが妨げられているかどうかを検討することができ
る。
て結果的にこれらの位置での正確な濃度値を識別するこ
とが不可能になり診断上関係ある何らかの情報が失われ
てしまったかどうか、を判断することが彼にとって可能
となろう。この評価は実施される試験の種類によって及
び予想される診断によって左右される。
に関係あるような場合、オペレータ又は放射線技師は、
表示に表れる飽和の原因である第1の転換は不十分であ
ると判断してよく、また、診断上関係する位置での飽和
を避けるよう第1の転換曲線の変更を決めてもよい。
の転換に相当する表検索での値を変更することにより履
行される。
い画素値が単一の縦座標値にマッピングされ且つ個々の
画素値間の変化が検出できるように、前記第1の変換転
換の勾配を増すことにより行うことができる。
となり、そこでは第1の転換に従う画素値の変換により
得られた像が表示されよう。画素値はその後、飽和画素
の横座標値に少なくとも一つの不連続がある第2の転換
で処理される。第2の転換に従って転換された画素値の
表示に生ずる像及び特に飽和画素の位置が検討される。
もし必要なら、第1の変換曲線は、診断上関係ある領域
の画素が飽和しないように修正される。次に、修正され
た「第2の変換曲線」が定められ、プロセスは再開され
る。この手順は、ディスプレー上に満足な像が得られる
まで繰り返される。最終的に、例えばフィルム上にハー
ドコピー像が得られるよう、表示された像をハードコピ
ー記録計に入力することができる。
についても実行可能であるということ及び従ってこの第
2の転換はラジオグラフィー像の画素値を表示可能濃度
値に転換するために使用されるということは、熟練した
当業者にとって明かとなろう。
ワークステーションで実行される: − 第1の信号対濃度転換を記憶する手段; − 前記第1の転換に従い変換されるとき飽和される
(単数・複数の)画素値に少なくとも一つの不連続があ
ることを除けば、実質上前記第1の転換と同一の第2の
信号対濃度転換を記憶する手段; − これらの転換の一つに従いディジタル像表示を転換
する手段; − 転換された像を表示する手段。
連続(複数の場合あり)の位置を前記ディジタル表示を
持つ信号値の範囲内で変え、且つ前記第2の転換が不連
続となる範囲幅を変えるために使われるユーザ用インタ
フェースを含んでよい。
さらに詳しい説明は、以下に下記図面を引用して述べ
る:
的に示す。
示せず)より送出されるX線に光刺激性リン光体スクリ
ーン(3) を露光させることにより(2) 前記スクリーンに
記憶された。その光刺激性蛍光スクリーンは、電気的消
去書込み可能読出し専用メモリ(EEPROM)(5) を組込んだ
カセット(4) に収めた。識別ステーション(6) では、種
々のデータ、例えば、患者識別データ(姓名、生年月
日)及び露光及び/又は処理パラメータのような信号処
理に関するデータがEEPROM(5) に書込まれた。
れた情報及び光刺激性リン光体スクリーンに記憶された
像が読取られた。読取り法は図2に説明する。記憶され
た像は、レーザ16で発射された刺激線でリン光体スク
リーンを走査して読取られた。刺激線は電流磁気偏向1
7によって主走査方向に偏向された。副走査は、リン光
体スクリーンを副走査方向18に搬送して実行された。
刺激された放射は、電気的像表示に変換するために光集
光器19によって光電子増倍管20の方へ向けられた。
次に、信号は抽出・保持回路21によって抽出され、対
数値log E (E露光値を表示)に変換されそしてアナログ
/ディジタル変換器22によって量子化された。ディジ
タル像信号23は、読取り装置(図1、図7)の画像処
理モジュールに送られそこで内部バッファに記憶され
た。
ため下見モニタ8にも入力され、これによりオペレータ
は実行された露光についての早期フィードバックが可能
となった。
らにワークステーション11及び関連した検閲コンソー
ル12に接続され、そこでオフライン処理が実施され
た。
クステーション及び関連検閲コンソールはバッファ9を
介して外部記録計10に接続された。
装置の処理ユニットでのオンラインか又はワークステー
ションでのオフラインで実行された画像処理は、この特
定の実施例においては多分解能表示で、より詳しくは画
像のピラミッド型多分解能表示で実行された。
ベルでの詳細像の連続したものと前記多分解能レベルの
最小値より低い分解能レベルでの残りの像とに分解する
ことにより得られた。各詳細像の画素数は各々の粗分解
能レベルでは減少する。
像は、それぞれ下記処理のK 回繰返しの結果として得ら
れた: a) 低パスフィルタを現在の繰返しに対応する概略像に
適用して次位の粗さのレベルでの概略像を計算し、どん
な方法にしろ、元の像を最初の繰返しの粗さで前記低パ
スフィルタへの入力として使いながら、その結果を空間
周波数帯域幅の低減に比例して準抽出すること; b) 詳細像を、現在の繰返しに対応する概略像と方法処
理サブa)に従って計算された次位の粗さの分解能レベル
での概略像(両像は後の方の像の適当な書き込みによっ
てレジスタに入れられる)との間の画素方式の差として
計算すること。従って残りの像は最後の繰返しによって
作られた概略像と等しい。
並びに他の処置例は欧州特開EP527525に記述さ
れている。
細像の画素値を与えるよう修正処置された。詳細像の画
素値は、引数値の増加につれ徐々に減少する勾配を持っ
た少なくとも一つの非線形単調増加奇数変換関数に従い
修正された。前述の変換関数の例もまた既に述べた欧州
特開527525に開示されている。
アルゴリズムを残りの像及び修正詳細像に適用すること
により計算された;ここで再構成アルゴリズムは、もし
それが残りの像に及び修正なしで詳細像に適用されれ
ば、前記元の像又はそれに接近した近似像が得られるよ
うなものである。
ルゴリズムは、最も粗い詳細像と残りの像からスタート
する次の処理をK 回繰り返すことにより計算された: − 画素に関して同一分解能レベルにある詳細像を先の
繰返しに対応するより粗い分解能レベルの近似像に付加
することにより、近似像を現在の分解能レベルで計算す
ること;ここで両像は、どんな方法にしろ、最初の繰返
しの粗さにある前記のより粗い近似像の代わりに残りの
像を使いながら、その後の像の適当な繰り返しによりレ
ジスタに入れられる。
527 525に開示されている。
を通さない「コリメーション」物質を使って彼らの患者
をX線に対する不必要な露出から防御している。その物
質は、診断上重要視されない患者のそうした部位を覆う
ようにX線束の通路に置かれる。この技法は患者の線量
を減らす上、さらに像の対象視野における散乱照射量を
制限するという利点がある。X線を通さない物質の投影
によって生ずる領域(シャドウ領域)は散乱照射によっ
てのみ露光される。
フイルム上で又は表示装置上で放射像を表示する際問題
を生ずることがある。そこではシャドウ領域は相対的に
明るく、そして無修正で表示されるとき、特に非露光領
域が比較的大きいとき、眩惑によって鋭敏な傷害の診断
をし損じるかも知れない。
定し、次いで、診断上関係あるゾーン外の画素が可視化
され、従って制御卓のスクリーン上でそのハードコピー
が観察されるか又は表示される時、前記像の部分より送
出される光が効果的に遮蔽されるような方法で、放射線
像の信号値を表示装置上のソフトコピーか又はフイルム
上のハードコピーとして表示される濃度値に変換するこ
とが提案された。
射線像内の診断上無関係な像のゾーン中で構成される画
素の電気信号は、0.5と2.5の間の濃度値に変換さ
れる。
の範囲内で診断上無意味な情報はもはや目に見えない。
れる光によりもたらされる眩惑の影響が低減されるとい
う事実に満足するが、一方で診断上無意味な部位にあっ
た情報の少なくとも痕跡を保持できればと望む。
放射線像の画素値を、診断上関係あるゾーン外の画素が
可視化又は再生され、従って制御卓のスクリーン上でそ
のハードコピーが観察されるか又は表示される時、前記
像の部分より送出される光が軽減される(それ故、無関
係ゾーンの平均輝度レベルが低減される)ような方法
で、変換することにより解決された。そうすることによ
り診断上無意味な領域の放射線像の部分についての解剖
学的情報は可視状態のまま残されることになる。
断上無意味なゾーンにある時より対象領域にある時に、
別々に処理することにより得られた。
内で別個の画素処理を可能にするためには、先ず最初に
どの画素が対象領域に属しているかを決めなければなら
ない。
が開発されている。欧州特開523771では、放射線
像の診断上関係あるゾーンを手動で詳細に描写すること
について言及された。
領域を自動的に決定する方法である。この方法は、欧州
特開610605の拡張部に記述された。
であるか否かについての情報がオーバレイ像によって得
られた。このオーバレイ像は、画素が対象領域に属して
いるか否かを識別する標識から成る低分解能2進画像で
ある。
らに表検索にかけられ、(1) 対象領域外の平均輝度及び
コントラストが低減され、且つ(2) この領域における未
処理像の情報が保持され且つ可視化されることを規定し
た変換を表示する。
は、識別転換によって変換され、それに反し、診断上関
係ある領域外の画素は、識別転換より上にある直線で表
すことができ且つ0と1の間の勾配を持つ変換関数に従
って変換され、その結果、診断上関係ある領域外の画素
は像の残りの部分のそれより高い平均濃度に変換され且
つこの像部分の像情報は失われない。
+(1-a)fmax で表され、ここでg(x,y ) は、転換後の画
素値であり、f(x,y)は変換前の画素値及び"a" は対象領
域、fmaxはgmaxに等しい。勾配"a" はゼロと1 の間の値
を持ち且つ診断上無関係なゾーンの画素値が変換され得
る最小濃度を決定する。
に分解され、次いで修正され最後に再生処理されるこの
特有の実施例において、関係ある像領域外の画素をさら
に変換するこの(表検索のように履行される)処理は、
画素は対象領域の部分であるか否かの情報を考慮しなが
ら、部分的に再生された像の画素に対して適用された。
後の情報は上述のオーバレイ像から利用できる情報が適
用される。
次いで全体的に再構成され且つ処理された像が得られる
よう欧州特開527525に記述された再生処理の残り
の工程にかけられた。
表示の前に階調(ぼかし)処理が施された。
使って得られたオーバレイ像の特殊性により、即ち、オ
ーバレイ像は低分解能2 進画像であるという事実によっ
て引起こされた障害となる人為的影響は、低パスフィル
タを応用して2進オーバレイを多重値オーバレイ像に転
換することにより避けられた。採用されたマッピング転
換の漸進的移行は、対象領域外の画素に対して講じられ
た。採用されたマッピング転換の勾配と交差は、2 進で
はなく多重値であるオーバレイ像における画素値により
制御された。採用されたマッピング転換は、コリメーシ
ョン物体のシャドウ領域で最大勾配を持ち、診断上関係
ある像領域での識別マッピングに等しい。この転換は数
式で次のように表すことができる: ここでf(x,y)は転換前の画素値であり、c(x,y)は、対象
領域内の画素に関してはゼロに等しく、コリメーション
シャドウ領域内部の画素に関しては63に等しく、"a" は
コリメーションシャドウ領域内のマッピング勾配を表
す。 − もし"a" が1 に等しければ、g(x,y)は像のどこにお
いてもf(x,y)と等しく、換言すれば、対象領域内又はコ
リメーションシャドウゾーン内の画素の間に何等差異は
存在しない。 − もし例えば"a" が1/3に等しければ;対象領域で
c(x,y)はゼロに、g(x,y)はf(x,y)に且つ対象領域外でc
(x,y)は63に等しく、故に g(x,y) =af(x,y) +(1-a)fmax となる。
り、c はゼロと63の間の値をとり従って漸進的移行が
両マッピング間で得られる。 − もし"a" がゼロに等しければ、g(x,y)は、コリメー
ションシャドウ領域においてfmaxと等しい。
なされた特別な処置はワークステーションと同様にオン
ラインで行える。
の表示の中に飽和画素を捜し出す方法の説明である。
によって得られたオリジナル像の信号値で占められてい
る像の範囲が、元の信号値より小さい表示装置のダイナ
ミックレンジにマッピングされた。次に、マッピングさ
れた信号値は、使用された表示装置(CRT 装置)に特有
の非線形転換曲線に従って、対応する濃度値に転換され
た。CRT 装置の非線形挙動に対する補正も適用した。図
3は、一般的に単調増加パスを有する信号対濃度転換の
説明図である。Pminより小さい画素値を有する画素は、
最小濃度値Dminにマッピングされるが、Pmaxより大きい
画素値を有する画素は、最大濃度値Dmaxにマッピングさ
れる。露光された光刺激性リン光体スクリーンから読取
られた像が上記転換に従って変換される時、及びその変
換像が表示用モニタにかけられる時、Pminより小さいか
又はPmaxより大きい画素値を有する画素は飽和される。
色調から成る)低濃度領域内又は(極めて暗い灰色の色
調から成る)非常に暗い領域内に位置しており、従って
表示可能範囲の限界にある濃度値を有する周囲領域にお
いては、飽和画素は不飽和画素と識別できなかった。
上の表示像でどの場所が飽和画素で占められそしてどこ
がそうでないかを見分けることは不可能であった。放射
線技師又は技術者にとって、いくつかの画素が表示に際
し飽和されたという事実は、これらの画素が診断上無関
係な像領域に位置していたという理由で容認できるかど
うか、あるいは又、診断上重要な情報が失われたかどう
か、を検討することが必然的に不可能であった。
発明の方法に従い−図4に示された第2の転換曲線に応
じてもう一度濃度値にマッピングされた。
された。より詳細には、この第2の転換曲線は、−図3
の転換に従って転換される時、それ故これらの不連続の
無い状態で−飽和が起こるであろう場所に位置した画素
に対応する横座標値に不連続が存在することを除けば、
図3の曲線の経路と同一であった。図4は、図3と同じ
曲線を示しているが、低信号レベル及び高信号レベルの
位置で濃度転換曲線に対する信号のパスは修正され、そ
の結果小さい信号値はモニタ上で表示し得る最大濃度値
に及び最高信号値は最小濃度にマッピングされることに
なる。
応する濃度値に変換することによって得られる像が表示
される。
とき飽和された画素の位置は、今やこれらの画素は周囲
の濃度とははるかに異なった濃度を持つため、目でみる
ことができた。事実、単調増加パスに従って画素が濃度
値にマッピングされたが、最小信号値より下の画素値は
最大濃度レベルに不連続に変換され、それ故、表示では
信号軸の低信号レベルの部分に位置する画素は最大表示
可能濃度にマッピングされ、そして低濃度の周囲領域内
で識別できるようになる。
号値は表示された像に使える最小濃度レベルに不連続的
にマッピングされ、従ってそれらは全般的に暗い灰色の
周囲領域内で「白色の」画素として現れる。
囲の限界に位置しない値の範囲(ウインドウとも呼ばれ
る)全域で不連続となった。
より、この不連続範囲のレベルとウインドウを、例えば
不連続範囲を飽和が起こる画素値の範囲で移動させるよ
うなやり方で変えることも可能である。この転換を適用
するとき、飽和画素の位置は表示装置上で可視化するこ
とができる。
分にあるか否かが放射線技師によって評価されれば、放
射線技師は第1の信号対濃度転換曲線に相当する表検索
の値の少なくともいくつかを修正でき、その結果診断上
重要な像領域においてもはや飽和は起こらないであろ
う。
濃度値に変換するということについてなされたものであ
る。熟練した当業者にとって、濃度値の代わりに強度値
が使われる場合に同じ方法が適用し得ることは明かであ
ろう。
である。
読取る装置の特定の実施例である。
Claims (13)
- 【請求項1】 第1の信号対濃度転換に従い放射線像の
ディジタル画素値が濃度値に転換されるとき飽和される
前記像における画素の位置決定方法において、 − 前記放射線像のディジタル画素値を前記第1 の転換
に従って変換されるとき飽和されるであろう前記像の画
素に対応する画素値に少なくとも一つの不連続を含んで
いることを除けば前記第1の転換と実質上同一である第
2の信号対濃度転換に従って濃度値に変換することと; − 前記第2の転換に従って変換された画素値を表示装
置に用いることとを特徴とする方法。 - 【請求項2】 前記第2の転換の不連続(複数の場合あ
り)の縦座標値と前記不連続に隣合う値の縦座標値との
間の差が前記表示装置のダイナミックレンジの範囲内で
可能な限り大きいことを特徴とする請求項1記載の方
法。 - 【請求項3】 不連続(複数の場合あり)の縦座標値が
前記表示装置のダイナミクレンジの最大値である請求項
1記載の方法。 - 【請求項4】 不連続(複数の場合あり)の縦座標値が
前記表示装置のダイナミックレンジの最小値である請求
項1記載の方法。 - 【請求項5】 前記転換が不連続挙動を有する横座標値
の範囲に対応する縦座標値が前記表示装置のダイナミッ
クレンジの最大値に等しい請求項1記載の方法。 - 【請求項6】 前記転換が不連続挙動を有する横座標値
の範囲に対応する縦座標値が前記表示装置のダイナミッ
クレンジの最小値に等しい請求項1記載の方法。 - 【請求項7】 前記ディジタル信号表示が、像を記憶し
ている光刺激性リン光体スクリーンを刺激照射で走査
し、刺激で放出された光を検出し且つ検出された光をデ
ィジタル信号表示に変換することにより得られる請求項
1乃至6のいずれ一に記載の方法。 - 【請求項8】 − 前記第1の転換は信号値を濃度値
に変換する単調増加又は減少関数によって表すことがで
き、それにより所定の最小値より小さい信号値は最小の
表示可能濃度値に等しい値に転換され且つ所定の最大値
を越える信号値は最大の表示可能濃度値に等しい値に転
換されることと− 前記第2の転換は前記最小値に第1
の不連続があり、前記最大値に第2の不連続があること
を除けば前記第1の転換と実質上同一の形状を有するこ
ととを特徴とする請求項1記載の方法。 - 【請求項9】 前記第1の不連続より小さい全ての画素
値の縦座標値は表示可能濃度値の最大値に等しく、且つ
前記第2の不連続より大きい全ての画素の縦座標値は表
示可能濃度値の最小値に等しい請求項8記載の方法。 - 【請求項10】 前記第2の転換に従って変換された信
号表示を前記表示装置に印加することによって得られた
表示像の不連続(複数の場合あり)の位置検討に応じて
前記第1の転換の形状が変化させられることを特徴とす
る請求項1乃至9のいずれか一記載の方法。 - 【請求項11】 − 第1の信号対濃度転換を記憶する
ための手段と; − 前記第1の転換に従って変換されるとき飽和される
画素値(複数の場合あり)に少なくとも一つの不連続が
あることを除けば前記第1の転換と実質上同一である第
2の信号対濃度転換を記憶する手段と; − これらの転換の一つに従ってディジタル像表示を転
換する手段と; − 転換された像を表示する手段とを含むことを特徴と
するディジタル像表示を処理し、且つ処理されたディジ
タル像表示に対応する像を表示するためのワークステー
ション。 - 【請求項12】 前記ディジタル表示の信号値の範囲内
で前記不連続(複数の場合あり)の位置を変えるために
採用されたユーザ用インターフェースを包含する請求項
11記載のワークステーション。 - 【請求項13】 前記第2の転換が不連続となる範囲の
幅を変化させる手段を包含する請求項11記載のワーク
ステーション。
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
EP93203277A EP0654761B1 (en) | 1993-11-23 | 1993-11-23 | Method and apparatus of locating saturated pixels in the display of a radiographic image |
DE93203277.4 | 1993-11-23 |
Publications (2)
Publication Number | Publication Date |
---|---|
JPH07222142A JPH07222142A (ja) | 1995-08-18 |
JP3538244B2 true JP3538244B2 (ja) | 2004-06-14 |
Family
ID=8214182
Family Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP30999094A Expired - Lifetime JP3538244B2 (ja) | 1993-11-23 | 1994-11-18 | ラジオグラフィー像の表示において飽和画素を捜し出す方法 |
JP31268694A Expired - Fee Related JP3662283B2 (ja) | 1993-11-23 | 1994-11-21 | ラジオグラフィー像の表示において診断上無関係な領域の可視化 |
Family Applications After (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP31268694A Expired - Fee Related JP3662283B2 (ja) | 1993-11-23 | 1994-11-21 | ラジオグラフィー像の表示において診断上無関係な領域の可視化 |
Country Status (4)
Country | Link |
---|---|
US (2) | US6249596B1 (ja) |
EP (1) | EP0654761B1 (ja) |
JP (2) | JP3538244B2 (ja) |
DE (2) | DE69308024T2 (ja) |
Families Citing this family (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS63177791A (ja) * | 1987-01-16 | 1988-07-21 | Kibun Kk | 固定化酵素又は固定化微生物の製造法 |
US5694490A (en) * | 1995-11-27 | 1997-12-02 | Sun Microsystems, Inc. | System and method for a simultaneous multi-band block-stop filter |
JP4083251B2 (ja) * | 1996-08-16 | 2008-04-30 | 富士フイルム株式会社 | 照射野外画像処理方法および装置 |
US6271536B1 (en) | 1997-10-08 | 2001-08-07 | Agfa-Gevaert | Radiographic image identification method |
FR2786293B1 (fr) * | 1998-11-24 | 2001-04-27 | Ge Medical Syst Sa | Procede de compensation de l'epaisseur d'un organe |
FR2796740B1 (fr) | 1999-07-19 | 2001-10-26 | Ge Medical Syst Sa | Procede et systeme de gestion de la saturation sur une image radiographique numerisee |
JP4005277B2 (ja) | 1999-09-06 | 2007-11-07 | 富士フイルム株式会社 | 画像処理装置、方法及び記録媒体 |
FR2803069B1 (fr) * | 1999-12-28 | 2002-12-13 | Ge Medical Syst Sa | Procede et systeme de compensation de l'epaisseur d'un organe |
US7072072B1 (en) * | 2000-05-02 | 2006-07-04 | Xerox Corporation | Color rendering optimized for text and line art |
EP1239415A1 (en) * | 2001-03-05 | 2002-09-11 | Ge Medical Systems Sa | Method and system of management of the dynamics of a digitized radiological image |
US7397964B2 (en) * | 2004-06-24 | 2008-07-08 | Apple Inc. | Gaussian blur approximation suitable for GPU |
US20060050983A1 (en) * | 2004-09-08 | 2006-03-09 | Everest Vit, Inc. | Method and apparatus for enhancing the contrast and clarity of an image captured by a remote viewing device |
JP5259175B2 (ja) * | 2007-12-26 | 2013-08-07 | 株式会社東芝 | 超音波診断装置、超音波画像処理装置及び超音波画像処理プログラム |
US7896229B2 (en) * | 2008-02-22 | 2011-03-01 | Apteryx, Inc. | Auto-distribution of scanned digital images based on standardized identifiers |
US20090271329A1 (en) * | 2008-04-25 | 2009-10-29 | Akiyoshi Mohri | Utility management system |
US8265369B2 (en) * | 2009-04-16 | 2012-09-11 | Apteryx, Inc. | Apparatus and method for virtual flaw removal from X-ray sensitive plates |
CN112635875B (zh) * | 2020-03-27 | 2023-05-26 | 苏州清陶新能源科技有限公司 | 一种电池用软壳、其制备方法和用途 |
CN114339447B (zh) | 2020-09-29 | 2023-03-21 | 北京字跳网络技术有限公司 | 图片转视频的方法、装置、设备及存储介质 |
Family Cites Families (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS56104645A (en) * | 1979-12-25 | 1981-08-20 | Fuji Photo Film Co Ltd | Radiation picture treating method and its device |
US4542459A (en) * | 1982-11-26 | 1985-09-17 | General Electric Company | Matched filter for x-ray hybrid subtraction |
JPS60156055A (ja) * | 1984-01-26 | 1985-08-16 | Fuji Photo Film Co Ltd | 放射線画像情報読取条件決定方法 |
US4803639A (en) * | 1986-02-25 | 1989-02-07 | General Electric Company | X-ray inspection system |
EP0352491B1 (en) * | 1988-06-30 | 1996-08-14 | Dainippon Screen Mfg. Co., Ltd. | Method of generating gradation correction curve for correcting gradation character of image |
EP0445450A1 (en) * | 1990-03-07 | 1991-09-11 | International Business Machines Corporation | Image processor mepping source pixel intensities to a limited range of display intensity values |
DE69214229T2 (de) * | 1991-08-14 | 1997-04-30 | Agfa Gevaert Nv | Verfahren und Vorrichtung zur Kontrastverbesserung von Bildern |
DE69217601T2 (de) * | 1991-09-27 | 1997-09-18 | Agfa Gevaert Nv | Verfahren zum Einsetzen von spezifischen Transformationskurven in Thermosublimationsdruckern |
DE69331719T2 (de) * | 1992-06-19 | 2002-10-24 | Agfa-Gevaert, Mortsel | Verfahren und Vorrichtung zur Geräuschunterdrückung |
EP0610603B1 (en) * | 1993-02-11 | 1999-09-08 | Agfa-Gevaert N.V. | Fast interactive off-line processing method for radiographic images |
US5471987A (en) * | 1993-03-30 | 1995-12-05 | Konica Corporation | Method of compressing a dynamic range for a radiation image |
US5447153A (en) * | 1993-07-02 | 1995-09-05 | Eastman Kodak Company | Real-time window/leveling on a radiographic workstation |
EP0639819B1 (en) * | 1993-08-16 | 1997-01-22 | Agfa-Gevaert N.V. | Method and apparatus for controlling the effectively applied irradiation amount |
-
1993
- 1993-11-23 DE DE69308024T patent/DE69308024T2/de not_active Expired - Lifetime
- 1993-11-23 EP EP93203277A patent/EP0654761B1/en not_active Expired - Lifetime
-
1994
- 1994-11-02 DE DE69420516T patent/DE69420516T2/de not_active Expired - Lifetime
- 1994-11-08 US US08/335,915 patent/US6249596B1/en not_active Expired - Lifetime
- 1994-11-08 US US08/335,917 patent/US6356652B1/en not_active Expired - Lifetime
- 1994-11-18 JP JP30999094A patent/JP3538244B2/ja not_active Expired - Lifetime
- 1994-11-21 JP JP31268694A patent/JP3662283B2/ja not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
EP0654761B1 (en) | 1997-02-05 |
EP0654761A1 (en) | 1995-05-24 |
DE69420516D1 (de) | 1999-10-14 |
JPH07226882A (ja) | 1995-08-22 |
JPH07222142A (ja) | 1995-08-18 |
US6249596B1 (en) | 2001-06-19 |
JP3662283B2 (ja) | 2005-06-22 |
DE69308024D1 (de) | 1997-03-20 |
DE69308024T2 (de) | 1997-08-14 |
US6356652B1 (en) | 2002-03-12 |
DE69420516T2 (de) | 2000-05-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP3538244B2 (ja) | ラジオグラフィー像の表示において飽和画素を捜し出す方法 | |
US4747052A (en) | Radiation image processing | |
US4317179A (en) | Method and apparatus for processing a radiographic image | |
JP3188491B2 (ja) | X線記録のダイナミック圧縮方法及びその装置 | |
EP0466907B1 (en) | Tone-scale generation method and apparatus for digital x-ray images | |
EP0032237B1 (en) | Method of and apparatus for processing a radiographic image | |
US8009892B2 (en) | X-ray image processing system | |
EP0523771B1 (en) | Processing method in radiation image recording systems | |
US6370265B1 (en) | Method for generating gray scale transfer functions for use in displaying a digital radiogram | |
US5978443A (en) | Automated removal of background regions from radiographic images | |
JP3290483B2 (ja) | 統計的画素サンプリングによってデジタル影像を処理してヒストグラムを作製する方法 | |
US4802093A (en) | X-ray image-processing apparatus utilizing grayscale transformation | |
EP0490532B1 (en) | Radiographic image processing apparatus and method | |
JP2005109790A (ja) | 医用画像処理装置 | |
JPH0792588A (ja) | 光刺激性蛍光スクリーンに記録された放射線画像を読み取るための装置の感度をモニターする方法 | |
JP4083251B2 (ja) | 照射野外画像処理方法および装置 | |
US6344858B1 (en) | Method of evaluating image processing performed on a radiographic image | |
JP2681482B2 (ja) | 放射線画像処理装置 | |
EP0654762B1 (en) | Visualisation of diagnostically irrelevant zones in a radiographic image | |
EP0527524B1 (en) | Image quality enhancing method | |
JP3197952B2 (ja) | オリジナル画像信号によって表現される画像の増大した可視画像を生成する方法 | |
JPH0572624B2 (ja) | ||
JP3181125B2 (ja) | 再生に関連がある信号レンジの決定方法 | |
JPH0516859B2 (ja) | ||
JP2542183B2 (ja) | X線画像処理方法および装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
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: 20040312 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20040319 |
|
R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313113 |
|
R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090326 Year of fee payment: 5 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100326 Year of fee payment: 6 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100326 Year of fee payment: 6 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110326 Year of fee payment: 7 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120326 Year of fee payment: 8 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120326 Year of fee payment: 8 |
|
S531 | Written request for registration of change of domicile |
Free format text: JAPANESE INTERMEDIATE CODE: R313531 |
|
S533 | Written request for registration of change of name |
Free format text: JAPANESE INTERMEDIATE CODE: R313533 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120326 Year of fee payment: 8 |
|
R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130326 Year of fee payment: 9 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130326 Year of fee payment: 9 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140326 Year of fee payment: 10 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
EXPY | Cancellation because of completion of term |