Sample program in Python / Spyder

Sample program in Python / Spider

# -*- coding: utf-8 -*-

"""

Python 3.7.11 [MSC v.1927 64 bit (AMD64)]

Created on Mon May 23 09:16:09 2022

/*****************************************************************************/

/* */

/* Template of a Python 3.7.11 [MSC v.1927 64 bit (AMD64)] program as a */

/* Spyder 3.3.6 project for some function calls of the 64 bit Geodetic */

/* Development Kit GeoDLL. */

/* The template can be extended for other GeoDLL function calls if */

/* necessary. */

/*---------------------------------------------------------------------------*/

/* */

/* Required external files: */

/* installation from GeoDLL in directory 'C:\\geodll\geodll\' */

/* 1. "C:\\geodll\geodll\GEODLL64.DLL" */

/* 2. "C:\\geodll\geodll\GEODLLBN.BIN" */

/* 3. Python installation: https://www.anaconda.com/products/distribution */

/* 4. Spyder installation: https://anaconda.org/anaconda/spyder */

/* */

/* Online user manual with all function declarations: */

/* https://www.killetsoft.de/h_geodll_e/handbuch_e.htm in English language */

/* https://www.killetsoft.de/h_geodll_e/handbuch_d.htm in German language */

/* */

/* Included help files with all function declarations see GeoDLL download: */

/* GeoDLL_e.chm in English language */

/* GeoDLL_e.chm in German language */

/* */

/* Author of this source: */

/* Dipl.-Ing. Andreas Neukirchen M.A. */

/* Ziviltechnikerkanzlei Baurat h.c. Dipl.-Ing. Andreas Neukirchen M.A. */

/* A-2531 Gaaden, Skodagasse 14 */

/* a.neukirchen@neuk.at */

/* */

/* Author of GeoDLL: */

/* Dipl.-Ing. Fred Killet */

/* Killet GeoSoftware Ing.-GbR */

/* https://www.killetsoft.de */

/* */

/*****************************************************************************/

"""

import ctypes

#------------------------------------------------------------------------------

# used class prototype

#------------------------------------------------------------------------------

class clsgeopoint:

def __init__(self,x,y,z,epsg):

self.X_dbl = ctypes.c_double(x)

self.Y_dbl = ctypes.c_double(y)

self.H_dbl = ctypes.c_double(z)

self.EPSG = ctypes.c_ushort(epsg)

self.X_float = self.X_dbl.value

self.Y_float = self.Y_dbl.value

self.H_float = self.H_dbl.value

#------------------------------------------------------------------------------

# used functions prototypes

#------------------------------------------------------------------------------

def c_str(string):

#Create ctypes char * from a Python string.

return ctypes.c_char_p(string.encode('utf-8'))

def c_dbl(float2double):

#Create ctypes double from a Python float.

return ctypes.c_double(float2double)

def c_ushort(int2ushort):

#Create ctypes ushort from a Python int.

return ctypes.c_ushort(int2ushort)

#------------------------------------------------------------------------------

# loading GeoDLL environment

#------------------------------------------------------------------------------

geodll = ctypes.CDLL("C:\\geodll\geodll\GEODLL64.DLL")

#print(geodll)

# Set GeoDLL internal language to English

geodll.setlanguage(ctypes.c_ushort(2))

# With no unlockkey set no_unlockkey = 1 otherwise 0

no_unlockkey = 1

# Calling setunlockcode()

# Set the original unlock parameters you got from Killetsoft here:

unlockcode = c_str("123456789-987654321")

username = c_str("MicroModia GmbH")

ppszError = c_str("")

if no_unlockkey != 1:

geodll.setunlockcode.argtypes = [ctypes.c_char_p,ctypes.c_char_p]

nRet = geodll.setunlockcode(unlockcode,username)

# Explanation

print("\n")

print("You can use all features of GeoDLL in your own application initially free.")

print("However, the number of function calls per program run is limited to 100")

print("geodetic calculations. Also in the free test version a message appears")

print("on the screen or is written to the EventLog and a tone sequence is sounded.")

print("\n")

print("After WIN10 PRO DLL installation:")

print("load DLL: geodll = ctypes.CDLL('C:\\geodll\geodll\GEODLL64.DLL')\n")

print("\n")

print("By purchasing a licence the restrictions will be canceled. For unlocking")

print("you get an user name and an unlock code from KilletSoft.")

print("\n")

print("The GeoDLL function groups required by your application must be unlocked")

print("in the initialization part of your program with the function setunlockcode().")

