Add PROJ projection (#5317)

* Rectifying 9.0.0 changes

* Explicitly add proj to ISIS deps

* Fix field output in qview advanced tracking tool

* Initial addition of PROJ into ISIS

* Added changelog entry

* Near finished state of the PROJ plugin

* Exposed using PROJ strings in cam2map

* Fixed radii acquisition in IProj

* Adjusted how radii is acquired from proj and remove equator check from XYRange

* Clean up headers in base projection classes

* Fix reading existing ISIS mapping groups from geotiffs

* Fixed changelog

* Get ellipsoid information from proj crs if no radii in mapping group

* Added IProj tests

* Remove forced planetographic lat conversion in IProj

* Update tests to resemble ISIS equirectangular projection tests

* Update cam2map and IProj to extract necessary parameters and error reporting from proj

* Updated useproj option exclusions

* Initial stub in for useproj example

* Added proj example jpegs

* Fixed cam2map xml spacing

* Finished up cam2map example

* Added thumbnail assets

Author: acpaquette

CHANGELOG.md

@@ -53,6 +53,7 @@ ctest FunctionalTestJigsawApollo to validate this output. [#5710](https://github
 - Added the ability to disable individual options for parameters in ISIS GUIs [#5849](https://github.com/DOI-USGS/ISIS3/pull/5849)
 - Added PAD or SHRINK options to crop for crops that extend beyond the source image [#5843](https://github.com/DOI-USGS/ISIS3/pull/5843)
 - Added std:: namespace for isinf, fixes build errors for some versions of c++
+- Adds PROJ into ISIS, and exposes the capability with a new class called IProj. [#5317](https://github.com/DOI-USGS/ISIS3/pull/5317)
 
 ### Changed
 - Removed Arm dependency on xalan-c, as it does not build for now on conda-forge. This requires turning off doc building on Arm. Also changed some variables to avoid name clashes on Arm with clang 16. [#5802](https://github.com/DOI-USGS/ISIS3/pull/5802)

isis/CMakeLists.txt

@@ -305,6 +305,7 @@ find_package(NN                REQUIRED)
 find_package(OpenCV    3.1.0   REQUIRED)
 find_package(PCL               REQUIRED)
 find_package(Protobuf          REQUIRED CONFIG)
+find_package(PROJ              REQUIRED)
 find_package(Qwt       6       REQUIRED)
 find_package(SuperLU   4.3     REQUIRED)
 find_package(SpiceQL           REQUIRED)

isis/src/base/apps/cam2map/assets/images/CAM2MAP_USEPROJ_GUI.jpeg

@@ -11,6 +11,8 @@
 #include "Target.h"
 #include "TProjection.h"
 
+#include <proj.h>
+
 using namespace std;
 
 namespace Isis {
@@ -31,7 +33,61 @@ namespace Isis {
 
     // Get the map projection file provided by the user
     Pvl userMap;
-    userMap.read(ui.GetFileName("MAP"));
+    if (ui.GetBoolean("USEPROJ")) {
+      PvlGroup mappingGroup("Mapping");
+      mappingGroup.addKeyword(PvlKeyword("ProjectionName", "IProj"));
+      mappingGroup.addKeyword(PvlKeyword("ProjStr", ui.GetAsString("PROJString")));
+      QString projStr = mappingGroup["ProjStr"][0];
+      PJ_CONTEXT *projContext = proj_context_create();
+      proj_log_level(projContext, PJ_LOG_ERROR);
+
+      /* Create a projection. */
+      PJ *outputProj = proj_create(projContext, projStr.toStdString().c_str());
+
+      if (!outputProj) {
+        proj_context_destroy(projContext);
+        QString msg = "Unable to create projection from [" + projStr + "].";
+        throw IException(IException::User, msg, _FILEINFO_);
+      }
+
+      // Read the radii information from the PROJ projection
+      PJ *ellipsoid = proj_get_ellipsoid(projContext, outputProj);
+
+      if (ellipsoid == nullptr) {
+        proj_destroy(outputProj);
+        proj_context_destroy(projContext);
+        QString msg = "Unable to create ellipsoid from [" + projStr + "]";
+        throw IException(IException::User, msg, _FILEINFO_);
+      }
+
+      double equatorialRadius;
+      double polarRadius;
+
+      int res = proj_ellipsoid_get_parameters(projContext, ellipsoid, 
+                                              &equatorialRadius,
+                                              &polarRadius,
+                                              nullptr,
+                                              nullptr);
+
+      proj_destroy(ellipsoid);
+      proj_destroy(outputProj);
+      proj_context_destroy(projContext);
+
+      if (res == 0) {
+        QString msg = "Unable to get ellipsoid information from [" + projStr + "]";
+        throw IException(IException::User, msg, _FILEINFO_);
+      }
+
+      mappingGroup.addKeyword(PvlKeyword("EquitorialRadius", toString(equatorialRadius, 15)));
+      mappingGroup.addKeyword(PvlKeyword("PolarRadius", toString(polarRadius, 15)));
+      mappingGroup.addKeyword(PvlKeyword("LatitudeType", "Planetographic"));
+      mappingGroup.addKeyword(PvlKeyword("LongitudeDirection", "PositiveEast"));
+      mappingGroup.addKeyword(PvlKeyword("LongitudeDomain", "180", "degrees"));
+      userMap.addGroup(mappingGroup);
+    }
+    else {
+      userMap.read(ui.GetFileName("MAP"));
+    }
     PvlGroup &userGrp = userMap.findGroup("Mapping", Pvl::Traverse);
 
     cam2map(&icube, userMap, userGrp, ui, log);