def run():
sf = shapelib.open(
os.getenv("SWFP_DATADIR") +
"/mapnik_render/world_boundaries/world_boundaries_m.shp")
d = dbflib.DBFFile(
os.getenv("SWFP_DATADIR") +
"/mapnik_render/world_boundaries/world_boundaries_m.dbf")
num_shapes = sf.info()[0]
assert num_shapes == d.record_count()
swedish_polygons = 0
for idx in xrange(num_shapes):
obj = sf.read_object(idx)
rec = d.read_record(idx)
if rec['CNTRY_NAME'] == 'Sweden':
#print "Sweden: ",obj.vertices()
swedish_polygons += 1
assert len(obj.vertices()) == 1
out = []
for vert in obj.vertices()[0]:
cd = prj.inverse(mapnik.Coord(vert[1], vert[0]))
#print "lat: %s, lon: %s,"%(cd.y,cd.x)
out.append(mapper.format_lfv(cd.y, cd.x))
print "Swedpol:", " - ".join(out)
print "Swedish polygons: %d" % (swedish_polygons, )
def run():
sf=shapelib.open(os.getenv("SWFP_DATADIR")+"/mapnik_render/world_boundaries/world_boundaries_m.shp")
d=dbflib.DBFFile(os.getenv("SWFP_DATADIR")+"/mapnik_render/world_boundaries/world_boundaries_m.dbf")
num_shapes=sf.info()[0]
assert num_shapes==d.record_count()
swedish_polygons=0
for idx in xrange(num_shapes):
obj=sf.read_object(idx)
rec=d.read_record(idx)
if rec['CNTRY_NAME']=='Sweden':
#print "Sweden: ",obj.vertices()
swedish_polygons+=1
assert len(obj.vertices())==1
out=[]
for vert in obj.vertices()[0]:
cd=prj.inverse(mapnik.Coord(vert[1],vert[0]))
#print "lat: %s, lon: %s,"%(cd.y,cd.x)
out.append(mapper.format_lfv(cd.y,cd.x))
print "Swedpol:", " - ".join(out)
print "Swedish polygons: %d"%(swedish_polygons,)
supname,=re.match(".*/(SUP_\d+_\d+.pdf)",sup).groups()
if limiter!=None and limiter(supname)==False:
continue
areas.extend(extract_single_sup(base+sup,sup,supname,opening_ours))
if "http://www.lfv.se/AIP/AIP%20Sup/SUP_52_10.pdf" in opening_ours:
opening_ours.remove("http://www.lfv.se/AIP/AIP%20Sup/SUP_52_10.pdf")
if not opening_ours:
opening_ours=['unknown']
return areas,list(opening_ours)[0]
if __name__=='__main__':
def limiter(x):
if len(sys.argv)<=1: return True
if x.count(sys.argv[1]):
return True
return False
if len(sys.argv)<=1:
limiter=None
areas,opening_hours=parse_all_sups(limiter)
f=open("aipsup-out.txt","w")
for area in sorted(areas,key=lambda x:x['name']):
t="%s: coords: %s "%(area['name'],"-".join(mapper.format_lfv(*mapper.from_str(c)) for c in area['points']))
print t
f.write(repr(t)+"\n")
print "Opening_hours url:",opening_hours
f.write("Hours:"+opening_hours)
f.close()
def ats(self):
try:
#waypoints=meta.Session.query(Waypoint).filter(sa.and_(
# Waypoint.user==tripuser(),Waypoint.trip==session['current_trip'])).order_by(Waypoint.ordering).all()
#c.waypoints=[]
self.standard_prep(c)
#print "Routes:",c.route
def break_subtrips(routes):
out = []
T = 0.0
for i, rt in enumerate(routes):
out.append(rt)
if rt.time_hours:
T += rt.time_hours
if rt.b.stay or i == len(routes) - 1:
if len(out):
yield dict(T=T), out
T = 0.0
out = []
def format_cruise(tas):
if tas > 999: tas = 999
if tas < 0: tas = 0
return "N%04d" % (tas, )
def format_alt(alt):
try:
alt = alt.upper().strip()
if alt.startswith("FL"):
ialt = int(float(alt[2:].strip()))
return "F%03d" % (ialt, )
ialt = int(float(alt)) / 100
print "parsed alt %s" % (repr(alt, )), "as", ialt
return "A%03d" % (ialt, )
except Exception:
raise AtsException(
"Bad altitude specification for some leg: <%s>" %
(alt))
c.atstrips = []
last_fuel_left = None
nr_persons = None
for meta, routes in break_subtrips(c.route):
print "===============New subtrip..............."
spaces = set()
fir_whenposname = []
accum_time = 0
#print "broke ruote",meta
if len(routes) == 0: continue
at = dict()
at['T'] = meta['T']
waypoints = [routes[0].a]
for rt in routes:
waypoints.append(rt.b)
wps = []
stay = routes[0].a.stay
dep_ad = "ZZZZ"
dep_ad_name = waypoints[0].waypoint
dep_ad_coords = mapper.format_lfv_ats(
*mapper.from_str(waypoints[0].pos))
dest_ad = "ZZZZ"
dest_ad_name = waypoints[-1].waypoint
dest_ad_coords = mapper.format_lfv_ats(
*mapper.from_str(waypoints[-1].pos))
extra_remarks = []
lastwppos = None
lastaltspeed = None
for i, wp in enumerate(waypoints):
print "Subtrip:", i, wp.waypoint
at['T'] = meta['T']
lat, lon = mapper.from_str(wp.pos)
if lastwppos:
assert i >= 1
curpos = (lat, lon)
crossing1 = airspace.get_fir_crossing(
lastwppos, curpos)
for sub in routes[i - 1].subs:
if crossing1:
posa,posb=mapper.merc2latlon(sub.subposa,13),\
mapper.merc2latlon(sub.subposb,13)
crossing = airspace.get_fir_crossing(
posa, posb)
if crossing:
fir, enterpos = crossing
bearing, along = mapper.bearing_and_distance(
posa, enterpos)
if sub.gs > 1e-6:
curtime = accum_time + along / sub.gs
fir_whenposname.append(
(curtime, enterpos, fir['icao']))
if sub.time != None:
accum_time += sub.time
else:
accum_time = 9999
for space in get_any_space_on_line(lastwppos, curpos):
spaces.add((space['name'],
space.get('floor', "<Unknown>"),
space.get('ceiling', "<Unknown>")))
lastwppos = (lat, lon)
symbolicpos = None
airport = None
if i == 0 or i == len(waypoints) - 1:
for ad in airspace.get_airfields(lat, lon, 11):
if not ad['icao'].upper() in ['ZZZZ', 'ESVF']:
airport = ad
symbolicpos = ad['icao'].upper()
if i == 0:
dep_ad = ad['icao'].upper()
if i == len(waypoints) - 1:
dest_ad = ad['icao'].upper()
break
else:
for sigp in airspace.get_sigpoints(lat, lon, 11):
if sigp['kind'] == "sig. point":
if len(sigp['name']) == 5:
sigfound = sigp
symbolicpos = sigp['name']
break
if symbolicpos == None:
symbolicpos = mapper.format_lfv_ats(lat, lon)
if i < len(routes):
altspeed = (format_alt(routes[i].altitude),
format_cruise(routes[i].tas))
if lastaltspeed != None:
if lastaltspeed != altspeed:
alt, speed = altspeed
symbolicpos += "/" + speed + alt
lastaltspeed = altspeed
wps.append(
dict(name=wp.waypoint,
airport=airport,
symbolicpos="DCT " + symbolicpos,
exactpos=mapper.format_lfv(lat, lon),
decimalpos="%.5f,%.5f" % (lat, lon)))
for when, pos, fir in fir_whenposname:
hour, minute = divmod(int(60 * when), 60)
extra_remarks.append("EET/%s%02d%02d" %
(fir, hour, minute))
if dep_ad == "ZZZZ":
extra_remarks.append(
u"DEP/%s %s" %
(dep_ad_coords, strip_accents(dep_ad_name.upper())))
if dest_ad == "ZZZZ":
extra_remarks.append(
u"DEST/%s %s" %
(dest_ad_coords, strip_accents(dest_ad_name.upper())))
if stay.date_of_flight.strip():
dof = stay.date_of_flight.replace("-", "").strip()
if len(dof) == 8 and dof.startswith("20"):
dof = dof[2:]
else:
dof = routes[0].depart_dt.strftime("%y%m%d")
print "dof:", dof
if len(dof) != 6:
raise AtsException(
u"ATS flight plans need takeoff date for all takeoffs!"
)
else:
extra_remarks.append(u"DOF/%s" % (dof, ))
if stay and stay.nr_persons:
nr_persons = stay.nr_persons
else:
if nr_persons == None:
raise AtsException(
u"You must enter the the number of persons who will be flying!"
)
tas = routes[0].tas
altitude = routes[0].altitude
at['wps'] = wps
if any(rt.time_hours == None for rt in routes):
raise AtsException(
"TAS is less than headwind for some part of trip.")
enroute_time = sum(rt.time_hours for rt in routes)
fuel = last_fuel_left
if stay and stay.fuel:
fuel = stay.fuel
if not c.ac:
raise AtsException(
"You must choose an aircraft type for this journey to be able to create an ATS flight plan"
)
if c.ac.cruise_burn > 1e-3 and fuel:
endurance = float(fuel) / float(c.ac.cruise_burn)
else:
endurance = 0.0
if endurance <= 0.0:
if fuel == None:
raise AtsException(
"Enter a value for 'Fuel at takeoff'!")
else:
raise AtsException(
"You do not have enough fuel for the entire journey! This means your endurance would be 0 or negative for one or more legs. Add a fuel stop, shorten the journey, or bring more fuel!"
)
if not c.user.realname:
raise AtsException(
"You should enter your name under profile settings, for use as the name of the commander in the flight plan"
)
phonenr = ""
if c.user.phonenr:
phonenr = c.user.phonenr
fir_whenposname.sort()
def eqp(x, s):
x = "".join(re.findall('[A-Z]', x.upper()))
if len(x) == 0:
return s
return x
dummy = u"""
FPL-SEVLI-VG
-ULAC/L-V/C
-EFKG1330
-N0075F065 DCT 5959N02016E DCT 5949N01936E DCT 5929N01818E DCT 5927N01742E
-ZZZZ0130
-DEST/5927N01742E FRÖLUNDA RMK/BORDER CROSSING 40MIN AFTER TAKEOFF DOF/101002 ORGN/ESSAZPZX
-E/0300 P/2
A/R W
C/ANDERS MUSIKKA +4670123123"""
atsfplan = u"""
(FPL-%(acreg)s-%(flight_rules)s%(type_of_flight)s
-%(actype)s/%(turbulence_category)s-%(equipment)s/%(transponder)s
-%(dep_ad)s%(eobt)s
-%(cruise_speed)s%(level)s %(route)s DCT
-%(dest_ad)s%(ete)s
-%(extra_remarks)s
-E/%(endurance)s P/%(nr_passengers)s
A/%(markings)s%(extra_equipment)s
C/%(commander)s %(phonenr)s)""" % (dict(
acreg=c.ac.aircraft.replace("-", ""),
actype=c.ac.atstype,
turbulence_category='L',
flight_rules='V',
type_of_flight='G',
equipment=eqp(c.ac.com_nav_equipment, 'V'),
transponder=eqp(c.ac.transponder_equipment, 'C'),
extra_equipment=u" %s" %
(c.ac.extra_equipment, ) if c.ac.extra_equipment else "",
dep_ad=dep_ad,
eobt=routes[0].depart_dt.strftime("%H%M"),
cruise_speed=format_cruise(tas),
level=format_alt(altitude),
route=(" ".join("%s" % (w['symbolicpos'], )
for w in wps[1:-1])),
dest_ad=dest_ad,
ete=lfvclockfmt(enroute_time),
extra_remarks=" ".join(extra_remarks),
endurance=lfvclockfmt(endurance),
nr_passengers=nr_persons,
markings=c.ac.markings,
commander=strip_accents(c.user.realname if c.user.realname
else u"UNKNOWN").replace(" ", ""),
phonenr=c.user.phonenr if c.user.phonenr else ""))
at['atsfplan'] = atsfplan.strip()
#print "Adding atstrip:",atsfplan
at['spacesummary'] = spaces
last_fuel_left = routes[-1].accum_fuel_left
c.atstrips.append(at)
c.atstrips = [at for at in c.atstrips if len(at['wps'])]
#response.headers['Content-Type'] = 'application/xml'
return render('/ats.mako')
except AtsException, ats:
redirect(
h.url_for(controller='flightplan',
action="index",
flash=unicode(ats)))
def extract_airfields(filtericao=lambda x:True,purge=True):
#print getxml("/AIP/AD/AD 1/ES_AD_1_1_en.pdf")
ads=[]
p=Parser("/AIP/AD/AD 1/ES_AD_1_1_en.pdf")
points=dict()
startpage=None
for pagenr in xrange(p.get_num_pages()):
page=p.parse_page_to_items(pagenr)
if page.count("Aerodrome directory"):
startpage=pagenr
break
if startpage==None:
raise Exception("Couldn't find aerodrome directory in file")
#print "Startpage: %d"%(startpage,)
#nochartf=open("nochart.txt","w")
for pagenr in xrange(startpage,p.get_num_pages()):
row_y=[]
page=p.parse_page_to_items(pagenr)
allines=[x for x in (page.get_lines(page.get_partially_in_rect(0,0,15,100))) if x.strip()]
for item,next in zip(allines,allines[1:]+[""]):
#print "item:",item
m=re.match(ur"^\s*[A-ZÅÄÖ]{3,}(?:/.*)?\b.*",item)
if m:
#print "Candidate, next is:",next
if re.match(r"^\s*[A-Z]{4}\b.*",next):
#print "Matched:",item
#print "y1:",item.y1
row_y.append(item.y1)
for y1,y2 in zip(row_y,row_y[1:]+[100.0]):
#print "Extacting from y-range: %f-%f"%(y1,y2)
items=list(page.get_partially_in_rect(0,y1-0.25,5.0,y2+0.25,ysort=True))
if len(items)>=2:
#print "Extract items",items
ad=dict(name=unicode(items[0].text).strip(),
icao=unicode(items[1].text).strip()
)
#print "Icao:",ad['icao']
assert re.match(r"[A-Z]{4}",ad['icao'])
if not filtericao(ad): continue
if len(items)>=3:
#print "Coord?:",items[2].text
m=re.match(r".*(\d{6}N)\s*(\d{7}E).*",items[2].text)
if m:
lat,lon=m.groups()
ad['pos']=parse_coords(lat,lon)
#print "Items3:",items[3:]
elev=re.findall(r"(\d{1,5})\s*ft"," ".join(t.text for t in items[3:]))
#print "Elev:",elev
assert len(elev)==1
ad['elev']=int(elev[0])
ads.append(ad)
big_ad=set()
for ad in ads:
if not ad.has_key('pos'):
big_ad.add(ad['icao'])
for ad in ads:
icao=ad['icao']
if icao in big_ad:
if icao in ['ESIB','ESNY','ESCM','ESPE']:
continue
try:
p=Parser("/AIP/AD/AD 2/%s/ES_AD_2_%s_6_1_en.pdf"%(icao,icao))
except:
p=Parser("/AIP/AD/AD 2/%s/ES_AD_2_%s_6-1_en.pdf"%(icao,icao))
ad['aipvacurl']=p.get_url()
for pagenr in xrange(p.get_num_pages()):
page=p.parse_page_to_items(pagenr)
"""
for altline in exitlines:
m=re.match(r"(\w+)\s+(\d+N)\s*(\d+E.*)",altline)
if not m: continue
name,lat,lon=m.groups()
try:
coord=parse_coords(lat,lon)
except Exception:
continue
points.append(dict(name=name,pos=coord))
"""
for kind in xrange(2):
if kind==0:
hits=page.get_by_regex(r"H[Oo][Ll][Dd][Ii][Nn][Gg]")
kind="holding point"
if kind==1:
hits=page.get_by_regex(r"[Ee]ntry.*[Ee]xit.*point")
kind="entry/exit point"
if len(hits)==0: continue
for holdingheading in hits:
items=sorted(page.get_partially_in_rect(holdingheading.x1+2.0,holdingheading.y2+0.1,holdingheading.x1+0.5,100),
key=lambda x:x.y1)
items=[x for x in items if not x.text.startswith(" ")]
#print "Holding items:",items
for idx,item in enumerate(items):
print "Holding item",item
y1=item.y1
if idx==len(items)-1:
y2=100
else:
y2=items[idx+1].y1
items2=[x for x in page.get_partially_in_rect(item.x1+1,y1+0.3,item.x1+40,y2-0.1) if x.x1>=item.x1-0.25 and x.y1>=y1-0.05 and x.y1<y2-0.05]
s=(" ".join(page.get_lines(items2))).strip()
print "Holding lines:",repr(page.get_lines(items2))
#if s.startswith("ft Left/3"): #Special case for ESOK
# s,=re.match("ft Left/3.*?([A-Z]{4,}.*)",s).groups()
#m=re.match("ft Left/\d+.*?([A-Z]{4,}.*)",s)
#if m:
# s,=m.groups()
if s.startswith("LjUNG"): #Really strange problem with ESCF
s=s[0]+"J"+s[2:]
if s.lower().startswith("holding"):
sl=s.split(" ",1)
if len(sl)>1:
s=sl[1]
s=s.strip()
if kind=="entry/exit point" and s.startswith("HOLDING"):
continue #reached HOLDING-part of VAC
#Check for other headings
#Fixup strange formatting of points in some holding items: (whitespace between coord and 'E')
s=re.sub(ur"(\d+)\s*(N)\s*(\d+)\s*(E)",lambda x:"".join(x.groups()),s)
m=re.match(r"([A-Z]{2,}).*?(\d+N)\s*(\d+E).*",s)
if not m:
m=re.match(r".*?(\d+N)\s*(\d+E).*",s)
if not m:
continue
assert m
lat,lon=m.groups()
#skavsta
if icao=="ESKN":
if s.startswith(u"Hold north of T"):
name="NORTH"
elif s.startswith(u"Hold south of B"):
name="SOUTH"
else:
assert 0
#add more specials here
else:
continue
else:
name,lat,lon=m.groups()
try:
coord=parse_coords(lat,lon)
except Exception:
print "Couldn't parse:",lat,lon
continue
#print name,lat,lon,mapper.format_lfv(*mapper.from_str(coord))
if name.count("REMARK") or len(name)<=2:
print "Suspicious name: ",name
#sys.exit(1)
continue
points[icao+' '+name]=dict(name=icao+' '+name,icao=icao,pos=coord,kind=kind)
#for point in points.items():
# print point
#sys.exit(1)
def fixhex11(s):
out=[]
for c in s:
i=ord(c)
if i>=0x20:
out.append(c)
continue
if i in [0x9,0xa,0xd]:
out.append(c)
continue
out.append(' ')
return "".join(out)
for ad in ads:
icao=ad['icao']
if icao in big_ad:
#print "Parsing ",icao
p=Parser("/AIP/AD/AD 2/%s/ES_AD_2_%s_en.pdf"%(icao,icao),loadhook=fixhex11)
ad['aiptexturl']=p.get_url()
firstpage=p.parse_page_to_items(0)
te="\n".join(firstpage.get_all_lines())
#print te
coords=re.findall(r"ARP.*(\d{6}N)\s*(\d{7}E)",te)
if len(coords)>1:
raise Exception("First page of airport info (%s) does not contain exactly ONE set of coordinates"%(icao,))
if len(coords)==0:
print "Couldn't find coords for ",icao
#print "Coords:",coords
ad['pos']=parse_coords(*coords[0])
elev=re.findall(r"Elevation.*?(\d{1,5})\s*ft",te,re.DOTALL)
if len(elev)>1:
raise Exception("First page of airport info (%s) does not contain exactly ONE elevation in ft"%(icao,))
if len(elev)==0:
print "Couldn't find elev for ",icao
ad['elev']=int(elev[0])
freqs=[]
found=False
thrs=[]
#uprint("-------------------------------------")
for pagenr in xrange(p.get_num_pages()):
page=p.parse_page_to_items(pagenr)
#uprint("Looking on page %d"%(pagenr,))
if 0: #opening hours are no longer stored in a separate document for any airports. No need to detect which any more (since none are).
for item in page.get_by_regex(".*OPERATIONAL HOURS.*"):
lines=page.get_lines(page.get_partially_in_rect(0,item.y2+0.1,100,100))
for line in lines:
things=["ATS","Fuelling","Operating"]
if not line.count("AIP SUP"): continue
for thing in things:
if line.count(thing):
ad['aipsup']=True
for item in page.get_by_regex(".*\s*RUNWAY\s*PHYSICAL\s*CHARACTERISTICS\s*.*"):
#uprint("Physical char on page")
lines=page.get_lines(page.get_partially_in_rect(0,item.y2+0.1,100,100))
seen_end_rwy_text=False
for line,nextline in izip(lines,lines[1:]+[None]):
#uprint("MAtching: <%s>"%(line,))
if re.match(ur"AD\s+2.13",line): break
if line.count("Slope of"): break
if line.lower().count("end rwy:"): seen_end_rwy_text=True
if line.lower().count("bgn rwy:"): seen_end_rwy_text=True
m=re.match(ur".*(\d{6}\.\d+)[\s\(\)\*]*(N).*",line)
if not m:continue
m2=re.match(ur".*(\d{6,7}\.\d+)\s*[\s\(\)\*]*(E).*",nextline)
if not m2:continue
latd,n=m.groups()
lond,e=m2.groups()
assert n=="N"
assert e=="E"
lat=latd+n
lon=lond+e
rwytxts=page.get_lines(page.get_partially_in_rect(0,line.y1+0.05,12,nextline.y2-0.05))
uprint("Rwytxts:",rwytxts)
rwy=None
for rwytxt in rwytxts:
#uprint("lat,lon:%s,%s"%(lat,lon))
#uprint("rwytext:",rwytxt)
m=re.match(ur"\s*(\d{2}[LRCM]?)\b.*",rwytxt)
if m:
assert rwy==None
rwy=m.groups()[0]
if rwy==None and seen_end_rwy_text:
continue
print "Cur airport:",icao
already=False
assert rwy!=None
seen_end_rwy_text=False
for thr in thrs:
if thr['thr']==rwy:
raise Exception("Same runway twice on airfield:"+icao)
thrs.append(dict(pos=mapper.parse_coords(lat,lon),thr=rwy))
assert len(thrs)>=2
for pagenr in xrange(0,p.get_num_pages()):
page=p.parse_page_to_items(pagenr)
matches=page.get_by_regex(r".*ATS\s+COMMUNICATION\s+FACILITIES.*")
#print "Matches of ATS COMMUNICATION FACILITIES on page %d: %s"%(pagenr,matches)
if len(matches)>0:
commitem=matches[0]
curname=None
callsign=page.get_by_regex_in_rect(ur"Call\s*sign",0,commitem.y1,100,commitem.y2+8)[0]
for idx,item in enumerate(page.get_lines(page.get_partially_in_rect(callsign.x1-0.5,commitem.y1,100,100),fudge=0.3,order_fudge=15)):
if item.strip()=="":
curname=None
if re.match(".*RADIO\s+NAVIGATION\s+AND\s+LANDING\s+AIDS.*",item):
break
#print "Matching:",item
m=re.match(r"(.*?)\s*(\d{3}\.\d{1,3})\s*MHz.*",item)
#print "MHZ-match:",m
if not m: continue
#print "MHZ-match:",m.groups()
who,sfreq=m.groups()
freq=float(sfreq)
if abs(freq-121.5)<1e-4:
if who.strip():
curname=who
continue #Ignore emergency frequency, it is understood
if not who.strip():
if curname==None: continue
else:
curname=who
freqs.append((curname.strip().rstrip("/"),freq))
for pagenr in xrange(0,p.get_num_pages()):
page=p.parse_page_to_items(pagenr)
matches=page.get_by_regex(r".*ATS\s*AIRSPACE.*")
#print "Matches of ATS_AIRSPACE on page %d: %s"%(pagenr,matches)
if len(matches)>0:
heading=matches[0]
desigitem,=page.get_by_regex("Designation and lateral limits")
vertitem,=page.get_by_regex("Vertical limits")
airspaceclass,=page.get_by_regex("Airspace classification")
lastname=None
subspacelines=dict()
subspacealts=dict()
for idx,item in enumerate(page.get_lines(page.get_partially_in_rect(desigitem.x2+1,desigitem.y1,100,vertitem.y1-1))):
if item.count("ATS airspace not established"):
assert idx==0
break
if item.strip()=="": continue
m=re.match(r"(.*?)(\d{6}N\s+.*)",item)
if m:
name,coords=m.groups()
name=name.strip()
else:
name=item.strip()
coords=None
if name:
lastname=name
if coords:
subspacelines.setdefault(lastname,[]).append(coords)
assert lastname
lastname=None
#print "Spaces:",subspacelines
#print "ICAO",ad['icao']
#altlines=page.get_lines(page.get_partially_in_rect(vertitem.x2+1,vertitem.y1,100,airspaceclass.y1-0.2))
#print "Altlines:",altlines
subspacealts=dict()
subspacekeys=subspacelines.keys()
allaltlines=" ".join(page.get_lines(page.get_partially_in_rect(vertitem.x1+0.5,vertitem.y1+0.5,100,airspaceclass.y1-0.2)))
single_vertlim=False
totalts=list(mapper.parse_all_alts(allaltlines))
#print "totalts:",totalts
if len(totalts)==2:
single_vertlim=True
for subspacename in subspacekeys:
ceil=None
floor=None
subnames=[subspacename]
if subspacename.split(" ")[-1].strip() in ["TIA","TIZ","CTR","CTR/TIZ"]:
subnames.append(subspacename.split(" ")[-1].strip())
#print "Parsing alts for ",subspacename,subnames
try:
for nametry in subnames:
if single_vertlim: #there's only one subspace, parse all of vertical limits field for this single one.
items=[vertitem]
else:
items=page.get_by_regex_in_rect(nametry,vertitem.x2+1,vertitem.y1,100,airspaceclass.y1-0.2)
for item in items:
alts=[]
for line in page.get_lines(page.get_partially_in_rect(item.x1+0.5,item.y1+0.5,100,airspaceclass.y1-0.2)):
#print "Parsing:",line
line=line.replace(nametry,"").lower().strip()
parsed=list(mapper.parse_all_alts(line))
if len(parsed):
alts.append(mapper.altformat(*parsed[0]))
if len(alts)==2: break
if alts:
#print "alts:",alts
ceil,floor=alts
raise StopIteration
except StopIteration:
pass
assert ceil and floor
subspacealts[subspacename]=dict(ceil=ceil,floor=floor)
spaces=[]
for spacename in subspacelines.keys():
altspacename=spacename
#print "Altspacename: %s, subspacesalts: %s"%(altspacename,subspacealts)
space=dict(
name=spacename,
ceil=subspacealts[altspacename]['ceil'],
floor=subspacealts[altspacename]['floor'],
points=parse_coord_str(" ".join(subspacelines[spacename])),
freqs=list(set(freqs))
)
if True:
vs=[]
for p in space['points']:
x,y=mapper.latlon2merc(mapper.from_str(p),13)
vs.append(Vertex(int(x),int(y)))
p=Polygon(vvector(vs))
if p.calc_area()<=30*30:
pass#print space
pass#print "Area:",p.calc_area()
assert p.calc_area()>30*30
#print "Area: %f"%(p.calc_area(),)
spaces.append(space)
#print space
ad['spaces']=spaces
found=True
if found:
break
assert found
ad['runways']=rwy_constructor.get_rwys(thrs)
#Now find any ATS-airspace
chartblobnames=[]
for ad in ads:
icao=ad['icao']
if icao in big_ad:
parse_landing_chart.help_plc(ad,"/AIP/AD/AD 2/%s/ES_AD_2_%s_2-1_en.pdf"%(icao,icao),
icao,ad['pos'],"se",variant="")
parse_landing_chart.help_plc(ad,"/AIP/AD/AD 2/%s/ES_AD_2_%s_6-1_en.pdf"%(icao,icao),
icao,ad['pos'],"se",variant="vac")
parse_landing_chart.help_plc(ad,"/AIP/AD/AD 2/%s/ES_AD_2_%s_2-3_en.pdf"%(icao,icao),
icao,ad['pos'],"se",variant="parking")
#aip_text_documents.help_parse_doc(ad,"/AIP/AD/AD 2/%s/ES_AD_2_%s_6_1_en.pdf"%(icao,icao),
# icao,"se",title="General Information",category="general")
aip_text_documents.help_parse_doc(ad,"/AIP/AD/AD 2/%s/ES_AD_2_%s_en.pdf"%(icao,icao),
icao,"se",title="General Information",category="general")
#if purge:
# parse_landing_chart.purge_old(chartblobnames,country="se")
#sys.exit(1)
for extra in extra_airfields.extra_airfields:
if filtericao(extra):
ads.append(extra)
print
print
for k,v in sorted(points.items()):
print k,v,mapper.format_lfv(*mapper.from_str(v['pos']))
#print "Num points:",len(points)
origads=list(ads)
for flygkartan_id,name,lat,lon,dummy in csv.reader(open("fplan/extract/flygkartan.csv"),delimiter=";"):
found=None
lat=float(lat)
lon=float(lon)
if type(name)==str:
name=unicode(name,'utf8')
mercf=mapper.latlon2merc((lat,lon),13)
for a in origads:
merca=mapper.latlon2merc(mapper.from_str(a['pos']),13)
dist=math.sqrt((merca[0]-mercf[0])**2+(merca[1]-mercf[1])**2)
if dist<120:
found=a
break
if found:
found['flygkartan_id']=flygkartan_id
else:
d=dict(
icao='ZZZZ',
name=name,
pos=mapper.to_str((lat,lon)),
elev=int(get_terrain_elev((lat,lon))),
flygkartan_id=flygkartan_id)
if filtericao(d):
ads.append(d)
minor_ad_charts=extra_airfields.minor_ad_charts
for ad in ads:
if ad['name'].count(u"Långtora"):
ad['pos']=mapper.to_str(mapper.from_aviation_format("5944.83N01708.20E"))
if ad['name'] in minor_ad_charts:
charturl=minor_ad_charts[ad['name']]
arp=ad['pos']
if 'icao' in ad and ad['icao'].upper()!='ZZZZ':
icao=ad['icao'].upper()
else:
icao=ad['fake_icao']
parse_landing_chart.help_plc(ad,charturl,icao,arp,country='raw',variant="landing")
"""
assert icao!=None
lc=parse_landing_chart.parse_landing_chart(
charturl,
icao=icao,
arppos=arp,country="raw")
assert lc
if lc:
ad['adcharturl']=lc['url']
ad['adchart']=lc
"""
#print ads
for ad in ads:
print "%s: %s - %s (%s ft) (%s)"%(ad['icao'],ad['name'],ad['pos'],ad['elev'],ad.get('flygkartan_id','inte i flygkartan'))
for space in ad.get('spaces',[]):
for freq in space.get('freqs',[]):
print " ",freq
#if 'spaces' in ad:
# print " spaces: %s"%(ad['spaces'],)
#if 'aiptext' in ad:
# print "Aip texts:",ad['aiptext']
#else:
# print "No aiptext"
print "Points:"
for point in sorted(points.values(),key=lambda x:x['name']):
print point
f=codecs.open("extract_airfields.regress.txt","w",'utf8')
for ad in ads:
r=repr(ad)
d=md5.md5(r).hexdigest()
f.write("%s - %s - %s\n"%(ad['icao'],ad['name'],d))
f.close()
f=codecs.open("extract_airfields.regress-details.txt","w",'utf8')
for ad in ads:
r=repr(ad)
f.write(u"%s - %s - %s\n"%(ad['icao'],ad['name'],r))
f.close()
return ads,points.values()
areas.extend(extract_single_sup(base + sup, sup, supname,
opening_ours))
if "http://www.lfv.se/AIP/AIP%20Sup/SUP_52_10.pdf" in opening_ours:
opening_ours.remove("http://www.lfv.se/AIP/AIP%20Sup/SUP_52_10.pdf")
if not opening_ours:
opening_ours = ['unknown']
return areas, list(opening_ours)[0]
if __name__ == '__main__':
def limiter(x):
if len(sys.argv) <= 1: return True
if x.count(sys.argv[1]):
return True
return False
if len(sys.argv) <= 1:
limiter = None
areas, opening_hours = parse_all_sups(limiter)
f = open("aipsup-out.txt", "w")
for area in sorted(areas, key=lambda x: x['name']):
t = "%s: coords: %s " % (area['name'], "-".join(
mapper.format_lfv(*mapper.from_str(c)) for c in area['points']))
print t
f.write(repr(t) + "\n")
print "Opening_hours url:", opening_hours
f.write("Hours:" + opening_hours)
f.close()
title="General Information",
category="general",
)
# if purge:
# parse_landing_chart.purge_old(chartblobnames,country="se")
# sys.exit(1)
for extra in extra_airfields.extra_airfields:
if filtericao(extra):
ads.append(extra)
print
print
for k, v in sorted(points.items()):
print k, v, mapper.format_lfv(*mapper.from_str(v["pos"]))
# print "Num points:",len(points)
origads = list(ads)
for flygkartan_id, name, lat, lon, dummy in csv.reader(open("fplan/extract/flygkartan.csv"), delimiter=";"):
found = None
lat = float(lat)
lon = float(lon)
if type(name) == str:
name = unicode(name, "utf8")
mercf = mapper.latlon2merc((lat, lon), 13)
for a in origads:
merca = mapper.latlon2merc(mapper.from_str(a["pos"]), 13)
dist = math.sqrt((merca[0] - mercf[0]) ** 2 + (merca[1] - mercf[1]) ** 2)
if dist < 120:
def ats(self):
try:
#waypoints=meta.Session.query(Waypoint).filter(sa.and_(
# Waypoint.user==tripuser(),Waypoint.trip==session['current_trip'])).order_by(Waypoint.ordering).all()
#c.waypoints=[]
self.standard_prep(c)
#print "Routes:",c.route
def break_subtrips(routes):
out=[]
T=0.0
for i,rt in enumerate(routes):
out.append(rt)
if rt.time_hours:
T+=rt.time_hours
if rt.b.stay or i==len(routes)-1:
if len(out):
yield dict(T=T),out
T=0.0
out=[]
def format_cruise(tas):
if tas>999: tas=999
if tas<0: tas=0
return "N%04d"%(tas,)
def format_alt(alt):
try:
alt=alt.upper().strip()
if alt.startswith("FL"):
ialt=int(float(alt[2:].strip()))
return "F%03d"%(ialt,)
ialt=int(float(alt))/100
print "parsed alt %s"%(repr(alt,)),"as",ialt
return "A%03d"%(ialt,)
except Exception:
raise AtsException("Bad altitude specification for some leg: <%s>"%(alt))
c.atstrips=[]
last_fuel_left=None
nr_persons=None
for meta,routes in break_subtrips(c.route):
print "===============New subtrip..............."
spaces=set()
fir_whenposname=[]
accum_time=0
#print "broke ruote",meta
if len(routes)==0: continue
at=dict()
at['T']=meta['T']
waypoints=[routes[0].a]
for rt in routes:
waypoints.append(rt.b)
wps=[]
stay=routes[0].a.stay
dep_ad="ZZZZ"
dep_ad_name=waypoints[0].waypoint
dep_ad_coords=mapper.format_lfv_ats(*mapper.from_str(waypoints[0].pos))
dest_ad="ZZZZ"
dest_ad_name=waypoints[-1].waypoint
dest_ad_coords=mapper.format_lfv_ats(*mapper.from_str(waypoints[-1].pos))
extra_remarks=[]
lastwppos=None
lastaltspeed=None
for i,wp in enumerate(waypoints):
print "Subtrip:",i,wp.waypoint
at['T']=meta['T']
lat,lon=mapper.from_str(wp.pos)
if lastwppos:
assert i>=1
curpos=(lat,lon)
crossing1=airspace.get_fir_crossing(lastwppos,curpos)
for sub in routes[i-1].subs:
if crossing1:
posa,posb=mapper.merc2latlon(sub.subposa,13),\
mapper.merc2latlon(sub.subposb,13)
crossing=airspace.get_fir_crossing(posa,posb)
if crossing:
fir,enterpos=crossing
bearing,along=mapper.bearing_and_distance(posa,enterpos)
if sub.gs>1e-6:
curtime=accum_time+along/sub.gs
fir_whenposname.append((curtime,enterpos,fir['icao']))
if sub.time!=None:
accum_time+=sub.time
else:
accum_time=9999
for space in get_any_space_on_line(lastwppos,curpos):
spaces.add((space['name'],space.get('floor',"<Unknown>"),space.get('ceiling',"<Unknown>")))
lastwppos=(lat,lon)
symbolicpos=None
airport=None
if i==0 or i==len(waypoints)-1:
for ad in airspace.get_airfields(lat,lon,11):
if not ad['icao'].upper() in ['ZZZZ','ESVF']:
airport=ad
symbolicpos=ad['icao'].upper()
if i==0:
dep_ad=ad['icao'].upper()
if i==len(waypoints)-1:
dest_ad=ad['icao'].upper()
break
else:
for sigp in airspace.get_sigpoints(lat,lon,11):
if sigp['kind']=="sig. point":
if len(sigp['name'])==5:
sigfound=sigp
symbolicpos=sigp['name']
break
if symbolicpos==None:
symbolicpos=mapper.format_lfv_ats(lat,lon)
if i<len(routes):
altspeed=(format_alt(routes[i].altitude),format_cruise(routes[i].tas))
if lastaltspeed!=None:
if lastaltspeed!=altspeed:
alt,speed=altspeed
symbolicpos+="/"+speed+alt
lastaltspeed=altspeed
wps.append(dict(
name=wp.waypoint,
airport=airport,
symbolicpos="DCT "+symbolicpos,
exactpos=mapper.format_lfv(lat,lon),
decimalpos="%.5f,%.5f"%(lat,lon)
))
for when,pos,fir in fir_whenposname:
hour,minute=divmod(int(60*when),60)
extra_remarks.append("EET/%s%02d%02d"%(fir,hour,minute))
if dep_ad=="ZZZZ":
extra_remarks.append(u"DEP/%s %s"%(dep_ad_coords,strip_accents(dep_ad_name.upper())))
if dest_ad=="ZZZZ":
extra_remarks.append(u"DEST/%s %s"%(dest_ad_coords,strip_accents(dest_ad_name.upper())))
if stay.date_of_flight.strip():
dof=stay.date_of_flight.replace("-","").strip()
if len(dof)==8 and dof.startswith("20"):
dof=dof[2:]
else:
dof=routes[0].depart_dt.strftime("%y%m%d")
print "dof:",dof
if len(dof)!=6:
raise AtsException(u"ATS flight plans need takeoff date for all takeoffs!")
else:
extra_remarks.append(u"DOF/%s"%(dof,))
if stay and stay.nr_persons:
nr_persons=stay.nr_persons
else:
if nr_persons==None:
raise AtsException(u"You must enter the the number of persons who will be flying!")
tas=routes[0].tas
altitude=routes[0].altitude
at['wps']=wps
if any(rt.time_hours==None for rt in routes):
raise AtsException("TAS is less than headwind for some part of trip.")
enroute_time=sum(rt.time_hours for rt in routes)
fuel=last_fuel_left
if stay and stay.fuel:
fuel=stay.fuel
if not c.ac:
raise AtsException("You must choose an aircraft type for this journey to be able to create an ATS flight plan")
if c.ac.cruise_burn>1e-3 and fuel:
endurance=float(fuel)/float(c.ac.cruise_burn)
else:
endurance=0.0
if endurance<=0.0:
if fuel==None:
raise AtsException("Enter a value for 'Fuel at takeoff'!")
else:
raise AtsException("You do not have enough fuel for the entire journey! This means your endurance would be 0 or negative for one or more legs. Add a fuel stop, shorten the journey, or bring more fuel!")
if not c.user.realname:
raise AtsException("You should enter your name under profile settings, for use as the name of the commander in the flight plan")
phonenr=""
if c.user.phonenr:
phonenr=c.user.phonenr
fir_whenposname.sort()
def eqp(x,s):
x="".join(re.findall('[A-Z]',x.upper()))
if len(x)==0:
return s
return x
dummy=u"""
FPL-SEVLI-VG
-ULAC/L-V/C
-EFKG1330
-N0075F065 DCT 5959N02016E DCT 5949N01936E DCT 5929N01818E DCT 5927N01742E
-ZZZZ0130
-DEST/5927N01742E FRÖLUNDA RMK/BORDER CROSSING 40MIN AFTER TAKEOFF DOF/101002 ORGN/ESSAZPZX
-E/0300 P/2
A/R W
C/ANDERS MUSIKKA +4670123123"""
atsfplan=u"""
(FPL-%(acreg)s-%(flight_rules)s%(type_of_flight)s
-%(actype)s/%(turbulence_category)s-%(equipment)s/%(transponder)s
-%(dep_ad)s%(eobt)s
-%(cruise_speed)s%(level)s %(route)s DCT
-%(dest_ad)s%(ete)s
-%(extra_remarks)s
-E/%(endurance)s P/%(nr_passengers)s
A/%(markings)s%(extra_equipment)s
C/%(commander)s %(phonenr)s)"""%(dict(
acreg=c.ac.aircraft.replace("-",""),
actype=c.ac.atstype,
turbulence_category='L',
flight_rules='V',
type_of_flight='G',
equipment=eqp(c.ac.com_nav_equipment,'V'),
transponder=eqp(c.ac.transponder_equipment,'C'),
extra_equipment=u" %s"%(c.ac.extra_equipment,) if c.ac.extra_equipment else "",
dep_ad=dep_ad,
eobt=routes[0].depart_dt.strftime("%H%M"),
cruise_speed=format_cruise(tas),
level=format_alt(altitude),
route=(" ".join("%s"%(w['symbolicpos'],) for w in wps[1:-1])),
dest_ad=dest_ad,
ete=lfvclockfmt(enroute_time),
extra_remarks=" ".join(extra_remarks),
endurance=lfvclockfmt(endurance),
nr_passengers=nr_persons,
markings=c.ac.markings,
commander=strip_accents(c.user.realname if c.user.realname else u"UNKNOWN").replace(" ",""),
phonenr=c.user.phonenr if c.user.phonenr else ""))
at['atsfplan']=atsfplan.strip()
#print "Adding atstrip:",atsfplan
at['spacesummary']=spaces
last_fuel_left=routes[-1].accum_fuel_left
c.atstrips.append(at)
c.atstrips=[at for at in c.atstrips if len(at['wps'])]
#response.headers['Content-Type'] = 'application/xml'
return render('/ats.mako')
except AtsException,ats:
redirect(h.url_for(controller='flightplan',action="index",flash=unicode(ats)))