-
Notifications
You must be signed in to change notification settings - Fork 94
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
Reconstruct CRS from WKT only #460
Labels
Comments
5 tasks
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
Feature Request
Is your feature request related to a problem? Please describe.
For some areas,
area.crs.to_wkt()
returns a dictionary that contains only thecrs_wkt
key and nothing else. I don't know if this is CF-conform, but this is the responsibility of pyproj. For some discussion, see the comments on pytroll/satpy#2226.For example:
Currently, pyresample cannot construct an AreaDefinition using only this CRS and the corresponding coordinates. Therefore, Satpy's
satpy_cf_nc
reader cannot read the files written by itscf
writer. See pytroll/satpy#2226 for the Satpy perspective.The NetCDF header may look like:
where trying to load the area will fail with
Describe the solution you'd like
I would like that, if the grid mapping variable contains only the attribute
crs_wkt
and no other information, pyresample should do its best to construct an area using the CRS WKT with other available information.Describe any changes to existing user workflow
There shouldn't be any.
Additional context
My current workaround involves writing all lat/lons and using a
SwathDefinition
. This is inefficient, in particular since I'm then resampling back to the area I had when I was writing the NetCDF file to begin with.The text was updated successfully, but these errors were encountered: