wradlib Documentation - Bitbucket

More documents

Recommendations

Info

$LaTeX-Tutorium III - Bitbucket$

wradlib Documentation, Release 0.1.1 latitude coordinates based on WGS84 lonc : array of floats longitude coordinates based on WGS84 projstr : string proj.4 projection string. Can be conveniently created by using function wradlib.georef.create_projstr Returns output : a tuple of 2 arrays (x and y coordinates) Examples Gauss-Krueger Zone 2: “+proj=tmerc +lat_0=0 +lon_0=6 +k=1 +x_0=2500000 +y_0=0 +ellps=bessel +towgs84=598.1,73.7,418.2,0.202,0.045,-2.455,6.7 +units=m +no_defs” Gauss-Krueger Zone 3: “+proj=tmerc +lat_0=0 +lon_0=9 +k=1 +x_0=3500000 +y_0=0 +ellps=bessel +towgs84=598.1,73.7,418.2,0.202,0.045,-2.455,6.7 +units=m +no_defs” UTM Zone 51 on the Northern Hemishpere “+proj=utm +zone=51 +ellps=WGS84” UTM Zone 51 on the Southern Hemishpere “+proj=utm +zone=51 +ellps=WGS84 +south” >>> import wradlib.georef as georef >>> # This is Gauss-Krueger Zone 3 (aka DHDN 3 aka Germany Zone 3) >>> gk3 = create_projstr("gk", zone=3) >>> latc = [54.5, 55.5] >>> lonc = [9.5, 9.8] >>> gk3_x, gk3_y = georef.project(latc, lonc, gk3) 3.5.7 wradlib.georef.create_projstr wradlib.georef.create_projstr(projname, **kwargs) Conveniently supports the construction of proj.4 projection strings Currently, the following projection names (argument projname) are supported: “aeqd”: Azimuthal Equidistant needs the following keyword arguments: lat_0 (latitude at projection center), lon_0 (longitude at projection center), x_0 (false Easting, also known as x-offset), y_0 (false Northing, also known as y-offset) “gk” : Gauss-Krueger (for Germany) only needs keyword argument zone (number of the Gauss-Krueger strip) “utm” : Universal Transmercator needs keyword arguments zone (integer) and optionally hemisphere (accepted values: “south”, “north”) see Wikipedia entry for UTM zones. “dwd-radolan” : RADOLAN Composite Coordinate System no additional arguments needed. Polar stereographic projection used by the German Weather Service (DWD) for all Radar composite products. See the final report on the RADOLAN project (available at http://www.dwd.de/RADOLAN) for details. 46 Chapter 3. Library Reference
wradlib Documentation, Release 0.1.1 Parameters projname : string (proj.4 projection acronym) kwargs : depends on projname - see above! Returns output : string (a proj.4 projection string) Examples >>> # Gauss-Krueger 2nd strip >>> print create_projstr("gk", zone=2) >>> # UTM zone 51 (northern hemisphere) >>> print create_projstr("utm", zone=51) >>> # UTM zone 51 (southern hemisphere) >>> print create_projstr("utm", zone=51, hemisphere="south") 3.5.8 wradlib.georef.projected_bincoords_from_radarspecs wradlib.georef.projected_bincoords_from_radarspecs(r, az, sitecoords, projstr, range_res=None) Convenience function to compute projected bin coordinates directly from radar site coordinates and range/azimuth specs Parameters r : array array of ranges [m]; r defines the exterior boundaries of the range bins! (not the centroids). Thus, values must be positive! az : array sitecoords : tuple array of azimuth angles containing values between 0° and 360°. The angles are assumed to describe the pointing direction fo the main beam lobe! The first angle can start at any values, but make sure the array is sorted continuously positively clockwise and the angles are equidistant. An angle if 0 degree is pointing north. projstr : string proj.4 projection string range_res : float range resolution of radar measurement [m] in case it cannot be derived from r (single entry in r-array) 3.6 Interpolation Interpolation allows to transfer data from one set of locations to another. This includes for example: • interpolating the data from a polar grid to a cartesian grid or irregular points • interpolating point observations to a grid or a set of irregular points • filling missing values, e.g. filling clutters Nearest Nearest-neighbour interpolation in N dimensions. Continued on next page 3.6. Interpolation 47
Page 1: wradlib Documentation Release 0.1.1
Page 4 and 5: 5.3 Users . . . . . . . . . . . . .
Page 6 and 7: wradlib Documentation, Release 0.1.
Page 100 and 101:
wradlib Documentation, Release 0.1.
Page 102 and 103:
Page 104 and 105:
Page 106 and 107:
show all

wradlib Documentation - Bitbucket

Create successful ePaper yourself

Delete template?

Save as template?