特許 6061133 テストチャート及びその使用方法

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B2)

(11)【特許番号】6061133

(24)【登録日】2016年12月22日

(45)【発行日】2017年1月18日

(54)【発明の名称】テストチャート及びその使用方法

(51)【国際特許分類】

   G01M 11/02 20060101AFI20170106BHJP

   H04N 5/225 20060101ALI20170106BHJP

【ＦＩ】

   G01M11/02 A

   H04N5/225 F

【請求項の数】18

【全頁数】18

(21)【出願番号】特願2012-244379(P2012-244379)

(22)【出願日】2012年11月6日

(65)【公開番号】特開2014-92508(P2014-92508A)

(43)【公開日】2014年5月19日

【審査請求日】2015年10月20日

(73)【特許権者】

【識別番号】000115603

【氏名又は名称】リーダー電子株式会社

(74)【代理人】

【識別番号】100140109

【弁理士】

【氏名又は名称】小野 新次郎

(74)【代理人】

【識別番号】100075270

【弁理士】

【氏名又は名称】小林 泰

(74)【代理人】

【識別番号】100096013

【弁理士】

【氏名又は名称】富田 博行

(74)【代理人】

【識別番号】100092967

【弁理士】

【氏名又は名称】星野 修

(74)【代理人】

【識別番号】100120112

【弁理士】

【氏名又は名称】中西 基晴

(72)【発明者】

【氏名】王 暁冬

【審査官】 横尾 雅一

(56)【参考文献】

【文献】 特開昭６２−２３７８８９（ＪＰ，Ａ）

【文献】 特開２００６−２４５８９１（ＪＰ，Ａ）

【文献】 特開２０１１−１０１１００（ＪＰ，Ａ）

【文献】 吉田興夫、外２名，「ＩＴＥテストチャートを用いたテレビジョンシステムの評価（ＩＩ）」，テレビジョン学会誌，１９８４年，第３８巻第３号，第２４９頁−第２５７頁

(58)【調査した分野】（Int.Cl.，ＤＢ名）

Ｇ０１Ｍ １１／００−１１／０８

Ｈ０４Ｎ ５／２２２−５／２５７

Ｈ０４Ｎ １３／００−１７／０６

Ｇ０３Ｂ ３／００−４３／０２

ＪＳＴＰｌｕｓ／ＪＭＥＤＰｌｕｓ／ＪＳＴ７５８０（ＪＤｒｅａｍＩＩＩ）

(57)【特許請求の範囲】

【請求項1】

少なくとも１つのチャートを含むテストチャートであって、
前記チャートが、
同心円状に配列された複数のリングであって、該複数のリングが、周辺から中心に向かう方向において増大する空間周波数を有する、複数のリングと、
該複数のリングの外に配置された少なくとも１つの１対のマークであって、該１対のマークが、前記複数のリングの中心を通る横方向または縦方向において配置された、少なくとも１つの１対のマークと、
を含み、
前記少なくとも１つの１対のマークの各対は、空間解像度の測定に使用する測定エリアを定め、該測定エリアの幅は、対応する対のマークの幅により定められる、テストチャート。

【請求項2】

請求項１記載のテストチャートにおいて、前記チャートは、前記複数のリングの中心に前記黒の円形領域を含み、
該円形領域は、前記マークが配置された前記横方向または前記縦方向に直交する方向における前記マークの幅に関連する直径を有する、
テストチャート。

【請求項3】

請求項１記載のテストチャートにおいて、前記複数のリングの各々は、前記空間周波数を定める波形の一方のピークで第１の色となり他方のピークで第２の色となり、前記一方のピーク前記他方のピークとの間において前記第１の色と前記第２の色との間で色が変化する、テストチャート。

【請求項4】

請求項１記載のテストチャートにおいて、
前記テストチャートは、矩形の形状を有し、かつ複数の前記チャートを含み、
該複数の前記チャートは、
前記テストチャートの前記矩形の中央に配置された中心チャートと、
前記矩形の４隅の各々に配置された４つの周辺チャートと、
を含む、テストチャート。

【請求項5】

ａ）被測定対象の撮像装置からの、請求項１〜４のいずれかに記載の前記テストチャートのデジタル撮像信号を受けるステップと、
ｂ）前記デジタル撮像信号から、前記テストチャートの前記少なくとも１つのチャートのうちの１つのチャートの前記１対のマークを検出することにより、前記１つのチャートのある角度における測定エリアを定めるステップと、
ｃ）前記測定エリア内において低い空間周波数から高い空間周波数に向かって複数の測定ラインを定めるステップと、
ｄ）各前記測定ラインの各々において、前記デジタル撮像信号から、前記測定ラインに沿った一連のデジタル・データを取得し、該一連のデジタル・データが定める波形の周期および振幅の少なくとも１つに基づいて各前記測定ラインにおける空間解像度を求めるステップと、
ｅ）複数の前記測定ラインにおいてそれぞれ求めた複数の空間解像度を使用して、前記被測定対象の撮像装置の、前記測定エリアにおける空間解像度を求めるステップと、
を含む、測定方法。

【請求項6】

請求項５記載の測定方法において、
各前記測定ラインにおける空間解像度を求めるステップは、
前記波形のピークの検出を、前記低い空間周波数から開始するステップと、
検出した前記ピークのうちの、前記波形の平均値を超えない一方の極性のピークを検出するステップと、
前記平均値を超えない一方の極性のピークの直前の、前記平均値を超える一方の極性のピークを有する波形位置に対応する前記空間周波数に基づいて、各前記測定ラインにおける空間解像度を定めるステップと、
を含む、測定方法。

【請求項7】

請求項６記載の測定方法において、
前記測定エリアの空間解像度は、前記複数の測定ラインにおいて求めた前記複数の空間解像度の平均値によって定める、測定方法。

【請求項8】

請求項５記載の測定方法において、
各前記測定ラインにおける空間解像度を求めるステップは、
前記波形の周期の検出を、前記低い空間周波数から開始するステップと、
検出した前記周期が、直前に検出した前記周期よりも長いかどうか判定するステップと、
検出した前記周期が直前に検出した周期よりも長いとき、前記直前に検出した周期を有する波形位置に対応する前記空間周波数に基づいて、各前記測定ラインにおける空間解像度を定めるステップと、
を含む、測定方法。

【請求項9】

請求項５記載の測定方法において、
各前記測定ラインにおける空間解像度を求めるステップは、
前記波形の振幅の検出を、前記低い空間周波数から開始するステップと、
検出した前記振幅が、しきい値よりも小さいかどうか判定するステップと、
検出した前記振幅が前記しきい値より小さいとき、該検出した振幅の直前に検出した振幅を有する波形位置に対応する前記空間周波数に基づいて、各前記測定ラインにおける空間解像度を定めるステップと、
を含む、測定方法。

【請求項10】

請求項５記載の測定方法において、
前記測定エリアを、前記１つのチャートの中心を中心としてある角度回転させて、回転後測定エリアを生成するステップと、
前記回転後測定エリアに関して、前記ステップｃ、ｄ、ｅを実行するステップと、
をさらに含む、測定方法。

【請求項11】

請求項５記載の測定方法において、
前記測定エリアにおいて求めた前記空間解像度を表す値を、前記測定エリアの角度と関連付けて表示するステップ、
をさらに含む、測定方法。

【請求項12】

請求項１１記載の測定方法において、
前記空間解像度を表す値と前記測定エリアの角度とを極座標で表示する、測定方法。

【請求項13】

請求項１２記載の測定方法において、
前記極座標における極からの距離は、前記複数のリングの中心と、該距離が表す空間周波数を有する前記チャートの前記複数のリング状の位置との間の距離に関係する、測定方法。

【請求項14】

請求項１３記載の測定方法において、
前記極座標における極からの距離が、前記空間解像度に関係する、測定方法。

【請求項15】

請求項１４記載の測定方法において、
前記極座標における中心からの距離が、前記空間解像度を表す値を表す、測定方法。

【請求項16】

請求項１１記載の測定方法において、
前記空間解像度を表す値と前記測定エリアの角度とを直交座標で表示する、測定方法。

【請求項17】

請求項５から１６のいずれかに記載の方法をコンピュータに実行させるためのコンピュータ・プログラム。

【請求項18】

請求項５から１６のいずれかに記載の方法をコンピュータに実行させるためのコンピュータ読み取り可能命令を格納したコンピュータ読み取り可能記憶媒体。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、カメラなどの撮像装置のためのテストチャート、並びにこれの使用方法に関する。

【背景技術】

【0002】

従来、カメラの空間周波数を評価するために、サーキュラゾーンプレート（CZP： Circular Zone Plate)を使用することが知られている。この公知のCZPチャートを例えばビデオカメラで撮影してそのビデオカメラの映像信号を波形モニターなどの測定器で評価するとき、測定器の観測波形の変化が、ビデオカメラのレンズの特性によるものなのかあるいはビデオカメラ本体の特性によるものなのか区別することができない、という問題がある。この問題に対処するため、特開２０１１−１０１１００号公報では、この公知のCZPチャートを用いて作成したテストチャートを開示している。この開示されたテストチャートとして、カメラ等の水平解像度を評価するものと垂直解像度を評価する２つのチャートがある。水平解像度評価用のテストチャート（上記公報の図２参照）は、公知のCZPチャートを中心部から左右に二等分する垂直線で切断し、切断された左側の左辺と他方の右辺とを接合して作成され、これによって、空間周波数特性が中心部において最も高く、水平方向周辺部へ向かうにつれて低くなるようにしている。一方、垂直解像度評価用のテストチャート（上記公報の図３参照）は、公知のCZPチャートを使用し、CZPチャートを中心部から上下に二等分する水平線で切断し、切断された上側の上辺と他方の下辺とを接合して作成され、これにより、空間周波数特性が中心部において最も高く、垂直方向周辺部へ向かうにつれて低くなるようにしている。これら水平解像度評価用と垂直解像度評価用のテストチャートにより、空間周波数が中心部において最も高く、水平方向周辺部あるいは垂直方向周辺部へ向かうにつれて低くなるようにすることにより、テストチャートの観測波形の変化が、ビデオカメラのレンズと本体のいずれによるものかの判別をしやすくしている。

