Hi!
I have been tracking your amazing library, which seems pretty cool. Yet, I have some problems with the Corine Land Cover geotiff (link to download)
As they say, the EPSG is 3035, my final goal is to get some coords in WGS84 (EPSG:4326).
First issue.
If I open the dataset like this:
It seems to understand everything, (from 3035 to 4326), yet when I pull coords (e.g. for regridding purposes), I get kind of a wrong transformation (see image)
Latitude in Norway shouldn't be 59º but something like 70º approx. and longitude values are also kind of shifted (like 2-3º), so if I try to reproject to another WGS84 grid you can see both shifts:
I have tried to use get_coords() insted of get_wgs84_coords() and transform the results using a pyproj transformer, but the error kind of persist... Could you validate it please?
Hi! I have been tracking your amazing library, which seems pretty cool. Yet, I have some problems with the Corine Land Cover geotiff (link to download) As they say, the EPSG is 3035, my final goal is to get some coords in WGS84 (EPSG:4326). First issue. If I open the dataset like this:
I get an error (pyproj error, asking for a CRS) when trying to pull coords, e.g
corine.get_coords(500, 500)
Second issue, if I open manually referencing the EPSG like this
It seems to understand everything, (from 3035 to 4326), yet when I pull coords (e.g. for regridding purposes), I get kind of a wrong transformation (see image)
Latitude in Norway shouldn't be 59º but something like 70º approx. and longitude values are also kind of shifted (like 2-3º), so if I try to reproject to another WGS84 grid you can see both shifts:
I have tried to use
get_coords()
insted ofget_wgs84_coords()
and transform the results using a pyproj transformer, but the error kind of persist... Could you validate it please?