quatrope / ProperImage

Proper Image utils for astronomy.
BSD 3-Clause "New" or "Revised" License
26 stars 6 forks source link

Strange background flux around transients #78

Open Astro-Sean opened 2 years ago

Astro-Sean commented 2 years ago

Hi,

Getting nice results with your code - however I have a strange issue when recovering transient flux - the background around the transient is around 4x times less than the image.

I'm using the code snippet from the tutorial. Setting Beta=True gives these results, without it I get strange residuals of the host galaxy. Any helps would be much appreciated.

result = subtract( ref=ref, new=sci, smooth_psf=False, fitted_psf=True, align=False, iterative=False, beta=True, shift=False)

image
BrunoSanchez commented 2 years ago

Hi! This looks strange. It seems related to the background when the code measures the PSF. As if the PSF would tend to sit on a negative background or similar... What is the mask ourtput?

pythonPleiades commented 1 year ago

Hello. I'm having the same issue and mostly at galaxy cores and sometime at bright stars. Any help to minimize or remove those artifacts would be appreciated. Cheers.

BrunoSanchez commented 1 year ago

Hi, can you give me more details? What is the mask output? I am not actively working on the code for a long time, but would be nice to hear more about this problem.

Is it possible to get some reproducible case? can you tar images + script? Thanks

pythonPleiades commented 1 year ago

Hello Bruno. Thanks for the reply. Similar to Astro-Sean (previous message), I'm using the code from the tutorial: result = subtract( ref=ref, new=sci, smooth_psf=True, fitted_psf=True, align=False, iterative=False, beta=True, shift=False). I tried different subtraction options but similar results. Here the result with Scorr_diff from NGC 3556:

scorr_diff

Most of the time, I get dark square with a bright star inside, at center. Here the mask result:

mask

Here a copy of the Science and Reference files: NGC 3556_Reference.zip Thanks again for your help.