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; }