Skip to content

Instantly share code, notes, and snippets.

Last active April 16, 2024 21:39
Show Gist options
  • Save sgillies/7e5cd548110a5b4d45ac1a1d93cb17a3 to your computer and use it in GitHub Desktop.
Save sgillies/7e5cd548110a5b4d45ac1a1d93cb17a3 to your computer and use it in GitHub Desktop.
Advanced Rasterio features notebook
Display the source blob
Display the rendered blob
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Copy link

chelm commented Dec 18, 2017

Great write up @sgillies - especially cool is the lazy warping

Copy link

Seagor commented Dec 19, 2017


Copy link


Copy link

lbustelo commented Mar 7, 2018

I'm trying to run rasterio in my notebook setup and keep getting when I import the module. Any pointers? I'm running my notebooks inside a container built based on the official jupyter docker stacks. I basically installing the python-numpy gdal-bin libgdal-dev on top of that image.

ImportError                               Traceback (most recent call last)
<ipython-input-1-34708b297916> in <module>()
----> 1 import rasterio

/opt/conda/lib/python3.5/site-packages/rasterio/ in <module>()
     13 import warnings
---> 15 from rasterio._base import (
     16     eval_window, window_shape, window_index, gdal_version)
     17 from rasterio.dtypes import (

ImportError: /usr/lib/ undefined symbol: sqlite3_column_table_name

Copy link

lbustelo commented Mar 7, 2018

Never mind the above... I had and old version of GDAL. After that, I pulled the --pre version of rasterio and was able to get going.

Copy link

renelikestacos commented Jul 6, 2018

by some reason i am getting error on from rasterio.vrt import WarpedVRT it shows me an import error.
ImportError: No module named 'rasterio.vrt'
I checked stackoverflow but couldn't find a solution by now. Anyone has the same issue?

Copy link

Thanks for sharing this wonderful write up. I noticed that AWS_access_key is required, does that mean the landsat data is requester-pay?

Copy link

mwkorver commented Dec 5, 2018

Copy link

Dehde commented Jan 18, 2019

First of all, thanks for the great introduction!
2 Questions:
1.What about GeoTIFFs that are striped instead of tiled? Is it still possible to access a subcrop efficiently?
2. Unfortunately the "with rasterio.Env(CPL_CURL_VERBOSE=True)" changes nothing in the output for me. I'm working on windows with rasterio.version = 1.0.8. Any ideas how I could achieve the verbose output? Been stuck on that for a bit..

Maybe interesting for other windows users: I got an SSL error when trying to access the landsat images. I got rid of it by exposing my SSL certificate with:
import os
os.environ["CURL_CA_BUNDLE"] = r"C:\path\to\curl-ca-bundle.crt"

Copy link

This is fantastic. I have been attempting to figure out how to manipulate LANDSAT and Sentinel data from AWS without downloading wholesale and without using packages like eo-learn and GEE (too hard to customize for my liking). Had seen these features in rasterio but never knew how they worked. This is a well-done write up, cheers.


Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment