I have a .UF file with the wrong latitude and longitude

141 views
Skip to first unread message

Eliton lima de figueiredo

unread,
Dec 2, 2015, 2:29:45 PM12/2/15
to Py-ART Users
Dear Py-ART users,

I have a .UF file with the wrong latitude and longitude (-28 775, -53.07) for the radar location. I have with the correct coordinates (-29, -54).
I have an attached figure showing the correction that I want to perform.

I would also like to correct or change the background map because it does not show the desired  political map.

sincerely

link UF file:
http://200.132.99.1/prev/temp/2014-10-19--03_20_23_00.mvol.UF

import numpy as np
import matplotlib.pyplot as plt
import pyart

# read in the file, create a RadarMapDisplay object
filename = '2014-10-19--03_20_23_00.mvol.UF'
radar = pyart.io.read(filename)
display = pyart.graph.RadarMapDisplay(radar)

# plot the second tilt
display.plot_ppi_map('corrected_reflectivity', 1, vmin=0, vmax=65,
                     min_lon=-65, max_lon=-45, min_lat=-40, max_lat=-23,
                     lon_lines=np.arange(-65, -45, 5),
                     lat_lines=np.arange(-40, -23, 5), resolution='h',
                     cmap='pyart_NWSRef')
plt.show()

Forum_2.png

Scott Collis

unread,
Dec 2, 2015, 2:48:39 PM12/2/15
to Eliton lima de figueiredo, Py-ART Users
Hey Elton,
First: nice data!

Second: you can edit the lat/lon yourself..

radar.latitude['data'] = np.array([mylat])
radar.longitude['data'] = np.array([mylon])

so:

import numpy as np
import matplotlib.pyplot as plt
import pyart

# read in the file, create a RadarMapDisplay object
filename = '2014-10-19--03_20_23_00.mvol.UF'
radar = pyart.io.read(filename)
radar.latitude['data'] = np.array([-29])
radar.longitude['data'] = np.array([-54])


display = pyart.graph.RadarMapDisplay(radar)

# plot the second tilt
display.plot_ppi_map('corrected_reflectivity', 1, vmin=0, vmax=65,
                     min_lon=-65, max_lon=-45, min_lat=-40, max_lat=-23,
                     lon_lines=np.arange(-65, -45, 5),
                     lat_lines=np.arange(-40, -23, 5), resolution='h',
                     cmap='pyart_NWSRef')
plt.show()


as for boundaries.. Check out:
http://matplotlib.org/basemap/users/geography.html

December 2, 2015 at 1:29 PM
--
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.

Eliton lima de figueiredo

unread,
Dec 2, 2015, 4:17:18 PM12/2/15
to Py-ART Users
worked, thank you very much for your help.

Eliton lima de figueiredo

unread,
May 17, 2016, 5:50:35 PM5/17/16
to Py-ART Users


Dear All,
I am trying to open a .mvol radar file using the function "read_gamic", but it is not working. I get several error messages. Is it possible to open a .mov file using "read_gamic"? If not, what alternative do I have?
Thanks in advance.

http://200.132.99.1/prev/temp/2014-10-19--03_20_23_00.mvol

filename = '2014-10-19--03_20_23_00.mvol'
radar = pyart.aux_io.read_gamic(filename)

Kai Muehlbauer

unread,
May 18, 2016, 2:37:53 AM5/18/16
to pyart...@googlegroups.com
Hi,

I opened an issue here https://github.com/ARM-DOE/pyart/issues/534.

An alternative would be to just open the mvol with an hdf5-reader (it's
actually an hdf5) and extract the needed data for adding to a
radar-object. I'm short of ideas if this is possible out of the box with
pyart, but surely Jonathan has some.

If this is really urgent, you might try wradlib-functions to read the
data and then import it into an pyart radar object or just visualize it.

Cheers,
Kai


Am 17.05.2016 um 23:50 schrieb Eliton lima de figueiredo:
--
Kai Muehlbauer
Meteorological Institute University of Bonn
Auf dem Huegel 20 | +49 228 739083
D-53121 Bonn | kai.mue...@uni-bonn.de
santiago.png

Kai Muehlbauer

unread,
May 19, 2016, 4:29:22 AM5/19/16
to pyart...@googlegroups.com
Hi,

I have a PR here which solves the reading problem of the OP as a first
step (https://github.com/ARM-DOE/pyart/pull/535).

The attached synthetic RHI is produced from the read radar-object.

Cheers,
Kai
santiago1.png
Reply all
Reply to author
Forward
0 new messages