JPH0737080A - Ｘ線画像の歪み補正方法及び装置 - Google Patents

Ｘ線画像の歪み補正方法及び装置

Info

Publication number
JPH0737080A
JPH0737080A JP6104060A JP10406094A JPH0737080A JP H0737080 A JPH0737080 A JP H0737080A JP 6104060 A JP6104060 A JP 6104060A JP 10406094 A JP10406094 A JP 10406094A JP H0737080 A JPH0737080 A JP H0737080A
Authority
JP
Japan
Prior art keywords
image
pixels
input image
value
output image
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
JP6104060A
Other languages
English (en)
Other versions
JP3672588B2 (ja
Inventor
Reiner Heinrich Dr Koppe
ハインリヒ コッペ ライナー
Erhard Paul Artur Klotz
パウル アルター クロッツ エアハルト
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.)
Koninklijke Philips NV
Original Assignee
Philips Electronics NV
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 Philips Electronics NV filed Critical Philips Electronics NV
Publication of JPH0737080A publication Critical patent/JPH0737080A/ja
Application granted granted Critical
Publication of JP3672588B2 publication Critical patent/JP3672588B2/ja
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Classifications

    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N19/00—Methods or arrangements for coding, decoding, compressing or decompressing digital video signals
    • H04N19/50—Methods or arrangements for coding, decoding, compressing or decompressing digital video signals using predictive coding
    • H04N19/503—Methods or arrangements for coding, decoding, compressing or decompressing digital video signals using predictive coding involving temporal prediction
    • H04N19/51—Motion estimation or motion compensation
    • H04N19/537—Motion estimation other than block-based
    • H04N19/54—Motion estimation other than block-based using feature points or meshes
    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T3/00—Geometric image transformations in the plane of the image
    • G—PHYSICS
    • G06—COMPUTING OR CALCULATING; COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T3/00—Geometric image transformations in the plane of the image
    • G06T3/18—Image warping, e.g. rearranging pixels individually

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Multimedia (AREA)
  • Signal Processing (AREA)
  • Image Processing (AREA)
  • Apparatus For Radiation Diagnosis (AREA)
  • Image Analysis (AREA)

Abstract

(57)【要約】 （修正有） 【目的】 Ｘ線照射により生成され、位置及び各々の入
力画像値により定義された画素により構成された歪んだ
入力画像を位置及び各々の出力画像値により定義された
画素により構成された実質的に歪みのない出力画像に変
換する。 【構成】 入力画像の４つの角の点は幾つかの画素に囲
まれた出力画像内の四角形（例えばＰ１．．．Ｐ４）を
定義する。これらの画素の各々に対して出力画像値は、
関連した入力画像画素の入力画像値から得られ、該出力
画像値は四角形の面領域全体に対する四角形により囲ま
れた画素の面領域の比に対応する。出力画像（入力画
像）の画素の角の点の位置は、入力画像（出力画像）の
角の点によって決定された多角形の面領域で決定され、
出力画像の画素内に位置する角の点の入力画像値が出力
画像の関連する出力の出力画像値を決定するために重み
付けされた和により用いられる。

Description

