JP4350468B2

JP4350468B2 - 画像処理方法および装置ならびにプログラム

Info

Publication number: JP4350468B2
Application number: JP2003315677A
Authority: JP
Inventors: 孝夫桑原
Original assignee: Fujifilm Corp
Current assignee: Fujifilm Corp
Priority date: 2003-09-08
Filing date: 2003-09-08
Publication date: 2009-10-21
Anticipated expiration: 2023-09-08
Also published as: US7657112B2; US20050053306A1; EP1513337A2; EP1513337B1; EP1513337A3; JP2005084902A

Description

本発明は画像処理方法および装置ならびにプログラムに関し、詳しくは、画像中の線状のノイズ成分を除去する画像処理方法および装置ならびにプログラムに関するものである。

従来より、多数の受光部を主走査方向に並べてなる光電変換素子列であるラインセンサに対して、画像担持体である原稿を上記主走査方向と交わる副走査方向に相対的に移動させながら、この原稿から発せられた光を上記ラインセンサで検出して、この原稿が担持する画像を示す画像データを取得する装置が知られている。このような装置においては、ラインセンサを構成する各受光部毎の特性の違い、例えば、各受光部毎のゲインや温度特性の違い等により、このラインセンサで検出して得られた画像データが示す画像中に、上記副走査方向に延びる線状のノイズ（以後、スジムラともいう）が現れることがある。また、ラインセンサへ加えられた衝撃により、このラインセンサで検出して得られた画像データが示す画像中に、上記主走査方向に延びるスジムラが現れることもある。そのため、例えば、このスジムラを含む画像を示す画像データをフィルタに通してこのスジムラ成分を抽出した後、このスジムラ成分を示す画像データを上記画像を示す画像データから差し引いてこのスジムラ成分を上記画像から除去する方法が提案されている（例えば、特許文献１参照）。
特開２００３−１５０９５４号公報

しかしながら、上記手法は、スジムラ成分を含む画像データ中に急激に値が変化する領域が存在すると、この画像データをフィルタに通してスジムラ成分を抽出するときに上記急激に画像データの値が変化する領域の近傍にリンギングが生じ、上記抽出したスジムラ成分を示す画像データ中にリンギングを示す成分が残ってしまうことがあり、このような場合には、上記スジムラ成分を除去して得られた画像中に上記リンギングを示す成分がスジムラのようになって現れる。また、上記スジムラ成分と線状画像情報成分とを分離して、線状画像情報成分を除去することなく画像中からスジムラ成分のみを除去したいという要請もある。

本発明は、上記事情に鑑みてなされたものであり、本発明の第１の目的は、画像中のスジムラ成分を抽出する処理におけるリンギングの発生を抑制し、これにより、スジムラ成分を除去して得られる画像の品質を高めることができる画像処理方法および装置ならびにプログラムを提供することを目的とするものであり、本発明の第２の目的は、スジムラ成分の除去に伴う画像中からの線状画像情報成分の消失を抑制することができる画像処理方法および装置ならびにプログラムを提供することを目的とするものである。

本発明の第１の画像処理方法は、画像データにフィルタ処理を施して、該画像データが示す画像中の１方向に延びる線状のノイズ成分を示す画像データを抽出し、該線状のノイズ成分を示す画像データを前記画像を示す画像データから差し引いてこの画像から線状のノイズ成分を除去する画像処理方法において、前記画像を示す画像データに前記フィルタ処理を施す前に、該画像データに対して、該画像中の前記１方向と直交する直交方向における画像データの値の急激な変化を緩やかにする前処理を施すことを特徴とするものである。

本発明の第２の画像処理方法は、多数の受光部を主走査方向に並べてなる光電変換素子列に対して、画像担持体を主走査方向と交わる副走査方向に相対的に移動させながら、この画像担持体から発せられた光を光電変換素子列で検出して得られたこの画像担持体が担持する画像を示す画像データに、フィルタ処理を施して、この画像データが示す画像中の副走査方向および／または前記主走査方向に延びる線状のノイズ成分を示す画像データを抽出し、前記線状のノイズ成分を示す画像データを前記画像を示す画像データから差し引いてこの画像から前記線状のノイズ成分を除去する画像処理方法において、前記画像を示す画像データに前記フィルタ処理を施す前に、この画像データに対して、前記画像中の副走査方向と直交する副走査直交方向、および／または前記画像中の主走査方向と直交する主走査直交方向における画像データの値の急激な変化を緩やかにする前処理を施すことを特徴とするものである。

本発明の画像処理装置は、多数の受光部を主走査方向に並べてなる光電変換素子列に対して、画像担持体を主走査方向と交わる副走査方向に相対的に移動させながら、この画像担持体から発せられた光を光電変換素子列で検出して得られたこの画像担持体が担持する画像を示す画像データにフィルタ処理を施して、この画像データが示す画像中の前記副走査方向および／または前記主走査方向に延びる線状のノイズ成分を示す画像データを抽出する線状ノイズ抽出手段と、この線状のノイズ成分を示す画像データを前記画像を示す画像データから差し引く差引手段とを備え、前記画像から前記線状のノイズ成分を除去する画像処理装置において、線状ノイズ抽出手段が前記画像を示す画像データに前記フィルタ処理を施す前に、この画像データに対して、前記画像中の副走査方向と直交する副走査直交方向、および／または前記画像中の主走査方向と直交する主走査直交方向における画像データの値の急激な変化を緩やかにする前処理手段を備えたことを特徴とするものである。

前記前処理手段は、この前処理手段が副走査直交方向における画像データの値の急激な変化を緩やかにするものである場合には、前記画像中の副走査直交方向において互いに隣り合うこの画像を示す画像データの差分値を示す画像データを得、この差分値を示す画像データに対して、この画像データのうちの所定の閾値を越える値を持つ画像データの値の絶対値を小さくする処理を施して得られた画像データを副走査直交方向に累積加算するものとしたり、前記画像中の副走査直交方向において互いに隣り合うこの画像を示す画像データの差分値を示す画像データを得、この差分値を示す画像データを副走査直交方向に周波数処理し、前記差分値を示す画像データに対して、前記周波数処理した画像データのうちの所定の閾値を越える値を持つ画像データに対応する前記差分値を示す画像データの値の絶対値を小さくする処理を施して得られた画像データを副走査直交方向に累積加算するものとすることができる。

