
カメラ、顕微鏡、人間の目、プロジェクターなどの光学システムの光伝達関数(OTF)は、そのイメージングコントラストのスケール依存的な記述です。その大きさは、調和強度パターンの画像コントラストです。空間周波数の関数として、一方、その複素引数は周期的なパターンの位相シフトを示します。光伝達関数は、光学技術者が、光学系が物体やシーンからの光を写真フィルム、検出器アレイ、網膜、スクリーン、あるいは単に光伝送チェーンの次の要素にどのように投影するかを記述するために使用します。
正式には、光伝達関数は点像分布関数(PSF、すなわち光学系のインパルス応答、点光源の像)のフーリエ変換として定義されます。フーリエ変換であるOTFは一般に複素数値ですが、中心に対して対称なPSFの場合には実数値となります。実際には、光伝達関数の大きさまたは絶対値で表される画像コントラストが最も重要です。この導出された関数は一般に変調伝達関数(MTF )と呼ばれます。
右側の画像は、パネル (a) と (d) に 2 つの異なる光学系の光伝達関数を示しています。前者は、円形瞳を持つ理想的な回折限界イメージング システムに対応します。その伝達関数は、空間周波数とともにほぼ徐々に減少し、回折限界に達します。この場合、500 サイクル/ミリメートル、つまり周期 2 μm です。このイメージング システムでは、この周期のような小さな周期的特徴が捉えられるため、その解像度は 2 μm であると言えます。[ 1 ]パネル (d) は、焦点がずれた光学系を示しています。これにより、回折限界イメージング システムと比較してコントラストが急激に低下します。コントラストは、約 250 サイクル/mm、つまり周期 4 μm 付近でゼロになっていることがわかります。これが、焦点がずれたシステム (e、f) の画像が、回折限界システム (b、c) の画像よりもぼやけている理由です。焦点が合っていないシステムでは、空間周波数が約 250 サイクル/mm 付近ではコントラストが非常に低いものの、回折限界である 500 サイクル/mm をわずかに下回る空間周波数では、理想的なシステムと同程度のコントラストが得られることに注意してください。パネル (f) の画像を詳しく観察すると、スポークターゲットの中心付近にあるスポーク密度が高い部分の画像が比較的鮮明であることがわかります。
光伝達関数[ 2 ](OTF)は点像分布関数(PSF)のフーリエ変換として定義されるため、一般に空間周波数の複素数値関数となります。特定の周期パターンの投影は、絶対値と複素引数がそれぞれ投影された投影の相対コントラストと並進に比例する複素数で表されます。

