Single input

Single input usage.

from geodesiclinestogis import GeodesicLine2Gisfile

lons_lats: input coordinates. (start longitude, start latitude, end longitude, end latitude)

lons_lats = (-3.6,40.5,-118.4,33.9)

Folder path to store output file and filename:

folderpath = '/tmp'

layername = "geodesicline"

Create object. You can pass two parameters: - antimeridian: [True, False] to solve antimeridian problem (default is True). - prints: [True, False] print output messages (default is True).

gtg = GeodesicLine2Gisfile()

Launch computations. You can pass two parameter: - lons_lats: input coords returned by gcComp. - km_pts: compute one point each n km (default is 20 km)

cd = gtg.gdlComp(lons_lats, km_pts=30)

Dump geodetic line coords to Linestring Feature and store in a GIS file.

Output formats: “ESRI Shapefile” (default), “GeoJSON”

# shapefile output
gtg.gdlToGisFile(cd, folderpath, layername)

# geojson output
gtg.gdlToGisFile(cd, folderpath, layername, fmt="GeoJSON")

Multiple input

Multiple input usage.

from geodesiclinestogis import GeodesicLine2Gisfile

data = [
        (-6.,37.,-145.,11.),
        (-150.,37.,140.,11.),
        (-6.,37.,120.,50.),
        (-3.6,40.5,-118.4,33.9),
        (-118.4,33.9,139.8,35.5),
        (-118.4,33.9,104.,1.35),
        (-118.4,33.9,151.,-33.9),
        (-20.4,33.9,178.,-33.9)
    ]

folderpath = "/tmp/geod_line"

layername = "geodesicline"

gtg = GeodesicLine2Gisfile()

gtg.gdlToGisFileMulti(data, folderpath, layername)

Numpy array (multiple) input

Numpy array input usage.

import numpy as np
from geodesiclinestogis import GeodesicLine2Gisfile

data_ = [
        (-6.,37.,-145.,11.),
        (-150.,37.,140.,11.),
        (-6.,37.,120.,50.),
        (-3.6,40.5,-118.4,33.9),
        (-118.4,33.9,139.8,35.5),
        (-118.4,33.9,104.,1.35),
        (-118.4,33.9,151.,-33.9),
        (-20.4,33.9,178.,-33.9)
    ]

data = np.array(data_)

folderpath = "/tmp/geod_line_np"

layername = "geodesicline"

gtg = GeodesicLine2Gisfile()

gtg.gdlToGisFileMulti(data, folderpath, layername)