GEOTRANSFORMER: Difference between revisions

Revision as of 11:17, 1 March 2016

Basic Information

Name of Program

GEOTRANSFORMER

Version-Date

February 2016

Description-Date

February 2016

Catchwords

coordinate transformation
coordinate system
Gauss-Krueger
Universal Transverse Mercator
Universal Polar Stereographic
European Terrestrial Reference System 1989 (ETRS89)
world geodetic system 1984 (WGS84)
european datum 1950 (ED50)
Potsdam datum (Bessel 1841)
Krassovski ellipsoid
NTv2 method
GNTRANS-WSV
BETA2007
Rijksdriehoeksmeting (RD, coordinate system of the netherlands)
universal transverse mercator (UTM)
European petrol service group (EPSG)

Short Description of Functionality

This program transforms coordinates between different coordinate systems. Currently implemented coordinate systems are:

Gauß-Krüger (stripe 2 to 5)
spherical coordinates (longitude, latitude)
Universal Transverse Mercator (UTM)
Rijksdatum (Netherlands)
local metric central projection

Currently implemented geodetic 7-parameter transformations:

European Terrestrial Reference System 1989 (ETRS89)
World Geodetic System 1984 (WGS84)
European Datum 1950 (ED50)
Federal Intstitute for cartography und geodesy, standard parameter
Federal Intstitute for cartography und geodesy, latitudes >52,3 degrees N
BAW C. Maushake,
Krassovsky standard (Wikipedia)
Krassovsky (WSA Stralsund)
Federal Intstitute for cartography und geodesy, latitudes between 50,3 degrees N and 52,3 degrees N
Federal Intstitute for cartography und geodesy, latitudes <50,3 Grad N

Additionally the NTv2 grid shift method is provided for datum transformation of DHDN90/STN to ETRS89 and back, based on GNTRANS-WSV or BETA2007.

You can run the software in batch mode. There is no limit in amount of data to work with, because files are processed sequentially. To process 25 Mio. pairs of coordinates (Gauß-Krüger/DHDN <-> UTM/ETRS89) the software uses 280.68 minutes on a common INTEL CPU including Input/Output.

Input-Files