多くの場合、コントラストの低下が最も重要であり、パターンの平行移動は無視できます。相対コントラストは、光伝達関数の絶対値で表され、この関数は一般的に変調伝達関数(MTF)と呼ばれます。その値は、空間周波数の関数として、対象物のコントラストが画像にどれだけ捉えられているかを示します。MTFは、空間周波数が1から0(回折限界)まで増加するにつれて減少する傾向がありますが、この関数は単調ではない場合が多いです。一方、パターンの平行移動も重要な場合、光伝達関数の複素引数は、一般に位相伝達関数(PhTF )と呼ばれる2番目の実数値関数として表すことができます。複素数値の光伝達関数は、これら2つの実数値関数の組み合わせとして考えることができます。
どこ
そしては複素引数関数を表し、は周期パターンの空間周波数である。一般にこれは各次元の空間周波数を持つベクトルであり、つまり周期的なパターンの方向も示しています。
焦点が適切に調整された光学系のインパルス応答は、焦点面に最大値を持つ三次元強度分布であり、検出器を軸方向に移動させながら画像スタックを記録することで測定できます。したがって、三次元光伝達関数は、インパルス応答の三次元フーリエ変換として定義できます。通常は一次元、あるいは場合によっては二次元断面のみが用いられますが、三次元光伝達関数を用いることで、構造化照明顕微鏡などの顕微鏡の理解を深めることができます。
伝達関数の定義により、点光源オブジェクトから検出された光の割合を示す必要があります。ただし、通常は検出された光の総量に対するコントラストが最も重要です。したがって、光伝達関数を検出強度で正規化するのが一般的です。。
一般的に、光伝達関数は、放射光のスペクトルや偏光、点光源の位置などの要因に依存します。例えば、画像のコントラストと解像度は通常、画像の中心で最適となり、視野の端に向かうにつれて低下します。大きな変化が生じる場合は、代表的な位置や色を用いて光伝達関数を計算することができます。
二値的な白黒ストライプパターンに基づいて伝達関数を定義する方が実際的な場合もある。等幅の白黒周期パターンの伝達関数は、コントラスト伝達関数(CTF)と呼ばれる。[ 3 ]
理想的なレンズシステムは、周期的なパターンをずらすことなく高コントラストの投影を実現するため、光伝達関数は変調伝達関数と同一になります。通常、コントラストは光学系の解像度によって決まる点で徐々にゼロに近づきます。例えば、可視光波長500nmで使用される、収差のない理想的なf / 4光学結像システムの場合、 光伝達関数は右図に示されています。
グラフから、コントラストが徐々に低下し、空間周波数500サイクル/ミリメートルでゼロになることがわかります。つまり、画像投影の光学解像度は1/500ミリメートルであり、これは2マイクロメートルの特徴サイズに相当します。500サイクル/ミリメートルを超えると、このイメージングシステムのコントラスト、ひいてはその変調伝達関数はゼロになります。したがって、この特定のイメージングデバイスでは、スポークは中心に向かってますますぼやけ、最終的には灰色で解像度の低い円盤状に融合します。
光学伝達関数は、対象空間または試料空間、観測角度、フィルム幅などの単位で表される場合や、理論上の最大値で正規化される場合があることに注意してください。単位間の変換は、通常、乗算または除算によって行われます。例えば、顕微鏡は通常、対象物を10~100倍に拡大し、一眼レフカメラは一般的に、5メートル離れた対象物を100~200倍に縮小します。
デジタルイメージングデバイスの解像度は、光学系だけでなく、画素数、特に画素間の距離によっても制限されます。ナイキスト・シャノン標本化定理で説明されているように、上記の例の光学解像度に合わせるには、各カラーチャンネルの画素を1マイクロメートル(1ミリメートルあたり500サイクルの半分の周期)で分離する必要があります。同じセンサーサイズで画素数を増やしても、より細かいディテールを解像度化することはできません。一方、画素間隔が1マイクロメートルより大きい場合、解像度は画素間の距離によって制限されます。さらに、エイリアシングによって画像の忠実度がさらに低下する可能性があります。
不完全で収差のある結像系は、次の図に示すような光伝達関数を持つ可能性がある。
理想的なレンズシステムでは、コントラストは空間周波数500サイクル/ミリメートルでゼロになります。しかし、それより低い空間周波数では、コントラストは前の例の理想的なシステムよりもかなり低くなります。実際、500サイクル/ミリメートルよりも低い空間周波数でも、コントラストがゼロになる場合がいくつかあります。これが、上の図に示されているスポーク画像の灰色の円形の帯を説明するものです。灰色の帯の間では、スポークが黒から白、または白から黒に反転しているように見えます。これはコントラスト反転と呼ばれ、光学伝達関数の実部の符号反転に直接関係しており、一部の周期的なパターンでは半周期のシフトとして表されます。
理想的なシステムと不完全なシステムの両方の解像度は2μm 、つまり500LP/mmであると主張できるかもしれないが、後者の例の画像の方が鮮明度が低いことは明らかである。知覚される品質により近い解像度の定義としては、最初のゼロが現れる空間周波数、すなわち10μm 、つまり100LP/mmを用いるべきだろう。解像度の定義は、たとえ理想的なイメージングシステムであっても、大きく異なる。より完全で明確な定義は、光伝達関数によって得られる。

