Commit 86d7c515 authored by Martin Reinecke's avatar Martin Reinecke
Browse files

Merge branch 'demo_2' into 'redesign'

added signal and residual plot to getting_started_2

See merge request ift/nifty-dev!84
parents c9a6a709 a5bffb79
...@@ -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])
...@@ -98,6 +98,10 @@ if __name__ == '__main__': ...@@ -98,6 +98,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)
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