回答編集履歴
2
誤記訂正
answer
CHANGED
|
@@ -10,9 +10,7 @@
|
|
|
10
10
|
|
|
11
11
|
以下、煩雑さを避けるために解析対象の画像が正方形であり、垂直方向の符号は数学でおなじみの座標系(複素平面)に合わせる前提とします。
|
|
12
12
|
|
|
13
|
-
|
|
13
|
+
画像の辺の画素数をWとすると、dftshift後のスペクトル配列spectrumの要素spectrum[W // 2 + re, W // 2 - im]は、複素数 c = re + 1j * im を用いて以下のような波の成分(スペクトル)となります。
|
|
14
|
-
|
|
15
|
-
複素数 c = re + 1j * im を用いて以下のような波の成分(スペクトル)となります。
|
|
16
14
|
|
|
17
15
|
・方向: cの偏角
|
|
18
16
|
・波長: W/|c| (単位は画素)
|
1
日本語が変なところを訂正
answer
CHANGED
|
@@ -17,7 +17,7 @@
|
|
|
17
17
|
・方向: cの偏角
|
|
18
18
|
・波長: W/|c| (単位は画素)
|
|
19
19
|
|
|
20
|
-
numpyには複素数の偏角を求めるangle、複素数の絶対値(複素平面原点からの距離)を
|
|
20
|
+
numpyには複素数の偏角を求めるangle、複素数の絶対値(複素平面原点からの距離)を求めるabsなどの関数があります。これらの機能を利用すると次のように集計できます。(本件では同じ波長範囲を集計すればよく、コードでは波長ではなく周波数で範囲をチェックしています)
|
|
21
21
|
|
|
22
22
|
```python
|
|
23
23
|
import numpy as np
|