【発明の詳細な説明】
【０００１】
【産業上の利用分野】本発明はＸ線照射により生成され
位置及び各々の入力画像値により定義された画素により
構成された歪んだ入力画像を、位置及び各々の出力画像
値により定義された画素により構成された実質的に歪み
のない出力画像に変換する方法に係わる。
【０００２】
【従来の技術】画素は、それの環境から得られた均一な
輝度値を割り当てられうる面領域の最小単位を意味する
と理解される。Ｘ線照射中のＸ線の吸収の測定値である
この輝度値を、以下に画像値と称する。歪んだ入力画像
の画素に関連した画像値を「入力画像値」と称し、一方
出力画像の画像値を「出力画像値」と称する。入力画像
値又は入力画像はＸ線照射により生成される。出力画像
値は変換処理により入力画像値から得られる。
【０００３】この種の方法は、例えばＴｅｈｒａｎｉ等
による（ＩＥＥＥ Ｃｏｍｐ．ｉｎＣａｒｄｉｏｌｏｇ
ｙ，Ｂｏｓｔｏｎ，Ｏｃｔｏｂｅｒ７−１０，１９８６
（１９８７），ｐｐ．６１５−６１８）の論文より知ら
れている。この方法においては、ある画像の画素は変換
後の他の画像内でもその強度及び形状を保持すると仮定
している。変換が実行され、それにより出力画像の全て
の画素に対して入力画像内のそれらの位置が計算され
る。この位置が入力画像内の画素の位置に正確に対応す
る場合、それの入力画像値は出力画像の画素に対する出
力画像値として割り当てられる。しかしながら、変換さ
れた画素の代わりに入力画像内の幾つかの（通常４つ
の）画素が部分的に囲われる場合、出力画像の画素は変
換された画素により部分的に囲われた入力画像の画素に
関連した入力画像値の重み付けされた和に対応する出力
画像値に割り当てられる。この所謂双一次（ｂｉｌｉｎ
ｅａｒ）の補間は歪み補正された出力画像の質を高める
が（出力画像の画素が変換された画素にもっとも近い入
力画像内の画素の入力画像値のみに割り当てられている
場合と比較して）、画質は、特に強い歪みのある場合
は、まだ満足できるものではない。
【０００４】
【発明が解決しようとする課題】本発明の目的は、出力
画像の画質を向上させる上記のような方法を提供するこ
とにある。
【０００５】
【課題を解決するための手段】第一の好ましい解決方法
は、入力画像内の出力画像の画素の角の点の位置が、入
力画像の角の点によって決められた多角形の面領域と同
様に決定され、出力画像の画素は関連した多角形により
囲われた入力画像面領域上の入力画像値の積分値に対応
するそれぞれの出力画像値を割り当てられる。
【０００６】第二の好ましい解決方法は、入力画像の画
素の角の点の位置は、出力画像内及び出力画像内の角の
点によって決められた多角形の面領域で決定され、出力
画像の画素内に位置する角の点の画素の入力画像値が出
力画像の関連する画素の出力画像値を決定するために重
み付けされた和により用いられる。本発明は、知られて
いる方法の基礎となっている仮定、すなわち、画素は変
換後でもその強度及び形状を保持するという仮定は画像
内の歪みが大きくなるにつれてより不正確になる単なる
近似であるという事実の認識に基づく。これは、強い歪
みのある場合は（例えば）ある画像内の正方形の画素は
変換後に他の画像内の菱形又は台形といった矩形にな
り；該矩形の面領域は入力画像内の関連した画素の面領
域から離れるためである。形状及び変換による画素の面
領域の変化に固有の不正確さは、変換された画像の位置
の決定は単一の点（例えば、画素の中心）のみに基づく
のではなく関連する画素の全ての角の点（すなわち正方
形の画素の場合四つの角の点）に基づく本発明により回
避される。
【０００７】これらの考察に基づいた解決方法は、第一
の解決方法が入力画像内のそれらの角の点の位置の出力
画像の画素に基づき、一方第二の解決方法が出力画像内
のそれらの角の点の位置の入力画像の画素に基づいてい
ることで区別される。第一の解決方法によれば、関連す
る出力画像画素に関連する出力画像値は直接得られるの
に対し、第二の解決方法によれば、幾つかの画像値の和
をとらなければならない（重み付けされた方法で）。故
に、第一の解決方法がより簡単で、好ましい。
【０００８】本発明による方法が実施され、それにより
始めに単一の画素のみが考慮され、続いて出力画像値又
は変換によりそれから得られる多角形がこの画素に対し
て得られ、更なる画素の処理が続き等々となる場合、各
画素に対して各々の他の画像内のそれの角の点（通常４
点）の位置が決定されなければならないため、比較的高
い計算負荷が要求される。この負荷は、入力画像の画像
値が出力画像のそれに割り当てられる前に、まず全ての
画素の角の点の位置が決定され、格納される本発明によ
る更なる実施例により減少できる。この解決方法は各角
の点は、画像の端に位置する画素の角の点を除いて幾つ
かの画素に属し；長方形の画素の場合、画像内に位置す
る角の点は４つの画素に属するという仮定に基づく。従
ってこれらの角の点の位置は一度だけ決定される必要が
あり、それにより計算は実質的に加速される。
【０００９】しかしながら、この解決方法によれば、入
力画像及び出力画像だけでなく、各々の他の画像内の一
の画像の角の点の位置もまた格納される必要がある。こ
のような格納の負荷は、画素の第一の行又は列に関係す
る角の点の位置が計算され、格納され、続いてこの行又
は列の画素及び他の画像内の対応する多角形の割り当て
がなされ、続いて第一の行又は列の画素のこれらの画素
に関連した角の点が用いられる一方で第一の行又は列に
隣接する画素の第二の行又は列の角の点の位置が決定さ
れ、その後で画素の第二の行又は列に対して画像値の割
り当てがなされる本発明により減少できる。この解決方
法は画素の二つの行（又は列）の角の点に対してのみ中
間的な格納が要求される。
【００１０】本発明による方法を実施する装置は、 ａ） Ｘ線照射により入力画像を形成するＸ線装置と、 ｂ） 入力画像を入力画像を構成する一連のデジタル画
像値に変換する画像変換器と、 ｃ） 入力画像値及び出力画像値を格納する格納装置
と、 ｄ） 他の画像内の角の点により決められた多角形の面
領域と同様に他の画像内の一の画像の角の点の位置を決
定する手段と、 ｅ） 他の画像内の多角形に関連した画素の画像値を割
り当てる手段とからなることを特徴とする。
【００１１】
【実施例】以下に図を参照しながら本発明を詳細に説明
する。図１に示すＸ線装置は、テーブルトップ２上に位
置する患者３を照射するＸ線源１からなる。かくして生
成されたＸ線レリーフは輝度増強出力画像がテレビジョ
ンカメラ５により収集される画像増強装置４の入口スク
リーンに印加される。画像増強装置４の入口スクリーン
は湾曲しているので、それにより画像増強装置から供給
された画像はクッション状の歪みを示す。更に、より少
ない歪みは地磁気及びビデオカメラ５により引き起こさ
れる。
【００１２】ビデオカメラ５により供給されたアナログ
映像信号はアナログーデジタル変換器６によってメモリ
ー１０内の連続したアドレスに格納されるデジタルデー
タ語の系列に変換される。この過程はビデオ制御器７に
より制御される。斯くして歪んだ入力画像を含むメモリ
ー１０はバスシステム９を介してマイクロコンピュータ
８に接続され、該マイクロコンピュータはこのバスシス
テムを介して更なるメモリー１１及び１２にもアクセス
可能である。メモリー１１はそれから変換により歪みを
除去される出力画像を格納する。デジタルーアナログ変
換器１３を介して図示されない表示ユニット（モニタ
ー）に印加されるために、出力画像値はビデオ制御ユニ
ット７を介してメモリー１１から読み出されうる。
【００１３】Ｘ線照射後には、メモリー１０は歪んだ入
力画像（Ｂi ）を含む。変換によるメモリー１１内の歪
みのない画像（Ｂo ）の形成を以下に図２及び３を参照
して詳細に説明する。図２に入力画像Ｂi 及び出力画像
Ｂo の画像マトリックスの一部を示す。該部分は３ｘ３
の領域を構成する。通常は、画素は正方形の形状を有す
る。しかしながら、原理的には他の形、例えば長方形又
は六角形もまた可能である。入力画像Ｂi の画素は、そ
れらの位置、すなわちそれらの座標ｕ，ｖ（ｕ（ｕ1 ，
ｕ2 ，ｕ3 ．．．）及びｖ（ｖ1 ，ｖ2 ，
ｖ3 ．．．））を特徴とする。同様に、座標ｘ（ｘ1 ，
ｘ2 ，ｘ3 ．．．）及びｙ（ｙ1 ，ｙ2 ，ｙ3 ．．．）
は出力画像内の画素の位置を特徴とする。出力画像内の
画素は線ｉ（ｉ1 ，ｉ2 ，ｉ3 ．．．）により水平方向
に境界をなし、垂直方向には線ｊ（ｊ1 ，ｊ2 ，
ｊ3 ．．．）により境界をなす。このような境界線はま
た入力画像Ｂi にも存在するが、本発明の方法に対して
示される必要はない。二つの線ｉ，ｊの交差の点は各々
の角の点を決める。ＮｘＮ画素からなる画像内には（Ｎ
＋１）ｘ（Ｎ＋１）の角の点がある。実際Ｎの総計は、
例えば５１２に達する。
【００１４】入力画像Ｂi 及び出力画像Ｂo 内の各々の
画素（又はその中心）の位置を特徴とするｕ，ｖ及び
ｘ，ｙに対する離散量は好ましくは整数値である。入力
画像Ｂ i 内の組み合わせｕ，ｖ又は出力画像Ｂo 内の組
み合わせｘ，ｙの各々はメモリー１０及びメモリー１１
内の所定のアドレスを有するメモリー位置に各々割り当
てられる。
【００１５】本発明の第一のステップは、出力画像の画
素に基づいて、入力画像の画素との空間的な関係を決定
することよりなる。この目的のために、知られている方
法は画素につき単一の点の位置、例えば中心を計算し、
該位置及び強度は変化しないままであると仮定するが、
一方本発明による方法では画素につき幾つかの点の位
置、すなわち全ての（４つの）角の点の位置を計算す
る。しかしながら、これらの計算の負荷はわずかに増加
するのみである。なぜならば画像の端の角の点を除いた
画素の各々の角の点はまた３つの更なる画素に属してい
るからである。ＮｘＮの計算の代わりに本発明では（Ｎ
＋１）ｘ（Ｎ＋１）の計算が要求される。
【００１６】図３を参照して説明される方法に対して、
先行するプログラム区画１０１中で入力画像内の線ｊ＝
ｊ1 上の出力画像の全ての角の点の座標ｕ，ｖが計算さ
れ、メモリー１２内の区画Ａに格納されると仮定してい
る。次のステップ１０２中で、座標ｕ，ｖは次の線（ｊ
＝ｊ2 ）上に位置する出力画像内の角の点として計算さ
れる。Ｔｅｈｒａｎｉ等による論文から知られているよ
うに、この計算は式（１）及び（２）に従って多項式に
よりなされる： ｕ ＝ ａ0 ＋ａ1 ｉ＋ａ2 ｊ＋ａ3 ｉｊ＋ａ4 ｉ2 ＋ａ5 ｊ2 （１） ｖ ＝ ｂ0 ＋ｂ1 ｉ＋ｂ2 ｊ＋ｂ3 ｉｊ＋ｂ4 ｉ2 ＋ｂ5 ｊ2 （２） 係数ａ0 ．．．ａ5 及びｂ0 ．．．ｂ5 は、格納され、
全ての後に続く画像の計算に用いるよう、ラスター様の
試験対象により予め決定される。強い歪みの場合、計算
は例えば高次の多項式によるなど、高い精度でなされな
ければならない。かくして計算された座標ｕ，ｖはメモ
リー１２内の区画Ｂに格納される。
【００１７】それで行ｙ＝ｙ1 内の画素の全ての角の点
の位置は格納され、どの出力画像値が出力画像の画素に
割り当てられなければならないかが計算されうる。出力
画像の画像値は以下にＢo （ｘ，ｙ）と称し、入力画像
の画像値はＢi （ｕ，ｖ）と称する。図３に座標ｘ2 ，
ｙ1 を有する出力画像画素の４つの角の点が入力画像内
の４つの点Ｐ1 ，Ｐ2 ，Ｐ3 ，Ｐ4 になることを示す。
４つの角の点を相互接続する出力画像内の４つの直線は
また入力画像内でも直線になり、それにより多角形、す
なわち４つの角の点の場合は矩形が得られるとまた仮定
する。この仮定は完全に正確ではないが（点Ｐ1 ．．．
Ｐ4 間を結ぶ線はまた曲線であり得る）、適切な近似を
なしうる。
【００１８】ブロック１０３の計算は出力画像値は四角
形Ｐ1 ．．．Ｐ4 により囲まれた入力画像画素の入力画
像値の重み付けされた和により形成される出力画像値で
あるという考えに基づく。入力画像値は関連する画像内
の四角形により囲われた面領域がより大きくなるに従っ
てより大きくなる重みと共に出力画像値に寄与する。こ
れらの面領域が図２に従って符号Ｐ1 ．．．Ｐ4 で示さ
れる場合、以下の式が出力画像値Ｂo （ｘ2 ，ｙ1 ）に
対して得られる： Ｂo （ｘ2 ，ｙ1 ） ＝ ｃ（Ｆ1 Ｂi （ｕ1 ，ｖ1 ）＋Ｆ2 Ｂi （ｕ2 ，ｖ1 ）＋Ｆ3 Ｂi （ｕ1 ，ｖ2 ）＋Ｆ4 Ｂi （ｕ2 ，ｖ2 ））／Ｆ （３） ここで、ｃは一画像の全ての画素に対して同一な適切に
選ばれた定数であり、Ｆは入力画像内の四角形
（Ｐ1 ．．．Ｐ4 ）の全面領域である（Ｆ＝Ｆ1 ＋Ｆ2
＋Ｆ3 ＋Ｆ4 ）。一般的に言って、出力画像値は四角形
の面領域にわたる入力画像値の積分値に対応する。
【００１９】この計算はステップ１０３中で行ｙ＝ｙ1
の全ての画素に対してなされ、かくして得られた出力画
像値Ｂo （ｘ，ｙ1 ）は出力画像Ｂo としてメモリー１
１内に格納される。続いて、角の点の最後の行（ｊ＝Ｎ
＋１）が既に処理されたかどうかをチェックする（１０
４）。そうでなければ、行数は１だけ加算され（例えば
現在の例としてｊ＝ｊ3 まで）、同様に画素行の数も
（ｙ＝ｙ2 （ステップ１０５））加算される。
【００２０】出力画像の角の点の新たな行に対して、入
力画像内のそれらの位置ｕ，ｖはステップ１０２中で再
び計算される。行ｊ＝ｊ1 に対して計算された値は斯く
して計算された値によりメモリー１２の区画Ａ内に上書
きされる。このように、メモリー１２の区画Ａ及びＢは
行ｙ＝ｙ2 の画素を定義する行ｊ＝ｊ2 及びｊ＝ｊ3に
対する角の点の座標ｕ，ｖを格納する。関連した出力画
像値はステップ１０３中で計算され、その後チェック１
０４が再び実施される。ループ１０２．．．１０５がも
う一度完了したときに、ステップ１０２中で計算された
値ｕ，ｖは区画Ｂ及びＡに格納されている値を交互に消
去するよう上書きされる。チェックにより最後の行（ｊ
＝Ｎ＋１）が既に処理されたことが判明した場合、出力
画像の全ての画像値が計算され、プログラムは終了する
（１０６）。
【００２１】角の点の位置はまた行順ではなく列順に計
算されうる（ｉ＝ｉ1 ，ｉ2 等々）。この場合、列（ｘ
＝ｘ1 ）に対する出力画像値はステップ１０３で連続的
に計算され、それからステップ１０４中で角の点の最終
列が既に処理されたかどうか（ｉ≧Ｎ？）チェックされ
ねばならない。計算負荷は理論的には一度出力画像の全
ての角の点に対して座標ｕ，ｖを計算することにより、
及びこれらの座標を格納する（アドレスリストの形で）
ことにより減少される。しかしながら、特に複数のＸ線
画像が異なる放射源位置から画像増強装置４に入射する
Ｘ線により形成される断層合成（ｔｏｍｏｓｙｎｔｈｅ
ｓｉｓ）の場合のように異なる歪みを示す入力画像が処
理される場合、これは実質的に要求される格納容量を増
加させる。故に、各照射に対する角の点の位置の更新さ
れた計算が好ましく、何故ならばその場合は唯一の多項
式の係数の組ａ0．．．ａ5 及びｂ0 ．．．ｂ5 が各入
力画像に対して格納される必要があるからである（式
（１）及び（２）参照）。
【００２２】第二の解決方法を図４、５を参照して以下
に説明する。図４に示すように、この解決方法は入力画
像の角の点及びこれらの出力画像の角の点が最初に決定
される位置に基づく。入力画像の４つの角の点は幾つか
の画素（この場合４つ）に囲まれた出力画像内の四角形
（例えばＰ1 ．．．Ｐ4 ）を定義する。これらの画素の
各々に対して出力画像値（の一部）は、関連した入力画
像画素の入力画像値から得られ、該出力画像値（の一
部）は四角形の面領域全体に対する四角形により囲まれ
た画素の面領域の比に対応する。
【００２３】図４と関連した図５を参照して説明される
図４の変換方法に対して、メモリー１２内の区画Ａ内で
出力画像内の関連した座標ｘ，ｙの位置は入力画像の角
の点（例えばｊ＝ｊ2 ）の第一の行に対して計算され格
納される（ブロック２０１）と仮定している。第一のテ
ップ２０２中で、出力画像内の関連した座標ｘ，ｙは入
力画像内の次の行（ｊ＝ｊ3 ）に対して計算される。ｕ
がｘで、ｖがｙで置き換えられる場合、この計算は式
（１）及び（２）に従って多項式によりなされうるが、
逆関数が関係する故に、多項式の係数ａ0 ．．．ａ5 及
びｂ0 ．．．ｂ5は第一の解決方法に従った方法により
得られた多項式の係数と対応しない。斯くして計算され
たｘ，ｙはメモリー１２の区画Ｂに書き込まれる（図
１）。
【００２４】それでメモリー１２は画素の一行を共に定
義する角の点の二行の位置を格納する。続いて、出力画
像内の四角形により囲まれた画素に対して、出力画像値
ｄＢo（の一部）は入力画像内の関連した画素の入力画
像値から得られる（ブロック２０３）。図４に従って、
例えば座標ｘ＝ｘ2 又はｘ＝ｘ3 及びｙ＝ｙ1 又はｙ＝
ｙ2 を有する出力画像画素は入力画像値Ｂi （ｕ2 ，ｖ
2 ）により影響される。例えば、出力画素ｘ2 ，ｙ2 に
対して： ｄＢo （ｘ2 ，ｙ2 ）＝ｃ．Ｆ4 ／Ｆ0 ．Ｂi （ｕ2 ，ｖ2 ） （４） が得られる。ここで、ｃは定数、Ｆ0 は出力画像内の点
Ｐ1 ．．．Ｐ4 により定義された四角形の全面領域、Ｆ
4 は四角形に囲まれた画素ｘ2 ，ｙ2 の部分面領域であ
る。
【００２５】しかしながら、図４に示したように、出力
画像値Ｂo （ｘ2 ，ｙ2 ）は入力画像値Ｂi （ｕ2 ，ｖ
2 ）のみならず、入力画像内の画素ｕ2 ，ｖ2 の下に
及び／又は隣接して位置する画素に関連した入力画像値
Ｂi （ｕ1 ，ｖ2 ）、Ｂi （ｕ1 ，ｖ3 ）、Ｂ
i （ｕ2 ，ｖ3 ）にも依存する。従ってステップ２０４
中で値ｄＢo はメモリー１１内の関連した出力画像画素
に対してメモリー１１内に格納された出力画像値に加え
られる。
【００２６】角の点の最後の行がまだ処理されない場合
（チェック２０５）、値ｊ及びｖは１増加し（ｊ＝
ｊ4 、ｖ＝ｖ3 （ステップ３０６））、ステップ１０２
中の出力画像内のこの行の角の点の座標ｘ，ｙが計算さ
れる。計算された値はメモリー１２内の区画Ａ内に最後
から２番目（ｊ＝ｊ2 ）として格納された座標に上書き
される。続いてステップ２０３の計算処理は次の行（ｖ
＝ｖ2 ）に対して繰り返され、斯くして計算された値は
メモリー１１に既に格納された値に加えられる。プログ
ラムループ２０４から２０６が再び完了した場合、新た
に計算された座標値はメモリー１２内の区画Ｂ内に最後
から２番目（ｊ＝ｊ3 ）の座標に上書きされる。ステッ
プ２０５中で入力画像の画素の全ての行がこのような方
法で扱われることが確立した場合、この方法は終了す
る。
【００２７】出力画像値は第二の解決方法によれば加算
又は累積によりステップ順にのみ形成され、一方完全な
画像値は図２及び３による方法による各ステップ（１０
３）中に形成されるので、第一の方法が好ましい。明ら
かに、第一の方法と類似で、出力画像内の角の点の座標
は行順に決定されうる代わりに列順に決定されうる。更
に、対応したより大きな格納容量を用いることにより、
後に続く歪み除去操作に用いるために、計算された座標
は一度計算され、アドレスリストの形で格納されうる。
【図面の簡単な説明】
【図１】本発明を実施する装置のブロック系統図であ
る。
【図２】第一の解決方法に対する入力画像及び出力画像
の詳細を示す図である。
【図３】対応するフローチャートを示す図である。
【図４】第二の解決方法に対する入力画像及び出力画像
の詳細を示す図である。
【図５】対応するフローチャートを示す図である。
【符号の説明】
１ Ｘ線源 ２ テーブルトップ ３ 患者 ４ 画像増強装置 ５ テレビジョンカメラ ６ アナログデジタル変換器 ７ ビデオ制御器 ８ マイクロコンピュータ ９ バスシステム １０、１１、１２ メモリー １３ デジタルアナログ変換器 Ｂi 入力画像 Ｂo 出力画像
───────────────────────────────────────────────────── フロントページの続き (51)Int.Cl.6 識別記号 庁内整理番号 ＦＩ 技術表示箇所 9287−5Ｌ Ｇ０６Ｆ 15/62 ３９０ Ａ

