We should preserve wcs when manipulating so_map object. This line
|
noise.data = np.random.randn(size[0], size[1]) * rms_uKarcmin_T / np.sqrt(pixArea) |
removes
wcs information : we must use
noise.data[:] to only modify the
numpy array and keep
wcs unchanged