Skip to main content

Libraries and command-line utilities for geospatial data processing/analysis

Project description


# pygeotools Libraries and command line tools for geospatial data processing and analysis

## Features - Resample/warp rasters to common resolution/extent/projection - Many functions to handle rasters with NoData gaps using [NumPy masked arrays]( - Point data coordinate transformations, sampling, and interpolation routines (e.g., arrays of xyz points) - Common raster filtering operations

### Libraries [pygeotools/lib](./pygeotools/lib) - geolib - Coordinate transformations, raster to vector, vector to raster - malib - NumPy Masked Array operations, DEMStack class - warplib - On-the-fly GDAL warp operations for abitrary number of input datasets - iolib - File input/output, wrappers for GDAL I/O, masked array write to disk - timelib - Time conversions, extract timestamps from filenames, useful for raster time series analysis - filtlib - Collection of filters for 2D masked arrays (Gauss, rolling median, high pass, etc.)

### Command-line utilities (run with no arguments for usage) - - Warp arbitrary rasters to common res/extent/proj - - Create a “stack” of input rasters (a raster time series object) and compute stats - - Clip and mask an input raster using a polygon shapefile - - Apply mask from one raster to another - - Apply various filters available in filtlib - - Remove rows/cols containing only NoData from raster margins - - Replace NoData value - - Automatically determine projection for input lat/lon or raster - - Create polygon shapefile of input raster footprints - - Merge shapefiles - …

## Examples

### Warping multiple datasets to common grid, computing difference, writing out ` from pygeotools.lib import iolib, warplib, malib fn1 = 'raster1.tif' fn2 = 'raster2.tif' ds_list = warplib.memwarp_multi_fn([fn1, fn2], res='max', extent='intersection', t_srs='first', r='cubic') r1 = iolib.ds_getma(ds_list[0]) r2 = iolib.ds_getma(ds_list[1]) rdiff = r1 - r2 malib.print_stats(rdiff) out_fn = 'raster_diff.tif' iolib.writeGTiff(rdiff, out_fn, ds_list[0]) ` or, from the command line…

Warp all to match raster1.tif projection with common intersection and largest pixel size: -tr max -te intersection -t_srs first raster1.tif raster2.tif raster3.tif

Create version of raster1.tif that matches resolution, extent, and projection of raster2.tif: -tr raster2.tif -te raster2.tif -t_srs raster2.tif raster1.tif

Reproject and clip to user-defined extent, preserving original resolution of each input raster: -tr source -te ‘439090 5285360 458630 5306450’ -t_srs EPSG:32610 raster1.tif raster2.tif

### Creating a time series “stack” object: ` from pygeotools.lib import malib fn_list = ['20080101_dem.tif', '20090101_dem.tif', '20100101_dem.tif'] s = malib.DEMStack(fn_list, res='min', extent='union') #Stack standard deviation s.stack_std #Stack linear trend s.stack_trend ` or, from the command line… -tr ‘min’ -te ‘union’ 20*.tif

## Documentation

## Installation

Install the latest release from PyPI:

pip install pygeotools

Note: by default, this will deploy executable scripts in /usr/local/bin

### Building from source

Clone the repository and install:

git clone pip install -e pygeotools

The -e flag (“editable mode”, setuptools “develop mode”) will allow you to modify source code and immediately see changes.

### Core requirements - [GDAL/OGR]( - [NumPy]( - [SciPy](

### Optional requirements (needed for some functionality) - [matplotlib]( - [NASA Ames Stereo Pipeline (ASP)](

## Disclaimer

This originated as a personal repo that I am slowly cleaning up and distributing. There are some useful things that work very well, other things that were hastily written for a one-off task several years ago, and some confusing things that were never finished.

Contributions, bug reports, and general feedback are all welcome. My time is limited, I have some bad habits, and I could really use some help. Thanks in advance.

This was all originally developed for Python 2.X, but should now also work with Python 3.X thanks to [@dlilien](

Some of this functionality now exists in the excellent, mature, well-supported [rasterio]( Eventually, I will integrate rasterio API calls where appropriate.

## License

This project is licensed under the terms of the MIT License.

Project details

Download files

Download the file for your platform. If you're not sure which to choose, learn more about installing packages.

Source Distribution

pygeotools-1.1.1.tar.gz (98.3 kB view hashes)

Uploaded Source

Supported by

AWS AWS Cloud computing and Security Sponsor Datadog Datadog Monitoring Fastly Fastly CDN Google Google Download Analytics Microsoft Microsoft PSF Sponsor Pingdom Pingdom Monitoring Sentry Sentry Error logging StatusPage StatusPage Status page