rasterio.transform.from_origin(west, north, xsize, ysize) Return an Affine transformation given upper left and pixel sizes. https://rasterio.readthedocs.io/en/latest/api/rasterio.transform.html#raster

You can also use json: import json poly = json.loads(obj.job_loc.geojson)

rasterio.mask.mask returns a two element tuple - element 0 is the array, element 1 is the transform. Try something like (untested): # Mask Raster masked_raster, transform = rasterio.mask.mask(file, po

GDALs Sentinel2 driver exposes the data as subdatasets, with (sort of) one for each resolution - https://gdal.org/drivers/raster/sentinel2.html#level-2a To access the subdataset that contains the 10m

Have you tried keeping your windows all the same size (including the edge windows) and specifying `boundless=True` and a `fill_value=some appropriate value` to fill out your windows in your ` `dataset

rasterio.transform.array_bounds should do the trick - https://rasterio.readthedocs.io/en/latest/api/rasterio.transform.html#rasterio.transform.array_bounds

Try going the other way - multiply the resolution by a scale factor (i.e. 3 in your case) and divide the dimensions - see example gist

No. Exactly as I show in my code example. No need to write it out again, as you are just updating the georeferencing with the dataset opened in update mode ("r+"). The dataset is saved automatically w

To override the projection with rasterio: import rasterio from rasterio.crs import CRS in_tif = 'path/to.tif' with rasterio.open(in_tif, 'r+') as src: src.crs = CRS.from_epsg(32629)

Is your matrix of values a numpy array? You can easily create a GeoTIFF in rasterio if you have an array of data and the basic metadata (georeferencing, dimensions, band count, data type). https://ras

How did you install rasterio? If you installed via Anaconda / conda using the default channel, your gdal may lack bigtiff support. https://github.com/ContinuumIO/anaconda-issues/issues/9887

You could use a context manager to clean it up automatically: from contextlib import contextmanager import rasterio from rasterio import MemoryFile @contextmanager def mem_raster(data, **profile): wit

I've put a couple of methods in my answer to your SO question, by expanding the block windows and also for arbitrary windows (unrelated to internal blocks/tiles).

I stacked some singleband rasters into a multiband raster using `rio stack` with a compress=deflate creation option. The compressed output filesize was approx 1.4x the sum of the uncompressed input fi

My GDAL_DATA env var wasn't set. Explanation: I was running the code within the PyCharm IDE that supposedly supports conda envs but doesn't actually *activate* them properly. There are init scripts in

Hi all. Congrats on the 1.0.0 release! Using rasterio 1.0.0 (installed from conda-forge in a fresh conda env). I'm having trouble copying the CRS from on dataset to another. The output dataset CRS is