【発明の概要】

【発明が解決しようとする課題】

【0003】

しかし、特開２０１１−１０１１００号公報に開示されたテストチャートでは、テストチャートの中心から水平方向に延びたラインと、中心から垂直方向に延びたラインの２つのラインにおいてしか、空間解像度を評価することができない。

【0004】

したがって、本発明の目的は、撮像装置の空間解像度を、撮像装置の視野内のより多くの部分において判定することができるテストチャートを提供することである。
また、本発明の目的は、上記テストチャートを使用して撮像装置の空間解像度を測定する測定方法を提供することである。

【0005】

本発明のさらに別の目的は、上記測定方法を実施するためのコンピュータ読み取り可能記憶媒体及びコンピュータ・プログラムを提供することである。

【課題を解決するための手段】

【0006】

本発明の１実施形態にしたがい提供するテストチャートは、少なくとも１つのチャートを含むテストチャートであって、前記チャートが、同心円状に配列された複数のリングであって、該複数のリングが、周辺から中心に向かう方向において増大する空間周波数を有する、複数のリングと、該複数のリングの外に配置された少なくとも１つの１対のマークであって、前記１対のマークが、前記複数のリングの中心を通る横方向または縦方向において配置された、少なくとも１つの１対のマークと、を含む。

【0007】

また、本発明の別の実施形態にしたがい提供する測定方法は、ａ）被測定対象の撮像装置からの、請求項１〜５のいずれかに記載の前記テストチャートのデジタル撮像信号を受けるステップと、ｂ）前記デジタル撮像信号から、前記テストチャートの前記少なくとも１つのチャートのうちの１つのチャートの前記１対のマークを検出することにより、前記１つのチャートのある角度における測定エリアを定めるステップと、ｃ）前記測定エリア内において低い空間周波数から高い空間周波数に向かって複数の測定ラインを定めるステップと、ｄ）各前記測定ラインの各々において、前記デジタル撮像信号から、前記測定ラインに沿った一連のデジタル・データを取得し、該一連のデジタル・データが定める波形の周期および振幅の少なくとも１つに基づいて各前記測定ラインにおける空間解像度を求めるステップと、ｅ）複数の前記測定ラインにおいてそれぞれ求めた複数の空間解像度を使用して、前記被測定対象の撮像装置の、前記測定エリアにおける空間解像度を求めるステップと、を含む。

【0008】

また、本発明のさらに別の実施形態にしたがって、上記のテストチャートおよび測定方法に関して種々のより具体的な特徴を提供する。

【図面の簡単な説明】

【0009】

【図1】図１は、本発明の１実施形態によるテストチャート１００の全体を示す平面図である。

【図2】図２（Ａ）は、図１のテストチャートのうちの中央チャートの上側半分を拡大して示す図であり、図２（Ｂ）は、中央チャートの一部をさらに拡大した図であり、図２（Ｃ）は、図２（Ｂ）の拡大部分が定める輝度の波形を示す図。

【図3】図３は、本発明の１実施形態のテストチャートを使用して撮像装置の空間解像度を評価する測定システムを示す図。

【図4】図４は、図３の測定装置で実行する、撮像装置の空間解像度を評価するための測定フローを示すフローチャート。

【図5】図５は、図４の測定フローのうちの空間解像度測定ステップを詳細に示す１実施形態の測定処理を示すフローチャート。

【図6】図６（Ａ）は、静止画のうちの複数の測定ラインにおける画像を示し、図６（Ｂ）は、それら複数の測定ラインのうちの異なった２つの測定ラインの画像が定める輝度の波形を示す図。

【図7】図７は、図６の波形図をさらに詳細に示す図。

【図8】図８（Ａ）は、図５の測定処理において得られるある測定ライン内のデジタル・データの値を示し、図８（Ｂ）は、この測定ライン内の一部のデジタル・データの値をより詳細に示す図。

【図9】図９（Ａ）、（Ｂ）は、図４の測定フローの結果表示ステップにより表示される、指定角度間隔が５度の場合の測定結果の表示例を示す図。

【図10】図１０（Ａ）、（Ｂ）は、指定角度間隔が５度の場合のさらに別の測定結果表示例を示す図。

【図11】図１１（Ａ）、（Ｂ）は、指定角度間隔が１度の場合の図９（Ａ）、（Ｂ）と同様の測定結果表示例を示す図である。

【図12】図１２（Ａ）、（Ｂ）は、指定角度間隔が１度の場合の図１０（Ａ）、（Ｂ）と同様の測定結果表示例を示す図である。

【図13】図１３は、図４の解像度測定ステップの別の実施形態の測定処理を示すフローチャート。

【図14】図１４は、図１３の測定処理において利用する周期を説明するための波形図。

【図15】図１５は、図４の解像度測定ステップのさらに別の実施形態の測定処理を示すフローチャート。

【発明を実施するための形態】

【0010】