Claims (5)

    【特許請求の範囲】
  1. 【請求項１】 Ｘ線照射により生成され位置（ｕ，ｖ）
    及び各々の入力画像値（Ｂi （ｕ，ｖ））により定義さ
    れた画素により構成された歪んだ入力画像（Ｂi ）を、
    位置（ｘ，ｙ）及び各々の出力画像値（Ｂo （ｘ，
    ｙ））により定義された画素により構成された実質的に
    歪みのない出力画像（Ｂo ）に変換する方法であって、
    入力画像（Ｂi ）内の出力画像（Ｂo ）の画素の角の点
    （ｉ，ｊ）の位置（ｕ，ｖ）が、入力画像の角の点によ
    って決められた多角形の面領域（Ｆ 1 ．．．Ｆ4 ）と同
    様に決定され、出力画像の画素は関連した多角形により
    囲われた入力画像面領域上の入力画像値の積分値に対応
    するそれぞれの出力画像値（Ｂo （ｘ，ｙ））を割り当
    てられることを特徴とする方法。
  2. 【請求項２】 Ｘ線照射により生成され位置（ｕ，ｖ）
    及び各々の入力画像値（Ｂi （ｕ，ｖ））により定義さ
    れた画素により構成された歪んだ入力画像（Ｂi ）を、
    位置及び各々の出力画像値により定義された画素により
    構成された実質的に歪みのない出力画像（Ｂo ）に変換
    する方法であって、入力画像の画素の角の点（ｉ，ｊ）
    の位置（ｘ，ｙ）が、出力画像内及び出力画像内の角の
    点によって決められた多角形の面領域（Ｆ1 ．．．
    Ｆ4 ）で決定され、出力画像の画素内に位置する角の点
    の画素の入力画像値（例えば、Ｂi （ｕ2 ，ｖ2 ））が
    出力画像の関連する画素の出力画像値Ｂo （ｘ2 ，
    ｙ2 ）を決定するために重み付けされた和により用いら
    れることを特徴とする方法。
  3. 【請求項３】 入力画像及び出力画像の画像値が相互に
    割り当てられる前に、まず全ての画素の角の点の位置が
    決定され、格納されることを特徴とする請求項１記載の
    方法。
  4. 【請求項４】 画素の第一の行又は列に属する角の点の
    位置が計算され、格納され、続いてこの行又は列の画素
    及び他の画像内の対応する多角形の割り当てがなされ、
    続いて第一の行又は列の画素のこれらの画素に関連した
    角の点が用いられる一方で第一の行又は列に隣接する画
    素の第二の行又は列の角の点の位置が決定され、その後
    で画素の第二の行又は列に対して画像値の割り当てがな
    されることを特徴とする請求項１記載の方法。
  5. 【請求項５】ａ） Ｘ線照射により入力画像を形成する
    Ｘ線装置（１、４）と、 ｂ） 入力画像を入力画像（Ｂi ）を構成する一連のデ
    ジタル画像値に変換する画像変換器と、 ｃ） 入力画像値及び出力画像値を格納する格納装置
    （１０、１１）と、 ｄ） 他の画像内の角の点により決められた多角形の面
    領域（Ｆ1 ．．．Ｆ4 ）と同様に他の画像内の一の画像
    （Ｂo ）の角の点の位置（例えばｕ，ｖ）を決定する手
    段（７；１０２；２０２）と、 ｅ） 他の画像内の多角形に関連した画素の画像値を割
    り当てる手段（７；１０３；２０３）とからなることを
    特徴とする請求項１記載の方法を実施する装置。
