Sorry for starting a new thread but i wasn't sure if the experts were still watching an old thread.
Eden, Richard & Jean kindly helped take two functions (psf & pixel response) and convolved them (by oversampling), create beautiful graphics and downsample to give what a CCD would detect. Eden also gave a very nice integration function that allowed the half total integral distance to be determined.
I have taken this further and now allowed for non-point objects - modelled real objects as gaussian for now. Eden commented on the need for recentring when performing convolution, but it appears if the convolution involves 3 functions then this is not needed. Is this correct or is my convolution wrong.
Also, the plot at the end for fractional volume radius displays 2 plots, one is the raw (unconvolved) psf, and the other is the total convolution. For some reason the raw psf (red curve) changes when the pixel value changes which it shouldn't do. I have spent a week trying to see why with no success.
Jason