光学系、特に光学収差は、必ずしも回転対称ではありません。そのため、周期が同じであっても、向きが異なる周期パターンは、コントラストが異なって画像化される可能性があります。したがって、光伝達関数または変調伝達関数は、一般に2次元関数となります。以下の図は、回転対称性を持たない収差であるトレフォイルを持つ光学系について、先に述べた理想的なシステムと不完全なシステムの2次元等価物を示しています。
光学的伝達関数は必ずしも実数値とは限りません。周期パターンは、システムの収差に応じて任意の量だけシフトする可能性があります。これは一般的に、非回転対称収差の場合に当てはまります。上の図の表面プロットの色相は位相を示しています。回転対称収差の場合、位相は0またはπのいずれかであり、したがって伝達関数は実数値であるのに対し、非回転対称収差の場合、伝達関数は虚数成分を持ち、位相は連続的に変化することがわかります。
カメラシステムに関して一般的に用いられる光学解像度は、画像内のピクセル数、ひいては微細なディテールを表示する可能性のみを表すのに対し、伝達関数は、空間周波数の変動パターンに応じて隣接するピクセルが黒から白に変化する能力、すなわち、コントラストが最大か最小かにかかわらず、微細なディテールを表示する実際の能力を表します。高空間周波数で「ロールオフ」する光学伝達関数で再現された画像は、日常的な言葉で言えば「ぼやけて」見えます。
1920×1080ピクセルの現在の高解像度(HD)ビデオシステムを例にとると、ナイキストの定理によれば、理想的なシステムでは、合計1920本の白黒交互線を完全に(真の黒から白への遷移で)解像できるはずであり、これは画像幅あたり1920/2=960ラインペア、または画像幅あたり960サイクルの空間周波数と呼ばれます(単位角度あたりまたは単位ミリメートルあたりのサイクルで定義することもできますが、カメラを扱う場合は一般的に不明瞭で、望遠鏡などに適しています)。実際には、これは当てはまらず、ナイキスト周波数に近い空間周波数は一般的に振幅が減少して再現されるため、細かいディテールは見えるものの、コントラストが大幅に低下します。このことから、例えば、後述するオーバーサンプリングを用いたフィルムスキャナーから得られた標準画質のテレビ映像は、変調伝達関数が劣るカメラで撮影された高精細映像よりも鮮明に見える可能性があるという興味深い観察結果が得られる。この2つの映像には、しばしば見落とされがちな興味深い違いがある。前者はある一定の点までは細部まで十分なコントラストを備えているが、それ以降は細かいディテールがほとんどないのに対し、後者はより細かいディテールは含まれているものの、コントラストが低下しているため、全体的に劣っているように見えるのである。

一般的に画像は平面、つまり二次元と考えられがちですが、撮像システムは原理的には測定可能な三次元の強度分布を画像空間に生成します。例えば、二次元センサーを移動させることで三次元の強度分布を捉えることができます。点光源の画像もまた三次元(3D)の強度分布であり、3D点像分布関数で表すことができます。例として、右の図は広視野顕微鏡(a)と共焦点顕微鏡(c)の物体空間における3D点像分布関数を示しています。開口数1.49の同じ対物レンズを使用していますが、共焦点顕微鏡の点像分布関数は横方向(x、y)と軸方向(z)の両方においてよりコンパクトであることがわかります。したがって、共焦点顕微鏡の解像度は、三次元すべてにおいて広視野顕微鏡よりも優れていると結論づけることができます。
3次元光伝達関数は、3D点像分布関数の3次元フーリエ変換として計算できます。その色分けされた大きさは、それぞれパネル(a)と(c)に示されている点像分布関数に対応するパネル(b)と(d)にプロットされています。広視野顕微鏡の伝達関数は、 3次元すべてにおいて共焦点顕微鏡の半分のサポートを持ち、広視野顕微鏡の解像度が低いことを以前に確認しました。z軸に沿って、x = y = 0の場合、原点以外では伝達関数がゼロになることに注意してください。この欠落した円錐は、広視野顕微鏡を使用した光学的断層撮影を妨げるよく知られた問題です。[ 4 ]
焦点面における2次元光伝達関数は、3D光伝達関数をz軸に沿って積分することで計算できます。広視野顕微鏡(b)の3D伝達関数はz ≠ 0のときz軸上でゼロですが、その積分である2D光伝達はx = y = 0で最大値に達します。これは、3D光伝達関数が原点x = y = z = 0で発散するためだけに可能です。3D光伝達関数のz軸に沿った関数値はディラックのデルタ関数に対応します。
ほとんどの光学設計ソフトウェアには、レンズ設計の光学伝達関数または変調伝達関数を計算する機能があります。ここに示す例のような理想的なシステムは、Julia、GNU Octave、Matlabなどのソフトウェアを使用して数値的に簡単に計算でき、特定のケースでは解析的に計算することもできます。光学伝達関数は、次の2つのアプローチに従って計算できます。[ 5 ]
数学的には、どちらの手法も同等である。数値計算は通常、フーリエ変換を用いるのが最も効率的であるが、解析計算は自己相関法を用いる方が扱いやすい場合がある。
光伝達関数は点像分布関数のフーリエ変換であり、点像分布関数は逆フーリエ変換された瞳関数の絶対値の二乗であるため、光伝達関数は瞳関数から直接計算することもできます。畳み込み定理から、光伝達関数は実際には瞳関数の自己相関であることがわかります。[ 5 ]
円形開口部を持つ理想的な光学系の瞳関数は、半径1の円盤である。したがって、このようなシステムの光伝達関数は、距離1の2つの同一の円盤の交差領域から幾何学的に計算できる。、 どこ空間周波数は、最高送信周波数で正規化されます。[ 2 ]一般に、光伝達関数は、最大値が1になるように正規化されます。結果として得られる面積は、。
交差する部分の面積は、2つの同一の円弧の面積の合計として計算できます。、 どこは円弧の角度です。等式を用いてそして面積の式は次のように書き換えることができます。したがって、正規化された光伝達関数は次のように表される。
一次元光伝達関数は、線幅関数の離散フーリエ変換として計算できます。このデータは空間周波数データに対してグラフ化されます。この場合、MTF対空間周波数曲線に6次多項式をフィッティングして傾向を示します。50%カットオフ周波数を決定することで、対応する空間周波数が得られます。このようにして、このデータから被試験ユニットの最適焦点のおおよその位置が決定されます。

