UV image simulation from Galex

This is an algorithm to generate frames from Galex image.


You make a table (dimentions 4x1000) for Poisson statistic. x-axis is average photons at a pixel P_ij and Y-axis is the number of photons. X-axis values range from P_ij = 0.001, 0.002,... 1.000, in steps of 0.001 (1000 values). Y-axis has 4 values, which are equal to probabilities of number of photons incident = 1, 2, 3, 4.

For example: at x = P_ij = 0.1, the four values (corresponding to probability of 1,2,3 and 4 photons incident on that pixel) are 0.1, 0.035, 0.002, and 0.0005.

For a given frame_k, for each pixel

1. Find a random number R equally distributed between 0 and 1.


2. If R > 0.1, then number of photons for that pixel = 0. If R > 0.035, 0.002 and 0.0005, then number of photons incident on that pixel are 1, 2, 3.