質問編集履歴

1

コードの追加

2021/01/04 07:00

投稿

shocora
shocora

スコア8

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
  ![イメージ説明](a693a4a05eeef9ac81a1fea03dfcad48.png)
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
+ ```