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

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

16.4. MANIPULATING THE MASK 195<br />

and .CUBE – these fields contains the vast bulk of the data in a PDS, so removing them greatly<br />

decreases the PDS size and saves disk and memory space.<br />

CIA> raster_pds = prune_pds( raster_pds )<br />

CIA> save, file=’raster_pds.dat’, raster_pds, /verb<br />

% SAVE: Portable (XDR) SAVE/RESTORE file.<br />

% SAVE: Saved variable: raster_pds.<br />

raster pds is a raster PDS which is originally more than 2 MB in size. prune pds reduces<br />

it to about 0.5 MB.<br />

CIA> $ls -la raster_pds.dat<br />

-rw-r--r-- 1 mdelaney ssamr 500568 Oct 14 13:17 raster_pds.dat<br />

16.4 Manipulating the MASK<br />

This section describes CIA MASK manipulation routines. For more details on the MASK in<br />

CIA data structures and a list of the possible MASK values see Section 15.2.2.18. Also take a<br />

look at Section 2.3.4 for how to configure the MASK handling.<br />

• ia put mask routine to modify the MASK in a CIA data structure or IDL array. To work<br />

with ordinary IDL arrays you can use put mask.<br />

• ia get mask routine to extract the MASK from a CIA data structure.<br />

16.4.1 Extracting the MASK from CIA data structures<br />

Follow the procedure below to extract and examine a MASK from a CIA data structure.<br />

1. Suppose we have a CIA structure. The first step is to extract the entire MASK from the<br />

structure.<br />

• If our structure was an SCD this can be done as:<br />

CIA> mask = scd_get( ’mask’, ’CSCD143006010110_97092611592803’ )<br />

• If it is a PDS:<br />

CIA> mask = pds.mask<br />

2. Now, to look at the dead pixels use ia get mask to extract and tviso to display.<br />

CIA> dead_mask = ia_get_mask( ’dead’, mask )<br />

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

If you look at the values in dead mask you will see that they are either 0 or 1. However<br />

the values in original mask can be 0,1,2,4,6,8,16,32,64,128, each value having a separate<br />

meaning (see the table in Section 15.2.2.18). ia get mask translates these values to a simple<br />

0 or 1. In our example above, we asked for the pixels masked as dead to be translated, so<br />

all the pixels with a value of 2 in original mask became 1 in dead mask.

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

Saved successfully!

Ooh no, something went wrong!