線幅広がり関数(LSF)のフーリエ変換は、以下の式では解析的に求めることができません。
したがって、フーリエ変換は離散フーリエ変換を用いて数値的に近似される。[ 6 ]
どこ
次に、MTFを空間周波数に対してプロットし、そのグラフからこのテストに関するすべての関連データを読み取ることができます。
顕微鏡などで見られるような高開口数では、光を運ぶ場のベクトル性を考慮することが重要です。デカルト座標軸に対応する3つの独立した成分に波を分解することで、各成分の点像分布関数を計算し、それらを組み合わせてベクトル点像分布関数を求めることができます。同様に、ベクトル光伝達関数は( [ 7 ])および([ 8 ] )に示すように決定できます。
光伝達関数は、光学システムの設計に役立つだけでなく、製造されたシステムの特性評価にも有用である。
光伝達関数は、光学系のインパルス応答のフーリエ変換として定義され、点像分布関数とも呼ばれます。したがって、光伝達関数は、まず点光源の像を取得し、その像に2次元離散フーリエ変換を適用することで容易に得られます。このような点光源としては、例えば、ピンホールのあるスクリーンの背後にある明るい光、蛍光体や金属の微小球、あるいは単にスクリーン上に描かれた点などが挙げられます。点像分布関数を用いた光伝達関数の計算は汎用性が高く、点光源の様々な位置や波長スペクトルに対してこの手順を繰り返すことで、空間的に変化する収差や色収差を持つ光学系を完全に特性評価することができます。
収差が空間的に不変であると仮定できる場合、線やエッジなどの代替パターンを使用して光伝達関数を決定できます。対応する伝達関数は、それぞれ線状拡散関数とエッジ拡散関数と呼ばれます。このような拡張オブジェクトは画像内のより多くのピクセルを照射するため、信号対雑音比が大きくなり、測定精度が向上します。この場合、光伝達関数は画像の2次元離散フーリエ変換を計算し、拡張オブジェクトのフーリエ変換で割ることによって求められます。一般的には、線または白黒のエッジが使用されます。
原点を通る直線の2次元フーリエ変換は、その直線に直交し、かつ原点を通る直線になります。したがって、除数は1次元を除いてすべてゼロとなるため、光伝達関数は単一の線広がり関数(LSF)を用いて1次元についてのみ決定できます。必要に応じて、様々な角度の直線を用いて測定を繰り返すことで、2次元光伝達関数を決定できます。
線広がり関数は、2つの異なる方法で求めることができます。スリットテストターゲットによって得られる理想的な線近似から直接求める方法と、次のサブセクションで説明するエッジ広がり関数から導出する方法です。
エッジの 2 次元フーリエ変換も、エッジに直交する単一の線上でのみ非ゼロになります。この関数は、エッジ スプレッド ファンクション(ESF) と呼ばれることもあります。[ 9 ] [ 10 ]ただし、この線上の値は原点からの距離に反比例します。この手法で取得される測定画像はカメラの広い領域を照らしますが、これは主に低空間周波数での精度に有利です。ライン スプレッド ファンクションと同様に、各測定では光学伝達関数の 1 つの軸しか決定されないため、光学系が回転対称であると仮定できない場合は、繰り返し測定する必要があります。

