-
Notifications
You must be signed in to change notification settings - Fork 17
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
Option to use fixed angle step in Rainbow reader #75
Comments
Thanks @egouden for raising this. Since xradar 0.0.9 we have angle reindexing implemented. Please have a look here for details here https://docs.openradarscience.org/projects/xradar/en/latest/notebooks/angle_reindexing.html If you know the layout of your scans, you can read the data like this:
reindex_angle=dict(tolerance=1.0, start_angle=0, stop_angle=360, angle_res=1.0, direction=1)
ds = xr.open_dataset(rainbow_file, group="sweep_0", engine="rainbow", reindex_angle=reindex_angle) So for ray at 0.5deg center-beam the algorithm would take the nearest ray which has a maximum distance of 1.0deg. Update: python example code, minor fixes |
Thank you @kmuehlbauer for the suggestion. I am wondering if an automatic option based on the file metadata would be interesting. |
That would be specific to that backend, but we might have the chance to use |
I think reading the nominal angles should be the default for all readers for several reasons:
The actual angles (due to imperfection of the radar system) could be read as optional/mandatory metadata (ODIM2.4). What do you think? |
@egouden I fully agree to present the dataset in nominal angles. That's exactly whats I'd make this either nominal or accurate angles, not both. If a user really want to have both |
Yeah - I agree that setting some reasonable defaults/adding a bit more logic on attempting to calculate some of the fields required would be great! This is a similar practice to reading in terabytes of climate model output where the exact latitude/longitude can vary a bit. When using |
In Rainbow, the radar measurements can be sampled for a given fixed angle step (fixselect>Anglestep)
It would be nice to have the option to use this as regular nominal coordinates instead of the actual azimuth angles.
The text was updated successfully, but these errors were encountered: