'02013' : 'Aleutians East, AK' ,
'02016' : 'Aleutians West, AK' ,
'02020' : 'Anchorage, AK' ,
-'02230' : 'Skagway Municipality' ,
'02050' : 'Bethel, AK' ,
'02060' : 'Bristol Bay, AK' ,
'02068' : 'Denali, AK' ,
'02070' : 'Dillingham, AK' ,
'02090' : 'Fairbanks North Star, AK' ,
'02100' : 'Haines, AK' ,
-'02105' : 'Hoonah-Angoon Census Area' ,
+'02105' : 'Hoonah-Angoon Census Area, AK' ,
'02110' : 'Juneau, AK' ,
'02122' : 'Kenai Peninsula, AK' ,
'02130' : 'Ketchikan Gateway, AK' ,
'02150' : 'Kodiak Island, AK' ,
+'02158' : 'Kusilvak Census Area, AK' ,
'02164' : 'Lake and Peninsula, AK' ,
'02170' : 'Matanuska-Susitna, AK' ,
'02180' : 'Nome, AK' ,
'02185' : 'North Slope, AK' ,
'02188' : 'Northwest Arctic, AK' ,
-'02195' : 'Petersburg Census Area' ,
-'02198' : 'Prince of Wales-Hyder Census Area' ,
+'02195' : 'Petersburg, AK' ,
+'02198' : 'Prince of Wales-Hyder Census Area, AK' ,
'02201' : 'Prince of Wales-Outer Ketchikan, AK' ,
'02220' : 'Sitka, AK' ,
+'02230' : 'Skagway, AK' ,
'02232' : 'Skagway-Hoonah-Angoon, AK' ,
'02240' : 'Southeast Fairbanks, AK' ,
'02261' : 'Valdez-Cordova, AK' ,
'02270' : 'Wade Hampton, AK' ,
-'02275' : 'Wrangell City and Borough' ,
+'02275' : 'Wrangell City and Borough, AK' ,
'02280' : 'Wrangell-Petersburg, AK' ,
'02282' : 'Yakutat, AK' ,
'02290' : 'Yukon-Koyukuk, AK' ,
'17093' : 'Kendall, IL' ,
'17095' : 'Knox, IL' ,
'17097' : 'Lake, IL' ,
-'17099' : 'La Salle, IL' ,
+'17099' : 'LaSalle, IL' ,
'17101' : 'Lawrence, IL' ,
'17103' : 'Lee, IL' ,
'17105' : 'Livingston, IL' ,
'22053' : 'Jefferson Davis, LA' ,
'22055' : 'Lafayette, LA' ,
'22057' : 'Lafourche, LA' ,
-'22059' : 'La Salle, LA' ,
+'22059' : 'LaSalle, LA' ,
'22061' : 'Lincoln, LA' ,
'22063' : 'Livingston, LA' ,
'22065' : 'Madison, LA' ,
'24043' : 'Washington, MD' ,
'24045' : 'Wicomico, MD' ,
'24047' : 'Worcester, MD' ,
-'24510' : 'Baltimore city, MD' ,
+'24510' : 'Baltimore City, MD' ,
'25000' : 'MASSACHUSETTS' ,
'25001' : 'Barnstable, MA' ,
'25003' : 'Berkshire, MA' ,
'29225' : 'Webster, MO' ,
'29227' : 'Worth, MO' ,
'29229' : 'Wright, MO' ,
-'29510' : 'St. Louis city, MO' ,
+'29510' : 'St. Louis City, MO' ,
'30000' : 'MONTANA' ,
'30001' : 'Beaverhead, MT' ,
'30003' : 'Big Horn, MT' ,
'32029' : 'Storey, NV' ,
'32031' : 'Washoe, NV' ,
'32033' : 'White Pine, NV' ,
-'32510' : 'Carson City city, NV' ,
+'32510' : 'Carson City, NV' ,
'33000' : 'NEW HAMPSHIRE' ,
'33001' : 'Belknap, NH' ,
'33003' : 'Carroll, NH' ,
'46097' : 'Miner, SD' ,
'46099' : 'Minnehaha, SD' ,
'46101' : 'Moody, SD' ,
+'46102' : 'Oglala Lakota, SD' ,
'46103' : 'Pennington, SD' ,
'46105' : 'Perkins, SD' ,
'46107' : 'Potter, SD' ,
'56043' : 'Washakie, WY' ,
'56045' : 'Weston, WY' ,
'60000' : 'AMERICAN SAMOA',
-'60010' : 'Eastern district, AS' ,
-'60020' : 'Manua district, AS' ,
-'60030' : 'Rose island, AS' ,
-'60040' : 'Swains island, AS' ,
-'60050' : 'Western district, AS' ,
+'60010' : 'Eastern District, AS' ,
+'60020' : 'Manu\'a District, AS' ,
+'60030' : 'Rose Island, AS' ,
+'60040' : 'Swains Island, AS' ,
+'60050' : 'Western District, AS' ,
'66000' : 'GUAM' ,
'66010' : 'Guam, GU' ,
'69000' : 'COMMONWEALTH OF THE NORTHERN MARIANA ISLANDS' ,
-'69085' : 'Commonwealth of the Northern Mariana Islands, MP' ,
-'69100' : 'Commonwealth of the Northern Mariana Islands, MP' ,
-'69110' : 'Commonwealth of the Northern Mariana Islands, MP' ,
-'69120' : 'Commonwealth of the Northern Mariana Islands, MP' ,
-'72000' : 'PUERTO RICO',
-'72001' : 'Adjuntas, PR',
-'72003' : 'Aguada, PR',
-'72005' : 'Aguadilla, PR',
-'72007' : 'Aguas Buenas, PR',
-'72009' : 'Aibonito, PR',
-'72011' : 'Anasco, PR',
-'72013' : 'Arecibo, PR',
-'72015' : 'Arroyo, PR',
-'72017' : 'Barceloneta, PR',
-'72019' : 'Barranquitas, PR',
-'72021' : 'Bayamon, PR',
-'72023' : 'Cabo Rojo, PR',
-'72025' : 'Caguas, PR',
-'72027' : 'Camuy, PR',
-'72029' : 'Canovanas, PR',
-'72031' : 'Carolina, PR',
-'72033' : 'Catano, PR',
-'72035' : 'Cayey, PR',
-'72037' : 'Ceiba, PR',
-'72039' : 'Ciales, PR',
-'72041' : 'Cidra, PR',
-'72043' : 'Coamo, PR',
-'72045' : 'Comerio, PR',
-'72047' : 'Corozal, PR',
-'72049' : 'Culebra, PR',
-'72051' : 'Dorado, PR',
-'72053' : 'Fajardo, PR',
-'72054' : 'Florida, PR',
-'72055' : 'Guanica, PR',
-'72057' : 'Guayama, PR',
-'72059' : 'Guayanilla, PR',
-'72061' : 'Guaynabo, PR',
-'72063' : 'Gurabo, PR',
-'72065' : 'Hatillo, PR',
-'72067' : 'Hormigueros, PR',
-'72069' : 'Humacao, PR',
-'72071' : 'Isabela, PR',
-'72073' : 'Jayuya, PR',
-'72075' : 'Juana Diaz, PR',
-'72077' : 'Juncos, PR',
-'72079' : 'Lajas, PR',
-'72081' : 'Lares, PR',
-'72083' : 'Las Marias, PR',
-'72085' : 'Las Piedras, PR',
-'72087' : 'Loiza, PR',
-'72089' : 'Luquillo, PR',
-'72091' : 'Manati, PR',
-'72093' : 'Maricao, PR',
-'72095' : 'Maunabo, PR',
-'72097' : 'Mayaguez, PR',
-'72099' : 'Moca, PR',
-'72101' : 'Morovis, PR',
-'72103' : 'Naguabo, PR',
-'72105' : 'Naranjito, PR',
-'72107' : 'Orocovis, PR',
-'72109' : 'Patillas, PR',
-'72111' : 'Penuelas, PR',
-'72113' : 'Ponce, PR',
-'72115' : 'Quebradillas, PR',
-'72117' : 'Rincon, PR',
-'72119' : 'Rio Grande, PR',
-'72121' : 'Sabana Grande, PR',
-'72123' : 'Salinas, PR',
-'72125' : 'San German, PR',
-'72127' : 'San Juan, PR',
-'72129' : 'San Lorenzo, PR',
-'72131' : 'San Sebastian, PR',
-'72133' : 'Santa Isabel, PR',
-'72135' : 'Toa Alta, PR',
-'72137' : 'Toa Baja, PR',
-'72139' : 'Trujillo Alto, PR',
-'72141' : 'Utuado, PR',
-'72143' : 'Vega Alta, PR',
-'72145' : 'Vega Baja, PR',
-'72147' : 'Vieques, PR',
-'72149' : 'Villalba, PR',
-'72151' : 'Yabucoa, PR',
-'72153' : 'Yauco, PR',
+'69085' : 'Northern Islands, MP' ,
+'69100' : 'Rota, MP' ,
+'69110' : 'Saipan, MP' ,
+'69120' : 'Tinian, MP' ,
+'72000' : 'PUERTO RICO' ,
+'72001' : 'Adjuntas, PR' ,
+'72003' : 'Aguada, PR' ,
+'72005' : 'Aguadilla, PR' ,
+'72007' : 'Aguas Buenas, PR' ,
+'72009' : 'Aibonito, PR' ,
+'72011' : 'Anasco, PR' ,
+'72013' : 'Arecibo, PR' ,
+'72015' : 'Arroyo, PR' ,
+'72017' : 'Barceloneta, PR' ,
+'72019' : 'Barranquitas, PR' ,
+'72021' : 'Bayamon, PR' ,
+'72023' : 'Cabo Rojo, PR' ,
+'72025' : 'Caguas, PR' ,
+'72027' : 'Camuy, PR' ,
+'72029' : 'Canovanas, PR' ,
+'72031' : 'Carolina, PR' ,
+'72033' : 'Catano, PR' ,
+'72035' : 'Cayey, PR' ,
+'72037' : 'Ceiba, PR' ,
+'72039' : 'Ciales, PR' ,
+'72041' : 'Cidra, PR' ,
+'72043' : 'Coamo, PR' ,
+'72045' : 'Comerio, PR' ,
+'72047' : 'Corozal, PR' ,
+'72049' : 'Culebra, PR' ,
+'72051' : 'Dorado, PR' ,
+'72053' : 'Fajardo, PR' ,
+'72054' : 'Florida, PR' ,
+'72055' : 'Guanica, PR' ,
+'72057' : 'Guayama, PR' ,
+'72059' : 'Guayanilla, PR' ,
+'72061' : 'Guaynabo, PR' ,
+'72063' : 'Gurabo, PR' ,
+'72065' : 'Hatillo, PR' ,
+'72067' : 'Hormigueros, PR' ,
+'72069' : 'Humacao, PR' ,
+'72071' : 'Isabela, PR' ,
+'72073' : 'Jayuya, PR' ,
+'72075' : 'Juana Diaz, PR' ,
+'72077' : 'Juncos, PR' ,
+'72079' : 'Lajas, PR' ,
+'72081' : 'Lares, PR' ,
+'72083' : 'Las Marias, PR' ,
+'72085' : 'Las Piedras, PR' ,
+'72087' : 'Loiza, PR' ,
+'72089' : 'Luquillo, PR' ,
+'72091' : 'Manati, PR' ,
+'72093' : 'Maricao, PR' ,
+'72095' : 'Maunabo, PR' ,
+'72097' : 'Mayaguez, PR' ,
+'72099' : 'Moca, PR' ,
+'72101' : 'Morovis, PR' ,
+'72103' : 'Naguabo, PR' ,
+'72105' : 'Naranjito, PR' ,
+'72107' : 'Orocovis, PR' ,
+'72109' : 'Patillas, PR' ,
+'72111' : 'Penuelas, PR' ,
+'72113' : 'Ponce, PR' ,
+'72115' : 'Quebradillas, PR' ,
+'72117' : 'Rincon, PR' ,
+'72119' : 'Rio Grande, PR' ,
+'72121' : 'Sabana Grande, PR' ,
+'72123' : 'Salinas, PR' ,
+'72125' : 'San German, PR' ,
+'72127' : 'San Juan, PR' ,
+'72129' : 'San Lorenzo, PR' ,
+'72131' : 'San Sebastian, PR' ,
+'72133' : 'Santa Isabel, PR' ,
+'72135' : 'Toa Alta, PR' ,
+'72137' : 'Toa Baja, PR' ,
+'72139' : 'Trujillo Alto, PR' ,
+'72141' : 'Utuado, PR' ,
+'72143' : 'Vega Alta, PR' ,
+'72145' : 'Vega Baja, PR' ,
+'72147' : 'Vieques, PR' ,
+'72149' : 'Villalba, PR' ,
+'72151' : 'Yabucoa, PR' ,
+'72153' : 'Yauco, PR' ,
'78000' : 'VIRGIN ISLANDS' ,
'78010' : 'St. Croix, VI' ,
'78020' : 'St. John, VI' ,
return tags
if fips not in fipscodes:
- raise KeyError, 'missing FIPS code', fips
+ raise KeyError('missing FIPS code', fips)
state, statecode, isocode = fipscodes[fips]
county_fips_code = fips + "" + countyfp
poLayer = poDS.GetLayer( 0 )
+ fieldNameList = []
+ layerDefinition = poLayer.GetLayerDefn()
+ for i in range(layerDefinition.GetFieldCount()):
+ fieldNameList.append(layerDefinition.GetFieldDefn(i).GetName())
+ # sys.stderr.write(",".join(fieldNameList))
+
poLayer.ResetReading()
ret = []
# WAY ID
tags[iSource + ":way_id"] = int( poFeature.GetField("TLID") )
- # FEATURE IDENTIFICATION
+ # FEATURE IDENTIFICATION
mtfcc = poFeature.GetField("MTFCC");
if mtfcc != None:
- if mtfcc == "L4010": #Pipeline
- tags["man_made"] = "pipeline"
- if mtfcc == "L4020": #Powerline
- tags["power"] = "line"
- if mtfcc == "L4031": #Aerial Tramway/Ski Lift
- tags["aerialway"] = "cable_car"
- if mtfcc == "L4110": #Fence Line
- tags["barrier"] = "fence"
- if mtfcc == "L4125": #Cliff/Escarpment
- tags["natural"] = "cliff"
- if mtfcc == "L4165": #Ferry Crossing
- tags["route"] = "ferry"
- if mtfcc == "R1011": #Railroad Feature (Main, Spur, or Yard)
- tags["railway"] = "rail"
- ttyp = poFeature.GetField("TTYP")
- if ttyp != None:
- if ttyp == "S":
- tags["service"] = "spur"
- if ttyp == "Y":
- tags["service"] = "yard"
- tags["tiger:ttyp"] = ttyp
- if mtfcc == "R1051": #Carline, Streetcar Track, Monorail, Other Mass Transit Rail)
- tags["railway"] = "light_rail"
- if mtfcc == "R1052": #Cog Rail Line, Incline Rail Line, Tram
- tags["railway"] = "incline"
- if mtfcc == "S1100":
- tags["highway"] = "primary"
- if mtfcc == "S1200":
- tags["highway"] = "secondary"
- if mtfcc == "S1400":
- tags["highway"] = "residential"
- if mtfcc == "S1500":
- tags["highway"] = "track"
- if mtfcc == "S1630": #Ramp
- tags["highway"] = "motorway_link"
- if mtfcc == "S1640": #Service Drive usually along a limited access highway
- tags["highway"] = "service"
- if mtfcc == "S1710": #Walkway/Pedestrian Trail
- tags["highway"] = "path"
- if mtfcc == "S1720":
- tags["highway"] = "steps"
- if mtfcc == "S1730": #Alley
- tags["highway"] = "service"
- tags["service"] = "alley"
- if mtfcc == "S1740": #Private Road for service vehicles (logging, oil, fields, ranches, etc.)
- tags["highway"] = "service"
- tags["access"] = "private"
- if mtfcc == "S1750": #Private Driveway
- tags["highway"] = "service"
- tags["access"] = "private"
- tags["service"] = "driveway"
- if mtfcc == "S1780": #Parking Lot Road
- tags["highway"] = "service"
- tags["service"] = "parking_aisle"
- if mtfcc == "S1820": #Bike Path or Trail
- tags["highway"] = "cycleway"
- if mtfcc == "S1830": #Bridle Path
- tags["highway"] = "bridleway"
- tags["tiger:mtfcc"] = mtfcc
+ if mtfcc == "L4010": #Pipeline
+ tags["man_made"] = "pipeline"
+ if mtfcc == "L4020": #Powerline
+ tags["power"] = "line"
+ if mtfcc == "L4031": #Aerial Tramway/Ski Lift
+ tags["aerialway"] = "cable_car"
+ if mtfcc == "L4110": #Fence Line
+ tags["barrier"] = "fence"
+ if mtfcc == "L4125": #Cliff/Escarpment
+ tags["natural"] = "cliff"
+ if mtfcc == "L4165": #Ferry Crossing
+ tags["route"] = "ferry"
+ if mtfcc == "R1011": #Railroad Feature (Main, Spur, or Yard)
+ tags["railway"] = "rail"
+ ttyp = poFeature.GetField("TTYP")
+ if ttyp != None:
+ if ttyp == "S":
+ tags["service"] = "spur"
+ if ttyp == "Y":
+ tags["service"] = "yard"
+ tags["tiger:ttyp"] = ttyp
+ if mtfcc == "R1051": #Carline, Streetcar Track, Monorail, Other Mass Transit Rail)
+ tags["railway"] = "light_rail"
+ if mtfcc == "R1052": #Cog Rail Line, Incline Rail Line, Tram
+ tags["railway"] = "incline"
+ if mtfcc == "S1100":
+ tags["highway"] = "primary"
+ if mtfcc == "S1200":
+ tags["highway"] = "secondary"
+ if mtfcc == "S1400":
+ tags["highway"] = "residential"
+ if mtfcc == "S1500":
+ tags["highway"] = "track"
+ if mtfcc == "S1630": #Ramp
+ tags["highway"] = "motorway_link"
+ if mtfcc == "S1640": #Service Drive usually along a limited access highway
+ tags["highway"] = "service"
+ if mtfcc == "S1710": #Walkway/Pedestrian Trail
+ tags["highway"] = "path"
+ if mtfcc == "S1720":
+ tags["highway"] = "steps"
+ if mtfcc == "S1730": #Alley
+ tags["highway"] = "service"
+ tags["service"] = "alley"
+ if mtfcc == "S1740": #Private Road for service vehicles (logging, oil, fields, ranches, etc.)
+ tags["highway"] = "service"
+ tags["access"] = "private"
+ if mtfcc == "S1750": #Private Driveway
+ tags["highway"] = "service"
+ tags["access"] = "private"
+ tags["service"] = "driveway"
+ if mtfcc == "S1780": #Parking Lot Road
+ tags["highway"] = "service"
+ tags["service"] = "parking_aisle"
+ if mtfcc == "S1820": #Bike Path or Trail
+ tags["highway"] = "cycleway"
+ if mtfcc == "S1830": #Bridle Path
+ tags["highway"] = "bridleway"
+ tags["tiger:mtfcc"] = mtfcc
# FEATURE NAME
if poFeature.GetField("FULLNAME"):
name = poFeature.GetField( "FULLNAME" )
tags["name"] = name
- #Attempt to guess highway grade
- if name[0:2] == "I-":
- tags["highway"] = "motorway"
- if name[0:3] == "US ":
- tags["highway"] = "primary"
- if name[0:3] == "US-":
- tags["highway"] = "primary"
- if name[0:3] == "Hwy":
- if tags["highway"] != "primary":
- tags["highway"] = "secondary"
+ #Attempt to guess highway grade
+ if name[0:2] == "I-":
+ tags["highway"] = "motorway"
+ if name[0:3] == "US ":
+ tags["highway"] = "primary"
+ if name[0:3] == "US-":
+ tags["highway"] = "primary"
+ if name[0:3] == "Hwy":
+ if tags["highway"] != "primary":
+ tags["highway"] = "secondary"
- divroad = poFeature.GetField("DIVROAD")
- if divroad != None:
- if divroad == "Y" and "highway" in tags and tags["highway"] == "residential":
- tags["highway"] = "tertiary"
- tags["tiger:separated"] = divroad
+ # TIGER 2017 no longer contains this field
+ if 'DIVROAD' in fieldNameList:
+ divroad = poFeature.GetField("DIVROAD")
+ if divroad != None:
+ if divroad == "Y" and "highway" in tags and tags["highway"] == "residential":
+ tags["highway"] = "tertiary"
+ tags["tiger:separated"] = divroad
statefp = poFeature.GetField("STATEFP")
countyfp = poFeature.GetField("COUNTYFP")
# ====================================
projcs_wkt = \
"""GEOGCS["GCS_North_American_1983",
- DATUM["D_North_American_1983",
- SPHEROID["GRS_1980",6378137,298.257222101]],
- PRIMEM["Greenwich",0],
- UNIT["Degree",0.017453292519943295]]"""
+ DATUM["D_North_American_1983",
+ SPHEROID["GRS_1980",6378137,298.257222101]],
+ PRIMEM["Greenwich",0],
+ UNIT["Degree",0.017453292519943295]]"""
from_proj = osr.SpatialReference()
from_proj.ImportFromWkt( projcs_wkt )
pointid, (lat, lon) = nodelist[ round_point( point ) ]
if first:
first = False
- else:
- #The approximate number of feet in one degree of longitute
+ else:
+ #The approximate number of feet in one degree of longitute
lrad = math.radians(lat)
lon_feet = 365527.822 * math.cos(lrad) - 306.75853 * math.cos(3 * lrad) + 0.3937 * math.cos(5 * lrad)
- distance += math.sqrt(((lat - previous[0])*lat_feet)**2 + ((lon - previous[1])*lon_feet)**2)
- previous = (lat, lon)
+ distance += math.sqrt(((lat - previous[0])*lat_feet)**2 + ((lon - previous[1])*lon_feet)**2)
+ previous = (lat, lon)
return distance
def addressways(waylist, nodelist, first_id):
distance = float(address_distance)
ret = []
- for waykey, segments in waylist.iteritems():
+ for waykey, segments in waylist.items():
waykey = dict(waykey)
rsegments = []
lsegments = []
rsegment = []
lastpoint = None
- #Don't pull back the ends of very short ways too much
- seglength = length(segment, nodelist)
- if seglength < float(address_pullback) * 3.0:
- pullback = seglength / 3.0
- else:
+ #Don't pull back the ends of very short ways too much
+ seglength = length(segment, nodelist)
+ if seglength < float(address_pullback) * 3.0:
+ pullback = seglength / 3.0
+ else:
pullback = float(address_pullback)
if "tiger:lfromadd" in waykey:
lfromadd = waykey["tiger:lfromadd"]
if "tiger:rtoadd" in waykey:
rtoadd = waykey["tiger:rtoadd"]
else:
- rtoadd = None
+ rtoadd = None
if rfromadd != None and rtoadd != None:
right = True
- else:
- right = False
+ else:
+ right = False
if lfromadd != None and ltoadd != None:
left = True
- else:
- left = False
+ else:
+ left = False
if left or right:
- first = True
+ first = True
firstpointid, firstpoint = nodelist[ round_point( segment[0] ) ]
finalpointid, finalpoint = nodelist[ round_point( segment[len(segment) - 1] ) ]
for point in segment:
pointid, (lat, lon) = nodelist[ round_point( point ) ]
- #The approximate number of feet in one degree of longitute
+ #The approximate number of feet in one degree of longitute
lrad = math.radians(lat)
lon_feet = 365527.822 * math.cos(lrad) - 306.75853 * math.cos(3 * lrad) + 0.3937 * math.cos(5 * lrad)
#Calculate the points of the offset ways
if lastpoint != None:
- #Skip points too close to start
- if math.sqrt((lat * lat_feet - firstpoint[0] * lat_feet)**2 + (lon * lon_feet - firstpoint[1] * lon_feet)**2) < pullback:
- #Preserve very short ways (but will be rendered backwards)
- if pointid != finalpointid:
- continue
- #Skip points too close to end
- if math.sqrt((lat * lat_feet - finalpoint[0] * lat_feet)**2 + (lon * lon_feet - finalpoint[1] * lon_feet)**2) < pullback:
- #Preserve very short ways (but will be rendered backwards)
- if (pointid != firstpointid) and (pointid != finalpointid):
- continue
+ #Skip points too close to start
+ if math.sqrt((lat * lat_feet - firstpoint[0] * lat_feet)**2 + (lon * lon_feet - firstpoint[1] * lon_feet)**2) < pullback:
+ #Preserve very short ways (but will be rendered backwards)
+ if pointid != finalpointid:
+ continue
+ #Skip points too close to end
+ if math.sqrt((lat * lat_feet - finalpoint[0] * lat_feet)**2 + (lon * lon_feet - finalpoint[1] * lon_feet)**2) < pullback:
+ #Preserve very short ways (but will be rendered backwards)
+ if (pointid != firstpointid) and (pointid != finalpointid):
+ continue
X = (lon - lastpoint[1]) * lon_feet
- Y = (lat - lastpoint[0]) * lat_feet
+ Y = (lat - lastpoint[0]) * lat_feet
if Y != 0:
- theta = math.pi/2 - math.atan( X / Y)
- Xp = math.sin(theta) * distance
- Yp = math.cos(theta) * distance
+ theta = math.pi/2 - math.atan( X / Y)
+ Xp = math.sin(theta) * distance
+ Yp = math.cos(theta) * distance
else:
Xp = 0
- if X > 0:
+ if X > 0:
Yp = -distance
- else:
+ else:
Yp = distance
- if Y > 0:
- Xp = -Xp
- else:
- Yp = -Yp
-
- if first:
- first = False
- dX = - (Yp * (pullback / distance)) / lon_feet #Pull back the first point
- dY = (Xp * (pullback / distance)) / lat_feet
- if left:
+ if Y > 0:
+ Xp = -Xp
+ else:
+ Yp = -Yp
+
+ if first:
+ first = False
+ dX = - (Yp * (pullback / distance)) / lon_feet #Pull back the first point
+ dY = (Xp * (pullback / distance)) / lat_feet
+ if left:
lpoint = (lastpoint[0] + (Yp / lat_feet) - dY, lastpoint[1] + (Xp / lon_feet) - dX)
lsegment.append( (id, lpoint) )
- id += 1
- if right:
+ id += 1
+ if right:
rpoint = (lastpoint[0] - (Yp / lat_feet) - dY, lastpoint[1] - (Xp / lon_feet) - dX)
rsegment.append( (id, rpoint) )
- id += 1
+ id += 1
- else:
- #round the curves
- if delta[1] != 0:
- theta = abs(math.atan(delta[0] / delta[1]))
- else:
- theta = math.pi / 2
- if Xp != 0:
- theta = theta - abs(math.atan(Yp / Xp))
- else: theta = theta - math.pi / 2
- r = 1 + abs(math.tan(theta/2))
- if left:
- lpoint = (lastpoint[0] + (Yp + delta[0]) * r / (lat_feet * 2), lastpoint[1] + (Xp + delta[1]) * r / (lon_feet * 2))
+ else:
+ #round the curves
+ if delta[1] != 0:
+ theta = abs(math.atan(delta[0] / delta[1]))
+ else:
+ theta = math.pi / 2
+ if Xp != 0:
+ theta = theta - abs(math.atan(Yp / Xp))
+ else: theta = theta - math.pi / 2
+ r = 1 + abs(math.tan(theta/2))
+ if left:
+ lpoint = (lastpoint[0] + (Yp + delta[0]) * r / (lat_feet * 2), lastpoint[1] + (Xp + delta[1]) * r / (lon_feet * 2))
lsegment.append( (id, lpoint) )
id += 1
- if right:
+ if right:
rpoint = (lastpoint[0] - (Yp + delta[0]) * r / (lat_feet * 2), lastpoint[1] - (Xp + delta[1]) * r / (lon_feet * 2))
-
+
rsegment.append( (id, rpoint) )
id += 1
#Add in the last node
- dX = - (Yp * (pullback / distance)) / lon_feet
- dY = (Xp * (pullback / distance)) / lat_feet
- if left:
+ dX = - (Yp * (pullback / distance)) / lon_feet
+ dY = (Xp * (pullback / distance)) / lat_feet
+ if left:
lpoint = (lastpoint[0] + (Yp + delta[0]) / (lat_feet * 2) + dY, lastpoint[1] + (Xp + delta[1]) / (lon_feet * 2) + dX )
lsegment.append( (id, lpoint) )
id += 1
- if right:
+ if right:
rpoint = (lastpoint[0] - Yp / lat_feet + dY, lastpoint[1] - Xp / lon_feet + dX)
rsegment.append( (id, rpoint) )
id += 1
#Generate the tags for ways and nodes
- rtags = []
- ltags = []
- tags = []
+ rtags = []
+ ltags = []
+ tags = []
zipr = ''
zipl = ''
name = ''
county = ''
- if "tiger:zip_right" in waykey:
- zipr = waykey["tiger:zip_right"]
+ if "tiger:zip_right" in waykey:
+ zipr = waykey["tiger:zip_right"]
rtags.append( "<tag k=\"addr:postcode\" v=\"%s\" />" % zipr )
- if "tiger:zip_left" in waykey:
- zipl = waykey["tiger:zip_left"]
+ if "tiger:zip_left" in waykey:
+ zipl = waykey["tiger:zip_left"]
ltags.append( "<tag k=\"addr:postcode\" v=\"%s\" />" % zipl )
if "name" in waykey:
name = waykey["name"]
tags.append( "<tag k=\"addr:street\" v=\"%s\" />" % name )
- if "is_in:state" in waykey:
- state = waykey["is_in:state"]
+ if "is_in:state" in waykey:
+ state = waykey["is_in:state"]
tags.append( "<tag k=\"addr:state\" v=\"%s\" />" % state )
- if "tiger:county" in waykey:
- county = waykey["tiger:county"]
+ if "tiger:county" in waykey:
+ county = waykey["tiger:county"]
tags.append( "<tag k=\"addr:county\" v=\"%s\" />" % county )
- if "is_in:country_code" in waykey:
- country = waykey["is_in:country_code"]
+ if "is_in:country_code" in waykey:
+ country = waykey["is_in:country_code"]
tags.append( "<tag k=\"addr:country\" v=\"%s\" />" % country )
- if "tiger:separated" in waykey:
- separated = waykey["tiger:separated"]
- else:
- separated = "N"
- ltags.extend(tags)
- rtags.extend(tags)
+ if "tiger:separated" in waykey:
+ separated = waykey["tiger:separated"]
+ else:
+ separated = "N"
+ ltags.extend(tags)
+ rtags.extend(tags)
#Write the nodes of the offset ways
- if right:
+ if right:
rlinestring = [];
for i, point in rsegment:
rlinestring.append( "%f %f" % (point[1], point[0]) )
- if left:
+ if left:
llinestring = [];
for i, point in lsegment:
llinestring.append( "%f %f" % (point[1], point[0]) )
- if right:
+ if right:
rsegments.append( rsegment )
- if left:
+ if left:
lsegments.append( lsegment )
- rtofromint = right #Do the addresses convert to integers?
- ltofromint = left #Do the addresses convert to integers?
- if right:
- try: rfromint = int(rfromadd)
- except:
- print("Non integer address: %s" % rfromadd)
- rtofromint = False
- try: rtoint = int(rtoadd)
- except:
- print("Non integer address: %s" % rtoadd)
- rtofromint = False
- if left:
- try: lfromint = int(lfromadd)
- except:
- print("Non integer address: %s" % lfromadd)
- ltofromint = False
- try: ltoint = int(ltoadd)
- except:
- print("Non integer address: %s" % ltoadd)
- ltofromint = False
- import_guid = time.strftime( '%Y%m%d%H%M%S' )
- if right:
- id += 1
+ rtofromint = right #Do the addresses convert to integers?
+ ltofromint = left #Do the addresses convert to integers?
+ if right:
+ try: rfromint = int(rfromadd)
+ except:
+ print("Non integer address: %s" % rfromadd)
+ rtofromint = False
+ try: rtoint = int(rtoadd)
+ except:
+ print("Non integer address: %s" % rtoadd)
+ rtofromint = False
+ if left:
+ try: lfromint = int(lfromadd)
+ except:
+ print("Non integer address: %s" % lfromadd)
+ ltofromint = False
+ try: ltoint = int(ltoadd)
+ except:
+ print("Non integer address: %s" % ltoadd)
+ ltofromint = False
+ import_guid = time.strftime( '%Y%m%d%H%M%S' )
+ if right:
+ id += 1
interpolationtype = "";
- if rtofromint:
+ if rtofromint:
if (rfromint % 2) == 0 and (rtoint % 2) == 0:
- if separated == "Y": #Doesn't matter if there is another side
+ if separated == "Y": #Doesn't matter if there is another side
# ret.append( "<tag k=\"addr:interpolation\" v=\"even\" />" )
interpolationtype = "even";
- elif ltofromint and (lfromint % 2) == 1 and (ltoint % 2) == 1:
+ elif ltofromint and (lfromint % 2) == 1 and (ltoint % 2) == 1:
interpolationtype = "even";
# ret.append( "<tag k=\"addr:interpolation\" v=\"even\" />" )
- else:
+ else:
interpolationtype = "all";
# ret.append( "<tag k=\"addr:interpolation\" v=\"all\" />" )
elif (rfromint % 2) == 1 and (rtoint % 2) == 1:
- if separated == "Y": #Doesn't matter if there is another side
+ if separated == "Y": #Doesn't matter if there is another side
interpolationtype = "odd";
# ret.append( "<tag k=\"addr:interpolation\" v=\"odd\" />" )
- elif ltofromint and (lfromint % 2) == 0 and (ltoint % 2) == 0:
+ elif ltofromint and (lfromint % 2) == 0 and (ltoint % 2) == 0:
interpolationtype = "odd";
# ret.append( "<tag k=\"addr:interpolation\" v=\"odd\" />" )
- else:
+ else:
interpolationtype = "all";
# ret.append( "<tag k=\"addr:interpolation\" v=\"all\" />" )
- else:
+ else:
interpolationtype = "all";
# ret.append( "<tag k=\"addr:interpolation\" v=\"all\" />" )
- else:
+ else:
interpolationtype = "all";
# ret.append( "<tag k=\"addr:interpolation\" v=\"all\" />" )
-# ret.extend(rtags)
+# ret.extend(rtags)
# ret.append( "<tag k=\"source\" v=\"%s_import_v%s_%s\" />" % (iSource, VERSION, import_guid) )
# ret.append( "<tag k=\"attribution\" v=\"%s\" />" % (iAttrib) )
# ret.append( "</way>" )
-
- ret.append( "select tigger_create_interpolation(ST_GeomFromText('LINESTRING(%s)',4326), '%s', '%s', '%s', '%s', '%s', '%s');" %
+ # call new tiger_line_import function to save the lines in the DB.
+ ret.append( "select tiger_line_import(ST_GeomFromText('LINESTRING(%s)',4326), '%s', '%s', '%s', '%s', '%s', '%s');" %
( ",".join(rlinestring), rfromadd.replace("'", "''"), rtoadd.replace("'", "''"), interpolationtype.replace("'", "''"), name.replace("'", "''"), county.replace("'", "''"), zipr.replace("'", "''") ) )
- if left:
- id += 1
- if ltofromint:
+ if left:
+ id += 1
+ if ltofromint:
if (lfromint % 2) == 0 and (ltoint % 2) == 0:
- if separated == "Y":
+ if separated == "Y":
interpolationtype = "even";
- elif rtofromint and (rfromint % 2) == 1 and (rtoint % 2) == 1:
+ elif rtofromint and (rfromint % 2) == 1 and (rtoint % 2) == 1:
interpolationtype = "even";
- else:
+ else:
interpolationtype = "all";
elif (lfromint % 2) == 1 and (ltoint % 2) == 1:
- if separated == "Y":
+ if separated == "Y":
interpolationtype = "odd";
- elif rtofromint and (rfromint %2 ) == 0 and (rtoint % 2) == 0:
+ elif rtofromint and (rfromint %2 ) == 0 and (rtoint % 2) == 0:
interpolationtype = "odd";
- else:
+ else:
interpolationtype = "all";
- else:
+ else:
interpolationtype = "all";
- else:
+ else:
interpolationtype = "all";
- ret.append( "select tigger_create_interpolation(ST_GeomFromText('LINESTRING(%s)',4326), '%s', '%s', '%s', '%s', '%s', '%s');" %
+ ret.append( "select tiger_line_import(ST_GeomFromText('LINESTRING(%s)',4326), '%s', '%s', '%s', '%s', '%s', '%s');" %
( ",".join(llinestring), lfromadd.replace("'", "''"), ltoadd.replace("'", "''"), interpolationtype.replace("'", "''"), name.replace("'", "''"), county.replace("'", "''"), zipl.replace("'", "''") ) )
return ret
#Group by iSource:way_id
for geom, tags in parsed_gisdata:
way_key = tags.copy()
- way_key = ( way_key[iSource + ':way_id'], tuple( [(k,v) for k,v in way_key.iteritems()] ) )
+ way_key = ( way_key[iSource + ':way_id'], tuple( [(k,v) for k,v in way_key.items()] ) )
if way_key not in waylist:
waylist[way_key] = []
waylist[way_key].append( geom )
ret = {}
- for (way_id, way_key), segments in waylist.iteritems():
+ for (way_id, way_key), segments in waylist.items():
if way_id != blank_way_id:
ret[way_key] = glom_all( segments )
import_guid = time.strftime( '%Y%m%d%H%M%S' )
- print "parsing shpfile"
+ print("parsing shpfile")
parsed_features = parse_shp_for_osm( shp_filename )
- print "compiling nodelist"
+ print("compiling nodelist")
i, nodelist = compile_nodelist( parsed_features )
- print "compiling waylist"
+ print("compiling waylist")
waylist = compile_waylist( parsed_features, blank_way_id )
filenumber = 1
objectcount = 0
seen = {}
- print "preparing address ways"
+ print("preparing address ways")
ret = addressways(waylist, nodelist, i)
osm_filename = "%s%d.osm" % (base_filename, filenumber)
- print "writing %s" %osm_filename
+ print("writing %s" %osm_filename)
fp = open( osm_filename, "w" )
fp.write( "\n".join( ret ) )
fp.close()
filenumber += 1
- print "constructing osm xml file"
+ print("constructing osm xml file")
ret = []
ret.append( "<?xml version='1.0' encoding='UTF-8'?>" )
ret.append( "<osm version='0.6' generator='shape_to_osm.py'>" )
- for waykey, segments in waylist.iteritems():
+ for waykey, segments in waylist.items():
for segment in segments:
- #write the nodes
+ #write the nodes
for point in segment:
id, (lat, lon) = nodelist[ round_point( point ) ]
if id not in seen:
seen[id] = True
- #write node
+ #write node
ret.append( " <node id='-%d' action='create' visible='true' lat='%f' lon='%f' >" % (id, lat, lon) )
ret.append( " </node>" )
objectcount += 1
pass
#print "Skipping node %d" %id
- #write the way
+ #write the way
ret.append( " <way id='-%d' action='create' visible='true'>" % i )
ids = [ nodelist[ round_point( point ) ][0] for point in segment ]
for id in ids:
count += 1
ret.append( " <nd ref='-%d' />" % id )
- if (count % Max_Waylength == 0) and (count != len(ids)): #Split the way
+ if (count % Max_Waylength == 0) and (count != len(ids)): #Split the way
for k, v in waykey:
ret.append( " <tag k=\"%s\" v=\"%s\" />" % (k, escape(str(v))) )
ret.append( " <tag k=\"source\" v=\"%s_import_v%s_%s\" />" % (iSource, VERSION, import_guid) )
ret.append( " </way>" )
objectcount += 1
- i += 1
+ i += 1
ret.append( " <way id='-%d' action='create' visible='true'>" % i )
ret.append( " <nd ref='-%d' />" % id )
i += 1
- if objectcount > maxNodes: #Write a file
+ if objectcount > maxNodes: #Write a file
ret.append( "</osm>" )
osm_filename = "%s%d.osm" % (base_filename, filenumber)
- print "writing %s" %osm_filename
+ print("writing %s" %osm_filename)
fp = open( osm_filename, "w" )
fp.write( "\n".join( ret ) )
fp.close()
ret.append( "</osm>" )
osm_filename = "%s%d.osm" % (base_filename, filenumber)
- print "writing %s" %osm_filename
+ print("writing %s" %osm_filename)
fp = open( osm_filename, "w" )
fp.write( "\n".join( ret ) )
fp.close()
if __name__ == '__main__':
import sys, os.path
if len(sys.argv) < 2:
- print "%s filename.shp [filename.osm]" % sys.argv[0]
+ print("%s filename.shp [filename.osm]" % sys.argv[0])
sys.exit()
shape = sys.argv[1]
if len(sys.argv) > 2:
else:
osm = shape[0:-4] + ".osm"
id = "1.shp"
- # Left over from massGIS unknown usage, but works fine hardcoded to "1.shp" which was the valu on a test of the actual mass data,
- #id = os.path.basename(shape).split("_")[-1]
+ # Left over from massGIS unknown usage, but works fine hardcoded to "1.shp" which was the valu on a test of the actual mass data,
+ #id = os.path.basename(shape).split("_")[-1]
shape_to_osm( shape, osm, id )
projects / nominatim.git / blobdiff