Geospatial Data

To use this notebook you must have a (free) FEO account, available here. You may then install the feo-client library, if necessary, and then use the package to log in. You can also run notebook in a Colab notebook via the badge below.

In [ ]:

%pip install "feo-client[geo]"

%pip install "feo-client[geo]"

In [ ]:

skip-execution

from feo.client.auth import login

login()

from feo.client.auth import login login()

Geospatial Data¶

Object-oriented interface for Vector Features and Geometries via the Python client

FEO comprises a wide variety of geospatial data, covering vector and raster data types. Vector types include node geometries, hydrobasins, exclusive economic zones and protected areas. Each of these exists as a Collection that can be browsed and searched via several methods.

In this notebook we will gather some vector data for a node, namely:

• The node geometry
• The protected areas intersecting the node geometry

NOTE: To take full advantage of geospatial support, please ensure you install the optional geo dependencies. The install command at the top of this notebook already includes these.

In [ ]:

from feo.client import Features, Node

from feo.client import Features, Node

In [ ]:

# Let's use Bali
node_name = "Bali"
nodes = Node.search(node_name)
try:
    bali = nodes[0]
    print(bali)
except IndexError:
    raise IndexError(f"No nodes found for {node_name}")

# Let's use Bali node_name = "Bali" nodes = Node.search(node_name) try: bali = nodes[0] print(bali) except IndexError: raise IndexError(f"No nodes found for {node_name}")

For a more detailed exploration of Nodes see the Nodes notebook.

In [ ]:

# Access the node's geometry
bali_geom = bali.geometry

# Access the node's geometry bali_geom = bali.geometry

In [ ]:

type(bali_geom)

type(bali_geom)

The geometry property is a Geometry class. This class exposes the geometry type and its (longitude, latitude) coordinates as properties.

In [ ]:

# The Bali geometry is a MultiPolygon; it comprises multiply polygon geometries
bali_geom.type

# The Bali geometry is a MultiPolygon; it comprises multiply polygon geometries bali_geom.type

In [ ]:

# The coordinates list is a nested list of coordinate pairs that define the geometry
bali_geom.coordinates[0]

# The coordinates list is a nested list of coordinate pairs that define the geometry bali_geom.coordinates[0]

The Geometry class proivides convenience methods for downstream purposes; you can convert it to a dictionary, GeoJSON or shapely geometry representation.

In [ ]:

# E.g. convert to GeoJSON string
print(bali_geom.to_geojson()[:100], "...}")  # sliced to limit output

# E.g. convert to GeoJSON string print(bali_geom.to_geojson()[:100], "...}") # sliced to limit output

A shapely geometry provides access to some core geospatial operations and, in notebooks, can be rendered as an image.

In [ ]:

# Create shapely geometry
bali_geom_shply = bali_geom.to_shape()
centroid = bali_geom_shply.centroid
print(
    "Bali centroid:", f"\tLON: {centroid.x:0.6f}", f"\tLAT: {centroid.y:0.6f}", sep="\n"
)

# Create shapely geometry bali_geom_shply = bali_geom.to_shape() centroid = bali_geom_shply.centroid print( "Bali centroid:", f"\tLON: {centroid.x:0.6f}", f"\tLAT: {centroid.y:0.6f}", sep="\n" )

In [ ]:

bali_geom_shply

bali_geom_shply

Search for other vector data¶

FEO provides access to a range of supplimentary geospatial datasets. In this example we search for Protected Areas within Bali. The FEO catalogue will continue to grow!

In [ ]:

# Search for Protected Areas
pa_feats = Features.search(collection_id="biodiversity-wdpa", geometry=bali_geom)

# Search for Protected Areas pa_feats = Features.search(collection_id="biodiversity-wdpa", geometry=bali_geom)

pa_feats is a Features objects; a collection of features, each with properties and a geometry

In [ ]:

# How many protected areas were returned?
len(pa_feats.features)

# How many protected areas were returned? len(pa_feats.features)

In [ ]:

# Examine a feature
ft = pa_feats.features[0]
ft.properties

# Examine a feature ft = pa_feats.features[0] ft.properties

In [ ]:

# Render its geometry
ft.geometry.to_shape()

# Render its geometry ft.geometry.to_shape()

If you have the geo dependencies installed you can convert the Features object to a GeoDataFrame

In [ ]:

gdf = pa_feats.to_geodataframe()
gdf

gdf = pa_feats.to_geodataframe() gdf

Visualise the data¶

Features and Geometry objects both have GeoJSON representations, allowing for easy interoperability with other geospatial tools. In this example, we use folium to plot an interactive map of Bali and its protected areas.

In [ ]:

%pip install folium

%pip install folium

In [ ]:

import folium
from folium.features import GeoJson

import folium from folium.features import GeoJson

In [ ]:

# Define the layer styles
bali_style = {"fillColor": "#f2ff00", "color": "#f2ff00"}
pa_style = {"fillColor": "#000000", "color": "#000000"}

# Define the layer styles bali_style = {"fillColor": "#f2ff00", "color": "#f2ff00"} pa_style = {"fillColor": "#000000", "color": "#000000"}

In [ ]:

bali_geom_shply = (
    bali_geom.to_shape()
)  # convert our Geometry object to a shapely object
map_center = bali_geom_shply.centroid  # access the centroid using the shape property

# Configure the map
map = folium.Map(location=(map_center.y, map_center.x), zoom_start=10)

# Add vector data to the map via the GeoJSON interface
map.add_child(
    GeoJson(bali_geom.to_geojson(), name="Bali", style_function=lambda x: bali_style)
)
map.add_child(
    GeoJson(
        pa_feats.to_geojson(),
        name="World Database on Protected Areas",
        style_function=lambda x: pa_style,
        popup=folium.GeoJsonPopup(fields=["name", "designation"]),
    )
)
folium.LayerControl().add_to(map)
map

bali_geom_shply = ( bali_geom.to_shape() ) # convert our Geometry object to a shapely object map_center = bali_geom_shply.centroid # access the centroid using the shape property # Configure the map map = folium.Map(location=(map_center.y, map_center.x), zoom_start=10) # Add vector data to the map via the GeoJSON interface map.add_child( GeoJson(bali_geom.to_geojson(), name="Bali", style_function=lambda x: bali_style) ) map.add_child( GeoJson( pa_feats.to_geojson(), name="World Database on Protected Areas", style_function=lambda x: pa_style, popup=folium.GeoJsonPopup(fields=["name", "designation"]), ) ) folium.LayerControl().add_to(map) map

In [ ]: