24.12.2013 Views

wradlib Documentation - Bitbucket

wradlib Documentation - Bitbucket

wradlib Documentation - Bitbucket

SHOW MORE
SHOW LESS

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

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!