| Re: How to create synthetic images of stars with gaussian psf in IDL. [message #92987 is a reply to message #92969] |
Tue, 05 April 2016 07:18   |
Jeremy Bailin
Messages: 618
Registered: April 2008
|
Senior Member |
|
|
On Monday, April 4, 2016 at 10:59:09 AM UTC-4, Sonu Tabitha wrote:
> On Monday, April 4, 2016 at 1:13:27 PM UTC+5:30, Sonu Tabitha wrote:
>> I want to create a synthetic image to resemble a stellar field (say covering an area of 255 by 255 pixels, with a total of about 4 stars having Gaussian PSF). The rest of the image has to be filled with constant background and noise, as present in the entire image, but devoid of stars. I want to create this to test a source extraction algorithm that I have developed. I am a beginner in IDL. Can you please help me out in coding?
>
> Thanks Jeremy! But I am very new to IDL and I am not so familiar with its functions. Can you please give me a sample code?
1) RANDOMU takes as its first argument a random seed, and then the dimensions of the output array. So, for example, if you wanted to get a 256 x 512 random array, you could do:
seed = 1L ; set it to something so your results are repeatable
noise_image = RANDOMU(seed, 256, 512)
By default RANDOMU gives you a uniform random deviate, but you probably want a Poisson distribution, which requires giving the POISSON keyword with the mean expected value. For example, if your noise level is 2.5 electrons, you might say:
RANDOMU(seed, 256, 512, POISSON=2.5)
2) PSF_GAUSSIAN takes a few options. The ones you care about are:
NPIXEL: The dimensions of the output array. So you could use [256, 512] to make an array the same size as above.
CENTROID: The position on the iamge of the center of the PSF. For example [10.5, 186] for a center at x=10.5, y=186 (note they don't have to be at pixel centers).
FWHM or ST_DEV: Width of the Gaussian -- you can specify it either using the FWHM or Gaussian sigma. For an isotropic PSF, just give one value (e.g. ST_DEV=3.0 for a 3-pixel standard deviation), or for an anisotropic PSF you can give a separate x- and y-width (e.g. FWHM=[4.0, 5.5] for a PSF that is slightly elongated in the y-direction).
e.g. psf1 = PSF_GAUSSIAN(NPIXEL=[256,512], CENTROID=[10.5,186], ST_DEV=3.0)
3) Add them.
image = noise_image + psf1
-Jeremy.
|
|
|
|
comp.lang.idl-pvwave archive
Members
Search
Help
Login
Home
