Документ взят из кэша поисковой машины. Адрес оригинального документа : http://xmm.vilspa.esa.es/calibration/documentation/CALHB/node575.html
Дата изменения: Tue Sep 20 19:45:20 2011
Дата индексирования: Mon Apr 11 03:12:00 2016
Кодировка:
Поисковые слова: http www.badastronomy.com phpbb index.php

Calibration Access and Data Handbook

Procedure

This routine is used by the SAS.
The PSF is stored as a radial distribution in the CCF component PSF1DRB. The routine reads the PSF from the PSF1DRB CCF component, where the PSF with the closest CFRR to the input countFrameRateRatio is selected. The radial width of the PSF is stretched by the factor f$_{stretch}$ by multiplying the radial bin limits with f$_{stretch}$:

$$RMIN \rightarrow RMIN\cdot f_{stretch}\\ RMAX \rightarrow RMAX\cdot f_{stretch}\\$$

The stretching factor f$_{stretch}$ is calculated for the specified field position (rawX, rawY) from a 3rd order polynomial.

$$f_{stretch}= a_0 + a_1 Y + a_2 Y^2 + a_3 Y^3 \\ a_4 X + a_5 X Y + a_6 X Y^2 +\\ a_7 X^2 + a_8 X^2Y + \\ a_9 X^3 \\$$ (35)

where $X = (rawX-1024.5)^i and Y=(rawY-1024.5).$
The 10 polynomial coefficients $a_0, ..., a_9$ are listed as keywords in the extension headers:

keyword	variable name
	in equ. 35
PARX0Y0	$a_0$
PARX0Y1	$a_1$
PARX0Y2	$a_2$
PARX0Y3	$a_3$
PARX1Y0	$a_4$
PARX1Y1	$a_5$
PARX1Y2	$a_6$
PARX2Y0	$a_7$
PARX2Y1	$a_8$
PARX3Y0	$a_9$

The PSF map is generated from the radial distribution in the following steps:

1. identify PSF1DRB extension associated with FilterId
2. read polynomial coefficients in header
3. calculate the stretching factor $f_{strech}$ for the specified field position (rawX, rawY)
4. read binary table extension
5. apply stretching factor $f_{strech}$
6. initialize PSF map as specified by input variables:
   (nPixelsNegX+nPixelsPosX+1, nPixelsNegY+nPixelsPosY+1)
7. for each pixel of the PSF map
   • calculate the minimum distance r$_{min}$ of the pixel to the PSF centre (rawX, rawY)
   • calculate the maximum distance r$_{max}$ of the pixel to the PSF centre (rawX, rawY)
   • calculate the normalization factor f$_{norm}$:
     
     
     $$f_{norm}={1\over [(RMAX-RMIN)^2\cdot\pi] \cdot xRes \cdot yRes}$$
     
     

   • set the PSF pixel value to the interpolated intensity value at the radius $r=r_{min}+(r_{max}-r_{min})/2$
   • normalize PSF pixel value with the normalization factor f$_{norm}$
8. normalize PSF map to 1