pygeoutils 0.1.4

A set of utilities for manipulating (Geo)JSON and GeoTIFF data.

Features

PyGeoUtils is a part of Hydrodata software stack and provides utilities for manipulating (Geo)JSON and GeoTIFF data:

• json2geodf: For converting (Geo)JSON objects to GroPandas dataframe.

• arcgis2geojson: For converting ESRIGeoJSON objects to standard GeoJSON format.

• gtiff2xarray: For converting (Geo)TIFF objects to xarray datasets.

All these function handle all necessary CRS transformations. Moreover, requests for additional functionalities can be submitted via issue tracker.

Installation

You can install pygeoutils using pip after installing libgdal on your system (for example, in Ubuntu run sudo apt install libgdal-dev):

$ pip install pygeoutils

Alternatively, pygeoutils can be installed from the conda-forge repository using Conda:

$ conda install -c conda-forge pygeoutils

Quickstart

To demonstrate capabilities of PyGeoUtils lets use PyGeoOGC to access National Wetlands Inventory from WMS, and FEMA National Flood Hazard via WFS, then convert the outpus to GeoDataFrame and xarray.Dataset using PyGeoUtils.

import pygeoutils as geoutils
from pygeoogc import WFS, WMS
from shapely.geometry import Polygon


geometry =  Polygon(
    [
        [-118.72, 34.118],
        [-118.31, 34.118],
        [-118.31, 34.518],
        [-118.72, 34.518],
        [-118.72, 34.118],
    ]
)

url_wms = "https://www.fws.gov/wetlands/arcgis/services/Wetlands_Raster/ImageServer/WMSServer"
wms = WMS(
    url_wms,
    layers="0",
    outformat="image/tiff",
    crs="epsg:3857",
)
r_dict = wms.getmap_bybox(
    geometry.bounds,
    1e3,
    box_crs="epsg:4326",
)
wetlands = geoutils.gtiff2xarray(r_dict, geometry, "epsg:4326")

url_wfs = "https://hazards.fema.gov/gis/nfhl/services/public/NFHL/MapServer/WFSServer"
wfs = WFS(
    url_wfs,
    layer="public_NFHL:Base_Flood_Elevations",
    outformat="esrigeojson",
    crs="epsg:4269",
)
r = wfs.getfeature_bybox(geometry.bounds, box_crs="epsg:4326")
flood = geoutils.json2geodf(r.json(), "epsg:4269", "epsg:4326")

Contributing

Contirbutions are very welcomed. Please read CODE_OF_CONDUCT.rst and CONTRIBUTING.rst files for instructions.