(57)【要約】【課題】簡易な方法で撮像装置の撮像素子の感度分布を測定し、撮像装置を測定器として用いて画像表示装置のキャリブレーションを行う方法を提供する。【解決手段】撮像素子の感度分布を測定する測定装置であって、被写体の像が形成される撮像素子上の位置、又は、被写体の像を撮像素子上に形成する合焦距離を変更可能な光学系の合焦距離が異なる、２つの撮像条件で前記撮像素子により被写体の像を撮像して得られる第１の画像データ及び第２の画像データを取得する取得手段と、前記第１の画像データと前記第２の画像データとに基づき前記撮像素子の感度分布を算出する算出手段と、を備える測定装置。【選択図】図１９
2015119344.tif 000002
【特許請求の範囲】
【請求項１】
撮像素子の感度分布を測定する測定装置であって、
被写体の像が形成される撮像素子上の位置、又は、被写体の像を撮像素子上に形成する合焦距離を変更可能な光学系の合焦距離が異なる、２つの撮像条件で前記撮像素子により被写体の像を撮像して得られる第１の画像データ及び第２の画像データを取得する取得手段と、
前記第１の画像データと前記第２の画像データとに基づき前記撮像素子の感度分布を算出する算出手段と、
を備える測定装置。
【請求項２】
前記算出手段は、前記第１の画像データと前記第２の画像データとの撮像条件の差異と、前記第１の画像データと、に基づき、撮像条件が前記第２の画像データの撮像条件と同じである第３の画像データを生成し、前記第２の画像データと前記第３の画像データとの差異に基づき前記撮像素子の感度分布を算出する請求項１に記載の測定装置。
【請求項３】
前記第１の画像データの撮像条件における前記被写体の像が形成される撮像素子上の位置は、前記第２の画像データの撮像条件における前記被写体の像が形成される撮像素子上の位置に対し、第１のずれ量だけ相違しており、
前記算出手段は、前記第１の画像データに対し、前記第１のずれ量に基づく変換処理を行うことにより、前記被写体の像が形成される撮像素子上の位置が前記第２の画像データの前記被写体の像が形成される撮像素子上の位置と同じである第３の画像データを生成する請求項２に記載の測定装置。
【請求項４】
前記第１の画像データと前記第２の画像データとは、前記撮像素子を備えた撮像装置により、同一の被写体を、被写体に対する撮像装置の位置を異ならせて撮像して得られた画像データである請求項３に記載の測定装置。
【請求項５】
前記第１のずれ量は、並行移動、回転、及び拡大縮小の組み合わせで表される量であり、
前記算出手段は、前記第１の画像データに対し、前記第１のずれ量を表す並行移動、回転、及び拡大縮小の組み合わせに対応するアフィン変換処理を行うことにより、前記第３の画像データを生成する請求項３又は４に記載の測定装置。
【請求項６】
前記算出手段は、前記第１の画像データと前記第２の画像データとを解析することにより前記第１のずれ量の情報を取得し、取得した情報に基づき前記変換処理を行う請求項３〜５のいずれか１項に記載の測定装置。
【請求項７】
前記算出手段は、画像データを構成する各画素について、前記第２の画像データの画素値と前記第３の画像データの画素値との比を求め、当該画素値同士の比に基づき、前記撮像素子を構成する各画素について、その画素の感度と、その画素に対し前記第１のずれ量だけ相違する位置にある画素の感度と、の比である第１の感度比を算出し、前記第１の感度比に基づき前記撮像素子の感度分布を算出する請求項３〜６のいずれか１項に記載の測定装置。
【請求項８】
前記算出手段は、前記第１の感度比を順次積算して前記撮像素子の前記第１のずれ量ずつ離れた画素の感度を算出することにより第１の離散的な感度分布を算出し、前記第１の離散的な感度分布に対し補間処理を行うことにより前記撮像素子の全ての画素の感度分布を算出する請求項７に記載の測定装置。
【請求項９】
前記取得手段は、前記第１の画像データ及び前記第２の画像データの撮像条件に対し被写体の像が形成される撮像素子上の位置が異なる第４の画像データを更に取得し、
前記第４の画像データの撮像条件における前記被写体の像が形成される撮像素子上の位置は、前記第２の画像データの撮像条件における前記被写体の像が形成される撮像素子上の位置に対し、第２のずれ量だけ相違しており、
前記第２のずれ量の大きさは前記第１のずれ量の大きさより小さく、
前記算出手段は、前記第４の画像データに対し、前記第２のずれ量に基づく変換処理を行うことにより、前記被写体の像が形成される撮像素子上の位置が前記第２の画像データの前記被写体の像が形成される撮像素子上の位置と同じである第５の画像データを生成し、画像データを構成する各画素について、前記第２の画像データの画素値と前記第５の画像データの画素値との比を求め、当該画素値同士の比に基づき、前記撮像素子を構成する各画素について、その画素の感度と、その画素に対し前記第２のずれ量だけ相違する位置にある画素の感度と、の比である第２の感度比を算出し、前記第１の感度比と前記第２の感度比とに基づき前記撮像素子の感度分布を算出する請求項７又は８に記載の測定装置。
【請求項１０】
前記算出手段は、前記第１の感度比を順次積算して前記撮像素子の前記第１のずれ量ずつ離れた画素の感度を算出することにより第１の離散的な感度分布を算出し、前記第１の離散的な感度分布における各画素の感度に前記第２の感度比を順次積算して前記撮像素子の前記第１のずれ量だけ離れた画素の間にある前記第２のずれ量ずつ離れた画素の感度を算出することにより第２の離散的な感度分布を算出し、前記第１の離散的な感度分布及び前記第２の離散的な感度分布を合わせて得られる第３の離散的な感度分布に対し補間処理を行うことにより前記撮像素子の全ての画素の感度分布を算出する請求項９に記載の測定装置。
【請求項１１】
前記第１の画像データの撮像条件と前記第２の画像データの撮像条件とは、前記合焦距離が異なり、
前記算出手段は、前記合焦距離の差と、第１の画像データと、に基づき、前記合焦距離が前記第２の画像データの合焦距離と同じである第３の画像データを生成する請求項２に記載の測定装置。
【請求項１２】
前記第１の画像データの撮像条件では、合焦距離が光学系と被写体との距離に等しく、
前記被写体は、画像表示装置に表示された、撮像素子の感度分布の測定用のパターン画像であり、当該パターン画像は、背景に背景と異なる色の点が配置されて構成された画像であり、
前記算出手段は、前記合焦距離の差に対応する前記光学系の点像分布関数と、前記第１の画像データの前記点に対応する位置の画素の画素値と、に基づき前記第３の画像データを生成する請求項１１に記載の測定装置。
【請求項１３】
前記算出手段は、画像データを構成する各画素について、前記第２の画像データの画素値と前記第３の画像データの画素値との比を求め、当該画素値同士の比に基づき、前記撮像素子を構成する各画素の感度分布を算出する請求項１１又は１２に記載の測定装置。
【請求項１４】
前記取得手段は、背景における点の配置パターンが異なる複数のパターン画像について前記第１の画像データと前記第２の画像データを取得し、
前記算出手段は、前記複数のパターン画像の各々について取得された前記第１の画像データから前記第３の画像データを生成し、前記複数のパターン画像の各々について取得された前記第２の画像データの画素値と前記第３の画像データの画素値との比を求め、当該画素値同士の比に基づき、前記撮像素子を構成する画素のうち各画像における点の位置に対応する画素及びその周囲の画素の感度分布を算出し、前記複数の画像の各々について算出された前記感度分布を繋ぎ合わせることにより前記撮像素子の全ての画素の感度分布を
算出する請求項１２又は１３に記載の測定装置。
【請求項１５】
請求項１〜１４のいずれか１項に記載の測定装置と、前記撮像素子と、を備えた撮像装置。
【請求項１６】
画像表示装置のキャリブレーションを行うキャリブレーション装置又は画像表示装置に対し前記測定装置により測定された撮像素子の感度分布の情報を送信する送信手段を更に備えた請求項１５に記載の撮像装置。
【請求項１７】
前記送信手段は、画像表示装置に表示されたキャリブレーション用のパターン画像を前記撮像素子により撮像して得られた画像データを更に送信する請求項１６に記載の撮像装置。
【請求項１８】
請求項１〜１４のいずれか１項に記載の測定装置と、
前記撮像素子と、
前記測定装置により測定された前記撮像素子の感度分布の情報に基づき、画像表示装置に表示されたキャリブレーション用のパターン画像を前記撮像素子により撮像して得られた画像データから、前記撮像素子の感度分布の影響を除去する補正を行う補正手段と、
前記補正手段により補正された画像データに基づき前記画像表示装置のキャリブレーションを行うキャリブレーション手段と、
を備える撮像装置。
【請求項１９】
請求項１〜１４のいずれか１項に記載の測定装置により測定された撮像素子の感度分布の情報を取得する第１の取得手段と、
画像表示装置に表示されたキャリブレーション用のパターン画像を前記撮像素子を備えた撮像装置により撮像して得られた画像データを取得する第２の取得手段と、
前記感度分布の情報に基づき、前記画像データから前記撮像素子の感度分布の影響を除去する補正を行う補正手段と、
前記補正手段により補正された画像データに基づき前記画像表示装置のキャリブレーションを行うキャリブレーション手段と、
を備える画像表示装置のキャリブレーション装置。
【請求項２０】
請求項１９に記載のキャリブレーション装置を備えた画像表示装置。
【請求項２１】
撮像素子の感度分布を測定する測定装置の制御方法であって、
被写体の像が形成される撮像素子上の位置、又は、被写体の像を撮像素子上に形成する合焦距離を変更可能な光学系の合焦距離が異なる、２つの撮像条件で前記撮像素子により被写体の像を撮像して得られる第１の画像データ及び第２の画像データを取得する取得工程と、
前記第１の画像データと前記第２の画像データとに基づき前記撮像素子の感度分布を算出する算出工程と、
を有する測定装置の制御方法。
【請求項２２】
前記算出工程では、前記第１の画像データと前記第２の画像データとの撮像条件の差異と、前記第１の画像データと、に基づき、撮像条件が前記第２の画像データの撮像条件と同じである第３の画像データを生成し、前記第２の画像データと前記第３の画像データとの差異に基づき前記撮像素子の感度分布を算出する請求項２１に記載の測定装置の制御方法。
【請求項２３】
請求項２１又は２２に記載の測定装置の制御方法により測定された撮像素子の感度分布
の情報を取得する第１の取得工程と、
画像表示装置に表示されたキャリブレーション用のパターン画像を前記撮像素子を備えた撮像装置により撮像して得られた画像データを取得する第２の取得工程と、
前記感度分布の情報に基づき、前記画像データから前記撮像素子の感度分布の影響を除去する補正を行う補正工程と、
前記補正工程により補正された画像データに基づき前記画像表示装置のキャリブレーションを行うキャリブレーション工程と、
を有する画像表示装置のキャリブレーション方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、撮像素子の感度分布の測定装置及びその制御方法、画像表示装置のキャリブレーション装置及びその制御方法に関するものである。
【背景技術】
【０００２】
画像表示装置のキャリブレーションでは、画面の輝度や色度をキャリブレーション用の測定器で測定し、測定される輝度や色度が目標値になるように、バックライトの発光制御や画像処理を行う。キャリブレーション用の測定器は、画面の表示領域の一部である測定対象領域の輝度や色度を測定するものが一般的である。しかし、キャリブレーション用の測定器は、ユーザが用意できない場合もある。
【０００３】
また、画像表示装置の輝度や色度は、画面内で均一でない。これはバックライトの発光特性のばらつきや経時劣化等による。キャリブレーション用の測定器を用いたキャリブレーションによって測定対象領域の輝度や色度を目標値に調整しても、画面全体の輝度や色度にむらがある場合、画面全体の輝度や色度を目標値に調整することは難しい。
【０００４】
そこで、広く一般に普及しているデジタルカメラやカメラ付スマートフォンなどの撮像素子を用いた撮像装置を、キャリブレーション用の測定器として用いることが考えられる。
【０００５】
撮像素子の感度は画素毎にばらつく場合がある。撮像装置で画像表示装置の画面を撮像した場合、撮像素子の感度分布の影響を反映した画像データが得られることになる。そのため、この画像データを用いて画像表示装置のキャリブレーションを精度良く行うためには、画像データから撮像素子の感度分布の影響を除去し、画面の表示輝度分布を正確に反映した画像データになるよう補正する必要がある。そのためには、撮像素子の感度分布の情報が必要になる。撮像素子の感度分布の情報は、例えば、撮像素子に光が均一の照度で照射する状態で撮像を行って得られる画像データに基づき取得することができる。この方法は、均一性の保証された基準面光源をユーザが用意する必要があり、簡便な方法とは言えない。
【０００６】
特許文献１には、基準面光源を用いずに、撮像素子の感度分布の情報を取得する方法が記載されている。特許文献１には、レンズ交換可能なカメラ本体に取り付け可能な筒状の補正用部材であって、入射光が入射する側の端面に光拡散部材を備え、筒の内面に光反射防止処理が施された補正用部材が開示されている。特許文献１では、この補正用部材をカメラ本体に取り付けることで、基準面光源を用いることなく、撮像素子に光が均一の照度で照射する状態を作り出すことを図っている。
【先行技術文献】
【特許文献】
【０００７】
【特許文献１】特開平１０−２００７９６号公報
【発明の概要】
【発明が解決しようとする課題】
【０００８】
しかしながら、特許文献１の技術も、ユーザが特別な補正用部材を用意する必要があり、かつ撮像装置がレンズ交換可能なカメラに限定されるという点で、やはり簡便な方法とは言えない。また、撮像素子に照射する光の照度を精度良く均一にするためには、補正用
部材に入射する光が均一である必要があり、暗室や照明が厳密に管理された室内等の環境が必要になるため、その点でも簡便な方法とは言えない。
【０００９】
そこで本発明は、簡易な方法で撮像装置の撮像素子の感度分布を測定し、撮像装置を測定器として用いて画像表示装置のキャリブレーションを行う方法を提供することを目的とする。
【課題を解決するための手段】
【００１０】
本発明は、撮像素子の感度分布を測定する測定装置であって、
被写体の像が形成される撮像素子上の位置、又は、被写体の像を撮像素子上に形成する合焦距離を変更可能な光学系の合焦距離が異なる、２つの撮像条件で前記撮像素子により被写体の像を撮像して得られる第１の画像データ及び第２の画像データを取得する取得手段と、
前記第１の画像データと前記第２の画像データとに基づき前記撮像素子の感度分布を算出する算出手段と、
を備える測定装置である。
【００１１】
本発明は、撮像素子の感度分布を測定する測定装置の制御方法であって、
被写体の像が形成される撮像素子上の位置、又は、被写体の像を撮像素子上に形成する合焦距離を変更可能な光学系の合焦距離が異なる、２つの撮像条件で前記撮像素子により被写体の像を撮像して得られる第１の画像データ及び第２の画像データを取得する取得工程と、
前記第１の画像データと前記第２の画像データとに基づき前記撮像素子の感度分布を算出する算出工程と、
を有する測定装置の制御方法である。
【発明の効果】
【００１２】
本発明により、簡易な方法で撮像装置の撮像素子の感度分布を測定し、撮像装置を測定器として用いて画像表示装置のキャリブレーションを行う方法を提供することができる。
【図面の簡単な説明】
【００１３】
【図１】撮像装置の内部構造の一例の略図を示す図
【図２】撮像装置によって表示装置の画面を撮像する様子を示す図
【図３】実施例１において測定用パターンを撮像装置により撮像する様子を示す図
【図４】撮像装置を移動させて測定用パターンを撮像する様子を示す図
【図５】パターンの撮像データから撮像素子の感度分布を推定する方法を示す図
【図６】パターンの撮像データから撮像素子の感度分布を推定する方法を示す図
【図７】撮像素子の感度分布を推定する処理を示すフローチャート
【図８】キャリブレーション用のパターンを撮像装置により撮像する様子を示す図
【図９】全白パターンの撮像データから画面の輝度むらを推定する方法を示す図
【図１０】撮像素子の感度分布の推定方法を示す図
【図１１】測定用パターンを撮像装置により撮像する様子を示す図。
【図１２】全白パターンの撮像データから画面の輝度むらを推定する方法を示す図
【図１３】実施例２において測定用パターンを撮像装置により撮像する様子を示す図
【図１４】ピントをぼかした状態で測定用パターンを撮像する様子を示す図
【図１５】ドットの配置を切り替えた測定用パターンを撮像する様子を示す図
【図１６】撮像素子の感度分布を推定する方法を示す図
【図１７】撮像素子の感度分布を推定する処理を示すフローチャート
【図１８】撮像素子の感度分布を推定する処理を示すフローチャート
【図１９】感度分布測定装置の構成を示す図
【図２０】感度分布測定装置の構成を示す図
【図２１】キャリブレーション装置の構成を示す図
【図２２】感度分布測定装置及びキャリブレーション装置の構成を示す図
【図２３】キャリブレーション装置の構成を示す図
【発明を実施するための形態】
【００１４】
（実施例１）
本発明に係る、撮像素子の感度分布を測定する測定装置の実施例１を説明する。実施例１では、被写体に対する撮像素子の相対位置を変化させることにより、撮像素子上に投影する光分布に変化を与える。具体的には、被写体の像が形成される撮像素子上の位置に変化を与える。それに伴う撮像素子出力の変化を比較することにより、撮像素子の感度分布を推定する。本実施例では、被写体は、撮像素子感度分布測定用のパターン画像を表示している表示装置の画面である。
【００１５】
図１に、撮像装置（カメラ）の内部構造の一例の略図を示す。
図１において、撮像装置１は、撮像素子１−１、レンズ等の光学系１−２、筐体１−３、ファインダー１−４を有する。
【００１６】
図２に、撮像装置１によって表示装置の画面の像を撮像素子に結像投影している様子を示す。図２において、撮像装置１は、表示装置２の画面２−１の像を、光学系１−２によって撮像素子１−１上に結像投影している。
【００１７】
図１９（Ａ）に、本実施例の感度分布測定装置１０の機能構成を示す。感度分布測定装置１０は、取得部１１、生成部１２、及び感度分布算出部１３を有する。
【００１８】
取得部１１は、被写体の像が形成される撮像素子上の位置が異なる２つの撮像条件で撮像素子１４により被写体を撮像して得られる第１の画像データ及び第２の画像データを取得する。本実施例では、第１の画像データの撮像条件における被写体の像が形成される撮像素子上の位置は、第２の画像データの撮像条件における被写体の像が形成される撮像素子上の位置に対し、第１のずれ量だけ相違しているものとする。
【００１９】
なお、図１９（Ａ）では取得部１１が撮像素子１４から直接、画像データを取得するように描かれているが、取得部１１が画像データを取得する方法はこれに限らない。例えば、撮像素子１４を備えるカメラから有線又は無線の通信手段により画像データを取得しても良いし、画像データが記録（記憶）されたメモリカード等のデータ格納メディアから読み取ることで画像データを取得しても良い。また、図１９（Ｂ）に示すように、撮像素子１４を備えるカメラ１５が、感度分布測定装置１０を備えており、カメラ１５内部で撮像素子１４から感度分布測定装置１０へ画像データが受け渡される構成でも良い。このような構成は、例えば、カメラ付きスマートフォンに感度分布測定装置１０の機能を有するハードウェアが備わる構成により実現できる。図１９（Ｂ）の構成は、或いは、ＣＰＵが実行することにより感度分布測定装置１０の機能が実現されるプログラム、ソフトウェア、アプリケーションとして実現できる。
【００２０】
生成部１２は、第１の画像データ及び第２の画像データの撮像条件の差異と、第１の画像データと、に基づき、撮像条件が第２の画像データの撮像条件と同じである第３の画像データを生成する。本実施例では、生成部１２は、第１の画像データに対し、第１のずれ量に基づく変換処理を行うことにより、被写体の像が形成される撮像素子上の位置が第２の画像データの被写体の像が形成される撮像素子上の位置と同じである第３の画像データを生成する。変換処理については後述する。
【００２１】
感度分布算出部１３は、第２の画像データと第３の画像データとの差異に基づき、撮像素子の感度分布を算出する。感度分布の算出方法については後述する。
【００２２】
感度分布測定装置１０は、算出した撮像素子の感度分布の情報を、キャリブレーション対象となるモニタ１７や、モニタ１７のキャリブレーションを行うキャリブレーション装置１８へ出力する。なお、図１９（Ｂ）に示すように、メモリカード１６等のデータ格納メディアに感度分布の情報を書き込み、メモリカード１６を用いてモニタ１７やキャリブレーション装置１８へ感度分布の情報を入力する構成が考えられる。また、図２０（Ａ）に示すように、有線又は無線で通信可能な送信部１９により感度分布の情報をモニタ１７やキャリブレーション装置１８へ送信する構成も考えられる。図２０（Ａ）では、カメラ１５が、撮像素子１４、感度分布測定装置１０、及び送信部１９を備えた構成を示すが、図２０（Ｂ）に示すように、撮像素子１４を備えるカメラ１５と、送信部１９を有する感度分布測定装置１０と、が別体構成でも良い。
【００２３】
図３に、画像表示装置（モニタ）の画面に撮像素子感度分布測定用のパターン画像を表示し、この表示された測定用のパターン画像を撮像装置によって撮像している様子を示す。この測定用パターンは、背景に背景と異なる色の点が配置された構成された画像であり、本実施例では背景は白色、点（ドット）は黒色である。図３（ａ）は表示装置及び撮像装置の配置図を示す。図３（ｂ）は撮像素子上に結像投影された像を示す。図３（ｃ）は撮像素子によって取得される画像データのうち、図３（ｂ）のｘ軸上の画像データｇ（ｘ）を示す（第２の画像データ）。図３（ｄ）は撮像素子のｘ軸上の感度分布ｓ（ｘ）を示す。
【００２４】
図３（ａ）に示すように、撮像装置は表示装置の画面の中央と正対して配置されているとする。従って、図３（ｂ）に示すように、画面中央に表示されているドットの像が撮像素子上の中央の位置に形成されるように、画面全体（パターン画像）の像が結像されている。このとき、ｘ軸上の画像データは、図３（ｃ）に示すように、パターン画像のドットに対応する部分（暗部である黒点）のデータ値は低く、それ以外の部分（明部である白背景）のデータ値は高くなっている。また、図３（ｄ）に示すように、撮像素子の感度は、ばらつきのために、ｘ軸上の位置ｘ０より右側（＋ｘ方向）で位置ｘ０より左側（−ｘ方向）より高くなっているものとする。この感度ばらつきを反映して、図３（ｃ）に示すｘ軸上の画像データの値も、位置ｘ０より右側の方が位置ｘ０より左側よりやや高くなっている。
【００２５】
図４に、表示装置はそのままで、撮像装置だけをｘ軸方向に−Ａだけ横移動させて、同一の被写体（測定用パターン画像）を撮像している様子を示す。図４（ｅ）は表示装置及び撮像装置の配置図を示す。図４（ｆ）は撮像素子上に結像投影された像を示す。図４（ｇ）は撮像素子によって取得される画像データのうち、上記図４（ｆ）のｘ軸上の画像データｇ’（ｘ）を示す（第１の画像データ）。
【００２６】
図４（ｆ）の撮像素子上に結像投影された像は、図３（ｂ）と比較して、全体にｘ軸方向にａだけ移動している。すなわち、第１の画像データの撮像条件における被写体の像が形成される撮像素子上の位置は、第２の画像データの撮像条件における被写体の像が形成される撮像素子上の位置に対し、第１のずれ量ａだけ相違している。ａは一般にはＡとは等しくない。光学系による結像投影が縮小投影又は拡大投影であればａはＡとは絶対値の異なる値となる。また、正立結像か倒立結像かによってaとＡの符号の異同は変わる。本
実施例では、説明を簡単にするため、撮像素子上の第１のずれ量ａと撮像装置の移動量Ａとの関係がａ＝−Ａである場合を例に説明するが、以下の説明は、必ずしもａ＝−Ａでなくとも成立する。
【００２７】
図４（ｇ）の第１の画像データｇ’（ｘ）においては、パターン画像のドットに対応する部分は、図３（ｃ）の第２の画像データｇ（ｘ）に対し、ｘ軸上でａだけ移動している。しかし、得られる画像データにおける撮像素子の感度分布（感度ばらつき）の影響の現れ方（本実施例では感度が若干上がる位置ｘ０）は、撮像を行う位置（撮像装置の移動）によらない。そのため、図４（ｇ）の第１の画像データｇ’（ｘ）において値が高くなる位置ｘ０は、図３（ｃ）の第２の画像データｇ（ｘ）に対し、移動していない。
【００２８】
このようにして得られた図３（ｃ）の第２の画像データｇ（ｘ）と図４（ｇ）の第１の画像データｇ’（ｘ）とを比較することにより、撮像素子の感度分布を推定することができる。推定方法を図５、図６を参照して説明する。
【００２９】
図５（ｃ）は、図３（ｃ）の画像データ、すなわち、撮像装置をｘ軸方向にＡだけ動かす前に撮像素子によって取得された第２の画像データｇ（ｘ）である。
第２の画像データｇ（ｘ）は、被写体の像を撮像素子上に形成する撮像装置の光学系により撮像素子上に結像された像の光の強度分布ｇdisplay（ｘ）と、撮像素子のｘ軸上の
感度分布ｓ（ｘ）と、の積ｇdisplay（ｘ）・ｓ（ｘ）である。
【数１】
2015119344.tif 000003
一般的には、光学系により結像された像の光の強度分布には、物体側（表示装置の画面）の輝度分布の他に、光学系固有の光学特性が影響するが、簡単のためここでは考えないものとする。従って、ｇdisplay（ｘ）は表示装置の画面のＸ軸上の輝度分布Ｂ（Ｘ）を
そのまま反映する強度分布になっているものとする。光学系固有の光学特性の影響については、別途、本実施例の末尾において述べる。
【００３０】
図５（ｇ）は、図４（ｇ）の画像データ、すなわち、撮像装置をｘ軸方向にＡだけ動かした後に撮像素子によって取得された第１の画像データｇ’（ｘ）である。
第１の画像データｇ’（ｘ）は、光学系により撮像素子上に結像された、ａだけｘ軸方向にずれた像の光の強度分布ｇdisplay（ｘ−ａ）と、撮像素子のｘ軸上の感度分布ｓ（
ｘ）と、の積ｇdisplay（ｘ−ａ）・ｓ（ｘ）である。
【数２】
2015119344.tif 000004
【００３１】
図５（ｈ）は、図５（ｇ）の第１の画像データｇ’（ｘ）を解析して移動量ａを求め、第１の画像データｇ’（ｘ）を逆方向の移動量−ａだけシフトさせて得た第３の画像データｇ’（ｘ＋ａ）である。この第３の画像データｇ’（ｘ＋ａ）は、生成部１２が、第１の画像データｇ’（ｘ）に対し、第１のずれ量ａに基づく変換処理（−ａのシフト）を行って生成した画像データである。この第３の画像データｇ’（ｘ＋ａ）は、被写体の像が形成される撮像素子上の位置が第２の画像データｇ（ｘ）の被写体の像が形成される撮像素子上の位置と同じであるような画像データを第１の画像データｇ’（ｘ）から生成したものである。画像データｇ’（ｘ＋ａ）を移動前推測画像データと称する。画像データの解析は、具体的にはパターン画像のドットに対応する画像データ（すなわち画素値が落ち込んでいる部分）のｘ座標を読み取ることにより行う。生成部１２がこの画像データの解析処理を行う。例えば、中央のドットに対応して画像データの値が小さくなっている部分のｘ座標から、移動量ａを推定できる。
【００３２】
撮像素子に感度むらがなければ、第３の画像データｇ’（ｘ＋ａ）と第２の画像データｇ（ｘ）とはよく一致するはずである。しかし、撮像素子に感度むらがある場合、ｇ’（
ｘ＋ａ）とｇ（ｘ）とで移動前後で感度むらの影響が出る位置が、第１のずれ量（移動量と移動方向）に応じてずれる。そのため、ｇ’（ｘ＋ａ）とｇ（ｘ）との間には感度むらに基づくずれが生じる。ｇ’（ｘ＋ａ）は、簡単な数学的変形、
【数３】
2015119344.tif 000005
によって、撮像素子上に形成された画像表示装置に表示された画像の像における光の強度分布ｇdisplay（ｘ）と、第１のずれ量に応じて移動した感度分布ｓ（ｘ＋ａ）と、の積
と等価であることがわかる。
【００３３】
図５（ｉ）は、図５（ｃ）の第２の画像データｇ（ｘ）と図５（ｈ）の第３の画像データｇ’（ｘ＋ａ）とを同一座標軸上に描いて比較しやすくした図である。
これら２つの画像データｇ（ｘ）及びｇ’（ｘ＋ａ）を、画面に表示されたパターン画像に対応する光の像の強度分布と、撮像素子の感度分布と、に分解して考える。式１と式３から、これら２つの画像データの間で、パターン画像に対応する光の像の強度分布ｇdisplay（ｘ）は共通で、撮像素子の感度分布だけがｓ（ｘ）とｓ（ｘ＋ａ）とで異なって
いることがわかる。
【００３４】
図５（ｉ）では、パターン画像に対応する光の像の強度分布ｇdisplay（ｘ）に由来す
る、表示面のパターン画像の白色部分の画素値やパターン画像のドットに対応する部分の位置、落ち込み量等は、２つの画像データｇ（ｘ）とｇ’（ｘ＋ａ）とで一致する。このことは、表示装置の画面における輝度むらの有無によらない。画面に輝度むらが有る場合には、輝度むらも含めて、両者は一致する。
【００３５】
しかし、撮像素子の感度分布に由来する部分は、２つの画像データｇ（ｘ）とｇ’（ｘ＋ａ）とでずれが生じるため、重ならなくなる。具体的には、図５（ｃ）の第２の画像データｇ（ｘ）では、撮像素子の感度分布の影響がｓ（ｘ）で画像データに現れるため、画素値が高くなる位置はｘ０である。一方、図５（ｈ）の第３の画像データｇ’（ｘ＋ａ）では、撮像素子の感度分布の影響はｘ軸方向にａだけシフトしたｓ（ｘ＋ａ）で画像データに表れるため、画素値が高くなる位置はｘ０−ａとなる。
【００３６】
このように、表示装置に表示された同一のパターン画像を、撮像装置だけを移動させて異なる位置から撮像した場合、表示装置の画面の輝度むらと、撮像素子の感度むら（感度分布）とは、撮像により得られる画像データにおいて異なる現れ方をする。すなわち、画面の輝度むらの影響は撮像装置の移動前後で現れる位置が変化するが、撮像素子の感度むらの影響は撮像装置の移動前後で現れる位置が変化しない。この違いを利用することによって、表示装置の画面の輝度むらの情報と撮像素子の感度分布の情報とを分離することができる。撮像素子の感度分布の情報のみを抽出することができれば、撮像素子の感度分布に基づいて撮像データにおける感度分布の影響を除去する補正をすることができる。従って、デジタルカメラやカメラ付スマートフォン等を測定器として用いて、精度良くキャリブレーションを行うことができ、均一性の保証された基準面光源が必要なくなる。
【００３７】
次に撮像素子の感度分布の推定方法を図６を用いて説明する。
図６（ｉ）は、図５（ｉ）に示した、２つの画像データｇ（ｘ）とｇ’（ｘ＋ａ）とを重ねあわせて比較した図である。図６（ｊ）は、ｇ’（ｘ＋ａ）をｇ（ｘ）で除算した値ｇ’（ｘ＋ａ）／ｇ（ｘ）を示す図である。２つの画像データの画素値同士の比ｇ’（ｘ＋ａ）／ｇ（ｘ）は、簡単な数学的変形、
【数４】
2015119344.tif 000006
によって、撮像素子におけるｘ軸方向にａ（第１のずれ量）だけ離れた画素間の感度比（第１の感度比）と等価であることがわかる。感度分布算出部１３が、第３の画像データｇ’（ｘ＋ａ）と第２の画像データｇ（ｘ）との比ｇ’（ｘ＋ａ）／ｇ（ｘ）に基づき、第１の感度比ｓ（ｘ＋ａ）／ｓ（ｘ）を算出する。
【００３８】
図６（ｋ）は、ｇ’（ｘ＋ａ）／ｇ（ｘ）の値を用いて求められた撮像素子の感度分布である。感度分布の求め方は、まずある位置ｘ１の画素の感度ｓ（ｘ１）を基準値として定め、そこからａずつ離れた画素の感度を順次求める。具体的には、位置ｘ１の画素の感度ｓ（ｘ１）を基準値Ｃとし、ｓ（ｘ１＋ａ）、ｓ（ｘ１＋２ａ）・・・を順次求める。ここでｇ’（ｘ＋ａ）／ｇ（ｘ）＝ｆ（ｘ）と書くと、上記の式４の関係から、
【数５】
2015119344.tif 000007
の関係が成り立つので、
【数６】
2015119344.tif 000008
【００３９】
このように、第１の感度比を順次積算して撮像素子の第１のずれ量ａずつ離れた画素の感度を算出することにより、感度分布ｓ（ｘ）の離散的な値を求めることができる（第１の離散的な感度分布）。求められた第１の離散的な感度分布に対し補間処理を行うことにより、撮像素子の全画素に対する感度分布ｓ（ｘ）を求めることができる。感度分布算出部１３が、第１の感度比に基づく第１の離散的な感度分布の算出及び補間処理による全画素の感度分布の算出を行う。
【００４０】
図７に、図３〜図６で説明した、撮像素子の感度分布を求める方法を示すフローチャートである。各ステップに併記した（１）〜（８）の符号は、図３〜図６中の（１）〜（８）の符号に対応している。
【００４１】
ステップＳ１では、表示装置（モニタ１７）の画面に測定用パターン画像が表示される。測定用パターン画像は、図３（ａ）に示したような白地背景に黒のドットが配置された画像とする。この処理は、ユーザが手動でモニタ１７やモニタ１７に画像を出力する画像出力装置（不図示）を操作して実行しても良い。或いは、感度分布測定装置１０又はカメラ１５が、モニタ１７や画像出力装置に指示することで実行されても良い。
【００４２】
ステップＳ２では、撮像装置（カメラ１５）によって、表示装置の画面に表示された撮像素子の感度分布測定用パターン画像の撮像が行われる。これにより図３（ｂ）、図３（ｃ）で示したような第２の画像データｇ（ｘ）が得られる。この処理は、ユーザが手動でカメラ１５を操作して実行しても良いし、感度分布測定装置１０がカメラ１５に指示することで実行されても良い。
【００４３】
ステップＳ３では、撮像装置が横移動させられた状態で、再度、表示装置の画面に表示された同じ感度分布測定用パターン画像の撮像が行われる。これにより図４（ｆ）、図４（ｇ）で示したような第１の画像データｇ’（ｘ）が得られる。この処理は、ユーザが手動でカメラ１５を操作して実行しても良いし、感度分布測定装置１０がカメラ１５に指示することで実行されても良い。
【００４４】
ステップＳ４では、感度分布測定装置１０の取得部１１は、第１の画像データｇ’（ｘ）と第２の画像データｇ（ｘ）を取得し、生成部１２がこれらの画像データを解析する。そして、図５（ｇ）に示したように、パターン画像の像が形成される撮像素子上の位置の移動量ａ（第１のずれ量）を算出する。
【００４５】
ステップＳ５では、感度分布測定装置１０の生成部１２は、図５（ｈ）に示したように、算出した移動量ａと、第１の画像データｇ’（ｘ）とから、第１のずれ量ａに基づく変換処理により、第３の画像データｇ’（ｘ＋ａ）（移動前推測画像）を生成する。更に、図６（ｊ）に示したように、感度分布測定装置１０の感度分布算出部１３は、第２の画像データｇ（ｘ）と第３の画像データｇ’（ｘ＋ａ）との比ｇ’（ｘ＋ａ）／ｇ（ｘ）に基づき、撮像素子のａだけ離れた画素間の感度比ｓ（ｘ＋ａ）／ｓ（ｘ）を求める。
【００４６】
ステップＳ６では、感度分布測定装置１０の感度分布算出部１３は、図６（ｋ）に示したように、第１の感度比ｓ（ｘ＋ａ）／ｓ（ｘ）を順次積算して、撮像素子の第１のずれ量ａずつ離れた画素の第１の離散的な感度分布を算出する。更に、感度分布算出部１３は、第１の離散的な感度分布に対し補間処理を行うことにより撮像素子の全画素の感度分布ｓ（ｘ）を算出する。
【００４７】
以上のように、撮像素子の感度分布を測定（推定）することにより、この撮像素子を備える撮像装置（カメラ１５）をキャリブレーション用の測定器として用いて、表示装置（モニタ１７）のキャリブレーションを行うことができる。撮像装置を用いた表示装置のキャリブレーションについて説明する。
【００４８】
図２１は、本実施例に係るキャリブレーション装置の機能構成を示す図である。本実施例のキャリブレーション装置１８は、感度分布取得部２０、画像データ取得部２１、補正部２２、及びキャリブレーション部２３を有する。
【００４９】
感度分布取得部２０は、上述した感度分布測定装置１０により測定された撮像素子の感度分布の情報を取得する。図２１では、感度分布取得部２０は感度分布測定装置１０から直接、撮像素子の感度分布の情報を取得するように描いているが、感度分布取得部２０が感度分布の情報を取得する方法はこれに限らない。図１９（Ｂ）に示すように、メモリカード１６のようなデータ格納メディアから読み出して取得しても良いし、図２０（Ａ）、
図２０（Ｂ）に示すように、有線又は無線の通信手段により取得しても良い。また、図２２（Ａ）、図２２（Ｂ）に示すように、カメラ１５が感度分布測定装置１０を備えている場合は、感度分布取得部２０はカメラ１５からデータ格納メディアや通信手段により感度分布の情報を取得する構成とすることができる。また、図２３（Ａ）に示すように、カメラ或いはカメラ機能を有するキャリブレーションシステムが、感度分布測定装置１０とキャリブレーション装置１８を備えている場合は、感度分布取得部２０は内部的に感度分布の情報を取得する構成とすることができる。
【００５０】
画像データ取得部２１は、モニタ１７に表示されたキャリブレーション用のパターン画像を、撮像素子１４を備えたカメラ１５により撮像して得られた画像データを取得する。図２１では、画像データ取得部２１はカメラ１５から直接、画像データを取得するように描いているが、画像データ取得部２１が画像データを取得する方法はこれに限らない。上記の感度分布取得部２０が感度分布の情報を取得する場合と同様、画像データ取得部２１は、メモリカードのようなメディア、有線又は無線の通信手段、内部的な処理により画像データを取得することができる。また、図２２（Ａ）、図２２（Ｂ）に示すように、カメラ１５が感度分布測定装置１０と撮像素子１４を備えている場合、カメラ１５は、感度分布の情報と画像データとを共通の情報伝達経路で取得する構成とすることもできる。図２２（Ｂ）の例では、送信部１９は、感度分布測定装置１０から出力される感度分布の情報と、撮像素子１４から出力される画像データと、をキャリブレーション装置１８へ出力する。
【００５１】
補正部２２は、感度分布取得部２０が取得した撮像素子１４の感度分布の情報に基づき、画像データ取得部２１が取得した画像データから撮像素子１４の感度分布の影響を除去する補正を行う。
【００５２】
キャリブレーション部２３は、補正部２２により補正された画像データに基づき、モニタ１７のキャリブレーションを行う。
【００５３】
なお、キャリブレーション装置１８は、上記の各機能を実現するプログラム（ソフトウェア）をコンピュータが実行することにより実現することができる。或いは、上記の各機能をコンピュータに追加する機能拡張ボードのようなハードウェアにより実現することができる。キャリブレーション装置１８は、図２３（Ａ）に示すように、カメラ１５と一体のキャリブレーションシステムとして構成することもできる。この場合、キャリブレーションシステムは、カメラ付スマートフォンやノート型ＰＣに上記の各機能を実行するアプリケーションをインストールすることにより実現することができる。また、図２３（Ｂ）に示すようにモニタ１７がキャリブレーション装置１８を備える構成とすることもできる。
【００５４】
キャリブレーションの具体的な方法について説明する。
図８に、表示装置の画面にキャリブレーション用のパターン画像が表示され、撮像装置によってパターン画像が撮像されている様子を示す。この例では、パターン画像は、全画面で最大輝度の白色で構成される。図８（ｌ）は表示装置及び撮像装置の配置図である。図８（ｍ）は撮像素子上に結像投影された像を示す。図８（ｎ）は撮像素子によって取得される画像データのうち、図８（ｍ）のｘ軸上の画像データｇ（ｘ）を示す。図８（ｏ）は撮像素子のｘ軸上の感度分布ｓ（ｘ）を示す。図８（ｐ）は表示装置の画面のＸ軸上の輝度分布Ｂ（Ｘ）を示す。
【００５５】
図８（ｌ）に示すように、撮像装置は表示装置の画面の中央と正対して配置されている。図８（ｍ）に示すように、撮像素子上には全白パターンを表示している画面全体の像が結像されている。このとき、ｘ軸上の画像データは、図８（ｎ）に示すように、表示装置
の画面のＸ軸上の輝度分布Ｂ（Ｘ）と撮像素子のｘ軸上の感度分布ｓ（ｘ）とを乗じた画像データとなっている。ここでは、表示装置の画面の輝度分布が均一である場合を例に説明する。この場合、図８（ｎ）の画像データは、図８（ｏ）の撮像素子の感度分布に比例した形となる。
【００５６】
図９に、表示装置の画面に表示された全白パターンを撮像した画像データから、画面の輝度むらを推定する方法を示す。図９（ｎ）は、図８（ｎ）に示した、画面に全白パターンが表示されているときに撮像素子によって取得されたｘ軸上の画像データｇ（ｘ）を示す。図９（ｋ）は、図６（ｋ）で示した、本実施例の感度分布測定装置１０によって測定された撮像素子のｘ軸上の感度分布ｓ（ｘ）を示す。図９（ｑ）は、図９（ｎ）の画像データｇ（ｘ）を、図９（ｋ）の撮像素子のｘ軸上の感度分布ｓ（ｘ）で除算することで補正された画像データを示す。
【００５７】
先に述べたように、撮像素子によって取得されたｘ軸上の画像データｇ（ｘ）は、画面の実際の輝度分布ｇdisplay（ｘ）と撮像素子の感度分布ｓ（ｘ）との積と考えられる。
従って、撮像素子によって取得されたｘ軸上の画像データｇ（ｘ）を撮像素子の感度分布ｓ（ｘ）で除算することで、画面の実際の輝度分布ｇdisplay（ｘ）が求められる。ｇdisplay（ｘ）は、表示装置の画面のＸ軸上の実際の輝度分布Ｂ（Ｘ）を表すと考えられるので、これを用いて表示装置のキャリブレーションを行うことができる。
【００５８】
上記の説明では、表示装置の画面の輝度分布Ｂ（Ｘ）は均一としたため、得られたｇdisplay（ｘ）は実際の輝度分布を反映して均一となっている。表示装置の画面の輝度分布
Ｂ（Ｘ）が均一でなくむらがある場合は、上記の計算で得られるｇdisplay（ｘ）は、実
際の画面に生じている輝度むらが反映された分布となる。
【００５９】
従って、表示装置の画面の輝度分布に、バックライト光源の劣化等による輝度むらが生じている場合には、上記のように求められたｇdisplay（ｘ）に基づいて、その輝度むら
を精度良く推定することができる。キャリブレーション装置１８は、推定した表示装置の画面の輝度分布に生じている輝度むらに基づき、入力画像信号を補正するための情報や、バックライト光源の制御信号を補正するための情報を、キャリブレーション結果として出力する。キャリブレーション装置１８は、キャリブレーション結果を、モニタ１７や、モニタ１７に画像信号を出力する画像出力装置（不図示）へ出力する。
【００６０】
以上説明した方法では、感度分布測定装置は、撮像素子の離散的な感度分布ｓ（ｘ）を求め、それに対し補間処理を行って撮像素子の全画素の感度分布ｓ（ｘ）を求めた。しかし、離散データの間隔が粗いと、補間処理による誤差が大きくなる。この離散データの間隔を細かくするためには、移動量ａ（被写体の像が形成される撮像素子上の位置のずれ量）の値を小さくすれば良い。すなわち、撮像装置の移動量Ａの大きさを小さくすれば良い。しかし、移動量ａの値が小さくなるほど、離散的な感度分布を求める際の感度比を順次積算する計算のステップ数が多くなり、補間処理による誤差は小さくなるものの、積算による誤差が増えてしまう。
【００６１】
これに対する改善方法の一例を図１０に示す。図１０（ｋ）は、図６（ｋ）で示した方法により算出した撮像素子の離散的な感度分布（第１の離散的な感度分布）を示している。図１０（ｒ）は、より高精度に撮像素子の感度分布を算出する方法を示す。
【００６２】
図１０（ｒ）に示す方法では、移動量ａ（被写体の像が形成される撮像素子上の位置のずれ量）を段階的に変えて、すなわち撮像装置の移動量Ａを段階的に変えて撮像された複数の画像データを取得する。そして、移動量（撮像条件の差異）が大きい画像データに基づき求められた感度比の順次計算と、移動量が小さい画像データに基づき求められた感度
比の順次計算とを組み合わせる。これにより、順次積算の回数の増大を抑制しつつ、離散間隔を縮小して補間処理による誤差も低減し、以て精度を改善する。
【００６３】
例えば、取得部１１は、移動量ＡでＸ軸方向に撮像装置を移動して撮像して得られた画像データｇ１（ｘ）（第１の画像データ）と、移動量ＢでＸ軸方向に撮像装置を移動して撮像して得られた画像データｇ２（ｘ）（第４の画像データ）と、を取得する。移動量Ａに対応する撮像素子上の像形成位置の間隔をａ（第１のずれ量）、移動量Ｂに対応する撮像素子上の像形成位置の間隔をｂ（第２のずれ量）とする。ここでは、Ａ＞Ｂ、ａ＞ｂとする。特に、図１０（ｒ）に示すように、ａ＝４ｂとするが、この関係は説明のための一例であり、この例に限られない。
【００６４】
生成部１２は、ｇ１（ｘ）とｇ２（ｘ）のそれぞれについて画像解析等により撮像装置の移動量Ａ，Ｂに対応する撮像素子上の像形成位置の間隔ａ、ｂ（第１のずれ量、第２のずれ量）を求める。生成部１２は、第１の画像データｇ１（ｘ）、第４の画像データｇ２（ｘ）に対し、これら第１のずれ量、第２のずれ量に基づく変換処理を行うことにより、第３の画像データｇ１（ｘ＋ａ）、第５の画像データｇ２（ｘ＋ｂ）を生成する。そして、式４と同様に、移動前の撮像により得られた第２の画像データｇ（ｘ）との比を算出し、撮像素子上で第１のずれ量ａだけ離れた画素同士の第１の感度比ｆ１（ｘ）、第２のずれ量ｂだけ離れた画素同士の第２の感度比ｆ２（ｘ）を求める。
【数７】
2015119344.tif 000009
上記の式より、
【数８】
2015119344.tif 000010
の関係が成り立つので、例えば位置ｘ１＋ａ＋２ｂにおける感度は、
【数９】
2015119344.tif 000011
のように求められ、感度比の積算は３回で済む。しかし、離散間隔の細かい移動量Ｂによる第４の画像データｇ２（ｘ）しか取得していなければ、
【数１０】
2015119344.tif 000012
のように、第２の感度比ｆ２（ｘ）を６回積算をする必要がある。
【００６５】
上記説明では、第１の画像データｇ’（ｘ）、ｇ１（ｘ）、第４の画像データｇ２（ｘ）は、撮像装置を横方向（ｘ軸方向）に移動させて第２の画像データｇ（ｘ）と同じ被写体を撮像された画像データである例を示したが、これは一例にすぎない。撮像条件を異ならせるために撮像装置を移動させる方向はｘ方向に限られるわけではない。縦方向であるｙ方向や斜め方向であってもよい。いずれの方向の移動であっても、生成部１２は、撮像装置の移動後の撮像により取得された画像データを解析して、撮像素子上での像形成位置のずれ量（移動方向と移動量）を求めることができる。さらに、生成部１２は、ずれ量に基づき、移動前に撮像された画像データと同じ撮像条件（撮像素子上での像形成位置が同じ）である画像データ（移動前推測画像データ）を生成することができる。感度分布算出部１３は、移動前推測画像データを、撮像装置の移動前の撮像により取得された画像データと比較することにより、撮像素子上でずれ量（移動量と移動方向）に対応する位置関係にある画素間の撮像素子の感度比を算出する。感度分布算出部１３は、この感度比に基づき、撮像素子の感度分布ｓ（ｘ）を算出する。
【００６６】
さらに一般には、第１の画像データｇ’（ｘ）、ｇ１（ｘ）、第４の画像データｇ２（ｘ）は、撮像装置を、ｘ軸方向やｙ軸方向への並行移動だけではなく、回転や拡大縮小を含む一般的な移動を行って撮像された画像データであっても良い。この場合第１のずれ量は、並行移動、回転、及び拡大縮小の組み合わせで表される量であり、生成部１２は、移動後の撮像装置による撮像により得られた画像データに対し、前記組み合わせに対応するアフィン変換処理を行うことで移動前推測画像データを生成する。生成部１２は、撮像装置の移動後の撮像により取得された画像データを解析して移動に対応するアフィン変換の変換行列を求め、更にその逆変換行列を求めて、移動後の撮像により取得された画像データに適用することにより、移動前推測画像データを求める。感度分布算出部１３は、移動前推測画像データを、撮像装置の移動前の撮像により取得された画像データと比較することにより、撮像素子の感度分布ｓ（ｘ）を算出する。
【００６７】
以上の説明は、先に述べたように、簡単のため、光学系により撮像素子上に結像された像の光の強度分布における、撮像装置の光学系固有の光学特性の影響を考慮していなかった。以下、光学系固有の光学特性の影響を考慮する場合について説明する。
【００６８】
図１１に、表示装置の画面に撮像素子感度分布測定用のパターン画像（この例では、白地背景に黒色のドットが配置された画像）が表示され、撮像装置によって撮像している様子を示す。これは先に述べた図３に対応する図である。図３においては光学系固有の光学特性の影響を省略していたが、図１１においては光学系固有の光学特性の影響を図１１（ｄ’）に示している。図１１（ａ）〜図１１（ｄ）は図３（ａ）〜図３（ｄ）に対応する。図１１（ｄ’）は、光学系固有の光学特性として、光学系に起因する光量低下（口径蝕等による周辺光量の低下）の特性ｏ（ｘ）を示す。
【００６９】
撮像素子によって取得される画像データには、この光学系固有の光学特性ｏ（ｘ）の影響と、撮像素子の感度分布ｓ（ｘ）の影響とが、重畳して現れる。すなわち、
【数１１】
2015119344.tif 000013
【００７０】
このため、図１１（ｃ）に示す撮像によって得られた画像データにおけるｘ軸上の画像データｇ（ｘ）は、周辺部が低下した形状となっている。
【００７１】
次に、表示装置はそのままで、撮像装置だけをｘ軸方向に−Ａだけ移動させて同じ測定
用パターン画像表示を撮像した場合を考える。なお、図示が煩雑になるのを避けるため、図４、図５に対応する図は省略する。
【００７２】
表示装置の画面のＸ軸上の輝度分布Ｂ（Ｘ）をそのまま反映する強度分布ｇdisplay（
ｘ）とする。この場合、図４、図５等で説明したように、移動後の撮像により得られる画像データにおいて、輝度分布Ｂ（Ｘ）の影響は、全体にａだけ移動してｇdisplay（ｘ−
ａ）の形で現れる。
【００７３】
しかし、撮像装置を移動させても、光学系と撮像素子との相対的な位置関係は変わらないため、移動後の撮像により得られる画像データにおいて、光学系固有の光学特性の影響は、ｏ（ｘ）の形で現れ、移動前と変わらない。以上のことから、撮像素子上に結像投影される像は、ｇdisplay（ｘ−ａ）とｏ（ｘ）との積ｇdisplay（ｘ−ａ）・ｏ（ｘ）となる。また、図４，図５等で説明したように、撮像装置を移動させても、移動後の撮像により得られる画像データにおいて、撮像素子の感度分布の影響は、ｓ（ｘ）の形で現れ、移動前と変わらない。従って、撮像装置の移動後に撮像素子によって取得されるｘ軸上の画像データｇ’（ｘ）は、画面の光の強度分布ｇdisplay（ｘ−ａ）と、光学系固有の光学
特性ｏ（ｘ）と、撮像素子のｘ軸上の感度分布ｓ（ｘ）と、の積となる。
【数１２】
2015119344.tif 000014
【００７４】
このように、表示装置に表示された同一の画像を、撮像装置だけを移動させて異なる位置から撮像した場合、光学系固有の光学特性の影響と撮像素子の感度分布の影響とは、撮像により得られる画像データにおいて、移動前と同様の現れ方をする。一方、表示装置の表示輝度分布は、撮像装置の移動後の撮像により得られる画像データにおける現れ方が移動前と異なる。この違いを利用することによって、表示装置の画面の輝度分布の情報と、感度分布と光学特性の積ｏ（ｘ）・ｓ（ｘ）の情報と、は分離することができる。以下、感度分布と光学特性の積ｏ（ｘ）・ｓ（ｘ）を撮像特性ｐ（ｘ）と称する。
【００７５】
この考えに基づき、前述の図５、図６に示した方法によって、移動量ａだけ離れた撮像素子の画素間での撮像特性ｐ（ｘ）の比ｐ（ｘ＋ａ）／ｐ（ｘ）を求めることができる。この撮像特性比を図６（ｋ）や図１０（ｒ）と同様に順次積算し、更に適切な補間処理を行うことにより、撮像素子の全画素について、撮像特性ｐ（ｘ）の値が求められる。
【００７６】
図１２に、表示装置の画面に表示された全白パターンを撮像して得られた画像データから、画面の輝度分布を推定する方法を示す。図１２（ｎ）は、画面に全白パターンが表示されているときに撮像素子によって取得されたｘ軸上の画像データｇ（ｘ）である。図１２（ｋ）は、上述の方法によって計算された撮像特性ｐ（ｘ）（光学系固有の光学特性ｏ（ｘ）と撮像素子のｘ軸上の感度分布ｓ（ｘ）との積ｏ（ｘ）・ｓ（ｘ））を示す。
【００７７】
図１２（ｑ）は、図１２（ｎ）の画像データｇ（ｘ）を、図１２（ｋ）の撮像特性ｐ（ｘ）で除算することで補正された画像データである。先に述べたように、画像データｇ（ｘ）は、表示装置の画面の輝度分布ｇdisplay（ｘ）と、撮像特性ｐ（ｘ）と、の積と考
えられる。従って、この除算の結果、画像データから撮像装置の撮像特性の影響を除去して、表示装置の画面の輝度分布ｇdisplay（ｘ）が求められる。すなわち、光学系固有の
光学特性を考慮した場合でも、光学系固有の光学特性を考慮しない場合と同様の計算により、表示装置の画面の輝度分布を、撮像装置による撮像データから推定することができる。従って、光学系固有の光学特性を考慮した場合でも、本実施例の方法により、表示装置のキャリブレーションに必要な、画面の輝度分布ｇdisplay（ｘ）が得られることがわか
る。
【００７８】
このように、撮像装置の光学系固有の光学特性の情報が未知でも、表示装置のキャリブレーションに必要な画面の輝度分布情報を得ることができる。従って、撮像装置をキャリブレーション用の測定器として用いて表示装置のキャリブレーションを行うことができる。なお、撮像素子の感度分布ｓ（ｘ）のみの情報が必要な場合は、光学系固有の光学特性ｏ（ｘ）の情報に基づき、上述の方法で得られた撮像特性ｐ（ｘ）をｏ（ｘ）で除算してｓ（ｘ）を推定することができる。
【００７９】
（実施例２）
本発明に係る、撮像素子の感度分布を測定する測定装置の実施例２を説明する。実施例２では、撮像装置の光学系に対する撮像素子の相対位置を変化させることにより、撮像素子上に投影する光分布に変化を与える。撮像装置の光学系は、キャリブレーション用のパターン画像を表示している表示装置の画面を撮像素子上に結像投影する光学系であり、被写体の像を撮像素子上に形成する。この光学系は、合焦距離を変更可能である。光学系に対する撮像素子の相対位置を変化させることは、合焦距離を変化させることである。この撮像条件の変化に伴う撮像素子出力の変化を比較することにより、撮像素子の感度分布を推定する。本実施例では、被写体は、撮像素子の感度分布測定用のパターン画像を表示している表示装置の画面である。
【００８０】
本実施例の感度分布測定装置１０の構成は実施例１と同様である。取得部１１が取得する第１の画像データ及び第２の画像データ、生成部１２による第３の画像データの生成方法、感度分布算出部１３による撮像素子１４の感度分布の算出方法が実施例１と異なる。以下、相違点を中心に説明する。
【００８１】
図１３に、表示装置の画面に撮像素子の感度分布測定用のパターン画像（ここでは、黒地背景に白色のドットが配置された画像）が表示された状態で、当該パターン画像を撮像装置によって撮像している様子を示す。このとき撮像装置のピントは表示装置の画面に合っているとする。すなわち、光学系と被写体である画面との距離と合焦距離とが等しい撮像条件で撮像が行われる。図１３（ａ）は表示装置及び撮像装置の配置図である。図１３（ｂ）は、撮像素子上に結像投影された像である。図１３（ｃ）は、撮像素子によって取得される画像データのうち、上記図１３（ｂ）の中央のドット近傍の画像データｇ（ｘ）（第１の画像データ）である。図１３（ｄ）は、撮像素子のｘ軸上の感度分布ｓ（ｘ）を示す。
【００８２】
図１３（ａ）に示すように、撮像装置が表示装置の画面の中央と正対して配置されている。そのため、図１３（ｂ）に示すように、撮像素子上には、画面中央に表示されているドットの像が撮像素子上の中央の位置に形成されるように、画面全体（パターン画像）の像が結像されている。中央ドット近傍の画像データにおいて、撮像装置のピントは表示装置の画面に合わされている。そのため、図１３（ｃ）に示すように、中央ドット（明部である白点）に対応する部分の画像データ値が鋭く高くなっており、それ以外の部分（暗部である黒背景）のデータ値が低く、ほぼ０になっている。
【００８３】
図１４に、表示装置はそのままで、撮像装置の光学系だけが撮像素子から僅かに離れる方向（又は近づく方向）に動かされた状態で、同一の被写体（測定用パターン画像を表示している画像表示装置の画面）を撮像している様子を示す。すなわち図１４では、光学系と被写体との距離と合焦距離とが等しくない撮像条件でピントをぼかして撮像が行われている。図１４（ｅ）は表示装置及び撮像装置の配置図である。図１４（ｆ）は撮像素子上に結像投影された像を示す。図１４（ｇ）は撮像素子によって取得される画像データのうち、図１４（ｇ）の中央ドット近傍の画像データｇ’（ｘ）を示す（第２の画像データ）。図１４（ｄ）は撮像素子のｘ軸上の感度分布ｓ（ｘ）を示す。
【００８４】
図１４（ｆ）の撮像素子上に結像投影された像は、ピントが合っていないために図１３（ｂ）よりもぼけており、第２の画像データｇ’（ｘ）のパターン画像のドットに対応する部分は、図１３（ｃ）と比較して広がって大きくなっている。
【００８５】
図１４（ｇ）の画像データにおいては、パターン画像のドットに対応する部分の画素値は図１３（ｃ）の場合と比べて低くなり、また画素値の分布は周辺に裾を引いて広がった形状になっている。図１４（ｇ）において、実線は、第１の画像データとは異なる撮像条件で実際に撮像して得られた第２の画像データｇ’（ｘ）を示す。破線は、ピントが合った状態の第１の画像データと、光学系（レンズ）の特性と、ピントぼかし量（合焦距離の差）と、に基づき生成部１２が生成した、ピントをずらした状態で得られる画像データの予想データｇ（ｘ０）・ｌ（ｘ−ｘ０）である。この予想データは、合焦距離が第２の画像データの合焦距離と同じである画像データを第１の画像データの画素値を用いて生成した第３の画像データである。
【００８６】
ここで、ｇ（ｘ０）は、図１３（ｃ）で示した、ピントが合っている状態で取得された第１の画像データの位置ｘ＝ｘ０におけるドットに対応するピーク部分における値である。ここでは、パターン画像の中央のドットを例に説明しているので、ピーク値はｘ＝０のところの値である。ｌ（ｘ−ｘ０）は、ピントをずらした時の像のボケ方を示す点像分布関数で、光学系の固有の特性と合焦距離の差によって決まる。点像分布関数は、ピントが合っている状態のピーク値で規格化されている。生成部１２は、この点像分布関数に、ピントが合っている状態で実際に撮像して得られた第１の画像データのピーク値ｇ（ｘ０）を乗じることによりスケーリングを行って、第３の画像データ（予想データ）を生成する。
【００８７】
生成されたピントを外した画像データの予想データｇ（ｘ０）・ｌ（ｘ−ｘ０）と、実際にピントを外して撮像により得られた画像データｇ’（ｘ）とは、撮像素子の感度分布が均一であれば、一致するはずである。
【００８８】
しかし、この例では、図１４（ｄ）に示すように、撮像素子の感度が、ｘ軸上の右の方（ｘ＝ｘｃの近傍より大きい方）で高くなっている。これを反映して、図１４（ｇ）に示すように、画像データの裾を引いて広がっている部分のｘ＝ｘｃ近傍において、実線のデータ（ぼけた像を撮像した第２の画像データ）が破線のデータ（点像分布関数から予想した第３の画像データ）よりやや高くなっている。また、撮像素子の画素毎にも細かい感度のばらつきがある。そのため、予想データｇ（ｘ０）・ｌ（ｘ−ｘ０）は滑らかな曲線になっているのに対し、実際に得られる画像データｇ’（ｘ）はばらつきの影響で細かな凹凸を有している。
【００８９】
このように、ピントをぼかした状態で撮像して得られた第２の画像データと、第１の画像データ及び点像分布関数から予想した第３の画像データとで、裾を引いて広がっている部分を比較することで、撮像素子の感度分布の情報が得られる。
【００９０】
本実施例においては、撮像素子上に結像するパターン画像のドットの像が、光学系のピントをずらしたときにボケる範囲内で、撮像素子の感度分布の情報が得られる。すなわち、上記の方法で得られるのは、撮像素子の局所的な感度分布の情報である。撮像素子の全域の感度分布の情報を得るために、パターン画像に複数のドットを高密度に配置することが考えられる。あるいは以下に示すように、ドットの配置パターンが異なる複数のパターン画像について、ピントが合っている状態とずれている状態との画像データを取得する方法がある。
【００９１】
図１５に、ドットの配置パターンが図１３，図１４の場合と異なるパターン画像に切り替えてピントが合っている状態とずれている状態との画像データを取得する様子を示す。図１５（ｈ）は表示装置及び撮像装置の配置図を示す。図１５（ｉ）は、ピントが合っている状態で撮像素子によって取得される画像データのうち、画面中央近傍の画像データｇ（ｘ）を示す（第１の画像データ）。図１５（ｊ）は、ピントがずれている状態で撮像素子によって取得される画像データのうち、画面中央近傍の画像データｇ’（ｘ）を示す（第２の画像データ）。図１５（ｄ）は撮像素子のｘ軸上の感度分布ｓ（ｘ）を示す。
【００９２】
図１５（ｈ）に示すように、表示装置の画面に表示されているパターン画像のドットの配置パターンが図１３（ａ）、図１４（ｅ）に示すパターン画像のドット配置パターンから変更されている。図１５（ｈ）に示すパターン画像は、図１３（ａ）、図１４（ｅ）に示すパターン画像のドットとドットの間の位置にドットを配置したものである。図１５（ｈ）に示すパターン画像では、図１３（ａ）、図１４（ｅ）に示すパターン画像のドットの位置にはドットは存在しない。このようなマーカパターンを表示させた状態で上述した計算を行うことで、より広範囲の撮像素子の感度分布の情報が得られるようになる。
【００９３】
図１５（ｉ）は、ピントが合っている状態で撮像素子によって取得される画像データのうち、画面の縦方向中央近傍（ｘ軸近傍）の画像データｇ（ｘ）である（第１の画像データ）。図１５（ｄ）に示すように、撮像素子の感度が、ｘ軸上の右の方（ｘ＝ｘｃ近傍）で高くなっている。これを反映して、第１の画像データｇ（ｘ）において、右側のドット（ｘ＝ｘ＋１）に対応する部分のピーク値ｇ（ｘ＋１）が左側のドット（ｘ＝ｘ−１）に対応する部分のピーク値ｇ（ｘ−１）より高くなっている。
【００９４】
図１５（ｊ）の、ピントがずれている状態で撮像素子によって取得される画像データｇ’（ｘ）（第２の画像データ）においても、右側のドットに対応する部分のピーク値が左側のドットに対応する部分のピーク値より高くなっている。左右の２つのドットの各々に対応する予想データｇ（ｘ−１）・ｌ（ｘ−ｘ−１）、ｇ（ｘ＋１）・ｌ（ｘ−ｘ＋１）も、同様に右側のドットに対応する部分のピーク値が左側のドットに対応する部分のピーク値より高くなっている。これは、ピントが合っている状態で得られた画像データ（第１の画像データ）ｇ（ｘ）における、右側のドットに対応する部分のピーク値ｇ（ｘ＋１）が、左側のドットに対応する部分のピーク値ｇ（ｘ−１）より大きいためである。
【００９５】
図１６に基づき、上記のようにして得られた画像データから撮像素子の感度分布を算出する方法を説明する。
図１６（ｇ）は、図１４（ｇ）と同様、画面に表示された図１４（ｅ）のパターン画像をピントをずらして実際に撮像して得られた第２の画像データと、ピントを合わせて撮像して得られた第１の画像データから生成した予想データ（第３の画像データ）と、を示す。図１６（ｊ）は、図１５（ｊ）と同様、画面に表示された図１５（ｈ）のパターン画像をピントをずらして実際に撮像して得られた第２の画像データと、ピントを合わせて撮像して得られた第１の画像データから生成した予想データ（第３の画像データ）と、を示す。図１６（ｋ）は、図１６（ｇ）に示す２つの画像データ及び図１６（ｊ）に示す２つの画像データのそれぞれにおいて、実際の撮像により得られた第２の画像データを予想データ（第３の画像データ）で除算して得られた撮像素子の局所感度分布を示す。図１６（ｌ）は、局所感度分布を繋ぎ合わせて得た、撮像素子の全体の感度分布である。
【００９６】
図１６（ｋ）において、予想データｇ（ｘ０）・ｌ（ｘ−ｘ０）と第２の画像データｇ’（ｘ）との差は、撮像素子を構成する画素のうち位置ｘ０のドットの像が結像している位置の画素及びその周囲の画素の局所的な感度分布（むらやばらつき）によるものである。従って、実際の撮像により得られた第２の画像データを予想データ（第３の画像データ）で除算すれば、そのドットの像の位置の画素及びその周囲の画素の撮像素子の局所的な
感度分布が得られる。これを、局所感度分布と称する。予想データは、ｇ（ｘ０）・ｌ（ｘ−ｘ０）、ｇ（ｘ−１）・ｌ（ｘ−ｘ−１）、ｇ（ｘ＋１）・ｌ（ｘ−ｘ＋１）である。ｇ（ｘ０）、ｇ（ｘ−１）、ｇ（ｘ＋１）には、各ドットに対応する位置ｘ＝ｘ０、ｘ＝ｘ−１、ｘ＝ｘ＋１における撮像素子の画素の感度の情報が含まれる。そのため、第２の画像データを予想データで除算して得られる、各ドットに対応する位置の画素の撮像素子の局所感度分布は、各ドットに対応する位置の画素の撮像素子の感度を１とする相対感度になっている。
【００９７】
図１６（ｌ）に示すように、感度分布算出部１３は、図１６（ｋ）で得られた局所感度分布を繋ぎ合わせて、撮像素子の全体感度分布を得る。繋ぎ合わせる際に、それぞれの局所感度が相対感度であることを考慮して、適宜スケーリングしながら繋ぎ合わせていく必要がある。本実施例では、第１の画像データ及び第２の画像データは、隣接する２つのドットに対応する局所感度分布の各々に、同じ画素に対応する感度の情報が含まれるように、ドット配置が設定されたパターン画像を撮像して得られた画像データであるとする。すなわち、複数のパターン画像は、隣り合う２つのドットに対応する局所感度分布にオーバーラップが生じるように、ドット配置が設定された画像であるとする。感度分布算出部１３は、隣接するドットに対応する局所感度分布同士で、同じ画素に対応する感度の値が同じ値になるよう、局所感度分布にスケーリングを行って繋ぎ合わせていく。
以上のようにして、撮像素子の感度分布ｓ（ｘ）を得ることができる。
【００９８】
図１７及び図１８は、図１３〜図１６で説明した、撮像素子の感度分布を求める方法を示すフローチャートである。各ステップの横に併記した（１）〜（９）の符号は、図１３〜図１６中の（１）〜（９）の符号に対応している。
【００９９】
ステップＳ１１では、表示装置（モニタ１７）の画面に測定用パターン画像が表示される。測定用パターン画像は、図１３（ａ）に示すような、黒地背景に白のドットの配置された画像とする。この処理は、ユーザが手動でモニタ１７やモニタ１７に画像を出力する画像出力装置（不図示）を操作して実行しても良い。或いは、感度分布測定装置１０又はカメラ１５が、モニタ１７や画像出力装置に指示することで実行されても良い。
【０１００】
ステップＳ１２では、撮像装置（カメラ１５）によって、表示装置の画面に表示された撮像素子の感度分布測定用パターン画像の撮像が行われる。この撮像は、被写体であるパターン画像と撮像装置の光学系の合焦距離とが等しい撮像条件で行われ、図１３（ｃ）に示すような第１の画像データｇ（ｘ）が得られる。この処理は、ユーザが手動でカメラ１５を操作して実行しても良いし、感度分布測定装置１０がカメラ１５に指示することで実行されても良い。
【０１０１】
ステップＳ１３では、撮像装置のピントが所定量だけずらされた状態で、再度、同じ感度分布測定用パターン画像の撮像が行われる。これにより図１４（ｇ）に示すような第２の画像データｇ’（ｘ）が得られる。この処理は、ユーザが手動でカメラ１５を操作して実行しても良いし、感度分布測定装置１０がカメラ１５に指示することで実行されても良い。
【０１０２】
ステップＳ１４では、表示装置の画面に測定用パターン画像が表示される。この測定用パターン画像は、図１５（ｈ）に示すように、ステップＳ１１で表示された測定用パターン画像（図１３（ａ）や図１４（ｅ）に示す）とはドットの配置が異なる測定用パターン画像である。この処理は、ユーザが手動でモニタ１７やモニタ１７に画像を出力する画像出力装置（不図示）を操作して実行しても良い。或いは、感度分布測定装置１０又はカメラ１５が、モニタ１７や画像出力装置に指示することで実行されても良い。
【０１０３】
ステップＳ１５では、撮像装置によって、表示装置の画面に表示された撮像素子の感度分布測定用パターン画像の撮像が行われる。この撮像は、被写体であるパターン画像と撮像装置の光学系の合焦距離とが等しい撮像条件で行われ、図１５（ｉ）に示すような第１の画像データｇ（ｘ）が得られる。この処理は、ユーザが手動でカメラ１５を操作して実行しても良いし、感度分布測定装置１０がカメラ１５に指示することで実行されても良い。
【０１０４】
ステップＳ１６では、撮像装置のピントが所定量だけずらされた状態で、再度、同じ感度分布測定用パターン画像の撮像が行われる。これにより図１５（ｊ）に示すような第２の画像データｇ’（ｘ）が得られる。この処理は、ユーザが手動でカメラ１５を操作して実行しても良いし、感度分布測定装置１０がカメラ１５に指示することで実行されても良い。
【０１０５】
ステップＳ１７では、感度分布測定装置１０の生成部１２は、ステップＳ１３、Ｓ１６のピント外し状態に対応する点像分布関数ｌ（ｘ）と、ステップＳ１２、Ｓ１５で実際に撮像により得られた第１の画像データｇ（ｘ）と、から予想データを生成する。予想データは第３の画像データである。ここでは、生成部１２は、図１４（ｇ）に示すように、ｇ（ｘ０）・ｌ（ｘ−ｘ０）、ｇ（ｘ−１）・ｌ（ｘ−ｘ−１）、ｇ（ｘ＋１）・ｌ（ｘ−ｘ＋１）を算出する。
【０１０６】
ステップＳ１８では、感度分布測定装置１０の感度分布算出部１３は、ステップＳ１３、Ｓ１６で得られたピントをずらした状態で実際に撮像して得られた第２の画像データｇ’（ｘ）とステップＳ１７で求めた第３の画像データ（予想データ）との比を計算する。これにより、局所感度分布が算出される。
【０１０７】
ステップＳ１９では、感度分布算出部１３は、ステップＳ１８で求めた局所感度分布を適宜スケーリングして繋ぎ合わせることにより、撮像素子の全体感度分布ｓ（ｘ）を算出する。
【０１０８】
以上のように、撮像素子の感度分布を測定することにより、この撮像素子を備える撮像装置をキャリブレーション用の測定器として用いて、表示装置のキャリブレーションを行うことができる。キャリブレーションの方法は実施例１と同様である。すなわち、モニタ１７に表示されたキャリブレーション用パターン画像を撮像素子１４で撮像して得られた画像データを、算出した撮像素子１４の感度分布を用いて補正し、補正後の画像データを用いてモニタ１７のキャリブレーションを行う。
【０１０９】
本実施例では、実施例１の方法とは異なり、予想データを算出する際に用いたｇ（ｘ０）、ｇ（ｘ−１）、ｇ（ｘ＋１）の中に光学系固有の光学特性の情報も含まれている。そのため、上記の方法によって得られる感度分布ｓ（ｘ）は、光学系固有の光学特性ｏ（ｘ）の情報を含まない。光学系固有の光学特性ｏ（ｘ）の情報を別途取得することにより、光学系固有の光学特性ｏ（ｘ）をも考慮してより高精度に表示装置のキャリブレーションを行うことも可能である。
【０１１０】
以上、本発明の発明を実施するための形態について実施例を用いて説明したが、実施の形態は実施例に限るものではない。記憶装置に記録されたプログラムを読み込み実行することで前述した実施形態の機能を実現するシステムや装置のコンピュータ（又はＣＰＵ、ＭＰＵ等のデバイス）によっても、本発明を実施することができる。また、例えば、記憶装置に記録されたプログラムを読み込み実行することで前述した実施形態の機能を実現するシステムや装置のコンピュータによって実行されるステップからなる方法によっても、本発明を実施することができる。この目的のために、上記プログラムは、例えば、ネット
ワークを通じて、又は、上記記憶装置となり得る様々なタイプの記録媒体（つまり、非一時的にデータを保持するコンピュータ読取可能な記録媒体）から、上記コンピュータに提供される。したがって、上記コンピュータ（ＣＰＵ、ＭＰＵ等のデバイスを含む）、上記方法、上記プログラム（プログラムコード、プログラムプロダクトを含む）、上記プログラムを非一時的に保持するコンピュータ読取可能な記録媒体は、いずれも本発明の範疇に含まれる。
【符号の説明】
【０１１１】
１０：感度分布測定装置、１１：取得部、１２：生成部、１３：感度分布算出部
【図１】
2015119344.tif 000015
【図２】
2015119344.tif 000016
【図３】
2015119344.tif 000017
【図４】
2015119344.tif 000018
【図５】
2015119344.tif 000019
【図６】
2015119344.tif 000020
【図７】
2015119344.tif 000021
【図８】
2015119344.tif 000022
【図９】
2015119344.tif 000023
【図１０】
2015119344.tif 000024
【図１１】
2015119344.tif 000025
【図１２】
2015119344.tif 000026
【図１３】
2015119344.tif 000027
【図１４】
2015119344.tif 000028
【図１５】
2015119344.tif 000029
【図１６】
2015119344.tif 000030
【図１７】
2015119344.tif 000031
【図１８】
2015119344.tif 000032
【図１９】
2015119344.tif 000033
【図２０】
2015119344.tif 000034
【図２１】
2015119344.tif 000035
【図２２】
2015119344.tif 000036
【図２３】
2015119344.tif 000037