search

OpenDrift / opendrift

Open source framework for ocean trajectory modelling
https://opendrift.github.io
GNU General Public License v2.0
245 stars 120 forks source link

Parsing netCDF CF projection with pyproj.CRS.from_cf instead of using… #1087

Closed knutfrode closed 1 year ago

knutfrode commented 1 year ago

… selfmade method in reader_netCDF_CF_generic

knutfrode commented 1 year ago

@gauteh Dette ser ut til å fikse biffen, og fungerer fint på eksisterande datasett frå thredds.met.no, samt på test-datasett utan proj4-streng: https://thredds.met.no/thredds/dodsC/metusers/eivinds/noproj4_testfiles/meps_lagged_6_h_subset_2 _5km_20230201T00Z.nc

Men kanskje litt bakstreversk å generere proj4 streng til slutt?

knutfrode commented 1 year ago

Ja, kan se ut til at det ikkje er optimalt å konvertere CRS-objektet til proj4-streng:

tests/plotting/test_norkyst.py: 2 warnings
  /opt/conda/envs/opendrift/lib/python3.8/site-packages/pyproj/crs/crs.py:1286: UserWarning: 
You will likely lose important projection information when converting to a PROJ string from another format.
See: https://proj.org/faq.html#what-is-the-best-format-for-describing-coordinate-reference-systems
    proj = self._crs.to_proj4(version=version)

https://proj.org/faq.html#what-is-the-best-format-for-describing-coordinate-reference-systems Men det er slik OpenDrift/readers er bygd pr i dag, so kan kanskje leve med dette til 2.0