All Questions
10
questions
2
votes
0
answers
278
views
Speeding up WarpedVRT in rasterio (change CRS + resample)
Given a source GeoTIFF of some unknown resolution (better than the target), what is the quickest way to re-project it in-memory to the desired GSD? Masking in the input data should be preserved.
...
- 151
1
vote
1
answer
431
views
Inconsistent results between gdalWarp (with cutline) and rasterio.mask
Cropping with GDAL cutline vs Rasterio.mask with crop=True results in arrays with different sizes (by one row). I've tested the different Rasterio options but no luck. I'm hoping to match the GDAL ...
- 429
0
votes
0
answers
412
views
How to find patches with data in Raster and then how to extract those patches?
I have a single band raster file as drawn in the labeled image attached.
The size of raster is around 4k pixel in each x and y axis. Large part of this raster is empty. Only in small portions of the ...
- 109
3
votes
1
answer
751
views
GeoTIFF is stretched sideway when translating CRS to EPSG:4326 from EPSG:32654
I want to convert a GeoTIFF's coordinate reference system to EPSG:4326 from EPSG:32654. But the output image is stretched to the side.
I have tried two different methods (gdalwarp and rasterio), but ...
- 41
-2
votes
1
answer
604
views
Using Python gdalwarp/rasterio to read gif file and world file, re-project to another picture [closed]
I want to re-project a radar image (https://radar.weather.gov/Conus/RadarImg/latest_radaronly.gif) which is in NAD83/EPSG4326 to WGS84/Pseudo-Mercator/EPSG3857. There is a world file associated with ...
- 1
0
votes
1
answer
1k
views
Create VRT from multiple hdf files with Rasterio/GDAL
I have several Landsat 8 Scenes, each consisting of 11 bands. I want to mosaic them by creating a VRT and afterwards exporting the mosaic as TIFF again. I tried using GDAL via Python. But maybe it ...
- 685
1
vote
0
answers
141
views
Rasterio 0-360 lon inputfile to -180:180lon WarpedVRT output
I have used rasterio to open netcdf files and by defining a WarpedVRT extract a subarea that can be read with the vrt.read() command.
It works like a charm except when the longitude grid is 0:360 ...
- 11
1
vote
0
answers
119
views
Getting subset with same cordinates for different rasters using gdal.translate?
I am a new into Python. I have processed 1000 NDVI rasters covering larger area. I need to burn those NDVI raster with Water mask. The water mask came for whole world I need to subset AOI and also ...
- 57
1
vote
1
answer
2k
views
Transforming a image to a shape defined by image corners in earth coordinates - Python
I have a source satellite image. It is a rectangular shape x by y. I also have the stored lon/lat coordinates of the corners and the center point in the xml meta file. Is there a way with the gdal or ...
- 43
3
votes
0
answers
537
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 ...
- 398