makesurface 0.0.1dev

Create vector datasets from raster surfaces

A few basic surface creation routines using (primarily) rasterio, fiona, and shapely

Installation

1. Install dependencies:

• rasterio

• fiona

• shapely

• scipy

2. Clone this repo:

git clone [e-mail]:mapbox/make-surface.git

3. CD into make-surface then pip install -e .

Alternatively, install straight from GitHub:

pip install -e [e-mail]:mapbox/make-surface.git#egg=makesurface

Usage

Takes an input raster, and converts into a stacked shapefile. Sort of like gdal polygonize with more control. Also prints out a cartocss template for stylizing (one style for each class).

Turns this:

Into this:

makesurface <input raster> <output shapefile> [OPTIONS]

Arguments

• <input raster> Input single-band raster to class and vectorize

• <output shapefile> Output shapefile to create - output will be “stacked”, with lower values including areas of higher values

Options

• --classfile <TEXT> - Filepath to one-line CSV of manually defined break values;

• --classes <INTEGER> (default = 10) - Number of sections to vectorize into, OR ‘all’ to vectorize into rounded input values; Ignored if classfile specified;

• --weight <FLOAT> (default = 1) - parameter to weight classification type; 0 = quantile breaks, 1 = equal interval, anywhere between 0 and 1 weights the classification between the two; Ignored if classfile specified;

• --smoothing <FLOAT> (default = no smoothing) - optional parameter that defines the upsampling and smoothing prior to classification; If none, no upsampling / smoothing is performed;

• --nodata <ANY NUMBER OR "min"> (default = nodata from metadata) - Manually defined nodata value - can be any number or “min”;

• --carto (default = none) - Flag to include stdout printing of cartoCSS for each class;

• --grib2 (default = none) - Flag to handle .grib2 files with bounds of 0 to 360;

• --help Show this message and exit;