また、前記前処理手段は、この前処理手段が主走査直交方向における画像データの値の急激な変化を緩やかにするものである場合には、前記画像中の主走査直交方向において互いに隣り合うこの画像を示す画像データの差分値を示す画像データを得、この差分値を示す画像データに対して、この画像データのうちの所定の閾値を越える値を持つ画像データの値の絶対値を小さくする処理を施して得られた画像データを主走査直交方向に累積加算するものとしたり、前記画像中の主走査直交方向において互いに隣り合うこの画像を示す画像データの差分値を示す画像データを得、この差分値を示す画像データを主走査直交方向に周波数処理し、前記差分値を示す画像データに対して、前記周波数処理した画像データのうちの所定の閾値を越える値を持つ画像データに対応する前記差分値を示す画像データの値の絶対値を小さくする処理を施して得られた画像データを主走査直交方向に累積加算するものとすることもできる。

本発明の画像処理をコンピュータに実行させるためのプログラムは、多数の受光部を主走査方向に並べてなる光電変換素子列に対して、画像担持体を主走査方向と交わる副走査方向に相対的に移動させながら、この画像担持体から発せられた光を光電変換素子列で検出して得られたこの画像担持体が担持する画像を示す画像データにフィルタ処理を施して、この画像データが示す画像中の副走査方向および／または主走査方向に延びる線状のノイズ成分を示す画像データを抽出する手順と、この線状のノイズ成分を示す画像データを前記画像を示す画像データから差し引いてこの画像から前記線状のノイズ成分を除去する手順とをコンピュータに実行させるためのプログラムにおいて、前記画像を示す画像データに前記フィルタ処理を施す前に、この画像データに対して、前記画像中の副走査方向と直交する副走査直交方向、および／または前記画像中の主走査方向と直交する主走査直交方向における画像データの値の急激な変化を緩やかにする前処理を行なわせる手順を実行させることを特徴とするものである。

なお、前記前処理は、前記画像データが示す画像中の副走査方向に延びる線状のノイズ成分を除去する場合に、前記画像データに対して、前記副走査直交方向における画像データの値の急激な変化を緩やかにする前処理を施し、前記画像データが示す画像中の主走査方向に延びる線状のノイズ成分を除去する場合に、前記主走査直交方向における画像データの値の急激な変化を緩やかにするものである。

前記画像データの値の急激な変化を緩やかにするとは、前記画像データに含まれる線状のノイズ成分が失われないように、かつ、この画像データに前記線状のノイズ成分を抽出する処理を施したときのリンギングの発生が抑制されるように、上記画像データの値の急激な変化を緩やかにすることを意味するものである。

前記画像データの値の絶対値を小さくする処理は、この画像データに前記線状のノイズ成分を抽出する処理を施したときのリンギングの発生が抑制されるように、画像データの値の絶対値を小さくするものである。

前記画像データを累積加算するとは、前記差分値を示す画像データから前記画像が再生されるように累積加算することを意味する。

前記所定の閾値は、前記線状のノイズ成分を抽出する処理を施したときにリンギングが大きく生じる領域を示す画像データと他の領域を示す画像データとを分離するように定められるものである。

前記所定の閾値を越える値とは、所定の閾値が負の値であるときにはこの閾値より小さな値を、所定の閾値が正の値であるときにはこの閾値より大きな値を意味するものである。なお、上記所定の閾値としては、正の値を持つ閾値と負の値を持つ閾値とを定めることが望ましい。

前記周波数処理は移動平均とすることができる。

本発明の第３の画像処理方法は、多数の受光部を主走査方向に並べてなる光電変換素子列に対して、画像担持体を前記主走査方向と交わる副走査方向に相対的に移動させながら、該画像担持体から発せられた光を前記光電変換素子列で検出して得られたこの画像担持体が担持する画像を示す画像データに、フィルタ処理を施して、該画像データが示す画像中の前記副走査方向および／または前記主走査方向に延びる線状のノイズ成分を示す画像データを抽出し、前記線状のノイズ成分を示す各画像データの値の絶対値を全体的に小さくした画像データを得、該画像データを前記画像を示す画像データから差し引いて前記画像から前記線状のノイズ成分を除去することを特徴とするものである。

本発明の他の画像処理装置は、多数の受光部を主走査方向に並べてなる光電変換素子列に対して、画像担持体を前記主走査方向と交わる副走査方向に相対的に移動させながら、該画像担持体から発せられた光を前記光電変換素子列で検出して得られたこの画像担持体が担持する画像を示す画像データにフィルタ処理を施して、該画像データが示す画像中の前記副走査方向および／または前記主走査方向に延びる線状のノイズ成分を示す画像データを抽出する線状ノイズ抽出手段と、前記線状のノイズ成分を示す各画像データの値の絶対値を全体的に小さくした画像データを得るデータ変換手段と、該画像データを前記画像を示す画像データから差し引く差引手段とを備え、前記画像から前記線状のノイズ成分を除去するものであることを特徴とするものである。

本発明の画像処理をコンピュータに実行させるための他のプログラムは、多数の受光部を主走査方向に並べてなる光電変換素子列に対して、画像担持体を前記主走査方向と交わる副走査方向に相対的に移動させながら、該画像担持体から発せられた光を前記光電変換素子列で検出して得られたこの画像担持体が担持する画像を示す画像データにフィルタ処理を施して、該画像データが示す画像中の前記副走査方向および／または前記主走査方向に延びる線状のノイズ成分を示す画像データを抽出する手順と、該線状のノイズ成分を示す各画像データの値の絶対値を全体的に小さくした画像データを得る手順と、該画像データを前記画像を示す画像データから差し引いて前記画像から前記線状のノイズ成分を除去する手順とをコンピュータに実行させることを特徴とするものである。

前記「各画像データの値の絶対値を全体的に小さくする」とは、少なくとも各画像データの値の互いの大小関係を維持したまま上記各画像データの値の絶対値を小さくすることを意味する。

前記光電変換素子列は、多数の受光部が直線状に１列に並ぶ１次元のラインセンサとしたり、あるいは多数の受光部が直線状に複数列並ぶ２次元のエリアセンサとすることができる。