JP10406094A 1993-05-19 1994-05-18 Ｘ線画像の歪み補正方法及び装置 Expired - Fee Related JP3672588B2 (ja)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
DE4316847A DE4316847A1 (de) 1993-05-19 1993-05-19 Verfahren zum Entzerren von Röntgenaufnahmen und Anordnung zur Durchführung des Verfahrens
DE4316847:7 1993-05-19

Publications (2)

Publication Number Publication Date
JPH0737080A true JPH0737080A (ja) 1995-02-07
JP3672588B2 JP3672588B2 (ja) 2005-07-20

Family

ID=6488529

Family Applications (1)

Application Number Title Priority Date Filing Date
JP10406094A Expired - Fee Related JP3672588B2 (ja) 1993-05-19 1994-05-18 Ｘ線画像の歪み補正方法及び装置

Country Status (4)

Country Link
US (1) US5671297A (ja)
EP (1) EP0625762B1 (ja)
JP (1) JP3672588B2 (ja)
DE (2) DE4316847A1 (ja)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6146579A (en) * 1997-05-27 2000-11-14 Sumitomo Chemical Co., Ltd. Process for producing thermoplastic resin hollow molded article
JP2001078074A (ja) * 1999-09-01 2001-03-23 Hitachi Ltd デジタルカメラ

