Local Emission

Description

This is the process to generate local emission for each LROC NAC image. It requires the DTM to have already been generated. This is a circuitous process that involves ISIS and SPC and ends with an ASCII image.

Procedure

Overview

ISIS to Emission sub-block

• NOTE: We generate the emission data first since this shell script first initializes the LROC NAC EDR, which will be used by incidence and I/F. The shell script here generates phoE.IMG. ISISe.sh

• Generating the IMG takes awhile, so we have a wrapper to run all the images at once. NOTE: This wrapper requires directory setup. Each image has its own folder; the folder name is the last 5 digits of the NAC, plus the L or R for left/right. The wrapper should be placed in the directory containing all the image folders.

Convert LRO NAC emission image to SPC sub-block

Replace with correct SUMFILES

• Copy the SUMFILE from the topography SPC directory to the “working” directory in the Cube generation directory. Nothing crazy here, but you will need to change the filename from NXXXXX.SUM to NEXXXXX.SUM (for emission) or NPXXXXX.SUM (for phase). For I/F, you can leave the filename the same.
• No need to change the contents of the SUMFILE, though if you want to keep the standard SUMFILE format, the first line (the SUMFILE name) should match the new filename. I don’t change the first line since it doesn’t affect the results.

Extract Pixels in ROI

Convert to Local Emission

• I have a program that needs a lot of setup for each image and bigmap. The base version is dumpMapEmiss.f, but I create a new file for each location on the Moon. I am only giving the overall idea here…
• Since LRO is a line scan, the angle between the surface and the s/c is always east or west, with no north or south component. Update dumpMapEmiss.f to include various parameters that solve for the local emission angle for each image. Since the emission image in ASCII format is essentially a “scalar” value (i.e. doesn’t include east or west), determining the parameters is a manual process that requires the user to determine if the s/c was east or west of the western portion of the bigmap. You also have to set a parameter which checks for an inflection point in the emission image in ASCII format. This step is required because for some of the images, the s/c to topography vector flips from east to west, or vice versa. The flip occurs because the image is mostly nadir pointing.
• The program uses the DN from the emission image in ASCII format, along with the user defined east/west, and does vector addition with the normal of the topography at each pixel to determine the local emission angle. The program starts at the western edge of the bigmap and works east, checking for an inflection point in the emission image in ASCII format. Other special situations require additional parameters, such as when an image strip is diagonal with respect to the bigmap, such that the image strip only covers a portion of the western edge of the bigmap. When this happens it confuses the check for an inflection point, so an exception must be explicitly defined. There is almost certainly a more effective way to program the check for the inflection point so it doesn’t require so much manual massaging, but this is what I came up with in the allotted time.