JPH0927902A - 画像データ補間演算方法および装置 - Google Patents
画像データ補間演算方法および装置Info
- Publication number
- JPH0927902A JPH0927902A JP7177015A JP17701595A JPH0927902A JP H0927902 A JPH0927902 A JP H0927902A JP 7177015 A JP7177015 A JP 7177015A JP 17701595 A JP17701595 A JP 17701595A JP H0927902 A JPH0927902 A JP H0927902A
- Authority
- JP
- Japan
- Prior art keywords
- image data
- interpolation
- original image
- enlargement ratio
- 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.)
- Withdrawn
Links
Classifications
-
- 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/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)
- Image Processing (AREA)
- Editing Of Facsimile Originals (AREA)
- Radiography Using Non-Light Waves (AREA)
Abstract
(57)【要約】
【課題】 原画像データを所定の拡大率により拡大する
際に行う補間演算処理において、拡大率に拘わらず画像
再生装置において再生される補間画像のレスポンスを略
等しくする。 【解決手段】 原画像データＤに対して補間演算処理を
施し、原画像データＤとは間隔が異なる補間画像データ
を求め、補間画像データを画像を再生する感熱記録装置
28に送り、ハードコピー29として再生する。この際、Cu
bic スプライン、Ｂスプラインなどの補間方法により、
原画像に対する拡大率に応じて原画像データＤに異なる
補間演算処理を行う。これにより拡大率に拘わらず補間
画像のレスポンスが略等しく揃えられる。
際に行う補間演算処理において、拡大率に拘わらず画像
再生装置において再生される補間画像のレスポンスを略
等しくする。 【解決手段】 原画像データＤに対して補間演算処理を
施し、原画像データＤとは間隔が異なる補間画像データ
を求め、補間画像データを画像を再生する感熱記録装置
28に送り、ハードコピー29として再生する。この際、Cu
bic スプライン、Ｂスプラインなどの補間方法により、
原画像に対する拡大率に応じて原画像データＤに異なる
補間演算処理を行う。これにより拡大率に拘わらず補間
画像のレスポンスが略等しく揃えられる。
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からなる。
【００１７】ここで、拡大率による補間画像データのレ
スポンス特性は、図２に示すようなものとなる。図２に
おいて、拡大率ａ，ｂ，ｃに拘わらずレスポンスを揃え
たい空間周波数をαとし、この空間周波数αに対応する
補間画像データのＭＴＦ（Modulation Transfer Functi
on，周波数依存特性）をそれぞれｍ，ｍ′，ｍ″とする
と、空間周波数αにおけるレスポンスを揃えるために
は、拡大率ｂのレスポンスが拡大率ａのレスポンスに対
してｍ／ｍ′となるような関数の１次元補間を、また、
拡大率ｃのレスポンスが拡大率ａのレスポンスに対して
ｍ／ｍ″となるような関数の１次元補間を画像データＤ
に施す必要がある。そしてこのために、補間演算部25に
おいては、補間画像のレスポンスが各拡大率ａ，ｂ，ｃ
において略等しくなるような補間演算処理が画像データ
Ｄに対して施される。
スポンス特性は、図２に示すようなものとなる。図２に
おいて、拡大率ａ，ｂ，ｃに拘わらずレスポンスを揃え
たい空間周波数をαとし、この空間周波数αに対応する
補間画像データのＭＴＦ（Modulation Transfer Functi
on，周波数依存特性）をそれぞれｍ，ｍ′，ｍ″とする
と、空間周波数αにおけるレスポンスを揃えるために
は、拡大率ｂのレスポンスが拡大率ａのレスポンスに対
してｍ／ｍ′となるような関数の１次元補間を、また、
拡大率ｃのレスポンスが拡大率ａのレスポンスに対して
ｍ／ｍ″となるような関数の１次元補間を画像データＤ
に施す必要がある。そしてこのために、補間演算部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 ）を表すものである。
Ｂk ｘ＋Ｃk ）を表すものである。
【００２４】また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として再生される。
【００３９】このようにして再生された補間画像は、拡
大率によるレスポンスの差が補償され、再生画像のレス
ポンスが拡大率に拘わらず略等しく揃えられる。これに
より、拡大率によりレスポンスが異なることを防止する
ことができる。
大率によるレスポンスの差が補償され、再生画像のレス
ポンスが拡大率に拘わらず略等しく揃えられる。これに
より、拡大率によりレスポンスが異なることを防止する
ことができる。
【００４０】なお、上述した実施の形態においては、拡
大率ａの場合にＢスプライン補間演算を、拡大率ｂの場
合にCubic スプライン補間演算を施すようにしている
が、画像データＤに対してＢスプライン補間およびCubi
c スプライン補間の双方を用いて、拡大率に応じてそれ
ぞれの補間係数に重み付けをして補間演算を行うように
してもよい。すなわち、下記の式(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+1 ＋｛α・ａk+2 ＋（１−α）・ｂk+2 ｝Ｙk+2 (21) となる。
大率ａの場合にＢスプライン補間演算を、拡大率ｂの場
合にCubic スプライン補間演算を施すようにしている
が、画像データＤに対してＢスプライン補間およびCubi
c スプライン補間の双方を用いて、拡大率に応じてそれ
ぞれの補間係数に重み付けをして補間演算を行うように
してもよい。すなわち、下記の式(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+1 ＋｛α・ａk+2 ＋（１−α）・ｂk+2 ｝Ｙk+2 (21) となる。
【００４１】そして、拡大率ａの場合には、Ｂスプライ
ン補間の補間係数の重み付けが大きくなるように、重み
係数ｔの値を小さくして補間演算を行い、拡大率ｂの場
合には、Cubic スプライン補間の補間係数の重み付けが
大きくなるように、重み係数ｔの値を大きくして補間演
算を行う。このように、Cubic スプライン補間およびＢ
スプライン補間のそれぞれの補間係数を重み付けして、
補間演算処理を行うことにより、拡大率に拘わらず再生
された補間画像のレスポンスを等しくすることができ
る。
ン補間の補間係数の重み付けが大きくなるように、重み
係数ｔの値を小さくして補間演算を行い、拡大率ｂの場
合には、Cubic スプライン補間の補間係数の重み付けが
大きくなるように、重み係数ｔの値を大きくして補間演
算を行う。このように、Cubic スプライン補間およびＢ
スプライン補間のそれぞれの補間係数を重み付けして、
補間演算処理を行うことにより、拡大率に拘わらず再生
された補間画像のレスポンスを等しくすることができ
る。
【００４２】さらに、原画像データに対して拡大率に応
じた補間演算処理を施す際に、例えば、特開昭55-16347
2 号公報に記載されたボケマスク処理などの周波数処理
を施し、これにより拡大率に応じた補間画像のレスポン
スを調整するようにしてもよい。
じた補間演算処理を施す際に、例えば、特開昭55-16347
2 号公報に記載されたボケマスク処理などの周波数処理
を施し、これにより拡大率に応じた補間画像のレスポン
スを調整するようにしてもよい。
【００４３】これは、拡大率に応じて原画像データＤor
g に対して、下記の演算処理 Ｄ′＝Ｄorg ＋β（Ｄorg −Ｄus） に相当するフィルタリング処理を行うものである。なお
βは強調係数、Ｄusは各画素に関する超低空間周波数に
対応するボケマスクデータで、マスクサイズをＮ、その
マスクサイズ内の画素の画素値の合計値をΣとすると、
Ｄus＝Σ／Ｎである。
g に対して、下記の演算処理 Ｄ′＝Ｄ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） と置き代えるフィルタリング処理がなされる。
【００４５】このように、拡大率に応じて原画像データ
に補間演算処理を施す際に、拡大率に拘わらず補間画像
のレスポンスが略等しくなるような周波数処理を施すこ
とにより、上述した実施の形態と同様に、拡大率による
補間画像のレスポンスの差が補償され、再生画像のレス
ポンスが拡大率に拘わらず略等しく揃えられる。これに
より、拡大率によりレスポンスが異なることを防止する
ことができる。
に補間演算処理を施す際に、拡大率に拘わらず補間画像
のレスポンスが略等しくなるような周波数処理を施すこ
とにより、上述した実施の形態と同様に、拡大率による
補間画像のレスポンスの差が補償され、再生画像のレス
ポンスが拡大率に拘わらず略等しく揃えられる。これに
より、拡大率によりレスポンスが異なることを防止する
ことができる。
【００４６】また、本実施の形態においては、Cubic ス
プライン補間演算を行うための上記式（６）および式
（７）の代わりに、パラメータαを付加した下記の式
（６′）および（７′） ｆk ′（Ｘk ）＝α（Ｙk+1 −Ｙk-1 ）／（Ｘk+1 −Ｘk-1 ） （６′） ｆk ′（Ｘk+1 ）＝α（Ｙk+2 −Ｙk ）／（Ｘk+2 −Ｘk ） （７′） を用い、パラメータαを種々変更して補間演算処理を行
うようにしてもよい。
プライン補間演算を行うための上記式（６）および式
（７）の代わりに、パラメータαを付加した下記の式
（６′）および（７′） ｆk ′（Ｘk ）＝α（Ｙk+1 −Ｙk-1 ）／（Ｘk+1 −Ｘk-1 ） （６′） ｆk ′（Ｘk+1 ）＝α（Ｙk+2 −Ｙk ）／（Ｘk+2 −Ｘk ） （７′） を用い、パラメータαを種々変更して補間演算処理を行
うようにしてもよい。
【００４７】すなわち、従来のCubic スプライン補間演
算では式（６′）および（７′）におけるパラメータα
が、「１」に固定されていたためこのCubic スプライン
補間演算を単独で、あるいはＢスプライン補間演算を単
独で行った場合には、常に１つのレスポンスの２次画像
しか得られず、また前述の特開平2-278478号に示したCu
bic スプライン補間演算とＢスプライン補間演算とを組
み合わせたスプライン補間演算ではCubic スプライン補
間演算で得られるレスポンスとＢスプライン補間演算で
得られるレスポンスとの範囲内でしかレスポンスの調整
はできないが、パラメータαを変更できるようにし、例
えば、パラメータαを「１」より大きく設定することに
より従来のCubic スプライン補間演算で得られるレスポ
ンスより高いレスポンスの２次画像を得ることができ、
このパラメータαを大きくするにしたがってレスポンス
を高くすることができる。また、パラメータαを「１」
より小さく設定することにより従来のCubic スプライン
補間演算で得られるレスポンスより低いレスポンスの２
次画像を得ることができ、このパラメータαを小さくす
るにしたがってレスポンスを低くすることができる。し
たがって、このパラメータαの値を任意に変更して、所
望のレスポンスの補間画像データを得ることができる。
算では式（６′）および（７′）におけるパラメータα
が、「１」に固定されていたためこのCubic スプライン
補間演算を単独で、あるいはＢスプライン補間演算を単
独で行った場合には、常に１つのレスポンスの２次画像
しか得られず、また前述の特開平2-278478号に示したCu
bic スプライン補間演算とＢスプライン補間演算とを組
み合わせたスプライン補間演算ではCubic スプライン補
間演算で得られるレスポンスとＢスプライン補間演算で
得られるレスポンスとの範囲内でしかレスポンスの調整
はできないが、パラメータαを変更できるようにし、例
えば、パラメータαを「１」より大きく設定することに
より従来のCubic スプライン補間演算で得られるレスポ
ンスより高いレスポンスの２次画像を得ることができ、
このパラメータαを大きくするにしたがってレスポンス
を高くすることができる。また、パラメータαを「１」
より小さく設定することにより従来のCubic スプライン
補間演算で得られるレスポンスより低いレスポンスの２
次画像を得ることができ、このパラメータαを小さくす
るにしたがってレスポンスを低くすることができる。し
たがって、このパラメータαの値を任意に変更して、所
望のレスポンスの補間画像データを得ることができる。
【００４８】このように、拡大率に応じて原画像データ
に対して補間演算を施す際に、レスポンスを調整するた
めのパラメータαを設定することにより、上述した実施
の形態と同様に、拡大率によるレスポンスの差が補償さ
れ、再生画像のレスポンスが拡大率に拘わらず略等しく
揃えられる。これにより、拡大率によりレスポンスが異
なることを防止することができる
に対して補間演算を施す際に、レスポンスを調整するた
めのパラメータαを設定することにより、上述した実施
の形態と同様に、拡大率によるレスポンスの差が補償さ
れ、再生画像のレスポンスが拡大率に拘わらず略等しく
揃えられる。これにより、拡大率によりレスポンスが異
なることを防止することができる
【図１】本発明の方法を実施する装置の一例を示す概略
構成図
構成図
【図２】拡大率による補間画像のレスポンス特性を示す
グラフ
グラフ
【図３】等間隔の周期でサンプリングされた一方向に配
列されたサンプリング点（画素）の原画像データから補
間画像データを求める方法を説明するグラフ
列されたサンプリング点（画素）の原画像データから補
間画像データを求める方法を説明するグラフ
10 蓄積性蛍光体シート 11 シート搬送手段 12 励起光源 13 レーザビーム 14 回転多面鏡 15 輝尽発光光 16 集光体 17 フォトマルチプライヤー 23 バッファメモリ 24 メモリ 25 補間演算処理部 26 画像処理部 27 画像処理装置 28 感熱記録装置 29 ハードコピー
─────────────────────────────────────────────────────
【手続補正書】
【提出日】平成８年１０月２日
【手続補正１】
【補正対象書類名】明細書
【補正対象項目名】００２３
【補正方法】変更
【補正内容】
【００２３】 ｆ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 スプライ
ン補間演算は厳密には第２階微分係数の連続条件も含む
ものであるが、式が複雑化するため上述したように簡略
化して用いるのが一般的である。
【手続補正２】
【補正対象書類名】明細書
【補正対象項目名】００２６
【補正方法】変更
【補正内容】
【００２６】 ｆ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 について整理すると、下記式（９）で表すことがで
きる。
【手続補正３】
【補正対象書類名】明細書
【補正対象項目名】００３４
【補正方法】変更
【補正内容】
【００３４】 ｆ（Ｘ）＝ｆ(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）および（14）〜（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）および（14）〜（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）で表すことがで
きる。
【手続補正４】
【補正対象書類名】明細書
【補正対象項目名】００４０
【補正方法】変更
【補正内容】
【００４０】なお、上述した実施の形態においては、拡
大率ａの場合にＢスプライン補間演算を、拡大率ｂの場
合にCubic スプライン補間演算を施すようにしている
が、画像データＤに対してＢスプライン補間およびCubi
c スプライン補間の双方を用いて、拡大率に応じてそれ
ぞれの補間係数に重み付けをして補間演算を行うように
してもよい。すなわち、下記の式(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+1 ＋｛α・ａk+2 ＋（１−α）・ｂk+2 ｝Ｙk+2 (21) となる。
大率ａの場合にＢスプライン補間演算を、拡大率ｂの場
合にCubic スプライン補間演算を施すようにしている
が、画像データＤに対してＢスプライン補間およびCubi
c スプライン補間の双方を用いて、拡大率に応じてそれ
ぞれの補間係数に重み付けをして補間演算を行うように
してもよい。すなわち、下記の式(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+1 ＋｛α・ａk+2 ＋（１−α）・ｂk+2 ｝Ｙk+2 (21) となる。
【手続補正５】
【補正対象書類名】明細書
【補正対象項目名】００４１
【補正方法】変更
【補正内容】
【００４１】そして、拡大率ａの場合には、Ｂスプライ
ン補間の補間係数の重み付けが大きくなるように、重み
係数αの値を小さくして補間演算を行い、拡大率ｂの場
合には、Cubic スプライン補間の補間係数の重み付けが
大きくなるように、重み係数αの値を大きくして補間演
算を行う。このように、Cubic スプライン補間およびＢ
スプライン補間のそれぞれの補間係数を重み付けして、
補間演算処理を行うことにより、拡大率に拘わらず再生
された補間画像のレスポンスを等しくすることができ
る。
ン補間の補間係数の重み付けが大きくなるように、重み
係数αの値を小さくして補間演算を行い、拡大率ｂの場
合には、Cubic スプライン補間の補間係数の重み付けが
大きくなるように、重み係数αの値を大きくして補間演
算を行う。このように、Cubic スプライン補間およびＢ
スプライン補間のそれぞれの補間係数を重み付けして、
補間演算処理を行うことにより、拡大率に拘わらず再生
された補間画像のレスポンスを等しくすることができ
る。
Claims (2)
- 【請求項１】 所定の間隔で格子状に配された画素のそ
れぞれについての画素値を示す原画像データについて、
該原画像データにより表される原画像に対する拡大率に
応じた補間演算処理により該原画像データとは間隔が異
なる補間画像データを求め、該補間画像データを画像再
生装置に送る画像データの補間演算方法において、 前記拡大率に応じて前記原画像データに異なる補間演算
処理を行うことにより、前記画像再生装置に表示される
前記補間画像データのレスポンスを前記拡大率に拘わら
ず略等しくすることを特徴とする画像データ補間演算方
法。 - 【請求項２】 所定の間隔で格子状に配された画素のそ
れぞれについての画素値を示す原画像データについて、
該原画像データにより表される原画像に対する拡大率に
応じた補間演算処理手段により該原画像データとは間隔
が異なる補間画像データを求め、該補間画像データを画
像再生装置に送る画像データの補間演算装置において、 前記補間演算処理手段が、前記拡大率に応じて前記原画
像データに異なる補間演算処理を行うことにより、前記
画像再生装置に表示される前記補間画像データのレスポ
ンスを前記拡大率に拘わらず略等しくする手段であるこ
とを特徴とする画像データ補間演算装置。
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP7177015A JPH0927902A (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 |
|---|---|---|---|
| JP7177015A JPH0927902A (ja) | 1995-07-13 | 1995-07-13 | 画像データ補間演算方法および装置 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JPH0927902A true JPH0927902A (ja) | 1997-01-28 |
Family
ID=16023679
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP7177015A Withdrawn JPH0927902A (ja) | 1995-07-13 | 1995-07-13 | 画像データ補間演算方法および装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JPH0927902A (ja) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO1999041703A1 (en) * | 1998-02-17 | 1999-08-19 | Matsushita Electric Industrial Co., Ltd. | Method and apparatus for interpolated pixel generation |
-
1995
- 1995-07-13 JP JP7177015A patent/JPH0927902A/ja not_active Withdrawn
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO1999041703A1 (en) * | 1998-02-17 | 1999-08-19 | Matsushita Electric Industrial Co., Ltd. | Method and apparatus for interpolated pixel generation |
| AU750460B2 (en) * | 1998-02-17 | 2002-07-18 | Matsushita Electric Industrial Co., Ltd. | Method and apparatus for interpolated pixel generation |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US5454044A (en) | Apparatus for enhancing image data using a monotonously decreasing function | |
| US4747052A (en) | Radiation image processing | |
| US5905822A (en) | Interpolating operation method and apparatus for image signals | |
| JP3196864B2 (ja) | 画像のダイナミックレンジ圧縮処理方法 | |
| JP3599435B2 (ja) | 画像データ補間演算方法および装置 | |
| JP2849964B2 (ja) | 画像処理方法および装置 | |
| JP3494764B2 (ja) | 画像データの補間演算方法および装置 | |
| US5726766A (en) | Interpolating operation method and apparatus for image signals | |
| JP2952418B2 (ja) | 放射線画像信号生成方法 | |
| JPH06274614A (ja) | 画像処理方法 | |
| JPH0927902A (ja) | 画像データ補間演算方法および装置 | |
| JPH0927901A (ja) | 画像データ補間演算方法および装置 | |
| JPH0951429A (ja) | 画像データ補間演算方法および装置 | |
| JPH0950516A (ja) | 画像データ補間演算方法および装置 | |
| JP3399684B2 (ja) | 画像処理方法および装置 | |
| JPH01237887A (ja) | 放射線画像の粒状雑音低減方法 | |
| US5237176A (en) | Method of energy subtraction for radiation image and apparatus for carrying out the method | |
| JPH0993426A (ja) | 画像データの補間演算方法およびその方法を使用した装置 | |
| JPH09160149A (ja) | 画像処理方法および装置 | |
| JP3604201B2 (ja) | 画像データのスプライン補間演算方法およびその方法を使用した装置 | |
| US7366358B2 (en) | Method of and system for generating image signal | |
| JPH0951430A (ja) | 画像データ補間演算方法および装置 | |
| JP3205416B2 (ja) | 画像データ補正方法 | |
| JPH05232608A (ja) | 放射線画像の重ね合せ処理方法および装置 | |
| JPH01106275A (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 |