Commit 05611b92 authored by Russell Luke's avatar Russell Luke
Browse files

OrbitalTomography started!

parent 95620cc7
This diff is collapsed.
import numpy as np
from scipy.io import loadmat
import matplotlib.pyplot as plt
inp=loadmat('../../../InputData/OrbitalTomog/coronen_homo1_fourier_noise15.mat')
imshow(inp['I2'])
# Berechnung und Speicherung als png-Datei des Supports
# Autokorrelation
threshold_autocorr=0.1;
support0=np.fft.fftshift(np.fft.ifft2(np.fft.ifftshift(inp['I2'])));
cond2=(support0 < threshold_autocorr*np.amax(support0)).astype(uint);
# Anfangsbedingungen
support=support0;
F0=inp['I2']**0.5;
ph_init = 2*np.pi*np.random.rand(F0.shape);
# ph_init = np.angle(np.fft.fft2(ph_init));
u0 = F0 * np.exp(1j.*ph_init);
previous = np.fft.fftshift(np.fft.ifft2(np.fft.ifftshift(u0)));
plt.figure()
plt.subplot(121)
plt.imshow(support0.real)
plt.colorbar()
plt.subplot(122)
plt.imshow(cond2)
plt.colorbar()
\ No newline at end of file
This diff is collapsed.
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment