vconrado/simple_geo

Python Client API for Web Services (WFS and WTSS)

Python Client API for Web Services (WFS and WTSS)

Building and installing simple_geo.py from source

1. In the shell, type

  git clone https://github.com/vconrado/simple_geo.py.git
  cd simple_geo.py/src
  pip install .

Using simple_geo.py

Retrieving Features

from simple_geo import simple_geo as sgeo 

# Connecting to servers
b = sgeo(wfs="http://localhost:8080/geoserver-esensing/",
        wtss="http://localhost:7654")

# Retrieving the list of all available features in the service
ft_list = b.list_features()
print(cv_list)

# Retrieving the metadata of a given feature
ft_scheme = b.describe_feature("esensing:focos_bra_2016")
print(cv_scheme)

# Retrieving the collection for a given feature
fc, fc_metadata = b.feature_collection("esensing:focos_bra_2016")

# Retrieving a selected elements for a given feature
fc, fc_metadata = b.feature_collection("esensing:focos_bra_2016",
                                       attributes=("id", "municipio", "timestamp", "regiao"),
                                       within="POLYGON((-49.515628859948507 -19.394602563415745,-48.020567850467053 -19.610579617637825,-48.354439522883652 -21.052347219666608,-49.849500507163917 -20.836369963642884,-49.515628859948507 -19.394602563415745))",
                                       filter=["satelite_referencia='true'", "timestamp>='2016-01-01'",
                                               "timestamp<'2016-02-01'"],
                                       sort_by=("regiao", "municipio"),
                                       max_features=10)
print(fc)

# Retrieving collection length of selected elements for a given feature
fc_len = b.feature_collection_len("esensing:focos_bra_2016",
                                  within="POLYGON((-49.515628859948507 -19.394602563415745,-48.020567850467053 -19.610579617637825,-48.354439522883652 -21.052347219666608,-49.849500507163917 -20.836369963642884,-49.515628859948507 -19.394602563415745))",
                                  filter=["satelite_referencia='true'", "timestamp>='2016-01-01'",
                                          "timestamp<'2016-02-01'"])
print(fc_len)

Retrieving Coverages

from simple_geo import simple_geo as sgeo

# Connecting to servers
b = sgeo(wfs="http://localhost:8080/geoserver-esensing/",
        wtss="http://localhost:7654")

# Retrieving the list of all available coverages in the service
cv_list = b.list_coverages()
print(cv_list)

# Retrieving the metadata of a given coverage
cv_scheme = b.describe_coverage("climatologia")
print(cv_scheme)

# Retrieving the time series for a given location
ts, ts_metadata = b.time_series("climatologia", ("precipitation", "temperature", "humidity"), -12, -54)
print(ts)
print(ts_metadata)