右図に示すように、オペレータは、黒体によって背面から照明されたナイフエッジテストターゲット画像の端を囲むボックス領域を定義します。ボックス領域は、フレーム全体の領域の約10%に定義されます。画像ピクセルデータは、2次元配列(ピクセル強度とピクセル位置)に変換されます。配列内の各ラインの振幅(ピクセル強度)は正規化され、平均化されます。これにより、エッジスプレッド関数が得られます。
どこ
線広がり関数は、エッジ広がり関数の1階微分と同一であり、 [ 11 ]数値的手法を用いて微分されます。エッジ広がり関数を測定する方が実用的である場合は、次のように線広がり関数を決定できます。
通常、ESFは離散的な点でのみ既知であるため、LSFは有限差分法を用いて数値的に近似される。
どこ:
「鮮明度」は、白黒の線が交互に並んだグリッドパターンで評価されることが多いが、厳密には、黒から白への正弦波変化(通常のパターンをぼかしたバージョン)を用いて測定すべきである。方形波パターン(単純な白黒の線)を用いる場合、エイリアシングのリスクが高まるだけでなく、方形波の基本波成分が方形波自体の振幅よりも大きい(高調波成分によってピーク振幅が減少する)という事実にも注意する必要がある。そのため、方形波テストチャートでは、実際よりも高い空間周波数の解像度を示す、楽観的な結果となる。方形波の結果は、「コントラスト伝達関数」(CTF)と呼ばれることもある。
実際には、多くの要因によって再現画像にかなりのぼやけが生じ、ナイキスト周波数よりわずかに低い空間周波数のパターンは見えなくなる場合があり、最も細かいパターンは白黒ではなく、グレーの濃淡として「ぼやけて」見えることがあります。主な要因の一つは、通常、完璧な「ブリックウォール」光学フィルター(デジタルカメラやビデオカメラでは「位相板」または特定のぼかし特性を持つレンズとして実現されることが多い)を作ることが不可能であることです。このようなフィルターは、ディスプレイのナイキスト周波数を超える空間周波数を除去することでエイリアシングを低減するために必要です。
カメラなどのデジタル画像システムで理論上のシャープネスに実際に近づく唯一の方法は、カメラセンサーのピクセル数を最終画像のサンプル数よりも多くし、ナイキスト周波数を超える高周波をカットオフしてエイリアシングを回避しつつ、その周波数まで比較的フラットなMTFを維持する特殊なデジタル処理を使用して「ダウンコンバート」または「補間」することです。このアプローチは、1970年代にフライングスポットスキャナ、そして後にCCDラインスキャナが開発された際に初めて採用されました。これらは必要なピクセル数よりも多くサンプリングしてからダウンコンバートしたため、映画はビデオカメラで撮影された他の素材よりもテレビで常にシャープに見えました。理論的に正しい補間またはダウンコンバートの唯一の方法は、強力な処理を必要とする2次元sin( x )/ x重み付け関数との畳み込みによって実現される急峻なローパス空間フィルタを使用することです。実際には、処理要件を軽減するために、これに対するさまざまな数学的近似が使用されています。これらの近似は現在、ビデオ編集システムやPhotoshopなどの画像処理プログラムで広く実装されています。
高コントラストMTFの標準画質ビデオがオーバーサンプリングによってのみ実現可能であるのと同様に、理論上の最高の鮮明度を持つHDテレビは、はるかに高い解像度のカメラから始めてデジタルフィルタリングを行うことによってのみ実現可能です。映画が4K、さらには映画館向けに8Kビデオで撮影されるようになった現在、HDTVで最高の画質を得られるのは、より高い基準で撮影された映画や素材に限られると予想されます。カメラで使用されるピクセル数をどれだけ増やしても、完璧な光学空間フィルタが存在しない限り、このことは常に当てはまります。同様に、5メガピクセルの静止画カメラで得られた5メガピクセルの画像は、同等の品質の10メガピクセルの静止画カメラからダウンコンバージョンして得られた5メガピクセルの画像よりも鮮明になることは決してありません。高コントラストMTFを維持するという問題のため、BBCのような放送局は長い間、標準画質テレビを維持しつつ、より多くのピクセルで撮影・視聴することで画質を向上させることを検討していました(ただし、前述のように、このようなシステムは印象的ではありますが、最終的には、減衰はするものの真のHD視聴の効果を高める非常に細かいディテールが欠けてしまいます)。
デジタルカメラやビデオカメラにおけるもう一つの重要な要素は、レンズの解像度です。レンズは1920本の水平線を「解像」できると言われますが、これは黒から白への完全な変調で解像できるという意味ではありません。「変調伝達関数」(位相を無視した光学伝達関数の振幅を表す用語)は、レンズの性能を真に測定するものであり、振幅と空間周波数の関係を示すグラフで表されます。
レンズの絞りによる回折もMTFを制限します。レンズの絞りを小さくすると通常は収差が減り、MTFの平坦性が向上しますが、どのレンズとイメージセンサーのサイズにも最適な絞り値があり、それを超えると回折によって光がイメージセンサー全体に拡散し、解像度が低下します。これはプレートカメラや35mm フィルムの時代にはほとんど問題になりませんでしたが、一部のデジタルカメラ、特にビデオカメラで使用されている非常に小さなフォーマットのセンサーでは克服できない制限となっています。第一世代のHDコンシューマーカムコーダーは1/4インチセンサーを使用しており、約f4より小さい絞り値では解像度が制限され始めます。プロ用ビデオカメラでもほとんどが2/3インチセンサーを使用しており、フィルムフォーマットでは標準と考えられていたf16付近の絞り値の使用は禁止されています。一部のカメラ(ペンタックスK10D など)には「MTF自動露出」モードがあり、絞り値の選択が最大シャープネスになるように最適化されています。通常、これは絞り値範囲の中間あたりを意味します。[ 12 ]
近年、低照度感度と浅い被写界深度効果へのニーズの高まりから、大判デジタル一眼レフカメラの使用が主流になりつつあります。その結果、映画やテレビ番組制作者の中には、その「映画的な」表現力に魅力を感じ、プロ仕様のHDビデオカメラよりもこうしたカメラを好む人も現れています。理論的には、16メガピクセルや21メガピクセルのセンサーを搭載したカメラは、カメラ内部でのダウンコンバージョンとデジタルフィルタリングによるエイリアシング除去によって、ほぼ完璧なシャープネスを実現できる可能性があります。こうしたカメラは非常に優れた結果を生み出し、デジタルフィルタリングを用いた大判ダウンコンバージョンによる、エイリアシングのないフラットなMTFを実現する標準的な手法として、ビデオ制作の分野をリードしているようです。
光学効果により、コントラストが最適ではなくなり、ディスプレイのナイキスト周波数に達する前にゼロに近づくことがあります。光学的なコントラスト低下は、表示または後続処理の前に空間周波数を選択的にデジタル増幅することで部分的に回復できます。より高度なデジタル画像復元手順も存在しますが、ウィーナー逆畳み込みアルゴリズムはその簡便性と効率性からよく使用されます。この手法は画像の空間スペクトル成分を乗算するため、エイリアシングなどによるノイズやエラーも増幅します。したがって、信号対雑音比が十分に高い高品質の録画にのみ有効です。
一般的に、点光源の像を表す点像分布関数は、波長(色)や視野角(点光源の横方向の位置)などの要因にも依存します。このような変化が十分に緩やかな場合、光学系は一連の光伝達関数によって特徴付けられます。しかし、点光源の像が急激に変化する場合、光伝達関数は光学系を正確に記述しません。このような不正確さは、適切な波長または視野位置における光伝達関数の集合によって軽減できる場合が多くあります。ただし、ライトフィールドカメラなどの一部のイメージングシステムでは、より複雑な特性評価が必要になる場合があります。