JPH0927901A - 画像データ補間演算方法および装置 - Google Patents
画像データ補間演算方法および装置Info
- Publication number
- JPH0927901A JPH0927901A JP7177014A JP17701495A JPH0927901A JP H0927901 A JPH0927901 A JP H0927901A JP 7177014 A JP7177014 A JP 7177014A JP 17701495 A JP17701495 A JP 17701495A JP H0927901 A JPH0927901 A JP H0927901A
- Authority
- JP
- Japan
- Prior art keywords
- image data
- directions
- image
- interpolation
- interpolation calculation
- 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.)
- Withdrawn
Links
- 238000000034 method Methods 0.000 title claims abstract description 20
- 230000004044 response Effects 0.000 claims abstract description 41
- 238000004364 calculation method Methods 0.000 claims description 62
- 230000008569 process Effects 0.000 claims description 10
- 230000005855 radiation Effects 0.000 description 15
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 description 14
- 230000006870 function Effects 0.000 description 12
- 230000005284 excitation Effects 0.000 description 5
- 230000008859 change Effects 0.000 description 2
- 238000001914 filtration Methods 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
- 230000005540 biological transmission Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000012888 cubic function Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000010894 electron beam technology Methods 0.000 description 1
- 230000001747 exhibiting effect Effects 0.000 description 1
- 230000001678 irradiating effect Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000000465 moulding Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
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/4007—Scaling of whole images or parts thereof, e.g. expanding or contracting based on interpolation, e.g. bilinear interpolation
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Radiography Using Non-Light Waves (AREA)
- Image Processing (AREA)
- Editing Of Facsimile Originals (AREA)
Abstract
(57)【要約】
【課題】 原画像データを所定の拡大率により拡大する
際に行う補間演算処理において、画像再生装置において
再生される画像のＸ、Ｙ方向の鮮鋭度を略一定に揃え
る。 【解決手段】 所定の間隔でＸ、Ｙ方向に格子状に配さ
れた画素のそれぞれについての画素値を示す原画像デー
タＤに対して補間演算処理を施し、原画像データＤとは
間隔が異なる補間画像データを求め、補間画像データ
を、Ｘ、Ｙ走査により画像を再生する感熱記録装置28に
送り、ハードコピー29として再生する。この際、感熱記
録装置28のＸ、Ｙ方向のレスポンスの差に応じて、Ｘ、
Ｙ方向についてCubic スプライン、Ｂスプラインなどの
補間方法により、互いに異なる補間演算処理を行う。こ
れにより感熱記録装置28におけるＸ、Ｙ方向の最終的な
レスポンスが互いに略等しく揃えられる。
際に行う補間演算処理において、画像再生装置において
再生される画像のＸ、Ｙ方向の鮮鋭度を略一定に揃え
る。 【解決手段】 所定の間隔でＸ、Ｙ方向に格子状に配さ
れた画素のそれぞれについての画素値を示す原画像デー
タＤに対して補間演算処理を施し、原画像データＤとは
間隔が異なる補間画像データを求め、補間画像データ
を、Ｘ、Ｙ走査により画像を再生する感熱記録装置28に
送り、ハードコピー29として再生する。この際、感熱記
録装置28のＸ、Ｙ方向のレスポンスの差に応じて、Ｘ、
Ｙ方向についてCubic スプライン、Ｂスプラインなどの
補間方法により、互いに異なる補間演算処理を行う。こ
れにより感熱記録装置28におけるＸ、Ｙ方向の最終的な
レスポンスが互いに略等しく揃えられる。
Description
【０００１】
【発明の属する技術分野】本発明は画像データの補間演
算方法および装置に関し、とくに詳細には再生画像の
Ｘ、Ｙ方向のレスポンスを互いに略等しく揃えられるよ
うにした画像データの補間演算方法および装置に関する
ものである。
算方法および装置に関し、とくに詳細には再生画像の
Ｘ、Ｙ方向のレスポンスを互いに略等しく揃えられるよ
うにした画像データの補間演算方法および装置に関する
ものである。
【０００２】
【従来の技術】ある種の蛍光体に放射線（Ｘ線、α線、
β線、γ線、電子線、紫外線等）を照射すると、この放
射線のエネルギーの一部がその蛍光体中に蓄積され、そ
の後該蛍光体に可視光等の励起光を照射すると、蓄積さ
れたエネルギーに応じて蛍光体が輝尽発光を示す。この
ような性質を示す蛍光体を蓄積性蛍光体と言う。
β線、γ線、電子線、紫外線等）を照射すると、この放
射線のエネルギーの一部がその蛍光体中に蓄積され、そ
の後該蛍光体に可視光等の励起光を照射すると、蓄積さ
れたエネルギーに応じて蛍光体が輝尽発光を示す。この
ような性質を示す蛍光体を蓄積性蛍光体と言う。
【０００３】この蓄積性蛍光体を利用して、人体等の被
写体の放射線画像情報を一旦蓄積性蛍光体のシート（以
下、蓄積性蛍光体シートと称する）に記録し、該シート
を励起光で走査して輝尽発光させ、この輝尽発光光を光
電的に読み取って画像信号を得、この画像信号を処理し
て診断適性の良い被写体の放射線画像を得る放射線画像
情報記録再生システムが提案されている（例えば特開昭
55−12429 号、同55−116340号、同55−163472号、同56
−11395 号、同56−104645号など）。
写体の放射線画像情報を一旦蓄積性蛍光体のシート（以
下、蓄積性蛍光体シートと称する）に記録し、該シート
を励起光で走査して輝尽発光させ、この輝尽発光光を光
電的に読み取って画像信号を得、この画像信号を処理し
て診断適性の良い被写体の放射線画像を得る放射線画像
情報記録再生システムが提案されている（例えば特開昭
55−12429 号、同55−116340号、同55−163472号、同56
−11395 号、同56−104645号など）。
【０００４】上記のように画像信号を得てこの画像信号
に基づいて可視画像を再生するシステムにおいて、その
可視画像のうち観察対象となる関心領域をより詳細に観
察したいとき、その領域を拡大して再生することがあ
る。この可視画像は、原画像を読み取って得られた原画
像データに対して所定の補間演算処理を施して原画像デ
ータ数とは異なるデータ数の２次元的な画像データであ
る補間画像データを求め、この補間画像データに基づい
た可視画像の再生を行うことによって得ることができ
る。
に基づいて可視画像を再生するシステムにおいて、その
可視画像のうち観察対象となる関心領域をより詳細に観
察したいとき、その領域を拡大して再生することがあ
る。この可視画像は、原画像を読み取って得られた原画
像データに対して所定の補間演算処理を施して原画像デ
ータ数とは異なるデータ数の２次元的な画像データであ
る補間画像データを求め、この補間画像データに基づい
た可視画像の再生を行うことによって得ることができ
る。
【０００５】このような画像データの補間法としては、
線形補間による方法の他、２次あるいは３次のスプライ
ン補間関数を用いる方法など種々の方法が提案されてい
る。例えば、３次のスプライン補間関数を用いるCubic
スプライン補間演算は、元のサンプル点（画素）を通る
ことと、その第１階微分係数が各区間間で連続すること
が必要とされており、比較的鮮鋭度の高いシャープな２
次画像（補間により得られる画像）を再生するためのも
のである。また、Cubic スプライン補間演算に対して比
較的鮮鋭度の低い滑らかな２次画像を再生するための補
間画像データを得るＢスプライン補間演算も知られてい
る。このように２次画像を高い鮮鋭度でシャープに再生
したい場合は、Cubic スプライン補間演算を用い、低い
鮮鋭度で滑らかに再生したい場合はＢスプライン補間演
算を用いればよい。
線形補間による方法の他、２次あるいは３次のスプライ
ン補間関数を用いる方法など種々の方法が提案されてい
る。例えば、３次のスプライン補間関数を用いるCubic
スプライン補間演算は、元のサンプル点（画素）を通る
ことと、その第１階微分係数が各区間間で連続すること
が必要とされており、比較的鮮鋭度の高いシャープな２
次画像（補間により得られる画像）を再生するためのも
のである。また、Cubic スプライン補間演算に対して比
較的鮮鋭度の低い滑らかな２次画像を再生するための補
間画像データを得るＢスプライン補間演算も知られてい
る。このように２次画像を高い鮮鋭度でシャープに再生
したい場合は、Cubic スプライン補間演算を用い、低い
鮮鋭度で滑らかに再生したい場合はＢスプライン補間演
算を用いればよい。
【０００６】
【発明が解決しようとする課題】しかしながら、補間画
像データを再生するための装置が、Ｘ、Ｙ走査により画
像を再生するもの、例えばＣＲＴ表示装置、サーマルヘ
ッドを用いる感熱記録装置、光走査記録装置等である
と、Ｘ、Ｙ方向のレスポンスが互いに異なり、方向によ
って被写体の見え方が変わってしまうという不具合が従
来より認められていた。この不具合は、画像再生装置の
Ｘ、Ｙ方向のレスポンスが互いに異なっていることに起
因していると考えられる。
像データを再生するための装置が、Ｘ、Ｙ走査により画
像を再生するもの、例えばＣＲＴ表示装置、サーマルヘ
ッドを用いる感熱記録装置、光走査記録装置等である
と、Ｘ、Ｙ方向のレスポンスが互いに異なり、方向によ
って被写体の見え方が変わってしまうという不具合が従
来より認められていた。この不具合は、画像再生装置の
Ｘ、Ｙ方向のレスポンスが互いに異なっていることに起
因していると考えられる。
【０００７】この場合、再生画像の鮮鋭度を強調して、
レスポンスが低い方の鮮鋭度を十分となるようにするこ
とも可能であるが、そのようにすると、レスポンスが高
い方向について過度の周波数強調がなされてアーチファ
クト（偽画像）が発生することがある。
レスポンスが低い方の鮮鋭度を十分となるようにするこ
とも可能であるが、そのようにすると、レスポンスが高
い方向について過度の周波数強調がなされてアーチファ
クト（偽画像）が発生することがある。
【０００８】本発明は上記の事情に鑑み、補間演算処理
を行った上で、画像再生装置において再生される画像の
Ｘ、Ｙ方向のレスポンスを略一定に揃えることができる
画像データの補間演算方法および装置を提供することを
目的とするものである。
を行った上で、画像再生装置において再生される画像の
Ｘ、Ｙ方向のレスポンスを略一定に揃えることができる
画像データの補間演算方法および装置を提供することを
目的とするものである。
【０００９】
【課題を解決するための手段】本発明による画像データ
補間演算方法および装置は、所定の間隔でＸ、Ｙ方向に
格子状に配された画素のそれぞれについての画素値を示
す原画像データについて、補間演算処理により該原画像
データとは間隔が異なる補間画像データを求め、該補間
画像データを、Ｘ、Ｙ走査により画像を再生する画像再
生装置に送る画像データの補間演算方法および装置にお
いて、前記画像再生装置のＸ、Ｙ方向のレスポンスの差
に応じて、該Ｘ、Ｙ方向について互いに異なる補間演算
処理を行うことにより、前記画像再生装置において再生
される画像のＸ、Ｙ方向の最終的なレスポンスを互いに
略等しく揃えることを特徴とするものである。
補間演算方法および装置は、所定の間隔でＸ、Ｙ方向に
格子状に配された画素のそれぞれについての画素値を示
す原画像データについて、補間演算処理により該原画像
データとは間隔が異なる補間画像データを求め、該補間
画像データを、Ｘ、Ｙ走査により画像を再生する画像再
生装置に送る画像データの補間演算方法および装置にお
いて、前記画像再生装置のＸ、Ｙ方向のレスポンスの差
に応じて、該Ｘ、Ｙ方向について互いに異なる補間演算
処理を行うことにより、前記画像再生装置において再生
される画像のＸ、Ｙ方向の最終的なレスポンスを互いに
略等しく揃えることを特徴とするものである。
【００１０】
【作用および発明の効果】本発明による補間画像データ
演算方法および装置は、画像再生装置のＸ、Ｙ方向につ
いて互いに異なる補間演算を行うことにより、画像再生
装置のＸ、Ｙ方向のレスポンスの差を補償するようにし
たため、再生画像のＸ、Ｙ方向の最終的なレスポンスを
互いに略等しく揃えることができる。これにより、Ｘ、
Ｙ方向のレスポンスのバランスが異なることを防止する
ことができる。また、Ｘ、Ｙ方向のうちの一方向につい
て過度の周波数強調処理をかけてアーチファクトを発生
させてしまうことも防止することができる。
演算方法および装置は、画像再生装置のＸ、Ｙ方向につ
いて互いに異なる補間演算を行うことにより、画像再生
装置のＸ、Ｙ方向のレスポンスの差を補償するようにし
たため、再生画像のＸ、Ｙ方向の最終的なレスポンスを
互いに略等しく揃えることができる。これにより、Ｘ、
Ｙ方向のレスポンスのバランスが異なることを防止する
ことができる。また、Ｘ、Ｙ方向のうちの一方向につい
て過度の周波数強調処理をかけてアーチファクトを発生
させてしまうことも防止することができる。
【００１１】なお本発明は、前述した蓄積性蛍光体シー
トを用いる放射線画像情報記録再生システムにおいて適
用されれば、診断適性が極めて高い放射線画像を再生で
きるので特に好ましいが、このような放射線画像情報記
録再生システム以外に適用された場合でも、上記の通り
の作用、効果を奏することは勿論である。
トを用いる放射線画像情報記録再生システムにおいて適
用されれば、診断適性が極めて高い放射線画像を再生で
きるので特に好ましいが、このような放射線画像情報記
録再生システム以外に適用された場合でも、上記の通り
の作用、効果を奏することは勿論である。
【００１２】
【発明の実施の形態】以下、図面に示す実施の形態に基
づいて本発明を詳細に説明する。図１は本発明の補間演
算を行う装置を内包した放射線画像情報記録再生システ
ムの一例を示すものである。このシステムは、前述の蓄
積性蛍光体シートに記録された放射線画像情報を読み取
り、その読み取った放射線画像を一例として感熱記録装
置により再生記録するものである。
づいて本発明を詳細に説明する。図１は本発明の補間演
算を行う装置を内包した放射線画像情報記録再生システ
ムの一例を示すものである。このシステムは、前述の蓄
積性蛍光体シートに記録された放射線画像情報を読み取
り、その読み取った放射線画像を一例として感熱記録装
置により再生記録するものである。
【００１３】例えばＸ線等の放射線が人体等の被写体を
介して照射されることによりこの被写体の透過放射線画
像情報を蓄積記録した蓄積性蛍光体シート10は、エンド
レスベルト等のシート搬送手段11により、副走査のため
に矢印Ｙ方向に搬送される。半導体レーザ等の励起光源
12から射出された励起光（読取光）としてのレーザビー
ム13は、高速回転する回転多面鏡14によって反射偏向さ
れ、通常ｆ・θレンズからなる走査レンズ18によって集
束され、ミラー19で反射して蓄積性蛍光体シート10上を
上記副走査方向Ｙと略直角な矢印Ｘ方向に主走査する。
介して照射されることによりこの被写体の透過放射線画
像情報を蓄積記録した蓄積性蛍光体シート10は、エンド
レスベルト等のシート搬送手段11により、副走査のため
に矢印Ｙ方向に搬送される。半導体レーザ等の励起光源
12から射出された励起光（読取光）としてのレーザビー
ム13は、高速回転する回転多面鏡14によって反射偏向さ
れ、通常ｆ・θレンズからなる走査レンズ18によって集
束され、ミラー19で反射して蓄積性蛍光体シート10上を
上記副走査方向Ｙと略直角な矢印Ｘ方向に主走査する。
【００１４】こうしてレーザビーム13が照射されたシー
ト10の箇所からは、蓄積記録されている放射線画像情報
に応じた光量の輝尽発光光15が発散され、この輝尽発光
光15は集光体16によって集光され、光検出器としてのフ
ォトマルチプライヤー（光電子増倍管）17によって光電
的に検出される。
ト10の箇所からは、蓄積記録されている放射線画像情報
に応じた光量の輝尽発光光15が発散され、この輝尽発光
光15は集光体16によって集光され、光検出器としてのフ
ォトマルチプライヤー（光電子増倍管）17によって光電
的に検出される。
【００１５】上記集光体16はアクリル板等の導光性材料
を成形して作られたものであり、直線状をなす入射端面
16ａが蓄積性蛍光体シート10上のビーム走査線に沿って
延びるように配され、円環状に形成された出射端面16ｂ
に上記フォトマルチプライヤー17の受光面が結合されて
いる。上記入射端面16ａから集光体16内に入射した輝尽
発光光15は、該集光体16の内部を全反射を繰り返して進
み、出射端面16ｂから出射してフォトマルチプライヤー
17に受光され、前記放射線画像情報を担持する輝尽発光
光15の光量がこのフォトマルチプライヤー17によって検
出される。
を成形して作られたものであり、直線状をなす入射端面
16ａが蓄積性蛍光体シート10上のビーム走査線に沿って
延びるように配され、円環状に形成された出射端面16ｂ
に上記フォトマルチプライヤー17の受光面が結合されて
いる。上記入射端面16ａから集光体16内に入射した輝尽
発光光15は、該集光体16の内部を全反射を繰り返して進
み、出射端面16ｂから出射してフォトマルチプライヤー
17に受光され、前記放射線画像情報を担持する輝尽発光
光15の光量がこのフォトマルチプライヤー17によって検
出される。
【００１６】フォトマルチプライヤー17のアナログ出力
信号（画像信号）Ｓは対数増幅器20によって増幅され、
Ａ／Ｄ変換器21において所定の収録スケールファクター
でデジタル化される。こうして得られた、２次元画像を
担持するデジタルの画像データＤは、画像処理装置27に
おいて所定の拡大率により拡大するために、補間演算処
理が施される。補間演算処理が施された画像データＤ
は、階調処理、周波数強調処理等の画像処理を受ける。
この画像処理を受けた画像データＤ′は、サーマルヘッ
ド28ａを用いる感熱記録装置28に入力され、感熱記録装
置28において画像データＤ′が示す２次元画像がハード
コピー29に再生記録される。
信号（画像信号）Ｓは対数増幅器20によって増幅され、
Ａ／Ｄ変換器21において所定の収録スケールファクター
でデジタル化される。こうして得られた、２次元画像を
担持するデジタルの画像データＤは、画像処理装置27に
おいて所定の拡大率により拡大するために、補間演算処
理が施される。補間演算処理が施された画像データＤ
は、階調処理、周波数強調処理等の画像処理を受ける。
この画像処理を受けた画像データＤ′は、サーマルヘッ
ド28ａを用いる感熱記録装置28に入力され、感熱記録装
置28において画像データＤ′が示す２次元画像がハード
コピー29に再生記録される。
【００１７】次に、上記画像処理装置27における補間演
算処理について詳しく説明する。この画像処理装置27
は、１画像分の画素値を示しているデジタル画像データ
を一時的に格納しておくバッファメモリ23、後述する補
間係数等を記憶しておくメモリ24、補間演算部25および
画像処理部26からなる。
算処理について詳しく説明する。この画像処理装置27
は、１画像分の画素値を示しているデジタル画像データ
を一時的に格納しておくバッファメモリ23、後述する補
間係数等を記憶しておくメモリ24、補間演算部25および
画像処理部26からなる。
【００１８】ここで、感熱記録装置28のＸ方向（主走査
方向）におけるレスポンス特性ａとＹ方向（副走査方
向）におけるレスポンス特性ｂとは、図２に示すような
ものとなる。図２において、Ｘ、Ｙ方向で主にレスポン
スを揃えたい空間周波数をαとし、この空間周波数αに
対応するＸ、Ｙ方向のＭＴＦ（Modulation Transfer Fu
nction，周波数依存特性）をそれぞれｍ、ｍ′とする
と、Ｘ、Ｙ方向のレスポンスを揃えるためには、Ｘ方向
に対してＹ方向のレスポンスがｍ／ｍ′となるような関
数の１次元補間をＸ、Ｙ方向にかける必要がある。そし
てこのために、補間演算部25においては、補間画像デー
タのレスポンスがＸ、Ｙ方向においてそれぞれ異なるよ
うな補間演算処理が画像データＤに対して施される。
方向）におけるレスポンス特性ａとＹ方向（副走査方
向）におけるレスポンス特性ｂとは、図２に示すような
ものとなる。図２において、Ｘ、Ｙ方向で主にレスポン
スを揃えたい空間周波数をαとし、この空間周波数αに
対応するＸ、Ｙ方向のＭＴＦ（Modulation Transfer Fu
nction，周波数依存特性）をそれぞれｍ、ｍ′とする
と、Ｘ、Ｙ方向のレスポンスを揃えるためには、Ｘ方向
に対してＹ方向のレスポンスがｍ／ｍ′となるような関
数の１次元補間をＸ、Ｙ方向にかける必要がある。そし
てこのために、補間演算部25においては、補間画像デー
タのレスポンスがＸ、Ｙ方向においてそれぞれ異なるよ
うな補間演算処理が画像データＤに対して施される。
【００１９】本実施の形態においては、図２に示すよう
にＸ方向のレスポンスがＹ方向のレスポンスよりも高い
ため、Ｘ方向の画像データＤに対してＢスプライン補間
演算を施して、低い鮮鋭度により滑らかな補間を行い、
Ｙ方向の画像データＤに対してCubic スプライン補間演
算により高い鮮鋭度によりシャープな補間を行う。
にＸ方向のレスポンスがＹ方向のレスポンスよりも高い
ため、Ｘ方向の画像データＤに対してＢスプライン補間
演算を施して、低い鮮鋭度により滑らかな補間を行い、
Ｙ方向の画像データＤに対してCubic スプライン補間演
算により高い鮮鋭度によりシャープな補間を行う。
【００２０】以下Cubic スプライン補間演算およびＢス
プライン補間演算について詳細に説明する。
プライン補間演算について詳細に説明する。
【００２１】原画像からデジタル的に読み取って得られ
た、連続する画素Ｘk-2 ，Ｘk-1 ，Ｘk ，Ｘk+1 ，Ｘ
k+2 ，…の画像データ（原画像データ）を図３に示すよ
うにそれぞれＹk-2 ，Ｙk-1 ，Ｙk ，Ｙk+1 ，Ｙk+2 ，
…とする。ここで、３次のスプライン補間関数は、各区
間Ｘk-2 〜Ｘk-1 ，Ｘk-1 〜Ｘk ，Ｘk 〜Ｘk+1 ，Ｘ
k+1 〜Ｘk+2 ごとにそれぞれ設定され、各区間に対応す
るスプライン補間関数をｆk-2 ，ｆk-1 ，ｆk ，
ｆk+1 ，ｆk+2 とする。この補間関数はいずれも各区間
の位置を変数とする３次関数である。
た、連続する画素Ｘk-2 ，Ｘk-1 ，Ｘk ，Ｘk+1 ，Ｘ
k+2 ，…の画像データ（原画像データ）を図３に示すよ
うにそれぞれＹk-2 ，Ｙk-1 ，Ｙk ，Ｙk+1 ，Ｙk+2 ，
…とする。ここで、３次のスプライン補間関数は、各区
間Ｘk-2 〜Ｘk-1 ，Ｘk-1 〜Ｘk ，Ｘk 〜Ｘk+1 ，Ｘ
k+1 〜Ｘk+2 ごとにそれぞれ設定され、各区間に対応す
るスプライン補間関数をｆk-2 ，ｆk-1 ，ｆk ，
ｆk+1 ，ｆk+2 とする。この補間関数はいずれも各区間
の位置を変数とする３次関数である。
【００２２】ここでまず、補間しようとする点（以下、
補間点という）Ｘp が区間Ｘk 〜Ｘk+1 の範囲にある場
合について説明する。なお、区間Ｘk 〜Ｘk+1 に対応す
るスプライン補間関数ｆk は下記式（１）で表される。
補間点という）Ｘp が区間Ｘk 〜Ｘk+1 の範囲にある場
合について説明する。なお、区間Ｘk 〜Ｘk+1 に対応す
るスプライン補間関数ｆk は下記式（１）で表される。
【００２３】 ｆk （ｘ）＝Ａk ｘ3 ＋Ｂk ｘ2 ＋Ｃk ｘ＋Ｄk （１） Cubic スプライン補間演算においては、スプライン補間
関数ｆk は元のサンプル点（画素）を通ることと、その
第１階微分係数が各区間間で連続することが必要とさ
れ、これらの条件から下記式（２）〜（５）を満たす必
要がある。
関数ｆk は元のサンプル点（画素）を通ることと、その
第１階微分係数が各区間間で連続することが必要とさ
れ、これらの条件から下記式（２）〜（５）を満たす必
要がある。
【００２４】 ｆk （Ｘk ）＝Ｙk （２） ｆk （Ｘk+1 ）＝Ｙk+1 （３） ｆk ′（Ｘk ）＝ｆk-1 ′（Ｘk ） （４） ｆk ′（Ｘk+1 ）＝ｆk+1 ′（Ｘk+1 ） （５） なお、ｆk ′は関数ｆk の第１階微分（３Ａk ｘ2 ＋２
Ｂk ｘ＋Ｃk ）を表すものである。また、Cubic スプラ
イン補間演算は厳密には第２階微分係数の連続条件も含
むものであるが、式が複雑化するため上述したように簡
略化して用いるのが一般的である。
Ｂk ｘ＋Ｃk ）を表すものである。また、Cubic スプラ
イン補間演算は厳密には第２階微分係数の連続条件も含
むものであるが、式が複雑化するため上述したように簡
略化して用いるのが一般的である。
【００２５】またCubic スプライン補間演算において
は、画素Ｘk における第１階微分係数が、その画素Ｘk
の前後の画素であるＸk-1 とＸk+1 とについて、これら
の画像データＹk-1 、Ｙk+1 の勾配（Ｙk+1 −Ｙk-1 ）
／（Ｘk+1 −Ｘk-1 ）に一致することが条件であるか
ら、下記式（６）を満たす必要がある。
は、画素Ｘk における第１階微分係数が、その画素Ｘk
の前後の画素であるＸk-1 とＸk+1 とについて、これら
の画像データＹk-1 、Ｙk+1 の勾配（Ｙk+1 −Ｙk-1 ）
／（Ｘk+1 −Ｘk-1 ）に一致することが条件であるか
ら、下記式（６）を満たす必要がある。
【００２６】 ｆk ′（Ｘk ）＝（Ｙk+1 −Ｙk-1 ）／（Ｘk+1 −Ｘk-1 ） （６） 同様に、画素Ｘk+1 における第１階微分係数が、その画
素Ｘk+1 の前後の画素であるＸk とＸk+2 とについて、
これらの画像データＹk 、Ｙk+2 の勾配（Ｙk+2 −
Ｙk ）／（Ｘk+2 −Ｘk ）に一致することが条件である
から、下記式（７）を満たす必要がある。
素Ｘk+1 の前後の画素であるＸk とＸk+2 とについて、
これらの画像データＹk 、Ｙk+2 の勾配（Ｙk+2 −
Ｙk ）／（Ｘk+2 −Ｘk ）に一致することが条件である
から、下記式（７）を満たす必要がある。
【００２７】 ｆk ′（Ｘk+1 ）＝（Ｙk+2 −Ｙk ）／（Ｘk+2 −Ｘk ） （７） ここで、各区間Ｘk-2 〜Ｘk-1 ，Ｘk-1 〜Ｘk ，Ｘk 〜
Ｘk+1 ，Ｘk+1 〜Ｘk+2 の間隔（格子間隔という）を１
とし、画素Ｘk からの画素Ｘk+1 方向への補間点Ｘp の
位置をｔ（０≦ｔ≦１）とすれば、式（１）〜（４）お
よび（６）、（７）より、 ｆk （０）＝Ｄk ＝Ｙk ｆk （１）＝Ａk ＋Ｂk ＋Ｃk ＋Ｄk ＝Ｙk+1 ｆk ′（０）＝Ｃk ＝（Ｙk+1 −Ｙk-1 ）／２ ｆk ′（１）＝３Ａk ＋２Ｂk ＋Ｃk ＝（Ｙk+2 −
Ｙk ）／２ したがって、 Ａk ＝（Ｙk+2 −３Ｙk+1 ＋３Ｙk −Ｙk-1 ）／２ Ｂk ＝（−Ｙk+2 ＋４Ｙk+1 −５Ｙk ＋２Ｙk-1 ）／２ Ｃk ＝（Ｙk+1 −Ｙk-1 ）／２ Ｄk ＝Ｙk なお、スプライン補間関数ｆk （ｘ）は上述の通り、Ｘ
＝ｔなる変数変換をしているため、 ｆk （ｘ）＝ｆk （ｔ） となる。よって、補間点Ｘp における補間画像データＹ
p は、 Ｙp ＝ｆk （ｔ）＝Ａk ｔ3 ＋Ｂk ｔ2 ＋Ｃk ｔ＋Ｄk （８） で表すことができる。ここで上記各係数Ａk ，Ｂk ，Ｃ
k ，Ｄk を式（８）に代入すると、 Ｙp ＝｛（Ｙk+2 −３Ｙk+1 ＋３Ｙk −Ｙk-1 ）／２｝ｔ3 ＋｛（−Ｙk+2 ＋４Ｙk+1 −５Ｙk ＋２Ｙk-1 ）／２｝ｔ2 ＋｛（Ｙk+1 −Ｙk-1 ）／２｝ｔ＋Ｙk となり、これを画像データＹk-1 ，Ｙk ，Ｙk+1 ，Ｙ
k+2 について整理すると、下記式（９）で表すことがで
きる。
Ｘk+1 ，Ｘk+1 〜Ｘk+2 の間隔（格子間隔という）を１
とし、画素Ｘk からの画素Ｘk+1 方向への補間点Ｘp の
位置をｔ（０≦ｔ≦１）とすれば、式（１）〜（４）お
よび（６）、（７）より、 ｆk （０）＝Ｄk ＝Ｙk ｆk （１）＝Ａk ＋Ｂk ＋Ｃk ＋Ｄk ＝Ｙk+1 ｆk ′（０）＝Ｃk ＝（Ｙk+1 −Ｙk-1 ）／２ ｆk ′（１）＝３Ａk ＋２Ｂk ＋Ｃk ＝（Ｙk+2 −
Ｙk ）／２ したがって、 Ａk ＝（Ｙk+2 −３Ｙk+1 ＋３Ｙk −Ｙk-1 ）／２ Ｂk ＝（−Ｙk+2 ＋４Ｙk+1 −５Ｙk ＋２Ｙk-1 ）／２ Ｃk ＝（Ｙk+1 −Ｙk-1 ）／２ Ｄk ＝Ｙk なお、スプライン補間関数ｆk （ｘ）は上述の通り、Ｘ
＝ｔなる変数変換をしているため、 ｆk （ｘ）＝ｆk （ｔ） となる。よって、補間点Ｘp における補間画像データＹ
p は、 Ｙp ＝ｆk （ｔ）＝Ａk ｔ3 ＋Ｂk ｔ2 ＋Ｃk ｔ＋Ｄk （８） で表すことができる。ここで上記各係数Ａk ，Ｂk ，Ｃ
k ，Ｄk を式（８）に代入すると、 Ｙp ＝｛（Ｙk+2 −３Ｙk+1 ＋３Ｙk −Ｙk-1 ）／２｝ｔ3 ＋｛（−Ｙk+2 ＋４Ｙk+1 −５Ｙk ＋２Ｙk-1 ）／２｝ｔ2 ＋｛（Ｙk+1 −Ｙk-1 ）／２｝ｔ＋Ｙk となり、これを画像データＹk-1 ，Ｙk ，Ｙk+1 ，Ｙ
k+2 について整理すると、下記式（９）で表すことがで
きる。
【００２８】 Ｙp ＝｛（−ｔ3 ＋２ｔ2 −ｔ）／２｝Ｙk-1 ＋｛（３ｔ3 −５ｔ2 ＋２）／２｝Ｙk ＋｛（−３ｔ3 ＋４ｔ2 ＋ｔ）／２｝Ｙk+1 ＋｛（ｔ3 −ｔ2 ）／２｝Ｙk+2 （９） と表すことができる。
【００２９】ここで、原画像データＹk-1 、Ｙk 、Ｙ
k+1 、Ｙk+2 の各係数を補間係数ａk- 1 、ａk 、
ａk+1 、ａk+2 と称する。すなわち、式（９）における
原画像データＹk-1 、Ｙk 、Ｙk+1 、Ｙk+2 にそれぞれ
対応する補間係数ａk-1 、ａk 、ａk+1 、ａk+2 は、 ａk-1 ＝（−ｔ3 ＋２ｔ2 −ｔ）／２ ａk ＝（３ｔ3 −５ｔ2 ＋２）／２ ａk+1 ＝（−３ｔ3 ＋４ｔ2 ＋ｔ）／２ ａk+2 ＝（ｔ3 −ｔ2 ）／２ となる。
k+1 、Ｙk+2 の各係数を補間係数ａk- 1 、ａk 、
ａk+1 、ａk+2 と称する。すなわち、式（９）における
原画像データＹk-1 、Ｙk 、Ｙk+1 、Ｙk+2 にそれぞれ
対応する補間係数ａk-1 、ａk 、ａk+1 、ａk+2 は、 ａk-1 ＝（−ｔ3 ＋２ｔ2 −ｔ）／２ ａk ＝（３ｔ3 −５ｔ2 ＋２）／２ ａk+1 ＝（−３ｔ3 ＋４ｔ2 ＋ｔ）／２ ａk+2 ＝（ｔ3 −ｔ2 ）／２ となる。
【００３０】以上の演算を各区間Ｘk-2 〜Ｘk-1 ，Ｘ
k-1 〜Ｘk ，Ｘk 〜Ｘk+1 ，Ｘk+1 〜Ｘk+2 について繰
り返すことにより、原画像データの全体について原画像
データとは間隔の異なる補間画像データを求めることが
できる。
k-1 〜Ｘk ，Ｘk 〜Ｘk+1 ，Ｘk+1 〜Ｘk+2 について繰
り返すことにより、原画像データの全体について原画像
データとは間隔の異なる補間画像データを求めることが
できる。
【００３１】ところで、上記Cubic スプライン補間演算
は、前述したように元のサンプル点（画素）を通ること
と、その第１階微分係数が各区間間で連続することが必
要とされていて、比較的鮮鋭度の高いシャープな２次画
像（補間により得られる画像）を再生するための補間画
像データを得るものであるが、一方、比較的鮮鋭度の低
い滑らかな２次画像を再生するための補間画像データを
得るＢスプライン補間演算も知られている。このＢスプ
ライン補間演算は、元のサンプル点（画素）を通ること
は必要とされない代わりに、第１階微分係数および第２
階微分係数（ｆ″（Ｘ）で表す）が各区間間で連続する
ことが必要とされる。
は、前述したように元のサンプル点（画素）を通ること
と、その第１階微分係数が各区間間で連続することが必
要とされていて、比較的鮮鋭度の高いシャープな２次画
像（補間により得られる画像）を再生するための補間画
像データを得るものであるが、一方、比較的鮮鋭度の低
い滑らかな２次画像を再生するための補間画像データを
得るＢスプライン補間演算も知られている。このＢスプ
ライン補間演算は、元のサンプル点（画素）を通ること
は必要とされない代わりに、第１階微分係数および第２
階微分係数（ｆ″（Ｘ）で表す）が各区間間で連続する
ことが必要とされる。
【００３２】すなわち、 ｆk （ｘ）＝Ａk ｘ3 ＋Ｂk ｘ2 ＋Ｃk ｘ＋Ｄk （１） において、 ｆk ′（Ｘk ）＝ｆk-1 ′（Ｘk ） （10） ｆk ′（Ｘk+1 ）＝ｆk+1 ′（Ｘk+1 ） （11） ｆk ″（Ｘk ）＝ｆk-1 ″（Ｘk ） （12） ｆk ″（Ｘk+1 ）＝ｆk+1 ″（Ｘk+1 ） （13） が条件となる。ただし、画素Ｘk における第１階微分係
数が、その画素Ｘk の前後の画素であるＸk-1 とＸk+1
とについて、これらの画像データＹk-1 、Ｙk+1の勾配
（Ｙk+1 −Ｙk-1 ）／（Ｘk+1 −Ｘk-1 ）に一致するこ
とが条件であるから、下記式（14）を満たす必要があ
る。
数が、その画素Ｘk の前後の画素であるＸk-1 とＸk+1
とについて、これらの画像データＹk-1 、Ｙk+1の勾配
（Ｙk+1 −Ｙk-1 ）／（Ｘk+1 −Ｘk-1 ）に一致するこ
とが条件であるから、下記式（14）を満たす必要があ
る。
【００３３】 ｆk ′（Ｘk ）＝（Ｙk+1 −Ｙk-1 ）／（Ｘk+1 −Ｘk-1 ） （14） 同様に、画素Ｘk+1 における第１階微分係数が、その画
素Ｘk+1 の前後の画素であるＸk とＸk+2 とについて、
これらの画像データＹk 、Ｙk+2 の勾配（Ｙk+2 −
Ｙk ）／（Ｘk+2 −Ｘk ）に一致することが条件である
から、下記式（15）を満たす必要がある。
素Ｘk+1 の前後の画素であるＸk とＸk+2 とについて、
これらの画像データＹk 、Ｙk+2 の勾配（Ｙk+2 −
Ｙk ）／（Ｘk+2 −Ｘk ）に一致することが条件である
から、下記式（15）を満たす必要がある。
【００３４】 ｆk ′（Ｘk+1 ）＝（Ｙk+2 −Ｙk ）／（Ｘk+2 −Ｘk ） （15） また関数ｆ（Ｘ）は一般に下記式（16）に示すもので近
似される。
似される。
【００３５】 ｆ（Ｘ）＝ｆ(0) ＋ｆ′(0) Ｘ＋｛ｆ″(0)/２｝Ｘ2 （16） ここで、各区間Ｘk-2 〜Ｘk-1 ，Ｘk-1 〜Ｘk ，Ｘk 〜
Ｘk+1 ，Ｘk+1 〜Ｘk+2 の間隔（格子間隔という）を１
とし、画素Ｘk からの画素Ｘk+1 方向への補間点Ｘp の
位置をｔ（０≦ｔ≦１）とすれば、式（10）〜（13）お
よび（16）より、 ｆk ′（０）＝Ｃk ＝（Ｙk+1 −Ｙk-1 ）／２ ｆk ′（１）＝３Ａk ＋２Ｂk ＋Ｃk ＝（Ｙk+2 −
Ｙk ）／２ ｆk ″（０）＝Ｙk+1 −２Ｙk ＋Ｙk-1 ＝２Ｂ したがって、 Ａk ＝（Ｙk+2 −３Ｙk+1 ＋３Ｙk −Ｙk-1 ）／６ Ｂk ＝（Ｙk+1 −２Ｙk ＋Ｙk-1 ）／２ Ｃk ＝（Ｙk+1 −Ｙk-1 ）／２ ここで、Ｄk は未知のため、 Ｄk ＝（Ｄ1 Ｙk+2 ＋Ｄ2 Ｙk+1 ＋Ｄ3 Ｙk ＋Ｄ4 Ｙ
k-1 ）／６ とおく。また、スプライン補間関数ｆk （ｘ）は上述の
通り、Ｘ＝ｔなる変数変換をしているため、 ｆk （ｘ）＝ｆk （ｔ） となる。よって、 ｆk （ｔ）＝｛（Ｙk+2 −３Ｙk+1 ＋３Ｙk −Ｙk-1 ）／６｝ｔ3 ＋｛（Ｙk+1 −２Ｙk ＋Ｙk-1 ）／２｝ｔ2 ＋｛（Ｙk+1 −Ｙk-1 ）／２｝ｔ ＋（Ｄ1 Ｙk+2 ＋Ｄ2 Ｙk+1 ＋Ｄ3 Ｙk ＋Ｄ4 Ｙk-1 ）／６ となり、これを画像データＹk-1 ，Ｙk ，Ｙk+1 ，Ｙ
k+2 について整理すると、下記式（17）で表すことがで
きる。
Ｘk+1 ，Ｘk+1 〜Ｘk+2 の間隔（格子間隔という）を１
とし、画素Ｘk からの画素Ｘk+1 方向への補間点Ｘp の
位置をｔ（０≦ｔ≦１）とすれば、式（10）〜（13）お
よび（16）より、 ｆk ′（０）＝Ｃk ＝（Ｙk+1 −Ｙk-1 ）／２ ｆk ′（１）＝３Ａk ＋２Ｂk ＋Ｃk ＝（Ｙk+2 −
Ｙk ）／２ ｆk ″（０）＝Ｙk+1 −２Ｙk ＋Ｙk-1 ＝２Ｂ したがって、 Ａk ＝（Ｙk+2 −３Ｙk+1 ＋３Ｙk −Ｙk-1 ）／６ Ｂk ＝（Ｙk+1 −２Ｙk ＋Ｙk-1 ）／２ Ｃk ＝（Ｙk+1 −Ｙk-1 ）／２ ここで、Ｄk は未知のため、 Ｄk ＝（Ｄ1 Ｙk+2 ＋Ｄ2 Ｙk+1 ＋Ｄ3 Ｙk ＋Ｄ4 Ｙ
k-1 ）／６ とおく。また、スプライン補間関数ｆk （ｘ）は上述の
通り、Ｘ＝ｔなる変数変換をしているため、 ｆk （ｘ）＝ｆk （ｔ） となる。よって、 ｆk （ｔ）＝｛（Ｙk+2 −３Ｙk+1 ＋３Ｙk −Ｙk-1 ）／６｝ｔ3 ＋｛（Ｙk+1 −２Ｙk ＋Ｙk-1 ）／２｝ｔ2 ＋｛（Ｙk+1 −Ｙk-1 ）／２｝ｔ ＋（Ｄ1 Ｙk+2 ＋Ｄ2 Ｙk+1 ＋Ｄ3 Ｙk ＋Ｄ4 Ｙk-1 ）／６ となり、これを画像データＹk-1 ，Ｙk ，Ｙk+1 ，Ｙ
k+2 について整理すると、下記式（17）で表すことがで
きる。
【００３６】 ｆk （ｔ）＝｛（−ｔ3 ＋３ｔ2 −３ｔ＋Ｄ4 ）／６｝Ｙk-1 ＋｛（３ｔ3 −６ｔ2 ＋Ｄ3 ）／６｝Ｙk ＋｛（−３ｔ3 ＋３ｔ2 ＋３ｔ＋Ｄ2 ）／６｝Ｙk+1 ＋｛（ｔ3 ＋Ｄ1 ）／６｝Ｙk+2 （17） ここで、ｔ＝１とおけば、 ｆk （１）＝｛（Ｄ4 −１）／６｝Ｙk-1 ＋｛（Ｄ3 −３）／６｝Ｙk ＋｛（Ｄ2 ＋３）／６｝Ｙk+1 ＋｛（Ｄ1 ＋１）／６｝Ｙk+2 次に区間Ｘk+1 〜Ｘk+2 についての式（17）は、 ｆk+1 （ｔ）＝｛（−ｔ3 ＋３ｔ2 −３ｔ＋Ｄ4 ）／６｝Ｙk ＋｛（３ｔ3 −６ｔ2 ＋Ｄ3 ）／６｝Ｙk+1 ＋｛（−３ｔ3 ＋３ｔ2 ＋３ｔ＋Ｄ2 ）／６｝Ｙk+2 ＋｛（ｔ3 ＋Ｄ1 ）／６｝Ｙk+3 （18） ここで、ｔ＝０とおけば、 ｆk+1 （０）＝（Ｄ4 ／６）Ｙk ＋（Ｄ3 ／６）Ｙk+1
＋（Ｄ2 ／６）Ｙk+2 ＋（Ｄ1 ／６）Ｙk+3 連続性の条件（ｆk （１）＝ｆk+1 （０））、および各
原画像データに対応する係数同士が等しいという条件に
より、Ｄ4 −１＝０，Ｄ3 −３＝Ｄ4 ，Ｄ2 ＋３＝
Ｄ3 ，Ｄ1 ＋１＝Ｄ2 ，Ｄ1 ＝０、となり、したがっ
て、 Ｄk ＝（Ｙk+1 ＋４Ｙk ＋Ｙk-1 ）／６ となる。よって、 Ｙp ＝ｆk （ｔ）＝｛（−ｔ3 ＋３ｔ2 −３ｔ＋１）／６｝Ｙk-1 ＋｛（３ｔ3 −６ｔ2 ＋４）／６｝Ｙk ＋｛（−３ｔ3 ＋３ｔ2 ＋３ｔ＋１）／６｝Ｙk+1 ＋｛ｔ3 ／６｝Ｙk+2 （19） したがって、原画像データＹk-1 、Ｙk 、Ｙk+1 、Ｙ
k+2 にそれぞれ対応する補間係数ｂk-1 、ｂk 、
ｂk+1 、ｂk+2 は、 ｂk-1 ＝（−ｔ3 ＋３ｔ2 −３ｔ＋１）／６ ｂk ＝（３ｔ3 −６ｔ2 ＋４）／６ ｂk+1 ＝（−３ｔ3 ＋３ｔ2 ＋３ｔ＋１）／６ ｂk+2 ＝ｔ3 ／６ となる。
＋（Ｄ2 ／６）Ｙk+2 ＋（Ｄ1 ／６）Ｙk+3 連続性の条件（ｆk （１）＝ｆk+1 （０））、および各
原画像データに対応する係数同士が等しいという条件に
より、Ｄ4 −１＝０，Ｄ3 −３＝Ｄ4 ，Ｄ2 ＋３＝
Ｄ3 ，Ｄ1 ＋１＝Ｄ2 ，Ｄ1 ＝０、となり、したがっ
て、 Ｄk ＝（Ｙk+1 ＋４Ｙk ＋Ｙk-1 ）／６ となる。よって、 Ｙp ＝ｆk （ｔ）＝｛（−ｔ3 ＋３ｔ2 −３ｔ＋１）／６｝Ｙk-1 ＋｛（３ｔ3 −６ｔ2 ＋４）／６｝Ｙk ＋｛（−３ｔ3 ＋３ｔ2 ＋３ｔ＋１）／６｝Ｙk+1 ＋｛ｔ3 ／６｝Ｙk+2 （19） したがって、原画像データＹk-1 、Ｙk 、Ｙk+1 、Ｙ
k+2 にそれぞれ対応する補間係数ｂk-1 、ｂk 、
ｂk+1 、ｂk+2 は、 ｂk-1 ＝（−ｔ3 ＋３ｔ2 −３ｔ＋１）／６ ｂk ＝（３ｔ3 −６ｔ2 ＋４）／６ ｂk+1 ＝（−３ｔ3 ＋３ｔ2 ＋３ｔ＋１）／６ ｂk+2 ＝ｔ3 ／６ となる。
【００３７】以上の演算を各区間Ｘk-2 〜Ｘk-1 ，Ｘ
k-1 〜Ｘk ，Ｘk 〜Ｘk+1 ，Ｘk+1 〜Ｘk+2 について繰
り返すことにより、原画像データの全体について原画像
データとは間隔の異なる補間画像データを求めることが
できる。
k-1 〜Ｘk ，Ｘk 〜Ｘk+1 ，Ｘk+1 〜Ｘk+2 について繰
り返すことにより、原画像データの全体について原画像
データとは間隔の異なる補間画像データを求めることが
できる。
【００３８】このようにして求められた補間係数
ａk-1 、ａk 、ａk+1 、ａk+2 およびｂk-1 、ｂk 、ｂ
k+1 、ｂk+2 は、メモリ24に記憶される。そして、画像
データＤが補間演算部25に入力されると、補間演算部25
は、メモリ24から補間係数を読み出し、この補間係数に
基づいて、画像データＤの補間演算を行う。すなわち、
画像データＤのＸ方向に対してＢスプライン補間演算を
施して低い鮮鋭度により滑らかな補間を行い、Ｙ方向に
対してCubic スプライン補間演算により高い鮮鋭度によ
りシャープな補間を行う。
ａk-1 、ａk 、ａk+1 、ａk+2 およびｂk-1 、ｂk 、ｂ
k+1 、ｂk+2 は、メモリ24に記憶される。そして、画像
データＤが補間演算部25に入力されると、補間演算部25
は、メモリ24から補間係数を読み出し、この補間係数に
基づいて、画像データＤの補間演算を行う。すなわち、
画像データＤのＸ方向に対してＢスプライン補間演算を
施して低い鮮鋭度により滑らかな補間を行い、Ｙ方向に
対してCubic スプライン補間演算により高い鮮鋭度によ
りシャープな補間を行う。
【００３９】このようにして補間演算処理が施されて、
所定の拡大率に拡大された画像データは、画像処理部26
に入力され、周波数処理などの画像処理が施され、さら
に画像処理が施された画像データＤ′は感熱記録装置28
によりハードコピー29として再生される。
所定の拡大率に拡大された画像データは、画像処理部26
に入力され、周波数処理などの画像処理が施され、さら
に画像処理が施された画像データＤ′は感熱記録装置28
によりハードコピー29として再生される。
【００４０】このようにして再生された放射線画像は、
感熱記録装置28のＸ、Ｙ方向のレスポンスの差が補償さ
れ、再生画像のＸ、Ｙ方向の最終的なレスポンスが互い
に略等しく揃えらる。これにより、Ｘ、Ｙ方向のレスポ
ンスのバランスが異なることを防止することができ、ま
た、Ｘ、Ｙ方向のうちの一方向について過度の周波数強
調処理をかけてアーチファクトを発生させてしまうこと
も防止することができる。
感熱記録装置28のＸ、Ｙ方向のレスポンスの差が補償さ
れ、再生画像のＸ、Ｙ方向の最終的なレスポンスが互い
に略等しく揃えらる。これにより、Ｘ、Ｙ方向のレスポ
ンスのバランスが異なることを防止することができ、ま
た、Ｘ、Ｙ方向のうちの一方向について過度の周波数強
調処理をかけてアーチファクトを発生させてしまうこと
も防止することができる。
【００４１】なお、上述した実施の形態においては、原
画像データのＸ方向にＢスプライン補間演算を、Ｙ方向
にCubic スプライン補間演算を施すようにしているが、
Ｘ、Ｙ方向のそれぞれについて、Ｂスプライン補間およ
びCubic スプライン補間の双方をそれぞれの補間係数に
重み付けをして補間演算を行うようにしてもよい。すな
わち、下記の式(20) Ｆ ＝α・Ａ＋（１−α）・Ｂ (20) 但し、Ｆ：新たな補間係数 Ａ，Ｂ：補間係数 α：０より小さい範囲および／または１より大きい範囲
を含む重み係数 において、ＡをCubic スプライン補間の補間係数、Ｂを
Ｂスプライン補間の補間係数とし、原画像データのＸ、
Ｙ方向のそれぞれの最終的なレスポンスが略等しくなる
ように、重み係数αの値を変更して補間演算を行う。例
えば、原画像データＹk-1 、Ｙk 、Ｙk+1 、Ｙk+2 に対
するCubic スプライン補間係数をａk-1 、ａk 、
ａk+1 、ａk+2 、Ｂスプライン補間係数をｂk-1 、
ｂk 、ｂk+1 、ｂk+2としたとき、補間値Ｓ′は、 Ｓ′＝｛α・ａk-1 ＋（１−α）・ｂk-1 ｝Ｙk-1 ＋｛α・ａk ＋（１−α）・ｂk ｝Ｙk ＋｛α・ａk+1 ＋（１−α）・ｂk+1 ｝Ｙｋ＋１ ＋｛α・ａｋ＋２ ＋（１−α）・ｂk+2 ｝Ｙk+2 (21) となる。
画像データのＸ方向にＢスプライン補間演算を、Ｙ方向
にCubic スプライン補間演算を施すようにしているが、
Ｘ、Ｙ方向のそれぞれについて、Ｂスプライン補間およ
びCubic スプライン補間の双方をそれぞれの補間係数に
重み付けをして補間演算を行うようにしてもよい。すな
わち、下記の式(20) Ｆ ＝α・Ａ＋（１−α）・Ｂ (20) 但し、Ｆ：新たな補間係数 Ａ，Ｂ：補間係数 α：０より小さい範囲および／または１より大きい範囲
を含む重み係数 において、ＡをCubic スプライン補間の補間係数、Ｂを
Ｂスプライン補間の補間係数とし、原画像データのＸ、
Ｙ方向のそれぞれの最終的なレスポンスが略等しくなる
ように、重み係数αの値を変更して補間演算を行う。例
えば、原画像データＹk-1 、Ｙk 、Ｙk+1 、Ｙk+2 に対
するCubic スプライン補間係数をａk-1 、ａk 、
ａk+1 、ａk+2 、Ｂスプライン補間係数をｂk-1 、
ｂk 、ｂk+1 、ｂk+2としたとき、補間値Ｓ′は、 Ｓ′＝｛α・ａk-1 ＋（１−α）・ｂk-1 ｝Ｙk-1 ＋｛α・ａk ＋（１−α）・ｂk ｝Ｙk ＋｛α・ａk+1 ＋（１−α）・ｂk+1 ｝Ｙｋ＋１ ＋｛α・ａｋ＋２ ＋（１−α）・ｂk+2 ｝Ｙk+2 (21) となる。
【００４２】このように、Ｘ、Ｙ方向のそれぞれの補間
係数を重み付けして、補間演算処理を行うことにより、
Ｘ、Ｙ方向のレスポンスを等しくするようにレスポンス
を微調整することができる。
係数を重み付けして、補間演算処理を行うことにより、
Ｘ、Ｙ方向のレスポンスを等しくするようにレスポンス
を微調整することができる。
【００４３】さらに、Ｘ、Ｙ方向の原画像データに対し
て補間演算処理を施す際に、例えば、特開昭55-163472
号公報に記載されたボケマスク処理などの周波数処理を
施し、これによりＸ、Ｙ方向のそれぞれのレスポンスを
調整するようにしてもよい。
て補間演算処理を施す際に、例えば、特開昭55-163472
号公報に記載されたボケマスク処理などの周波数処理を
施し、これによりＸ、Ｙ方向のそれぞれのレスポンスを
調整するようにしてもよい。
【００４４】これは、Ｘ、Ｙ方向にそれぞれ延びる１走
査ラインごとに、そのラインを構成している各画素に関
する原画像データＤorg に対して、下記の演算処理 Ｄ′＝Ｄorg ＋β（Ｄorg −Ｄus） に相当するフィルタリング処理を行うものである。なお
βは強調係数、Ｄusは各画素に関する超低空間周波数に
対応するボケマスクデータで、マスクサイズをＮ、その
マスクサイズ内の画素の画素値の合計値をΣとすると、
Ｄus＝Σ／Ｎである。
査ラインごとに、そのラインを構成している各画素に関
する原画像データＤorg に対して、下記の演算処理 Ｄ′＝Ｄorg ＋β（Ｄorg −Ｄus） に相当するフィルタリング処理を行うものである。なお
βは強調係数、Ｄusは各画素に関する超低空間周波数に
対応するボケマスクデータで、マスクサイズをＮ、その
マスクサイズ内の画素の画素値の合計値をΣとすると、
Ｄus＝Σ／Ｎである。
【００４５】例えば、Ｎ＝３とし、Ｘ方向に並ぶ３つの
画素の原画像データをｘ1 、ｘ2 、ｘ3 とすると、中央
の画素の原画像データｘ2 については、 Ｄ′＝ｘ2 ＋β｛ｘ2 −（ｘ1 ＋ｘ2 ＋ｘ3 ）／３｝ ＝（−β／３）ｘ1 ＋（１＋２β／３）ｘ2 ＋（−β／３）ｘ3 であるから、 ｘ2 ′＝ａｘ1 ＋ｂｘ2 ＋ｃｘ3 （ｂ＝１＋２β／３、ａ＝ｃ＝−β／３） ……（22） と置き代えるフィルタリング処理がなされる。
画素の原画像データをｘ1 、ｘ2 、ｘ3 とすると、中央
の画素の原画像データｘ2 については、 Ｄ′＝ｘ2 ＋β｛ｘ2 −（ｘ1 ＋ｘ2 ＋ｘ3 ）／３｝ ＝（−β／３）ｘ1 ＋（１＋２β／３）ｘ2 ＋（−β／３）ｘ3 であるから、 ｘ2 ′＝ａｘ1 ＋ｂｘ2 ＋ｃｘ3 （ｂ＝１＋２β／３、ａ＝ｃ＝−β／３） ……（22） と置き代えるフィルタリング処理がなされる。
【００４６】このように、原画像データのＸ、Ｙ方向の
それぞれに対して補間演算処理を施す際に、Ｘ、Ｙ方向
のレスポンスを調整するための周波数処理を施すことに
より、上述した実施の形態と同様に、感熱記録装置28の
Ｘ、Ｙ方向のレスポンスの差が補償され、再生画像の
Ｘ、Ｙ方向の最終的なレスポンスが互いに略等しく揃え
られる。これにより、Ｘ、Ｙ方向のレスポンスのバラン
スが異なることを防止することができ、また、Ｘ、Ｙ方
向のうちの一方向について過度の周波数強調処理をかけ
てアーチファクトを発生させてしまうことも防止するこ
とができる。
それぞれに対して補間演算処理を施す際に、Ｘ、Ｙ方向
のレスポンスを調整するための周波数処理を施すことに
より、上述した実施の形態と同様に、感熱記録装置28の
Ｘ、Ｙ方向のレスポンスの差が補償され、再生画像の
Ｘ、Ｙ方向の最終的なレスポンスが互いに略等しく揃え
られる。これにより、Ｘ、Ｙ方向のレスポンスのバラン
スが異なることを防止することができ、また、Ｘ、Ｙ方
向のうちの一方向について過度の周波数強調処理をかけ
てアーチファクトを発生させてしまうことも防止するこ
とができる。
【００４７】また、上述した実施の形態においては、Cu
bic スプライン補間演算を行うための上記式（６）およ
び式（７）の代わりに、パラメータαを付加した下記の
式（６′）および（７′） ｆk ′（Ｘk ）＝α（Ｙk+1 −Ｙk-1 ）／（Ｘk+1 −Ｘk-1 ） （６′） ｆk ′（Ｘk+1 ）＝α（Ｙk+2 −Ｙk ）／（Ｘk+2 −Ｘk ） （７′） を用い、パラメータαを種々変更して補間演算処理を行
うようにしてもよい。
bic スプライン補間演算を行うための上記式（６）およ
び式（７）の代わりに、パラメータαを付加した下記の
式（６′）および（７′） ｆk ′（Ｘk ）＝α（Ｙk+1 −Ｙk-1 ）／（Ｘk+1 −Ｘk-1 ） （６′） ｆk ′（Ｘk+1 ）＝α（Ｙk+2 −Ｙk ）／（Ｘk+2 −Ｘk ） （７′） を用い、パラメータαを種々変更して補間演算処理を行
うようにしてもよい。
【００４８】すなわち、従来のCubic スプライン補間演
算では式（６′）および（７′）におけるパラメータα
が、「１」に固定されていたためこのCubic スプライン
補間演算を単独で、あるいはＢスプライン補間演算を単
独で行った場合には、常に１つの鮮鋭度の２次画像しか
得られず、また前述の特開平2-278478号に示したCubic
スプライン補間演算とＢスプライン補間演算とを組み合
わせたスプライン補間演算ではCubic スプライン補間演
算で得られる鮮鋭度とＢスプライン補間演算で得られる
鮮鋭度との範囲内でしか鮮鋭度の調整はできないが、パ
ラメータαを変更できるようにし、例えば、パラメータ
αを「１」より大きく設定することにより従来のCubic
スプライン補間演算で得られる鮮鋭度より高い鮮鋭度の
２次画像を得ることができ、このパラメータαを大きく
するにしたがって鮮鋭度を高くすることができる。ま
た、パラメータαを「１」より小さく設定することによ
り従来のCubic スプライン補間演算で得られる鮮鋭度よ
り低い鮮鋭度の２次画像を得ることができ、このパラメ
ータαを小さくするにしたがって鮮鋭度を低くすること
ができる。したがって、このパラメータαの値を任意に
変更して、所望の鮮鋭度の補間画像データを得ることが
できる。
算では式（６′）および（７′）におけるパラメータα
が、「１」に固定されていたためこのCubic スプライン
補間演算を単独で、あるいはＢスプライン補間演算を単
独で行った場合には、常に１つの鮮鋭度の２次画像しか
得られず、また前述の特開平2-278478号に示したCubic
スプライン補間演算とＢスプライン補間演算とを組み合
わせたスプライン補間演算ではCubic スプライン補間演
算で得られる鮮鋭度とＢスプライン補間演算で得られる
鮮鋭度との範囲内でしか鮮鋭度の調整はできないが、パ
ラメータαを変更できるようにし、例えば、パラメータ
αを「１」より大きく設定することにより従来のCubic
スプライン補間演算で得られる鮮鋭度より高い鮮鋭度の
２次画像を得ることができ、このパラメータαを大きく
するにしたがって鮮鋭度を高くすることができる。ま
た、パラメータαを「１」より小さく設定することによ
り従来のCubic スプライン補間演算で得られる鮮鋭度よ
り低い鮮鋭度の２次画像を得ることができ、このパラメ
ータαを小さくするにしたがって鮮鋭度を低くすること
ができる。したがって、このパラメータαの値を任意に
変更して、所望の鮮鋭度の補間画像データを得ることが
できる。
【００４９】このように、原画像データのＸ、Ｙ方向の
それぞれに対して補間演算を施す際に、Ｘ、Ｙ方向のレ
スポンスを調整するためのパラメータαを設定すること
により、上述した実施の形態と同様に、感熱記録装置28
のＸ、Ｙ方向のレスポンスの差が補償され、再生画像の
Ｘ、Ｙ方向の最終的なレスポンスが互いに略等しく揃え
られる。これにより、Ｘ、Ｙ方向のレスポンスのバラン
スが異なることを防止することができ、また、Ｘ、Ｙ方
向のうちの一方向について過度の周波数強調処理をかけ
てアーチファクトを発生させてしまうことも防止するこ
とができる
それぞれに対して補間演算を施す際に、Ｘ、Ｙ方向のレ
スポンスを調整するためのパラメータαを設定すること
により、上述した実施の形態と同様に、感熱記録装置28
のＸ、Ｙ方向のレスポンスの差が補償され、再生画像の
Ｘ、Ｙ方向の最終的なレスポンスが互いに略等しく揃え
られる。これにより、Ｘ、Ｙ方向のレスポンスのバラン
スが異なることを防止することができ、また、Ｘ、Ｙ方
向のうちの一方向について過度の周波数強調処理をかけ
てアーチファクトを発生させてしまうことも防止するこ
とができる
【図１】本発明の方法を実施する装置の一例を示す概略
構成図
構成図
【図２】本発明に関わる、画像再生装置のＸ、Ｙ各方向
のレスポンス特性を示すグラフ
のレスポンス特性を示すグラフ
【図３】等間隔の周期でサンプリングされた一方向に配
列されたサンプリング点（画素）の原画像データから補
間画像データを求める方法を説明するグラフ
列されたサンプリング点（画素）の原画像データから補
間画像データを求める方法を説明するグラフ
10 蓄積性蛍光体シート 11 シート搬送手段 12 励起光源 13 レーザビーム 14 回転多面鏡 15 輝尽発光光 16 集光体 17 フォトマルチプライヤー 23 バッファメモリ 24 メモリ 25 補間演算処理部 26 画像処理部 27 画像処理装置 28 感熱記録装置 29 ハードコピー
Claims (2)
- 【請求項１】 所定の間隔でＸ、Ｙ方向に格子状に配さ
れた画素のそれぞれについての画素値を示す原画像デー
タについて、補間演算処理により該原画像データとは間
隔が異なる補間画像データを求め、該補間画像データ
を、Ｘ、Ｙ走査により画像を再生する画像再生装置に送
る画像データの補間演算方法において、 前記画像再生装置のＸ、Ｙ方向のレスポンスの差に応じ
て、該Ｘ、Ｙ方向について互いに異なる補間演算処理を
行うことにより、前記画像再生装置において再生される
画像のＸ、Ｙ方向の最終的なレスポンスを互いに略等し
く揃えることを特徴とする画像データ補間演算方法。 - 【請求項２】 所定の間隔でＸ、Ｙ方向に格子状に配さ
れた画素のそれぞれについての画素値を示す原画像デー
タについて、補間演算処理手段により該原画像データと
は間隔が異なる補間画像データを求め、該補間画像デー
タを、Ｘ、Ｙ走査により画像を再生する画像再生装置に
送る画像データの補間演算装置において、 前記補間演算処理手段が、前記画像再生装置のＸ、Ｙ方
向のレスポンスの差に応じて、該Ｘ、Ｙ方向について互
いに異なる補間演算処理を行い、前記画像再生装置にお
いて再生される画像のＸ、Ｙ方向の最終的なレスポンス
を互いに略等しく揃える手段であることを特徴とする画
像データ補間演算装置。
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP7177014A JPH0927901A (ja) | 1995-07-13 | 1995-07-13 | 画像データ補間演算方法および装置 |
| US08/678,612 US5726766A (en) | 1995-07-13 | 1996-07-15 | Interpolating operation method and apparatus for image signals |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP7177014A JPH0927901A (ja) | 1995-07-13 | 1995-07-13 | 画像データ補間演算方法および装置 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JPH0927901A true JPH0927901A (ja) | 1997-01-28 |
Family
ID=16023662
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP7177014A Withdrawn JPH0927901A (ja) | 1995-07-13 | 1995-07-13 | 画像データ補間演算方法および装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JPH0927901A (ja) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6760489B1 (en) | 1998-04-06 | 2004-07-06 | Seiko Epson Corporation | Apparatus and method for image data interpolation and medium on which image data interpolation program is recorded |
-
1995
- 1995-07-13 JP JP7177014A patent/JPH0927901A/ja not_active Withdrawn
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6760489B1 (en) | 1998-04-06 | 2004-07-06 | Seiko Epson Corporation | Apparatus and method for image data interpolation and medium on which image data interpolation program is recorded |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US4747052A (en) | Radiation image processing | |
| US5454044A (en) | Apparatus for enhancing image data using a monotonously decreasing function | |
| US5319719A (en) | Processing apparatus for radiographic image signals | |
| JP3196864B2 (ja) | 画像のダイナミックレンジ圧縮処理方法 | |
| US5737101A (en) | Interpolating operation method and apparatus for image signals | |
| US5301107A (en) | Method for forming energy subtraction images | |
| JP3494764B2 (ja) | 画像データの補間演算方法および装置 | |
| US5726766A (en) | Interpolating operation method and apparatus for image signals | |
| JP2952418B2 (ja) | 放射線画像信号生成方法 | |
| JPH06274615A (ja) | 画像処理方法 | |
| JPH06274614A (ja) | 画像処理方法 | |
| JPH0927901A (ja) | 画像データ補間演算方法および装置 | |
| JPH0927902A (ja) | 画像データ補間演算方法および装置 | |
| JP3494787B2 (ja) | 画像データの補間演算方法およびその方法を使用した装置 | |
| US6587594B1 (en) | Moire eliminating filter and image processing method and apparatus using the filter | |
| JP3399684B2 (ja) | 画像処理方法および装置 | |
| JPH0950516A (ja) | 画像データ補間演算方法および装置 | |
| JPH01237887A (ja) | 放射線画像の粒状雑音低減方法 | |
| JPH09160149A (ja) | 画像処理方法および装置 | |
| JPH0951429A (ja) | 画像データ補間演算方法および装置 | |
| US5237176A (en) | Method of energy subtraction for radiation image and apparatus for carrying out the method | |
| JP3604201B2 (ja) | 画像データのスプライン補間演算方法およびその方法を使用した装置 | |
| US7366358B2 (en) | Method of and system for generating image signal | |
| JPH0247778A (ja) | 放射線画像処理装置 | |
| JPH0951430A (ja) | 画像データ補間演算方法および装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A300 | Application deemed to be withdrawn because no request for examination was validly filed |
Free format text: JAPANESE INTERMEDIATE CODE: A300 Effective date: 20021001 |