AlxisDataBlog-数学・物理学


MathJax = {
  chtml: {
    matchFontHeight: false
  },
  tex: {
    inlineMath: [['$', '$']]
  }
};






以前 OTF の説明で PSF のスペクトルを計算しました。

\begin{align}
\mathscr{F}\mathrm{PSF}=\dfrac{g\ast\tilde{g}}{\|g\|_2^2}
\end{align}

$g$ は開口関数をスケーリングしたものであり台有界です。したがって $\tilde{g}$ も台有界となり $g,\tilde{g}$ を畳込んだ $g\ast \tilde{g}$ も当然のことながら台有界となります。

\begin{align}
\{\sqrt{p^2+q^2}:(p,q)\in\mathrm{supp}(g\ast\tilde{g})\}
\end{align}

の最大値を遮断周波数と呼びます。物理学的には開口は有界なので、常に遮断周波数は存在します。試しにエアリーパターンの遮断周波数を考えてみましょう。開口関数は直径 $d_0$ の円板定義関数でしたね。このとき $g$ は直径 $1/(\lambda F)$ の円板定義関数になります。したがって $g\ast\tilde{g}$ の台は半径 $1/(\lambda F)$ の円板。このとき遮断周波数は $1/(\lambda F)$ です。


光学像の畳込みモデルでは遮断周波数は極めて重要な係数と考えられます。

\begin{align}
h=\mathrm{PSF}\ast f
\end{align}

畳込みモデルのフーリエ像は、畳込み定理により

\begin{align}
\mathscr{F}h=\mathrm{OTF}\mathscr{F}f
\end{align}

となります。遮断周波数を超えるスペクトル領域で、右辺は零である故に光学像の空間周波数である左辺も零になってしまいます。


零となってしまった周波数からは $\mathscr{F}f$ へ回復させることができないという問題が生じます。今回はここまで。次回からは、話題を変えて弊社独自のデコンボリューションの話をします。(written by J.Nishigata)


Copy Right 2020-現在 ALXIS DATA, Inc. All Right Reserved.

遮断周波数

AlxisDataBlog-数学・物理学


MathJax = {
  chtml: {
    matchFontHeight: false
  },
  tex: {
    inlineMath: [['$', '$']]
  }
};






任意の凸体 $K$ の定義関数を $k$ とします。$k$ による畳込み

\begin{align}
\int_{\mathbb{R}^2}k(x)f(y-x)\,dS
\end{align}

で許可されている運動は並進運動です。これを剛体運動まで拡張した

\begin{align}
\mathscr{P}f(g)=\int_{\mathbb{R}^2}(g^{-1}k)(x)f(x)\,dS
\end{align}

はポンペイユ変換と呼ばれています。ここに $g$ は剛体運動群の元です。ポンペイユ変換に関する以下の問いを考えてみましょう。


ポンペイユの順問題: $K$ が円であるならばポンペイユ変換は単射であるだろうか。


円板定義関数のフーリエ像は $J_1(\|\cdot\|)/\|\cdot\|$ であり、原点を中心とする円周状の零点をもちます。この円に対してポンペイユ変換の自明ではない零点が得られる故にポンペイユ変換は単射にはなりません。


ポンペイユの順問題は否定的に解かれました。改めて正しい命題を書いてみると「$K$ が円であるならばポンペイユ変換は単射ではない」となるでしょう。


この真の命題から対偶をとり「ポンペイユ変換が単射であるならば $K$ は円ではない」、さらに逆の命題にしてみましょう。


ポンペイユの逆問題: $K$ が円ではないならばポンペイユ変換は単射であるだろうか。


この問題は約 $100$ 年間ものあいだ未解決のまま残る、非常に有名な問題です。ただし、未解決とはいっても完全に手付かずの未解決というわけではありません。特殊な $K$ や制限された $f$ のクラスに対しては肯定的な結果が得られています。


今回はここまで。次回は特殊な場合として、$K$ が正方形の場合を考えてみましょう。(written by J.Nishigata)


Copy Right 2020-現在 ALXIS DATA, Inc. All Right Reserved.

デコンボリューションの問題設定