Skip to content
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

Can we use the data of Single Polarization Radar? #16

Open
nay-chi-18 opened this issue Feb 13, 2024 · 12 comments
Open

Can we use the data of Single Polarization Radar? #16

nay-chi-18 opened this issue Feb 13, 2024 · 12 comments

Comments

@nay-chi-18
Copy link

Hi @wolfidan,
Let me ask a question. I already installed Py-art and had run it on my machine by using your script and data. My testing was successful. But when I replace our radar data (kyaukphyu_N191710_E0933126_H0092_202311200617_EL+0008P1_S03_ppi_Z_450km_intensity.nc), pyart don't work. The system shows this error message " KeyError: 'range'. Our radar is single polarization type. My question is that can the data of single polarization radar use with Py-art?

@wolfidan
Copy link
Collaborator

Hi @nay-chi-18
No the fact that it is single-pol should not be an issue. What pyart reader did you use for your data? Do you know if your data is of cfradial type?

Thanks

@nay-chi-18
Copy link
Author

nay-chi-18 commented Feb 15, 2024 via email

@wolfidan
Copy link
Collaborator

wolfidan commented Feb 15, 2024

Thanks you @nay-chi-18 ,
That's perfect, I will check your data and let you know,
If you don't want to publicly share your data, feel free to send it to daniel.wolfensberger@meteoswiss.ch. I will test it and then delete it.
Daniel

@wolfidan
Copy link
Collaborator

wolfidan commented Feb 19, 2024

Hi @nay-chi-18

Thank you for the data, I could read it. I have a few constatations:

  • The data you provided is already gridded (Cartesian lat/lon coordinates, instead of polar angle/range coordinates)
  • The data lacks any metadata and does not follow any standard (CFRadial, ODIM or else)

Because the data is not polar the pyart.io.read reader will not work, as it is meant to be used on polar data (hence the missing "Range" in your data).

The reader for gridded CFRadial is pyart.io.read_grid, unfortunately it fails in your case because you have missing metadata, for example the keys origin_altitude, origin_latitude and origin_longitude are missing. They should indicate the coordinates of the radar. You can see the CFRadial standard here.

I must say that pyart is mostly designed to process polar data and most of its processing routines will not work on gridded data (attenuation correction, clutter filtering, KDP processing, ...). Do you also have access to the polar data of your radars? If yes please provide some to me and I will test them.

Thanks

@nay-chi-18
Copy link
Author

nay-chi-18 commented Feb 20, 2024 via email

@wolfidan
Copy link
Collaborator

Thanks @nay-chi-18 ,

I will check your files. Unfortunately, I had to ask for permission to access them. Could you grant access to me?

@nay-chi-18
Copy link
Author

nay-chi-18 commented Feb 21, 2024 via email

@wolfidan
Copy link
Collaborator

Thank you @nay-chi-18

The data seems to be in an unknown binary format. Can you provide information on how they are structured? Do you have a custom reader for them, in python, C or else?
Thank you

@nay-chi-18
Copy link
Author

nay-chi-18 commented Feb 22, 2024 via email

@wolfidan
Copy link
Collaborator

@nay-chi-18
Sorry I cannot see your attached files in the message

@nay-chi-18
Copy link
Author

nay-chi-18 commented Feb 27, 2024 via email

@wolfidan
Copy link
Collaborator

wolfidan commented Mar 5, 2024

@nay-chi-18

No sorry, there was no attachment, please send the files directly to daniel.wolfensberger@meteoswiss.ch.

Thank you

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants