Closed meteoswiss-mdr closed 3 weeks ago
@meteoswiss-mdr I agree completely with the recommendation. Thanks for submitting this. Adding labels.
So Just checking I understand @meteoswiss-mdr : x and y axes in km instead of lat/lon?
Yes, that is what I meant
Sent from my iPhone
On 13 Jul 2017, at 18:23, Scott notifications@github.com<mailto:notifications@github.com> wrote:
So Just checking I understand @meteoswiss-mdrhttps://github.com/meteoswiss-mdr : x and y axes in km instead of lat/lon?
— You are receiving this because you were mentioned. Reply to this email directly, view it on GitHubhttps://github.com/ARM-DOE/pyart/issues/597#issuecomment-315129417, or mute the threadhttps://github.com/notifications/unsubscribe-auth/APYptep7uo8CbGSqCF6DlM-0yoGiL-eyks5sNkR-gaJpZM4KAKPN.
I recently needed to do a CAPPI, so I have coded up a notebook that do just that: https://github.com/vlouf/radar_cappi
It wouldn't be a lot of work to add this to Py-ART, so I can do this if you want.
Strictly speaking, no interpolation is made, just looking for the nearest neighbour at each grid point for a given altitude.
My own code for this can be found here:
https://meteoswiss-mdr.github.io/pyrad/_modules/pyrad/graph/plots_vol.html#plot_cappi
I guess this feature hasn't been added to pyart yet? I need some help on this.
@mustafaalicelik Seems this one went under the radar, no pun intended, I'll see if I can get the pyrad implementation provided by meteoswiss into pyart.
Any improvement?
Can’t you just make a single level grid?
-sent from a mobile device-
On Feb 14, 2021, at 6:53 AM, malicelik notifications@github.com wrote:
Any improvement?
— You are receiving this because you commented. Reply to this email directly, view it on GitHub, or unsubscribe.
@mustafaalicelik I'm hoping to have a try at this sometime tomorrow or early next week.
Hi @hamidrizvi ,
Frankly I don't remember well, it was quite some time ago. 1/ It's to convert the azimuth angle (North is 0 degree and going clockwise) to mathematical angle (0 degree is x-axis and going counter-clockwise). 2/ Altitude is in meter I believe, you should use 2000, not 2. 3/ Nearest neighbour.
Has this CAPPI product already been added to pyart?
@gabrielcassol Hello! At the moment, it has not, still on our list for an addition before the next big release.
@gabrielcassol I would suggest having a look at this, https://github.com/syedhamidali/CAPPI-NETCDF/blob/main/simple_cappi.ipynb
You can make the data gridded, then plot it at whatever height you want. Here I have plotted Max-CAPPI in this example.
@syedhamidali - is there any interest in contributing this to the main Py-ART repo? If not, we can draft something from our side and ping you for review.
@mgrover1 Sounds good! I’ll get started on a script for plotting max_cappi with any projection. If that aligns with what you’re thinking, awesome. If not, just let me know what you’d like me to focus on instead.
It would be convenient to create a griddisplay class in line with the radardisplay class to be able to plot CAPPIs in km units.
In addition it may be convenient to be able to plot gridmapdisplay and radarmapdisplay in km units.