本発明の第１の画像処理方法は、画像データにフィルタ処理を施して、該画像データが示す画像中の１方向に延びる線状のノイズ成分を示す画像データを抽出し、線状のノイズ成分を示す画像データを前記画像を示す画像データから差し引いてこの画像から線状のノイズ成分を除去する画像処理方法において、前記画像を示す画像データに前記フィルタ処理を施す前に、該画像データに対して、該画像中の前記１方向と直交する直交方向における画像データの値の急激な変化を緩やかにする前処理を施すようにしたので、上記線状のノイズ成分であるスジムラ成分を抽出する処理を施すときのリンギングの発生を抑制して上記スジムラ成分を示す画像データを抽出することができ、スジムラ成分を除去して得られる画像への上記リンギングを示す成分の混入を抑制できるので、上記画像の品質を高めることができる。これとともに、スジムラ成分の除去に伴う画像中からの線状画像情報成分の消失を抑制することもできる。

本発明の第２の画像処理方法、および画像処理装置ならびに画像処理プログラムによれば、上記光電変換素子列で検出して得られた画像担持体が担持する画像を示す画像データにフィルタ処理を施してこの画像中の副走査方向および／または主走査方向に延びる線状のノイズ成分であるスジムラ成分を抽出する前に、上記画像データに対して、この画像中の副走査直交方向および／または主走査直交方向における画像データの値の急激な変化を緩やかにする前処理を施すようにしたので、上記スジムラ成分を抽出する処理を施すときのリンギングの発生を抑制して上記スジムラ成分を示す画像データを抽出することができ、スジムラ成分を除去して得られる画像への上記リンギングを示す成分の混入を抑制できるので、上記画像の品質を高めることができる。これとともに、スジムラ成分の除去に伴う画像中からの線状画像情報成分の消失を抑制することもできる。

また、前処理手段を、この前処理手段が副走査直交方向における画像データの値の急激な変化を緩やかにするものである場合に、画像中の副走査直交方向において互いに隣り合うこの画像を示す画像データの差分値を示す画像データを得、この差分値を示す画像データに対して、この画像データのうちの所定の閾値を越える値を持つ画像データの値の絶対値を小さくする処理を施して得られた画像データを副走査直交方向に累積加算するものとしたり、あるいは、前記前処理を、画像中の副走査直交方向において互いに隣り合うこの画像を示す画像データの差分値を示す画像データを得、この差分値を示す画像データを副走査直交方向に周波数処理し、上記差分値を示す画像データに対して、周波数処理した画像データのうちの所定の閾値を越える値を持つ画像データに対応する上記差分値を示す画像データの値の絶対値を小さくする処理を施して得られた画像データを副走査直交方向に累積加算するものとすれば、上記閾値を設定することにより、前記スジムラ成分を抽出する処理によってリンギングが生じ易い領域を示す画像データを抽出することができ、また、上記差分値を示す画像データの値の絶対値を小さくする処理により、上記リンギングが大きく生じる領域を示す画像データの値の変化を緩やかにすることができるので、上記スジムラ成分を抽出する処理を施す際のリンギングの発生をより確実に抑制して上記スジムラ成分を示す画像データを抽出することができる。これにより、上記スジムラ成分を除去して得られる画像への上記リンギングを示す成分の混入をより確実に抑制でき、上記画像の品質を高めることができる。これとともに、スジムラ成分の除去に伴う画像中からの線状画像情報成分の消失をより確実に抑制することもできる。

また、前処理手段を、この前処理手段が主走査直交方向における画像データの値の急激な変化を緩やかにするものである場合に、画像中の主走査直交方向において互いに隣り合うこの画像を示す画像データの差分値を示す画像データを得、この差分値を示す画像データに対して、この画像データのうちの所定の閾値を越える値を持つ画像データの値の絶対値を小さくする処理を施して得られた画像データを主走査直交方向に累積加算するものとしたり、あるいは、前記前処理を、画像中の主走査直交方向において互いに隣り合うこの画像を示す画像データの差分値を示す画像データを得、この差分値を示す画像データを主走査直交方向に周波数処理し、上記差分値を示す画像データに対して、周波数処理した画像データのうちの所定の閾値を越える値を持つ画像データに対応する上記差分値を示す画像データの値の絶対値を小さくする処理を施して得られた画像データを主走査直交方向に累積加算するものとすれば、上記と同様に、上記スジムラ成分を除去して得られる画像への上記リンギングを示す成分の混入をより確実に抑制でき、上記画像の品質を高めることができる。これとともに、スジムラ成分の除去に伴う画像中からの線状画像情報成分の消失をより確実に抑制することもできる。

なお、上記絶対値を小さくする処理が施された画像データに対応する画像中の領域においては、この処理によりスジムラ成分を示す画像情報が消滅したり不正確なものとなったりするので、この領域においては上記スジムラ成分の除去を正確に実施することはできないが、これらの領域は画像データの値が急激に変化する領域なので、上記領域中に除去されずに残ったスジムラ成分に起因する画像品質の劣化は極めて少なく、このノイズ成分が上記画像の品質を高める効果を妨げることはない。

本発明の第３の画像処理方法、および他の画像処理装置ならびに他の画像処理プログラムによれば、上記光電変換素子列で検出して得られた画像担持体が担持する画像を示す画像データにフィルタ処理を施して、上記画像データが示す画像中の副走査方向および／または主走査方向に延びる線状のノイズ成分を示す画像データを抽出し、上記線状のノイズ成分を示す各画像データの値の絶対値を全体的に小さくした画像データを得、上記画像データを上記画像を示す画像データから差し引いてこの画像から線状のノイズ成分を除去するようにしたので、スジムラ成分を除去して得られる画像への上記リンギングを示す成分の混入を抑制できるとともに、スジムラ成分の除去に伴う画像中からの線状画像情報成分の消失を抑制することもでき、上記画像の品質を高めることができる。

以下、本発明の実施の形態について、図面を用いて説明する。図１は本発明実施の形態による画像処理装置の概略構成を示す図、図２（ａ）はラインセンサで検出され取得された画像データで表わされる画像を示す図であり、図２（ｂ）はラインセンサで検出され取得された画像データの値を示す図である。

図示の本発明の実施の形態による画像処理装置１００は、多数の受光部１１を主走査方向（図１中の矢印Ｘ′方向）に並べてなる光電変換素子列であるラインセンサ１０に対して、画像担持体である原稿１を主走査方向と交わる副走査方向（図１中の矢印Ｙ′方向）に相対的に移動させながら、照明されてこの原稿１から発せられた光をラインセンサ１０で検出して得られたこの原稿１が担持する画像である原画像２０Ｅを示す画像データに対して画像処理を施すものである。この画像処理装置１００は、上記原画像２０Ｅを示す画像データにフィルタ処理を施して、この画像データが示す原画像２０Ｅ中の副走査方向（図２中の矢印Ｙ方向）に延びる線状のノイズ成分（以後、スジムラ成分２１Ｅという）を示す画像データを抽出する線状ノイズ抽出手段３０と、上記線状ノイズ抽出手段３０が原画像２０Ｅを示す画像データにフィルタ処理を施す前に、この画像データに対して、原画像２０Ｅ中の上記副走査方向と直交する副走査直交方向（図２中の矢印Ｘ方向）における画像データの値の急激な変化を緩やかにする前処理を実施する前処理手段３１と、上記スジムラ成分２１Ｅを示す画像データを原画像２０Ｅを示す画像データから差し引く差引手段４０とを備えている。上記原画像２０Ｅ中の副走査方向（図２中の矢印Ｙ方向）は、上記ラインセンサ１０に対して原稿１を相対的に移動させる副走査方向（図１中の矢印Ｙ′方向）に対応する方向である。

上記線状ノイズ抽出手段３０は、上記原画像２０Ｅを示す画像データにフィルタ処理を施して、この画像データが示す原画像２０Ｅ中の主走査方向に延びるスジムラ成分を示す画像データを抽出することもでき、上記前処理手段３１は、線状ノイズ抽出手段３０が原画像２０Ｅを示す画像データにフィルタ処理を施す前に、この画像データに対して、原画像２０Ｅ中の上記主走査方向と直交する主走査直交方向における画像データの値の急激な変化を緩やかにする前処理を実施することもできる。上記原画像２０Ｅ中の主走査方向（図示は省略）は、上記ラインセンサ１０において多数の受光部１１が並べられている主走査方向（図１中の矢印Ｘ′方向）に対応する方向である。

なお、上記スジムラ成分２１Ｅは、主に、ラインセンサ１０を構成する各受光部毎の特性の違いに起因して生じるものであり、実際には、原画像２０Ｅ中に副走査方向へ延びる大小様々な多数のスジムラが現れることがあるが、図２（ａ）においては上記原画像２０Ｅ中の２本の大きなスジムラ成分２１Ｅのみを示した。また、図２（ｂ）に示す画像データの値は、原画像２０Ｅにおいて上記副走査直交方向に延びる線ｈ−ｈ′上に位置する画像データの値を示している。

また、上記原画像２０Ｅを示す画像データは、画像中の副走査直交方向において画像データの値が急激に変化する領域を有するものである。

後に詳しく説明するが、上記前処理手段３１は、原画像２０Ｅ中の上記副走査直交方向において互いに隣り合う上記原画像２０Ｅを示す画像データＤｅの差分値を示す画像データを得、上記差分値を示す画像データに対して、これらの画像データのうちの所定の閾値を越える値を持つ画像データの値の絶対値を小さくする処理を施して得られた画像データを副走査直交方向に累積加算したり、あるいは、原画像２０Ｅ中の副走査直交方向において互いに隣り合う原画像２０Ｅを示す画像データＤｅの差分値を示す画像データを得、これらの差分値を示す画像データを副走査直交方向に移動平均し、上記差分値を示す画像データに対して、上記移動平均した画像データのうちの所定の閾値を越える値を持つ画像データに対応する上記差分値を示す画像データの値の絶対値を小さくする処理を施して得られた画像データを副走査直交方向に累積加算するものである。また、上記前処理手段３１は、原画像２０Ｅ中の上記主走査直交方向において互いに隣り合う上記原画像２０Ｅを示す画像データＤｅの差分値を示す画像データを得、上記差分値を示す画像データに対して、これらの画像データのうちの所定の閾値を越える値を持つ画像データの値の絶対値を小さくする処理を施して得られた画像データを主走査直交方向に累積加算したり、あるいは、原画像２０Ｅ中の主走査直交方向において互いに隣り合う原画像２０Ｅを示す画像データＤｅの差分値を示す画像データを得、これらの差分値を示す画像データを主走査直交方向に移動平均し、上記差分値を示す画像データに対して、上記移動平均した画像データのうちの所定の閾値を越える値を持つ画像データに対応する上記差分値を示す画像データの値の絶対値を小さくする処理を施して得られた画像データを主走査直交方向に累積加算することもできる。

なお、上記移動平均した画像データのうちの所定の閾値を越える値を持つ画像データと、この画像データに対応する上記差分値を示す画像データとは上記原画像２０Ｅ中の同一位置を示すものである。また、上記移動平均は周波数処理の１例として示したものである。

次に、上記画像処理装置の作用について説明する。

以下、上記線状ノイズ抽出手段でのスジムラ成分を抽出する処理で生じるリンギングを、前処理手段３１による前処理で抑制する作用について説明する。

まず始めに、比較例として、原画像２０Ｅ中の副走査方向に延びるスジムラ成分２１Ｅを示す画像データを、この原画像２０Ｅを示す画像データに対して上記前処理を施すことなく抽出する場合について説明する。図３は前処理を実施することなく線状ノイズ抽出手段でスジムラ成分を抽出する様子を示す図であり、図３（ａ−１）は上記原画像を示す図、図３（ａ―２）は上記原画像を示す画像データの値を示す図、図３（ｂ−１）は原画像を示す画像データにフィルタ処理を施して得られた画像データで表わされる画像を示す図、図３（ｂ―２）は上記原画像を示す画像データにフィルタ処理を施して得られた画像データの値を示す図、図３（ｃ−１）上記原画像から抽出されたスジムラ成分を示す図、図３（ｃ―２）は上記原画像から抽出されたスジムラ成分を表す画像データの値を示す図である。なお、上記各画像データの値を示す図は、各画像中の副走査直交方向に延びるｈ−ｈ′線上に位置する画像データの値を示すものである。

上記原画像２０Ｅを示す画像データＤｅに対して上記の前処理を実施することなく、フィルタ処理である上記副走査直交方向へのローパスフィルタ処理を施して低周波成分からなる画像２０Ｆを示す画像データＤｆを得る。このフィルタ処理が施された画像２０Ｆはボケて広がる画像となり、画像データＤｅの値が上記副走査直交方向に急激に変化する原画像２０Ｅ中の領域Ｒｅに対応する上記画像２０Ｆ中の領域Ｒｆの近傍を示す画像データＤｆの値が、上記画像データＤｅの値から大きく変化する（図３（ｂ−１）および図３（ｂ―２）参照）。このような画像データＤｆを、原画像２０Ｅを示す画像データＤｅから差し引いて高周波成分からなる画像データＤｇを取得すると、この画像データＤｇが示す画像２０Ｇは、スジムラ成分２１Ｇの他に上記領域Ｒｅに対応する領域Ｒｇの近傍に上記画像データＤｅの値と画像データＤｆの値との差を示す上記副走査直交方向に波打つ画像（リンギングが生じた画像）となってしまう（図３（ｃ−１）および図３（ｃ―２）参照）。

このような場合には、上記フィルタ処理を実施する前に、上記前処理手段３１による前処理を実施することで、上記原画像２０Ｅからスジムラを示す画像データを上記リンギングの発生なく抽出することができる。以下、前処理を実施してスジムラ成分を示す画像データを抽出する場合について説明する。図４は画像データに前処理を施す様子を示す図であり、図４（ａ）は上記原画像を表す画像データを示す図、図４（ｂ）は上記画像データの差分値を示す画像データ、図４（ｃ）は上記差分値を示す画像データを移動平均して得られた画像データを示す図、図４（ｄ）は上記移動平均した画像データのうちの所定の閾値以上の値を持つ画像データに対応する上記差分値を示す画像データの値を０にして得られた画像データを示す図である。図５は上記差分値を示す画像データの値を０にした画像データを累積加算した画像データを示す図であり、図５（ａ）上記累積加算した画像データで表わされる画像を示す図、図５（ｂ）は累積加算した画像データを示す図である。

上記前処理手段３１で実施する前処理は、始めに、上記原画像２０Ｅ中の副走査直交方向において互いに隣り合うこの原画像２０Ｅを示す画像データＤｅの差分値を示す画像データＳｅを得る。具体的には、例えば、図４（ａ）、（ｂ）に示すように、画像データＤｅ中の互いに隣り合う画像データである画像データＤｅ１、Ｄｅ２について、画像データＤｅ２から画像データＤｅ１を差し引いた値である差分値Ｕ１を得、この差分値Ｕ１を原画像２０Ｅ中において画像データＤｅ１が示す位置に対応する画像データＳｅ１の値とする。また、画像データＤｅ中の互いに隣り合う画像データである画像データＤｅ２、Ｄｅ３について、画像データＤｅ中の画像データＤｅ３から画像データＤｅ２を差し引いた値である差分値Ｕ２を得、この差分値Ｕ２を原画像２０Ｅ中において画像データＤｅ２が示す位置に対応する画像データＳｅ２の値とする。このようにして上記差分値の全体を示す画像データＳｅを得る。

次に、上記差分値を示す画像データＳｅを副走査直交方向に移動平均し、移動平均した画像データＨｅを得る。具体的には、例えば、画像データＳｅ１、Ｓｅ２、Ｓｅ３の３つの画像データの値を加算平均した値Ｖ２を得、この値Ｖ２を原画像２０Ｅ中において上記画像データＳｅ２が示す位置に対応する画像データＨｅ２の値とする。このようにして上記移動平均した全体を示す画像データＨｅを得る。

つづいて、上記差分値を示す画像データＳｅに対して、上記移動平均した画像データＨｅのうちの所定の閾値±Ｋｈを越える値を持つ画像データに対応するこの画像データＳｅ中の画像データの値の絶対値を小さくする処理を施して画像データＴｅを得る。具体的には、例えば、上記正の値を持つ閾値＋Ｋｈを越える値を持つ画像データである、画像データＨｅ２、画像データＨｅ３が示す原画像２０Ｅ中の位置に対応する画像データＴｅ２、画像データＴｅ３の値を、例えば０にする。このようにして上記画像データＴｅを得る。すなわち、画像データＴｅは、画像データＤｅのうちの副走査直交方向に値の変化が激しい領域Ｊ１、Ｊ２に対応する画像データの値を０にしたものである。

なお、上記閾値±Ｋｈは、画像２０Ｅ中の画像データＤｅの値が急激に変化する領域である上記スジムラ成分を抽出する処理を施したときにリンギングが大きい領域Ｊ１、Ｊ２を示す画像データと他の領域を示す画像データとを分離できるように定められたものである。また、上記閾値の設定によってスジムラ成分と線状画像情報成分との分離もなされる。

次に、画像データＴｅを上記副走査直交方向に累積加算して画像データＤｅ′を得る。具体的には、例えば、画像データＴｅを、副走査直交方向の一方の側（ｈ−ｈ′線におけるｈの側）から順次加算して（累積加算して）求められた画像データＤｅ１′の値に画像データＴｅ１の値を加算して画像データＤｅ２′の値を求め、画像データＤｅ２′の値に画像データＴｅ２の値を加算して画像データＤｅ３′の値を求めるという演算を繰り返して、上記累積加算した全体を示す画像データＤｅ′を得る。

ここで得られた画像データＤｅ′は、上記原画像２０Ｅ中の副走査直交方向における画像データＤｅの値の急激な変化を緩やかにしたものである。図５（ａ）に上記画像データＤｅ′が表す画像２０Ｅ′を示す。また、図５（ｂ）に示す画像データＤｅ′の値中の領域Ｊ１、Ｊ２が、上記画像データＤｅの値の急激な変化を緩やかにした領域である。

なお、上記前処理において移動平均の処理を実施しない場合であっても、上記原画像２０Ｅ中の副走査直交方向における画像データＤｅの値の急激な変化を緩やかにすることができる。すなわち、図４（ｂ）に示すように、上記差分値を示す画像データＳｅに対して、この画像データのうちの所定の閾値±Ｋｓを越える値を持つ画像データについて、その値の絶対値を小さくする処理を施して得られた画像データを上記と同様に副走査直交方向に累積加算して、画像中の副走査直交方向における画像データの値の急激な変化を緩やかにした画像データを作成するようにしてもよい。

なお、上記閾値は、画像データＳｅについて閾値を設定する場合と、画像データＨｅについて閾値を設定する場合とでは、上記スジムラ成分を抽出する処理でリンギングが大きい領域と他の領域とを分離する分離性能が異なる。すなわち、一般に画像データが示す成分はスジムラ成分（高周波成分）より低周波成分が多い。移動平均は高周波成分を減衰させるため、スジムラ成分を分離する閾値を小さくすることができ、画像データ中からスジムラ成分を分離させる分離性能が良くなる。

最後に、この画像データＤｅ′に対して副走査直交方向にフィルタ処理（例えばローパスフィルタ処理）を施して上記スジムラ成分２１Ｅ′を示す画像データを抽出する。図６は画像データにフィルタ処理を施してスジムラ成分を抽出する様子を示す図であり、図６（ａ−１）は前処理が施された画像データで表わされる画像を示す図、図６（ａ−２）は前処理が施された画像データの値を示す図、図６（ｂ−１）は前処理が施された画像データにフィルタ処理を施して得られた画像データで表わされる画像を示す図、図６（ｂ−２）は前処理が施された画像データにフィルタ処理を施して得られた画像データの値を示す図、図６（ｃ−１）は抽出されたスジムラ成分を示す図、図６（ｃ−２）は抽出されたスジムラ成分を表す画像データの値を示す図である。なお、上記各画像データの値を示す図は、各画像中の副走査直交方向に延びるｈ−ｈ′線上に位置する画像データの値を示すものである。

線状ノイズ抽出手段３０により、上記画像２０Ｅ′を示す画像データＤｅ′に上記フィルタ処理である副走査直交方向（図中矢印Ｘ方向）へのローパスフィルタ処理を施して低周波成分からなる画像２０Ｂを示す画像データＤｂを得、上記画像データＤｅ′から画像データＤｂを差し引いて高周波成分からなる画像２０Ｃを示す画像データＤｃを取得する。この画像２０Ｃがスジムラ成分２１Ｅ′を示す画像となる。なお、上記スジムラ成分を示す画像データを抽出する処理は、画像データＤｅ′に対して上記ローパスフィルタ処理を行なった後、さらに副走査方向にもローパスフィルタ処理を行なって、低周波成分からなる画像２０Ｂ′を示す画像データＤｂ′を取得するようにしてもよい。

上記のように、前処理で副走査直交方向における画像データの値の変化が緩やかになった上記画像データＤｅ′に対してスジムラ成分を抽出する処理を実施してスジムラ成分を示す画像データを抽出することにより、上記前処理を実施しない画像データＤｅに対してフィルタ処理を実施する比較例に比して、リンギングの発生を抑えてスジムラ成分を示す画像データを抽出することができる。

次に、差引手段により画像からスジムラ成分を除去する。図７は画像からスジムラ成分を除去する様子を示す図であり、図７（ａ−１）は上記原画像を示す図、図７（ａ−２）は上記原画像を表す画像データの値示す図、図７（ｂ−１）は抽出されたスジムラ成分を示す図、図７（ｂ−２）は抽出されたスジムラ成分を表す画像データの値を示す図、図７（ｃ−１）は上記原画像からスジムラ成分を除去した画像を示す図、図７（ｃ−２）は上記原画像を示す画像データからスジムラ成分を除去した画像データの値を示す図である。なお、上記各画像データの値を示す図は、各画像中の副走査直交方向に延びるｈ−ｈ′線上に位置する画像データの値を示すものである。

差引手段４０により、原画像２０Ｅを示す画像データＤｅからスジムラ成分２１Ｅ、すなわち原画像２０Ｅ中の高周波成分を示す画像２０Ｃを示す画像データＤｃを差し引いて、原画像２０Ｅからスジムラ成分２１Ｅを除去した画像２０Ｄを示す画像データＤｄを得る。上記のことにより、スジムラ成分の除去に伴う画像中からの線状画像情報成分の消失を抑制することができる。

このとき、例えば、画像データＤｅから差し引くスジムラ成分２１Ｅを示す画像データＤｃの各値をそのまま差し引くと、スジムラ成分２１Ｅを含む高周波成分が画像データＤｅから必要以上に差し引かれて、原画像２０Ｅの全体的な質感と、スジムラ成分を除去した画像２０Ｄの全体的な質感とに差が生じることがある。したがって、原画像２０Ｅからスジムラ成分２１を除去する際には、上記画像２０Ｄの全体的な質感の変化が少なくなるように、上記抽出されたスジムラ成分を示す画像データＤｃの各値に、例えば一定の係数を掛け合わて、画像データＤｃの各値の大きさを変更した後、上記原画像２０Ｅを示す画像データＤｅから上記のようにして作成されたスジムラ成分を示す画像データを差し引いて、原画像２０Ｅからスジムラ成分を除去することが望ましい。

上記質感の差をフーリエ空間において説明すると以下のようになる。図８は画像の２次元スペクトルを示す図であり、図８（ａ）は原画像の２次元スペクトルを示す図、図８（ｂ）はスジムラ成分が除去された画像の２次元スペクトルを示す図である。

図８（ａ）に示すように上記スジムラ成分２１Ｅを含む原画像２０Ｅを２次元フーリエ変換して得られた２次元スペクトル２５Ｅには、上記スジムラ成分２１Ｅが、副走査方向（図中矢印Ｙ方向）おいては低周波成分に、かつ副走査直交方向（図中矢印Ｘ方向）においては高周波数成分を示す領域２６Ｅに現れる。例えば、画像データＤｅから差し引くスジムラ成分２１Ｅを含む画像データＤｃの各値をそのまま差し引くと、ここで得られた画像を２次元フーリエ変換して得られた２次元スペクトル２５Ｅ′中の上記領域２６Ｅに対応する領域２６Ｅ´には、上記スジムラ成分２１Ｅを含む上記スペクトル成分が殆どなくなってしまう。すなわち、上記スジムラ成分が除去された画像において、２次元スペクトル２５Ｅ′の領域２６Ｅ´で示される特定のスペクトル成分が極端に少なくなってしまう。そのため、上記画像２０Ｄの質感と原画像２０Ｅの質感とに差が生じることがある。

また、上記画像データに対して施すフィルタ処理を、高速で行なうため、図９に示すように、図中フィルタＦｅを、互いに移動平均する範囲が異なるフィルタＦｅ１、Ｆｅ２、Ｆｅ３、Ｆｅ４、Ｆｅ５で近似し、上記画像データをフィルタＦｅに通す処理を、この画像データをそれぞれのフィルタＦｅ１、Ｆｅ２、Ｆｅ３、Ｆｅ４、Ｆｅ５に通す個別の処理に分けて実行することにより、上記画像データをフィルタＦｅに通す処理と同等な処理をより少ない演算量で実施することができる。すなわち、上記画像データをフィルタＦｅに通す際の乗算処理を、この画像データを各フィルタＦｅ１、Ｆｅ２、Ｆｅ３、Ｆｅ４、Ｆｅ５に通す際の加算処理で代替することにより全体の演算量を少なくすることができる。

なお、上記前処理は、上記原画像を示す画像データを、スジムラが延びる副走査方向における高周波数成分あるいは低周波数成分を低減するフィルタ等を通した後に実施するようにしてもよい。

