ミス発見

Webでも見れる「Principles of Computerized Tomographic Imaging」にミス発見。

http://www.slaney.org/pct/pct-toc.html

Chapter.3、p.83、(96)式が間違ってる。

g(n\alpha)=\left{ \begin{array}{ll}\frac{1}{8\alpha^2},&\textrm{$n=0$} \\ 0, &\textrm{$n$ is even}\\(\frac{\alpha}{\pi \alpha sin n\alpha})^2, &\textrm{$n$ is odd.} \end{array}\right

正しくは、

g(n\alpha)=\left{ \begin{array}{ll}\frac{1}{8\alpha^2},&\textrm{$n=0$} \\ 0, &\textrm{$n$ is even}\\-\frac{1}{2}(\frac{1}{\pi sin n\alpha})^2, &\textrm{$n$ is odd.} \end{array}\right

となる。
そのままプログラム組んでたので結果がおかしかった。
よく考えるとカーネルが全部正になるわけないじゃん。これじゃ平滑化フィルタの一種だ。という話。
まあまだオレのプログラム自体はバグが残ってて結果がおかしいわけだが。