search

DOI-USGS / ISIS3

Integrated Software for Imagers and Spectrometers v3. ISIS3 is a digital image processing software package to manipulate imagery collected by current and past NASA and International planetary missions.
https://isis.astrogeology.usgs.gov
Other
198 stars 167 forks source link

"Observation Mode" in jigsaw is not working. #1392

Closed ascbot closed 5 years ago

ascbot commented 5 years ago

Author Name: Kenneth Edmundson (Kenneth Edmundson)

Original Assignee: Kenneth Edmundson

Reported by Megan Henriksen of ASU. Problem occurred when processing LRO NAC image left/right pairs as single observations. Megan's email is below...

On Wed, Jun 13, 2012 at 3:49 PM, Megan Henriksen megan.henriksen@asu.edu wrote:

Hi,

Thank you so much for your help last time.

I have been regularly encountering an error when I try to use jigsaw with observations=yes on several different mosaics, although they will converge with varying degrees of success with observations=no. I have been using isis 3.4.0 and checking the control networks with cnetcheck for islands, etc. I don't know if it is relevant, but I have also been using the new LROC NAC test kernels for spiceinit.

This is an example of the command I have been trying to use, although the bundle adjustments converge with a variety of camsolve, spsolve, and sigma parameters and with outlier_rejection on or off:
jigsaw fromlist=phil_level1_cubes.lis cnet=phil_qnet.cnet onet=phil_jigsaw.net observations=yes method=sparse outlier_rejection=no camsolve=angles twist=yes spsolve=position spacecraft_position_sigma=100 camera_angles_sigma=0.05 file_prefix=phil_

Jigsaw does this:
jigsaw: Loading Control Points...
100% Processed
jigsaw: Setting input images...
100% Processed
Validating network...
starting iteration 1

matrix block retrieval failure at column 3, row 3Total # of block columns: 6Total # of blocks: 8Object = jigsaw

then it fails with this error:
Group = Error
    Program = jigsaw
    Class   = "USER ERROR"
    Code    = 2
    Message = "Unable to bundle adjust network [ap17_tl_qnet.net]"
    File    = jigsaw.cpp
    Line    = 232
  End_Group

  Group = Error
    Program = jigsaw
    Class   = "PROGRAMMER ERROR"
    Code    = 3
    Message = "CHOLMOD: Failed to load Triplet matrix"
    File    = BundleAdjust.cpp
    Line    = 2256
  End_Group

Emerson Speyerer informed me that some of you have access to /ser/estore/; if you are interested in seeing my data it is located here:
/ser/estore/lroc/work/scratch/mhenriksen/2012-06-07_copernicus_controlled_mosaic
/ser/estore/lroc/work/scratch/mhenriksen/2012-06-07_philolaus_controlled_mosaic
/ser/estore/lroc/work/scratch/mhenriksen/2012-06-12_Apollo17_Taurus_Littrow

If you can't get to it through estore and you would still like to see it, please let me know. I know some of the LROC people here are going to the workshop there in about two weeks, so I could probably send it with them.

Thank you very much for your help. Treating the nac pairs as a single observation is something we are very interested in doing, so it would be great if you could help me get this to work. If you need anything else, please let me know.

Thank you,
Megan

Steps to reproduce:

Ella created two small LRO NAC test sets. When setting "observation mode" to true, the error message reported by Megan above is reproduced.

ascbot commented 5 years ago

Original Redmine Comment Author Name: Kenneth Edmundson (Kenneth Edmundson) Original Date: 2012-06-22T17:48:55Z

IMPACT The impacted programs in Isis include jigsaw. Anticipated impact to other parts of Isis is none.

ascbot commented 5 years ago

Original Redmine Comment Author Name: Ella (Contractor) Lee (Ella (Contractor) Lee) Original Date: 2012-06-28T23:30:31Z

I ran findimageoverlaps, autoseed, and pointreg on the 2 pairs, so all the images are linked together. I ran jigsaw with "observation=no" and "observation=yes," and got different results. The residual error is higher when I set observation to yes. Note the pairs are offset from each other by more than 200 pixels.

