Working_with_Lunar_Reconnaissance_Orbiter_MiniRF_Data

Working with Lunar Reconnaissance Orbiter MiniRF Data ¶

---

• Working with Lunar Reconnaissance Orbiter MiniRF Data
• ISIS3 Processing of MiniRF S-zoom data
  • S-Zoom (LSZ)
    • Data Acquisition
    • Data Ingestion
    • spiceinit (Load NAIF SPICE kernels)
  • Generate Stoke Parameter Images
    • S1
    • S2
    • S3
    • S4
  • Replace invalid pixels along the vertical edge beyond the lunar data for the S1 file
  • Stack the masked S1_NoZ output with S2, S3 and S4
  • Replace invalid pixels for stacked S2, S3 and S4
  • Reduce the 4-band stacked, cleaned cube to approximately 15 meters/pixel (from the original 30 meters/pixel)
  • Create a maptemplate to map project the stacked and reduced file
  • Map project the stacked and reduced file
    • CPR

ISIS3 Processing of MiniRF S-zoom data ¶

---

S-Zoom (LSZ) ¶

---

Data Acquisition ¶

Data Ingestion ¶

Note:
* Both the .LBL and .IMG must be in the same directory together
* The LBL and IMG extensions must be in upper or lower case consistent with each other (.LBL/.IMG or .lbl/.img).

Example command:

mrf2isis from=LSZ_04485_1CD_XKU_87S205_V1.LBL to=LSZ_04485_1CD_XKU_87S205_V1_lev1.cub

The application spiceinit will add the appropriate SPICE information to the ISIS3 image cube. Successful spiceinit allows for further geometric processing of the data.

spiceinit (Load NAIF SPICE kernels) ¶

SPICE (Spacecraft & Planetary ephemeredes, Instrument C-matrix and Event kernels)

Note:
For the moon, ISIS3 currently defaults to the global LOLA for the local topographic shapemodel.

Example command:
Specify the "Smithed" SPK kernels for LRO (not a default in spiceinit)

spiceinit from=LSZ_04485_1CD_XKU_87S205_V1_lev1.cub spksmithed=true

---

Generate Stoke Parameter Images ¶

S1 ¶

fx f1=LSZ_04485_1CD_XKU_87S205_V1_lev1.cub+1  f2=LSZ_04485_1CD_XKU_87S205_V1_lev1.cub+2  
   equation="f1 + f2" 
   to=LSZ_04485_1CD_XKU_87S205_V1_S1.cub

S2 ¶

fx f1=LSZ_04485_1CD_XKU_87S205_V1_lev1.cub+1  f2=LSZ_04485_1CD_XKU_87S205_V1_lev1.cub+2  
   equation="f1 - f2" 
   to=LSZ_04485_1CD_XKU_87S205_V1_S2.cub

S3 ¶

fx f1=LSZ_04485_1CD_XKU_87S205_V1_lev1.cub+3 
   equation="f1 * 2.0" 
   to=LSZ_04485_1CD_XKU_87S205_V1_S3.cub

S4 ¶

fx f1=LSZ_04485_1CD_XKU_87S205_V1_lev1.cub+4 
   equation="f1 * -2.0" 
   to=LSZ_04485_1CD_XKU_87S205_V1_S4.cub

---

Replace invalid pixels along the vertical edge beyond the lunar data for the S1 file ¶

The output pixels will be replaced with NULL
Example Command:

mask from=LSZ_04485_1CD_XKU_87S205_V1_S1.cub to=LSZ_04485_1CD_XKU_87S205_V1_S1_NoZ.cub 
     min=-20.0 max=0.0 preserve=outside spixels=NULL

Stack the masked S1_NoZ output with S2, S3 and S4 ¶

Example Commands:

ls LSZ_04485_1CD_XKU_87S205_V1_S1_NoZ.cub > stack.lis
ls LSZ_04485_1CD_XKU_87S205_V1_S2.cub >> stack.lis
ls LSZ_04485_1CD_XKU_87S205_V1_S3.cub >> stack.lis
ls LSZ_04485_1CD_XKU_87S205_V1_S4.cub >> stack.lis


cubeit fromlist=stack.lis to=LSZ_04485_1CD_XKU_87S205_V1_STACK.cub

Replace invalid pixels for stacked S2, S3 and S4 ¶

The NULL values in the S1_NoZ will applied to the remaining bands

bandtrim from=LSZ_04485_1CD_XKU_87S205_V1_STACK.cub to=LSZ_04485_1CD_XKU_87S205_V1_STACK-TRIM.cub

Reduce the 4-band stacked, cleaned cube to approximately 15 meters/pixel (from the original 30 meters/pixel) ¶

Reduce will average the input pixels which results in an improved output instead of specifying the reduced resolution within a map template for cam2map

reduce from=LSZ_04485_1CD_XKU_87S205_V1_STACK-TRIM.cub to=LSZ_04485_1CD_XKU_87S205_V1_STACK-TRIM_15m.cub 
       lscale=2 sscale=2

Create a maptemplate to map project the stacked and reduced file ¶

Refer to: Learning_About_Map_Projections

Example map template for MiniRF (to apply on a polar image in this example):

Group = Mapping
ProjectionName = POLARSTEREOGRAPHIC
CenterLongitude = 0.0
TargetName = MOON
EquatorialRadius = 1737400.0
PolarRadius = 1737400.0
LatitudeType = Planetocentric
LongitudeDirection = PositiveEast
LongitudeDomain = 360
PixelResolution = 14.806323449292
Scale = 2048.0
CenterLatitude = -90.0
End_Group

Map project the stacked and reduced file ¶

For a mosaic, remember to make sure all images are projected to the same projection, clat, clon and map scale/resolution.
For resolution, specify the scale in meters/pixel in cam2map.

Example cam2map command:

cam2map from=LSZ_04485_1CD_XKU_87S205_V1_STACK-TRIM_15m.cub to=LSZ_04485_1CD_XKU_87S205_V1_STACK-TRIM_15m_map.cub 
        pixres=camera resolution=2048 warpalgorithm=forwardpatch patchsize=5 map=predefined.map

To create a mosaic, refer to: Power Tip: Making Mosaics Learning_About_Map_Projections

CPR ¶

Create an 8-bit CPR cube from the Stokes-1 and Stokes-4 bands of the stacked cube or the multi-band mosaic:

fx f1=LSZ_04485_1CD_XKU_87S205_V1_STACK-TRIM_15m_map.cub+1 
   f2=LSZ_04485_1CD_XKU_87S205_V1_STACK-TRIM_15m_map.cub+4 
   to=LSZ_04485_1CD_XKU_87S205_V1_STACK-TRIM_15m_map_CPR.cub+UnsignedByte+0.0:2.53 
   equation="(f1-f2)/f1+f2)" mode=cubes