Families Citing this family (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE19542355A1 (de) * 1995-11-14 1997-05-15 Hans Thilo Richter Verfahren und Vorrichtung zur Erzeugung eines Weitwinkelbildes
US6038548A (en) * 1997-11-26 2000-03-14 International Business Machines Corporation System and method for conducting electronic commerce in a computer network using a cashier desk payment framework
US7067414B1 (en) 1999-09-01 2006-06-27 Micron Technology, Inc. Low k interlevel dielectric layer fabrication methods
US6975755B1 (en) 1999-11-25 2005-12-13 Canon Kabushiki Kaisha Image processing method and apparatus
GB2359884B (en) * 1999-11-25 2004-06-30 Canon Kk Image processing method and apparatus
US6671349B1 (en) 2000-11-13 2003-12-30 Olganix Corporation Tomosynthesis system and registration method
WO2002089059A2 (de) * 2001-04-25 2002-11-07 Siemens Aktiengesellschaft Bildverarbeitungsverfahren
DE10141085A1 (de) * 2001-08-22 2003-05-15 Siemens Ag Röntgendiagnostikeinrichtung mit einer Bildvorverarbeitung und Verfahren zur Entzerrung von Bildern
US7302111B2 (en) * 2001-09-12 2007-11-27 Micronic Laser Systems A.B. Graphics engine for high precision lithography
US6707998B2 (en) 2002-05-14 2004-03-16 Eastman Kodak Company Method and system for correcting non-symmetric distortion in an image
US7327390B2 (en) * 2003-02-04 2008-02-05 Eastman Kodak Company Method for determining image correction parameters
DE10352411B4 (de) 2003-11-10 2007-02-22 Yxlon International Security Gmbh Verfahren zur Entzerrung eines Röntgenbildes eines Gepäckstücks
US20050226535A1 (en) * 2004-04-13 2005-10-13 Srinka Ghosh Method and system for rectilinearizing an image of a microarray having a non-rectilinear feature arrangement
US20050275733A1 (en) * 2004-06-10 2005-12-15 Philip Chao Method and apparatus of rendering a video image by polynomial evaluation
US9229970B2 (en) 2009-12-07 2016-01-05 International Business Machines Corporation Methods to minimize communication in a cluster database system

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE3843232A1 (de) * 1988-12-22 1990-06-28 Philips Patentverwaltung Schaltungsanordnung zur geometrischen bildtransformation
NL9000766A (nl) * 1990-04-02 1991-11-01 Koninkl Philips Electronics Nv Inrichting voor geometrische correctie van een vertekend beeld.
US5173948A (en) * 1991-03-29 1992-12-22 The Grass Valley Group, Inc. Video image mapping system
JP2791231B2 (ja) * 1991-04-17 1998-08-27 三菱電機株式会社 Ｘ線画像診断装置の糸巻歪補正方法
GB2256989B (en) * 1991-06-21 1995-02-08 Sony Broadcast & Communication Video image capture apparatus

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6146579A (en) * 1997-05-27 2000-11-14 Sumitomo Chemical Co., Ltd. Process for producing thermoplastic resin hollow molded article
JP2001078074A (ja) * 1999-09-01 2001-03-23 Hitachi Ltd デジタルカメラ

Also Published As

Publication number Publication date
JP3672588B2 (ja) 2005-07-20
DE4316847A1 (de) 1994-11-24
US5671297A (en) 1997-09-23
DE59409916D1 (de) 2001-11-29
EP0625762B1 (de) 2001-10-24
EP0625762A1 (de) 1994-11-23

Similar Documents

Publication Publication Date Title
JP3672588B2 (ja) Ｘ線画像の歪み補正方法及び装置
EP0450718B1 (en) Apparatus for geometrical correction of a distorted image
US4870692A (en) Device for and method of subtraction imaging with reduction of motion artefacts
JPH01126779A (ja) 補間方法
JP2905059B2 (ja) 色値処理方法および処理装置
US4682301A (en) Digital filter for processing two-dimensional digital image
EP0743617B1 (en) Image processing method and apparatus thereof
JPS6160066A (ja) デイジタル・イメージ記録投影処理方法
CN101027679A (zh) 表达通用的二维空间变换的系统和方法
GB2220540A (en) Device for the digital processing of images to obtain special geometrical effects
JP3240235B2 (ja) 放射線画像の表示方法
US5881163A (en) X-ray examination apparatus with an imaging arrangement having a plurality of image sensors
JP3298732B2 (ja) 放射線画像の部分の表示方法
JPS5831594B2 (ja) 画像表示装置
US5144687A (en) Image processing apparatus including spatial shift variant filter
EP0618719B1 (en) X-ray examination apparatus with an imaging arrangement having a plurality of image sensors
JPH09307789A (ja) 画像処理装置
US5473653A (en) Method of generating slice images, and device for carrying out the method
US5220428A (en) Digital video effects apparatus for image transposition
JP2905060B2 (ja) 色値処理方法および装置
JPH06348835A (ja) 歪んだ入力画像を変換する方法及び装置
US4939758A (en) Apparatus for the generation of tomographic X-ray images
JP2003037769A (ja) 信号処理装置、撮像装置、及び信号処理方法
JPH1153513A (ja) 画像内の領域設定方法及び装置
JP2001134780A (ja) ３次元対象を２次元画像によってイメージングする方法

Legal Events

Date Code Title Description
A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20041022

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20041026

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20050125

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20050420

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

LAPS Cancellation because of no payment of annual fees