/**
* Calculates the lat/lon of the ControlNet.
*
* @param incubes The filename of the list of cubes in the ControlNet
* @param cnet The filename of the ControlNet
*/
void setControlPointLatLon(SerialNumberList &snl, ControlNet &cnet) {
CubeManager manager;
manager.SetNumOpenCubes(50); //Should keep memory usage to around 1GB
Progress progress;
progress.SetText("Calculating Lat/Lon");
progress.SetMaximumSteps(cnet.GetNumPoints());
progress.CheckStatus();
for (int cp = 0; cp < cnet.GetNumPoints(); cp++) {
ControlPoint *point = cnet.GetPoint(cp);
ControlMeasure *cm = point->GetRefMeasure();
Cube *cube = manager.OpenCube(snl.FileName(cm->GetCubeSerialNumber()));
try {
cube->camera()->SetImage(cm->GetSample(), cm->GetLine());
g_surfacePoints[point->GetId()] = cube->camera()->GetSurfacePoint();
}
catch (IException &e) {
QString msg = "Unable to create camera for cube file [";
msg += snl.FileName(cm->GetCubeSerialNumber()) + "]";
throw IException(e, IException::Unknown, msg, _FILEINFO_);
}
cube = NULL; //Do not delete, manager still has ownership
progress.CheckStatus();
}
manager.CleanCubes();
}
/**
* Removes control points not in the lat/lon range provided in the unput
* parameters.
*
* @param outNet The output control net being removed from
* @param noLanLonPoint The keyword recording all of the control points removed
* due to the provided lat/lon range
* @param noLanLonPoint The keyword recording all of the control points removed
* due to the inability to calculate the lat/lon for that
* point
*/
void ExtractLatLonRange( ControlNet & outNet, PvlKeyword & nonLatLonPoints,
PvlKeyword & cannotGenerateLatLonPoints, map<iString,iString> sn2filename ) {
if( outNet.Size() == 0 ) { return; }
UserInterface &ui = Application::GetUserInterface();
// Get the lat/lon and fix the range for the internal 0/360
double minlat = ui.GetDouble("MINLAT");
double maxlat = ui.GetDouble("MAXLAT");
double minlon = ui.GetDouble("MINLON");
if( minlon < 0.0 ) { minlon += 360; }
double maxlon = ui.GetDouble("MAXLON");
if( maxlon < 0.0 ) { minlon += 360; }
bool useNetwork = ui.GetBoolean("USENETWORK");
Progress progress;
progress.SetText("Calculating lat/lon");
progress.SetMaximumSteps(outNet.Size());
progress.CheckStatus();
CubeManager manager;
manager.SetNumOpenCubes( 50 ); //Should keep memory usage to around 1GB
for( int cp = outNet.Size()-1; cp >= 0; cp --) {
progress.CheckStatus();
// If the Contorl Network takes priority, use it
double pointLat = outNet[cp].UniversalLatitude();
double pointLon = outNet[cp].UniversalLongitude();
bool useControlNet = useNetwork && pointLat > -1000 && pointLon > -1000;
if( outNet[cp].Type() == Isis::ControlPoint::Ground || useControlNet ) {
if( NotInLatLonRange( outNet[cp].UniversalLatitude(),
outNet[cp].UniversalLongitude(),
minlat, maxlat, minlon, maxlon ) ) {
nonLatLonPoints += outNet[cp].Id();
outNet.Delete( cp );
}
}
/**
* If the lat/lon cannot be determined from the point, then we need to calculate
* lat/lon on our own
*/
else if( ui.WasEntered("FROMLIST") ) {
// Find a cube in the Control Point to get the lat/lon from
int cm = 0;
iString sn = "";
double lat = 0.0;
double lon = 0.0;
double radius = 0.0;
// First check the reference Measure
if( outNet[cp].HasReference() ) {
cm = outNet[cp].ReferenceIndex();
if( !sn2filename[outNet[cp][cm].CubeSerialNumber()].empty() ) {
sn = outNet[cp][cm].CubeSerialNumber();
}
}
// Search for other Control Measures if needed
if( sn.empty() ) {
// Find the Serial Number if it exists
for( int cm = 0; (cm < outNet[cp].Size()) && sn.empty(); cm ++ ) {
if( !sn2filename[outNet[cp][cm].CubeSerialNumber()].empty() ) {
sn = outNet[cp][cm].CubeSerialNumber();
}
}
}
// Connot fine a cube to get the lat/lon from
if( sn.empty() ) {
cannotGenerateLatLonPoints += outNet[cp].Id();
outNet.Delete( cp );
}
// Calculate the lat/lon and check for validity
else {
bool remove = false;
Cube *cube = manager.OpenCube( sn2filename[sn] );
Camera *camera = cube->Camera();
if (camera == NULL) {
try {
Projection *projection = ProjectionFactory::Create( (*(cube->Label())) );
if(!projection->SetCoordinate(outNet[cp][cm].Sample(),outNet[cp][cm].Line())) {
nonLatLonPoints += outNet[cp].Id();
remove = true;
}
lat = projection->Latitude();
lon = projection->Longitude();
radius = projection->LocalRadius();
delete projection;
projection = NULL;
} catch ( iException &e ) {
remove = true;
e.Clear();
}
}
else {
if(!camera->SetImage(outNet[cp][cm].Sample(),outNet[cp][cm].Line())) {
nonLatLonPoints += outNet[cp].Id();
remove = true;
}
lat = camera->UniversalLatitude();
lon = camera->UniversalLongitude();
radius = camera->LocalRadius();
camera = NULL;
}
cube = NULL;
if( remove || NotInLatLonRange( lat, lon, minlat, maxlat, minlon, maxlon ) ) {
nonLatLonPoints += outNet[cp].Id();
outNet.Delete( cp );
}
else { // Add the reference lat/lon/radius to the Control Point
outNet[cp].SetUniversalGround( lat, lon, radius );
}
}
}
else {
cannotGenerateLatLonPoints += outNet[cp].Id();
outNet.Delete( cp );
}
}
manager.CleanCubes();
}
void IsisMain() {
// Get user interface
UserInterface &ui = Application::GetUserInterface();
bool register_ignored = ui.GetBoolean("REGISTERIGNOREDONLY");
// Open the files list in a SerialNumberList for
// reference by SerialNumber
SerialNumberList files(ui.GetFilename("FILES"));
// Create a ControlNet from the input file
ControlNet inNet(ui.GetFilename("CNET"));
// Create an AutoReg from the template file
Pvl pvl(ui.GetFilename("TEMPLATE"));
AutoReg *ar = AutoRegFactory::Create(pvl);
// Create the output ControlNet
ControlNet outNet;
outNet.SetType(inNet.Type());
outNet.SetUserName(Application::UserName());
outNet.SetDescription(inNet.Description());
outNet.SetCreatedDate(iTime::CurrentLocalTime());
outNet.SetTarget(inNet.Target());
outNet.SetNetworkId(inNet.NetworkId());
Progress progress;
progress.SetMaximumSteps(inNet.Size());
progress.CheckStatus();
int ignored=0, unmeasured=0, registered=0, unregistered=0, validated=0;
CubeManager cubeMgr;
cubeMgr.SetNumOpenCubes(50);
// Register the points and create a new
// ControlNet containing the refined measurements
for (int i=0; i<inNet.Size(); i++) {
ControlPoint &inPoint = inNet[i];
ControlPoint outPoint;
outPoint.SetType(inPoint.Type());
outPoint.SetId(inPoint.Id());
outPoint.SetUniversalGround(inPoint.UniversalLatitude(), inPoint.UniversalLongitude(), inPoint.Radius());
outPoint.SetHeld(inPoint.Held());
outPoint.SetIgnore(inPoint.Ignore());
// CHECK TO SEE IF THE CONTROL POINT SHOULD BE REGISTERED
// "Ignore" point and we are not registering ignored
if (inPoint.Ignore() && !register_ignored){
ignored++;
// add "Ignored" to network only if indicated
if (ui.GetBoolean("OUTPUTIGNORED")) {
// only include appropriate control measures
for (int j = 0; j < inPoint.Size(); j++) {
if (inPoint[j].IsMeasured()){
outPoint.Add(inPoint[j]);
}
else{
unmeasured++;
if (ui.GetBoolean("OUTPUTUNMEASURED")){
outPoint.Add(inPoint[j]);
}
}
}
// only add this point if OUTPUTIGNORED
outNet.Add(outPoint);
}
// go to next control point
continue;
}
// Not "Ignore" point (i.e. "valid") and we are only registering "Ignored"
else if (!inPoint.Ignore() && register_ignored) {
// add all "valid" points to network
// only include appropriate control measures
for (int j = 0; j < inPoint.Size(); j++) {
if (inPoint[j].IsMeasured()){
outPoint.Add(inPoint[j]);
}
else{
unmeasured++;
if (ui.GetBoolean("OUTPUTUNMEASURED")) {
outPoint.Add(inPoint[j]);
}
}
}
// add this point since it is not ignored
outNet.Add(outPoint);
// go to next control point
continue;
}
// "Ignore" point or "valid" point to be registered
else { // if ( (inPoint.Ignore() && register_ignored) || (!inPoint.Ignore() && !register_ignored ) ) {
if (inPoint.Ignore()) { outPoint.SetIgnore(false); }
ControlMeasure &patternCM = inPoint[inPoint.ReferenceIndex()];
Cube &patternCube = *cubeMgr.OpenCube(files.Filename(patternCM.CubeSerialNumber()));
ar->PatternChip()->TackCube(patternCM.Sample(), patternCM.Line());
ar->PatternChip()->Load(patternCube);
if (patternCM.IsValidated()) validated++;
if (!patternCM.IsMeasured()) continue;
if(!patternCM.IsReference()) {
patternCM.SetReference(true);
patternCM.SetChooserName("Application pointreg");
patternCM.SetDateTime();
}
outPoint.Add(patternCM);
// reset goodMeasureCount for this point before looping measures
int goodMeasureCount = 0;
// Register all the unvalidated measurements
for (int j = 0; j < inPoint.Size(); j++) {
// don't register the reference, go to next measure
if (j == inPoint.ReferenceIndex()){
if (!inPoint[j].Ignore()) goodMeasureCount++;
continue;
}
// if the measurement is valid, keep it as is and go to next measure
if (inPoint[j].IsValidated()) {
validated++;
outPoint.Add(inPoint[j]);
if (!inPoint[j].Ignore()) goodMeasureCount++;
continue;
}
// if the point is unmeasured, add to output only if necessary and go to next measure
if (!inPoint[j].IsMeasured()) {
unmeasured++;
if (ui.GetBoolean("OUTPUTUNMEASURED")) {
outPoint.Add(inPoint[j]);
}
continue;
}
ControlMeasure searchCM = inPoint[j];
// refresh pattern cube pointer to ensure it stays valid
Cube &patternCube = *cubeMgr.OpenCube(files.Filename(patternCM.CubeSerialNumber()));
Cube &searchCube = *cubeMgr.OpenCube(files.Filename(searchCM.CubeSerialNumber()));
ar->SearchChip()->TackCube(searchCM.Sample(), searchCM.Line());
try {
ar->SearchChip()->Load(searchCube,*(ar->PatternChip()),patternCube);
// If the measurements were correctly registered
// Write them to the new ControlNet
AutoReg::RegisterStatus res = ar->Register();
double score1, score2;
ar->ZScores(score1, score2);
searchCM.SetZScores(score1, score2);
if(res == AutoReg::Success) {
registered++;
searchCM.SetType(ControlMeasure::Automatic);
searchCM.SetError(searchCM.Sample() - ar->CubeSample(), searchCM.Line() - ar->CubeLine());
searchCM.SetCoordinate(ar->CubeSample(),ar->CubeLine());
searchCM.SetGoodnessOfFit(ar->GoodnessOfFit());
searchCM.SetChooserName("Application pointreg");
searchCM.SetDateTime();
searchCM.SetIgnore(false);
outPoint.Add(searchCM);
goodMeasureCount++;
}
// Else use the original marked as "Estimated"
else {
unregistered++;
searchCM.SetType(ControlMeasure::Estimated);
if(res == AutoReg::FitChipToleranceNotMet) {
searchCM.SetError(inPoint[j].Sample() - ar->CubeSample(), inPoint[j].Line() - ar->CubeLine());
searchCM.SetGoodnessOfFit(ar->GoodnessOfFit());
}
searchCM.SetChooserName("Application pointreg");
searchCM.SetDateTime();
searchCM.SetIgnore(true);
outPoint.Add(searchCM);
}
} catch (iException &e) {
e.Clear();
unregistered++;
searchCM.SetType(ControlMeasure::Estimated);
searchCM.SetChooserName("Application pointreg");
searchCM.SetDateTime();
searchCM.SetIgnore(true);
outPoint.Add(searchCM);
}
}
// Jeff Anderson put in this test (Dec 2, 2008) to allow for control
// points to be good so long as at least two measure could be
// registered. When a measure can't be registered to the reference then
// that measure is set to be ignored where in the past the whole point
// was ignored
if (goodMeasureCount < 2) {
if (!outPoint.Held() && outPoint.Type() != ControlPoint::Ground) {
outPoint.SetIgnore(true);
}
}
// Otherwise, ignore=false. This is already set at the beginning of the registration process
// Check to see if the control point has now been assigned
// to "ignore". If not, add it to the network. If so, only
// add it to the output if the OUTPUTIGNORED parameter is selected
// 2008-11-14 Jeannie Walldren
if (!outPoint.Ignore()) {
outNet.Add(outPoint);
}
else{
ignored++;
if (ui.GetBoolean("OUTPUTIGNORED")) outNet.Add(outPoint);
}
}
progress.CheckStatus();
}
// If flatfile was entered, create the flatfile
// The flatfile is comma seperated and can be imported into an excel
// spreadsheet
if (ui.WasEntered("FLATFILE")) {
string fFile = Filename(ui.GetFilename("FLATFILE")).Expanded();
ofstream os;
os.open(fFile.c_str(),ios::out);
os << "PointId,OriginalMeasurementSample,OriginalMeasurementLine," <<
"RegisteredMeasurementSample,RegisteredMeasurementLine,SampleDifference," <<
"LineDifference,ZScoreMin,ZScoreMax,GoodnessOfFit" << endl;
os << NULL8 << endl;
for (int i=0; i<outNet.Size(); i++) {
// get point from output control net and its
// corresponding point from input control net
ControlPoint outPoint = outNet[i];
ControlPoint *inPoint = inNet.Find(outPoint.Id());
if (outPoint.Ignore()) continue;
for (int i = 0; i<outPoint.Size();i++) {
// get measure and find its corresponding measure from input net
ControlMeasure cmTrans = outPoint[i];
ControlMeasure cmOrig = (*inPoint)[cmTrans.CubeSerialNumber()];
double inSamp = cmOrig.Sample();
double inLine = cmOrig.Line();
double outSamp = cmTrans.Sample();
double outLine = cmTrans.Line();
double sampErr = cmTrans.SampleError();
double lineErr = cmTrans.LineError();
double zScoreMin = cmTrans.GetZScoreMin();
if (fabs(zScoreMin) <= DBL_EPSILON || zScoreMin == NULL8) zScoreMin = 0;
double zScoreMax = cmTrans.GetZScoreMax();
if (fabs(zScoreMax) <= DBL_EPSILON || zScoreMax == NULL8) zScoreMax = 0;
double goodnessOfFit = cmTrans.GoodnessOfFit();
if (fabs(goodnessOfFit) <= DBL_EPSILON || goodnessOfFit == NULL8) goodnessOfFit = 0;
string pointId = outPoint.Id();
os << pointId << "," << inSamp << ","
<< inLine << "," << outSamp << ","
<< outLine << "," << sampErr << ","
<< lineErr << "," << zScoreMin << ","
<< zScoreMax << "," << goodnessOfFit << endl;
}
}
}
PvlGroup pLog("Points");
pLog+=PvlKeyword("Ignored", ignored);
Application::Log(pLog);
PvlGroup mLog("Measures");
mLog+=PvlKeyword("Validated", validated);
mLog+=PvlKeyword("Registered", registered);
mLog+=PvlKeyword("Unregistered", unregistered);
mLog+=PvlKeyword("Unmeasured", unmeasured);
Application::Log(mLog);
// Log Registration Statistics
Pvl arPvl = ar->RegistrationStatistics();
for(int i = 0; i < arPvl.Groups(); i++) {
Application::Log(arPvl.Group(i));
}
// add the auto registration information to print.prt
PvlGroup autoRegTemplate = ar->RegTemplate();
Application::Log(autoRegTemplate);
outNet.Write(ui.GetFilename("TO"));
delete ar;
}
