Hello, Is there any recipe on how to (1.) georeference single sweep (!) data,(2.) transform it to cartesian coordinates and then (3.) to a numpy array? Basically what happened here ( Recipe #2: Reading and visualizing an ODIM_H5 polar volume — wradlib ) under [8], but for files that include only on…

Hi @atmoboran , this is the simplest approach to create cartesian representation of the polar data. Note that this is just Nearest neigbour interpolation. More sophisticated interpolation might be used for proper handling. Upper image shows polar source data, lower image cartesian gridded data. im…

[image] kmuehlbauer: Azimuthal Equidistant Projection Hey Kai, thank you very much! This worked well. But now I ran into a subsequent issue: I would like to have the cartesian data on some coordinate system/ projection (e.g. EPSG 32632), so it is comparable with other data (like the CAPPI f…

Hi Boran @atmoboran , I’ve updated the above example. HTH, Kai

Hi Kai, thanks so much! I’ve also done this with the default WG84 projection in lon lat coordinates: # georeference swp = swp.wrl.georef.georeference(crs = wrl.georef.get_default_projection()) # default = WGS 84 in lon lat # stack source data st = swp.stack(onedim=("azimuth", "range")) src = np.…

@atmoboran I’ve added stacked source data and made the code coyp&paste-able. All in all this looks OK, but you do not have an equal-area projection anymore. Is this just for display or is this used for further processing? In that case I’d suggest to stick with an equal area projection. HTH, Kai

@kmuehlbauer Hello again, just want to make sure: after georeferencing the coordinate system is on ground, i.e. the given distances from radar refer to the distance on the ground and not along the beam?! Best regards, Boran

Hi @atmoboran , Yes, after georeferencing the created x,y,z coordinates are in Azimuthal Equidistant (default) or any other given projection with the radar at it’s center. The slant range (along ray) is still available in the range coordinate. HTH, Kai

Georeferencing single sweep data and converting to numpy array

Python

Topic		Replies	Views
Generate Geotiff from Rainbow file with wradlib Retrievals wradlib	10	529	March 7, 2025
Facing difficultyin converting Coordinate system from x,y to lat,lon Retrievals wradlib	7	331	July 5, 2023
Problems with georeferencing Corrections/Data Quality	2	75	April 12, 2025
Spherical to cartesian, interpolation and back to spherical. How? Python wradlib	2	232	March 8, 2024
Conversion from Lat/Lon to Coordinates Python wradlib	6	1170	February 17, 2023

Georeferencing single sweep data and converting to numpy array

Related topics