JavaScript disabled You need to activate JavaScript in order to use OpenProject! Learn more
added beam_height_ft_doviak, which uses a slightly different formular than collier. The advantage of this formula is that it is invertible for theta (which may be of interest for visibility studies). Also, it is deriveable using simple trigonometry, whereas the derivation of the Collier-formula is rather obscure.
Added literature reference to beam_height_ft
fixed wrong treatment of large and small overflow sectors.
added calls to interp_atten and bisectReferenceAttenuation. Everything not tested. bisectReferenceAttenuation is also missing a reduction of b (for those rare cases where this is necessary), but this might also be addressed by another function before the bisection.
removed some of the mind-bogglingness of multidimensionality by reshaping the original data to something three-dimensional ('the rest', azimuth, range).
added interp_atten. Function to interpolate reference attenuation to invalid beams. Not yet tested.
made sector_filter work with arbitrary sector sizes (=dynamic origin).
Added testing routines for atten module
minor change in tutorial
resized some tutorial images and changed some tutorial text
removed the oversized tutorial images from the readDX tutorial
updated the io- and the vis-tutorial georef_example.py: multiply the radar-range array by 1000, due to SI-convention
changes in wrl.ipol.interpolate_polar: 1. the mask array is optional, if the data array is masked 2. in terms of recommended efficiency ONLY masked values will be interpolated 3. in consequence for recalculation of linear interpolation (values at the rim can't...
added convenient interpolation function for polar data: ipol.interpolate_polar
Added first (example of a) testing file.
Added preliminary code for the new (and almost complete) implementation of the Kraemer algorithm.
Renamed some variables in atten.py for better understanding: k -> pia (because it integrates all foregoing k's) kn -> k (the explicit indication that it stands for the specific attenuation of bin n is now (first rename) unessential
Attention for users of georef: changed units for earth radius in polar2latlon so this function expects ranges in SI-units([m]) This affects the functions polar2polyvert, polar2centroids and projected_bincoords_from_radarspecs Further inout parameter re (earth radius) is removed in functions, which do not...
Attention for users of georef.projected_bincoords_from_radarspecs: changed units for range-array from [km] to [m] (SI-units) will change that in the near future for the other functions, too
typo in the reference of ipol.linear documentation
Added an Error Message, in case of assigning a range gate array with less than two values to georef._get_range_resolution()
additional parameter (range_resolution) for functions polar2centroids and projected_bincoords_from_radarspecs, just for the case the sitecoords shall be calculated from one range, the resolution cannot be derived from the range array (_get_range_resolution)
changed io.readDX in a way that the elevations and azimuths are returned anyway as an attrs-dictionary. By the way the DX data is now immidiately converted to dBZ data in the readDX function. Just can't imagine that anyone needs the RVP6-units. And if so it can be easily converted back.
Added contstrained attenuation correction routine. Extended iterations to reduce parameter b as well, if first iteration of a won't fulfill all constraints.
remove old cluttercor image
changed cluttercorrectio tutorial added the figure closure in vis
fixed bug that occurs when reading a great many of netCDF-files (the have to be closed explizitly)
added PolarPlot class and polar_plot2 function. class PolarPlot is a an attempt to refactor the plotting capabilities of the polar_plot function into something that provides more flexibility when creating custom plots while avoiding large function signatures and an...
added arguments classes and extend to vis.polar_plot
<classes> is a sequence of class boundaries. Colors are assigned according to the class boundaries and the colormap argument. It is recommended to use this plotting mode because the behaviour of automatic stretching is quite unpredictable....
further edits on tutorial_reading_dx
editing the tutorial_reading_dx
several modifications in the documentation concerning
index: added some introductory text and an image (we can replace that) in order to create a nice welcome
getting started: much more details in order to make things easier for Python newbies
tutorials: minor edits...
Removed autogenerated files from version control
made vis.polar_plot behave nicely when fig keyword argument is set but not axpos argument.
added clutter tutorial
merged Revisions 146 and 148
added a draft for a new map plot function for grids in vis.py ...but is not operational, yet.
limited the range array in read_EDGE_netcdf according to range_lim
dummy push
io.read_EDGE_netcdf: changed all range units to meters (SI-convention) added the range_lim default and the unit-specification to the docstring
io.read_EDGE_netcdf: set an range_lim default value of 200km avoid errors in case range_lim is greater than the maximum range of the data
io.read_EDGE_netcdf: included the toctree added an range limitation asinput argument rounding of theta0 instead of cuting the decimals rolling the azimuths instead of sorting added max_range to the returned attributes
renamed SUB2dBZ to read_EDGE_netcdf and modified and enhanced the function. It still needs to be clarified whether the data array needs to flipped.
modified project function in georef.py so it can also do inverse projection
minor edit in georef_example.py
added convenience function projected_bincoords_from_radarspecs to documentation minor edit in projected_bincoords_from_radarspecs (now returns ravelled coordinate arrays)
added the missing import
added io.SUB2dBZ for reading philipinean s-band radar data
added tutorials and some minor edits
added colormap argument to vis.polar_plot
correctAttenuationHJ in atten.py processes even if the input gateset contains NaNs; in this case all corresponding k(attenuation)-beams are set to NaN, too; but the remaining PIA will be calculated;
if this approves, the code for the other to correction functions should also be...
atten.py: moved the NaN-check before first calculations and memory allocations
change in make.bat: deleting source/generated before each built
Fixed bug that occurred if only one point is passed to the IDW Interpolator
removed superfluous calculations estimated speedup is roughly by a factor of 1.7
removed superfluous call to np.where
added _z2rEnhanced_mdcorr function same result as z2rEnhanced but roughly 50 times faster.
Removed the trafo.rvp2r function: for clarity such processchains should be executed in a more basic way e.g. zr.z2r(trafo.idecibel(trafo.rvp2dBZ(x))), in order to avoid undue code and import errors (caused by endless cross-referencing).
Fixed a refencing bug when nnearest equals 1 in idw
Added convenience function ipol.interpolate
This function should not only offer a uniform interface for all interpolation classes. Particularly, it allows to efficiently deal with missing values in the interpolation source values (i.e. np.nan values) if your input value array contains missing values...
changed warning message
Submitted version including Maik's additions.
Added short absract for ERAD2012
Made all correctAttenuation functions behave the same for attenuation overflows
Adjusted overflow handling so that it fits the documentation.
a few readability and logical corrections. - k is defined as np.zeros, so there is no need to zero out any subset of k directly afterwards. - If there are *any* NaNs in the gateset we cannot continue. So let's say that in the code and throw a proper exception instead of returning something useless....
changed class.rst which should enable *some* kind of class documentation. Seems to work although it generates a host of strange error messages during build.
Fixed some documentation bugs
added **kwargs to polar_plot This way additional keywords can be passed to pcolormesh like vmin and vmax.
removed jsmath-part. jsmath is a large download and needs its fonts installed on the *user's* system, which most people won't do just to read our documentation. So I'll abandon this road for now and stick with LaTeX.
changed dtype of output array to float32 in ipol.py Idw
minor edits in ipol docstrings
Temporarily added the current dependencies to pyproj and basemap until we will have found an agreement on the dependy policy.
import of utils in __init__.py
added util to __init__.py
minor edit in util.py docstrings
completing function aggregate_in_time several minro edits in docstrings of georef.py
reworked togrid so that multidimensional data will be handled correctly. Now more than one image may be gridded simultaneously.
added new moduile util.py
added convenience function projected_bincoords_from_radarspecs in georef.py
the check if the input array is nan, now don't ends in an exception raise, but in a Zero-attenuation array in the output in the next version i will check if it is better to return a nan array for attenuation
minro edits in index.rst, ipol.rst, verify.rst, verify.py
minor edit
still experimental map plotting in vis.py and vis_example.py, please do not use, yet
added a class to verify.py which extracts radar bins in the neighbourhood of a rain gage, added a corresponding example workflow in verify_example.py
bugfix in georef
adding a basemap plot of polar radar data to vis and vis_example - still needs refinement. This is just the outline. This plot type allows polar data to be plotted on arbitrary map projections including shapefile overlay.
some refinement of the example for georef
lots of changes in georef and the corresponding example: introduced a new approach for computing the polygon vertices and a more explicit method for computing the centroids of the bins. But both are founded on polar2latlon
added a new function polar2polyvert in georef.py which does not use centroids to compute the polygon vertices
added an example file georef_example.py which should contain a workflow for georeferencing radar observations only from the location of the radar and its azimuth and range specs
added projection examples for UTM
added necessary imports
added the project function to georef.py, I hope the bad merging problem is now resolved...
merging went wrong, trying to recommit thomas revision no 83 and then late the new changes from me
merge conflict - removed conflicting function for while being
added funciton to convert lat/lon coordinates to native map coordinates (projection)
changed the rotation properties of the polar plot
substituted the dictionary for attenuation coefficients by direct parameterization.
some docstring changes
Also available in: Atom