| {"features": [{"geometry": {"coordinates": [[[-111.73527526855469, 41.995094299316406], [-111.65931701660156, 41.99627685546875], [-111.6587142944336, 41.9921875], [-111.65888977050781, 41.95676803588867], [-111.67082977294922, 41.91230010986328], [-111.67332458496094, 41.905494689941406], [-111.67088317871094, 41.90049362182617], [-111.66474914550781, 41.893211364746094], [-111.6506576538086, 41.875465393066406], [-111.64759826660156, 41.87091827392578], [-111.64640808105469, 41.86273956298828], [-111.64334869384766, 41.858192443847656], [-111.63720703125, 41.85499572753906], [-111.633544921875, 41.847267150878906], [-111.63053894042969, 41.83409118652344], [-111.6330337524414, 41.82728576660156], [-111.63983154296875, 41.8227653503418], [-111.6484603881836, 41.82188034057617], [-111.66077423095703, 41.82327651977539], [-111.6712417602539, 41.82330322265625], [-111.67618560791016, 41.82013702392578], [-111.68803405761719, 41.78792953491211], [-111.69361114501953, 41.77931594848633], [-111.70162200927734, 41. |
| >>> from shapely.ops import unary_union | |
| >>> from shapely.geometry import shape | |
| >>> unary_union([shape(g) for g, v in shapes(red, mask=(red < 250))]) | |
| <shapely.geometry.multipolygon.MultiPolygon object at 0x107c1b5d0> | |
| >>> unary_union([shape(g) for g, v in shapes(red, mask=(red < 10))]) | |
| <shapely.geometry.multipolygon.MultiPolygon object at 0x105844310> | |
| >>> import shapely | |
| >>> print shapely.geos.geos_version | |
| (3, 4, 2) |
| $ rio insp ~/Desktop/logo.png Rasterio 0.17.1 Interactive Inspector (Python 2.7.9) | |
| Type "src.meta", "src.read_band(1)", or "help(src)" for more information. | |
| >>> red = src.read(1) | |
| >>> red | |
| array([[255, 255, 255, ..., 255, 255, 255], | |
| [255, 255, 255, ..., 255, 255, 255], | |
| [255, 255, 255, ..., 255, 255, 255], | |
| ..., | |
| [255, 255, 255, ..., 255, 255, 255], | |
| [255, 255, 255, ..., 255, 255, 255], |
| import numpy | |
| import rasterio | |
| from rasterio import Affine as A | |
| from rasterio.warp import reproject, RESAMPLING | |
| with rasterio.open('rasterio/tests/data/RGB.byte.tif') as src: | |
| src_transform = src.affine | |
| # Zoom out by a factor of 2 from the center of the source | |
| # dataset. The destination transform is the product of the |
| # Rasterio is a C extension, need to apt-get gcc and g++. | |
| sudo apt-get update | |
| sudo apt-get install curl gcc g++ | |
| # Download miniconda. | |
| curl "http://mapbox-cloudless-testing.s3.amazonaws.com/landsat8/setup/Miniconda-3.7.0-Linux-x86_64.sh" -o "Miniconda-3.7.0-Linux-x86_64.sh" | |
| chmod +x Miniconda-3.7.0-Linux-x86_64.sh | |
| ./Miniconda-3.7.0-Linux-x86_64.sh -b | |
| # Put conda on path. |
A Rasterio dataset's index() method gives you the row, col index of the pixel containing the point x, y (in the dataset's CRS units).
Define a 1x1 pixel read window starting at that index (see the tiny_window() function below) and use it to read an ndarray. The one in this case has shape (3, 1, 1): 3 bands, 1 row, and 1 column. The values of the 3 bands at the point x, y for this case are 28, 29, and 27.
| import fiona | |
| with fiona.open('-') as src: | |
| print src.bounds |
OS X 10.7+ Python 2.7 (2.7.9 would be best because it includes pip, a package installer)