def make_KML_average_columns(averages_dict): kml = KML() aves = {} for point, bucket in averages_dict.items(): aves[point] = {"vel": [], "azm": []} vel = Mean() + point.velocity azm = Mean() + point.azimuth for en in ens: vel += en.velocity azm += en.azimuth aves[point]["vel"] = vel aves[point]["azm"] = azm f = open("../data/column_averages.xml", 'w') f.write( "latitude, longitude, velocity_mean, v_std, v_n, azimuth_mean, a_std, a_n\n" ) mkstr = lambda x: ("%s" % x).replace(", ", "<br/>").replace("=", ": ") for point, ave in aves.items(): vel = ave["vel"] azm = ave["azm"] text = "<h2>Velocity:</h2><p>%s</p>" % mkstr(vel) text += "<h2>Azimuth:</h2><p>%s</p>" % mkstr(azm) desc = description_from_text(text) pt = Placemark(int(azm), point.point, desc) kml.addPlacemark(pt) line = list(point.point) line.append(vel()) line.append(azm()) f.write(",".join(line) + "\n") f.close() kml.output("../data/column_averages.kml")
def make_KML_average_columns(averages_dict): kml = KML() aves = {} for point, bucket in averages_dict.items(): aves[point] = {"vel": [], "azm": []} vel = Mean() + point.velocity azm = Mean() + point.azimuth for en in ens: vel += en.velocity azm += en.azimuth aves[point]["vel"] = vel aves[point]["azm"] = azm f = open("../data/column_averages.xml",'w') f.write("latitude, longitude, velocity_mean, v_std, v_n, azimuth_mean, a_std, a_n\n") mkstr = lambda x: ("%s" % x).replace(", ","<br/>").replace("=",": ") for point, ave in aves.items(): vel = ave["vel"] azm = ave["azm"] text = "<h2>Velocity:</h2><p>%s</p>" % mkstr(vel) text += "<h2>Azimuth:</h2><p>%s</p>" % mkstr(azm) desc = description_from_text(text) pt = Placemark(int(azm), point.point, desc) kml.addPlacemark(pt) line = list(point.point) line.append(vel()) line.append(azm()) f.write(",".join(line) + "\n") f.close() kml.output("../data/column_averages.kml")
def make_KML_at_depths(averages_dict, depths=(3,6,9)): kmls = {} for depth in depths: kmls[depth] = KML() for point, aves in averages_dict.items(): vel = aves[depth]["vel"] azimuth = aves[depth]["azm"] desc = description_from_text("Velocity: %0.3f<br/>Azimuth: %0.3f" % (velocity, azimuth)) pt = Placemark(int(azimuth), point.point, desc) kmls[depth].addPlacemark(pt) for key in kmls.keys(): kmls[key].output("../data/depth%i.kml" % key)
def make_KML_at_depths(averages_dict, depths=(3, 6, 9)): kmls = {} for depth in depths: kmls[depth] = KML() for point, aves in averages_dict.items(): vel = aves[depth]["vel"] azimuth = aves[depth]["azm"] desc = description_from_text("Velocity: %0.3f<br/>Azimuth: %0.3f" % (velocity, azimuth)) pt = Placemark(int(azimuth), point.point, desc) kmls[depth].addPlacemark(pt) for key in kmls.keys(): kmls[key].output("../data/depth%i.kml" % key)