initial commit/should be working

- use -w file[.kml] to write to file

2 files changed, 174 insertions, 0 deletions

diff --git a/gk2kml.py b/gk2kml.py
new file mode 100755
index 0000000..37d3521
--- /dev/null
+++ b/gk2kml.py
@@ -0,0 +1,114 @@
+#!/usr/bin/env python
+#
+# gk2kml
+#
+# converts Gauss-Krueger coordinates to WGS84
+# and also output a .kml so we can easily
+# load it up in Google Earth without copy-pasting
+#
+# 2022-04-27 Oliver Meckmann
+
+import argparse
+from textwrap import dedent
+from sys import exit
+
+from pyproj import CRS # Coordinate Reference System
+from pyproj import Transformer
+
+# https://spatialreference.org/ref/epsg/wgs-84/
+# https://en.wikipedia.org/wiki/World_Geodetic_System
+# https://epsg.io/?q=Gauss-Kruger%20kind%3APROJCRS
+
+def run_argparser():
+    argparser = argparse.ArgumentParser(\
+            formatter_class=argparse.RawTextHelpFormatter,
+            description=dedent('''
+                Convert Gauss-Kru(e)ger Coordinates to WGS84
+                '''))
+    argparser.add_argument(
+        '-z',
+        '--zone',
+        default='3',
+        help=dedent('''
+            Input Gauss-Kruger Zone (Default: 3)
+            2 = 3deg Zone 2: Germany
+            3 = 3deg Zone 3: Germany onshore between 7°30'E and 10°30'E
+            4 = 3deg Zone 4: Germany onshore between 10°30'E and 13°30'E
+            5 = 3deg Zone 5: Germany onshore east of 13°30'E
+            '''),
+        nargs='?',
+    )
+    argparser.add_argument(
+        '-w',
+        '--write',
+        help='Write to file [filename] (.kml gets attached)',
+        type=str,
+        nargs='?'
+    )
+    argparser.add_argument(
+        'coord',
+        metavar='coord',
+        type=str,
+        help=dedent('''
+            [easting northing] ([Rechtswert Hochwert])
+            '''),
+        nargs=2)
+
+    args = argparser.parse_args()
+    return args
+
+def main():
+    args = run_argparser()
+
+    if len(args.coord[0]) < 6 or len(args.coord[1]) < 6:
+        print('Error: coordinates too short/inaccurate!')
+        exit()
+
+    easting = args.coord[0]
+    northing = args.coord[1]
+    # sometimes we don't have THAT accurate coords, we then fill up with a 0
+    if len(easting) == 6:
+        easting += '0'
+    if len(northing) == 6:
+        northing += '0'
+
+    # set zone
+    if args.zone == '2':
+        src_crs = 'EPSG:5682'
+    elif args.zone == '3':
+        src_crs = 'EPSG:5683'
+    elif args.zone == '4':
+        src_crs = 'EPSG:5684'
+    elif args.zone == '5':
+        src_crs = 'EPSG:5685'
+    else:
+        print('Error: Wrong zone!')
+        exit()
+
+    transformer = Transformer.from_crs(src_crs, 'EPSG:4326')
+    result = transformer.transform(easting, northing)
+    
+    if(str(result[0]) != 'inf'):
+        print(str(result[0]) + ', ' + str(result[1]))
+
+        if args.write:
+            print('Writing to kml file: ' + args.write + '.kml')
+            template = open('template.kml', 'r', encoding='utf-8')
+            outf = open(args.write + '.kml', 'w', encoding='utf-8', errors='ignore')
+            line = template.readline()
+
+            while line:
+                line = line.replace('$NAME', args.write)
+                line = line.replace('$LON', str(result[1]))
+                line = line.replace('$LAT', str(result[0]))
+                #print(line, end='')
+                outf.write(line)
+                line = template.readline()
+
+            outf.close()
+
+    else:
+        print('Error: no valid result! Check input coords!')
+
+if __name__ == '__main__':
+    main()
diff --git a/template.kml b/template.kml
new file mode 100644
index 0000000..9469872
--- /dev/null
+++ b/template.kml
@@ -0,0 +1,60 @@
+<?xml version="1.0" encoding="UTF-8"?>
+<kml xmlns="http://www.opengis.net/kml/2.2" xmlns:gx="http://www.google.com/kml/ext/2.2" xmlns:kml="http://www.opengis.net/kml/2.2" xmlns:atom="http://www.w3.org/2005/Atom">
+<Document>
+	<name>$NAME.kml</name>
+	<Style id="s_ylw-pushpin_hl">
+		<IconStyle>
+			<scale>1.4</scale>
+			<Icon>
+				<href>http://maps.google.com/mapfiles/kml/shapes/open-diamond.png</href>
+			</Icon>
+		</IconStyle>
+		<LabelStyle>
+			<scale>0.3</scale>
+		</LabelStyle>
+		<ListStyle>
+		</ListStyle>
+	</Style>
+	<StyleMap id="m_ylw-pushpin">
+		<Pair>
+			<key>normal</key>
+			<styleUrl>#s_ylw-pushpin</styleUrl>
+		</Pair>
+		<Pair>
+			<key>highlight</key>
+			<styleUrl>#s_ylw-pushpin_hl</styleUrl>
+		</Pair>
+	</StyleMap>
+	<Style id="s_ylw-pushpin">
+		<IconStyle>
+			<scale>1.2</scale>
+			<Icon>
+				<href>http://maps.google.com/mapfiles/kml/shapes/open-diamond.png</href>
+			</Icon>
+		</IconStyle>
+		<LabelStyle>
+			<scale>0.3</scale>
+		</LabelStyle>
+		<ListStyle>
+		</ListStyle>
+	</Style>
+	<Placemark>
+		<name>$NAME</name>
+        <description>$NAME</description>
+		<LookAt>
+			<longitude>$LON</longitude>
+			<latitude>$LAT</latitude>
+			<altitude>0</altitude>
+			<heading>0</heading>
+			<tilt>0</tilt>
+			<range>999.9998585735252</range>
+			<gx:altitudeMode>relativeToSeaFloor</gx:altitudeMode>
+		</LookAt>
+		<styleUrl>#m_ylw-pushpin</styleUrl>
+		<Point>
+			<gx:drawOrder>1</gx:drawOrder>
+            <coordinates>$LON,$LAT</coordinates>
+		</Point>
+	</Placemark>
+</Document>
+</kml>