|
Get PROJ.4 representation of CRS object
pyproj.CRS has a to_proj4 method that should be able to do that: https://pyproj4.github.io/pyproj/stable/examples.html
pyproj.CRS has a to_proj4 method that should be able to do that: https://pyproj4.github.io/pyproj/stable/examples.html
|
By
Alan Snow
· #680
·
|
|
Rasterio 1.2b1
I think this is the problem: https://pyproj4.github.io/pyproj/stable/gotchas.html#internal-proj-error-sqlite-error-on-select
I think this is the problem: https://pyproj4.github.io/pyproj/stable/gotchas.html#internal-proj-error-sqlite-error-on-select
|
By
Alan Snow
· #674
·
|
|
Rasterio 1.2b1
The new wheel now includes the proj.db, which is close to the size of grids that were removed. Also, PROJ requires more dependencies, so that is also a likely factor. It might work, but I wouldn't rec
The new wheel now includes the proj.db, which is close to the size of grids that were removed. Also, PROJ requires more dependencies, so that is also a likely factor. It might work, but I wouldn't rec
|
By
Alan Snow
· #669
·
|
|
Rasterio 1.2b1
Thanks for working on this Sean! I verified that PROJ_NETWORK=ON works with a basic test: rio_geom.py: from rasterio.warp import transform_geom geometry = [ { "type": "Polygon", "coordinates": [ [ [-9
Thanks for working on this Sean! I verified that PROJ_NETWORK=ON works with a basic test: rio_geom.py: from rasterio.warp import transform_geom geometry = [ { "type": "Polygon", "coordinates": [ [ [-9
|
By
Alan Snow
· #665
·
|
|
opening file with forced CRS
This should be fixed in 1.1.5 IIRC: #1248
This should be fixed in 1.1.5 IIRC: #1248
|
By
Alan Snow
· #628
·
|
|
Best way to sample all data in a raster
I believe that this is what you are looking for: https://rasterio.readthedocs.io/en/latest/topics/resampling.html
I believe that this is what you are looking for: https://rasterio.readthedocs.io/en/latest/topics/resampling.html
|
By
Alan Snow
· #623
·
|
|
Request for comment: rasterio wheels on PyPI and GDAL 3.1
My first thought is that if you do make the switch, I think it would make sense to do it in the 1.2 release. > I have heard from those who would rather stick with 2.4.4 for a while yet due to increase
My first thought is that if you do make the switch, I think it would make sense to do it in the 1.2 release. > I have heard from those who would rather stick with 2.4.4 for a while yet due to increase
|
By
Alan Snow
· #606
·
|
|
can resampling not populate cell if source has nodata?
Here is an example using rioxarray: https://corteva.github.io/rioxarray/stable/examples/interpolate_na.html It uses scipy griddata.
Here is an example using rioxarray: https://corteva.github.io/rioxarray/stable/examples/interpolate_na.html It uses scipy griddata.
|
By
Alan Snow
· #547
·
|
|
Interconverting between pyproj and rasterio CRS objects
Add some history for reference: https://rasterio.groups.io/g/dev/topic/switch_from_proj_4_format_to/28991934
Add some history for reference: https://rasterio.groups.io/g/dev/topic/switch_from_proj_4_format_to/28991934
|
By
Alan Snow
· #445
·
|
|
Interconverting between pyproj and rasterio CRS objects
Here are instructions for rasterio to geopandas: https://gis.stackexchange.com/a/350231/144357 For geopandas to rasterio, you will need something like: if LooseVersion(rasterio.__gdal_version__) > Loo
Here are instructions for rasterio to geopandas: https://gis.stackexchange.com/a/350231/144357 For geopandas to rasterio, you will need something like: if LooseVersion(rasterio.__gdal_version__) > Loo
|
By
Alan Snow
· #441
·
|
|
Creating rasterio dataset without IO
Hi Nic, I think you are looking for rioxarray. An example of what you want to do is here: https://corteva.github.io/rioxarray/stable/examples/clip_geom.html Best, Alan
Hi Nic, I think you are looking for rioxarray. An example of what you want to do is here: https://corteva.github.io/rioxarray/stable/examples/clip_geom.html Best, Alan
|
By
Alan Snow
· #435
·
|
|
Geodatabase, FileGeodatabase, and Esri Geodatabase
Hello, It appears the GeoPackage format can support vector and raster data in GDAL similar to the FIleGDB format. https://gdal.org/drivers/raster/gpkg.html https://gdal.org/drivers/vector/gpkg.html No
Hello, It appears the GeoPackage format can support vector and raster data in GDAL similar to the FIleGDB format. https://gdal.org/drivers/raster/gpkg.html https://gdal.org/drivers/vector/gpkg.html No
|
By
Alan Snow
· #363
·
|
|
Flaw in rasterio design?
Hi Amine, rasterio is definitely a great tool. But, if you are looking for an interface to it similar to pandas, you may be interested in rioxarray. It wraps the rasterio code and gives you an xarray
Hi Amine, rasterio is definitely a great tool. But, if you are looking for an interface to it similar to pandas, you may be interested in rioxarray. It wraps the rasterio code and gives you an xarray
|
By
Alan Snow
· #335
·
|
|
Numpy error when masking a Landsat image with a polygon [Rasterio 1.0.28]
After looking at the docs, I think you need to remove `.read(1)` as that returns a numpy array and `mask` expects a dataset: https://rasterio.readthedocs.io/en/stable/api/rasterio.mask.html#rasterio.m
After looking at the docs, I think you need to remove `.read(1)` as that returns a numpy array and `mask` expects a dataset: https://rasterio.readthedocs.io/en/stable/api/rasterio.mask.html#rasterio.m
|
By
Alan Snow
· #311
·
|
|
Numpy error when masking a Landsat image with a polygon [Rasterio 1.0.28]
I wonder if you are looking for `from rasterio.features import geometry_mask`? https://rasterio.readthedocs.io/en/stable/api/rasterio.features.html#rasterio.features.geometry_mask Example code: https:
I wonder if you are looking for `from rasterio.features import geometry_mask`? https://rasterio.readthedocs.io/en/stable/api/rasterio.features.html#rasterio.features.geometry_mask Example code: https:
|
By
Alan Snow
· #310
·
|
|
multi-dimensional support
I think this would be quite useful. i made a start at it using subdatasets in rioxarray here: https://github.com/corteva/rioxarray/pull/33 But, it is missing 1 dimensional variables such as time at th
I think this would be quite useful. i made a start at it using subdatasets in rioxarray here: https://github.com/corteva/rioxarray/pull/33 But, it is missing 1 dimensional variables such as time at th
|
By
Alan Snow
· #289
·
|
|
Geotiff max/min coordinates question
I am not entirely sure what you are referring to, but I am assuming you want the bounds of the raster in lat/lon? If so, then `EPSG:4326` is a pretty standard projection for lat/lon coordinates. So, y
I am not entirely sure what you are referring to, but I am assuming you want the bounds of the raster in lat/lon? If so, then `EPSG:4326` is a pretty standard projection for lat/lon coordinates. So, y
|
By
Alan Snow
· #285
·
|
|
problem with window sizes when parallelizing a funciton
You are indeed correct Sean. Thanks for catching that and providing the correct answer! It seems I skimmed through the original message too fast.
You are indeed correct Sean. Thanks for catching that and providing the correct answer! It seems I skimmed through the original message too fast.
|
By
Alan Snow
· #222
·
|
|
problem with window sizes when parallelizing a funciton
Hi Javier, The blocksize is not a dynamic attribute. It represents the blocksize of the file on disk. So, if you want to change the blocksize, you will need to write it to a new file with a new blocks
Hi Javier, The blocksize is not a dynamic attribute. It represents the blocksize of the file on disk. So, if you want to change the blocksize, you will need to write it to a new file with a new blocks
|
By
Alan Snow
· #220
·
|
|
Is it possible to create a rasterio object without exporting a file?
Hi Leah, If you would like to try the rioxarray as an xarray extension, it is now available on pypi. An example of what you would like to do is here: https://corteva.github.io/rioxarray/html/examples/
Hi Leah, If you would like to try the rioxarray as an xarray extension, it is now available on pypi. An example of what you would like to do is here: https://corteva.github.io/rioxarray/html/examples/
|
By
Alan Snow
· #198
·
|