Skip to content

Instantly share code, notes, and snippets.

@dbuscombe-usgs
Created September 9, 2015 18:27
Show Gist options
  • Save dbuscombe-usgs/9c1626c8523f79f59379 to your computer and use it in GitHub Desktop.
Save dbuscombe-usgs/9c1626c8523f79f59379 to your computer and use it in GitHub Desktop.
creating a trackline kml in PyHum using non-smoothed easting/northings
from scipy.io import loadmat
import simplekml
import pyproj
import numpy as np
my_meta_file = '/home/dbuscombe/pyhum_test/testclipmeta.mat'
cs2cs_args="epsg:26949"
metadat = loadmat(my_meta_file)
try:
trans = pyproj.Proj(init=cs2cs_args)
except:
trans = pyproj.Proj(cs2cs_args.lstrip(), inverse=True)
lon, lat = trans(np.squeeze(metadat['e']), np.squeeze(metadat['n']), inverse=True)
kml = simplekml.Kml()
ls = kml.newlinestring(name='trackline')
ls.coords = zip(lon,lat)
ls.extrude = 1
ls.altitudemode = simplekml.AltitudeMode.relativetoground
ls.style.linestyle.width = 5
ls.style.linestyle.color = simplekml.Color.red
kml.save('my_new_trackline.kml')
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment