JP2004134861A - Resolution evaluation method, resolution evaluation program, and optical apparatus - Google Patents
Resolution evaluation method, resolution evaluation program, and optical apparatus Download PDFInfo
- Publication number
- JP2004134861A JP2004134861A JP2002294967A JP2002294967A JP2004134861A JP 2004134861 A JP2004134861 A JP 2004134861A JP 2002294967 A JP2002294967 A JP 2002294967A JP 2002294967 A JP2002294967 A JP 2002294967A JP 2004134861 A JP2004134861 A JP 2004134861A
- Authority
- JP
- Japan
- Prior art keywords
- lines
- edge
- resolution
- resolution evaluation
- luminance
- 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.)
- Pending
Links
Images
Landscapes
- Facsimiles In General (AREA)
- Testing, Inspecting, Measuring Of Stereoscopic Televisions And Televisions (AREA)
Abstract
Description
【0001】
【発明の属する技術分野】
本発明は、電子カメラやスキャナ等の光学機器の解像度評価方法に関する。
【0002】
【従来の技術】
従来、電子カメラやスキャナ等の光学機器の第1の解像度評価方法として、低輝度部分と高輝度部分とが相互に配置された縞状の原稿を撮影する方法が知られている(例えば、特許文献1参照。)。より具体的には、縞の幅を相互に違える複数の原稿を撮影し、このデータを解析することにより解像度評価が行われている。この解像度は、縞の幅が狭くなればなる程、低輝度部分と高輝度部分との区別がつきにくくなるため解像度が落ちることになる。
【0003】
第2の解像度評価方法としては、低輝度部分と高輝度部分との境界部分(以下、エッジという)を有する原稿を予め用意し、この原稿を撮影する方法が知られている(例えば、非特許文献1参照。)。
【特許文献1】
特開1999−313186号公報
【非特許文献1】
”Photography−Electronic still picture cameras−resolution measurements”, Draft International Standards, 米国, International Organization for Standarization, 1997.11.18, ISO/TC42N 4231
【0004】
【発明が解決しようとする課題】
しかしながら、上述した第1の解像度評価方法では、複数の原稿を一度に撮影するか、複数回撮影しなければならなかった。複数の原稿を一度に撮影する場合は、1つの画像内で部分的な解像度の比較が困難であり、複数回撮像する場合は原稿の配置を幾度も変更しなけらばならず手間がかかるという問題があった。
【0005】
かかる問題は、第2の解像度評価方法では発生せず、第1の解像度評価方法よりも第2の解像度評価方法の方が一般的である。しかしながら、この第2の解像度評価方法では、光学機器の機械的振動により測定データが不正確になるという問題があった。
本発明は、かかる従来の問題を解決するためになされたもので、光学機器の機械的振動による影響を受けにくい解像度評価方法を提供することを目的とする。
【0006】
【課題を解決するための手段】
請求項1に記載の解像度評価方法は、原稿に描かれた模様のエッジを選択するステップと、前記模様のエッジ近傍における輝度の分布を複数ライン求めるステップと、この複数ラインのライン毎に、エッジ座標を求めるステップと、前記複数ラインのうち少なくとも1つのラインにおける前記エッジ座標でのエッジの傾きを求めるステップと、前記複数ラインのライン毎に、前記ライン状に計測された輝度の分布から所定間隔毎に計測データを選択し、相隣接する前記計測データの差分を求め、離散した差分データを得るステップと、この離散した差分データをフーリエ変換し複数の変調伝達関数を求めるステップと、この複数の変調伝達関数に基づき解像度を評価するステップとを有することを特徴とする。
【0007】
請求項2に記載の解像度評価方法は、原稿に描かれた模様のエッジを選択するステップと、前記模様のエッジ近傍における輝度の分布を複数ライン求めるステップと、この複数ラインのライン毎に、エッジ座標を求めるステップと、前記複数ラインのうち少なくとも1つのラインにおける前記エッジ座標でのエッジの傾きを求めるステップと、前記複数ラインのライン毎に、前記ライン状に計測された輝度の分布から所定間隔毎に計測データを選択し、相隣接する前記計測データの差分を求め、離散した差分データを得るステップと、前記複数ラインのライン毎に、前記離散した差分データに近似曲線をフィッティングさせるステップと、この複数の近似曲線をフーリエ変換し変調伝達関数を求めるステップと、この複数の変調伝達関数に基づき解像度を評価するステップとを有することを特徴とする。
【0008】
請求項3に記載の解像度評価方法は、請求項1または請求項2に記載の解像度評価方法において、前記エッジ座標を算出するために、前記被写体のエッジ近傍における低輝度部分の平均値と高輝度部分の平均値とを求め、この低輝度部分の平均値と高輝度部分の平均値とからエッジ閾値を求め、前記計測データを前記ラインの一方から検査して前記輝度がエッジ閾値を跨る位置を検出することを特徴とする。
【0009】
請求項4に記載の解像度評価方法は、請求項1ないし請求項3のいずれか1項に記載の解像度評価方法において、前記解像度の評価を、前記複数の変調伝達関数の平均値に基づいて行うことを特徴とする。
請求項5に記載の解像度評価方法は、請求項1ないし請求項3のいずれか1項に記載の解像度評価方法において、前記解像度の評価を、評価する周波数毎の前記複数の変調伝達関数の最大値に基づいて行うことを特徴とする。
【0010】
請求項6に記載の解像度評価方法は、請求項1ないし請求項3のいずれか1項に記載の解像度評価方法において、前記解像度の評価を、評価する周波数毎の前記複数の変調伝達関数の最小値に基づいて行うことを特徴とする。
請求項7に記載の解像度評価プログラムは、原稿に描かれた模様のエッジを選択する手順と、前記模様のエッジ近傍における輝度の分布を複数ライン求める手順と、この複数ラインのライン毎に、エッジ座標を求める手順と、前記複数ラインのうち少なくとも1つのラインにおける前記エッジ座標でのエッジの傾きを求める手順と、前記複数ラインのライン毎に、前記ライン状に計測された輝度の分布から所定間隔毎に計測データを選択し、相隣接する前記計測データの差分を求め、離散した差分データを得る手順と、この離散した差分データをフーリエ変換し複数の変調伝達関数を求める手順と、この複数の変調伝達関数に基づき解像度を評価する手順とをコンピュータに実行させる。
【0011】
請求項8に記載の解像度評価プログラムは、原稿に描かれた模様のエッジを選択する手順と、前記模様のエッジ近傍における輝度の分布を複数ライン求める手順と、この複数ラインのライン毎に、エッジ座標を求める手順と、前記複数ラインのうち少なくとも1つのラインにおける前記エッジ座標でのエッジの傾きを求めるステップと、前記複数ラインのライン毎に、前記ライン状に計測された輝度の分布から所定間隔毎に計測データを選択し、相隣接する前記計測データの差分を求め、離散した差分データを得る手順と、前記複数ラインのライン毎に、前記離散した差分データに近似曲線をフィッティングさせる手順と、この複数の近似曲線をフーリエ変換し変調伝達関数を求める手順と、この複数の変調伝達関数に基づき解像度を評価する手順とをコンピュータに実行させる。
【0012】
請求項9に記載の解像度評価プログラムは、請求項7または請求項8に記載の解像度評価プログラムにおいて、前記エッジ座標を算出するために、前記被写体のエッジ近傍における低輝度部分の平均値と高輝度部分の平均値とを求め、この低輝度部分の平均値と高輝度部分の平均値とからエッジ閾値を求め、前記計測データを前記ラインの一方から検査して前記輝度がエッジ閾値を跨る位置を検出する処理をコンピュータに実行させる。
【0013】
請求項10に記載の解像度評価プログラムは、請求項7ないし請求項9のいずれか1項に記載の解像度評価プログラムにおいて、前記解像度の評価を、前記複数の変調伝達関数の平均値に基づいて行う処理をコンピュータに実行させる。
請求項11に記載の解像度評価プログラムは、請求項7ないし請求項9のいずれか1項に記載の解像度評価プログラムにおいて、前記解像度の評価を、評価する周波数毎の前記複数の変調伝達関数の最大値に基づいて行う処理をコンピュータに実行させる。
【0014】
請求項12に記載の解像度評価プログラムは、請求項7ないし請求項9のいずれか1項に記載の解像度評価プログラムにおいて、前記解像度の評価を、評価する周波数毎の前記複数の変調伝達関数の最小値に基づいて行う処理をコンピュータに実行させる。
請求項13に記載の光学機器は、請求項7ないし請求項12のいずれか1項に記載の解像度評価プログラムを実行可能であることを特徴とする。
【0015】
【発明の実施の形態】
以下、本発明を図面を用いて詳細に説明する。
図1は、本発明の解像度評価方法の一実施形態の流れを示している。図2は、図1の解像度評価のためのデータ採取方法を示している。
図2に示すように、データ採取をするために、暗幕1の前に模様が描かれた計測用の原稿3が配置されている。また、原稿3の斜め前方には、2つのランプ5が配置されている。さらに、原稿3の正面には、光学機器である電子カメラ7が配置されている。そして、ランプ5の電子カメラ7側には、ランプ5からの光が直接電子カメラ7に行かないように遮光板9が設けられている。
【0016】
本発明の解像度評価方法は、図1に示すように、先ず、ステップS1において、撮影準備をしてから、電子カメラ7により一様な高輝度部分と低輝度部分とからなる原稿3の撮影を行う。ここで、原稿3は、エッジが含まれているものであれば良い。
次いで、ステップS2において、撮影した画像から、エッジを含むように、図3の点線で示すようなエリアを指定する。この指定は、例えば電子カメラ7にノート型電子計算機を接続し、専用ソフトウェア上で行う。
【0017】
次いで、ステップS3において、高輝度部分の計測を行い平均値を算出する。
平均値を算出するとは、グレースケールで例えば8ビット(256階調)分、輝度を分割した場合の輝度の平均値を算出する。
次いで、ステップS4において、低輝度部分の計測を行い、高輝度部分と同様に平均値を算出する。
【0018】
次いで、ステップS5において、低輝度部分の平均値と高輝度部分の平均値とから高輝度部分と低輝度部分との境界であると仮定されるエッジ閾値を求める。
具体的には、高輝度部分の平均値と低輝度部分の平均値との平均をさらに求めることで算出する。
次いで、ステップS6において、図3に点線で示した解像度算出の対象となるエリア内から図4に示すようなエッジ近傍におけるライン状の輝度データ(以下、ラインデータという)を抜き出し、メモリーに記録する。この実施形態では、30本のラインL1〜L30のラインデータを記録することにする。
【0019】
次いで、ステップS7において、ラインデータをラインの一方から所定間隔毎に検索して輝度がエッジ閾値を跨る位置を検出する。そして、輝度がエッジ閾値を跨る前後の位置の座標とその座標での輝度とから線形補間により、エッジ座標を求める。
次いで、ステップS8において、30本のラインL1〜L30の全てについてエッジ座標を求めたか否かが判断される。求め終わっていない場合には、ステップS7に戻る。一方、求め終わった場合には、ステップS9に進む。
【0020】
次いで、ステップS9において、周期の算出を行う。周期とは、ラインL1のエッジ座標をE1とした場合に、エッジ座標が(E1+1)となるまでのライン数を指し、このようにして求めたライン数は傾きに相当する。なお、エッジが全く傾いていなければこのような算出はできないことになるが、実際にエッジが全く傾いていないということは考えられないため、常に周期の算出が行われ、算出結果を得ることができる。なお、以下では、周期として求めたライン数が25である例を示す。
【0021】
次いで、ステップS10において、図5に示すようなラインデータから所定間隔毎にデータを選択し、相隣接する計測データの差分を求め、離散した差分データを算出する。なお、図5のグラフの縦軸は、255段階に階調を分けたときの階調値を示し、横軸は、ラインデータのライン上での位置を示している。
【0022】
次いで、ステップS11において、ステップS10において求めた離散した差分データに対し、最小二乗法(LSM: Least Squares Method)により図6に示すような近似曲線をフィッティングさせる。この図6のグラフの縦軸は、隣接画素との階調差を示し、横軸はラインデータのライン上での位置を示している。
次いで、ステップS12において、ステップS9で算出した傾きにあわせたライン数分の近似曲線を測定範囲として採用する。
【0023】
次いで、ステップS13において、最小二乗法により求めた近似曲線をフーリエ変換する。このフーリエ変換は、離散フーリエ変換(DFT: Discrete Fourier Transform、以下DFTという)により行う。
次いで、ステップS14において、ステップS13で求めた算出結果をメモリーに記録する。
【0024】
次いで、ステップS15において、フーリエ変換したライン数が所定数(ステップS9で求めた周期)の25本になったか否かが判断される。25本に満たなければ、ステップS10に戻り、次のラインで同様の処理を行う。一方、ライン数が25本に達した場合には、ステップS16に進む。
ステップS16では、ステップS9で求めたエッジの傾きに関するデータやステップS13で行ったフーリエ変換結果に基づき、フーリエ変換後の25本のラインの平均値を求める。そして、求めた25本のラインの平均値のライン形状より解像度を推定し、解像度評価を行う。図8に平均値を求める際の例を示す。なお、図8のグラフの縦軸は、DFT後の応答を示し、横軸は、周波数を示している。また、グラフ中、点線で示されるデータは、説明を簡単にするために、25本のラインのうち、ラインL1〜L4を示しており、実線は、これらの平均値を示している。
【0025】
なお、上述したような解像度評価を行うためのプログラムは、記録媒体に記録され、電子カメラ7内に備えられている。
この実施形態の解像度評価方法では、原稿3に描かれた模様のエッジ近傍における輝度の分布を30ライン求め、このライン毎にフーリエ変換を行い、変調伝達関数(MTF: Modulation Transfer Function 以下、MTFという)を求め、最後にライン毎のMTFの平均値を算出するようにしたので、より狭い範囲で、より多数のラインをフーリエ変換すれば正確度や精密度を高くすることができる。ひいては、電子カメラ7の機械的振動による影響を受けにくくすることができる。
【0026】
さらに、この実施形態の解像度評価方法では、原稿3に描かれた模様のエッジ近傍における輝度の分布を30ライン求め、またエッジ座標E1でのエッジの傾きを計算するので、撮影時の電子カメラ7の合焦精度や手ぶれの度合い等を含んだ電子カメラ7側の悪影響を定量的に把握することができる。そのため、撮影条件を変えて同一シーンを撮影した場合に、比較的焦点の合っている画像であるか、或いは比較的手ぶれの起こっていない画像であるかといったことを、解像度を確認することで的確に判断することができる。また、比較的解像度が低いディスプレイでは、少しでも良い解像度で観察することが観察者に望まれるので、この場合には特に効果的である。
【0027】
この実施形態の解像度評価プログラムでは、原稿3に描かれた模様のエッジ近傍における輝度の分布を30ライン求め、このライン毎にフーリエ変換を行い、MTFを求め、最後にライン毎のMTFの平均値を算出するように構成したので、より狭い範囲で、より多数のラインをフーリエ変換させるように、このプログラムを実行させれば、正確度や精密度を高くすることができる。
【0028】
なお、上述した実施形態では、フーリエ変換としてDFTを用いた例について説明したが、高速フーリエ変換(FFT: Fast Fourier Transformation、以下FFTという)を用いても良い。
また、上述した実施形態では、解像度の評価を、複数のMTFの平均値を算出して行った例について説明したが、評価する周波数毎に複数のMTFの中から最大値を有するものを使用するようにしても良い。または、評価する周波数毎に複数のMTFの中から最小値を有するものを使用するようにしても良い。要するに、複数のMTFに対し、相互に連関する数学的処理を施し、的確なMTFを算出することができるのならば如何なる数学的処理を施しても良い。
【0029】
また、上述した実施形態では、光学機器として電子カメラ7を用いて説明したが、解像度評価を行うために原稿3を撮影することができる光学機器なら如何なるものでも良い。例えば、スキャナーや顕微鏡でも良い。
【0030】
【発明の効果】
以上述べたように、本発明の解像度評価方法では、原稿に描かれた模様のエッジ近傍における輝度の分布を複数ライン求め、この複数ラインのライン毎にフーリエ変換し、MTFを求め、最後にライン毎のMTFの平均値を算出するようにしたので、より狭い範囲で、より多数のラインをフーリエ変換すれば正確度や精密度を高くすることができる。ひいては、光学機器の機械的振動による影響を受けにくくすることができる。
【0031】
また、本発明の解像度評価プログラムでは、原稿に描かれた模様のエッジ近傍における輝度の分布を複数ライン求め、この複数のラインのライン毎にフーリエ変換を行い、MTFを求め、最後にライン毎のMTFの平均値を算出するように構成したので、より狭い範囲で、より多数のラインをフーリエ変換させるようにこのプログラムを実行させれば、正確度や精密度を高くすることができる。
【0032】
さらに、本発明の光学機器では、原稿に描かれた模様のエッジ近傍における輝度の分布を複数ライン求め、この複数のラインのライン毎にフーリエ変換を行い、MTFを求め、最後にライン毎のMTFの平均値を算出するプログラムを実行可能であるので、より狭い範囲で、より多数のラインをフーリエ変換すれば、正確度や精密度を高くすることができる。ひいては、光学機器の機械的振動による影響を受けにくくすることができる。
【図面の簡単な説明】
【図1】本発明の解像度評価方法の一実施形態を示す流れ図である。
【図2】解像度評価の対象となる原稿を撮影するときの撮影方法を示す説明図である。
【図3】解像度評価の対象となる部分を示す説明図である。
【図4】解像度評価の対象となるエッジ部分のラインを示す説明図である。
【図5】ラインデータの階調変化を示す説明図である。
【図6】ラインデータの所定間隔毎の階調差に対しフィッティングをかけたとき説明図である。
【図7】フーリエ変換する範囲を示す説明図である。
【図8】フーリエ変換結果を示す説明図である。
【符号の説明】
3 原稿
7 電子カメラ
E1 エッジ座標
L1〜L30 ライン[0001]
TECHNICAL FIELD OF THE INVENTION
The present invention relates to a method for evaluating the resolution of an optical device such as an electronic camera or a scanner.
[0002]
[Prior art]
2. Description of the Related Art Conventionally, as a first resolution evaluation method of an optical device such as an electronic camera or a scanner, there is known a method of photographing a striped document in which a low-brightness portion and a high-brightness portion are mutually arranged (for example, see Patent Reference 1). More specifically, resolution evaluation is performed by photographing a plurality of originals having different stripe widths and analyzing the data. In this resolution, as the width of the stripe becomes narrower, it becomes more difficult to distinguish between a low luminance portion and a high luminance portion, so that the resolution decreases.
[0003]
As a second resolution evaluation method, there is known a method in which a document having a boundary portion (hereinafter, referred to as an edge) between a low-luminance portion and a high-luminance portion is prepared in advance, and this document is photographed (for example, Non-Patent Document 1). Reference 1).
[Patent Document 1]
JP-A-1999-313186 [Non-Patent Document 1]
"Photography-Electronic still picture cameras-resolution measurements", Draft International Standards, USA, International Organization for International Standards.
[0004]
[Problems to be solved by the invention]
However, in the first resolution evaluation method described above, a plurality of originals have to be photographed at once or a plurality of times. When photographing a plurality of documents at once, it is difficult to compare partial resolutions within one image, and when photographing a plurality of documents, the arrangement of the documents must be changed many times, which is troublesome. There was a problem.
[0005]
Such a problem does not occur in the second resolution evaluation method, and the second resolution evaluation method is more general than the first resolution evaluation method. However, the second resolution evaluation method has a problem that the measurement data becomes inaccurate due to mechanical vibration of the optical device.
The present invention has been made to solve such a conventional problem, and has as its object to provide a resolution evaluation method which is hardly affected by mechanical vibration of an optical device.
[0006]
[Means for Solving the Problems]
2. The resolution evaluation method according to
[0007]
3. The resolution evaluation method according to
[0008]
According to a third aspect of the present invention, in the resolution evaluation method according to the first or second aspect, in order to calculate the edge coordinates, an average value and a high luminance of a low luminance portion near an edge of the subject are calculated. Find the average value of the portion, find the edge threshold from the average value of the low luminance portion and the average value of the high luminance portion, inspect the measurement data from one of the lines and determine the position where the luminance crosses the edge threshold. It is characterized by detecting.
[0009]
A resolution evaluation method according to a fourth aspect is the resolution evaluation method according to any one of the first to third aspects, wherein the evaluation of the resolution is performed based on an average value of the plurality of modulation transfer functions. It is characterized by the following.
A resolution evaluation method according to
[0010]
A resolution evaluation method according to
The resolution evaluation program according to claim 7, further comprising: a step of selecting an edge of a pattern drawn on the document; a step of obtaining a plurality of lines of luminance distribution in the vicinity of the edge of the pattern; A step of obtaining coordinates, a step of obtaining an inclination of an edge at the edge coordinates in at least one of the plurality of lines, and a predetermined interval from a distribution of luminance measured in a line shape for each of the plurality of lines. Selecting measurement data for each, obtaining a difference between the adjacent measurement data, obtaining discrete difference data, a procedure for obtaining a plurality of modulation transfer functions by Fourier-transforming the discrete difference data, And a procedure for evaluating the resolution based on the modulation transfer function.
[0011]
9. The resolution evaluation program according to claim 8, further comprising: a step of selecting an edge of a pattern drawn on the document; a step of obtaining a plurality of lines of a luminance distribution in the vicinity of the edge of the pattern; A step of obtaining coordinates, a step of obtaining an inclination of an edge at the edge coordinates in at least one of the plurality of lines, and a predetermined interval from the distribution of the luminance measured in the form of a line for each of the plurality of lines. Selecting measurement data for each, determining the difference between adjacent measurement data, obtaining discrete difference data, and for each of the plurality of lines, fitting an approximate curve to the discrete difference data, Fourier transforming the plurality of approximate curves to obtain a modulation transfer function, and evaluating the resolution based on the plurality of modulation transfer functions To perform the procedure on your computer.
[0012]
The resolution evaluation program according to
[0013]
A resolution evaluation program according to a tenth aspect is the resolution evaluation program according to any one of the seventh to ninth aspects, wherein the evaluation of the resolution is performed based on an average value of the plurality of modulation transfer functions. Causes the computer to execute the processing.
A resolution evaluation program according to
[0014]
According to a twelfth aspect of the present invention, in the resolution evaluation program according to any one of the seventh to ninth aspects, the resolution evaluation is performed by determining a minimum of the plurality of modulation transfer functions for each frequency to be evaluated. Causes the computer to execute processing based on the value.
According to a thirteenth aspect of the present invention, an optical apparatus is capable of executing the resolution evaluation program according to any one of the seventh to twelfth aspects.
[0015]
BEST MODE FOR CARRYING OUT THE INVENTION
Hereinafter, the present invention will be described in detail with reference to the drawings.
FIG. 1 shows the flow of an embodiment of the resolution evaluation method of the present invention. FIG. 2 shows a data collection method for the resolution evaluation of FIG.
As shown in FIG. 2, in order to collect data, a
[0016]
In the resolution evaluation method of the present invention, as shown in FIG. 1, first, in step S1, after preparing for photographing, the electronic camera 7 photographs the original 3 composed of a uniform high-luminance portion and a low-luminance portion. Do. Here, the
Next, in step S2, an area as indicated by a dotted line in FIG. 3 is designated from the captured image so as to include an edge. This designation is performed, for example, by connecting a notebook computer to the electronic camera 7 and using dedicated software.
[0017]
Next, in step S3, a high-luminance portion is measured and an average value is calculated.
Calculating the average value means calculating the average value of the luminance when the luminance is divided by, for example, 8 bits (256 gradations) in the gray scale.
Next, in step S4, a low luminance portion is measured, and an average value is calculated similarly to the high luminance portion.
[0018]
Next, in step S5, an edge threshold value assumed to be a boundary between the high luminance portion and the low luminance portion is obtained from the average value of the low luminance portion and the average value of the high luminance portion.
Specifically, it is calculated by further obtaining the average of the average value of the high luminance portion and the average value of the low luminance portion.
Next, in step S6, line-shaped luminance data (hereinafter, referred to as line data) near the edge as shown in FIG. 4 is extracted from the area for resolution calculation indicated by the dotted line in FIG. 3 and recorded in the memory. . In this embodiment, line data of 30 lines L1 to L30 is recorded.
[0019]
Next, in step S7, line data is searched from one of the lines at predetermined intervals to detect a position where the luminance crosses the edge threshold. Then, edge coordinates are obtained by linear interpolation from the coordinates of the positions before and after the luminance crosses the edge threshold value and the luminance at the coordinates.
Next, in step S8, it is determined whether or not edge coordinates have been obtained for all of the 30 lines L1 to L30. If not, the process returns to step S7. On the other hand, if the calculation has been completed, the process proceeds to step S9.
[0020]
Next, in step S9, a cycle is calculated. The period refers to the number of lines until the edge coordinates become (E1 + 1) when the edge coordinates of the line L1 are E1, and the number of lines obtained in this manner corresponds to the slope. Note that such calculation cannot be performed if the edge is not inclined at all, but since it is not considered that the edge is not actually inclined at all, the period is always calculated and the calculation result may be obtained. it can. Hereinafter, an example in which the number of lines obtained as a cycle is 25 will be described.
[0021]
Next, in step S10, data is selected at predetermined intervals from the line data as shown in FIG. 5, a difference between adjacent measurement data is obtained, and discrete difference data is calculated. Note that the vertical axis of the graph in FIG. 5 indicates the gray scale value when the gray scale is divided into 255 levels, and the horizontal axis indicates the position of the line data on the line.
[0022]
Next, in step S11, an approximation curve as shown in FIG. 6 is fitted to the discrete difference data obtained in step S10 by a least squares method (LSM: Last Squares Method). The vertical axis of the graph of FIG. 6 indicates the gradation difference between adjacent pixels, and the horizontal axis indicates the position of the line data on the line.
Next, in step S12, an approximate curve for the number of lines corresponding to the inclination calculated in step S9 is adopted as a measurement range.
[0023]
Next, in step S13, the approximate curve obtained by the least square method is subjected to Fourier transform. The Fourier transform is performed by a discrete Fourier transform (DFT).
Next, in step S14, the calculation result obtained in step S13 is recorded in the memory.
[0024]
Next, in step S15, it is determined whether or not the number of Fourier-transformed lines has reached a predetermined number (the cycle obtained in step S9) of 25 lines. If the number is less than 25, the process returns to step S10, and the same processing is performed on the next line. On the other hand, if the number of lines has reached 25, the process proceeds to step S16.
In step S16, an average value of the 25 lines after the Fourier transform is obtained based on the data on the inclination of the edge obtained in step S9 and the result of the Fourier transform performed in step S13. Then, the resolution is estimated from the obtained line shape of the average value of the 25 lines, and the resolution is evaluated. FIG. 8 shows an example of calculating the average value. The vertical axis of the graph in FIG. 8 indicates the response after DFT, and the horizontal axis indicates the frequency. Further, in the graph, data indicated by a dotted line indicates lines L1 to L4 among the 25 lines for simplicity of explanation, and a solid line indicates an average value thereof.
[0025]
Note that a program for performing the above-described resolution evaluation is recorded on a recording medium and provided in the electronic camera 7.
In the resolution evaluation method according to this embodiment, 30 lines of luminance distribution near the edge of the pattern drawn on the
[0026]
Further, in the resolution evaluation method of this embodiment, 30 lines of luminance distribution near the edge of the pattern drawn on the
[0027]
In the resolution evaluation program of this embodiment, 30 lines of luminance distribution near the edge of the pattern drawn on the
[0028]
In the above-described embodiment, an example in which DFT is used as Fourier transform has been described, but fast Fourier transform (FFT: Fast Fourier Transformation, hereinafter referred to as FFT) may be used.
Also, in the above-described embodiment, an example has been described in which the evaluation of the resolution is performed by calculating the average value of a plurality of MTFs. However, the one having the maximum value among the plurality of MTFs is used for each frequency to be evaluated. You may do it. Alternatively, one having a minimum value among a plurality of MTFs for each frequency to be evaluated may be used. In short, any mathematical processing may be performed on a plurality of MTFs as long as the mathematical processing associated with each other can be performed and an accurate MTF can be calculated.
[0029]
Further, in the above-described embodiment, the electronic camera 7 has been described as the optical device. However, any optical device that can photograph the
[0030]
【The invention's effect】
As described above, in the resolution evaluation method of the present invention, a plurality of luminance distributions in the vicinity of the edge of a pattern drawn on a document are obtained, a Fourier transform is performed for each of the plurality of lines, an MTF is obtained, and finally a line is obtained. Since the average value of the MTF for each is calculated, the accuracy and precision can be increased by performing a Fourier transform on a larger number of lines in a narrower range. As a result, the optical device can be made less susceptible to mechanical vibration.
[0031]
In the resolution evaluation program according to the present invention, a plurality of luminance distributions in the vicinity of edges of a pattern drawn on a document are obtained, a Fourier transform is performed for each of the plurality of lines, an MTF is obtained, and finally, an MTF is obtained for each line. Since the average value of the MTF is configured to be calculated, the accuracy and precision can be increased by executing this program so as to perform a Fourier transform on a larger number of lines in a narrower range.
[0032]
Further, in the optical apparatus of the present invention, a plurality of luminance distributions in the vicinity of the edge of the pattern drawn on the document are obtained, a Fourier transform is performed for each of the plurality of lines, an MTF is obtained, and finally, an MTF for each line is obtained. Since the program for calculating the average value of can be executed, the accuracy and precision can be increased by performing a Fourier transform on a larger number of lines in a narrower range. As a result, the optical device can be made less susceptible to mechanical vibration.
[Brief description of the drawings]
FIG. 1 is a flowchart showing one embodiment of a resolution evaluation method of the present invention.
FIG. 2 is an explanatory diagram showing a photographing method when photographing a document to be evaluated for resolution.
FIG. 3 is an explanatory diagram showing a portion to be evaluated for resolution.
FIG. 4 is an explanatory diagram showing a line of an edge portion to be evaluated for resolution.
FIG. 5 is an explanatory diagram showing a gradation change of line data.
FIG. 6 is an explanatory diagram when fitting is performed on a gradation difference at predetermined intervals of line data.
FIG. 7 is an explanatory diagram showing a range for performing Fourier transform.
FIG. 8 is an explanatory diagram showing a Fourier transform result.
[Explanation of symbols]
3 Original 7 Electronic camera E1 Edge coordinates L1 to L30 Line
Claims (13)
前記模様のエッジ近傍における輝度の分布を複数ライン求めるステップと、
この複数ラインのライン毎に、エッジ座標を求めるステップと、
前記複数ラインのうち少なくとも1つのラインにおける前記エッジ座標でのエッジの傾きを求めるステップと、
前記複数ラインのライン毎に、前記ライン状に計測された輝度の分布から所定間隔毎に計測データを選択し、相隣接する前記計測データの差分を求め、離散した差分データを得るステップと、
この離散した差分データをフーリエ変換し複数の変調伝達関数を求めるステップと、
この複数の変調伝達関数に基づき解像度を評価するステップと、
を有することを特徴とする解像度評価方法。Selecting an edge of a pattern drawn on the manuscript;
Obtaining a plurality of lines of luminance distribution near the edge of the pattern;
Obtaining edge coordinates for each of the plurality of lines;
Determining the slope of the edge at the edge coordinates in at least one of the plurality of lines;
For each of the plurality of lines, select measurement data at predetermined intervals from the distribution of the luminance measured in the line shape, determine the difference between adjacent measurement data, obtain discrete difference data,
Fourier transforming the discrete difference data to obtain a plurality of modulation transfer functions;
Evaluating a resolution based on the plurality of modulation transfer functions;
A resolution evaluation method comprising:
前記模様のエッジ近傍における輝度の分布を複数ライン求めるステップと、
この複数ラインのライン毎に、エッジ座標を求めるステップと、
前記複数ラインのうち少なくとも1つのラインにおける前記エッジ座標でのエッジの傾きを求めるステップと、
前記複数ラインのライン毎に、前記ライン状に計測された輝度の分布から所定間隔毎に計測データを選択し、相隣接する前記計測データの差分を求め、離散した差分データを得るステップと、
前記複数ラインのライン毎に、前記離散した差分データに近似曲線をフィッティングさせるステップと、
この複数の近似曲線をフーリエ変換し変調伝達関数を求めるステップと、
この複数の変調伝達関数に基づき解像度を評価するステップと、
を有することを特徴とする解像度評価方法。Selecting an edge of a pattern drawn on the manuscript;
Obtaining a plurality of lines of luminance distribution near the edge of the pattern;
Obtaining edge coordinates for each of the plurality of lines;
Determining the slope of the edge at the edge coordinates in at least one of the plurality of lines;
For each of the plurality of lines, select measurement data at predetermined intervals from the distribution of the luminance measured in the line shape, determine the difference between adjacent measurement data, obtain discrete difference data,
Fitting an approximate curve to the discrete difference data for each of the plurality of lines;
Fourier transforming the plurality of approximate curves to obtain a modulation transfer function;
Evaluating a resolution based on the plurality of modulation transfer functions;
A resolution evaluation method comprising:
前記解像度の評価を、前記複数の変調伝達関数の平均値に基づいて行うことを特徴とする解像度評価方法。The resolution evaluation method according to any one of claims 1 to 3,
A resolution evaluation method, wherein the evaluation of the resolution is performed based on an average value of the plurality of modulation transfer functions.
前記解像度の評価を、評価する周波数毎の前記複数の変調伝達関数の最大値に基づいて行うことを特徴とする解像度評価方法。The resolution evaluation method according to any one of claims 1 to 3,
A resolution evaluation method, wherein the evaluation of the resolution is performed based on a maximum value of the plurality of modulation transfer functions for each frequency to be evaluated.
前記解像度の評価を、評価する周波数毎の前記複数の変調伝達関数の最小値に基づいて行うことを特徴とする解像度評価方法。The resolution evaluation method according to any one of claims 1 to 3,
A resolution evaluation method, wherein the evaluation of the resolution is performed based on a minimum value of the plurality of modulation transfer functions for each frequency to be evaluated.
前記模様のエッジ近傍における輝度の分布を複数ライン求める手順と、
この複数ラインのライン毎に、エッジ座標を求める手順と、
前記複数ラインのうち少なくとも1つのラインにおける前記エッジ座標でのエッジの傾きを求める手順と、
前記複数ラインのライン毎に、前記ライン状に計測された輝度の分布から所定間隔毎に計測データを選択し、相隣接する前記計測データの差分を求め、離散した差分データを得る手順と、
この離散した差分データをフーリエ変換し複数の変調伝達関数を求める手順と、
この複数の変調伝達関数に基づき解像度を評価する手順と、
をコンピュータに実行させるための解像度評価プログラム。Steps to select the edges of the pattern drawn on the manuscript,
A procedure for obtaining a plurality of lines of the distribution of luminance near the edge of the pattern;
A procedure for obtaining edge coordinates for each of the plurality of lines;
A step of obtaining an inclination of an edge at the edge coordinates in at least one of the plurality of lines;
For each line of the plurality of lines, select measurement data at predetermined intervals from the distribution of the luminance measured in the line shape, determine the difference between adjacent measurement data, a procedure to obtain discrete difference data,
Fourier transforming the discrete difference data to obtain a plurality of modulation transfer functions;
A step of evaluating the resolution based on the plurality of modulation transfer functions;
A resolution evaluation program for causing a computer to execute.
前記模様のエッジ近傍における輝度の分布を複数ライン求める手順と、
この複数ラインのライン毎に、エッジ座標を求める手順と、
前記複数ラインのうち少なくとも1つのラインにおける前記エッジ座標でのエッジの傾きを求めるステップと、
前記複数ラインのライン毎に、前記ライン状に計測された輝度の分布から所定間隔毎に計測データを選択し、相隣接する前記計測データの差分を求め、離散した差分データを得る手順と、
前記複数ラインのライン毎に、前記離散した差分データに近似曲線をフィッティングさせる手順と、
この複数の近似曲線をフーリエ変換し変調伝達関数を求める手順と、
この複数の変調伝達関数に基づき解像度を評価する手順と、
をコンピュータに実行させるための解像度評価プログラム。Steps to select the edges of the pattern drawn on the manuscript,
A procedure for obtaining a plurality of lines of the distribution of luminance near the edge of the pattern;
A procedure for obtaining edge coordinates for each of the plurality of lines;
Determining the slope of the edge at the edge coordinates in at least one of the plurality of lines;
For each line of the plurality of lines, select measurement data at predetermined intervals from the distribution of the luminance measured in the line shape, determine the difference between adjacent measurement data, a procedure to obtain discrete difference data,
A procedure for fitting an approximate curve to the discrete difference data for each of the plurality of lines,
Fourier transforming the plurality of approximate curves to obtain a modulation transfer function;
A step of evaluating the resolution based on the plurality of modulation transfer functions;
A resolution evaluation program for causing a computer to execute.
前記エッジ座標を算出するために、前記被写体のエッジ近傍における低輝度部分の平均値と高輝度部分の平均値とを求め、この低輝度部分の平均値と高輝度部分の平均値とからエッジ閾値を求め、前記計測データを前記ラインの一方から検査して前記輝度がエッジ閾値を跨る位置を検出する処理をコンピュータに実行させるための解像度評価プログラム。In the resolution evaluation program according to claim 7 or 8,
In order to calculate the edge coordinates, an average value of a low-luminance portion and an average value of a high-luminance portion in the vicinity of the edge of the subject are obtained, and an edge threshold is calculated from the average value of the low-luminance portion and the average value of the high-luminance portion. A resolution evaluation program for causing a computer to execute a process of detecting the measured data from one of the lines and detecting a position where the luminance crosses an edge threshold.
前記解像度の評価を、前記複数の変調伝達関数の平均値に基づいて行う処理をコンピュータに実行させるための解像度評価プログラム。In the resolution evaluation program according to any one of claims 7 to 9,
A resolution evaluation program for causing a computer to execute a process of evaluating the resolution based on an average value of the plurality of modulation transfer functions.
前記解像度の評価を、評価する周波数毎の前記複数の変調伝達関数の最大値に基づいて行う処理をコンピュータに実行させるための解像度評価プログラム。In the resolution evaluation program according to any one of claims 7 to 9,
A resolution evaluation program for causing a computer to execute processing for evaluating the resolution based on the maximum value of the plurality of modulation transfer functions for each frequency to be evaluated.
前記解像度の評価を、評価する周波数毎の前記複数の変調伝達関数の最小値に基づいて行う処理をコンピュータに実行させるための解像度評価プログラム。In the resolution evaluation program according to any one of claims 7 to 9,
A resolution evaluation program for causing a computer to execute a process of evaluating the resolution based on a minimum value of the plurality of modulation transfer functions for each frequency to be evaluated.
ことを特徴とする光学機器。An optical apparatus capable of executing the resolution evaluation program according to claim 7.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2002294967A JP2004134861A (en) | 2002-10-08 | 2002-10-08 | Resolution evaluation method, resolution evaluation program, and optical apparatus |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2002294967A JP2004134861A (en) | 2002-10-08 | 2002-10-08 | Resolution evaluation method, resolution evaluation program, and optical apparatus |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2004134861A true JP2004134861A (en) | 2004-04-30 |
Family
ID=32285362
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2002294967A Pending JP2004134861A (en) | 2002-10-08 | 2002-10-08 | Resolution evaluation method, resolution evaluation program, and optical apparatus |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP2004134861A (en) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2009038431A (en) * | 2007-07-31 | 2009-02-19 | Panasonic Corp | Moving image resolution evaluation device of display device |
JP2011238292A (en) * | 2004-06-16 | 2011-11-24 | Microsoft Corp | Method and system for reducing effects of undesired signals in an infrared imaging system |
JP2019194577A (en) * | 2018-04-27 | 2019-11-07 | 日本放送協会 | Mtf measurement device and program |
CN111045287A (en) * | 2017-08-08 | 2020-04-21 | 广东弘景光电科技股份有限公司 | Detection system of wide-angle camera module |
JP2020071144A (en) * | 2018-10-31 | 2020-05-07 | 日本放送協会 | Mtf measuring device and program therefor |
JP2020071145A (en) * | 2018-10-31 | 2020-05-07 | 日本放送協会 | Mtf measuring device and program therefor |
CN116847067A (en) * | 2023-08-30 | 2023-10-03 | 荣耀终端有限公司 | Imaging system evaluation method and device |
-
2002
- 2002-10-08 JP JP2002294967A patent/JP2004134861A/en active Pending
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2011238292A (en) * | 2004-06-16 | 2011-11-24 | Microsoft Corp | Method and system for reducing effects of undesired signals in an infrared imaging system |
JP2009038431A (en) * | 2007-07-31 | 2009-02-19 | Panasonic Corp | Moving image resolution evaluation device of display device |
CN111045287A (en) * | 2017-08-08 | 2020-04-21 | 广东弘景光电科技股份有限公司 | Detection system of wide-angle camera module |
CN111045287B (en) * | 2017-08-08 | 2021-08-24 | 广东弘景光电科技股份有限公司 | Detection system of wide-angle camera module |
JP2019194577A (en) * | 2018-04-27 | 2019-11-07 | 日本放送協会 | Mtf measurement device and program |
JP7257231B2 (en) | 2018-04-27 | 2023-04-13 | 日本放送協会 | MTF measuring device and its program |
JP2020071144A (en) * | 2018-10-31 | 2020-05-07 | 日本放送協会 | Mtf measuring device and program therefor |
JP2020071145A (en) * | 2018-10-31 | 2020-05-07 | 日本放送協会 | Mtf measuring device and program therefor |
JP7190328B2 (en) | 2018-10-31 | 2022-12-15 | 日本放送協会 | MTF measuring device and its program |
JP7190327B2 (en) | 2018-10-31 | 2022-12-15 | 日本放送協会 | MTF measuring device and its program |
CN116847067A (en) * | 2023-08-30 | 2023-10-03 | 荣耀终端有限公司 | Imaging system evaluation method and device |
CN116847067B (en) * | 2023-08-30 | 2023-11-24 | 荣耀终端有限公司 | Imaging system evaluation method and device |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Abdelhamed et al. | A high-quality denoising dataset for smartphone cameras | |
JP2008500529A (en) | Method for characterizing a digital imaging system | |
JP4416795B2 (en) | Correction method | |
KR100961027B1 (en) | Method for detecting streaks in digital images | |
CN112033965B (en) | 3D arc surface defect detection method based on differential image analysis | |
CN111083365B (en) | Method and device for rapidly detecting optimal focal plane position | |
WO2009029365A1 (en) | Systems and methods for determination of a camera imperfection for an image | |
JP5659623B2 (en) | Exposure attribute setting method and computer-readable storage medium | |
CN110915193B (en) | Image processing system, server device, image processing method, and recording medium | |
US20200389573A1 (en) | Image processing system, image processing method and storage medium | |
TW201702600A (en) | Measurement device and measurement method for rotation of round body and non-transitory information readable medium | |
JP5705711B2 (en) | Crack detection method | |
JP2007309764A (en) | Mtf measuring instrument, method, and program | |
KR101072945B1 (en) | Apparatus for balancing flashed image color using depth estimation and method therefor | |
TW201326776A (en) | Lens test device and method | |
JP2004134861A (en) | Resolution evaluation method, resolution evaluation program, and optical apparatus | |
JP2008281481A (en) | Apparatus and method for measuring resolution | |
WO2018094823A1 (en) | Uniformity test method and system for display screen | |
Branca et al. | Texture MTF from images of natural scenes | |
Williams et al. | Low-frequency MTF estimation for digital imaging devices using slanted-edge analysis | |
CN112734721B (en) | Optical axis deflection angle detection method, device, equipment and medium | |
JP6897448B2 (en) | Line width estimation programs, devices, and methods | |
Boiangiu et al. | Target Validation and Image Color Calibration | |
CN113034553A (en) | Evaluation method of image registration algorithm, terminal and storage medium | |
van Zwanenberg et al. | Camera system performance derived from natural scenes |