I created a circle using geopandas and it returned a shapely polygon:
POLYGON: ((...))
I want this same polygon as a geojson object. I ran across this:
shapely.geometry.mapping(shapelyObject)
which returns this:
{'type': 'Polygon', 'coordinates': (((570909.9247264927, 125477.71811034005)...}
But when I try to map this in mapbox it does not show anything. I think maybe it is not fully a geojson object.
If you don't want to create this dict manually, you can also rely on geopandas
creating it:
In [1]: import shapely.geometry
In [2]: import geopandas
In [3]: shapely_polygon = shapely.geometry.Polygon([(0, 0), (0, 1), (1, 0)])
In [4]: geopandas.GeoSeries([shapely_polygon]).__geo_interface__
Out[4]:
{'bbox': (0.0, 0.0, 1.0, 1.0),
'features': [{'bbox': (0.0, 0.0, 1.0, 1.0),
'geometry': {'coordinates': (((0.0, 0.0),
(0.0, 1.0),
(1.0, 0.0),
(0.0, 0.0)),),
'type': 'Polygon'},
'id': '0',
'properties': {},
'type': 'Feature'}],
'type': 'FeatureCollection'}
(Note that this gives a FeatureCollection and not a single feature.)
Or to a string (or file):
In [4]: geopandas.GeoSeries([shapely_polygon]).to_json()
Out[4]: '{"features": [{"bbox": [0.0, 0.0, 1.0, 1.0], "geometry": {"coordinates": [[[0.0, 0.0], [0.0, 1.0], [1.0, 0.0], [0.0, 0.0]]], "type": "Polygon"}, "properties": {}, "id": "0", "type": "Feature"}], "bbox": [0.0, 0.0, 1.0, 1.0], "type": "FeatureCollection"}'
Something like this should do the trick:
features = [{'type': 'Feature', 'properties': {}, 'geometry': shapely.geometry.mapping(shapelyObject)}]
Now you can try to map features
in mapbox.
Hope this helps.
Reference:
https://gis.stackexchange.com/questions/213717/geometry-workflow-from-shapely-to-geojson
To write a standard geojson object using pandas you shall use the driver provided by fiona
as recommended in the documentation
gdf.to_file('path/to/file.geojson', driver='GeoJSON')
See import fiona; fiona.supported_drivers
for a list of fully supported drivers
Shapely returns a python dict where all the coordinates are in tuples. You need to convert to JSON in order for mapbox, etc... to properly accept it.
json.dumps(shapely.geometry.mapping(shapelyObject))
Use the driver provided by fiona:
data=shapefile.to_file("file.geojson",driver='GeoJSON')
data=geopandas.read_file("file.geojson")
data