print("The function gets the user name and the unlock code as parameters.")

print("\n")

print("How to do that I will show you here. Of course, these unlock parameters")

print("are fictitious and do not work. Therefore the error message and following")

print("tone sequence are correct in this case!")

print("\n")

print("nRet = geodll.setunlockcode(\"123456789-987654321\",\"MicroModia GmbH\");")

print("if(!= nRet)")

print(" geodll.geterrorcode(ppszError) with stringError = "" and ")

print("ppszError = ctypes.c_char_p(stringError.encode('utf-8'))")

print("strError = ppszError.value.decode('utf-8')\n")

if no_unlockkey != 1:

if nRet == 1:

print("GeoDLL functions unlocked!")

#print("Funktionsgruppe Koordinatentransformation freigeschaltet!")

print(username.value.decode("utf-8"))

else:

ppszError = c_str("")

#geodll.geterrorcode(ctypes.byref(ppszError))

#print("type of ppszError: " + str(type(ppszError)))

#print(ppszError)

strp = ppszError.value.decode("utf-8")

print("GeoDLL functions remain locked!")

print("Error Message from GeodDLL: " + strp)

# Wait for Key Press or Mouse Click.

input("Press Enter to continue...")

#------------------------------------------------------------------------------

# Coordinate Transformations with the function coordtrans()

#------------------------------------------------------------------------------

#Explanation

print("\n")

print("Now we start transforming coordinates. Here only the use of the GeoDLL")

print("function coordtrans() is demonstrated. But you can also use any other")

print("function for coordinate transformations available in GeoDLL. The required")

print("Coordinate Systems and Reference Systems you find sorted by countries in")

print("the User Manual GeoDLL_e.chm, chapter \"Coodinate Transformations\", in the")

print("file README.DOC or in the Online User Manual.")

print("\n")

# Gauss-Kueger, DHDN --> UTM, ETRS89

# origin

dCoordXQ = c_dbl(2531418.55) # double dCoordXQ

dCoordYQ = c_dbl(5695935.71) # double dCoordYQ

pszDummy = c_str("") # char *pszCoordQ

nCoordSysQ = c_ushort(2) # unsigned short Gauss-Krueger

nRefSysQ = c_ushort(17) # unsigned short DHDN

# destination

dCoordXZ = c_dbl(0.0) # double *dCoordXZ

dCoordYZ = c_dbl(0.0) # double *dCoordYZ

ppszDummy = c_str("") # char **pszCoordZ

nCoordSysZ = c_ushort(3) # unsigned short UTM

nRefSysZ = c_ushort(4) # unsigned short ETRS89

nStripZ = c_ushort(0) # unsigned short

nRet = geodll.coordtrans(\

dCoordXQ,\

dCoordYQ,\

pszDummy,\

nCoordSysQ,\

nRefSysQ,\

ctypes.byref(dCoordXZ),\

ctypes.byref(dCoordYZ),\

ctypes.byref(ppszDummy),\

nCoordSysZ,\

nRefSysZ,\

nStripZ)

# Output of the results.

if nRet == 1:

print("Gauss-Kueger, DHDN --> UTM, ETRS89")

str1 = "{0} {1:12.6f}".format("Gauss-Krueger Easting = ",dCoordXQ.value)

str2 = "{0} {1:12.6f}".format("Gauss-Krueger Northing = ",dCoordYQ.value)

str3 = "{0} {1:12.6f}".format("UTM Easting = ",dCoordXZ.value)

str4 = "{0} {1:12.6f}".format("UTM Easting = ",dCoordYZ.value)

print(str1)

print(str2)

print(str3)

print(str4)

else:

ppszError = c_str("")

geodll.geterrorcode(ctypes.byref(ppszError))

strp = ppszError.value.decode("utf-8")

print("GeoDLL coordinate transformation failed!")

print("Error: " + strp)

print("\n")

#input("Press Enter to continue...")

# UTM, ETRS89 -- > UTMRef (Centimeter Grid), ED50 (Europe)

# Note:

# GeoDLL is able to allocate memory for returnable PSZ parameters internally

# and makes them accessible to the calling program by global stored addresses.

# In the calling program then neither memory must allocated nor memory must

# released later. For more information see function setstringallocate().

dCoordXQ = dCoordXZ

dCoordYQ = dCoordYZ

pszDummy = c_str("") # char *pszCoordQ

nCoordSysQ = c_ushort(3) # unsigned short UTM

nRefSysQ = c_ushort(4) # unsigned short ETRS89

# destination

