Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

update to socetframesettings #8

Closed
thareUSGS opened this issue Feb 6, 2018 · 5 comments
Closed

update to socetframesettings #8

thareUSGS opened this issue Feb 6, 2018 · 5 comments
Assignees

Comments

@thareUSGS
Copy link

thareUSGS commented Feb 6, 2018

During the write of the ISD (keywords.lis) file we need to add two more keywords. These make sense to be towards the top of the file - perhaps just under "SENSOR_TYPE" and "SENDOR_MODE". The two keywords are:

SEMI_MAJOR_AXIS
ECCENTRICITY

Perhaps some pseudo code. I have no idea if this is the correct method to grab equatorial and polar radii from a level 0 or 1 cube.

    if (mapping.hasKeyword("EquatorialRadius") && mapping.hasKeyword("PolarRadius")) {
      m_equatorialRadius = new Distance(toDouble(mapping["EquatorialRadius"][0]),
                                        Distance::Meters);
      m_polarRadius = new Distance(toDouble(mapping["PolarRadius"][0]),
                                   Distance::Meters);
    }
    else {
      try {
        PvlGroup radiiGrp = Target::radiiGroup(mapping["TargetName"][0]);
        m_equatorialRadius = new Distance(toDouble(radiiGrp["EquatorialRadius"][0]),
                                          Distance::Meters);
        m_polarRadius = new Distance(toDouble(radiiGrp["PolarRadius"][0]),
                                     Distance::Meters);
      }
      catch (IException &e) {
        QString msg = "Unable to grab radii from cube.";
        throw IException(e, IException::Unknown, msg, _FILEINFO_);
      }
    }

     m_eccentricity = 1.0 -
                      (m_polarRadius * m_polarRadius) /
                      (m_equatorialRadius * m_equatorialRadius);
     m_eccentricity = sqrt(m_eccentricity);

so:
SEMI_MAJOR_AXIS m_equatorialRadius
ECCENTRICITY m_eccentricity

e.g. Mars
RECTIFICATION_TERMS
3.00000000000000e+04 0.00000000000000e+00 1.00000000000000e+00
1.00000000000000e+04 1.00000000000000e+00 0.00000000000000e+00
GROUND_ZERO 4.27444857125941e-01 -3.26805190989393e-01 0.00000000000000e+00
LOAD_PT 4.27444857125941e-01 -3.26805190989393e-01 0.00000000000000e+00
COORD_SYSTEM 1
IMAGE_MOTION 0
SENSOR_TYPE USGSAstroLineScanner
SENSOR_MODE UNKNOWN
SEMI_MAJOR_AXIS 3.39619000000000e+006
ECCENTRICITY 1.08339143554195e-001
...

@jessemapel
Copy link
Contributor

@thareUSGS Can you make a Redmine ticket for this and #9 in ISIS

@jessemapel jessemapel reopened this Mar 23, 2018
@jessemapel jessemapel self-assigned this Mar 23, 2018
@jessemapel
Copy link
Contributor

@thareUSGS
Copy link
Author

added redmine ticket: https://isis.astrogeology.usgs.gov/fixit/issues/5363

@thareUSGS
Copy link
Author

changes were made and can be viewed but the frame import in Socet Set or GXP will probably not use these keys. While it shouldn't hurt, they are deemed unnecessary for now for the framing cameras.

@jlaura
Copy link
Collaborator

jlaura commented Sep 5, 2018

Managed in redmine, external to this repo.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

3 participants