JP3599435B2 - 画像データ補間演算方法および装置 - Google Patents
画像データ補間演算方法および装置 Download PDFInfo
- Publication number
- JP3599435B2 JP3599435B2 JP20320795A JP20320795A JP3599435B2 JP 3599435 B2 JP3599435 B2 JP 3599435B2 JP 20320795 A JP20320795 A JP 20320795A JP 20320795 A JP20320795 A JP 20320795A JP 3599435 B2 JP3599435 B2 JP 3599435B2
- Authority
- JP
- Japan
- Prior art keywords
- interpolation
- image data
- original image
- original
- neighboring
- 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 - Fee Related
Links
- 238000004364 calculation method Methods 0.000 title claims description 58
- 238000000034 method Methods 0.000 claims description 18
- 230000008569 process Effects 0.000 claims description 8
- 230000005855 radiation Effects 0.000 description 12
- 238000010586 diagram Methods 0.000 description 10
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 description 9
- 229910052737 gold Inorganic materials 0.000 description 6
- 230000005284 excitation Effects 0.000 description 5
- 229910052774 Proactinium Inorganic materials 0.000 description 4
- 230000008859 change Effects 0.000 description 4
- 229910052745 lead Inorganic materials 0.000 description 4
- 239000000463 material Substances 0.000 description 3
- 230000007423 decrease Effects 0.000 description 2
- NIXOWILDQLNWCW-UHFFFAOYSA-N acrylic acid group Chemical group C(C=C)(=O)O NIXOWILDQLNWCW-UHFFFAOYSA-N 0.000 description 1
- 230000005260 alpha ray Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000005250 beta ray Effects 0.000 description 1
- 210000000988 bone and bone Anatomy 0.000 description 1
- 238000012888 cubic function Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000010894 electron beam technology Methods 0.000 description 1
- 230000005251 gamma ray Effects 0.000 description 1
- 238000000465 moulding Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 229910052709 silver Inorganic materials 0.000 description 1
- 239000004332 silver Substances 0.000 description 1
- -1 silver halide Chemical class 0.000 description 1
- 230000004936 stimulating effect Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T3/00—Geometric image transformations in the plane of the image
- G06T3/40—Scaling of whole images or parts thereof, e.g. expanding or contracting
- G06T3/403—Edge-driven scaling; Edge-based scaling
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
- H04N1/00—Scanning, transmission or reproduction of documents or the like, e.g. facsimile transmission; Details thereof
- H04N1/40—Picture signal circuits
- H04N1/40068—Modification of image resolution, i.e. determining the values of picture elements at new relative positions
Landscapes
- Engineering & Computer Science (AREA)
- Multimedia (AREA)
- Signal Processing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
- Editing Of Facsimile Originals (AREA)
- Radiography Using Non-Light Waves (AREA)
Description
【発明の属する技術分野】
本発明は画像データの補間演算方法および装置に関し、詳細には画像のうち、少なくとも濃度変化が急峻なエッジ部が斜め方向に延びる画像部分について適用する画像データの補間演算方法および装置に関するものである。
【0002】
【従来の技術】
従来より、放射線写真フイルムに記録された放射線画像を光電的に読み取って画像信号を得、この画像信号に適切な画像処理を施した後、画像を再生記録することが種々の分野で行われている。たとえば、後の画像処理に適合するように設計されたガンマ値の低いX線フィルムを用いてX線画像を記録し、このX線画像が記録されたフィルムからX線画像を読み取って電気信号に変換し、この電気信号(画像信号)に画像処理を施した後、コピー写真等に可視像として再生することにより、コントラスト、シャープネス、粒状性等の画質性能の良好な再生画像を得ることが行われている(特公昭61−5193 号参照)。
【0003】
また本願出願人により、放射線(X線,α線,β線,γ線,電子線,紫外線等)を照射すると、この放射線エネルギーの一部が蓄積され、その後可視光等の励起光を照射すると蓄積されたエネルギーに応じて輝尽発光を示す蓄積性蛍光体(輝尽性蛍光体)を利用して、人体等の被写体の放射線画像情報を一旦シート状の蓄積性蛍光体に記録し、この蓄積性蛍光体シートをレーザー光等の励起光で走査して輝尽発光光を生ぜしめ、得られた輝尽発光光を光電的に読み取って画像信号を得、この画像データに基づき被写体の放射線画像を写真感光材料等の記録材料、CRT等に可視像として出力させる放射線画像記録再生システムがすでに提案されている(特開昭55−12429号,同56−11395号,同55−163472 号,同56−104645 号,同55−116340 号等)。このシステムは、従来の銀塩写真を用いる放射線写真システムと比較して極めて広い放射線露出域にわたって画像を記録し得るという実用的な利点を有している。
【0004】
上記のように画像信号を得てこの画像信号に基づいて可視画像を再生するシステムにおいて、その可視画像のうち観察対象となる関心領域をより詳細に観察したいとき、その領域を拡大して再生することがある。この拡大画像は、原画像を読み取って得られた原画像データに対して所定の補間演算を施して原画像データ数とは異なるデータ数の2次的な画像データである補間画像データを求め、この補間画像データに基づいた可視画像の再生を行うことによって得ることができる。
【0005】
ところで、画像入出力装置の構成のし易さの観点から一般的に用いられている、各画像データを担持する画素が所定の間隔で縦横方向に正方形格子状に配列されて画像を形成するようにしたものにおいては、上記画像の拡大処理における補間演算は、その補間画像データを、補間して新たに設定しようとする画素(補間画素)の近傍4点の原画素における原画像データを線形補間することによって行っている。
【0006】
例えば、図8(A)に示すように、正方形格子状に配列された原画像の原画素P(○記号が記された点)について、その原画素Pが配列された間隔とは異なる間隔で配列される補間画素P′(×記号が記された点)の補間画像データを求めようとするときは、例えば補間画素P′0 の補間画像データについては以下の手順によって求める。
【0007】
補間画素P′0 を囲む近傍4点の原画像の原画素PA 、PB 、PC 、PD の画像データSA 、SB 、SC 、SD を用いる。
【0008】
ここで原画素PA 〜PB 間、PC 〜PD 間、PA 〜PC 間、PB 〜PD 間のピッチをそれぞれ1とし、補間画素P′0 の、原画素PA (PC )からのx軸方向(横方向)の距離がTx(図8(B)参照)、原画素PA (PB )からのy軸方向(縦方向)の距離がTyである場合、まず補間画素P′0 のx軸方向の位置に対応する補間画素P′m、P′nの補間画像データS′m、S′nを、下記の式(1)、(2)の線形補間の演算により求める。
【0009】
S′m=(1−Tx)SA +TxSB (17)
S′n=(1−Tx)SC +TxSD (18)
次いで、補間画素P′0 のy軸方向について補間画像データS′m、S′nを用いた下記の式(18)の線形補間の演算を行って、補間画像データS′0 を求める。
【0010】
S′0 =(1−Ty)S′m+TyS′n (19)
以上の演算を他の補間点P′についても同様に適用して、各補間画像データS′を求めることができる。
【0011】
また、画像データの補間法としては、上述した線形補間による方法の他、2次あるいは3次のスプライン補間関数を用いる方法など種々の方法が提案されている。例えば、3次のスプライン補間関数を用いるCubic スプライン補間演算は、元のサンプル点(画素)を通ることと、その第1階微分係数が各区間間で連続することが必要とされており、この条件に基づいて補間点近傍4点の原画像の画素に乗じる補間係数を算出し、この補間係数を補間点近傍4点の原画像の画素に乗じて補間画像データを得る方法である。このCubic スプライン補間演算は、比較的鮮鋭度の高いシャープな2次画像(補間により得られる画像)を再生するためのものである。また、Cubic スプライン補間演算に対して比較的鮮鋭度の低い滑らかな2次画像を再生するための補間画像データを得るBスプライン補間演算も知られている。このように2次画像を高い鮮鋭度でシャープに再生したい場合は、Cubic スプライン補間演算を用い、低い鮮鋭度で滑らかに再生したい場合はBスプライン補間演算を用いればよい。
【0012】
【発明が解決しようとする課題】
ところで、再生された可視画像においては、例えば放射線画像における骨部のような濃度(輝度)の変化が急峻なエッジ部もあり、このようなエッジ部を拡大することもある。
【0013】
しかし、このようなエッジ部が、上記原画像の画素が配列された正方形格子に対して斜め方向に延びる場合に、上述した式(16)〜(18)による線形補間、あるいはCubic スプライン補間、Bスプライン補間を行ったのでは、その斜め方向に延びるエッジ部の拡大画像は階段状の段付きが目立つものとなる。
【0014】
例えば、図9(A)に示すような斜め方向に延びるエッジ部を有する画像においては、微視的には図9(B)に示す高濃度の点(黒丸で示す)の領域と低濃度の点(白丸で示す)の領域との境界線(エッジ)が斜め方向に延びている部分について、前述した補間演算を適用して補間画像データを求めた場合、その補間画像データS′0 は、低濃度の原画像データSD にも依存するため、高濃度の画素PA 、PB 、PC よりも濃度が少し低下した中間濃度を示す画像データとなる(図10(B)参照)。このため、得られた補間画像データS′0 に基づいた拡大画像を再生した場合、エッジ部が図10(B)の破線で示すように階段状の段付きが拡大された画像となる。すなわち画像全体としては、図9(A)に示すように斜め方向に延びるエッジ部は拡大処理によって、図10(A)に示すように、そのエッジ部の階段状の段付きがそのまま拡大されることになる。
【0015】
このようなエッジ部の段付きは、そのエッジ部近傍を観察するうえで画像読影の障害となり、画像の診断性能を低下させる虞がある。
【0016】
本発明は上記事情に鑑み、斜め方向に延びるエッジ部について滑らかに、かつシャープにした拡大画像を得ることができる画像データの補間演算方法および装置を提供することを目的とするものである。
【0017】
【課題を解決するための手段】
本発明による画像データ補間演算方法および装置は、所定の間隔で縦横方向に格子状に配された原画素のそれぞれについての画素値を表す多数の原画像データについて、該原画素とは間隔が異なる補間画素における補間画像データを求める画像データの補間演算方法であって、該補間画素近傍の複数点の近傍原画素における近傍原画像データに基づいて該各近傍原画像データに乗じる補間係数をそれぞれ算出するとともに、該各補間係数を該近傍原画像データにそれぞれ乗じる補間演算処理を行うことにより、前記補間画像データを算出する画像データ補間演算方法において、前記原画像データにより表される原画像上の前記補間画素における濃度ベクトルを算出し、該濃度ベクトルに直交する線分と、前記各近傍原画素との距離を算出し、該濃度ベクトルおよび/または該距離が比較的大きいものとなる前記近傍原画素ほど、前記近傍原画像データに乗じる補間係数が小さくなるように該補間係数を補正し、該補正された補間係数に基づいて前記補間演算処理を行うことを特徴とするものである。
【0018】
ここで、補間画素近傍の複数点の近傍原画素とは、具体的には、原画素のうち補間画素に最も近い複数点(例えば4点)の画素のことをいう。
【0019】
また、補間係数は、線形のものの他、3次以上の高次のものや、各種の関数形式で表されたものを算出するようにしてもよい。
【0020】
【発明の効果】
原画像における濃度ベクトルは、濃度変化が大きい部分ほど大きくなるものであり、濃度ベクトルの大きさが比較的大きい場合、この濃度ベクトルに直交する線分は、画像中において濃度変化が急峻なエッジ部を表すものとなる。ここで、補間画素の補間画像データを求める際の、補間画素における濃度ベクトルの大きさと、この濃度ベクトルに直交する線分と補間画素近傍の原画素との距離について考えると、濃度ベクトルの大きさが大きい場合は補間画素はエッジ部にあり、濃度ベクトルに直交する線分はエッジを表すものと考えることができる。したがって、この線分からの距離が離れている原画素の原画像データの補間係数をそのまま用いると、エッジ部が画素の配列に対して斜め方向に延びている場合、エッジから離れた位置にあるデータを補間演算に用いることとなるため、補間画像データがエッジとは異なる濃度を有する画素に影響されてしまい、補間画像データはエッジ部の濃度とは異なるものとなる。このため得られた補間画像データに基づいて画像を再生した場合、斜めのエッジ部が原画像と同様の階段状の段付きが拡大された画像となってしまう。
【0021】
本発明はこの点に鑑みてなされたものであり、算出された補間係数をそのまま補間画素近傍の原画素における原画像データに乗じるのではなく、補間画素における濃度ベクトルの大きさおよび/または濃度ベクトルに直交する線分からの距離の大きさに応じて補正して原画像データに乗じるようにしたものである。すなわち、濃度ベクトルの大きさが大きいほど、および/または濃度ベクトルに直交する線分からの距離が大きいほど補間係数が小さくなるように補正するようにしたものである。したがって、エッジ部が斜めに延びている場合に、エッジ部にある補間画像データはエッジから離れた位置にある原画像データに依存する割合が低くなり、エッジ近傍すなわちエッジに沿った方向にある原画像データに依存する割合が高くなる。これにより、エッジ部にある補間画素はエッジ部と略同様の濃度に補間されることとなるため、原画像の階段状の段付きが拡大されることが無くなり、斜め方向に延びるエッジを階段状の段付きが生じることなく拡大することができる。
【0022】
なお、補間画素がエッジ部でない場合は、濃度ベクトルの大きさは比較的小さくなるため、原画像データが濃度ベクトルに直交する直線から比較的離れた位置にあっても、補間係数はそれほど大きく補正されない。したがって、補間画像データは従来の補間係数とほとんど変わらない補間係数により求められることとなる。
【0023】
また、予め画像におけるエッジ部とエッジ部以外とを判別している場合は、濃度ベクトルの大きさは考慮せず、濃度ベクトルと直交する線分からの距離のみに基づいて補正を行うようにしてもよい。
【0024】
【発明の実施の形態】
以下図面を参照して本発明の実施の形態について説明する。
【0025】
図1は本発明による画像データ補間演算装置を表す図である。図1に示すように本発明による画像データ補間演算装置1は、後述する原画像データを入力するための画像入力手段2から入力された原画像データに基づいて補間画素における補間係数を算出する補間係数算出手段3と、補間画素における原画像の濃度ベクトルを算出するとともに、この濃度ベクトルと直交する線分と補間に用いる原画素との距離を算出する手段4と、手段4において算出された濃度ベクトルの大きさおよび/または濃度ベクトルに直交する線分と補間に用いる原画素との距離の大きさに基づいて補間係数算出手段3において算出された補間係数を補正するための補正項を算出する補正項算出手段5と、補正項算出手段5において算出された補正項に基づいて補間係数を補正する補正手段6と、補正手段6において補正された補間係数により補間画像データを算出する補間手段7と、補間手段7において算出された補間画像データを可視像として再生するための画像出力手段8とからなるものである。
【0026】
画像入力手段2は、例えば図2に示すような放射線画像読取装置からなるものである。図2に示すように、放射線画像読取装置は、例えばX線等の放射線が人体等の被写体を介して照射されることによりこの被写体の透過放射線画像情報を蓄積記録した蓄積性蛍光体シート10は、エンドレスベルト等のシート搬送手段11により、副走査のために矢印Y方向に搬送される。半導体レーザ等の励起光源12から射出された励起光(読取光)としてのレーザビーム13は、高速回転する回転多面鏡14によって反射偏向され、通常f・θレンズからなる走査レンズ18によって集束され、ミラー19で反射して蓄積性蛍光体シート10上を上記副走査方向Yと略直角な矢印X方向に主走査する。
【0027】
こうしてレーザビーム13が照射されたシート10の箇所からは、蓄積記録されている放射線画像情報に応じた光量の輝尽発光光15が発散され、この輝尽発光光15は集光体16によって集光され、光検出器としてのフォトマルチプライヤー(光電子増倍管)17によって光電的に検出される。
【0028】
上記集光体16はアクリル板等の導光性材料を成形して作られたものであり、直線状をなす入射端面16aが蓄積性蛍光体シート10上のビーム走査線に沿って延びるように配され、円環状に形成された出射端面16bに上記フォトマルチプライヤー17の受光面が結合されている。上記入射端面16aから集光体16内に入射した輝尽発光光15は、該集光体16の内部を全反射を繰り返して進み、出射端面16bから出射してフォトマルチプライヤー17に受光され、前記放射線画像情報を担持する輝尽発光光15の光量がこのフォトマルチプライヤー17によって検出される。
【0029】
フォトマルチプライヤー17のアナログ出力信号(画像信号)Sは対数増幅器20によって増幅され、A/D変換器21において所定の収録スケールファクターでデジタル化される。こうして得られた、2次元画像を担持するデジタルの原画像データは、画像データ補間演算装置1に入力される。
【0030】
画像データ補間演算装置1においては、まず補間係数算出手段3において補間係数が算出される。以下、この補間係数の算出について説明する。なお、本実施の形態においては、3次のCubic スプライン補間演算により補間係数を算出するものとする。なお、本実施の形態において得られる原画像データは正方形格子状に配列されているが、ここでは説明のため1次元状に並べられた原画像データを例に挙げて説明する。
【0031】
原画像からデジタル的に読み取って得られた、連続する画素Xk−2 ,Xk−1 ,Xk ,Xk+1 ,Xk+2 ,…の原画像データを図3に示すようにそれぞれYk−2 ,Yk−1 ,Yk ,Yk+1 ,Yk+2 ,…とする。ここで、3次のスプライン補間関数は、各区間Xk−2 〜Xk−1 ,Xk−1 〜Xk ,Xk 〜Xk+1 ,Xk+1 〜Xk+2 ごとにそれぞれ設定され、各区間に対応するスプライン補間関数をfk−2 ,fk−1 ,fk ,fk+1 ,fk+2 とする。この補間関数はいずれも各区間の位置を変数とする3次関数である。
【0032】
ここでまず、補間しようとする点(以下、補間点という)Xp が区間Xk 〜Xk+1 の範囲にある場合について説明する。なお、区間Xk 〜Xk+1 に対応するスプライン補間関数fk は下記の式(1)で表される。
【0033】
fk (x)=Ak x3 +Bk x2 +Ck x+Dk (1)
Cubic スプライン補間演算においては、スプライン補間関数fk は元のサンプル点(画素)を通ることと、その第1階微分係数が各区間間で連続することが必要とされ、これらの条件から下記の式(2)〜(5)を満たす必要がある。
【0034】
fk (Xk )=Yk (2)
fk (Xk+1 )=Yk+1 (3)
fk ′(Xk )=fk−1 ′(Xk ) (4)
fk ′(Xk+1 )=fk+1 ′(Xk+1 ) (5)
なお、fk ′は関数fk の第1階微分(3Ak x2 +2Bk x+Ck )を表すものである。
【0035】
またCubic スプライン補間演算においては、画素Xk における第1階微分係数が、その画素Xk の前後の画素であるXk−1 とXk+1 とについて、これらの原画像データYk−1 、Yk+1 の勾配(Yk+1 −Yk−1 )/(Xk+1 −Xk−1 )に一致することが条件であるから、下記の式(6)を満たす必要がある。
【0036】
fk ′(Xk )=(Yk+1 −Yk−1 )/(Xk+1 −Xk−1 ) (6)
同様に、画素Xk+1 における第1階微分係数が、その画素Xk+1 の前後の画素であるXk とXk+2 とについて、これらの原画像データYk 、Yk+2 の勾配(Yk+2 −Yk )/(Xk+2 −Xk )に一致することが条件であるから、下記の式(7)を満たす必要がある。
【0037】
fk ′(Xk+1 )=(Yk+2 −Yk )/(Xk+2 −Xk ) (7)
ここで、各区間Xk−2 〜Xk−1 ,Xk−1 〜Xk ,Xk 〜Xk+1 ,Xk+1 〜Xk+2 の間隔(格子間隔という)を1とし、画素Xk からの画素Xk+1 方向への補間点Xp の位置をt(0≦t≦1)とすれば、式(1)〜(4)および(6)、(7)より、
fk (0)=Dk =Yk
fk (1)=Ak +Bk +Ck +Dk =Yk+1
fk ′(0)=Ck =(Yk+1 −Yk−1 )/2
fk ′(1)=3Ak +2Bk +Ck =(Yk+2 −Yk )/2
したがって、
Ak =(Yk+2 −3Yk+1 +3Yk −Yk−1 )/2
Bk =(−Yk+2 +4Yk+1 −5Yk +2Yk−1 )/2
Ck =(Yk+1 −Yk−1 )/2
Dk =Yk
なお、スプライン補間関数fk (x)は上述の通り、X=tなる変数変換をしているため、
fk (x)=fk (t)
となる。よって、補間点Xp における補間画像データYp は、
Yp =fk (t)=Ak t3 +Bk t2 +Ck t+Dk (8)
で表すことができる。ここで上記各係数Ak ,Bk ,Ck ,Dk を式(8)に代入すると、
Yp ={(Yk+2 −3Yk+1 +3Yk −Yk−1 )/2}t3
+{(−Yk+2 +4Yk+1 −5Yk +2Yk−1 )/2}t2
+{(Yk+1 −Yk−1 )/2}t+Yk
となり、これを原画像データYk−1 ,Yk ,Yk+1 ,Yk+2 について整理すると、下記の式(9)で表すことができる。
【0038】
Yp ={(−t3 +2t2 −t)/2}Yk−1
+{(3t3 −5t2 +2)/2}Yk
+{(−3t3 +4t2 +t)/2}Yk+1
+{(t3 −t2 )/2}Yk+2 (9)
と表すことができる。
【0039】
ここで、原画像データYk−1 ,Yk ,Yk+1 ,Yk+2 の各係数が補間係数ak−1 ,ak ,ak+1 ,ak+2 となる。すなわち、式(9)における原画像データYk−1 ,Yk ,Yk+1 ,Yk+2 にそれぞれ対応する補間係数ak−1 ,ak ,ak+1,ak+2 は、
ak−1 =(−t3 +2t2 −t)/2
ak =(3t3 −5t2 +2)/2
ak+1 =(−3t3 +4t2 +t)/2
ak+2 =(t3 −t2 )/2
となる。
【0040】
以上の演算を各区間Xk−2 〜Xk−1 ,Xk−1 〜Xk ,Xk 〜Xk+1 ,Xk+1 〜Xk+2 について繰り返すことにより、原画像データの全体についての補間係数を求めることができる。
【0041】
一方、濃度ベクトル算出および濃度ベクトルに直交する線分との距離を算出する手段4(以下手段4とする)においては、補間画素における濃度ベクトルと、この濃度ベクトルに直交する線分と補間演算に用いる原画像データを担持する原画素との距離が算出される。すなわち、図4に示すように、補間画像データを算出する補間画素P0 ′の近傍16個の画素を4つの領域a,b,c,dに分割し、各領域における原画像データの総和Wa,Wb,Wc,Wdをそれぞれ求め、この総和Wa,Wb,Wc,Wdについて、
濃度ベクトルPV =(Wb−Wa,Wd−Wc)
として、濃度ベクトルPV を算出する。
【0042】
あるいは図5に示すように、補間画像データを算出する補間画素P0 ′の近傍4画素Pa,Pb,Pc,Pdの原画像データSa,Sb,Sc,Sdについて、
濃度ベクトルPV =(Sb−Sa,Sd−Sc)
として、濃度ベクトルPV を算出するようにしてもよい。
【0043】
このようにして濃度ベクトルPV が算出されると、次いでこの濃度ベクトルPV に直交する線分lv が求められる。ここで、補間画素P0 ′の補間画像データS0 ′はこの補間画素P0 ′に隣接する4つの原画素Pa,Pb,Pc,Pdが担持する原画像データSa,Sb,Sc,Sdに前述した補間係数をそれぞれ乗じることにより算出されるが、本発明は、この4つの原画像データSa,Sb,Sc,Sdに乗じる補間係数Aa,Ab,Ac,Adを、濃度ベクトルPV の大きさおよび原画素Pa,Pb,Pc,Pdの線分lv からの距離に応じて補正する。この補正は以下のようにして行う。
【0044】
図6に示す原画素の座標値を(u,v)、座標(u,v)における補間係数をAu,v 、補間画素P0 ′の座標値を(dx,dy)、濃度ベクトルPV を(Px,Py)とする。そして、下記の式(10)により補正後の補間係数Au,v ′を求める。
【0045】
Au,v ′=Au,v /(k・f(PV ,t)+1) (10)
但し、k:定数
f(PV ,t):(濃度ベクトルPV )×(座標(u,v)の直交線分
lv からの距離)(補正項)
式(10)の計算は具体的には以下のようにして行う。
【0046】
まず、補正項算出手段5において、下記の式(11)により、手段4において求められた濃度ベクトル(Px,Py)の大きさ|PV |を求める。
【0047】
|PV |=√(Px2 +Py2 ) (11)
次いで下記の式(12)により、補間画素(dx,dy)を通り濃度ベクトルPV と直交する線分lv と、原画素(u,v)との距離tを求める。
【0048】
したがって、補正項f(PV ,t)は、
f(PV ,t)=|Px(u−dx)+Py(v−dy)| (13)となる。
【0049】
そして、補正手段6において、補正項算出手段5において算出された補正項f(PV ,t)を用いて補間係数の補正を行う。すなわち、式(13)を式(10)に代入して下記の式(10′)を求める。
【0050】
そして、この補間係数Au,v の補正を補間画素P0 ′の近傍の4つの原画素Pa,Pb,Pc,Pdについて行い、4つの補正補間係数Au,v ′を得、この求められた補正補間係数Au,v ′を正規化して最終的な補間係数
【0051】
【数1】
【0052】
を得る。
【0053】
このようにして得られた最終的な補間係数は、補間画像データS0 ′を求めるべき補間画素P0 ′における濃度ベクトルPV の傾きが大きく、かつ濃度ベクトルPV と直交する線分との距離が大きいほど、元の補間係数と比較して値が小さくなるように補正される。すなわち、原画像において濃度ベクトルPV が比較的大きくなるのは、濃度変化が急峻なエッジ部であり、この場合、濃度ベクトルPV に垂直な線分は原画像におけるエッジ部を示すものとなる。本発明においては、濃度ベクトルPV が比較的大きいエッジ部が、原画像の正方形格子状の画素に対して斜めに延びている場合において、原画像を拡大するための補間演算を行った際に、この拡大されたエッジ部に段付き形状が現れることを防止するためになされたものである。
【0054】
すなわち、濃度ベクトルPV の大きさが比較的大きい場合に、この濃度ベクトルPV に直交する線分lv から、離れた位置にある原画素の補間係数をそのまま補間演算に用いることは、エッジから離れた位置にある原画像データを補間演算に用いることとなるため、補間画像データS0 ′はエッジとは異なる濃度を有する画素に影響されてしまい、補間画像データS0 ′がエッジ部の濃度とは異なるものとなる。このため得られた補間画像データS0 ′に基づいて画像を再生した場合、斜めのエッジ部が原画像と同様の階段状の段付きが拡大された画像となってしまう。
【0055】
このため、本発明においては、算出された補間係数をそのまま補間画素P0 ′近傍の原画素における原画像データに乗じるのではなく、補間画素P0 ′における濃度ベクトルPV に直交する線分lv からの距離に応じて補正して原画像データに乗じるようにしたものである。すなわち、濃度ベクトルPV の大きさが大きいほど、また濃度ベクトルPV に直交する線分lv からの距離が大きいほど補間係数が小さくなるように補正するようにしたものである。したがって、エッジが斜めに延びている場合に、エッジ部にある補間画像データS0 ′はエッジから離れた位置にある原画像データに依存する割合が低くなり、エッジ近傍すなわちエッジに沿った方向にある原画像データに依存する割合が高くなる。これにより、エッジ部にある補間画素P0 ′はエッジ部と略同様の濃度に補間されることとなるため、原画像の階段状の段付きが拡大されることが無くなり、斜め方向に延びるエッジを階段状の段付きが生じることなく拡大することができる。
【0056】
なお、補間画素P0 ′がエッジ部でない場合は、濃度ベクトルPV の大きさは比較的小さいため、原画像データが濃度ベクトルPV に直交する線分lv から比較的離れた位置にあっても、補間係数の補正のされ方はそれほど大きくならない。したがって、補間画像データS0 ′は従来の補間係数とほとんど変わらない補間係数により求められることとなる。
【0057】
次いで、補間画像データ算出手段に7において、最終的な補間係数を用いて、補間画素P0 ′における補間画像データS0 ′を算出する。この補間画像データS0 ′の算出は上述した式(9)にしたがって行われる。すなわち、
但し、 Aa′,Ab′,Ac′,Ad′:補正補間係数
そして式(9′)に基づく演算を全ての補間画素について行い、拡大された原画像に対応する補間画像データS0 ′を得る。
【0058】
このようにして得られた補間画像データS0 ′は、CRTなどの画像出力手段8に入力され、可視像として表示される。画像出力手段8において再生された画像は、エッジが斜め方向に延びている部分についても、階段状の段付きが生じることが無く、滑らかでかつシャープなエッジとなる。
【0059】
なお、上述した実施の形態においては、式(10)にしたがって補間係数を補正するようにしているが、下記の式(14)にしたがって、補間係数を補正するようにしてもよい。
【0060】
Au,v ′=Au,v /(F(PV ,t)) (14)
ここで、F(PV ,t)は、図7に示すようにf(PV ,t)が増加するにしたがって、単調増加する関数である。
【0061】
また、上記実施の形態においては、図2に示すような放射線画像読取装置において得られた画像に対して補間演算を行うようにしているが、例えば予め記憶手段に記憶された原画像データに対して補間演算を行うようにしてもよいものである。
【0062】
さらに、上述した実施の形態においては、Cubic スプライン補間により補間演算を行うようにしているが、これに限定されるものではなく、Bスプライン補間により補間演算を行うようにしてもよい。このBスプライン補間は、Cubic スプライン補間とは異なり、元のサンプル点(画素)を通ることは必要とされない代わりに、第1階微分係数および第2階微分係数が各区間間で連続することが必要とされ、この条件の下に上記と同様の演算により、原画像データYk−1 ,Yk ,Yk+1 ,Yk+2 のBスプライン補間演算の補間係数bk−1 ,bk ,bk+1 ,bk+2 を求めればよい。
【0063】
また、原画像データに対してBスプライン補間およびCubic スプライン補間の双方を用いて、拡大率に応じてそれぞれの補間係数に重み付けをして補間演算を行うようにしてもよい。すなわち、下記の式(15)
F=t・A+(1−t)・B (15)
但し、F:補間値
A,B:補間係数
t:重み係数
において、AをCubic スプライン補間の補間係数、BをBスプライン補間の補間係数とし、重み係数tの値を変更して補間演算を行う。例えば、原画像データYk−1 ,Yk ,Yk+1 ,Yk+2 に対するCubic スプライン補間係数をak−1 ,ak ,ak+1 ,ak+2 、Bスプライン補間係数をbk−1 ,bk ,bk+1 ,bk+2 としたとき、補間値Fは、
F={t・ak−1 +(1−t)・bk−1 }Yk−1
+{t・ak +(1−t)・bk }Yk
+{t・ak+1 +(1−t)・bk+1 }Yk+1
+{t・ak+2 +(1−t)・bk+2 }Yk+2 (16)
となる。
【0064】
また、上記実施の形態においては、濃度ベクトルの大きさおよび濃度ベクトルに直交する線分と各近傍原画素との距離の大きさの双方に基づいて補正を行うようにしているが、画像におけるエッジ部とエッジ部でない部分を予め判別しておき、エッジ部の場合にのみ本発明による画像データ補間演算を行う場合には、濃度ベクトルの大きさ|PV |を固定値として、演算を簡略化するようにしてもよい。
【図面の簡単な説明】
【図1】本発明による画像データ補間演算装置の概略を示すブロック図
【図2】放射線画像読取装置を表す図
【図3】等間隔の周期でサンプリングされた一方向に配列されたサンプリング点(画素)の原画像データから補間画像データを求める方法を説明するグラフ
【図4】濃度ベクトルの算出を説明するための図
【図5】濃度ベクトルの算出を説明するための図
【図6】濃度ベクトル、濃度ベクトルに直交する線分および原画素のこの線分からの距離を表す図
【図7】補正項の他の例を表す図
【図8】原画像データを構成する画素と補間画像データを構成する画素とを示す図
【図9】原画像および原画像データを構成する画素を表す図
【図10】補間画像および補間画像データを構成する画素を表す図
【符号の説明】
1 画像データ補間演算装置
2 画像入力手段
3 補間係数算出手段
4 濃度ベクトル算出および直交線分との距離を算出する手段
5 補正項算出手段
6 補正手段
7 補間画像データ算出手段
8 画像出力手段
10 蓄積性蛍光体シート
11 シート搬送手段
12 励起光源
13 レーザビーム
14 回転多面鏡
15 輝尽発光光
16 集光体
17 フォトマルチプライヤー
Claims (4)
- 所定の間隔で縦横方向に格子状に配された原画素のそれぞれについての画素値を表す多数の原画像データについて、該原画素とは間隔が異なる補間画素における補間画像データを求める画像データの補間演算方法であって、該補間画素近傍の複数点の近傍原画素における近傍原画像データに基づいて該各近傍原画像データに乗じる補間係数をそれぞれ算出するとともに、該各補間係数を該近傍原画像データにそれぞれ乗じる補間演算処理を行うことにより、前記補間画像データを算出する画像データ補間演算方法において、
前記原画像データにより表される原画像上の前記補間画素における濃度ベクトルを算出し、
前記補間画素を通り、前記濃度ベクトルに直交する線分と、前記各近傍原画素との距離を算出し、
前記濃度ベクトルの大きさが大きいほど、かつ前記距離が大きいほど、前記各近傍画素における前記近傍原画像データに乗じる補間係数が小さくなるように該補間係数を補正し、
該補正された補間係数に基づいて前記補間演算処理を行うことを特徴とする画像データ補間演算方法。 - 所定の間隔で縦横方向に格子状に配された原画素のそれぞれについての画素値を表す多数の原画像データについて、該原画素とは間隔が異なる補間画素における補間画像データを求める画像データの補間演算方法であって、該補間画素近傍の複数点の近傍原画素における近傍原画像データに基づいて該各近傍原画像データに乗じる補間係数をそれぞれ算出するとともに、該各補間係数を該近傍原画像データにそれぞれ乗じる補間演算処理を行うことにより、前記補間画像データを算出する画像データ補間演算方法において、
前記原画像データにより表される原画像上の前記補間画素における濃度ベクトルを算出し、
前記補間画素を通り、前記濃度ベクトルに直交する線分と、前記各近傍原画素との距離を算出し、
前記原画像データにより表される原画像において、エッジ部とエッジ部以外とを判別し、
前記エッジ部においては、前記距離が大きいほど、前記各近傍画素における前記近傍原画像データに乗じる補間係数が小さくなるように該補間係数を補正し、
該補正された補間係数に基づいて前記補間演算処理を行うことを特徴とする画像データ補間演算方法。 - 所定の間隔で縦横方向に格子状に配された原画素のそれぞれについての画素値を表す多数の原画像データについて、該原画素とは間隔が異なる補間画素における補間画像データを求める画像データの補間演算装置であって、該補間画素近傍の複数点の近傍原画素における近傍原画像データに基づいて該各近傍原画像データに乗じる補間係数をそれぞれ算出するとともに、該各補間係数を該近傍原画像データにそれぞれ乗じる補間演算処理を行うことにより、前記補間画像データを算出する補間演算処理手段を備えた画像データ補間演算装置において、
前記原画像データにより表される原画像上の前記補間画素における濃度ベクトルを算出する濃度ベクトル算出手段と、
前記補間画素を通り、前記濃度ベクトルに直交する線分と、前記各近傍原画素との距離を算出する距離算出手段と、
前記濃度ベクトルの大きさが大きいほど、かつ前記距離が大きいほど、前記各近傍画素における前記近傍原画像データに乗じる補間係数が小さくなるよう該補間係数を補正する補正手段とをさらに備え、
前記補間演算処理手段が、該補正された補間係数に基づいて前記補間画像データを算出する手段であることを特徴とする画像データ補間演算装置。 - 所定の間隔で縦横方向に格子状に配された原画素のそれぞれについての画素値を表す多数の原画像データについて、該原画素とは間隔が異なる補間画素における補間画像データを求める画像データの補間演算装置であって、該補間画素近傍の複数点の近傍原画素における近傍原画像データに基づいて該各近傍原画像データに乗じる補間係数をそれぞれ算出するとともに、該各補間係数を該近傍原画像データにそれぞれ乗じる補間演算処理を行うことにより、前記補間画像データを算出する補間演算処理手段を備えた画像データ補間演算装置において、
前記原画像データにより表される原画像上の前記補間画素における濃度ベクトルを算出する濃度ベクトル算出手段と、
前記補間画素を通り、前記濃度ベクトルに直交する線分と、前記各近傍原画素との距離を算出する距離算出手段と、
前記原画像データにより表される原画像において、エッジ部とエッジ部以外とを判別する判別手段と、
前記エッジ部においては、前記距離が大きいほど、前記各近傍画素における前記近傍原画像データに乗じる補間係数が小さくなるように該補間係数を補正する補正手段とをさらに備え、
前記補間演算処理手段が、該補正された補間係数に基づいて前記補間画像データを算出する手段であることを特徴とする画像データ補間演算装置。
Priority Applications (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP20320795A JP3599435B2 (ja) | 1995-08-09 | 1995-08-09 | 画像データ補間演算方法および装置 |
US08/686,030 US5737101A (en) | 1995-08-09 | 1996-07-25 | Interpolating operation method and apparatus for image signals |
DE69618398T DE69618398T2 (de) | 1995-08-09 | 1996-08-08 | Verfahren und Vorrichtung zum Interpolieren von Bildsignalen |
EP96112833A EP0758179B1 (en) | 1995-08-09 | 1996-08-08 | Interpolating operation method and apparatus for image signals |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP20320795A JP3599435B2 (ja) | 1995-08-09 | 1995-08-09 | 画像データ補間演算方法および装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
JPH0951427A JPH0951427A (ja) | 1997-02-18 |
JP3599435B2 true JP3599435B2 (ja) | 2004-12-08 |
Family
ID=16470247
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP20320795A Expired - Fee Related JP3599435B2 (ja) | 1995-08-09 | 1995-08-09 | 画像データ補間演算方法および装置 |
Country Status (4)
Country | Link |
---|---|
US (1) | US5737101A (ja) |
EP (1) | EP0758179B1 (ja) |
JP (1) | JP3599435B2 (ja) |
DE (1) | DE69618398T2 (ja) |
Families Citing this family (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH08185149A (ja) * | 1995-01-05 | 1996-07-16 | Fujitsu Ltd | 画像処理における階調色獲得方法および装置 |
JP3706189B2 (ja) * | 1996-03-15 | 2005-10-12 | 富士写真フイルム株式会社 | 画像データの補間処理方法およびその方法を使用した装置 |
JP3753197B2 (ja) | 1996-03-28 | 2006-03-08 | 富士写真フイルム株式会社 | 画像データの補間演算方法およびその方法を実施する装置 |
JP3763901B2 (ja) * | 1996-10-17 | 2006-04-05 | 富士写真フイルム株式会社 | 画像情報読取装置 |
US6631216B2 (en) * | 1997-07-31 | 2003-10-07 | Canon Kabushiki Kaisha | Image interpolation apparatus |
EP0921494B1 (en) * | 1997-12-03 | 2003-10-15 | Fuji Photo Film Co., Ltd. | Interpolating operation method and apparatus for image signals |
US6252576B1 (en) | 1998-08-06 | 2001-06-26 | In-System Design, Inc. | Hardware-efficient system for hybrid-bilinear image scaling |
JP2001008037A (ja) * | 1999-06-23 | 2001-01-12 | Nec Corp | 画素補間方法および回路 |
JP2001013946A (ja) * | 1999-06-28 | 2001-01-19 | Alps Electric Co Ltd | 画像表示装置 |
KR100658660B1 (ko) * | 2000-08-23 | 2006-12-15 | 삼성전자주식회사 | 큐빅 컨벌루션 보간계수 생성방법 |
US6747630B2 (en) | 2002-07-31 | 2004-06-08 | Texas Instruments Incorporated | Method to up-sample frequency rich images without significant loss of image sharpness |
JP4144378B2 (ja) | 2003-02-28 | 2008-09-03 | ソニー株式会社 | 画像処理装置および方法、記録媒体、並びにプログラム |
JP4703504B2 (ja) | 2006-07-21 | 2011-06-15 | ソニー株式会社 | 画像処理装置、画像処理方法、及びプログラム |
US7800694B2 (en) * | 2006-08-31 | 2010-09-21 | Microsoft Corporation | Modular grid display |
KR101544914B1 (ko) * | 2009-02-12 | 2015-08-17 | 삼성전자주식회사 | 보간 픽셀을 생성하는 방법 |
US8264611B2 (en) * | 2010-03-27 | 2012-09-11 | Newport Media, Inc. | System and method for video downscaling |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5054100A (en) * | 1989-11-16 | 1991-10-01 | Eastman Kodak Company | Pixel interpolator with edge sharpening |
US5131057A (en) * | 1990-02-12 | 1992-07-14 | Wright State University | Method for video-to-printing image resolution conversion |
EP0517912A1 (en) * | 1990-12-31 | 1992-12-16 | E.I. Du Pont De Nemours And Company | Method of resizing an image, designing a filter therefor, and mapping the output image to the input image |
JP3143209B2 (ja) * | 1992-05-25 | 2001-03-07 | 株式会社リコー | 画像変倍装置 |
-
1995
- 1995-08-09 JP JP20320795A patent/JP3599435B2/ja not_active Expired - Fee Related
-
1996
- 1996-07-25 US US08/686,030 patent/US5737101A/en not_active Expired - Lifetime
- 1996-08-08 DE DE69618398T patent/DE69618398T2/de not_active Expired - Lifetime
- 1996-08-08 EP EP96112833A patent/EP0758179B1/en not_active Expired - Lifetime
Also Published As
Publication number | Publication date |
---|---|
DE69618398D1 (de) | 2002-02-14 |
JPH0951427A (ja) | 1997-02-18 |
US5737101A (en) | 1998-04-07 |
EP0758179A2 (en) | 1997-02-12 |
EP0758179A3 (en) | 1997-07-30 |
DE69618398T2 (de) | 2002-08-14 |
EP0758179B1 (en) | 2002-01-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP3599435B2 (ja) | 画像データ補間演算方法および装置 | |
EP0170270B1 (en) | Method of adjusting radiation image read-out conditions | |
EP0440166B1 (en) | Method for compressing dynamic ranges of images | |
EP0182423B1 (en) | Radiation image processing apparatus and method | |
JP3196864B2 (ja) | 画像のダイナミックレンジ圧縮処理方法 | |
JP2849964B2 (ja) | 画像処理方法および装置 | |
JPH06304159A (ja) | 画像の撮影***認識方法 | |
US5796865A (en) | Gradation correcting method and apparatus | |
JP2754068B2 (ja) | 放射線画像信号生成方法および放射線画像読取装置 | |
US5832142A (en) | Interpolating operation method and apparatus for image signals | |
EP1202555A2 (en) | Image signal generating method, apparatus and program | |
US5726766A (en) | Interpolating operation method and apparatus for image signals | |
JPH0950516A (ja) | 画像データ補間演算方法および装置 | |
JP3720081B2 (ja) | 階調補正方法および装置 | |
JP3604201B2 (ja) | 画像データのスプライン補間演算方法およびその方法を使用した装置 | |
JP3399684B2 (ja) | 画像処理方法および装置 | |
JPH0951429A (ja) | 画像データ補間演算方法および装置 | |
JPH09160149A (ja) | 画像処理方法および装置 | |
JPH0951430A (ja) | 画像データ補間演算方法および装置 | |
JPH04154387A (ja) | 放射線画像のエネルギーサブトラクション方法および装置 | |
JP3241488B2 (ja) | 画像処理方法 | |
JPH01194683A (ja) | 放射線画像処理方法 | |
JPH0786925B2 (ja) | 放射線画像処理方法および装置 | |
JPH02272681A (ja) | 放射線画像処理装置 | |
JPH0816767A (ja) | 画像拡大縮小方法および装置 |
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: 20040914 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20040914 |
|
R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313111 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20070924 Year of fee payment: 3 |
|
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: 20080924 Year of fee payment: 4 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20080924 Year of fee payment: 4 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090924 Year of fee payment: 5 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090924 Year of fee payment: 5 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100924 Year of fee payment: 6 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100924 Year of fee payment: 6 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110924 Year of fee payment: 7 |
|
LAPS | Cancellation because of no payment of annual fees |