next up previous
Next: An easy approach: the Up: FUSE & Optimal extraction Previous: Introduction

Description of optimal extraction

The purpose of this web page is to understand optimal extraction algorithm. I hope that after reading it, you will agree that it is worth using it for the data extraction. Hence, this web page will focus on the algorithm, and not on the reasons which lead us to it.

Here is the algorithm:

displaymath234

x and y are the coordinates of each pixel. x is corresponding to the wavelength.

In fact, behind this equation lay a simple idea: putting a weight on each pixels so the ones which are the most important prevail in case of low S/N ratio. Thus, we can rewrite the algorithm:

displaymath242

displaymath244

displaymath246

tex2html_wrap_inline248 is the weight for each pixel, and tex2html_wrap_inline250 is a normalization factor, allowing tex2html_wrap_inline252 to be conservative when the weights are changing.

tex2html_wrap_inline254 is the dataset coming from one aperture.
Example of dataset for segment LIF 1a:
tex2html_wrap232

tex2html_wrap_inline256 correspond to the background.

tex2html_wrap_inline258 is the probability for a photon with a given wavelength to come on that pixel. This knowledge is coming from a calibrated file named wgts*.fit. Note that tex2html_wrap_inline260 .
Example of the calibrated file tex2html_wrap_inline258 for segment LIF 1a:
tex2html_wrap276

tex2html_wrap_inline252 and tex2html_wrap_inline266 are the flux on a column x. tex2html_wrap_inline266 is an estimation of the flux and tex2html_wrap_inline252 is the flux obtained by optimal extraction. It means that before doing optimal extraction, you need to have an estimation of the flux, obtained here by using box extraction. Then, to increase the effectiveness, the algorithm is repeated several times, using the last flux results as an estimation for the new flux calculation. By default, optimal extraction is executed 3 times.



Sylvestre Lacour
Mon Nov 5 22:07:31 EST 2001