def main():
#sort out the command line stuff
parser = argparse.ArgumentParser(description=__progdesc__)
arga = get_args(parser)
results = parser.parse_args()
if arga.inosm == '':
print '{0} version {1}'.format(__progdesc__, __version__)
print 'OSM file name is required (-i)'
sys.exit(4)
if not os.path.exists(arga.inosm):
print '{0} version {1}'.format(__progdesc__, __version__)
print 'OSM file {0} cannot be found'.format(arga.inosm)
sys.exit(5)
if arga.verbose:
arga.quiet = False
if not arga.quiet:
print '{0} version {1}'.format(__progdesc__, __version__)
if arga.verbose:
print 'Input OSM file: {0}'.format(arga.inosm)
print 'Text ouput file: {0}'.format(arga.outtxt)
print 'Quiet mode: {0}'.format(arga.quiet)
print 'Verbose: {0}'.format(arga.verbose)
# parse the osm file
OSM = OSMHandler()
saxparser = make_parser()
saxparser.setContentHandler(OSM)
datasource = open(arga.inosm, "r")
saxparser.parse(datasource)
datasource.close
print 'Nodes {0}'.format(len(OSM.Nodes))
print 'Ways {0}'.format(len(OSM.Ways))
print 'Relations {0}'.format(len(OSM.Relations))
# open the database ready for updates
h, u, p, d = hupd(
) # get the host, username password and database from an import file so it's not hardcoded here
conn = mdb.connect(h, u, p, d)
croad = conn.cursor()
# this is a set of hoghway types I'm interested in
amenities = ('hospital', 'doctors', 'dentist')
# find the ways that are highways
for wid in OSM.Ways.keys():
way = OSM.Ways[wid]
if 'amenity' in way.Tags and way.Tags['amenity'] in amenities:
# get north, south east & west max
n = int(way.Nds[0])
north = OSM.Nodes[n].Lat
south = OSM.Nodes[n].Lat
east = OSM.Nodes[n].Lon
west = OSM.Nodes[n].Lon
for n in way.Nds:
if north < OSM.Nodes[int(n)].Lat:
north = OSM.Nodes[int(n)].Lat
if south > OSM.Nodes[int(n)].Lat:
south = OSM.Nodes[int(n)].Lat
if east < OSM.Nodes[int(n)].Lon:
east = OSM.Nodes[int(n)].Lon
if west > OSM.Nodes[int(n)].Lon:
west = OSM.Nodes[int(n)].Lon
croad.execute("INSERT INTO medicalfacility (osmid,north,south,east,west) VALUES(%s,%s,%s,%s,%s)",\
(way.WayID,north,south,east,west))
wayid = croad.lastrowid
for n in way.Nds:
lon = OSM.Nodes[int(n)].Lon
lat = OSM.Nodes[int(n)].Lat
croad.execute(
"INSERT INTO mpoints (medicalfacilityid,lon,lat) VALUES(%s,%s,%s)",
(wayid, lon, lat))
for tag in way.Tags:
croad.execute(
"INSERT INTO tags (medicalfacilityid,okey,oval) VALUES(%s,%s,%s)",
(wayid, tag, way.Tags[tag]))
for nid in OSM.Nodes.keys():
node = OSM.Nodes[nid]
if 'amenity' in node.Tags and node.Tags['amenity'] in amenities:
lat = node.Lat
lon = node.Lon
amenity = node.Tags['amenity']
if 'name' in node.Tags:
name = node.Tags['name']
else:
name = ''
croad.execute("INSERT INTO medicalfacility (osmid,north,south,east,west) VALUES(%s,%s,%s,%s,%s)",\
(node.NodeID,lat,lat,lon,lon))
mfid = croad.lastrowid
croad.execute("INSERT INTO mpoints (medicalfacilityid,lon,lat) VALUES(%s,%s,%s)",\
(mfid,lat,lon))
for tag in node.Tags:
croad.execute(
"INSERT INTO tags (medicalfacilityid,okey,oval) VALUES(%s,%s,%s)",
(node.NodeID, tag, node.Tags[tag]))
conn.close()
def main():
# sort out the command line stuff
parser = argparse.ArgumentParser(description=__progdesc__)
arga = get_args(parser)
results = parser.parse_args()
if arga.inosm == "":
print "{0} version {1}".format(__progdesc__, __version__)
print "OSM file name is required (-i)"
sys.exit(4)
if not os.path.exists(arga.inosm):
print "{0} version {1}".format(__progdesc__, __version__)
print "OSM file {0} cannot be found".format(arga.inosm)
sys.exit(5)
if arga.verbose:
arga.quiet = False
if not arga.quiet:
print "{0} version {1}".format(__progdesc__, __version__)
if arga.verbose:
print "Input OSM file: {0}".format(arga.inosm)
print "Text ouput file: {0}".format(arga.outtxt)
print "Quiet mode: {0}".format(arga.quiet)
print "Verbose: {0}".format(arga.verbose)
# parse the osm file
OSM = OSMHandler()
saxparser = make_parser()
saxparser.setContentHandler(OSM)
datasource = open(arga.inosm, "r")
saxparser.parse(datasource)
datasource.close
print "Nodes {0}".format(len(OSM.Nodes))
print "Ways {0}".format(len(OSM.Ways))
print "Relations {0}".format(len(OSM.Relations))
# open the database ready for updates
h, u, p, d = hupd() # get the host, username password and database from an import file so it's not hardcoded here
conn = mdb.connect(h, u, p, d)
croad = conn.cursor()
# this is a set of hoghway types I'm interested in
hwtype = (
"motorway",
"motorway_link",
"trunk",
"trunk_link",
"primary",
"primary_link",
"secondary",
"tertiary",
"residential",
"unclassified",
)
# find the ways that are highways
for wid in OSM.Ways.keys():
way = OSM.Ways[wid]
if "highway" in way.Tags and way.Tags["highway"] in hwtype:
# get north, south east & west max
n = int(way.Nds[0])
north = OSM.Nodes[n].Lat
south = OSM.Nodes[n].Lat
east = OSM.Nodes[n].Lon
west = OSM.Nodes[n].Lon
for n in way.Nds:
if north < OSM.Nodes[int(n)].Lat:
north = OSM.Nodes[int(n)].Lat
if south > OSM.Nodes[int(n)].Lat:
south = OSM.Nodes[int(n)].Lat
if east < OSM.Nodes[int(n)].Lon:
east = OSM.Nodes[int(n)].Lon
if west > OSM.Nodes[int(n)].Lon:
west = OSM.Nodes[int(n)].Lon
if "maxspeed" in way.Tags:
ms = way.Tags["maxspeed"]
else:
ms = ""
hw = way.Tags["highway"]
croad.execute(
"INSERT INTO road (osmid,north,south,east,west,highway,maxspeed) VALUES(%s,%s,%s,%s,%s,%s,%s)",
(way.WayID, north, south, east, west, hw, ms),
)
roadid = croad.lastrowid
for n in way.Nds:
lon = OSM.Nodes[int(n)].Lon
lat = OSM.Nodes[int(n)].Lat
croad.execute("INSERT INTO roadpoints (roadid,lon,lat) VALUES(%s,%s,%s)", (roadid, lon, lat))
conn.close()
def main():
#sort out the command line stuff
parser = argparse.ArgumentParser(description=__progdesc__)
arga=get_args(parser)
results = parser.parse_args()
if arga.inosm=='':
print '{0} version {1}'.format(__progdesc__,__version__)
print 'OSM file name is required (-i)'
sys.exit(4)
if not os.path.exists(arga.inosm):
print '{0} version {1}'.format(__progdesc__,__version__)
print 'OSM file {0} cannot be found'.format(arga.inosm)
sys.exit(5)
if arga.verbose:
arga.quiet=False
if not arga.quiet:
print '{0} version {1}'.format(__progdesc__,__version__)
if arga.verbose:
print 'Input OSM file: {0}'.format(arga.inosm)
print 'Text ouput file: {0}'.format(arga.outtxt)
print 'Quiet mode: {0}'.format(arga.quiet)
print 'Verbose: {0}'.format(arga.verbose)
# parse the osm file
OSM = OSMHandler()
saxparser = make_parser()
saxparser.setContentHandler(OSM)
datasource = open(arga.inosm,"r")
saxparser.parse(datasource)
datasource.close
print 'Nodes {0}'.format(len(OSM.Nodes))
print 'Ways {0}'.format(len(OSM.Ways))
print 'Relations {0}'.format(len(OSM.Relations))
# open the database ready for updates
h,u,p,d = hupd() # get the host, username password and database from an import file so it's not hardcoded here
conn = mdb.connect(h,u,p,d)
croad = conn.cursor()
# this is a set of hoghway types I'm interested in
amenities=('hospital','doctors','dentist')
# find the ways that are highways
for wid in OSM.Ways.keys():
way = OSM.Ways[wid]
if 'amenity' in way.Tags and way.Tags['amenity'] in amenities:
# get north, south east & west max
n=int(way.Nds[0])
north = OSM.Nodes[n].Lat
south = OSM.Nodes[n].Lat
east = OSM.Nodes[n].Lon
west = OSM.Nodes[n].Lon
for n in way.Nds:
if north < OSM.Nodes[int(n)].Lat:
north = OSM.Nodes[int(n)].Lat
if south > OSM.Nodes[int(n)].Lat:
south = OSM.Nodes[int(n)].Lat
if east < OSM.Nodes[int(n)].Lon:
east = OSM.Nodes[int(n)].Lon
if west > OSM.Nodes[int(n)].Lon:
west = OSM.Nodes[int(n)].Lon
croad.execute("INSERT INTO medicalfacility (osmid,north,south,east,west) VALUES(%s,%s,%s,%s,%s)",\
(way.WayID,north,south,east,west))
wayid=croad.lastrowid
for n in way.Nds:
lon=OSM.Nodes[int(n)].Lon
lat=OSM.Nodes[int(n)].Lat
croad.execute("INSERT INTO mpoints (medicalfacilityid,lon,lat) VALUES(%s,%s,%s)",(wayid,lon,lat))
for tag in way.Tags:
croad.execute("INSERT INTO tags (medicalfacilityid,okey,oval) VALUES(%s,%s,%s)",(wayid,tag,way.Tags[tag]))
for nid in OSM.Nodes.keys():
node = OSM.Nodes[nid]
if 'amenity' in node.Tags and node.Tags['amenity'] in amenities:
lat = node.Lat
lon = node.Lon
amenity = node.Tags['amenity']
if 'name' in node.Tags:
name = node.Tags['name']
else:
name=''
croad.execute("INSERT INTO medicalfacility (osmid,north,south,east,west) VALUES(%s,%s,%s,%s,%s)",\
(node.NodeID,lat,lat,lon,lon))
mfid = croad.lastrowid
croad.execute("INSERT INTO mpoints (medicalfacilityid,lon,lat) VALUES(%s,%s,%s)",\
(mfid,lat,lon))
for tag in node.Tags:
croad.execute("INSERT INTO tags (medicalfacilityid,okey,oval) VALUES(%s,%s,%s)",(node.NodeID,tag,node.Tags[tag]))
conn.close()