#cursor.execute(query, place_data)
except Exception as e:
print 'unable to insert %s: %s' % (row.ptol_id, e.message)
#connection.commit()
#cursor.close()
#connection.close()
session.commit()
session.close()
except Exception as e:
print 'unable to connect: %s' % (e.message, )
if __name__ == '__main__':
parser = argparse.ArgumentParser(
description='Load Ptolemy places to database.')
parser.add_argument('--sgdb', help='read from sgdb with given prefix')
parser.add_argument('--xlsx', help='xlsx to read from instead of sgdb')
args = parser.parse_args()
if args.sgdb:
places = common.read_places(args.sgdb)
elif args.xlsx:
places = common.read_places_xlsx(args.xlsx)
else:
sys.stderr.write('must specify one of --sgdb or --xlsx')
exit(1)
places = common.read_places_output_csv(args.csv)
#main(places)
kml.write('<Document>\n')
common.write_styles(kml)
for i in range(len(known)):
kml.write(' <Folder id="%s">\n' % (known.ix[i].ptol_id, ))
kml.write(' <name>%s</name>\n' % (known.ix[i].ptol_id, ))
ax = known.ix[i, XCOLS].values
ay = known.ix[i, YCOLS].values
alabel = known.ix[i].ptol_id
common.write_point(kml, 'red', ax, alabel)
common.write_point(kml, 'red', ay, alabel)
if verbose: print known.ix[i].ptol_id, ax, ay
points = [indices[i, j] for j in range(1, k)]
for m in range(len(points)):
j = points[m]
bx = known.ix[j, XCOLS].values
by = known.ix[j, YCOLS].values
blabel = known.ix[j].ptol_id
common.write_point(kml, 'yellow', bx, blabel)
common.write_point(kml, 'yellow', by, blabel)
common.write_line(kml, ax, bx, colors[m])
common.write_line(kml, ay, by, colors[m])
if verbose: print ' ', known.ix[j].ptol_id, bx, by
kml.write(' </Folder>\n')
kml.write('</Document>\n')
kml.write('</kml>\n')
if __name__ == '__main__':
known, _ = common.split_places(common.read_places())
write_neighbor_visualization_kml(known)
"""Compute unknown modern coordinates from known ones using the
given triangulation."""
colors = ['red', 'yellow', 'green']
with open('../Data/visualize_triangles.kml', 'wb') as kml:
kml.write('<?xml version="1.0" encoding="UTF-8"?>\n')
kml.write('<kml xmlns="http://www.opengis.net/kml/2.2">\n')
kml.write('<Document>\n')
common.write_styles(kml)
for simp in tri.simplices:
print simp
tri_name = 'Triangle_%s' % '_'.join(str(x) for x in simp)
kml.write(' <Folder id="%s">\n' % (tri_name, ))
kml.write(' <name>%s</name>\n' % (tri_name, ))
write_triangle(kml, tri_name, colors,
[(known.ix[x].ptol_lat, known.ix[x].ptol_lon)
for x in simp])
write_triangle(kml, tri_name, colors,
[(known.ix[x].modern_lat, known.ix[x].modern_lon)
for x in simp])
kml.write(' </Folder>\n')
kml.write('</Document>\n')
kml.write('</kml>\n')
if __name__ == '__main__':
places = common.read_places()
known, unknown = common.split_places(places)
points = known.loc[:, ['ptol_lat', 'ptol_lon']]
tri = Delaunay(points, furthest_site=False)
write_triangle_visualization_kml(tri, known, unknown)
def main(filename, model):
places = common.read_places()
known, unknown = common.split_places(places)
validate_each(known, model)
compute_errors(known)
known.to_csv(filename, encoding='cp1252')
except Exception as e:
print 'unable to insert %s: %s' % (row.ptol_id, e.message)
#connection.commit()
#cursor.close()
#connection.close()
session.commit()
session.close()
except Exception as e:
print 'unable to connect: %s' % (e.message, )
if __name__ == '__main__':
parser = argparse.ArgumentParser(
description='Load Ptolemy places to database.')
parser.add_argument('--sgdb', help='read from sgdb with given prefix')
parser.add_argument('--xlsx', help='xlsx to read from instead of sgdb')
args = parser.parse_args()
if args.sgdb:
places = common.read_places(args.sgdb)
elif args.xlsx:
places = common.read_places_xlsx(args.xlsx)
else:
sys.stderr.write('must specify one of --sgdb or --xlsx')
exit(1)
places = common.read_places_output_csv(args.csv)
#main(places)
common.write_line(kml, points[i], points[(i+1) % 3], colors[i])
def write_triangle_visualization_kml(tri, known, unknown):
"""Compute unknown modern coordinates from known ones using the
given triangulation."""
colors = ['red', 'yellow', 'green']
with open('../Data/visualize_triangles.kml', 'wb') as kml:
kml.write('<?xml version="1.0" encoding="UTF-8"?>\n')
kml.write('<kml xmlns="http://www.opengis.net/kml/2.2">\n')
kml.write('<Document>\n')
common.write_styles(kml)
for simp in tri.simplices:
print simp
tri_name = 'Triangle_%s' % '_'.join(str(x) for x in simp)
kml.write(' <Folder id="%s">\n' % (tri_name, ))
kml.write(' <name>%s</name>\n' % (tri_name, ))
write_triangle(kml, tri_name, colors, [(known.ix[x].ptol_lat, known.ix[x].ptol_lon) for x in simp])
write_triangle(kml, tri_name, colors, [(known.ix[x].modern_lat, known.ix[x].modern_lon) for x in simp])
kml.write(' </Folder>\n')
kml.write('</Document>\n')
kml.write('</kml>\n')
if __name__ == '__main__':
places = common.read_places()
known, unknown = common.split_places(places)
points = known.loc[:, ['ptol_lat','ptol_lon']]
tri = Delaunay(points, furthest_site=False)
write_triangle_visualization_kml(tri, known, unknown)
for i in range(len(known)):
kml.write(' <Folder id="%s">\n' % (known.ix[i].ptol_id,))
kml.write(" <name>%s</name>\n" % (known.ix[i].ptol_id,))
ax = known.ix[i, XCOLS].values
ay = known.ix[i, YCOLS].values
alabel = known.ix[i].ptol_id
common.write_point(kml, "red", ax, alabel)
common.write_point(kml, "red", ay, alabel)
if verbose:
print known.ix[i].ptol_id, ax, ay
points = [indices[i, j] for j in range(1, k)]
for m in range(len(points)):
j = points[m]
bx = known.ix[j, XCOLS].values
by = known.ix[j, YCOLS].values
blabel = known.ix[j].ptol_id
common.write_point(kml, "yellow", bx, blabel)
common.write_point(kml, "yellow", by, blabel)
common.write_line(kml, ax, bx, colors[m])
common.write_line(kml, ay, by, colors[m])
if verbose:
print " ", known.ix[j].ptol_id, bx, by
kml.write(" </Folder>\n")
kml.write("</Document>\n")
kml.write("</kml>\n")
if __name__ == "__main__":
known, _ = common.split_places(common.read_places())
write_neighbor_visualization_kml(known)
def main(filename, model):
places = common.read_places()
known, unknown = common.split_places(places)
validate_each(known, model)
compute_errors(known)
known.to_csv(filename, encoding='cp1252')