また、上記スジムラ成分を抽出する処理は、前処理が施された上記スジムラ成分を含む画像データをハイパスフィルタ、すなわち、上記副走査直交方向における高周波数成分を抽出するフィルタに通して行なうようにしてもよい。

また、上記差引手段により原画像からのスジムラ成分を除去する際の、上記スジムラ成分を示す画像データの各値の大きさの変更は、予め用意された変換テーブルや変換関数を用いて実施するようにしてもよい。

なお、上記実施の形態においては、スジムラ成分が画像中の副走査方向に延びるものである場合について説明したが、スジムラ成分が画像中の主走査方向に延びるものである場合、およびスジムラ成分が画像中の副走査方向および主走査方向に延びるものである場合についても、上記と同様に、リンギングの発生を抑制してスジムラ成分を示す画像データを抽出することができ、スジムラ成分を除去した画像を得る際に、上記リンギングを示す成分の画像への混入を抑制することができることができることは言うまでもない。

また、上記リンギングの発生を抑制してスジムラ成分を示す画像データを抽出する手法は、主走査方向あるいは副走査方向以外の方向に延びるスジムラ成分に対しても適用することができる。

なお、上記画像担持体あるいは原稿は、基板上に蓄積性蛍光体層を積層して作成された放射線像の記録に使用される放射線像変換パネルであってもよい（例えば、特開２００１−１３５９９号公報参照）。このような場合には、励起光の照射を受けて放射線像変換パネルから発生した輝尽発光光をラインセンサで検出することにより、上記放射線像変換パネルが担持する放射線像を示す画像データを得ることができる。ここで、上記ラインセンサは千鳥状に配置することもできる。

また、上記画像処理装置の機能を備えたプログラムをＣＤ−ＲＯＭ等の記録媒体に記録しコンピュータにインストールして、上記と同様の効果が得られるように機能させるようにしてもよい。

なお、上記原画像２０Ｅを示す画像データＤｅに前処理を施すことなくこの画像中からスジムラ成分を除去する手法として、以下のような手法が考えられる。

すなわち、図１０および図１１に示すように、線状ノイズ抽出手段３０により、上記スジムラ成分２１Ｅを含む原画像２０Ｅを示す画像データＤｅ（図１０（ａ）参照）にフィルタ処理を施して画像データＤｅから、予め定められた高周波成分の抽出を行なって上記スジムラ成分２１Ｅを示す画像データＷ１を抽出する（図１０（ｂ）参照）。

次に、例えば、画像データＤｅから上記スジムラ成分２１Ｅを示す画像データＷ１の各値をそのまま差し引くと、リンギングの成分およびスジムラ成分２１Ｅを含む高周波成分が画像データＤｅから必要以上に差し引かれて、原画像２０Ｅの全体的な質感とスジムラ成分が除去された画像の全体的な質感とに差が生じることがあるので、データ変換手段３５により、上記抽出されたスジムラ成分を示す画像データＷ１の各値に、例えば一定の係数を掛け合わて、画像データＷ１の各値の絶対値の大きさを一律に小さくした画像データＷ２（図１０（ｃ）参照）を得る。

その後、差引手段４０により上記原画像２０Ｅを示す画像データＤｅから上記のようにして作成されたスジムラ成分を示す画像データＷ２を差し引いて、原画像２０Ｅからスジムラ成分を除去した画像を示す画像データＷ３（図１０（ｄ）参照）を得る。これにより、原画像からスジムラ成分を除去する際の、リンギングの振幅、画像の全体的な質感の変化を少なくすることができる。これとともに、スジムラ成分の除去に伴う画像中からの線状画像情報成分の消失を抑制することもできる。

本発明の実施の形態による画像処理装置の概略構成を示す斜視図原画像および原画像を表す画像を示す図前処理を施すことなく線状ノイズ抽出手段でスジムラ成分を抽出する様子を示す図原画像を示す画像データに前処理を施す様子を示す図累積加算した画像データおよびこの画像データで表わされる画像を示す図前処理が施された画像データにフィルタ処理を施してスジムラ成分を抽出する様子を示す図原画像からスジムラ成分を除去する様子を示す図画像の２次元スペクトルを示す図重み付け移動平均を行う量子化したフィルタを示す図前処理を施すことなく画像中からスジムラ成分を除去する手法を示す図前処理を施すことなく画像中からスジムラ成分を除去する画像処理装置の概略構成を示す斜視図

符号の説明

１画像担持体
１０ラインセンサ
１１受光部
２０Ｅ原画像
２１Ｅスジムラ成分
３０線状ノイズ抽出手段
３１前処理手段
４０差引手段
１００画像処理装置

Claims

