Questions tagged [rasterio]
Clean and fast and geospatial raster I/O for Python programmers who use Numpy
0
votes
1answer
16 views
Rasterio.features.rastersize
I have to transform the list of polygons that is in vector form into a raster mask that should be geographically aligned with the NDVI image that I have. I was advised to use rasterio.features....
0
votes
0answers
21 views
Different values between original raster and cropped raster with rasterio.mask
I have a large sentinel-2 image and I want to clip it according polygons which are stored in a geojson file.
My code is:
df = gpd.read_file(os.path.join(geojson_folder,geojson_file))
df.head()
...
0
votes
0answers
18 views
Would like to randomly sample points from within a polygon using Python [duplicate]
I would like to achieve something similar to the ArcGIS Create Random Points tool using geopandas, shapely, or rasterio. I aim to have a polygon with some holes in it (removed by subtracting buffered ...
0
votes
1answer
29 views
Rasterio resolution warping via CLI
I'm trying to warp a raster to meet an exact resolution (0.0022500000000000007, -0.0022499999999999977). I've tried using the CLI to no avail.
The example in the docs is
rio warp input.tif output....
1
vote
1answer
63 views
How to mosaic images with different CRS with rasterio
I am doing mosaicing of 14 raster images (.img format) to completely cover a province.see original spatial location of separately opened tiles(). Mosiac operation is completed but 4 tiles are mosiaced ...
1
vote
0answers
75 views
Rasterizing shapefile with rasterio and gdal gives all cells 0 values [closed]
I am trying to rasterize my shapefile (20 polygons) with rasterio and fiona and turn the new raster into a mask for clipping.
I need to rasterize my all polygons in my shapefile and turn on mask ...
0
votes
0answers
65 views
How delete delete or mask nodata 0 value in the mask
I have raster (with 5 bands) and shapfile in the same projection.
I do this script for rastersize and do mask and clip raster or crop raster with mask but I need to have just clip raster without black ...
0
votes
0answers
11 views
Rasterio CLI warping multiple rasters like another raster in all minus one feature
I'm writing a bash script to warp all the rasters in a specific directory to match another raster. I'd like to warp all these rasters iteratively in to match the other raster in all was except shape ...
3
votes
0answers
67 views
Writing masked GeoTIFF in rasterio and visualizing in QGIS?
I have an interpolated raster, which is a masked array, that shows where the data is valid. I want to save it as masked TIFF with rasterio.
To test, here's the JSON file, and the code:
import ...
1
vote
0answers
23 views
Reading a RasterIO MemoryFile triggers a GTiff driver warning/error message
I have recently upgraded RasterIO from version 1.0.5 to 1.0.12. After this upgrade I now get CPLE_NotSupported in driver GTiff does not support open option MODE when reading a MemoryFile.
My code to ...
1
vote
2answers
159 views
Rasterio error: Input shapes do not overlap raster but reprojection not working
I'm trying to use rasterio (v1.0.13) and fiona to perform a raster clip on a geotiff using a geojson polygon. For reference, the clip works perfectly from the command line using GDAL:
gdalwarp -...
2
votes
1answer
72 views
Split GeoTIFF into multiple cells with Rasterio
I'd like to create training data for a machine learning pipeline. To do this I want to take a large raster and split it into multiple equal sized cells.
I have a very basic loop for generating the ...
0
votes
1answer
234 views
Reading a 4 band image with rasterio (python)
I am trying to view a tif satellite image which has 4 bands. I want to remove the last band (NIR) and view the RGB image only, so I am trying to split the NIR from the rest of the image. Here is my ...
1
vote
1answer
274 views
How to visualize multiband imagery using rasterio?
How do I plot a satellite image using rasterio? I'm able to plot the individual bands but not the entire image. Is there a way to combine the separate R, G, B numpy arrays to create a standard RGB ...
0
votes
1answer
68 views
Loading all rasters in a directory and naming the variables iteratively
Hi I'm working with prism data and I already have clipped (masked) rasters (*.bil files) of a specific area for each day of a year. I'd like to open all these rasters with rasterio, convert them into ...
1
vote
1answer
124 views
How to convert vector format (geopackage, shapefile) to raster format (geotiff) using rasterio without template raster
I'm looking to create a raster given just a geopackage (or shapefile, or other vector format). I don't have a template raster, and all I wish to do is burn the shapes into a raster using rasterio. ...
3
votes
1answer
59 views
How to use rasterio transform on all bands
I've imported a 10 band raster tif using rasterio as ds. I'm attempting to pad the image but keep getting errors with the transform.
print(ds.transform)
pilpadimg = rasterio.pad(ds.read(), transform=...
3
votes
1answer
88 views
Cannot upscale raster with rasterio
I am trying to upscale a raster using rasterio but I don't know why my process gets killed.
It can print "i am arrived here" but on the reproject() I have as terminal output:
killed
This is the ...
1
vote
1answer
36 views
cannot read raster of 14 MB
I would like to start using pure python to handle geospatial data since POSTGIS is not so optimized for rasters.
I starting using rasterio and with this code :
with rasterio.open(clc_r) as src:
...
3
votes
0answers
116 views
Using Rasterio to work with two rasters of different extent and resolution? [duplicate]
I have two spatial rasters of different extent and different cell size and I want to do a conditional statement with them. I decided to use rasterio for that. I use this code:
with rasterio.open('A....
1
vote
1answer
108 views
Using rasterio and matplotlib how can I create a window for an NDVI image to fit into?
I am using Rasterio, Numpy, and Matplotlib. I created a window with the dimensions of 800x600 and plotted the band 5 for a Landsat 8 image, that worked just fine. I also have the nir image so then I ...
0
votes
1answer
162 views
Using Rasterio, Matplotlib, and Numpy trying to calculate and plot NDVI of a raster image getting error: IndexError: list index out of range?
First I create two variables and store my raster images in them and I use MatPlotLib to plot a subset of the original "nirband" raster image. That works fine. I am getting an IndexError: list index ...
1
vote
1answer
162 views
Masking GeoTIFF file after GeoJSON through rasterio - “Input shapes do not overlap raster”
So my intention is quite simple:
I have a 4-Band TIFF file I want to crop after a polygon shape. However, I'm getting the following two exceptions
WindowError: windows do not intersect
ValueError: ...
1
vote
1answer
73 views
Reprojection of satellite images mosaic to CRS with latitude and longitudes
I downloaded MODIS product MOD13A1 for India. Then mosaiced the granules in GeoTIFF format using rasterio Python library. I have the appropriate shapefile to mask the mosaiced raster.
After mosaicing ...
1
vote
1answer
167 views
Creating subset of TIFF image and generating its NDVI in Python?
I have two Raster images, one from band4 with a B4 at the end and another with from band5 with B5 at the end. I want to subset the B5 raster to 800x600, then display it and save it as a GeoTiff. Then ...
0
votes
1answer
140 views
Rasterio showing single band images wrong
I'm using rasterio in Python at the moment, and I have no problem when I'm dealing with RGB satellite images. However, if I try using singleband Sentinel 2 images, or a singleband Sentinel 1 image, ...
1
vote
0answers
50 views
Summarising all the (NDVI) values pixelwise for a masked polygon
A .tif image has been created using the NDVI formula in QGIS. The resulting image shows a region, with varying grades of NDVI values over the image. A .shp file has then been created for an area of ...
1
vote
2answers
97 views
Unknown projection of GRIB file
I have issues extracting geometrical shaped data from a GRIB raster of soil moisture. The file is for example https://gimms.gsfc.nasa.gov/SMOS/jbolten/FAS/L03/20100113_20100115.as1.grib
Using either ...
2
votes
1answer
310 views
Finding pixel location in raster using coordinates?
I'm not a GIS person. I haven't found or at least identified a solution to my problem, so decided to ask.
I have a big rasters and a list of thousands of coordinates. I need to find the pixel ...
1
vote
1answer
57 views
Windowed writing in rasterio?
I have two large datasets, RP1 and RP2. I would like to subtract one from the other, however the datasets are extremely large and so exit out with a MemoryError when I try a direct read:
RP1 = ...
0
votes
1answer
29 views
Windowed read of a raster producing wrong projection?
I'm using rasterio to read in a .tif file and converting it to vectors using the rasterio.features module. For a particularly large .tif file, I thought I'd use the windowed read option, but for some ...
1
vote
0answers
729 views
Error installing Rasterio in PyCharm and in CMD Line
Error installing Rasterio in PyCharm (Settings > Project Interpreter)
Collecting rasterio
Using cached https://files.pythonhosted.org/packages/1a/d5/...
1
vote
0answers
27 views
Geologica data UDM(unusable data mask) have 8 bits per pixel
UNUSABLE DATA MASK FILE(UDM)
The unusable data mask file provides information on areas of unusable data within an image (e.g. cloud and
non-imaged areas).
The pixel size after orthorectification ...
0
votes
1answer
49 views
Generating a raster layer with intersecting line lengths as cell values
I have a polyline vector layer and a template raster dataset. For this polyline vector layer I would like to generate a new raster dataset based on the template raster dataset, but in which each cell ...
0
votes
1answer
96 views
Making subset of Sentinel-2 with rasterio?
I am looking for a solution to my problem with rasterio. I want to download just a specific part of the image from AWS based on lat lon coordinates. Rasterio offers a great tool (window) although it ...
3
votes
0answers
74 views
Why can't I transpose a raster across the pole with rasterio or gdal, and how could it be done?
I have a raster in e.g epsg:4326, it covers all world, from 180E to 180W, from 90S to 90N. I want to re-project a part of it to epsg:3031 and also change the grid size and resolution.
Here is a ...
0
votes
1answer
167 views
Rasterizing by attribute in rasterio
According to Rasterizing shapefiles with GDAL and Python? it seems that gdal allows rasterization of a layer according to a certain attribute
ds = gdal.RasterizeLayer(target_ds, [1], source_layer, ...
0
votes
0answers
207 views
Rasterio : rasterio.drivers() - 'module' object is not callable [duplicate]
I'm trying to pass a script on jupyter-hub with a specific kernel. With my last local installation it worked but now with my kernel I have this error:
TypeError ...
1
vote
0answers
116 views
Polygonizing raster cells in rasterio regardless of value?
I would like to polygonize every single cell of a raster dataset using rasterio. I have been experimenting with the following function based on this question:
def polygonize_raster(dataset):
# ...
1
vote
2answers
50 views
Unable to use GPS coordinates for GEOTiff map
I have been using a map to calculate some distances at sea. The map had a bit low resolution so I got another one with better resolution (both GEOTiff), but with the amounts of distances I need to ...
0
votes
1answer
91 views
How to convert masks generated by rasterio to channel last mask
The problem is about using rasterio to generate mask from GeoJSON file which contains geographic coordinates.
The function that I used is from:
from rasterio.tools.mask import mask
with rasterio....
5
votes
1answer
308 views
Open a raster file with read and write permissions in rasterio
I am trying to write a function which will help me construct a rasterio memory file and return the dataset reader object so I can then play around with it as needed. I am trying to avoid with ...
2
votes
1answer
314 views
rasterio mask method
I am using a simple ASCII raster file to test the rasterio.mask module
ncols 4
nrows 4
xllcorner 0
yllcorner 0
cellsize 0.5
nodata_value -9999
1 2 3 4
5 6 7 8
9 10 11 12
13 14 15 16
I generated the ...
1
vote
0answers
304 views
Writing numpy array to rasterio raster file generates empty file
I'm using Rasterio to read GeoTIFF from Sentinel 2 and calculating NDVI, but when I try to save as a new TIF file, the file turns to transparent and empty file. What am I missing?
That's the code:
...
1
vote
0answers
36 views
rasterio & cogeo
I am trying to use rasterio to read cloud-optimized geotiff stored in a gc bucket.
With method 1, GDAL makes use of the overviews and read is quick. Unfortunately, if the extent is not fully included ...
1
vote
0answers
46 views
rasterio - nodata per band [closed]
I am looking at a way to specify a different nodata value for each band in rasterio.
I tried to:
use the nodata key in the profile with a tuple or list of nodata values,
set the nodatavals attribute,...
0
votes
0answers
103 views
Clipping raster with different UTM zones using rasterio?
I have 2 Landsat images, with an overlapping area:
The one on the left is in WGS84 / UTM32N (EPSG:32632), the one on the right is in WGS84 / UTM33N (EPSG: 32633).
I have a polygon in WGS84 / UTM33N (...
2
votes
1answer
435 views
How can I superimpose a geopandas dataframe on a raster plot?
I have a geopandas dataframe countries with country polygons, and a raster dataset raster (read with rasterio). I've already aligned the two as follows:
countries = countries.to_crs(raster.crs.data)
...
2
votes
1answer
99 views
Reading a WMS/ raster map with rasterio
Can I read a raster map which has been constructed in mapserver as WMS with rasterio?
rasterio.open
('http://...&service=WMS&request=GetMap&...FORMAT=image/tiff')
I tried :
from ...
0
votes
1answer
99 views
Extracting pixels in multiple rasters based on the values of another raster's pixels for time series analysis
I'm working with 52 rasters of a specific region over a period of one year in python. I also have a mask raster of the same area that is a binary color (250 black pixels for disparate areas of ...