以下、本発明の種々の実施形態について、図面を参照して詳細に説明する。
図１は、本発明の１実施形態によるテストチャート１００の全体を示す平面図である。図示の通り、テストチャート１００は、矩形の形状を有し、１６：９アスペクト比（ＨＤＴＶ用）を有している。尚、被測定対象の撮像装置のアスペクト比に合わせてその他のアスペクト比にすることも可能である。図示のように、テストチャート１００は、矩形枠１０２の中心に配置した中央チャート１１０と、矩形枠１０２の４隅にそれぞれ配置した４つの周辺チャート１２０〜１８０とを有している。矩形枠２の縦のサイズと比べると、中央チャート１１０は、１／３程度のサイズであり、各周辺チャート１２０〜１８０は、１／５程度のサイズである。

【0011】

中央チャート１１０の上側半分を拡大した図２も参照して説明すると、図２（Ａ）に示すように、中央チャート１１０は、中央に配置した同心円状の複数のリング１１１と、このリングの周囲に配置した４つの位置検出マーク１１２，１１３，１１４、１１５と、リング１１１の中央にある中心円形領域１１６とを有する。同心円状の複数のリング１１１は、黒と白の間で輝度が連続的に変化することにより、周辺から中心に向かって空間周波数が高くなる同心円状の模様を定める。空間周波数は、３５０ＴＶ本（外側）から１１５０ＴＶ本（中心側）へと変化し、変化率は１０ＴＶ本／サイクルである。ＴＶ本とは画面の垂直方向の高さに等しい水平幅の中に、白黒の縦線が何本まで分離して見分けられるかを表す単位である。この際、白と黒の縦線をそれぞれ1本と数える。1本の走査線においては、白黒の縦線は輝度信号が矩形波（あるいは正弦波）の1サイクルに相当する。輝度信号の山が白、谷が黒として表示されるので、1サイクルはTV本で2本分の期間に相当する。図２（Ｂ）では、２サイクル分のリングをより詳細に示している。尚、図２（Ｂ）では、段階的にグレースケールが変化するように図示しているが、実際は、連続的に変化するグレースケールを使用している。図２（Ｃ）は、２サイクル分のリングが生成する輝度を、０〜１の値に正規化して示しており、輝度が関数ｓｉｎ^２ｘで連続的に変化している。

【0012】

中央チャート１１０に含まれた位置検出マーク１１２，１１３，１１４、１１５は、チャートの位置とその向きを自動で判定するために使用するマークである。また、マーク１１２と１１４はチャートの中心を通る水平方向（０度と１８０度の角度を示す）を示し（０度と１８０度の角度を示す）、そしてマーク１１３と１１５はチャートの中心を通る垂直方向を示す（９０度と２７０度を示す）。また、図２（Ａ）に示すように、中央チャート１１０に黒の中心円形領域１１６を含ませることによって、幅が２ｄの中央帯状部分１１８を定める。２ｄの幅は円形領域１１６の直径よりも小さくすることによって、中心側の端部を比較的直線に近い線にすることができ、これによって中心部でのリングの配置を回避して測定誤りの発生を減少させるようにする。この帯状部分１１８の領域からは、空間解像度を測定するための測定ラインとして使用する複数のラインを選択する。測定ラインとして複数のラインを使用する理由は、ライン毎の空間解像度の測定結果にバラツキが生ずるので、その影響を少なくするためである。

【0013】

周辺チャート１２０〜１８０は、互いに同じ構成のものであり、しかも中央チャート１１０と同様の構成を有している。各周辺チャートが中央チャートと異なる点は、周辺チャートの全体のサイズが小さく、また中心に配置された複数のリングが小さい点である。各周辺チャートの複数のリング（例えば周辺チャート１２０の複数のリング１２２）は、中央チャートと同じ３５０ＴＶ本（外側）〜１１５０ＴＶ本（中心側）の空間周波数を有するが、変化率は５０ＴＶ本／サイクルで、中央チャートと比べ１サイクル当たりの変化率は高い。

【0014】

以上に説明したテストチャート１においては、中心が黒の白黒配置であるが、この配置を反転させることもできる。
尚、以上では、テストチャートは、白黒のチャート、すなわち、第１の色である黒と第２の色である白との２つの色の間で階調すなわちグレースケールが変化する模様を構成するものとして説明した。白黒のチャートでは、赤（R）、緑（G）、青（B）の各コンポーネント（「チャンネル」とも呼ぶ）値は同じ値になるため、RGBのどれか一つのチャンネルのデータを空間周波数を定める輝度データとして使用すればよい。また、テストチャートには、白黒以外のその他の色を使用することもできる。例えば、第１と第２の色として黒白以外の任意の２つの色を用いることができる。下記の表は、それぞれの色を表現するときの、代表的な８色のRGBチャンネルの各々の値（８ビット表現）を示している。

【0015】

【表1】

【0016】

尚、表１は、上記代表的な８色を示しているに過ぎず、黒白以外の任意の２つの色として、その他の色を使用することも可能である。
黒白以外の任意の２つの色を用いるとき、正と負のピークの間ではそれら２つの色の間で色の階調を変化させる。これによって、色による解像度の分布の変化を解析することができる。また、階調の変化は、RGBチャンネルの３つのチャンネルの値を同時に変化させたり、３つのチャンネルのうちの２つのチャンネルを同時に変化させたり、あるいは１つのチャンネルのみを変化させたりすることができる。階調の変化のさせ方としては、例えばグラデーションが可能である。さらに、RGBチャンネルの３つのチャンネルを同相で変化させるだけでなく、互いに位相（例えば９０度等）をずらして変化させることもできる。階調の変化に位相差をもたせることにより、空間解像度の測定精度を高めることができる。

【0017】

以上に説明したテストチャート１００によれば、矩形枠２内に複数のチャート１１０，１２０〜１８０を配置することにより、視野内の異なった領域の空間解像度を互いに独立して評価することができる。また、各チャートが同心円状のリングを含むことにより、それら各領域において、複数の異なった方向の空間解像度を評価することができる。これにより、視野内の希望する任意の数の領域の空間解像度を、０度から３６０度の全ての方向において評価することができる。

【0018】

次に、図３を参照して、上記のテストチャート１を使用して撮像装置の空間解像度を評価する測定システム３００について詳細に説明する。図示の通り、測定システム３００は、図１に示したテストチャート１００を撮影する、被測定対象の撮像装置３１０と、この撮像装置３１０からの撮像出力信号を受ける測定装置３２０とを含んでいる。被測定対象の撮像装置３１０の例としては、デジタルカメラ、ビデオカメラ、その他の任意のカメラ機構を有するデバイス、例えばスマートフォン、携帯電話や、それらカメラ機構のモジュール等が含まれる。測定装置３２０は、測定専用の機器で構成したり、あるいはコンピュータのような汎用性のある機器で実現することができる。以下の説明では、測定装置をパーソナル・コンピュータＰＣで実施した例で説明する。このＰＣは、プロセッサ、入力装置、コンピュータ・プログラムを格納した記憶装置、ディスプレイ等の出力装置という公知の要素を備えている。尚、コンピュータ・プログラムは、コンピュータ読み取り可能記憶媒体に格納することもできる。

【0019】

図４を参照して、撮像装置３１０の空間解像度を評価するための測定装置３２０で実行するコンピュータ・プログラムで実施される測定フローを説明する。まず、最初のステップ４００で、撮像装置３１０から、テストチャート１００を撮影した結果の出力信号を受けて、ＰＣ内のバッファに置く。この出力信号は、例えば、カメラ映像の静止画ファイルであり、ＪＰＥＧ，ＢＩＴＭＡＰ，ＰＮＧ，ＴＩＦＦなどの任意のフォーマットが可能である。尚、静止画ファイルを受け取る方法としては、撮像装置３１０から直接受ける場合の他に、既に取得された静止画ファイルを、ネットワーク（例えばイーサネット）や、キャプチャボード及びＵＳＢなど記憶媒体から受けることもできる。次に、ステップ４０２において、受けた静止画における空間周波数の測定範囲を定める。すなわち、例えば５つのチャート１１０，１２０〜１８０のうちの中央チャート１１０を選択し、さらにこの選択した中央チャート１１０の空間解像度を測定する方向を角度で指定する。例えば、チャートの中心から右水平方向を示す０度の角度を、測定開始角度として指定し、さらに、測定する角度間隔を、０．０１度〜９０度の範囲内で指定する。尚、本実施形態では、４方向同時に測定するので、０度を測定開始角度として指定した場合、９０度、１８０度、２７０度での測定も指定したことになる。

【0020】

次に、ステップ４０４において、空間解像度の測定エリアを検出する。具体的には、上下左右４つの位置検出マーク１１２〜１１５を囲んでいる四角形の中で左から右までライン毎に走査し、白い四角形であるマーク１１２〜１１５の座標を確定する。次に、上の白い四角形すなわちマーク１１３の上側の左右の角の座標と、下の白い四角形すなわちマーク１１５の下側の左右の角の座標とで、縦方向の測定エリアを定める。同様に、左の白の四角形すなわちマーク１１４の左側の上下の角の座標と、右の四角形すなわちマーク１１２の右側の上下の角の座標とで、横方向の測定エリアを定める。次に、最初のリングのピークを求める。尚、空間解像度のより正確な測定のため、測定エリアは、複数の測定ラインを含むように定める。また、測定エリア内のどの測定ラインを使用するのかも定める。

【0021】

次に、ステップ４０６で、定めた測定ラインにおいて空間解像度の測定を行い、測定結果を格納する。この測定の詳細は、図５を参照して後述する。次に、ステップ４０８で、測定エリアでの測定が完了したかどうか判定し、完了していないときは、ステップ４１０で次の測定エリアを選択する。尚、ステップ４１０の処理の詳細については、後述する。次に、ステップ４０６の空間解像度測定を再開する。もしステップ４０８で測定が完了していると判定したとき、ステップ４１２で、格納された測定結果を表示し、このフローを終了する。

【0022】

