JP2004134861A - Resolution evaluation method, resolution evaluation program, and optical apparatus - Google Patents

Resolution evaluation method, resolution evaluation program, and optical apparatus Download PDF

Info

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
Application number
JP2002294967A
Other languages
Japanese (ja)
Inventor
Hiroshi Sano
佐野 央
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nikon Corp
Original Assignee
Nikon Corp
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Nikon Corp filed Critical Nikon Corp
Priority to JP2002294967A priority Critical patent/JP2004134861A/en
Publication of JP2004134861A publication Critical patent/JP2004134861A/en
Pending legal-status Critical Current

Links

Images

Landscapes

  • Facsimiles In General (AREA)
  • Testing, Inspecting, Measuring Of Stereoscopic Televisions And Televisions (AREA)

Abstract

<P>PROBLEM TO BE SOLVED: To provide a method for resolution evaluation of an optical apparatus such as an electronic camera and a scanner, which is hardly vulnerable to the mechanical vibration of the optical apparatus. <P>SOLUTION: The resolution evaluation method includes the steps of: selecting the edges of a pattern drawn on an original; obtaining the distribution of luminance values near the edges of the pattern by a plurality of lines; obtaining the edge coordinates of each of a plurality of the lines; obtaining the inclination of the edges at the edge coordinates of at least one line among the plurality of lines; selecting measured data from the distribution of the linearly measured luminance values at each prescribed interval, obtaining the difference between the measured data adjacent to each other and obtaining discrete difference data; applying Fourier transform to the discrete difference data to obtain a plurality of modulation transfer functions; and evaluating the resolution on the basis of the plurality of modulation transfer functions. <P>COPYRIGHT: (C)2004,JPO

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 claim 1, wherein the step of selecting an edge of the pattern drawn on the document, the step of obtaining a plurality of luminance distributions in the vicinity of the edge of the pattern, and the step of: Obtaining coordinates, obtaining an inclination of an edge at the edge coordinates in at least one of the plurality of lines, and determining, for each of the plurality of lines, a predetermined interval from a distribution of luminance measured in a line shape. Selecting measurement data for each, obtaining a difference between adjacent measurement data, obtaining discrete difference data, obtaining a plurality of modulation transfer functions by Fourier-transforming the discrete difference data; and Evaluating the resolution based on the modulation transfer function.
[0007]
3. The resolution evaluation method according to claim 2, wherein the step of selecting an edge of the pattern drawn on the document, the step of obtaining a plurality of lines of the luminance distribution near the edge of the pattern, Obtaining coordinates, obtaining an inclination of an edge at the edge coordinates in at least one of the plurality of lines, and determining, for each of the plurality of lines, a predetermined interval from a distribution of luminance measured in a line shape. Selecting measurement data for each, determining the difference between adjacent measurement data, obtaining discrete difference data, 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 Characterized by a step of evaluating can resolution.
[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 claim 5, wherein in the resolution evaluation method according to any one of claims 1 to 3, the evaluation of the resolution is performed by calculating the maximum of the plurality of modulation transfer functions for each frequency to be evaluated. It is characterized in that it is performed based on the value.
[0010]
A resolution evaluation method according to claim 6, wherein in the resolution evaluation method according to any one of claims 1 to 3, the evaluation of the resolution is performed by minimizing the plurality of modulation transfer functions for each frequency to be evaluated. It is characterized in that it is performed based on the value.
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 claim 9 is the resolution evaluation program according to claim 7 or 8, wherein the average value and the high luminance of a low luminance portion near the edge of the subject are calculated to calculate the edge coordinates. 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. Causes the computer to execute the detection processing.
[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 claim 11, wherein in the resolution evaluation program according to any one of claims 7 to 9, the resolution evaluation is performed by determining a maximum of the plurality of modulation transfer functions for each frequency to be evaluated. Causes the computer to execute processing based on the value.
[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 measurement document 3 on which a pattern is drawn is arranged in front of a blackout curtain 1. Two lamps 5 are arranged diagonally forward of the document 3. Further, an electronic camera 7 as an optical device is arranged in front of the document 3. A light-shielding plate 9 is provided on the electronic camera 7 side of the lamp 5 so that light from the lamp 5 does not directly go to the electronic camera 7.
[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 document 3 only needs to include an edge.
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 document 3 are obtained, Fourier transform is performed for each line, and a modulation transfer function (MTF: MTF). ), And finally calculate the average value of the MTF for each line. Therefore, 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 electronic camera 7 can be less affected by mechanical vibration.
[0026]
Further, in the resolution evaluation method of this embodiment, 30 lines of luminance distribution near the edge of the pattern drawn on the document 3 are obtained, and the inclination of the edge at the edge coordinate E1 is calculated. The adverse effects on the electronic camera 7 side, including the focusing accuracy and the degree of camera shake, can be quantitatively grasped. Therefore, when the same scene is shot under different shooting conditions, whether the image is relatively focused or relatively free from camera shake can be determined by checking the resolution. Can be determined. In addition, in a display having a relatively low resolution, it is particularly effective in this case because it is desired for the observer to observe the image with a little better resolution.
[0027]
In the resolution evaluation program of this embodiment, 30 lines of luminance distribution near the edge of the pattern drawn on the document 3 are obtained, Fourier transform is performed for each line, an MTF is obtained, and finally the average value of the MTF for each line is obtained. Is 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.
[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 document 3 for evaluating the resolution may be used. For example, a scanner or a microscope may be used.
[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:
請求項1または請求項2に記載の解像度評価方法において、前記エッジ座標を算出するために、前記被写体のエッジ近傍における低輝度部分の平均値と高輝度部分の平均値とを求め、この低輝度部分の平均値と高輝度部分の平均値とからエッジ閾値を求め、前記計測データを前記ラインの一方から検査して前記輝度がエッジ閾値を跨る位置を検出することを特徴とする解像度評価方法。3. The resolution evaluation method according to claim 1, wherein an average value of a low-luminance portion and an average value of a high-luminance portion in the vicinity of an edge of the subject are calculated to calculate the edge coordinates. A resolution evaluation method, wherein an edge threshold value is obtained from an average value of a portion and an average value of a high-brightness portion, and the measurement data is inspected from one of the lines to detect a position where the luminance crosses the edge threshold value. 請求項1ないし請求項3のいずれか1項に記載の解像度評価方法において、
前記解像度の評価を、前記複数の変調伝達関数の平均値に基づいて行うことを特徴とする解像度評価方法。
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.
請求項1ないし請求項3のいずれか1項に記載の解像度評価方法において、
前記解像度の評価を、評価する周波数毎の前記複数の変調伝達関数の最大値に基づいて行うことを特徴とする解像度評価方法。
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.
請求項1ないし請求項3のいずれか1項に記載の解像度評価方法において、
前記解像度の評価を、評価する周波数毎の前記複数の変調伝達関数の最小値に基づいて行うことを特徴とする解像度評価方法。
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.
請求項7または請求項8に記載の解像度評価プログラムにおいて、
前記エッジ座標を算出するために、前記被写体のエッジ近傍における低輝度部分の平均値と高輝度部分の平均値とを求め、この低輝度部分の平均値と高輝度部分の平均値とからエッジ閾値を求め、前記計測データを前記ラインの一方から検査して前記輝度がエッジ閾値を跨る位置を検出する処理をコンピュータに実行させるための解像度評価プログラム。
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.
請求項7ないし請求項9のいずれか1項に記載の解像度評価プログラムにおいて、
前記解像度の評価を、前記複数の変調伝達関数の平均値に基づいて行う処理をコンピュータに実行させるための解像度評価プログラム。
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.
請求項7ないし請求項9のいずれか1項に記載の解像度評価プログラムにおいて、
前記解像度の評価を、評価する周波数毎の前記複数の変調伝達関数の最大値に基づいて行う処理をコンピュータに実行させるための解像度評価プログラム。
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.
請求項7ないし請求項9のいずれか1項に記載の解像度評価プログラムにおいて、
前記解像度の評価を、評価する周波数毎の前記複数の変調伝達関数の最小値に基づいて行う処理をコンピュータに実行させるための解像度評価プログラム。
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.
請求項7ないし請求項12のいずれか1項に記載の解像度評価プログラムを実行可能である
ことを特徴とする光学機器。
An optical apparatus capable of executing the resolution evaluation program according to claim 7.
JP2002294967A 2002-10-08 2002-10-08 Resolution evaluation method, resolution evaluation program, and optical apparatus Pending JP2004134861A (en)

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)

* Cited by examiner, † Cited by third party
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

Cited By (12)

* Cited by examiner, † Cited by third party
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