dCoordXZ = c_dbl(0.0) # double *dCoordXZ

dCoordYZ = c_dbl(0.0) # double *dCoordYZ

ppszCoordZ = c_str("") # char **pszCoordZ

nCoordSysZ = c_ushort(661) # unsigned short UTMref

nRefSysZ = c_ushort(2) # unsigned short ED50

nStripZ = c_ushort(0) # unsigned short

nRet = geodll.coordtrans(\

dCoordXQ,\

dCoordYQ,\

pszDummy,\

nCoordSysQ,\

nRefSysQ,\

ctypes.byref(dCoordXZ),\

ctypes.byref(dCoordYZ),\

ctypes.byref(ppszCoordZ),\

nCoordSysZ,\

nRefSysZ,\

nStripZ)

# Output of the results.

if nRet == 1:

print("UTM, ETRS89 -- > UTMRef (Centimeter Grid), ED50 (Europe)")

print("UTMRef grid coordinate = ",ppszCoordZ.value.decode("utf-8"))

else:

ppszError = c_str("")

geodll.geterrorcode(ctypes.byref(ppszError))

strp = ppszError.value.decode("utf-8")

print("GeoDLL coordinate transformation failed!")

print("Error: " + strp)

print("\n")

#input("Press Enter to continue...")

# UTMRef (Centimeter Grid), ED50 (Europe) --> Geographic Coordinates, WGS84

dCoordXQ = c_dbl(0.0)

dCoordYQ = c_dbl(0.0)

pszCoordQ = c_str(ppszCoordZ.value.decode('utf-8')) # hard copy

nCoordSysQ = c_ushort(661) # unsigned short UTMref

nRefSysQ = c_ushort(2) # unsigned short ED50

# destination

dCoordXZ = c_dbl(0.0) # double *dCoordXZ

dCoordYZ = c_dbl(0.0) # double *dCoordYZ

ppszDummy = c_str("") # char **pszCoordZ

nCoordSysZ = c_ushort(6) # unsigned short Geo. Coord

nRefSysZ = c_ushort(10) # unsigned short WGS84

nStripZ = c_ushort(0) # unsigned short

nRet = geodll.coordtrans(\

dCoordXQ,\

dCoordYQ,\

pszCoordQ,\

nCoordSysQ,\

nRefSysQ,\

ctypes.byref(dCoordXZ),\

ctypes.byref(dCoordYZ),\

ctypes.byref(ppszDummy),\

nCoordSysZ,\

nRefSysZ,\

nStripZ)

# Output of the results.

if nRet == 1:

print("UTMRef (Centimeter Grid), ED50 (Europe) --> Geographic Coordinates, WGS84")

str3 = "{0} {1:12.9f}".format("Longitude = ",dCoordXZ.value)

str4 = "{0} {1:12.9f}".format("Latitude = ",dCoordYZ.value)

print(str3)

print(str4)

else:

ppszError = c_str("")

geodll.geterrorcode(ctypes.byref(ppszError))

strp = ppszError.value.decode("utf-8")

print("GeoDLL coordinate transformation failed!")

print("Error: " + strp)

print("\n")

input("Press Enter to continue...")

#------------------------------------------------------------------------------

# Coordinate Transformations with the function coordtransepsg()

#------------------------------------------------------------------------------

# Explanation

print("\n")

print("Coordinate Transformations with the function coordtransepsg():")

print("For simplicity, it is possible to perform all required settings for the")

print("parameters Coordinate System, Reference System, Meridian Strip, Meridian")

print("Strip Suffix and 2D/3D Modus by using so-called EPSG codes. EPSG is the")

print("acronym for \"European Petroleum Survey Group Geodesy\". This is a working")

print("group of the European oil and gas exploration companies. The EPSG and")

print("their successors OGP \"International Association of Oil and Gas Producers\"")

print("have made it to their task to build a system with globally unique")

print("identifiers (EPSG codes) for geodetic data, such as Coordinate Reference")

print("Systems, Reference Ellipsoids or Map Projections.")

print("\n")

print("The required EPSG Codes you find in the User Manual GeoDLL_e.chm, chapter")

print("\"Coordinate Transformations\", in the Online Manual or in the EPSG database.")

print("\n")

print("This EPSG Codes for Coordinate Reference Systems are used:")

print("4326: EPSG code for WGS 84 -- WGS84 - World Geodetic System 1984")

print("31254: EPSG code for MGI/Austria GK West")

print("\n")

# use class for conversion to ctypes

coordQ = clsgeopoint(11.3439685833333,47.2602099166667,0.0,4326) # origin coord

