質問編集履歴
1
コードの追加
test
CHANGED
|
File without changes
|
test
CHANGED
|
@@ -8,7 +8,7 @@
|
|
|
8
8
|
|
|
9
9
|
現在、ケプストラムを導くことはできたのですが、そこからPSFを推定するための方法が分からず詰まっています。ご教授していただけると助かります。
|
|
10
10
|
|
|
11
|
-
また、python3系での実装を考えています。
|
|
11
|
+
また、python3系での実装を考えています。また、ケプストラムは以下のコードで導きました。
|
|
12
12
|
|
|
13
13
|
|
|
14
14
|
|
|
@@ -23,3 +23,51 @@
|
|
|
23
23
|
求めたケプストラム
|
|
24
24
|
|
|
25
25
|
|
|
26
|
+
|
|
27
|
+
|
|
28
|
+
|
|
29
|
+
|
|
30
|
+
|
|
31
|
+
```python
|
|
32
|
+
|
|
33
|
+
import sys
|
|
34
|
+
|
|
35
|
+
import numpy as np
|
|
36
|
+
|
|
37
|
+
from matplotlib import pyplot as plt
|
|
38
|
+
|
|
39
|
+
from PIL import Image
|
|
40
|
+
|
|
41
|
+
from scipy import signal
|
|
42
|
+
|
|
43
|
+
import cv2
|
|
44
|
+
|
|
45
|
+
|
|
46
|
+
|
|
47
|
+
rekka_img = "bure1.PNG"
|
|
48
|
+
|
|
49
|
+
|
|
50
|
+
|
|
51
|
+
imput_img = cv2.imread(rekka_img,cv2.IMREAD_GRAYSCALE)
|
|
52
|
+
|
|
53
|
+
|
|
54
|
+
|
|
55
|
+
kepusutoramu = np.fft.fft2(imput_img)
|
|
56
|
+
|
|
57
|
+
kepusutoramu = np.log(np.abs(kepusutoramu))
|
|
58
|
+
|
|
59
|
+
kepusutoramu = np.fft.ifft2(kepusutoramu).real
|
|
60
|
+
|
|
61
|
+
kepusutoramu = np.fft.ifftshift(kepusutoramu)
|
|
62
|
+
|
|
63
|
+
|
|
64
|
+
|
|
65
|
+
kepusutoramu_c = np.clip((kepusutoramu-np.min(kepusutoramu)), 0, np.min(kepusutoramu)*-2.0)
|
|
66
|
+
|
|
67
|
+
plt.imshow((kepusutoramu_c), cmap = 'gray')
|
|
68
|
+
|
|
69
|
+
plt.title('rekka Cepstrum'), plt.xticks([]), plt.yticks([])
|
|
70
|
+
|
|
71
|
+
plt.show()
|
|
72
|
+
|
|
73
|
+
```
|