23.04.2015 Views

ISOCAM Interactive Analysis User's Manual Version 5.0 - ISO - ESA

ISOCAM Interactive Analysis User's Manual Version 5.0 - ISO - ESA

ISOCAM Interactive Analysis User's Manual Version 5.0 - ISO - ESA

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

262 CHAPTER 20. ADVANCED DATA CALIBRATION<br />

• Project the images contained in the FITS files and place the results in the file result.fits:<br />

CIA> spawn, !cia_exec+’/projection -i2,input -o result.fits’<br />

• Take a look at the result:<br />

CIA> result = readfits(’result.fits’, hdr)<br />

CIA> tviso, result[*, *, 0]<br />

20.15.4 Back projection<br />

The projection routines also allow the possibility to back project a raster MOSAIC to a set<br />

of EXPOSUREs. A real source will, due to the many effects present in the <strong><strong>ISO</strong>CAM</strong> data,<br />

have a somewhat different signal in each EXPOSURE. The purpose of the data calibration is to<br />

minimize or eliminate this difference. Knowledge of how the signal varies from EXPOSURE to<br />

EXPOSURE is not contained in the raster MOSAIC, so the back-projected EXPOSUREs will<br />

each contain an averaged source signal. Figures 20.7 and 20.8 illustrate this point. Note that<br />

this averaged source signal is really an idealized signal, i.e. it assumes that the same source has<br />

the same signal in each EXPOSURE. This assumption can be useful for testing the quality of<br />

the data analysis. The less the signal in the original EXPOSUREs deviates from the idealized<br />

signal in the back-projected EXPOSUREs the better the quality of the data calibration.<br />

The back projection is performed with the routine back project. You should of course only<br />

perform back projection on a fully calibrated PDS after creation of the raster MOSAIC:<br />

CIA> raster_scan, raster_pds<br />

CIA> help, raster_pds.image<br />

FLOAT = Array[32, 32, 8]<br />

CIA> back_project, raster_pds, images<br />

CIA> help, images<br />

IMAGES FLOAT = Array[32, 32, 8]<br />

You can use stat to check the quality of your calibration. The lower the RMS the less the<br />

averaged or idealized signal deviates from the real calibrated and corrected signal:<br />

CIA> stat, images - raster_pds.image<br />

-------------------------------------------<br />

Image dimensions: 32 32<br />

Number of frames: 8<br />

Total number of pixels: 8192<br />

-------------------------------------------<br />

Minimum Maximum Mean Median RMS<br />

-10.7506 14.4660 0.144570 0.00626230 0.902133

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!