MeteoIODoc 20240411.d3bdb3cb
Available coordinate systems

Geographic coordinates will be transparently and automatically converted to lat/lon and any other coordinate system that the client program uses. However, in order to do so, the input coordinate system must be specified. In order to output geolocalized data, the desired coordinate system must also be specified for the outputs (in the output section). This is done through the use of the COORDIN and COORDPARAM keys (see the documentation for each plugin).

There are two ways of supporting a given coordinate system: through the use of an adhoc implementation (that becomes part of MeteoIO) or through the use of an external library, Proj [ref: http://trac.osgeo.org/proj/]. The current internal implementations are the following (given by their keyword):

  • CH1903 or CH1903+ for coordinates in the Swiss Grid [ref: http://geomatics.ladetto.ch/ch1903_wgs84_de.pdf] (epsg codes, respectively 21781 and 2056)
  • UTM for UTM coordinates, the zone must be specified in the parameters, for example 31T [ref: http://www.oc.nps.edu/oc2902w/maps/utmups.pdf] (epsg codes as 32600+zoneNumber in the northern hemisphere or as 32700+zoneNumber in the southern hemisphere)
  • UPS for Universal Polar Stereographic coordinates (the zone, either N or S, must be specified in the parameters). [ref: J. Hager, J. Behensky, B. Drew, THE UNIVERSAL GRIDS: Universal Transverse Mercator (UTM) and Universal Polar Stereographic (UPS), 1989, Defense Mapping Agency, DMATM 8358.2] (epsg codes as 32661 for the north pole and 32761 for the south pole)
  • LOCAL for local coordinate system (using the horizontal and vertical distance from a reference point, see Coords::geo_distances for the available choice of distance algorithms)

Such an example of use is the following:

COORDSYS = UTM
COORDPARAM = 31T

On the other hand, when using the Proj library for handling the coordinate conversion, the EPSG codes of as illustrated below (21781 is the EPSG code for the CH1903 coordinate system. Such a code is 32767 at the maximum): the chosen projection must be specified (such codes can be found at http://spatialreference.org/ref/epsg/?page=1)

COORDSYS = PROJ
COORDPARAM = 21781