ratt-ru / pyxis

Python Extensions for Interferometry Scripting
3 stars 5 forks source link

im.wsclean with *wsclean-first-residual.fits and **wsclean.restored.fits #74

Open Ermiasabebe opened 8 years ago

Ermiasabebe commented 8 years ago

When I made image after simulation, I found an image with _wsclean-first-residual.fits besides the wsclean.residual.fits. The fluxes in the _wsclean-first-residual.fits* are almost half of that of the _wsclean.restored.fits. That means the fluxes in the initial skymodel will be shared between wsclean-first-residual.fits_ and *wsclean.restored.fits.

I expect that the fluxes in the wsclean.restored.fits has to be almost equal to the initial skymodel flux values, but in my case they are almost half of the fluxes of the initial sky model.

I just wonder how the *wsclean-first-residual.fit appears as I haven't seen such image before while using wsclean?

Ermiasabebe commented 8 years ago

Correction: The fluxes in the wsclean-first-residual.fits are almost equal to that of the wsclean.restored.fits.

modhurita commented 8 years ago

I too have got that first-residual.fits file as one of the outputs from wsclean, though I don't know what its significance is. Which wsclean are you using? The one specified in the 3C147 pipeline doesn't seem to be currently working for me:

mitra@elwood:~/PBandReduction$ /home/oms/src/wsclean-code/wsclean/build/wsclean -size 16384 16384 -scale 2.5asec -niter 10000 -threshold 1 -mgain 0.8 -name img16384_cal PBandData.MS -bash: /home/oms/src/wsclean-code/wsclean/build/wsclean: No such file or directory

Ermiasabebe commented 8 years ago

@modhurita. I am using wsclean version of 1.9. I am not sure which version 3C147 pipeline is using. The lwimager is too slow and it doesn't use "DATA" column to produce image though it is possible to copy the "DATA" column to "CORRECTED_DATA" and produce image.

SpheMakh commented 8 years ago

I usually just ignore that first-residual.fits file

Ermiasabebe commented 8 years ago

@SpheMakh: What I found strange is that the fluxes in the wsclean-first-residual.fits are almost equal to that of the wsclean.restored.fits. The fluxes in the initial skymodel are shared between these two images.