dCoordXQ = coordQ.X_dbl

dCoordYQ = coordQ.Y_dbl

dEllHgtQ = coordQ.H_dbl

epsgQ = coordQ.EPSG

sysorigin = "EPSG 4326 WGS84 World Geodetic System 1984"

# or you can do it by hand

#dCoordXQ = ctypes.c_double(11.3439685833333)

#dCoordYQ = ctypes.c_double(47.2602099166667)

#dEllHgtQ = ctypes.c_double(0.0)

#epsgQ = ctypes.c_ushort(4326)

coordZ = clsgeopoint(0.0,0.0,0.0,31254) # initialyze destination coord

dCoordXZ = coordZ.X_dbl

dCoordYZ = coordZ.Y_dbl

dEllHgtZ = coordQ.H_dbl

epsgZ = coordZ.EPSG

# you can do it by hand

#dCoordXZ = ctypes.c_double()

#dCoordYZ = ctypes.c_double()

#dEllHgtZ = ctypes.c_double()

#epsgZ = ctypes.c_ushort(31254)

sysdestination = "EPSG 31254 MGI/GK Austria West"

nRet = geodll.coordtransepsg(dCoordXQ,\

dCoordYQ,\

dEllHgtQ,\

epsgQ,\

ctypes.byref(dCoordXZ),\

ctypes.byref(dCoordYZ),\

ctypes.byref(dEllHgtZ),\

epsgZ)

#Output of the results.

if (nRet == 1):

print("origin: " + sysorigin + " --> ")

print("destination: " + sysdestination)

str1 = "WGS84 Longitude = {0:6.9f}".format(dCoordXQ.value)

str2 = "WGS84 Latitude = {0:6.9f}".format(dCoordYQ.value)

str3 = "MGI GK Austria West Easting = {0:12.3f}".format(dCoordXZ.value)

str4 = "MGI GK Austria West Northing = {0:12.3f}".format(dCoordYZ.value)

print(str1)

print(str2)

print(str3)

print(str4)

else:

ppszError = c_str("")

geodll.geterrorcode(ctypes.byref(ppszError))

strp = ppszError.value.decode("utf-8")

print("GeoDLL coordinate transformation failed!")

print("Error: " + strp)

input("Press Enter to continue...")

print("\n")

print("Coordinate Transformations with the function coordtransepsg():")

print("This EPSG Codes for Coordinate Reference Systems are used:")

print("31466: EPSG code for \"DHDN / 3-deg. Gauss-Krueger zone 2\"")

print(" 4647: EPSG code for \"ETRS89 / UTM zone 32N (zE-N)\"")

print("\n");

# Gauss-Kueger, DHDN --> UTM, ETRS89

# use class for conversion to ctypes

coordQ = clsgeopoint(2531418.55,5695935.71,0.0,31466) # origin coord

dCoordXQ = coordQ.X_dbl

dCoordYQ = coordQ.Y_dbl

dEllHgtQ = coordQ.H_dbl # double dEllHgtQ Ell. Height not usesed

epsgQ = coordQ.EPSG

# initialyze destination coord

coordZ = clsgeopoint(0.0,0.0,0.0,4647)

dCoordXZ = coordZ.X_dbl

dCoordYZ = coordZ.Y_dbl

dEllHgtZ = coordQ.H_dbl # double* dEllHgtZ Ell. Height not usesed

epsgZ = coordZ.EPSG

nRet = geodll.coordtransepsg(dCoordXQ,\

dCoordYQ,\

dEllHgtQ,\

epsgQ,\

ctypes.byref(dCoordXZ),\

ctypes.byref(dCoordYZ),\

ctypes.byref(dEllHgtZ),\

epsgZ)

#Output of the results.

if (nRet == 1):

print("Gauss-Kueger, DHDN --> UTM, ETRS89");

str1 = "Gauss-Krueger Easting = {0:12.6f}".format(dCoordXQ.value)

str2 = "Gauss-Krueger Northing = {0:12.6f}".format(dCoordYQ.value)

str3 = "UTM Easting = {0:12.6f}".format(dCoordXZ.value)

str4 = "UTM Northing = {0:12.6f}".format(dCoordYZ.value)

print(str1)

print(str2)

print(str3)

print(str4)

else:

ppszError = c_str("")

geodll.geterrorcode(ctypes.byref(ppszError))

strp = ppszError.value.decode("utf-8")

print("GeoDLL coordinate transformation failed!")

print("Error: " + strp)

input("Press Enter to continue...")