Hi Guys, Wradlib Version: 1.16.2 Python Version: 3.10.5 I am able to read and plot the graph with raa01-rw_10000-2208301350-dwd—bin.gz file with xy co-ordinates. The below code works fine. But I want your help to plot the graph with longitude and latitude. Could you please send me som…

Hi Muru, welcome to our openradar discourse and thanks for the question. I’ll move the discussion to Openradar until we have decided on the final structure. I’ve created a jupyter notebook gist with a code example for your use case here → radolan-rw.ipynb . One more word of clarification. In the no…

Hi Kai, Thank you for your support. It works very well. It prints like below, [Figure_lat_long_grid] <xarray.DataArray 'RW' ()> array(0., dtype=float32) Coordinates: y float64 -3.759e+06 x float64 37.83 lon float64 9.993 lat float64 54.9 Attributes: val…

Hi Muru, the RADOLAN RW data is taken verbatim from the datafile as uint16. The xarray reader is applying CF convention based scale_factor and add_offset (see at ds.RW.encoding) and let xarray do the CF decoding. Please follow-up with the netcdf-c attribute conventions for valid_min and valid_max. …

Hi Kai, Thank you, understood. So rainfall has occurred at this location. lon float64 9.993 lat float64 54.9 Attributes: valid_min: 0 valid_max: 4095 standard_name: rainfall_rate long_name: RW unit: mm h-1 But if you l…

Hi Muru, maybe I was too sparse with my explanations. The valid_min/valid_max are two dataset-wide values and not point-specific. The correct value in the above example is 0: <xarray.DataArray 'RW' ()> array(0., dtype=float32)

Hi Kai, Ohh Ok. def print_nearest_ds(lat, lon): # find nearest xy-grid point to a specific latlon-coordinate abslat = np.abs(dsg.lat - lat) abslon = np.abs(dsg.lon - lon) c = np.maximum(abslon, abslat) ([xloc], [yloc]) = np.where(c == np.min(c)) # Select index location at t…

Hi Muru, sorry that was my bad, I’ve made a mistake when adapting the solution from SO. I’ve fixed the notebook gist (see link above). For corrected code to determine the wanted pixel see below. There are two changes. First, y/lat is the first dim of our 2d array. Second, we need to use .isel to s…

Hi Kai, Thank you so much. It works like gem. I could now see the proper rainfall value [14.7]. <xarray.DataArray 'RW' ()> array(14.7, dtype=float32) Coordinates: y float64 -4.65e+06 x float64 -2.83e+05 lon float64 6.512 lat float64 47.19 Attributes: val…

Reading Radolan RW files

Openradar

Topic		Replies	Views
How can I able to access RADOLAN binary data for specific region for example Lübeck General	4	258	June 11, 2023
How to subset the RADOLAN data file based on the radarlocation? Data	12	363	December 23, 2022
Generate Geotiff from Rainbow file with wradlib Retrievals wradlib	10	551	March 7, 2025
RADOLAN scale factor for RW product? Corrections/Data Quality	0	147	June 7, 2023
Identifying Precipitation types using RADOLAN data (RW Product) Openradar wradlib	2	374	June 4, 2023

Reading Radolan RW files

Related topics