void IsisMain ()
{
stretch.ClearPairs();
for (int i=0; i<6; i++) {
gapCount[i] = 0;
suspectGapCount[i] = 0;
invalidCount[i] = 0;
lisCount[i] = 0;
hisCount[i] = 0;
validCount[i] = 0;
}
void TranslateHiriseEdrLabels (Filename &labelFile, Cube *);
void SaveHiriseCalibrationData (ProcessImportPds &process, Cube *,
Pvl &pdsLabel);
void SaveHiriseAncillaryData (ProcessImportPds &process, Cube *);
void FixDns8 (Buffer &buf);
void FixDns16 (Buffer &buf);
ProcessImportPds p;
Pvl pdsLabel;
UserInterface &ui = Application::GetUserInterface();
// Get the input filename and make sure it is a HiRISE EDR
Filename inFile = ui.GetFilename("FROM");
iString id;
bool projected;
try {
Pvl lab(inFile.Expanded());
id = (string) lab.FindKeyword ("DATA_SET_ID");
projected = lab.HasObject("IMAGE_MAP_PROJECTION");
}
catch (iException &e) {
string msg = "Unable to read [DATA_SET_ID] from input file [" +
inFile.Expanded() + "]";
throw iException::Message(iException::Io,msg, _FILEINFO_);
}
//Checks if in file is rdr
if( projected ) {
string msg = "[" + inFile.Name() + "] appears to be an rdr file.";
msg += " Use pds2isis.";
throw iException::Message(iException::User,msg, _FILEINFO_);
}
id.ConvertWhiteSpace();
id.Compress();
id.Trim(" ");
if (id != "MRO-M-HIRISE-2-EDR-V1.0") {
string msg = "Input file [" + inFile.Expanded() + "] does not appear to be " +
"in HiRISE EDR format. DATA_SET_ID is [" + id + "]";
throw iException::Message(iException::Io,msg, _FILEINFO_);
}
p.SetPdsFile (inFile.Expanded(), "", pdsLabel);
// Make sure the data we need for the BLOBs is saved by the Process
p.SaveFileHeader();
p.SaveDataPrefix();
p.SaveDataSuffix();
// Let the Process create the output file but override any commandline
// output bit type and min/max. It has to be 16bit for the rest of hi2isis
// to run.
// Setting the min/max to the 16 bit min/max keeps all the dns (including
// the 8 bit special pixels from changing their value when they are mapped
// to the 16 bit output.
CubeAttributeOutput &outAtt = ui.GetOutputAttribute("TO");
outAtt.PixelType (Isis::SignedWord);
outAtt.Minimum((double)VALID_MIN2);
outAtt.Maximum((double)VALID_MAX2);
Cube *ocube = p.SetOutputCube(ui.GetFilename("TO"), outAtt);
p.StartProcess ();
TranslateHiriseEdrLabels (inFile, ocube);
// Pull out the lookup table so we can apply it in the second pass
// and remove it from the labels.
// Add the UNLUTTED keyword to the instrument group so we know
// if the lut has been used to convert back to 14 bit data
PvlGroup &instgrp = ocube->GetGroup("Instrument");
PvlKeyword lutKey = instgrp["LookupTable"];
PvlSequence lutSeq;
lutSeq = lutKey;
// Set up the Stretch object with the info from the lookup table
// If the first entry is (0,0) then no lut was applied.
if ((lutKey.IsNull()) ||
(lutSeq.Size()==1 && lutSeq[0][0]=="0" && lutSeq[0][1]=="0")) {
stretch.AddPair(0.0, 0.0);
stretch.AddPair(65536.0, 65536.0);
instgrp.AddKeyword(PvlKeyword("Unlutted","TRUE"));
instgrp.DeleteKeyword ("LookupTable");
}
// The user wants it unlutted
else if (ui.GetBoolean("UNLUT")) {
for (int i=0; i<lutSeq.Size(); i++) {
stretch.AddPair(i, (((double)lutSeq[i][0] + (double)lutSeq[i][1]) / 2.0));
}
instgrp.AddKeyword(PvlKeyword("Unlutted","TRUE"));
instgrp.DeleteKeyword ("LookupTable");
}
// The user does not want the data unlutted
else {
stretch.AddPair(0.0, 0.0);
stretch.AddPair(65536.0, 65536.0);
instgrp.AddKeyword(PvlKeyword("Unlutted","FALSE"));
}
ocube->PutGroup(instgrp);
// Save the calibration and ancillary data as BLOBs. Both get run thru the
// lookup table just like the image data.
SaveHiriseCalibrationData (p, ocube, pdsLabel);
SaveHiriseAncillaryData (p, ocube);
// Save off the input bit type so we know how to process it on the
// second pass below.
Isis::PixelType inType = p.PixelType();
// All finished with the ImportPds object
p.EndProcess ();
// Make another pass thru the data using the output file in read/write mode
// This allows us to correct gaps, remap special pixels and accumulate some
// counts
lsbGap = ui.GetBoolean("LSBGAP");
ProcessByLine p2;
string ioFile = ui.GetFilename("TO");
CubeAttributeInput att;
p2.SetInputCube(ioFile, att, ReadWrite);
p2.Progress()->SetText("Converting special pixels");
section = 4;
p2.StartProcess((inType == Isis::UnsignedByte) ? FixDns8 : FixDns16);
p2.EndProcess();
// Log the results of the image conversion
PvlGroup results("Results");
results += PvlKeyword ("From", inFile.Expanded());
results += PvlKeyword ("CalibrationBufferGaps", gapCount[0]);
results += PvlKeyword ("CalibrationBufferLIS", lisCount[0]);
results += PvlKeyword ("CalibrationBufferHIS", hisCount[0]);
results += PvlKeyword ("CalibrationBufferPossibleGaps", suspectGapCount[0]);
results += PvlKeyword ("CalibrationBufferInvalid", invalidCount[0]);
results += PvlKeyword ("CalibrationBufferValid", validCount[0]);
results += PvlKeyword ("CalibrationImageGaps", gapCount[1]);
results += PvlKeyword ("CalibrationImageLIS", lisCount[1]);
results += PvlKeyword ("CalibrationImageHIS", hisCount[1]);
results += PvlKeyword ("CalibrationImagePossibleGaps", suspectGapCount[1]);
results += PvlKeyword ("CalibrationImageInvalid", invalidCount[1]);
results += PvlKeyword ("CalibrationImageValid", validCount[1]);
results += PvlKeyword ("CalibrationDarkGaps", gapCount[2]);
results += PvlKeyword ("CalibrationDarkLIS", lisCount[2]);
results += PvlKeyword ("CalibrationDarkHIS", hisCount[2]);
results += PvlKeyword ("CalibrationDarkPossibleGaps", suspectGapCount[2]);
results += PvlKeyword ("CalibrationDarkInvalid", invalidCount[2]);
results += PvlKeyword ("CalibrationDarkValid", validCount[2]);
results += PvlKeyword ("ObservationBufferGaps", gapCount[3]);
results += PvlKeyword ("ObservationBufferLIS", lisCount[3]);
results += PvlKeyword ("ObservationBufferHIS", hisCount[3]);
results += PvlKeyword ("ObservationBufferPossibleGaps", suspectGapCount[3]);
results += PvlKeyword ("ObservationBufferInvalid", invalidCount[3]);
results += PvlKeyword ("ObservationBufferValid", validCount[3]);
results += PvlKeyword ("ObservationImageGaps", gapCount[4]);
results += PvlKeyword ("ObservationImageLIS", lisCount[4]);
results += PvlKeyword ("ObservationImageHIS", hisCount[4]);
results += PvlKeyword ("ObservationImagePossibleGaps", suspectGapCount[4]);
results += PvlKeyword ("ObservationImageInvalid", invalidCount[4]);
results += PvlKeyword ("ObservationImageValid", validCount[4]);
results += PvlKeyword ("ObservationDarkGaps", gapCount[5]);
results += PvlKeyword ("ObservationDarkLIS", lisCount[5]);
results += PvlKeyword ("ObservationDarkHIS", hisCount[5]);
results += PvlKeyword ("ObservationDarkPossibleGaps", suspectGapCount[5]);
results += PvlKeyword ("ObservationDarkInvalid", invalidCount[5]);
results += PvlKeyword ("ObservationDarkValid", validCount[5]);
// Write the results to the log
Application::Log(results);
return;
}
// Construct a BLOb to contain the Hirise main line suffix and prefix data
void SaveHiriseAncillaryData (ProcessImportPds &process, Cube *ocube) {
vector<int> ConvertCalibrationPixels (int samples,
Isis::PixelType pixelType,
unsigned char *data);
// Setup a Table to hold the main image prefix/suffix data
TableField gap("GapFlag", TableField::Integer);
TableField line("LineNumber", TableField::Integer);
TableField buffer("BufferPixels", TableField::Integer, 12);
TableField dark("DarkPixels", TableField::Integer, 16);
TableRecord rec;
rec += gap;
rec += line;
rec += buffer;
rec += dark;
Table table("HiRISE Ancillary", rec);
table.SetAssociation (Table::Lines);
// Loop through all the prefix and suffix data and construct the table records
// In the case of HiRISE there is only one band so the outside vector
// only contains one entry. The inside vector contains nl entries.
vector<vector<char *> > pre = process.DataPrefix();
vector<vector<char *> > suf = process.DataSuffix();
vector<char *> prefix = pre.at(0);
vector<char *> suffix = suf.at(0);
Progress progress;
progress.SetText("Saving ancillary data");
progress.SetMaximumSteps(prefix.size());
progress.CheckStatus();
for (unsigned int l=0; l<prefix.size(); l++) {
unsigned char *linePrefix = (unsigned char *)(prefix[l]);
// Pull out the gap byte (in byte 0)
rec[0] = (int)(linePrefix[0]);
// Pull out the line number (bytes 3-5 3=MSB, 5=LSB)
int lineCounter = 0;
lineCounter += ((int)(linePrefix[3])) << 16;
lineCounter += ((int)(linePrefix[4])) << 8;
lineCounter += ((int)(linePrefix[5]));
rec[1] = lineCounter;
// Pull the 12 buffer pixels (same type as image data)
// from the image prefix area
linePrefix += 6;
section = 3;
rec[2] = ConvertCalibrationPixels (12, process.PixelType(),linePrefix);
linePrefix += 12 * SizeOf(process.PixelType());
// Pull the 16 dark pixels (same type as image data)
// from the image suffix area
unsigned char *lineSuffix = (unsigned char *)(suffix[l]);
section = 5;
rec[3] = ConvertCalibrationPixels (16, process.PixelType(),lineSuffix);
lineSuffix += 16 * SizeOf(process.PixelType());
// Add this record to the table
table += rec;
// Report the progress
progress.CheckStatus();
}
// Add the table to the output cube
ocube->Write(table);
}
void IsisMain (){
//initialize globals
summed = false;
summedOutput = NULL;
// Grab the file to import
ProcessImportPds p;
UserInterface &ui = Application::GetUserInterface();
Filename inFile = ui.GetFilename("FROM");
Filename out = ui.GetFilename("TO");
// Make sure it is a Galileo SSI image
Pvl lab(inFile.Expanded());
//Checks if in file is rdr
if( lab.HasObject("IMAGE_MAP_PROJECTION") ) {
string msg = "[" + inFile.Name() + "] appears to be an rdr file.";
msg += " Use pds2isis.";
throw iException::Message(iException::Io,msg, _FILEINFO_);
}
// data set id value must contain "SSI-2-REDR-V1.0"(valid SSI image)
// or "SSI-4-REDR-V1.0"(reconstructed from garbled SSI image)
string dataSetId;
dataSetId = (string)lab["DATA_SET_ID"];
try {
if (dataSetId.find("SSI-2-REDR-V1.0") == string::npos
&& dataSetId.find("SSI-4-REDR-V1.0") == string::npos) {
string msg = "Invalid DATA_SET_ID [" + dataSetId + "]";
throw iException::Message(iException::Pvl,msg,_FILEINFO_);
}
}
catch (iException &e) {
string msg = "Unable to read [DATA_SET_ID] from input file [" +
inFile.Expanded() + "]";
throw iException::Message(iException::Io,msg,_FILEINFO_);
}
// set summing mode
if(ui.GetString("FRAMEMODE") == "AUTO") {
double frameDuration = lab["FRAME_DURATION"];
// reconstructed images are 800x800 (i.e. not summed)
// even though they have frame duration of 2.333
// (which ordinarily indicates a summed image)
if (dataSetId.find("SSI-4-REDR-V1.0") != string::npos) {
summed = false;
}
else if (frameDuration > 2.0 && frameDuration < 3.0) {
summed = true;
}
// seti documentation implies valid frame duration values are 2.333, 8.667, 30.333, 60.667
// however some images have value 15.166 (see example 3700R.LBL)
else if (frameDuration > 15.0 && frameDuration < 16.0) {
summed = true;
}
}
else if(ui.GetString("FRAMEMODE") == "SUMMED") {
summed = true;
}
else {
summed = false;
}
Progress prog;
Pvl pdsLabel;
p.SetPdsFile(inFile.Expanded(),"",pdsLabel);
//Set up the output file
Cube *ocube;
if(!summed) {
ocube = p.SetOutputCube("TO");
p.StartProcess();
}
else {
summedOutput = new Cube();
summedOutput->SetDimensions(p.Samples()/2, p.Lines()/2, p.Bands());
summedOutput->SetPixelType(p.PixelType());
summedOutput->Create(ui.GetFilename("TO"));
p.StartProcess(TranslateData);
ocube = summedOutput;
}
TranslateLabels(pdsLabel, ocube);
p.EndProcess ();
if(summed) {
summedOutput->Close();
delete summedOutput;
}
return;
}