NOTES: ====== ORAC-DR will produce individual images (for each of the pol vectors) labelled _I, _P, _th, etc. It will also produce a "catalogue" file, .FIT: for example, gu20020101_123_I.sdf, gu20020101_123_P.sdf, gu20020101_123.FIT, etc. 1. Copy files to something more meaningful... > cp gu20020101_123.FIT smm1.FIT > cp gu20020101_123_I.sdf smm1_I.sdf 2. Bin up values in the catalogue file: > polbin IN - Input catalogue /@m87j_bin/ > smm1 <<------ smm1 BOX - Bin size /5/ > 3 METHOD - Binning method /'mean'/ > OUT - Output catalogue > smm1_bin <<------ smm1_bin 320 vectors written to the output catalogue. 3. Create a "sub"-catalogue which contains only select values, in this case for pixels where Intensity > 0 but Polarisation percent is < 0, and polarisation is at least 3-sigma > catselect CATIN - Input catalogue: /@smm1_avcube(,,1)/ > smm1_bin CATOUT - Output catalogue: > smm1_thr <<------ smm1_thr REJCAT - Rejects catalogue? /NO/ > SELTYP - Selection type: /''/ > e EXPR - Selection expression: > I>0 & P/dp > 3 & P<10 88 objects selected. 4. Display image (created by ORAC-DR) and overplot vectors: > display smm1_I MODE - Method to define the scaling limits /'per'/ > PERCENTILES - Percentiles for scaling /[10,90]/ > Data will be scaled from -0.0009998139 to 0.002607805. > lutgrey > polplot smm1_thr noclear frame=pixel vscale=8 Alignment has occurred within the PIXEL Domain. 88 vectors plotted. 5. Finally, create image from the smoothed catalogue file created earlier: > polimage IN - Input catalogue /@unpol_K_cat/ > smm1_thr COLDAT - Name of catalogue column holding data values /'U'/ > P OUT - Output image /@unpol_K_U/ > smm1_thr_P <<------ smm1_thr_P COLDAT - Name of catalogue column holding data values /'ang'/ > DP OUT - Output image /@smm1_thr_ang/ > @smm1_thr_dP COLDAT - Name of catalogue column holding data values /'z'/ > ang OUT - Output image /@dp/ > smm1_thr_ang > ls *_thr*.sdf smm1_thr_P.sdf smm1_thr_dP.sdf smm4_thr_ang.sdf smm1_thr_ang.sdf smm4_thr_P.sdf smm4_thr_dP.sdf