triangulated mesh (file type TICAD,gitter05.dat/bin)
digitized polylines (file type digi.gkk) or
digitized polylines (file type insel.dat) or
polygon file (file type poly.dat) or
polygone file to protect edges of a grid (file type nodes.save) oder
plot frames (file type frames.dat, center and size, or
description of boundary cells (file type rgz.dat) or
ASCII triple data, file with soundings (file type geom.dat) or
time series at single station (file type boewrt.dat) or
position description (file type geopos.dat) or
triangulated UNTRIM-VC mesh (file type untrim_grid.dat, version Vincenzo Casulli) or
grid of the numerical method DELFT3D (file type delft3d.grd) or
triangulated UNTRIM-BAW mesh (file type untrim_grid.dat, BAW version) or
system file for specific locations (file type location_grid.dat) or
file IPDS to initialize models area wide (file type ipds.dat) or
world files which describe the georeference of an Image (typical file extensions are .pngw, .jpgw, .gifw) or
mesh of profiles (file type profil05.bin) .

Output-Files

triangulated mesh (file type TICAD,gitter05.dat/bin)
digitized polylines (file type digi.gkk) or
digitized polylines (file type insel.dat) or
polygon file (file type poly.dat) or
polygone file to protect edges of a grid (file type nodes.save) oder
plot frames (file type frames.dat, center and size, or
description of boundary cells (file type rgz.dat) or
ASCII triple data, file with soundings (file type geom.dat) or
time series at single station (file type boewrt.dat) or
position description (file type geopos.dat) or
triangulated UNTRIM-VC mesh (file type untrim_grid.dat, version Vincenzo Casulli) or
grid of the numerical method DELFT3D (file type delft3d.grd) or
triangulated UNTRIM-BAW mesh (file type untrim_grid.dat, BAW version) or
system file for specific locations (file type location_grid.dat) or
file IPDS to initialize models area wide (file type ipds.dat) or
world files which describe the georeference of an Image (typical file extensions are .pngw, .jpgw, .gifw) or
mesh of profiles (file type profil05.bin) .

Methodology

Input and output list are opened at the same time. Every coordinate point will be read, transformed and written on the output list immediately. So the size of the input file is unlimited. Comments and other informations are copied "as is"!
The formatting of decimal coordinates depends on the actual projection of coordinates (geographic or map coordinates). File type 08 (ASCII triple data) allows user defined Output Format via Environment variable GEOMFMT.
The program first calculates spherical coordinates (WGS84) from input projected coordiantes. Then the new geodetic datum and projection is applied. In the datum cases DHDN or STN in input and/or output data, the NTv2 method with grid shift file *.gsb is provided as Default alternative. The program should be used with the actual NTv2 grid shift files provided by the GNTRANS-WSV System conform to this system. It can be used with the BETA2007.gsb file to get conformity with the ATKIS system.
on input in files of type 08 (geom.dat) the separators semicolon, tabulator, comma and colon are allowed. Additionally sperical coordinates are allowed to be written in nautical style, e.g. 4°17'16.62931"E 51°29'40.02423"N. Instead of characters (° ' ") also (d m s) or (^ ' ") may be used. The nautical format must include the information of hemisphere (N/S) and the orientation with respect to Greenwich (E/W).
on input in files of type 09 (boewrt.dat) or 10 (geopos.dat), 13 (untrim_grid.dat, BAW version) and 14 (location_grid.dat) the coordinate reference system is first changed to the input system, if an internal CRS Information is available. In files of type 01 (gitter05.dat), 02 (digi.gkk) and 08 (geom.dat) the( well known!) coordinate reference system of the input files can be specified through the comment "C CRS=#####" at the beginning of the file. "#####" here is an EPSG code which is supported by libgeodesy. In this case also a Transformation to the input system is done first. The input system must have a valid EPSG code too.
in input files of type 16 (world files) also the avarage pixel width and height transformation parameters are changed by assuming a virtual image of 800x800 pixels.
To all ASCII formats which are not NAMELIST based, a file header with the current coordinate information is added as comments.
UTM easting coordinates with 7 and 8 digits are recognized and changed into 6 digits representation. For recognition of valid 8 und 7 digits easting values both the number of digits and the actual UTM zone are used.

EPSG codes wich are supported by libgeodesy are:

Gauss-Krueger 3 degree / DHDN: EPSG 31466-31469 (EPSG = 31464 + stripe number)
Gauss-Krueger 3 degree / STN : EPSG 02398-02399
UTM / ETRS89 : EPSG 25828-25838 (EPSG = 25800 + zone number)
UTM / ED50 : EPSG 23028-23038 (EPSG = 23000 + zone number)
UTM / WGS84 : EPSG 32601-32660 (EPSG = 32600 + zone number)
UTM / WGS84 : EPSG 32701-32760 (EPSG = 32700 + zone number)
UPS / WGS84 : EPSG 32661 and 32761
RD Amersfoort New : EPSG 28992
Geographic ETRS89 : EPSG 4258
Geographic WGS84 : EPSG 4326
Geographic ED 50 : EPSG 4230
Geographic DHDN : EPSG 4314
Geographic Amersfoort : EPSG 4289

It is reccommended, to store the coordinate reference system as an EPSG code inside the basic datasets (raw data), if it is known. This can be done by the "magic comment" inside the most ASCII file types. Within file types, where storage is foreseen for an EPSG code, it should be placed there.

Actual datasets should be stored with UTM/ETRS89 coordinates of the valid zone or in geographic coordinates (ETRS89 or WGS84) . Transforming older datasets with Gauß-Krüger coordinates into UTM/ETRS89 should be done preferably with the NTv2 method after GNTRANS-WSV. The transformation grid is valid only for German area including coastal waters. Coordinates outside this area must not be transformed by NTv2 but should be transformed using a valid Helmert transformation parameter set.

Transforming coordinates is not a trivial task. Therefor the result of the transformation should be tested for correctness and usability.

Program(s) to run before this Program

all programs wich use above listed file formats

Program(s) to run after this Program

all programs wich use above listed file formats

Additional Information

Language

Fortran95

Additional software

NTv2 grids dhdn_to_etrs89_wsv_v1.gsb, stn_to_etrs89_wsv_v1.gsb and BETA2007.gsb

Original Version

G. Seiß (main program, file I/O, libgeodesy)

Maintenance

G. Seiß

Documentation/Literature

interactive dialog, self explaining
$PROGHOME/examples/geotransformer
Junkins, D.R. and Farley, S.A. (1995) National Transformation version 2 Users Guide, Geodetic Survey Division Geomatics Canada
Lott, Roger und OGP Geodesy Working Group (2015) Coordinate Conversions and Transformations including Formulas. Geomatics Guidance Note (IOGP Publication 373), 7. Part 2.
Bruijne, Arnoud de; van Buren, Joop; Kösters, Anton; van der Marel, Hans (2005): De geodetische referentiestelsels van Nederland. Geodetic reference frames in the Netherlands. Hg. v. Nederlandse Commissie voor Geodesie Netherlands Geodetic Commission. Delft. Download here.
Federal Fairway Authority (Eds.) (2012): Handlungsanweisung für die Transformation der Datenbestände der WSV in das System ETRS89/UTM. In German language, with contribution of Hendrik Hampe, Sudau, Gunther Braun, Cornelius Zschunke, Egon Feigel, Helga Panknin et al.

back to Program Descriptions

Overview

@@ Line 2: / Line 2: @@
 |name_de=GEOTRANSFORMER
 |name=GEOTRANSFORMER
-|version=January 2016
+|version=February 2016
-|version_descr=January 2016
+|version_descr=February 2016
 |catchwords=coordinate transformation<br/>
 coordinate system<br/>
@@ Line 88: / Line 88: @@
 * The formatting of decimal coordinates depends on the actual projection of coordinates (geographic or map coordinates). File type 08 (ASCII triple data) allows user defined Output Format via Environment variable GEOMFMT.
 * The program first calculates spherical coordinates (WGS84) from input projected coordiantes. Then the new geodetic datum and projection is applied. In the '''datum cases DHDN or STN in input and/or output data''', the NTv2 method with grid shift file *.gsb is provided as Default alternative. The program should be used with the actual NTv2 grid shift files provided by the GNTRANS-WSV System conform to this system. It can be used with the BETA2007.gsb file to get conformity with the ATKIS system.
-* on input in files of type 08 (geom.dat) the separators semicolon, comma and colon are allowed.
+* on input in files of type 08 (geom.dat) the separators semicolon, tabulator, comma and colon are allowed. Additionally sperical coordinates are allowed to be written in nautical style, e.g. 4°17'16.62931"E 51°29'40.02423"N. Instead of characters (° ' ") also (d m s) or (^ ' ") may be used. The nautical format must include the information of hemisphere (N/S) and the orientation with respect to Greenwich (E/W).
 * on input in files of type 09 ([[BOEWRT.DAT|boewrt.dat]]) or 10 ([[GEOPOS.DAT|geopos.dat]]), 13 ([[UNTRIM_GRID.DAT|untrim_grid.dat]], BAW version) and 14 ([[LOCATION_GRID.DAT|location_grid.dat]]) the coordinate reference system is first changed to the input system, if an internal CRS Information is available. In files of type  01 (gitter05.dat), 02 (digi.gkk) and 08 (geom.dat) the( well known!) coordinate reference system of the input files can be specified through the comment "C CRS=#####" at the beginning of the file. "#####" here is an EPSG code which is supported by libgeodesy. In this case also a Transformation to the input system is done first. The input system must have a valid EPSG code too.
 * in input files of type 16 (world files) also the avarage pixel width and height transformation parameters are changed by assuming a virtual image of 800x800 pixels.
@@ Line 100: / Line 100: @@
 * UTM / ED50               : EPSG 23028-23038 (EPSG = 23000 + zone number)
 * UTM / WGS84              : EPSG 32601-32660 (EPSG = 32600 + zone number)
+* UTM / WGS84              : EPSG 32701-32760 (EPSG = 32700 + zone number)
 * UPS / WGS84              : EPSG 32661 and 32761
 * RD Amersfoort New        : EPSG 28992
@@ Line 121: / Line 122: @@
 |contact_maintenance=[mailto:guntram.seiss@baw.de G. Seiß]
 |documentation=interactive dialog, self explaining<br/> $PROGHOME/examples/geotransformer <br/> Junkins, D.R. and Farley, S.A. (1995) '''N'''ational '''T'''ransformation '''v'''ersion '''2''' Users Guide, Geodetic Survey Division Geomatics Canada  <br/>
-Lott, Roger und OGP Geodesy Working Group (2015) Coordinate Conversions and Transformations including Formulas. [http://www.iogp.org/pubs/373-07-2.pdf Geomatics Guidance Note (IOGP Publication 373), 7. Part 2].
+Lott, Roger und OGP Geodesy Working Group (2015) Coordinate Conversions and Transformations including Formulas. [http://www.iogp.org/pubs/373-07-2.pdf Geomatics Guidance Note (IOGP Publication 373), 7. Part 2].<br/>
+Bruijne, Arnoud de; van Buren, Joop; Kösters, Anton; van der Marel, Hans (2005): De geodetische referentiestelsels van Nederland. Geodetic reference frames in the Netherlands. Hg. v. Nederlandse Commissie voor Geodesie Netherlands Geodetic Commission. Delft. [http://www.ncgeo.nl/index.php?option=com_k2&view=item&id=2361:gs-43-a-de-bruijne-de-geodetische-referentiestelsels-van-nederland&Itemid=178&lang=nl Download here]. <br/>
+Federal Fairway Authority (Eds.) (2012): Handlungsanweisung für die Transformation der Datenbestände der WSV in das System ETRS89/UTM. In German language, with contribution of Hendrik Hampe, Sudau, Gunther Braun, Cornelius Zschunke, Egon Feigel, Helga Panknin et al.
 }}