TEST1: observation=no Object = jigsaw IsisVersion = "3.4.1.0 alpha | 2012-05-21" ProgramVersion = 2011-12-20 ProgramPath = /work/projects/isis/latest/m00956/isis/bin ExecutionDateTime = 2012-06-28T12:54:48 HostName = astrovm2.wr.usgs.gov UserName = elee Description = "Improves camera pointing and a whole lot more!"

Group = UserParameters FROMLIST = all.lis CNET = tieall_ptregeditval.net ONET = tieall_ptregedit_jig1.net OBSERVATIONS = No RADIUS = YES UPDATE = No METHOD = SPARSE ERRORPROPAGATION = No OUTLIER_REJECTION = No REJECTION_MULTIPLIER = 3.0 MODEL1 = NONE SIGMA0 = 1.0e-10 MAXITS = 50 CKDEGREE = 2 SOLVEDEGREE = 2 CAMSOLVE = ANGLES TWIST = Yes SPSOLVE = NONE POINT_RADIUS_SIGMA = 100 CAMERA_ANGLES_SIGMA = 1 FILE_PREFIX = jig_no_obser BUNDLEOUT_TXT = yes OUTPUT_CSV = yes RESIDUALS_CSV = yes End_Group

Group = Iteration1 Sigma0 = 0.20972738082402 Observations = 666 Constrained_Point_Parameters = 159 Constrained_Image_Parameters = 12 Unknown_Parameters = 489 Degrees_of_Freedom = 348 Rejected_Measures = 0 End_Group

Group = Iteration2 Sigma0 = 0.17986467835764 Observations = 666 Constrained_Point_Parameters = 159 Constrained_Image_Parameters = 12 Unknown_Parameters = 489 Degrees_of_Freedom = 348 Rejected_Measures = 0 End_Group

Group = Iteration3 Sigma0 = 0.17985146831918 Observations = 666 Constrained_Point_Parameters = 159 Constrained_Image_Parameters = 12 Unknown_Parameters = 489 Degrees_of_Freedom = 348 Rejected_Measures = 0 End_Group

Group = "Iteration4: Final" Sigma0 = 0.17985146831217 Observations = 666 Constrained_Point_Parameters = 159 Constrained_Image_Parameters = 12 Unknown_Parameters = 489 Degrees_of_Freedom = 348 Rejected_Measures = 0 Converged = TRUE TotalElapsedTime = 0.14 End_Group

Group = JigsawResults Status = "Camera pointing NOT updated" End_Group

Group = Accounting ConnectTime = 00:00:01.0 CpuTime = 00:00:00.4 End_Group End_Object

TEST2: Observation=yes

Object = jigsaw IsisVersion = "3.4.1.0 alpha | 2012-05-21" ProgramVersion = 2011-12-20 ProgramPath = /work/projects/isis/latest/m00956/isis/bin ExecutionDateTime = 2012-06-28T12:55:46 HostName = astrovm2.wr.usgs.gov UserName = elee Description = "Improves camera pointing and a whole lot more!"

Group = UserParameters FROMLIST = all.lis CNET = tieall_ptregeditval.net ONET = tieall_ptregedit_jig2.net OBSERVATIONS = YES RADIUS = YES UPDATE = No METHOD = SPARSE ERRORPROPAGATION = No OUTLIER_REJECTION = No REJECTION_MULTIPLIER = 3.0 MODEL1 = NONE SIGMA0 = 1.0e-10 MAXITS = 50 CKDEGREE = 2 SOLVEDEGREE = 2 CAMSOLVE = ANGLES TWIST = Yes SPSOLVE = NONE POINT_RADIUS_SIGMA = 100 CAMERA_ANGLES_SIGMA = 1 FILE_PREFIX = jig2_obser BUNDLEOUT_TXT = yes OUTPUT_CSV = yes RESIDUALS_CSV = yes End_Group

Group = Iteration1 Sigma0 = 6.7663409152862 Observations = 666 Constrained_Point_Parameters = 159 Constrained_Image_Parameters = 6 Unknown_Parameters = 483 Degrees_of_Freedom = 348 Rejected_Measures = 0 End_Group