画像データに、該画像データが示す画像中の１方向に延びる線状のノイズ成分を示す画像データを抽出するフィルタ処理を施して、該線状のノイズ成分を示す画像データを前記画像を示す画像データから差し引いて前記画像から前記線状のノイズ成分を除去する画像処理方法において、
前記画像を示す画像データに前記フィルタ処理を施す前に、該画像データに対して、該画像中の前記１方向と直交する直交方向における画像データの値の急激な変化を緩やかにするように、前記画像中の前記直交方向において互いに隣り合う該画像を示す画像データの差分値を示す画像データに対して、該差分値を示す画像データのうちの所定の閾値を越える値を持つ画像データの値の絶対値を小さくする処理を施して得られた画像データを前記直交方向に累積加算する前処理を施すことを特徴とする画像処理方法。
多数の受光部を主走査方向に並べてなる光電変換素子列に対して、画像担持体を前記主走査方向と交わる副走査方向に相対的に移動させながら、該画像担持体から発せられた光を前記光電変換素子列で検出して得られたこの画像担持体が担持する画像を示す画像データに、該画像データが示す画像中の前記副走査方向に延びる線状のノイズ成分を示す画像データを抽出するフィルタ処理を施して、該線状のノイズ成分を示す画像データを前記画像を示す画像データから差し引いて前記画像から前記線状のノイズ成分を除去する画像処理方法において、
前記画像を示す画像データに前記フィルタ処理を施す前に、該画像データに対して、該画像中の前記副走査方向と直交する副走査直交方向における画像データの値の急激な変化を緩やかにするように、前記画像中の前記副走査直交方向において互いに隣り合う該画像を示す画像データの差分値を示す画像データに対して、該差分値を示す画像データのうちの所定の閾値を越える値を持つ画像データの値の絶対値を小さくする処理を施して得られた画像データを前記副走査直交方向に累積加算する前処理を施すことを特徴とする画像処理方法。
多数の受光部を主走査方向に並べてなる光電変換素子列に対して、画像担持体を前記主走査方向と交わる副走査方向に相対的に移動させながら、該画像担持体から発せられた光を前記光電変換素子列で検出して得られたこの画像担持体が担持する画像を示す画像データに、該画像データが示す画像中の前記主走査方向に延びる線状のノイズ成分を示す画像データを抽出するフィルタ処理を施して、該線状のノイズ成分を示す画像データを前記画像を示す画像データから差し引いて前記画像から前記線状のノイズ成分を除去する画像処理方法において、
前記画像を示す画像データに前記フィルタ処理を施す前に、該画像データに対して、該画像中の前記主走査方向と直交する主走査直交方向における画像データの値の急激な変化を緩やかにするように、前記画像中の前記主走査直交方向において互いに隣り合う該画像を示す画像データの差分値を示す画像データに対して、該差分値を示す画像データのうちの所定の閾値を越える値を持つ画像データの値の絶対値を小さくする処理を施して得られた画像データを前記主走査直交方向に累積加算する前処理を施すことを特徴とする画像処理方法。
多数の受光部を主走査方向に並べてなる光電変換素子列に対して、画像担持体を前記主走査方向と交わる副走査方向に相対的に移動させながら、該画像担持体から発せられた光を前記光電変換素子列で検出して得られたこの画像担持体が担持する画像を示す画像データに、該画像データが示す画像中の前記副走査方向に延びる線状のノイズ成分を示す画像データを抽出するフィルタ処理を施す線状ノイズ抽出手段と、該線状のノイズ成分を示す画像データを前記画像を示す画像データから差し引く差引手段とを備え、前記画像から前記線状のノイズ成分を除去する画像処理装置において、
前記線状ノイズ抽出手段が前記画像を示す画像データに前記フィルタ処理を施す前に、該画像データに対して、該画像中の前記副走査方向と直交する副走査直交方向における画像データの値の急激な変化を緩やかにするように、前記画像中の前記副走査直交方向において互いに隣り合う該画像を示す画像データの差分値を示す画像データに対して、該差分値を示す画像データのうちの所定の閾値を越える値を持つ画像データの値の絶対値を小さくする処理を施して得られた画像データを前記副走査直交方向に累積加算する前処理手段を備えたものであることを特徴とする画像処理装置。
前記所定の閾値を超える値を持つ画像データが、前記差分値を示す画像データを前記副走査直交方向に移動平均した画像データのうちの所定の閾値を超える値を持つ画像データに対応する前記差分値を示す画像データであることを特徴とする請求項４記載の画像処理装置。
多数の受光部を主走査方向に並べてなる光電変換素子列に対して、画像担持体を前記主走査方向と交わる副走査方向に相対的に移動させながら、該画像担持体から発せられた光を前記光電変換素子列で検出して得られたこの画像担持体が担持する画像を示す画像データに、該画像データが示す画像中の前記主走査方向に延びる線状のノイズ成分を示す画像データを抽出するフィルタ処理を施す線状ノイズ抽出手段と、該線状のノイズ成分を示す画像データを前記画像を示す画像データから差し引く差引手段とを備え、前記画像から前記線状のノイズ成分を除去する画像処理装置において、
前記線状ノイズ抽出手段が前記画像を示す画像データに前記フィルタ処理を施す前に、該画像データに対して、該画像中の前記主走査方向と直交する主走査直交方向における画像データの値の急激な変化を緩やかにするように、前記画像中の前記主走査直交方向において互いに隣り合う該画像を示す画像データの差分値を示す画像データに対して、該差分値を示す画像データのうちの所定の閾値を越える値を持つ画像データの値の絶対値を小さくする処理を施して得られた画像データを主走査直交方向に累積加算する前処理手段を備えたものであることを特徴とする画像処理装置。
前記所定の閾値を超える値を持つ画像データが、前記差分値を示す画像データを前記主走査直交方向に移動平均した画像データのうちの所定の閾値を超える値を持つ画像データに対応する前記差分値を示す画像データであることを特徴とする請求項６記載の画像処理装置。
多数の受光部を主走査方向に並べてなる光電変換素子列に対して、画像担持体を前記主走査方向と交わる副走査方向に相対的に移動させながら、該画像担持体から発せられた光を前記光電変換素子列で検出して得られたこの画像担持体が担持する画像を示す画像データに、該画像データが示す画像中の前記副走査方向に延びる線状のノイズ成分を示す画像データを抽出するフィルタ処理を施す手順と、該線状のノイズ成分を示す画像データを前記画像を示す画像データから差し引いて前記画像から前記線状のノイズ成分を除去する手順とをコンピュータに実行させるためのプログラムにおいて、
前記画像を示す画像データに前記フィルタ処理を施す前に、該画像データに対して、該画像中の前記副走査方向と直交する副走査直交方向における画像データの値の急激な変化を緩やかにするように、前記画像中の前記副走査直交方向において互いに隣り合う該画像を示す画像データの差分値を示す画像データに対して、該差分値を示す画像データのうちの所定の閾値を越える値を持つ画像データの値の絶対値を小さくする処理を施して得られた画像データを前記副走査直交方向に累積加算する前処理を行なわせる手順を実行させることを特徴とする画像処理をコンピュータに実行させるためのプログラム。
多数の受光部を主走査方向に並べてなる光電変換素子列に対して、画像担持体を前記主走査方向と交わる副走査方向に相対的に移動させながら、該画像担持体から発せられた光を前記光電変換素子列で検出して得られたこの画像担持体が担持する画像を示す画像データに、該画像データが示す画像中の前記主走査方向に延びる線状のノイズ成分を示す画像データを抽出するフィルタ処理を施す手順と、該線状のノイズ成分を示す画像データを前記画像を示す画像データから差し引いて前記画像から前記線状のノイズ成分を除去する手順とをコンピュータに実行させるためのプログラムにおいて、
前記画像を示す画像データに前記フィルタ処理を施す前に、該画像データに対して、該画像中の前記主走査方向と直交する主走査直交方向における画像データの値の急激な変化を緩やかにするように、前記主走査直交方向において互いに隣り合う該画像を示す画像データの差分値を示す画像データに対して、該差分値を示す画像データのうちの所定の閾値を越える値を持つ画像データの値の絶対値を小さくする処理を施して得られた画像データを前記主走査直交方向に累積加算する前処理を行なわせる手順を実行させることを特徴とする画像処理をコンピュータに実行させるためのプログラム。