| # Often when reading in a ShapeFile from Basemap, you'll get: "ValueError: readshapefile can only handle 2D shape types" | |
| # A trick can be to convert your geometry in your GeoPandas Dataframe and restoring the new flattened 2D geometry | |
| # series back into a shapefile and try again. | |
| # edit from http://stackoverflow.com/questions/33417764/basemap-readshapefile-valueerror | |
| from shapely.geometry import Polygon, MultiPolygon, shape, Point | |
| import geopandas as gp | |
| def convert_3D_2D(geometry): | |
| ''' |
This is for an update to the the SpatioTemporal Task View.
If you have any other suggestions, please leave a comment below - or raise an issue here: https://github.com/mdsumner/ctv-mdsumner/issues (or just email/tweet me).
(This is not the final presentation format).
| #!/usr/bin/env python | |
| # ndvi.py red.tif nir.tif output-ndvi.tif | |
| # Calculate NDVI (see Wikipedia). Assumes atmospheric correction. | |
| # (Although I use it without all the time for quick experiments.) | |
| import numpy as np | |
| from sys import argv | |
| from osgeo import gdal, gdalconst |