Commit a5bffb79 authored by Lukas Platz's avatar Lukas Platz

added signal and residual plot

parent 5d2241a3
...@@ -43,7 +43,7 @@ if __name__ == '__main__': ...@@ -43,7 +43,7 @@ if __name__ == '__main__':
# #
# # One dimensional regular grid with uniform exposure # # One dimensional regular grid with uniform exposure
# position_space = ift.RGSpace(1024) # position_space = ift.RGSpace(1024)
# exposure = np.ones(position_space.shape) # exposure = ift.Field.full(position_space, 1.)
# Two-dimensional regular grid with inhomogeneous exposure # Two-dimensional regular grid with inhomogeneous exposure
position_space = ift.RGSpace([512, 512]) position_space = ift.RGSpace([512, 512])
...@@ -99,6 +99,10 @@ if __name__ == '__main__': ...@@ -99,6 +99,10 @@ if __name__ == '__main__':
H, convergence = minimizer(H) H, convergence = minimizer(H)
# Plot results # Plot results
ift.plot(sky(H.position), title='Reconstruction') signal = sky(mock_position)
reconst = sky(H.position)
ift.plot(signal, title='Signal')
ift.plot(GR.adjoint(data), title='Data') ift.plot(GR.adjoint(data), title='Data')
ift.plot(reconst, title='Reconstruction')
ift.plot(reconst - signal, title='Residuals')
ift.plot_finish(name='getting_started_2.png', xsize=16, ysize=16) ift.plot_finish(name='getting_started_2.png', xsize=16, ysize=16)
Markdown is supported
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