Benjamin_Loison/Robust_image_source_identification_on_modern_smartphones
1
0
Fork 0
Code Issues 68 Pull Requests Packages Projects Releases Wiki Activity

Add titles

Browse Source
This commit is contained in:
Benjamin Loison 2024-05-14 02:12:54 +02:00
parent db2e19738e
commit 337a6eca84
Signed by: Benjamin_Loison
SSH Key Fingerprint: SHA256:BtnEgYTlHdOg1u+RmYcDE0mnfz1rhv5dSbQ2gyxW8B8
3 changed files with 13 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

1
datasets/raise/fft/fourier_to_image.py
View File

@ -1,5 +1,6 @@
from PIL import Image
import numpy as np
from scipy import fftpack
fft = np.array(Image.open('fft.png').convert('L'))# / 255
print(fft)

11
datasets/raise/fft/image_to_fourier.py
View File

@ -1,5 +1,6 @@
from PIL import Image
import numpy as np
from scipy import fftpack
image = np.array(Image.open('image_4.png').convert('L'))# / 255
print(image)
@ -8,8 +9,16 @@ print(image)
# fftpack.fftshift
fft1 = (fftpack.fft2(image))
height, width = fft1.shape
for x in range(width):
fft1[height // 2, x] = np.median(fft1[:, x])
# save the image
# np.log10
plt.imshow(np.log10(abs(fft1)))
plt.show()
#plt.imsave('ifft1.png', fft1)
#plt.imsave('ifft1.png', np.log10(abs(fft1)))

3
datasets/raise/fft/remove_period_patterns.py
View File

@ -47,8 +47,9 @@ fixedImage = interpolate_replace_nans(fft1, kernel)
figure, axes = plt.subplots(1, 2, sharex = True, sharey = True)
axes[0].set_title('Original FFT')
axes[0].imshow(np.log10(originalFft1))
axes[1].set_title('FFT with significant lines attenuated')
axes[1].imshow(np.log10(fixedImage))
#plt.imshow(np.log10(fixedImage))
plt.show()
#plt.imsave('fft.png', np.log10(fixedImage))