HARMONY-1859: Add WKT POINT support in harmony-py.

examples/tutorial.ipynb

@@ -325,7 +325,7 @@
    "id": "a22a34f6-36a5-4554-b248-1b3e599dc4c2",
    "metadata": {},
    "source": [
-    "Example of submitting a request with WKT spatial:"
+    "Example of submitting a request with WKT spatial. The supported WKT geometry types are listed at: https://harmony-py.readthedocs.io/en/latest/api.html#harmony.harmony.WKT"
    ]
   },
   {
@@ -384,6 +384,14 @@
     "    if filename.endswith(\"png\"):\n",
     "        helper.show_result(filename)"
    ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "09f132f3-6117-4ca0-b2b6-ae41c6b48981",
+   "metadata": {},
+   "outputs": [],
+   "source": []
   }
  ],
  "metadata": {

harmony/harmony.py

@@ -151,7 +151,18 @@ def __repr__(self) -> str:
 
 
 class WKT:
-    """The Well Known Text (WKT) representation of Spatial."""
+    """The Well Known Text (WKT) representation of Spatial.
+    Supported WKT geometry types are: POINT, MULTIPOINT, POLYGON, MULTIPOLYGON.
+
+    Example:
+        spatial=WKT('POINT(-40 10)')
+
+        spatial=WKT('MULTIPOINT((-77 38.9),(-40 10))')
+
+        spatial=WKT('POLYGON((-140 20, -50 20, -50 60, -140 60, -140 20))')
+
+        spatial=WKT('MULTIPOLYGON(((10 10, 20 20, 30 10, 10 10)),((40 40, 50 50, 60 40, 40 40)))')
+    """
 
     def __init__(self, wkt: str):
         """Constructs a WKT instance of spatial area.
@@ -256,7 +267,8 @@ class Request(BaseRequest):
     Args:
         collection: The CMR collection that should be queried
 
-        spatial: Bounding box spatial constraints on the data
+        spatial: Bounding box spatial constraints on the data or Well Known Text (WKT) string
+          describing the spatial constraints.
 
         temporal: Date/time constraints on the data provided as a dict mapping "start" and "stop"
           keys to corresponding start/stop datetime.datetime objects
@@ -607,6 +619,20 @@ def _http_method(self, request: BaseRequest) -> str:
         method = 'GET' if isinstance(request, CapabilitiesRequest) else 'POST'
         return method
 
+    def _wkt_to_edr_route(self, wkt_string: str) -> str:
+        """Returns the EDR route for the given WKT string."""
+        # Load the WKT string into a Shapely geometry object
+        geometry = loads(wkt_string)
+
+        if geometry.geom_type == 'Polygon' or geometry.geom_type == 'MultiPolygon':
+            return 'area'
+        elif geometry.geom_type == 'Point' or geometry.geom_type == 'MultiPoint':
+            return 'position'
+        elif geometry.geom_type == 'LineString' or geometry.geom_type == 'MultiLineString':
+            return 'trajectory'
+        else:
+            raise Exception(f"Unsupported geometry type: {geometry.geom_type}")
+
     def _submit_url(self, request: BaseRequest) -> str:
         """Constructs the URL for the request that is used to submit a new Harmony Job."""
         if isinstance(request, CapabilitiesRequest):
@@ -616,7 +642,7 @@ def _submit_url(self, request: BaseRequest) -> str:
                 f'{self.config.root_url}'
                 f'/ogc-api-edr/1.1.0/collections'
                 f'/{request.collection.id}'
-                f'/area'
+                f'/{self._wkt_to_edr_route(request.spatial.wkt)}'
             )
         else:
             return (