Group = Iteration2 Sigma0 = 8.2662868680761 Observations = 666 Constrained_Point_Parameters = 159 Constrained_Image_Parameters = 6 Unknown_Parameters = 483 Degrees_of_Freedom = 348 Rejected_Measures = 0 End_Group

Group = Iteration3 Sigma0 = 22.744787939855 Observations = 666 Constrained_Point_Parameters = 159 Constrained_Image_Parameters = 6 Unknown_Parameters = 483 Degrees_of_Freedom = 348 Rejected_Measures = 0 End_Group

Group = Iteration4 Sigma0 = 69.679826814376 Observations = 666 Constrained_Point_Parameters = 159 Constrained_Image_Parameters = 6 Unknown_Parameters = 483 Degrees_of_Freedom = 348 Rejected_Measures = 0 End_Group

Group = Iteration5 Sigma0 = 149.34696098823 Observations = 666 Constrained_Point_Parameters = 159 Constrained_Image_Parameters = 6 Unknown_Parameters = 483 Degrees_of_Freedom = 348 Rejected_Measures = 0 End_Group

Group = Iteration6 Sigma0 = 320.75689557784 Observations = 666 Constrained_Point_Parameters = 159 Constrained_Image_Parameters = 6 Unknown_Parameters = 483 Degrees_of_Freedom = 348 Rejected_Measures = 0 End_Group

Group = Iteration7 Sigma0 = 342.97296131283 Observations = 666 Constrained_Point_Parameters = 159 Constrained_Image_Parameters = 6 Unknown_Parameters = 483 Degrees_of_Freedom = 348 Rejected_Measures = 0 End_Group

Group = Iteration8 Sigma0 = 423.01899164469 Observations = 666 Constrained_Point_Parameters = 159 Constrained_Image_Parameters = 6 Unknown_Parameters = 483 Degrees_of_Freedom = 348 Rejected_Measures = 0 End_Group

Group = Iteration9 Sigma0 = 359.24298872106 Observations = 666 Constrained_Point_Parameters = 159 Constrained_Image_Parameters = 6 Unknown_Parameters = 483 Degrees_of_Freedom = 348 Rejected_Measures = 0 End_Group

Group = Iteration10 Sigma0 = 427.29667033109 Observations = 666 Constrained_Point_Parameters = 159 Constrained_Image_Parameters = 6 Unknown_Parameters = 483 Degrees_of_Freedom = 348 Rejected_Measures = 0 End_Group

Group = Iteration11 Sigma0 = 359.71656664335 Observations = 666 Constrained_Point_Parameters = 159 Constrained_Image_Parameters = 6 Unknown_Parameters = 483 Degrees_of_Freedom = 348 Rejected_Measures = 0 End_Group

Group = Iteration12 Sigma0 = 427.66286164845 Observations = 666 Constrained_Point_Parameters = 159 Constrained_Image_Parameters = 6 Unknown_Parameters = 483 Degrees_of_Freedom = 348 Rejected_Measures = 0 End_Group . . Oscillates between 427.7 and 359.76 . Group = Iteration48 Sigma0 = 427.70067027443 Observations = 666 Constrained_Point_Parameters = 159 Constrained_Image_Parameters = 6 Unknown_Parameters = 483 Degrees_of_Freedom = 348 Rejected_Measures = 0 End_Group

Group = Iteration49 Sigma0 = 359.76460473775 Observations = 666 Constrained_Point_Parameters = 159 Constrained_Image_Parameters = 6 Unknown_Parameters = 483 Degrees_of_Freedom = 348 Rejected_Measures = 0 End_Group

Group = Iteration50 Sigma0 = 427.70023851555 Observations = 666 Constrained_Point_Parameters = 159 Constrained_Image_Parameters = 6 Unknown_Parameters = 483 Degrees_of_Freedom = 348 Rejected_Measures = 0 End_Group

Group = JigsawResults Status = "Camera pointing NOT updated" End_Group

Group = Accounting ConnectTime = 00:00:02.0 CpuTime = 00:00:02.0 End_Group End_Object

My files are in this directory:

/work/users/elee/IsisSupport/m00956_jigsaw/lroc/Tie_between_all_four_images

Ella Mae