wradlib Documentation - Bitbucket
wradlib Documentation - Bitbucket
wradlib Documentation - Bitbucket
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
<strong>wradlib</strong> <strong>Documentation</strong>, Release 0.1.1<br />
latitude coordinates based on WGS84<br />
lonc : array of floats<br />
longitude coordinates based on WGS84<br />
projstr : string<br />
proj.4 projection string. Can be conveniently created by using function<br />
<strong>wradlib</strong>.georef.create_projstr<br />
Returns output : a tuple of 2 arrays (x and y coordinates)<br />
Examples<br />
Gauss-Krueger Zone 2: “+proj=tmerc +lat_0=0 +lon_0=6 +k=1 +x_0=2500000 +y_0=0 +ellps=bessel<br />
+towgs84=598.1,73.7,418.2,0.202,0.045,-2.455,6.7 +units=m +no_defs”<br />
Gauss-Krueger Zone 3: “+proj=tmerc +lat_0=0 +lon_0=9 +k=1 +x_0=3500000 +y_0=0 +ellps=bessel<br />
+towgs84=598.1,73.7,418.2,0.202,0.045,-2.455,6.7 +units=m +no_defs”<br />
UTM Zone 51 on the Northern Hemishpere “+proj=utm +zone=51 +ellps=WGS84”<br />
UTM Zone 51 on the Southern Hemishpere “+proj=utm +zone=51 +ellps=WGS84 +south”<br />
>>> import <strong>wradlib</strong>.georef as georef<br />
>>> # This is Gauss-Krueger Zone 3 (aka DHDN 3 aka Germany Zone 3)<br />
>>> gk3 = create_projstr("gk", zone=3)<br />
>>> latc = [54.5, 55.5]<br />
>>> lonc = [9.5, 9.8]<br />
>>> gk3_x, gk3_y = georef.project(latc, lonc, gk3)<br />
3.5.7 <strong>wradlib</strong>.georef.create_projstr<br />
<strong>wradlib</strong>.georef.create_projstr(projname, **kwargs)<br />
Conveniently supports the construction of proj.4 projection strings<br />
Currently, the following projection names (argument projname) are supported:<br />
“aeqd”: Azimuthal Equidistant<br />
needs the following keyword arguments: lat_0 (latitude at projection center), lon_0 (longitude at projection<br />
center), x_0 (false Easting, also known as x-offset), y_0 (false Northing, also known as y-offset)<br />
“gk” : Gauss-Krueger (for Germany)<br />
only needs keyword argument zone (number of the Gauss-Krueger strip)<br />
“utm” : Universal Transmercator<br />
needs keyword arguments zone (integer) and optionally hemisphere (accepted values: “south”, “north”) see<br />
Wikipedia entry for UTM zones.<br />
“dwd-radolan” : RADOLAN Composite Coordinate System<br />
no additional arguments needed.<br />
Polar stereographic projection used by the German Weather Service (DWD) for all Radar composite products.<br />
See the final report on the RADOLAN project (available at http://www.dwd.de/RADOLAN) for details.<br />
46 Chapter 3. Library Reference