Plot over satellite imagery?

[Mccrone, Paul J CIV FNMOC, N3]

Scott:

Thanks for all your effort on PyART.

Is there a way to use PyART to overlay reflectivity /etc/ over a satellite image?

I have HDF5 and netCDF4 files, each of which I can display using the h5py, netCDF4 , and Basemap python modules.

So If I have a h5py of nc4 object, can I overplot the PyART object over that image? Any ideas as to how to proceed?

Thank again.

Very Respectfully,
===========================
Paul J. McCrone, CCM
Meteorologist, FNMOC/N38DI
Paul.J....@navy.mil
Comm: 831-657-1503
===========================
U N C L A S S I F I E D For Official Use Only (FOUO)

[Scott Collis]

Yep!

Would strongly suggest using Cartopy for one.

Since you can pass the existing axes through you can plot the Sat data
first then the radar data

https://github.com/ARM-DOE/pyart/blob/master/pyart/graph/radarmapdisplay_cartopy.py#L113

BTW:

Strongly suggest you use Xarray for your data. Amazingly easy to do
cartopy pcolormeshes..

Anyone done what Paul is looking for?

I am interested enough to give this ago but it would take me some time
to get to it, because, work..

[Mccrone, Paul J CIV FNMOC, N3]

Ok, I will look at that. Thanks Scott!

Very interested to hear if anyone else has done this.

Very Respectfully,
===========================
Paul J. McCrone, CCM
Meteorologist, FNMOC/N38DI
Paul.J....@navy.mil
Comm: 831-657-1503
===========================
U N C L A S S I F I E D

--
You received this message because you are subscribed to the Google Groups "Py-ART Users" group.
To unsubscribe from this group and stop receiving emails from it, send an email to pyart-users...@googlegroups.com.
For more options, visit https://groups.google.com/d/optout.

[Martin Rugna]

Although I haven't done any overlay with satellite data, once I had to recover some methods of the basemap instance created by PyART (some simple modifications because the graphics defaults to CONUS).

Below is a part of the script I use, though I do not mess with projections because I didn't need to.

By the way, I'm switching to cartopy for satellite plotting so I'm really interested on what follows this conversation.

Hope this helps,

-Martin

display = pyart.graph.RadarMapDisplay(radar)

fig = plt.figure(figsize=(15, 10))

if var == 'TH' or var == 'total_power':

    display.plot_ppi_map(var, sweep=0, vmin=-10, vmax=70,

                         cmap=pyart.graph.cm.NWSRef, resolution='h',

                         embelish=False)

# elif something else not relevant for this

display.basemap.drawcoastlines(color='0.8', linewidth=1.25)

display.basemap.drawcountries(color='w', linewidth=1.25)

display.basemap.drawstates(color='0.8')

To unsubscribe from this group and stop receiving emails from it, send an email to pyart-users+unsubscribe@googlegroups.com.

For more options, visit https://groups.google.com/d/optout.

--
You received this message because you are subscribed to the Google Groups "Py-ART Users" group.

To unsubscribe from this group and stop receiving emails from it, send an email to pyart-users+unsubscribe@googlegroups.com.

[Scott Collis]

I just had a request from a collaborator on a DoE project for the same
thing.. I will be taking a look at this Friday.

[Scott Collis]

Okay!

This uses gridded data from Py-ART but the idea is the same..

https://github.com/scollis/notebooks/blob/master/sgp_goes_nexrad.ipynb

Enjoy and let me know if there are any questions.. Greatest challenge
was finding historic GOES data!

[Mccrone, Paul J CIV FNMOC, N3]

Nice!

Very Respectfully,
===========================
Paul J. McCrone, CCM
Meteorologist, FNMOC/N38DI
Paul.J....@navy.mil
Comm: 831-657-1503
===========================
U N C L A S S I F I E D For Official Use Only (FOUO)