次に、図５を参照して、空間解像度測定ステップ４０６の１実施形態の測定処理の詳細について説明する。図５の測定処理では、波形における正と負のピークの発生における異常を利用する。まずステップ５００において、測定バッファを初期化して、測定用パラメータ（ライン・カウンタ、検査データ・ポインタ、結果記録用バッファなど）を初期化する。次に、ステップ５０２において、４方向の測定エリアの各々内の測定ラインを選択する。測定エリア内での測定ラインの選択は、基本的には、左側から右側へ、上側から下側への順で行い、ライン・カウンタで記録する。各測定エリア内では、４０個のラインを選択してその平均値を測定値とする。次に、ステップ５０４で、選択した測定ラインのデータの中の高周波ノイズを、例えばガウスフィルタを使って除去する。具体的には、測定画像がＨＤ（１９２０ｘ１０８０）の場合、最高解像度の１周期は２ｄｏｔ(その周波数は７４．２５ＭＨｚの半分)なので、１周期２ｄｏｔ以下の周波数成分はノイズとみなして除去する。

【0023】

図６は、このようにして得られた画像（図６（Ａ））とその画像のデータの波形（図６（Ｂ））を示している。この図では、水平ライン５３５と水平ライン５４０を含む複数のラインの画像を示しており、下側の波形図には、水平ライン５３５と水平ライン５４０の輝度の変化を０〜１の範囲に正規化して示している。図７は、図６の波形図を拡大したものであり、水平ライン５３５と水平ライン５４０とでは、若干サンプリング点がずれていることが分かる。このため、上述した通り、本発明の１実施形態では、各測定エリア内において、複数の測定ラインで空間解像度を測定する。

【0024】

次に、ステップ５０６で、画像データの前処理を行うことによって、山と谷のピークのみをもつデータにする。具体的には、（１）画像データの波形のすべてのサイクルで山と谷のピークを特定し、（２）それぞれのサイクルで、山と谷の両ピークの間にあるデータのうち、立ち上がり途中のデータ、すなわち立ち上がり開始時の谷のピークから次の山のピークまでの間にあるデータは、すべてその立ち上がり開始時の谷のピークのデータ値に変換し、（３）それぞれのサイクルで、山と谷の両ピークの間にあるデータのうち、立ち下がり途中のデータ、すなわち、立ち下がり開始時の山のピークから次の谷のピークまでの間にあるデータは、すべてその立ち下がり開始時の山のピークのデータ値に変換する。次に、ステップ５０８において、前処理されたデータの平均値を求め、この平均値を、山と谷のピークを判定するためのしきい値と定め、またこのしきい値を「スライダ」の値とする。

【0025】

図８は、ステップ５０６と５０８の処理によって得られたデータを示している。詳細には、図８（Ａ）は、一例として０度と１８０度の２つの測定ラインを含む１つのラインにおける輝度データ１〜５６０のデジタル値を示しており、図８（Ｂ）は、図８（Ａ）において四角で囲んだ輝度データ１５７〜１７７を拡大して示している。尚、図中、本発明の理解のため、上記で求めた値をもつスライダを示している。図から分かるように、山と谷のピーク周辺以外のデータは含まれていない。

【0026】

次に、ステップ５１０において、測定ライン内の検査対象のデータをポイントする検査データ・ポインタの初期化して、検査を開始すべき最初のデータを指すようにする。具体的には、測定ラインの角度が０度の時は右端、９０度の時は上端、１８０度の時は左端、２７０度の時は下端のデータを指すようにする。例えば、図８（Ａ）の例では、検査データ・ポインタは、０度の測定ラインが選択されたときにはデジタル・データ５６０をポイントし、１８０度の測定ラインが選択されたときにはデジタル・データ１をポイントする。

【0027】

次に、ステップ５１２において、検査データ・ポインタが、チャートの中心、すなわち、選択された測定ラインの検査開始位置とは反対側に到達したかどうか判定する。この判定は、黒のデータがある一定数以上続くかどうかを検出することにより行う。到達したと判定したときには、選択された測定ラインでの検査が完了したと判断してステップ５２０に進む。到達していないと判定したときには、ステップ５１４に進む。

【0028】

ステップ５１４においては、ポインタが指したデータが波形の立下り部分に該当するかどうか判定する。立ち下がり部分の判定は、直前に検査したデータの値より小さいかどうか判定することにより行う。立ち下がり部分でない、すなわち直前に検査したデータと同じかあるいは大きい値であると判定したときには（図８（Ｂ）の例ではデータ８００や８０１）、ステップ５１６に進んで、検査データ・ポインタを中心方向に１つシフトさせてステップ５１２に戻る。一方、ステップ５１４において、立ち下がり部分であると判定したときには（図８（Ｂ）の例ではデータ８０２や８０４）、次のステップ５１７において、検査中のデータの値が、スライダ（しきい値を示す）より上であるかどうか判定する。スライダより上であると判定したとき（図８（Ｂ）の例ではデータ８０４）、スライダを超える谷が現れないこと（白黒がつぶれて判別できない）から空間解像度の上限を超えたと判定し、直前に検査したデータ（図８（Ｂ）の例ではデータ８０１）の位置に対応する空間解像度を上限値として判断し、その直前検査データのポインタの値を記録する。

【0029】

ステップ５１７でＹｅｓの場合に、それぞれ判定した空間解像度を、検査中の測定ラインでの測定結果として確定する。具体的には、ポインタの値をチャート上の位置に換算して、換算したチャート上位置に対応するTV本の本数を算出する。この算出したＴＶ本が空間解像度を表す。次に、ステップ５２２においてこの測定結果を指定のバッファに保存する。

【0030】

尚、本実施形態では、ステップ５１４での立ち下がり検出と、ステップ５１７でのスライダより上か否かの判定の組み合わせを使用しているが、他の方法も可能である。例えば、ステップ５１４で、立ち上がりを検出し、ステップ５１７で、スライダよりも下か否かの判定を行うこともできる。

【0031】

次に、ステップ５２４において、４方向（０度、９０度、１８０度、２７０度）の測定ライン全てについて測定を完了したかどうか判定し、完了していないときは、次の測定ラインを選択してステップ５１０に進み、その後ステップ５１２〜５２２のステップを繰り返す。一方、４方向の測定ライン全てに関して測定が完了したときには、ステップ５２６で、４つの測定エリアの各々内のすべての測定ラインでの測定を完了したかどうか判定し、完了していないときは、ステップ５０２に進んで、４つの測定エリア内の各々の次の測定ラインを選択する。例えば、図６の測定ライン５３５での測定が完了したが、測定ライン５４０の測定をしていないときには、測定ライン５４０を選択する。この後、上述のステップ５０４〜５２４を繰り返す。各測定エリア内の全ての測定ラインでの測定が完了したとき、ステップ５２８で、各測定エリア内の複数の測定エリアでの測定結果の平均値を求め、この平均値であるＴＶ本値を、該当の測定エリアにおける最終的な空間解像度として決定し、結果バッファに保存する。これにより、４方向における最終的な空間解像度が得られる。

【0032】

この後、図４のステップ４０８に進み、全ての測定エリアにおける測定が完了していないときには、ステップ４１０で次の測定エリアを選択する。すなわち、例えば、０度、９０度、１８０度、２７０度の測定エリアでの測定が完了したときには、指定された角度間隔だけ離れた次の４つの角度で測定エリアを選択する。例えば、指定された角度間隔が５度の場合には、５度、９５度、１８５度、２７５度の測定エリアが選択される。指定された角度間隔が１度のときには、１度、９１度、１８１度、２７１度の測定エリアが選択される。ここで、次の測定エリアの選択時には、測定エリアに対応する空間領域を定める座標を、チャートの中心を軸として指定角度間隔だけ空間的に回転させる演算を行い、そしてその結果としての回転後の空間領域の座標を次の測定エリアを定めるものとして得られる。尚、座標の空間回転は、良く知られて方法で可能であるが、例えば下記の式を使うことによって可能である。

【0033】

【数1】

【0034】

ここで、Ｆ（ｘ、ｙ）は元座標、Ｆ’（ｘ’，ｙ’）は目標座標、σは回転角度である。ＦとＦ’の原点は重なるので、Ｆ（０，０）＝Ｆ‘（０，０）である。
次に、ステップ４０８で、全ての測定エリアで測定が完了すると、ステップ４１２に進み、結果バッファ内に格納された測定結果を表示する。測定結果の１例（指定角度間隔＝５度）は下記の表に示す（中間の値は省略）。

【0035】

【表2】

【0036】

尚、ｘ座標は、空間解像度の上限位置を示すチャート上の位置を表すｘ軸（横軸）の座標であり、ｙ座標は、その位置のｙ軸（縦軸）の座標である。この測定結果は、上記の表形式で結果を表示する方法以外に、以下で説明するような種々の異なった形式で測定結果を表示することができる。

【0037】

図９（Ａ）は、横軸が測定エリアの角度、縦軸が空間解像度を表すＴＶ本の値を示す図である。図９（Ｂ）は、横軸が上記のｘ座標で、縦軸が上記のｙ座標である。したがって、図９（Ｂ）の結果表示は、その元となるチャートと空間的に対応している。

【0038】

次に、図１０（Ａ）は、極座標形式の表示であって、図９（Ｂ）の結果表示を「ｘ」のプロットとしてテストチャートに重ねて表示しており、空間解像度がチャート上のどの位置にあるのかを示している。したがって、３６０度に渡る空間解像度を、チャートに関連付けて把握することができる。また、図１０（Ｂ）も極座標形式での表示であって、ある角度でのＴＶ本の値を、その角度での極からの距離によって示している。したがって、極からの線の長さが空間解像度の大きさを表しており、中心は０ＴＶ本、白い円の円周部(縁)が１１５０ＴＶ本を示す。この表示形式でも、３６０度に渡る空間解像度の分布を把握を要にすることができる。

【0039】

図１１（Ａ）、図１１（Ｂ）、図１２（Ａ），図１２（Ｂ）は、図９（Ａ）、図９（Ｂ）、図１０（Ａ），図１０（Ｂ）とそれぞれ同様の図であって、指定角度間隔をより細かな１度にしたときの結果を示している。

【0040】

次に、図１３を参照して、図４の解像度測定ステップ４０６の別の実施形態である測定処理を説明する。この実施形態は、空間解像度の判定に波形における周期の長さを利用する。図示のように図１３の測定処理は、図５の測定処理と類似しており、異なっているのは、二値化ステップ１３０１を含み、また図５のステップ５０６とステップ５１２〜５２０の代わりに、ステップ１３００とステップ１３０２〜１３１０を含む点である。

【0041】

詳細には、ステップ１３００では、図５のステップ５０６の前処理に加え、山のピークから次の谷のピークまでは「谷データ」とし、谷のピークから次の山のピークまでは「山データ」とする。あるいは変わりの方法として、山のピークから次の谷のピークまでは「山データ」とし、谷のピークから次の山のピークまでを「谷データ」とすることもできる。いずれにしても、波形データを山データと谷データの２つに分ける。次に、ステップ１３０１で、この前処理されたデータに対し、スライダの値でスライスすることにより二値化し、これによって、図１４に示すような二値化波形を生成する。尚、点線の部分は、波形の１周期部分である。

【0042】

このような処理をした後、ステップ１３０２で、１周期を検出したかどうか判定する。一例として、最初に「谷データ」を検出したときは、１周期ではないため、ステップ１３０４に進み、ここで、検査データ・ポインタを中心側に１つ進めてステップ１３０２に戻る。ここで、再度１周期を検出したかどうか判定する。例えば、「山データ」を今回検出したときに、「谷データ」を前回検出しているため、１周期を検出したと判定する。このとき、ステップ１３０６に進み、ここで、１周期の検出が初回かどうか判定し、初回のときは、ステップ１３０８でその１周期の周期長をバッファに記録し、その後ステップ１３０４に進む。一方、ステップ１３０６において初回でないと判定したときは、ステップ１３１０で、これまでに検出した１周期、すなわちバッファに格納された周期長よりも長いかどうか判定する。長くない場合、すなわち短いか同じ長さの場合、今回検出した周期長をバッファに記録して周期長を更新した後、ステップ１３０４に進み、次の検査に移る。一方、ステップ１３１０において長いと判定したとき、今回検出した周期は、黒白がつぶれて区別できなくなったと判定し、直前の検査データ・ポインタの値から、チャート上の座標とこの座標の位置におけるＴＶ本の値を求める。この後、図５のステップ５２２以降のステップと同様のステップに進むことにより、最終的な測定値を生成することができる。

【0043】

また、図１５を参照して、図４の解像度測定ステップ４０６のさらに別の実施形態である測定処理を説明する。この実施形態は、空間解像度の判定に波形における振幅を利用する。図示のように図１５の測定処理は、図５の測定処理と類似しており、異なっているのは、図５のステップ５０８がなく、また、図５のステップ５０６及びステップ５１２〜５２０の代わりに、ステップ１５００及びステップ１５０２〜１５０６を含む点である。

【0044】

詳細には、ステップ１５００では、図５のステップ５０６の前処理に加え、図１３の処理と同様に、山のピークから次の谷のピークまでは「谷データ」とし、谷のピークから次の山のピークまでは「山データ」とする。あるいは変わりの方法として、山のピークから次の谷のピークまでは「山データ」とし、谷のピークから次の山のピークまでを「谷データ」とすることもできる。いずれにしても、波形データを山データと谷データの２つに分ける。このような処理をした後、ステップ１５０２において、１周期を検出したかどうか判定する。一例として、最初に「谷データ」を検出したときは、１周期ではないため、ステップ１５０４に進み、ここで、検査データ・ポインタを中心側に１つ進めてステップ１５０２に戻る。ここで、再度１周期を検出したかどうか判定する。例えば、「山データ」を今回検出したときに、「谷データ」を前回検出しているため、１周期を検出したと判定する。このとき、ステップ１５０６に進み、ここで、検出した１周期の振幅が、予め定めたしきい値より小さいかどうか判定する。しきい値より小さくない場合、白黒がつぶれずに検出されたと判断し、そしてステップ１５０４に進んで、次のデータに対する検査に移る。一方、しきい値より小さい場合、白黒がつぶれたと判断し、直前の検査データ・ポインタの値からチャート上の座標とこの座標の位置におけるＴＶ本の値を求める。この後、図５のステップ５２２以降のステップと同様のステップに進むことにより、最終的な測定値を生成することができる。

【0045】

以上に詳細に説明した本発明のテストチャートを使用すれば、空間的に異なった視野領域における空間解像度を各領域において３６０度内の異なった方向で評価することができ、これにより従来のものと比べ詳細な評価を可能にすることができる。また、本発明の測定方法によれば、上記テストチャートを使用して空間解像度を的確に評価し、その評価結果をユーザに分かりやすい形式で表示することができる。

【図3】

【図4】

【図5】

【図7】

【図9】

【図11】

【図13】

【図14】

【図15】

【図1】

【図2】

【図6】

【図8】

【図10】

【図12】