Handle generic __geo_interface__ objects via fiona and geopandas

Hacky attempt to handle non-geopandas objects (e.g.
shapely.geometry) that have a __geo_interface__ dictionary
property associated with it. Still assumes that geopandas
is installed (along with fiona), so not a perfect solution.
Also added another test with different geometry types.

Author: weiji14

pygmt/helpers/tempfile.py

@@ -127,6 +127,32 @@ def tempfile_from_geojson(geojson):
     """
     with GMTTempFile(suffix=".gmt") as tmpfile:
         os.remove(tmpfile.name)  # ensure file is deleted first
-        # Using geopandas.to_file to directly export to OGR_GMT format
-        geojson.to_file(filename=tmpfile.name, driver="OGR_GMT", mode="w")
+        ogrgmt_kwargs = dict(filename=tmpfile.name, driver="OGR_GMT", mode="w")
+        try:
+            # Using geopandas.to_file to directly export to OGR_GMT format
+            geojson.to_file(**ogrgmt_kwargs)
+        except AttributeError:
+            # pylint: disable=import-outside-toplevel
+            # Other 'geo' formats which implement __geo_interface__
+            import json
+
+            import fiona
+            import geopandas as gpd
+
+            with fiona.Env():
+                jsontext = json.dumps(geojson.__geo_interface__)
+                # Do Input/Output via Fiona virtual memory
+                with fiona.io.MemoryFile(file_or_bytes=jsontext.encode()) as memfile:
+                    geoseries = gpd.GeoSeries.from_file(filename=memfile)
+                    geoseries.to_file(**ogrgmt_kwargs)
+
+                #     with memfile.open(driver="GeoJSON") as collection:
+                #         # Get schema from GeoJSON
+                #         schema = collection.schema
+                # # Write to temporary OGR_GMT format file
+                # with fiona.open(
+                #     fp=tmpfile.name, mode="w", driver="OGR_GMT", schema=schema
+                # ) as ogrgmtfile:
+                #     ogrgmtfile.write(geojson.__geo_interface__)
+
         yield tmpfile.name

pygmt/tests/test_geopandas.py

@@ -46,3 +46,21 @@ def test_geopandas_info_geodataframe(gdf):
     """
     output = info(table=gdf, per_column=True)
     npt.assert_allclose(actual=output, desired=[0.0, 35.0, 0.0, 20.0])
+
+
+@pytest.mark.parametrize(
+    "geomtype,desired",
+    [
+        ("multipolygon", [0.0, 35.0, 0.0, 20.0]),
+        ("polygon", [20.0, 23.0, 10.0, 14.0]),
+        ("linestring", [20.0, 30.0, 15.0, 15.0]),
+    ],
+)
+def test_geopandas_info_shapely(gdf, geomtype, desired):
+    """
+    Check that info can return the bounding box region from a shapely.geometry
+    object that has a __geo_interface__ property.
+    """
+    geom = gdf.loc[geomtype].geometry
+    output = info(table=geom, per_column=True)
+    npt.assert_allclose(actual=output, desired=desired)