Hello,
I’m trying to work with Cordex NetCDF files (Europe).
The main, challenging, issue I’m facing is the rotation of coordinates used by Cordex.
Googling and trying for some days I found a Proj4 string that seems to work. At least it works in a python caller with pyproj lib.
crsROT = CRS("+proj=ob_tran +o_proj=longlat +o_lon_p=-162 +o_lat_p=39.25 +lon_0=180 +to_meter=0.01745329")
crs4326 = CRS('EPSG:4326')
x=feature.getCoordinate(0)
transformer = Transformer.from_crs(crs4326, crsROT)
p1 = transformer.transform(xf1], x[0])
The process is understandably very slow. Some improvement can be applied working with list of coordinates, but I hope FME can provide some better standard solutions.
I put the same Proj4 string in a PROJReprojector. Also considering that working with raster and then to clip the desired region is faster. The parameters I’m using are in the attached image.
Unfortunately, the results are completely different. Both applying the transformer to vector or Raster. For sure I’m using the transformer in the wrong way.
Attached also the correct and wrong representation (in color the correct one).
Suggestions?
TIA.
Roberto