CASSINI Procedures

To ingest images, you first need to setup the directory. The programs are (to JRW's knowledge) unique to CASSINI.

Directory Structure

Saturn

• CASSINI
  • DATA
  • UTILITIES
• Dione
• Mimas
• Etc

Contents of the CASSINI Directory

• DATA has all the kernels in it, with no subdirectories. Will also need cumindex.tab (Cumulative Image Index provided by the PDS) in this directory, as well as TABLE.TXT (used to convert images to 12 bit).
• UTILITIES has three files.

  • list.e - User inputs object name (i.e. DIONE, TETHYS, etc. all caps!!!) and the minimum size of the object in pixels. Use this to flush out low resolution images (list.e scales the size of the body) and make a list of images to download and an input file for the next utility. Source is here list.f

  • convert.e - Enter the object name as "OBJ.TXT" where OBJ is the object name you entered previously. Puts the DAT image files into IMAGEFILES and makes make_sumfiles.in. Source is here convert.f

  • make_sumfiles.e - Makes the SUMFILES. When it asks if you want to "Include Pole (y/n)" SAY NO!!! Bob toyed with the idea of putting the pole into the SUMFILE, but never implemented that. Source is here make_sumfiles.f

# Below is from Bob's email 29 July 2020
../CASSINI/UTILITIES/list.e
DIONE
100   <- minimum size in pixels
  [makes DIONE.TXT and b/w coverage plot]
  [makes ./NEW_IMAGES/download.b]

cd NEW_IMAGES
chmod +x download.b
./download.b

cd ..

../CASSINI/UTILITIES/convert.e     <- this is equivalent to PROCESS_IMG
DIONE.TXT

  [makes images and puts into IMAGEFILES]
  [makes make_sumfiles.in]

../CASSINI/UTILITIES/make_sumfiles.e
 0.  Create new SUMFILES & NOMINALS and ..
 1.  Update old NOMINALS and ..
 2.  Update old SUMFILES 
0
 Include Pole? (y/n)
n                                  <- ALWAYS SAY NO!

make_sumfiles.txt and INIT_LITHOS.TXT

Here is the top portion of make_sumfiles.txt with the camera model for the Cassini ISS. "608 IAPETUS" and "IAU_IAPETUS BFRAME" will need to be changed appropriately.

2                                                  NCAM
'CASSINI_ISS_NAC'
1024  1024                                         NPXA, NLNA
2002.703  512.5  512.5                             MMFLA, PX0A, LN0A
83.3333, 0, 0, 0, 83.3428, 0                       KMATA
 0, 8.28E-06, 5.45E-06, -19.67E-06                 DA
160, 4095                                          T1,T2
'CASSINI_ISS_WAC'
1024  1024                                         NPXB, NLNB
200.7761  512.5  512.5                             MMFLB, PX0B, LN0B
83.3333, 0, 0, 0, 83.3411, 0                       KMATB
 0, 60.88E-06, 5.28E-06, -71.86E-06                DB
160, 4095                                          T1,T2
1                                                  NSC
-82                                                CASSINI
1.d-1,1.d-1,1.d-1                                  sigV0
3.d-4,3.d-4,3.d-4                                  sigPT
608                                                IAPETUS
'IAU_IAPETUS'                                      BFRAME
'LT+S'                                             ABCORR
KERNELS
../CASSINI/DATA/naif0012.tls
../CASSINI/DATA/cas_iss_v10.ti
../CASSINI/DATA/cas_v40.tf
../CASSINI/DATA/cas_v42.tf
../CASSINI/DATA/cpck31Oct2017.tpc
../CASSINI/DATA/cas00172.tsc
../CASSINI/DATA/de421.bsp                                                                                        
../CASSINI/DATA/sat286.bsp                                                                                       
../CASSINI/DATA/051212AP_RE_90165_14363.bsp
../CASSINI/DATA/041014R_SCPSE_01066_04199.bsp
../CASSINI/DATA/041219R_SCPSE_04199_04247.bsp
<snip>

Example INIT_LITHOS.TXT is here INIT_LITHOS.TXT.

Potential Errors

When you run the output from list.e (i.e. NEW_IMAGES/download.b) you get an <Name>.IMG file for each image, even if there was an error in the download. If you don't notice there was an error and run convert.e, you get an error that looks like this ...

<rand> Iapetus$ ../CASSINI/UTILITIES/convert.e 
 INPUT LIST
IAPETUS.TXT
At line 118 of file convert.f (unit = 10, file = 'NEW_IMAGES/N1483150865_1.IMG')
Fortran runtime error: Non-existing record number

This happened to JRW because the directory structure on the PDS server is different from that in cumindex.tab. You'll have to re-download the files.