pyart icon indicating copy to clipboard operation
pyart copied to clipboard

Constant Altitude Plan Position Indicator (CAPPI)

Open meteoswiss-mdr opened this issue 7 years ago • 14 comments

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.

meteoswiss-mdr avatar Sep 19 '16 06:09 meteoswiss-mdr

@meteoswiss-mdr I agree completely with the recommendation. Thanks for submitting this. Adding labels.

jjhelmus avatar Sep 19 '16 15:09 jjhelmus

So Just checking I understand @meteoswiss-mdr : x and y axes in km instead of lat/lon?

scollis avatar Jul 13 '17 16:07 scollis

Yes, that is what I meant

Sent from my iPhone

On 13 Jul 2017, at 18:23, Scott <[email protected]mailto:[email protected]> 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.

meteoswiss-mdr avatar Jul 14 '17 07:07 meteoswiss-mdr

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.

vlouf avatar Jul 31 '19 06:07 vlouf

My own code for this can be found here:

https://meteoswiss-mdr.github.io/pyrad/_modules/pyrad/graph/plots_vol.html#plot_cappi

meteoswiss-mdr avatar Jul 31 '19 07:07 meteoswiss-mdr

I guess this feature hasn't been added to pyart yet? I need some help on this.

mustafaalicelik avatar Feb 12 '21 10:02 mustafaalicelik

@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.

zssherman avatar Feb 12 '21 15:02 zssherman

Any improvement?

mustafaalicelik avatar Feb 14 '21 12:02 mustafaalicelik

Can’t you just make a single level grid?

-sent from a mobile device-

On Feb 14, 2021, at 6:53 AM, malicelik [email protected] wrote:

 Any improvement?

— You are receiving this because you commented. Reply to this email directly, view it on GitHub, or unsubscribe.

scollis avatar Feb 14 '21 13:02 scollis

@mustafaalicelik I'm hoping to have a try at this sometime tomorrow or early next week.

zssherman avatar Feb 18 '21 22:02 zssherman

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.

vlouf avatar Mar 16 '21 23:03 vlouf

Has this CAPPI product already been added to pyart?

gabrielcassol avatar Nov 10 '22 13:11 gabrielcassol

@gabrielcassol Hello! At the moment, it has not, still on our list for an addition before the next big release.

zssherman avatar Nov 10 '22 15:11 zssherman

@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 avatar Nov 10 '22 17:11 syedhamidali