Repository navigation
CRS setting and copying #252
Description
I try to set the crs, and than to copy the raster. Turns out the self set crs, isn't copied.:
The functions for copying I use are write_crs and set_crs. Both give different results.
raster = open_rasterio(r"PATH\FILENAME.xyz") # No crs set
raster.rio.write_crs(32633) # crs is None
raster.rio.set_crs(32633) # result: CRS.from_epsg(32633)
crs = CRS("EPSG:32633")
raster.rio.write_crs(crs.to_wkt()) # None
raster.rio.set_crs(crs.to_wkt()) # result: CRS.from_epsg(32633)
# when copiing the raster
raster.rio.crs # CRS.from_epsg(32633)
r1 = raster.copy()
r1.rio.crs # NoneProblem description
The copy method only copies the original data, not the changed?
Expected Output
Environment Information
rioxarray (0.2.0) deps:
rasterio: 1.2.0
xarray: 0.16.2
GDAL: 3.1.4
Other python deps:
scipy: 1.6.0
pyproj: 3.0.0.post1
System:
python: 3.8.6 | packaged by conda-forge | (default, Jan 25 2021, 22:54:47) [MSC v.1916 64 bit (AMD64)]
executable: C:\Users\remote\anaconda3\envs\PROJEKTNAME\python.exe
machine: Windows-10-10.0.19041-SP0
Installation method
- conda-forge
Activity
snowman2 commentedon Feb 22, 2021
If you want to be able to copy and have the CRS persist, you need to use
rio.write_crs. Nothing can be done inrioxarrayif you userio.set_crs. Hopefully that gets you where you need to be.snowman2 commentedon Feb 22, 2021
I noticed you have
raster.rio.write_crs(crs.to_wkt()), so my previous comment was only partially helpful.You need to either:
inplace=Truesnowman2 commentedon Jun 2, 2022
For those finding this issue, I recommend reading the Getting Started: Introductory Information references for rioxarray. The Coordinare Reference System Management section is relevant for this issue.