Last active
June 10, 2022 12:25
-
-
Save vmarquar/6becfd9fa86605c21451830d46b7d4dc to your computer and use it in GitHub Desktop.
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
#### This scripts collects a bunch of helpful proj transformation strings and shows their | |
#### usage for coordinate systems that are mainly used in germany and their resprective federal states. | |
#### The script focuses on open source projects gdal / ogr and proj. | |
#### INFO: | |
## Die Transformation erfolgt immer nach dem Prinzip: | |
## (YRe, XHo)GK → (Br1, Ln1) Ellipsoid(GK) <-> (Br2, Ln2) Ellipsoid(ETRS89) → (YEa, XNo)UTM | |
## Damit werden zunächst alle projizierten Hoch- und Rechtswerte aus dem Gauß-Krüger-System in Lat/Long-Werte im Gauß-Krüger-System (geographisch) transformiert. | |
## Anschließend werden die geographischen Koordinaten mit GK-Bezugssystem in das UTM-Bezugssystem überführt (es blieben jedoch Lat/Long Werte). | |
## Im letzten Schritt werden die geopraphischen Koordinaten der UTM-"Kugel" in projizierte UTM-Easting und Northing Koordinaten umgewandelt. | |
#### Gute Zusammenfassung der NTV2-Transformations-Verfahren in Deutschland | |
## https://www.akgsoftware.de/docs/de/infravision/b56/akgcad/webhelp/topics/koordinaten_r_transformation.html#import_r_transformation__section_d55_1nv_vdb | |
# Bayern(BY KANU NTV2) | |
# Affin-Transformation (6 Parameter) | |
# Helmert (4 Parameter) | |
# Nordrhein-Westfalen (TRABBI) | |
# Rheinland-Pfalz (Affin) | |
# Hessen (LET) | |
# Hessen (HETA2010) | |
# Niedersachsen (GNTRANS_NI) | |
# Sachsen-Anhalt (NTv2LSBB_LSA) | |
# Sachsen-Anhalt (LSA_TRANS) | |
# Schleswig-Holstein (TrafoFF) | |
# Sachsen (NTv2_SN) | |
# Brandenburg (NTv2-BB) | |
# Thüringen (NTv2) | |
# Baden-Württemberg (BWTA2017) | |
# Meridiantransformation | |
# DB_Ref (GNTRANS Geo++) | |
# Derzeit öffentlich bereitgestellte NTV2-Gitter: | |
# Bundesland Gitterweite in B, L Gitterweite in m | |
# Deutschland BeTA2007 6' x 10' ca. 17.000 x 18.000 | |
# Thüringen 30'' x 60'' ca. 927 x 1230 | |
# Brandenburg 30'' ca. 927 x 615 | |
# Hamburg 250 | |
# Hessen 8'' x 12'' ca. 250 | |
# Sachsen 4'' x 6'' ca. 125 | |
# Bayern 1'' ca. 30 | |
# Baden-Württemberg 1,5'' x 2,5'' ca. 50 | |
#### TESTED NTV2 TRANSFORMATIONS | |
## DE-BY: | |
## EPSG:31468 to EPSG:25832 (high-precision, +- 1 cm, using BY_KANU.gsb grid) > tested and works | |
## EPSG:25832 to EPSG:31468 (high-precision, +- 1 cm, using BY_KANU.gsb grid) > tested and works | |
## DE-Gesamt | |
## EPSG:31466 to EPSG:25832 (lower-precision, +- 0.9 m, using BETA2007.gsb grid) > tested and works | |
## EPSG:31467 to EPSG:25832 (lower-precision, +- 0.9 m, using BETA2007.gsb grid) > tested and works | |
## EPSG:31468 to EPSG:25832 (lower-precision, +- 0.9 m, using BETA2007.gsb grid) > tested and works | |
#### TRANSFORMATIONS USING OGR2OGR | |
## DE-BY | |
## EPSG:31468 to EPSG:25832 | |
./ogr2ogr -s_srs "+proj=tmerc +lat_0=0 +lon_0=12 +k=1.000000 +x_0=4500000 +y_0=0 +ellps=bessel +datum=potsdam +nadgrids=$GRID_PATH/BY_KANU.gsb +units=m +no_defs" -t_srs EPSG:25832 -f "GPKG" -lco ENCODING=UTF-8 /Users/Valentin/Desktop/linienzug_UTM32.gpkg /Users/Valentin/Desktop/linienzug_GK4.gpkg linienzug_GK4 | |
./ogr2ogr -s_srs "+proj=tmerc +lat_0=0 +lon_0=12 +k=1.000000 +x_0=4500000 +y_0=0 +ellps=bessel +datum=potsdam +nadgrids=$GRID_PATH/BY_KANU.gsb +units=m +no_defs" -t_srs "+proj=utm +zone=32 +ellps=GRS80 +units=m +no_defs +nadgrids=@null" -a_srs EPSG:25832 -f "GPKG" -lco ENCODING=UTF-8 /Users/Valentin/Desktop/linienzug_UTM32.gpkg /Users/Valentin/Desktop/linienzug_GK4.gpkg linienzug_GK4 | |
## DE-BY | |
## EPSG:31468 to EPSG:25832 | |
## Abweichugen mit offiziellen Testdaten von ca. 2,25 mm | |
## beide untenstehenden Lösungen funktionieren exakt gleich | |
./ogr2ogr -s_srs "+proj=tmerc +lat_0=0 +lon_0=12 +k=1.000000 +x_0=4500000 +y_0=0 +ellps=bessel +datum=potsdam +nadgrids=$GRID_PATH/BY_KANU.gsb +units=m +no_defs" -t_srs "+proj=utm +zone=32 +ellps=GRS80 +units=m +no_defs +nadgrids=@null" -a_srs EPSG:25832 -f "GPKG" -lco ENCODING=UTF-8 /Users/Valentin/Documents/04_GIS/04_Projekte/2022-06-09_NTV2_Tests/VEKTOR/Umwandlung_GK4-zu_UTM32-cmd-line.gpkg /Users/Valentin/Documents/04_GIS/04_Projekte/2022-06-09_NTV2_Tests/VEKTOR/DE-BY_Testpunkte_als_GK4.gpkg DE-BY_Testpunkte_als_GK4 | |
./ogr2ogr -s_srs "+proj=tmerc +lat_0=0 +lon_0=12 +k=1.000000 +x_0=4500000 +y_0=0 +ellps=bessel +datum=potsdam +nadgrids=$GRID_PATH/BY_KANU.gsb +units=m +no_defs" -t_srs EPSG:25832 -f "GPKG" -lco ENCODING=UTF-8 /Users/Valentin/Documents/04_GIS/04_Projekte/2022-06-09_NTV2_Tests/VEKTOR/Umwandlung_GK4-zu_UTM32-cmd-line_v2.gpkg /Users/Valentin/Documents/04_GIS/04_Projekte/2022-06-09_NTV2_Tests/VEKTOR/DE-BY_Testpunkte_als_GK4.gpkg DE-BY_Testpunkte_als_GK4 | |
#### TRANSFORMATIONS USING THE DIRECT CS2CS COMMAND LINE UTILITY FROM PROJ | |
# 0) Set specific gdal version and set proj's verbositiy to high | |
cd /Applications/QGIS-LTR.app/Contents/MacOS/bin/ | |
./gdalinfo --version | |
export PROJ_DEBUG=on | |
export GRID_PATH=/Users/Valentin/.qgis3_ntv2_transformation_plugin/grids | |
# 1) Calculation Nr. 1 - GK4 to UTM 32 (Bavaria) | |
# Forward Transformation: GK4 to UTM32 mit offiziellem Testdatensatz --> erfolgreich 100% korrekt (siehe .csv Datei) | |
cs2cs +proj=tmerc +lat_0=0 +lon_0=12 +k=1.000000 +x_0=4500000 +y_0=0 +ellps=bessel +datum=potsdam +nadgrids=$GRID_PATH/BY_KANU.gsb +units=m +no_defs +to +proj=utm +zone=32 +ellps=GRS80 +units=m +no_defs +nadgrids=@null <<EOF | |
4285567.56 5551117.44 | |
4310409.60 5525223.81 | |
4310409.60 5551117.44 | |
4335251.64 5525223.81 | |
4335251.64 5551117.44 | |
4360093.68 5266287.43 | |
4360093.68 5292181.07 | |
4360093.68 5318074.70 | |
4360093.68 5343968.34 | |
4360093.68 5369861.98 | |
4360093.68 5499330.17 | |
4360093.68 5525223.81 | |
4360093.68 5551117.44 | |
4360093.68 5577011.08 | |
4384935.71 5292181.07 | |
4384935.71 5318074.70 | |
4384935.71 5343968.34 | |
4384935.71 5369861.98 | |
4384935.71 5447542.89 | |
4384935.71 5473436.53 | |
4384935.71 5499330.17 | |
4384935.71 5525223.81 | |
4384935.71 5551117.44 | |
4384935.71 5577011.08 | |
4409777.75 5266287.43 | |
4409777.75 5292181.07 | |
4409777.75 5318074.70 | |
4409777.75 5343968.34 | |
4409777.75 5369861.98 | |
4409777.75 5395755.62 | |
4409777.75 5421649.26 | |
4409777.75 5447542.89 | |
4409777.75 5473436.53 | |
4409777.75 5499330.17 | |
4409777.75 5525223.81 | |
4409777.75 5551117.44 | |
4409777.75 5577011.08 | |
4434619.79 5266287.43 | |
4434619.79 5292181.07 | |
4434619.79 5318074.70 | |
4434619.79 5343968.34 | |
4434619.79 5369861.98 | |
4434619.79 5395755.62 | |
4434619.79 5421649.26 | |
4434619.79 5447542.89 | |
4434619.79 5473436.53 | |
4434619.79 5499330.17 | |
4434619.79 5525223.81 | |
4434619.79 5551117.44 | |
4434619.79 5577011.08 | |
4459461.83 5266287.43 | |
4459461.83 5292181.07 | |
4459461.83 5318074.70 | |
4459461.83 5343968.34 | |
4459461.83 5369861.98 | |
4459461.83 5395755.62 | |
4459461.83 5421649.26 | |
4459461.83 5447542.89 | |
4459461.83 5473436.53 | |
4459461.83 5499330.17 | |
4459461.83 5525223.81 | |
4459461.83 5551117.44 | |
4459461.83 5577011.08 | |
4484303.87 5292181.07 | |
4484303.87 5318074.70 | |
4484303.87 5343968.34 | |
4484303.87 5369861.98 | |
4484303.87 5395755.62 | |
4484303.87 5421649.26 | |
4484303.87 5447542.89 | |
4484303.87 5473436.53 | |
4484303.87 5499330.17 | |
4484303.87 5525223.81 | |
4484303.87 5551117.44 | |
4484303.87 5577011.08 | |
4509145.91 5292181.07 | |
4509145.91 5318074.70 | |
4509145.91 5343968.34 | |
4509145.91 5369861.98 | |
4509145.91 5395755.62 | |
4509145.91 5421649.26 | |
4509145.91 5447542.89 | |
4509145.91 5473436.53 | |
4509145.91 5499330.17 | |
4509145.91 5525223.81 | |
4509145.91 5551117.44 | |
4533987.95 5292181.07 | |
4533987.95 5318074.70 | |
4533987.95 5343968.34 | |
4533987.95 5369861.98 | |
4533987.95 5395755.62 | |
4533987.95 5421649.26 | |
4533987.95 5447542.89 | |
4533987.95 5473436.53 | |
4533987.95 5499330.17 | |
4533987.95 5525223.81 | |
4558829.98 5292181.07 | |
4558829.98 5318074.70 | |
4558829.98 5343968.34 | |
4558829.98 5369861.98 | |
4558829.98 5395755.62 | |
4558829.98 5421649.26 | |
4558829.98 5447542.89 | |
4583672.02 5369861.98 | |
4583672.02 5395755.62 | |
4583672.02 5421649.26 | |
4583672.02 5447542.89 | |
4608514.06 5395755.62 | |
4608514.06 5421649.26 | |
4633356.10 5395755.62 | |
EOF | |
# 2) Calculation Nr. 2 - UTM32 to GK4 (Bavaria) | |
# Reverse Transformation: UTM32 zu GK4 mit offiziellem Testdatensatz --> 100% correct, see Testpunkte_Echtumstellung.csv for details | |
cs2cs +proj=utm +zone=32 +ellps=GRS80 +units=m +no_defs +nadgrids=@null +to +proj=tmerc +lat_0=0 +lon_0=12 +k=1.000000 +x_0=4500000 +y_0=0 +ellps=bessel +datum=potsdam +nadgrids=$GRID_PATH/BY_KANU.gsb +units=m +no_defs <<EOF | |
500274.08 5545039.34 | |
526110.77 5520184.44 | |
525073.88 5546035.85 | |
550913.89 5521177.54 | |
549876.74 5547032.18 | |
585895.11 5263580.2 | |
584893.92 5289440.06 | |
583889.23 5315299.78 | |
582880.73 5341159.12 | |
581868.45 5367018.58 | |
576753.92 5496312.61 | |
575720.28 5522170.91 | |
574683.09 5548028.88 | |
573642.39 5573886.82 | |
609705.16 5290402.18 | |
608700.22 5316265.42 | |
607691.66 5342128.44 | |
606679.36 5367991.36 | |
603620.64 5445578.77 | |
602593.97 5471440.86 | |
601563.72 5497302.7 | |
600530 5523164.33 | |
599492.65 5549025.73 | |
598451.74 5574887.03 | |
635521.34 5265497.95 | |
634519.82 5291364.65 | |
633514.64 5317231.37 | |
632505.89 5343097.71 | |
631493.37 5368964.25 | |
630477.26 5394830.38 | |
629457.54 5420696.33 | |
628434.17 5446562.13 | |
627407.33 5472427.69 | |
626376.87 5498293.07 | |
625342.98 5524157.93 | |
624305.38 5550022.78 | |
623264.25 5575887.39 | |
660339.25 5266456.99 | |
659337.73 5292327.37 | |
658332.5 5318197.62 | |
657323.6 5344067.7 | |
656310.86 5369937.46 | |
655294.51 5395807.04 | |
654274.54 5421676.58 | |
653251.1 5447545.71 | |
652223.95 5473414.8 | |
651193.24 5499283.37 | |
650159.03 5525151.83 | |
649121.23 5551019.9 | |
648080 5576887.92 | |
685160.73 5267416.12 | |
684159.11 5293290.24 | |
683153.71 5319164.13 | |
682144.4 5345037.65 | |
681131.59 5370910.98 | |
680114.98 5396784.05 | |
679094.91 5422656.9 | |
678071.17 5448529.63 | |
677043.77 5474401.98 | |
676012.74 5500274.03 | |
674978.31 5526145.88 | |
673940.4 5552017.4 | |
672898.91 5577888.71 | |
708983.99 5294253.37 | |
707978.22 5320130.67 | |
706968.76 5346007.81 | |
705955.62 5371884.64 | |
704938.88 5397761.22 | |
703918.48 5423637.54 | |
702894.42 5449513.57 | |
701866.92 5475389.49 | |
700835.65 5501264.96 | |
699801.05 5527140.03 | |
698762.7 5553015.01 | |
697721.03 5578889.76 | |
733812.41 5295216.65 | |
732806.29 5321097.54 | |
731796.49 5346978.14 | |
730783.03 5372858.52 | |
729765.96 5398738.62 | |
728745.33 5424618.43 | |
727721.04 5450497.91 | |
726693.15 5476377.2 | |
725661.76 5502255.93 | |
724626.77 5528134.54 | |
723588.43 5554012.95 | |
758644.12 5296180.26 | |
757637.66 5322064.72 | |
756627.77 5347948.88 | |
755613.94 5373832.81 | |
754596.56 5399716.35 | |
753575.57 5425599.6 | |
752550.91 5451482.54 | |
751522.75 5477365.28 | |
750491.03 5503247.48 | |
749455.89 5529129.33 | |
783479.35 5297144.14 | |
782472.57 5323032.24 | |
781462.09 5348919.9 | |
780448.15 5374807.31 | |
779430.42 5400694.32 | |
778409.1 5426581.11 | |
777384.26 5452467.36 | |
805285.73 5375782 | |
804267.72 5401672.45 | |
803246.06 5427562.75 | |
802220.91 5453452.53 | |
829108.23 5402651.06 | |
828086.29 5428544.58 | |
853952.06 5403629.88 | |
EOF | |
# 3) Calculation Nr. 3 - GK4 to UTM 32 (Bavaria, lower precision with BETA2007) | |
# | |
# Result: 279488.01 5654871.71 0.00 | |
# Prüfsumme (adv pdf) 279488.01 5654871.71 0.00 | |
cs2cs +proj=tmerc +lat_0=0 +lon_0=6 +k=1.000000 +x_0=2500000 +y_0=0 +ellps=bessel +units=m +nadgrids=$GRID_PATH/BETA2007.gsb +to +proj=utm +ellps=GRS80 +zone=32 +nadgrids=@null <<EOF | |
2490000.00 5652000.00 | |
EOF | |
# 4) Calculation Nr. 4 - GK2 to UTM 32 (Germany, using submeter-precition BETA2007) | |
# tested against official testing points provided by ADV Germany | |
cs2cs +proj=tmerc +lat_0=0 +lon_0=6 +k=1.000000 +x_0=2500000 +y_0=0 +ellps=bessel +units=m +nadgrids=$GRID_PATH/BETA2007.gsb +to +proj=utm +ellps=GRS80 +zone=32 +nadgrids=@null <<EOF | |
2598417.333 5930677.98 | |
2678509.792 5890320.495 | |
2601895.025 5761935.672 | |
2606412.76 5654454.412 | |
2604044.332 5596917.812 | |
2611430.565 5467270.624 | |
2590840.679 5466891.207 | |
2615145.447 5281966.148 | |
2643120.946 5551463.861 | |
2661073.96 6067930.994 | |
EOF | |
# Comparison with official testing points: | |
# UTM32-E UTM32-N deltaX deltaY | |
# 399340.60 5928794.18 -0,002 0,002 | |
# 477621.72 5885134.57 -0,002 0,003 | |
# 395783.50 5760119.72 0,003 0,005 | |
# 395892.87 5652585.90 0,005 0,005 | |
# 391195.03 5595215.13 0,000 0,002 | |
# 393381.12 5465427.35 -0,002 -0,001 | |
# 372799.65 5465865.76 0,002 0,005 | |
# 389829.27 5280195.60 0,002 -0,001 | |
# 428391.21 5548246.77 0,001 0,003 | |
# 467726.90 6063191.97 0,004 -0,004 | |
# 5) Calculation Nr. 5 - GK 3 to UTM 32 using BETA2007 | |
# Testing points --> last column is error in meters | |
# DE_DHDN_3GK3 3399371,19 5930724,531 ETRS89_UTM32 399340,6019 5928794,178 -0,0019 | |
# DE_DHDN_3GK3 3599586,686 5819391,66 ETRS89_UTM32 599474,9342 5817502,627 -0,0042 | |
# DE_DHDN_3GK3 3477684,063 5887048,677 ETRS89_UTM32 477621,7225 5885134,567 -0,0025 | |
# DE_DHDN_3GK3 3395815,327 5761982,907 ETRS89_UTM32 395783,4969 5760119,715 0,0031 | |
# DE_DHDN_3GK3 3395925,872 5654406,809 ETRS89_UTM32 395892,8652 5652585,895 -0,0052 | |
# DE_DHDN_3GK3 3391226,59 5597013,366 ETRS89_UTM32 391195,0301 5595215,128 -9,99999E-05 | |
# DE_DHDN_3GK3 3393414,08 5467174,397 ETRS89_UTM32 393381,1216 5465427,351 -0,0016 | |
# DE_DHDN_3GK3 3372824,499 5467612,907 ETRS89_UTM32 372799,6479 5465865,755 0,0021 | |
# DE_DHDN_3GK3 3389860,774 5281869,239 ETRS89_UTM32 389829,2676 5280195,601 0,0024 | |
# DE_DHDN_3GK3 3597874,422 5411397,512 ETRS89_UTM32 597759,8986 5409672,24 0,0014 | |
# DE_DHDN_3GK3 3428437,613 5550026,645 ETRS89_UTM32 428391,2092 5548246,767 0,0008 | |
# DE_DHDN_3GK3 3571307,006 5655705,338 ETRS89_UTM32 571204,5633 5653882,477 -0,0033 | |
# DE_DHDN_3GK3 3595392,782 6023387,96 ETRS89_UTM32 595286,0444 6021417,377 -0,0044 | |
# DE_DHDN_3GK3 3615881,001 5940351,728 ETRS89_UTM32 615764,364 5938413,819 -0,004 | |
cs2cs +proj=tmerc +lat_0=0 +lon_0=9 +k=1.000000 +x_0=3500000 +y_0=0 +ellps=bessel +datum=potsdam +nadgrids=$GRID_PATH/BETA2007.gsb +units=m +no_defs +to +proj=utm +ellps=GRS80 +zone=32 +nadgrids=@null <<EOF | |
3399371.19 5930724.531 | |
3599586.686 5819391.66 | |
3477684.063 5887048.677 | |
3395815.327 5761982.907 | |
3395925.872 5654406.809 | |
3391226.59 5597013.366 | |
3393414.08 5467174.397 | |
3372824.499 5467612.907 | |
3389860.774 5281869.239 | |
3597874.422 5411397.512 | |
3428437.613 5550026.645 | |
3571307.006 5655705.338 | |
3595392.782 6023387.96 | |
3615881.001 5940351.728 | |
EOF | |
# 6) Calculation Nr. 6 - GK 4 to UTM 32 using BETA2007 | |
# DE_DHDN_3GK4 4395886,919 5819485,694 ETRS89_UTM32 | |
# DE_DHDN_3GK4 4377657,795 5411879,84 ETRS89_UTM32 | |
# DE_DHDN_3GK4 4360897,154 5657085,679 ETRS89_UTM32 | |
# DE_DHDN_3GK4 4400271,506 6023480,198 ETRS89_UTM32 | |
# DE_DHDN_3GK4 4417225,999 5939654,081 ETRS89_UTM32 | |
# DE_DHDN_3GK4 4368411,664 5931484,902 ETRS89_UTM32 | |
# DE_DHDN_3GK4 4433598,022 5926410,007 ETRS89_UTM32 | |
# DE_DHDN_3GK4 4596699,815 5960328,297 ETRS89_UTM32 | |
# DE_DHDN_3GK4 4431385,772 5848536,122 ETRS89_UTM32 | |
# DE_DHDN_3GK4 4568999,834 5739119,061 ETRS89_UTM32 | |
# DE_DHDN_3GK4 4399252,521 5704414,901 ETRS89_UTM32 | |
# DE_DHDN_3GK4 4390237,958 5606306,172 ETRS89_UTM32 | |
# DE_DHDN_3GK4 4385715,06 5504412,339 ETRS89_UTM32 | |
# DE_DHDN_3GK4 4391285,797 5291109,755 ETRS89_UTM32 | |
# DE_DHDN_3GK4 4602499,566 5615431,38 ETRS89_UTM32 | |
# DE_DHDN_3GK4 4608008,856 5689725,987 ETRS89_UTM32 | |
# DE_DHDN_3GK4 4607459,254 5769472,054 ETRS89_UTM32 | |
# DE_DHDN_3GK4 4595844,51 5893520,179 ETRS89_UTM32 | |
# DE_DHDN_3GK4 4599624,347 5817537,418 ETRS89_UTM32 | |
# DE_DHDN_3GK4 4582806,458 5437104,667 ETRS89_UTM32 | |
# DE_DHDN_3GK4 4571363,305 5279411,441 ETRS89_UTM32 | |
cs2cs +proj=tmerc +lat_0=0 +lon_0=12 +k=1.000000 +x_0=4500000 +y_0=0 +ellps=bessel +datum=potsdam +nadgrids=$GRID_PATH/BETA2007.gsb +units=m +no_defs +to +proj=utm +ellps=GRS80 +zone=32 +nadgrids=@null <<EOF | |
4395886.919 5819485.694 | |
4377657.795 5411879.84 | |
4360897.154 5657085.679 | |
4400271.506 6023480.198 | |
4417225.999 5939654.081 | |
4368411.664 5931484.902 | |
4433598.022 5926410.007 | |
4596699.815 5960328.297 | |
4431385.772 5848536.122 | |
4568999.834 5739119.061 | |
4399252.521 5704414.901 | |
4390237.958 5606306.172 | |
4385715.06 5504412.339 | |
4391285.797 5291109.755 | |
4602499.566 5615431.38 | |
4608008.856 5689725.987 | |
4607459.254 5769472.054 | |
4595844.51 5893520.179 | |
4599624.347 5817537.418 | |
4582806.458 5437104.667 | |
4571363.305 5279411.441 | |
EOF | |
# 7) Calculation Nr. 7 - GK 3 to UTM 32 using HETA2010 (DE-HE) | |
cs2cs +proj=tmerc +lat_0=0 +lon_0=9 +k=1.000000 +x_0=3500000 +y_0=0 +ellps=bessel +datum=potsdam +nadgrids=$GRID_PATH/HETA2010.gsb +units=m +no_defs +to +proj=utm +ellps=GRS80 +zone=32 +nadgrids=@null <<EOF | |
3412360 5547670 | |
3426390 5586650 | |
3448730 5638860 | |
3465100 5660520 | |
3469090 5682630 | |
3501470 5709400 | |
3534060 5724990 | |
3544860 5703930 | |
3563050 5698900 | |
3586530 5673120 | |
3585270 5655170 | |
3573210 5635820 | |
3574280 5615880 | |
3561180 5584950 | |
3547700 5567250 | |
3536790 5551050 | |
3510650 5550390 | |
3503750 5539630 | |
3510930 5511800 | |
3508780 5486160 | |
3486920 5473190 | |
3468710 5487010 | |
3453620 5505140 | |
3452890 5527660 | |
3419570 5538230 | |
3444270 5545110 | |
3431550 5565790 | |
3483880 5674870 | |
3551160 5684760 | |
3476290 5498610 | |
3530000 5720000 | |
3480000 5690000 | |
3500000 5680000 | |
3580000 5660000 | |
3500000 5640000 | |
3480000 5620000 | |
3520000 5620000 | |
3440000 5610000 | |
3460000 5600000 | |
3500000 5600000 | |
3540000 5600000 | |
3570000 5600000 | |
3480000 5580000 | |
3520000 5580000 | |
3500000 5560000 | |
3420000 5550000 | |
3530000 5550000 | |
3500000 5520000 | |
3500000 5500000 | |
3490000 5480000 | |
EOF | |
### Below is a collection of helpful proj strings, not all of them are tested into detail. | |
### originally from: http://www.osmwms.de/html/leistungsbeschreibung.html | |
json_reference_doc = [{"epsg":4326, | |
"description":"WGS84, World Geodetic System", | |
"proj4_string":"+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs" | |
}, | |
{"epsg":900913, | |
"description":"Kugelprojektion des Suchmaschinenbetreibers Google", | |
"proj4_string":"+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +wktext +no_defs" | |
}, | |
{"epsg":3395, | |
"description":"Wie EPSG:900913,nur mit Ellipsoid", | |
"proj4_string":"+proj=merc +lon_0=0 +k=1.000000 +x_0=0 +y_0=0 +ellps=WGS84 +datum=WGS84 +units=m +no_defs" | |
}, | |
{"epsg":21781, | |
"description":"Schweiz (CH1903 / LV03)", | |
"proj4_string":"+proj=somerc +lat_0=46.95240555555556 +lon_0=7.439583333333333 +x_0=600000 +y_0=200000 +ellps=bessel +towgs84=674.374,15.056,405.346,0,0,0,0 +units=m +no_defs" | |
}, | |
{"epsg":25831, | |
"description":"UTM-Zone 31", | |
"proj4_string":"+proj=utm +zone=31 +ellps=GRS80 +units=m +no_defs" | |
}, | |
{"epsg":25832, | |
"description":"Lagestatus 489, UTM-Zone 32 (ST)", | |
"proj4_string":"+proj=utm +zone=32 +ellps=GRS80 +units=m +no_defs" | |
}, | |
{"epsg":25833, | |
"description":"Lagestatus 489, UTM-Zone 33 (ST)", | |
"proj4_string":"+proj=utm +zone=33 +ellps=GRS80 +units=m +no_defs" | |
}, | |
{"epsg":31466, | |
"description":"DHDN – Deutsches Hauptdreiecksnetz, Zone 2", | |
"proj4_string":"+proj=tmerc +lat_0=0 +lon_0=6 +k=1.000000 +x_0=2500000 +y_0=0 +ellps=bessel +datum=potsdam +nadgrids=BeTA2007 +units=m +no_defs" | |
}, | |
{"epsg":31467, | |
"description":"DHDN – Deutsches Hauptdreiecksnetz, Zone 3 ", | |
"proj4_string":"+proj=tmerc +lat_0=0 +lon_0=9 +k=1.000000 +x_0=3500000 +y_0=0 +ellps=bessel +datum=potsdam +nadgrids=BeTA2007 +units=m +no_defs" | |
}, | |
{"epsg":31468, | |
"description":"DHDN – Deutsches Hauptdreiecksnetz, Zone 4", | |
"proj4_string":"+proj=tmerc +lat_0=0 +lon_0=12 +k=1.000000 +x_0=4500000 +y_0=0 +ellps=bessel +datum=potsdam +nadgrids=BeTA2007 +units=m +no_defs" | |
}, | |
{"epsg":3398, | |
"description":"Lagestatus 110, 4.GK-Streifen", | |
"proj4_string":"+proj=tmerc +lat_0=0 +lon_0=12 +k=1.000000 +x_0=4500000 +y_0=0 +ellps=bessel +towgs84=612.4,77,440.2,-0.054,0.057,-2.797,2.55 +datum=potsdam +units=m +no_defs" | |
}, | |
{"epsg":3399, | |
"description":"Lagestatus 110, 5.GK-Streifen", | |
"proj4_string":"+proj=tmerc +lat_0=0 +lon_0=15 +k=1.000000 +x_0=5500000 +y_0=0 +ellps=bessel +towgs84=612.4,77,440.2,-0.054,0.057,-2.797,2.55 +datum=potsdam +units=m +no_defs" | |
}, | |
{"epsg":3396, | |
"description":"Lagestatus 120, 3.GK-Streifen", | |
"proj4_string":"+proj=tmerc +lat_0=0 +lon_0=9 +k=1.000000 +x_0=3500000 +y_0=0 +ellps=bessel +towgs84=599.4,72.4,419.2,-0.063,-0.022,-2.723,6.46 +datum=potsdam +units=m +no_defs" | |
}, | |
{"epsg":3397, | |
"description":"Lagestatus 120, 4.GK-Streifen", | |
"proj4_string":"+proj=tmerc +lat_0=0 +lon_0=12 +k=1.000000 +x_0=4500000 +y_0=0 +ellps=bessel +towgs84=599.4,72.4,419.2,-0.063,-0.022,-2.723,6.46 +datum=potsdam +units=m +no_defs" | |
}, | |
{"epsg":2397, | |
"description":"Lagestatus 150, 3.GK-Streifen", | |
"proj4_string":"+proj=tmerc +lat_0=0 +lon_0=9 +k=1.000000 +x_0=3500000 +y_0=0 +ellps=krass +towgs84=24.9,-126.4,-93.2,-0.063,-0.247,-0.041,1.01 +units=m +no_defs" | |
}, | |
{"epsg":2398, | |
"description":"Lagestatus 150, 4.GK-Streifen", | |
"proj4_string":"+proj=tmerc +lat_0=0 +lon_0=12 +k=1.000000 +x_0=4500000 +y_0=0 +ellps=krass +towgs84=24.9,-126.4,-93.2,-0.063,-0.247,-0.041,1.01 +units=m +no_defs" | |
}, | |
{"epsg":2399, | |
"description":"Lagestatus 150, 5.GK-Streifen", | |
"proj4_string":"+proj=tmerc +lat_0=0 +lon_0=15 +k=1.000000 +x_0=5500000 +y_0=0 +ellps=krass +towgs84=24.9,-126.4,-93.2,-0.063,-0.247,-0.041,1.01 +units=m +no_defs" | |
}, | |
{"epsg":2166, | |
"description":"Pulkovo 1942(83), 3.GK-Streifen", | |
"proj4_string":"+proj=tmerc +lat_0=0 +lon_0=9 +k=1.000000 +x_0=3500000 +y_0=0 +ellps=krass +towgs84=24,-123,-94,0.02,-0.25,-0.13,1.1 +units=m +no_defs" | |
}, | |
{"epsg":2167, | |
"description":"Pulkovo 1942(83),4.GK-Streifen", | |
"proj4_string":"+proj=tmerc +lat_0=0 +lon_0=12 +k=1.000000 +x_0=4500000 +y_0=0 +ellps=krass +towgs84=24,-123,-94,0.02,-0.25,-0.13,1.1 +units=m +no_defs" | |
}, | |
{"epsg":2168, | |
"description":"Pulkovo 1942(83),5.GK-Streifen", | |
"proj4_string":"+proj=tmerc +lat_0=0 +lon_0=15 +k=1.000000 +x_0=5500000 +y_0=0 +ellps=krass +towgs84=24,-123,-94,0.02,-0.25,-0.13,1.1 +units=m +no_defs" | |
}, | |
{"epsg":45832, | |
"description":"Lagestatus 489, (Rechtswert mit vorangestellter 2 -- Brandenburg)", | |
"proj4_string":"+proj=tmerc +lat_0=0 +lon_0=9 +k=0.9996 +x_0=2500000 +y_0=0 +ellps=GRS80 +units=m +no_defs" | |
}, | |
{"epsg":45833, | |
"description":"Lagestatus 489, (Rechtswert mit vorangestellter 3 -- Brandenburg)", | |
"proj4_string":"+proj=tmerc +lat_0=0 +lon_0=15 +k=0.9996 +x_0=3500000 +y_0=0 +ellps=GRS80 +units=m +no_defs" | |
}, | |
{"epsg":35831, | |
"description":"Lagestatus 489 (Rechtswert mit vorangestellter 31 ganz Deutschland außer BB, ST)", | |
"proj4_string":"+proj=tmerc +lat_0=0 +lon_0=3 +k=0.9996 +x_0=31500000 +y_0=0 +ellps=GRS80 +units=m +no_defs" | |
}, | |
{"epsg":35832, | |
"description":"Lagestatus 489 (Rechtswert mit vorangestellter 32 ganz Deutschland außer BB, ST)", | |
"proj4_string":"+proj=tmerc +lat_0=0 +lon_0=9 +k=0.9996 +x_0=32500000 +y_0=0 +ellps=GRS80 +units=m +no_defs" | |
}, | |
{"epsg":35833, | |
"description":"Lagestatus 489 (Rechtswert mit vorangestellter 33 ganz Deutschland außer BB, ST)", | |
"proj4_string":"+proj=tmerc +lat_0=0 +lon_0=15 +k=0.9996 +x_0=33500000 +y_0=0 +ellps=GRS80 +units=m +no_defs" | |
} | |
] |
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment