def invent_matrix(self,cksum,variant):
print "Variant:",variant
for ad in ec.get_airfields():
if not 'adcharts' in ad: continue
dbb=False
for adchart in ad['adcharts'].values():
if adchart['checksum']==cksum:
lat,lon=mapper.from_str(ad['pos'])
w,h=adchart['render_width'],adchart['render_height']
print "adpos:",ad['pos'],lat,lon
dbb=True
break
if dbb:break
else:
raise Exception("Can't find this chart in aipdata")
if variant.lower().count('vac'):
mercsizex=w
mercsizey=h
scale=7
else:
mercsizex=w
mercsizey=h
scale=1
print "lat,lon",lat,lon
m1=mapper.latlon2merc((lat,lon),13)
print "corner merc",m1
ns=[]
for offx,offy in [(0,0),
(mercsizex,0),
(0,mercsizey),
(mercsizex,mercsizey)]:
merc2=(m1[0]+(offx-mercsizex/2)*scale,m1[1]+(offy-mercsizey/2)*scale)
n=AirportMarker()
n.latitude,n.longitude=mapper.merc2latlon(merc2,13)
n.x,n.y=(offx,offy)
ns.append(n)
error,A,T=customproj.solve(ns)
matrix=list(A)+list(T)
print "Fake projection:",matrix
return matrix
def invent_matrix(self, cksum, variant):
print "Variant:", variant
for ad in ec.get_airfields():
if not 'adcharts' in ad: continue
dbb = False
for adchart in ad['adcharts'].values():
if adchart['checksum'] == cksum:
lat, lon = mapper.from_str(ad['pos'])
w, h = adchart['render_width'], adchart['render_height']
print "adpos:", ad['pos'], lat, lon
dbb = True
break
if dbb: break
else:
raise Exception("Can't find this chart in aipdata")
if variant.lower().count('vac'):
mercsizex = w
mercsizey = h
scale = 7
else:
mercsizex = w
mercsizey = h
scale = 1
print "lat,lon", lat, lon
m1 = mapper.latlon2merc((lat, lon), 13)
print "corner merc", m1
ns = []
for offx, offy in [(0, 0), (mercsizex, 0), (0, mercsizey),
(mercsizex, mercsizey)]:
merc2 = (m1[0] + (offx - mercsizex / 2) * scale,
m1[1] + (offy - mercsizey / 2) * scale)
n = AirportMarker()
n.latitude, n.longitude = mapper.merc2latlon(merc2, 13)
n.x, n.y = (offx, offy)
ns.append(n)
error, A, T = customproj.solve(ns)
matrix = list(A) + list(T)
print "Fake projection:", matrix
return matrix
def save(self):
print request.params
ad = request.params['ad']
chartobj = None
mapchecksum = request.params['mapchecksum']
for adobj in ec.get_airfields():
if adobj['name'] == ad:
bb = False
for adchart in adobj['adcharts'].values():
if adchart['checksum'] == mapchecksum:
chartobj = adchart
bb = True
break
if bb: break
else:
self.error("No such airport" + ad)
marks = dict()
for param, val in request.params.items():
if param in [
"save", "ad", 'mapchecksum', 'scroll_x', 'scroll_y',
'maptype', 'scale', 'north'
]:
continue
if param.startswith("del"):
continue
if param.startswith("set_"):
x, y = [int(v) for v in param.split("_")[1:]]
session['curadmarker'] = (x, y)
session.save()
continue
if param == "add":
xs = meta.Session.query(AirportMarker.x).filter(
sa.and_(AirportMarker.user == session['user'],
AirportMarker.airport == ad)).all()
if xs:
maxx = max(xs)[0] + 1
else:
maxx = 0
marks[(maxx, 0)] = dict(latitude=None,
longitude=None,
x=maxx,
y=0)
session['curadmarker'] = (maxx, 0)
session.save()
continue
sx, sy, attrib = re.match(ur"mark_(\d+)_(\d+)_(\w*)",
param).groups()
x = int(sx)
y = int(sy)
marks.setdefault((x, y), dict())[attrib] = val
thresholds = dict()
for rwy in adobj.get('runways', []):
ends = rwy['ends']
for end in ends:
thresholds[end['thr']] = mapper.from_str(end['pos'])
for param, val in request.params.items():
if param.startswith("del_"):
x, y = [int(v) for v in param.split("_")[1:]]
marks.pop((x, y))
continue
meta.Session.query(AirportMarker).filter(
sa.and_(AirportMarker.user == session['user'],
AirportMarker.airport == ad)).delete()
ms = []
arppos = mapper.from_str(adobj['pos'])
for (x, y), val in marks.items():
m = AirportMarker()
m.user = session['user']
m.airport = ad
m.mapchecksum = str(mapchecksum)
m.x = int(val['x'])
m.y = int(val['y'])
m.latitude, w1 = parselatlon(val['latitude'], arppos, thresholds,
0)
m.longitude, w2 = parselatlon(val['longitude'], arppos, thresholds,
1)
if w1 or w2:
m.weight = w1 + w2
else:
m.weigth = 1
meta.Session.add(m)
ms.append(m)
proj = meta.Session.query(AirportProjection).filter(
sa.and_(AirportProjection.user == session['user'],
AirportProjection.airport == ad,
AirportProjection.mapchecksum == str(mapchecksum))).one()
try:
proj.scale = float(request.params['scale'].strip())
except:
proj.scale = None
try:
proj.north = float(request.params['north'].strip())
except:
proj.north = None
def both_lat_lon(x):
return x.latitude and x.longitude
def neither_lat_lon(x):
return not x.latitude and not x.longitude
def just_lat(x):
return x.latitude and not x.longitude
def just_lon(x):
return not x.latitude and x.longitude
ms = [m for m in ms if not neither_lat_lon(m)]
"""
if (len(ms)==4 and
len([m for m in ms if just_lat(m)])==2 and
len([m for m in ms if just_lon(m)])==2):
extra=[]
for m in ms:
n=AirportMarker()
n.x=m.x
n.y=m.y
if just_lat(m):
n.latitude=m.latitude
n.x+=1000
extra.append(n)
if just_lon(m):
n.y+=1000
n.longitude=m.longitude
extra.append(n)
ms.extend(extra)
"""
if len(ms) == 1 and both_lat_lon(
ms[0]) and proj.scale and proj.north != None:
print "Scale/north triggered"
print "Adchart:", chartobj
if chartobj != None:
render_height = chartobj['render_height']
mark, = ms
pixelpos = (mark.x, mark.y)
mapsize = adchart.get('mapsize', (210, 297))
mapheight_meter = mapsize[1] / 1000.0 * proj.scale
mapheight_km = mapheight_meter / 1000.0
merc = mapper.latlon2merc((mark.latitude, mark.longitude), 17)
pixels = mapper.approx_scale(merc, 17, mapheight_km / 1.852)
newmerc = (merc[0], merc[1] - pixels)
northrad = proj.north / (180.0 / math.pi)
newpixelpos = (pixelpos[0] +
render_height * math.sin(northrad),
pixelpos[1] -
render_height * math.cos(northrad))
m = AirportMarker()
m.x = newpixelpos[0]
m.y = newpixelpos[1]
latlon = mapper.merc2latlon(newmerc, 17)
m.latitude = latlon[0]
m.longitude = latlon[1]
m.weight = 1
ms.append(m)
if len(ms) == 2 and all(both_lat_lon(x) for x in ms):
print "Have exactly two marks now"
mark1, mark2 = ms
lm1, lm2 = [
mapper.latlon2merc((mark.latitude, mark.longitude), 17)
for mark in [mark1, mark2]
]
ld = (lm2[0] - lm1[0], lm2[1] - lm1[1])
pd = (mark2.x - mark1.x, mark2.y - mark1.y)
lm3 = (lm1[0] - ld[1], lm1[1] + ld[0])
pm3 = (mark1.x - pd[1], mark1.y + pd[0])
llm3 = mapper.merc2latlon(lm3, 17)
m = AirportMarker()
m.x = pm3[0]
m.y = pm3[1]
m.latitude, w1 = llm3[0], 1
m.longitude, w2 = llm3[1], 1
ms.append(m)
print "delta pixels", pd
print "delta latlon", ld
print "extra end pixels", m.x, m.y
print "extra end latlon", m.latitude, m.longitude
eqns = 0
for m in ms:
if both_lat_lon(m): eqns += 2
elif just_lat(m): eqns += 1
elif just_lon(m): eqns += 1
try:
if eqns < 4: raise Exception("Unsolvable")
error, A, T = customproj.solve(ms)
matrix = list(A) + list(T)
if proj.matrix:
oldmatrix = list(proj.matrix)
newmatrix = list(A) + list(T)
diff = sum(abs(a - b) for a, b in zip(oldmatrix, newmatrix))
else:
diff = 1e30 #enough to trigger update
if diff > 1e-12:
proj.matrix = tuple(newmatrix)
proj.updated = datetime.utcnow().replace(microsecond=0)
except Exception, cause:
print "Couldn't solve projection equation %s" % (cause, )
proj.matrix = [1, 0, 0, 1, 0, 0]
proj.updated = datetime.utcnow().replace(microsecond=0)
meta.Session.add(proj)
def save(self):
print request.params
ad=request.params['ad']
chartobj=None
mapchecksum=request.params['mapchecksum']
for adobj in ec.get_airfields():
if adobj['name']==ad:
bb=False
for adchart in adobj['adcharts'].values():
if adchart['checksum']==mapchecksum:
chartobj=adchart
bb=True
break
if bb: break
else:
self.error("No such airport"+ad)
marks=dict()
for param,val in request.params.items():
if param in ["save","ad",'mapchecksum','scroll_x','scroll_y','maptype','scale','north']:
continue
if param.startswith("del"):
continue
if param.startswith("set_"):
x,y=[int(v) for v in param.split("_")[1:]]
session['curadmarker']=(x,y)
session.save()
continue
if param=="add":
xs=meta.Session.query(AirportMarker.x).filter(sa.and_(
AirportMarker.user==session['user'],
AirportMarker.airport==ad
)).all()
if xs:
maxx=max(xs)[0]+1
else:
maxx=0
marks[(maxx,0)]=dict(latitude=None,longitude=None,x=maxx,y=0)
session['curadmarker']=(maxx,0)
session.save()
continue
sx,sy,attrib=re.match(ur"mark_(\d+)_(\d+)_(\w*)",param).groups()
x=int(sx)
y=int(sy)
marks.setdefault((x,y),dict())[attrib]=val
thresholds=dict()
for rwy in adobj.get('runways',[]):
ends=rwy['ends']
for end in ends:
thresholds[end['thr']]=mapper.from_str(end['pos'])
for param,val in request.params.items():
if param.startswith("del_"):
x,y=[int(v) for v in param.split("_")[1:]]
marks.pop((x,y))
continue
meta.Session.query(AirportMarker).filter(sa.and_(
AirportMarker.user==session['user'],
AirportMarker.airport==ad)).delete()
ms=[]
arppos=mapper.from_str(adobj['pos'])
for (x,y),val in marks.items():
m=AirportMarker()
m.user=session['user']
m.airport=ad
m.mapchecksum=str(mapchecksum)
m.x=int(val['x'])
m.y=int(val['y'])
m.latitude,w1=parselatlon(val['latitude'],arppos,thresholds,0)
m.longitude,w2=parselatlon(val['longitude'],arppos,thresholds,1)
if w1 or w2:
m.weight=w1+w2
else:
m.weigth=1
meta.Session.add(m)
ms.append(m)
proj=meta.Session.query(AirportProjection).filter(sa.and_(
AirportProjection.user==session['user'],
AirportProjection.airport==ad,
AirportProjection.mapchecksum==str(mapchecksum))).one()
try:
proj.scale=float(request.params['scale'].strip())
except:
proj.scale=None
try:
proj.north=float(request.params['north'].strip())
except:
proj.north=None
def both_lat_lon(x):
return x.latitude and x.longitude
def neither_lat_lon(x):
return not x.latitude and not x.longitude
def just_lat(x):
return x.latitude and not x.longitude
def just_lon(x):
return not x.latitude and x.longitude
ms=[m for m in ms if not neither_lat_lon(m)]
"""
if (len(ms)==4 and
len([m for m in ms if just_lat(m)])==2 and
len([m for m in ms if just_lon(m)])==2):
extra=[]
for m in ms:
n=AirportMarker()
n.x=m.x
n.y=m.y
if just_lat(m):
n.latitude=m.latitude
n.x+=1000
extra.append(n)
if just_lon(m):
n.y+=1000
n.longitude=m.longitude
extra.append(n)
ms.extend(extra)
"""
if len(ms)==1 and both_lat_lon(ms[0]) and proj.scale and proj.north!=None:
print "Scale/north triggered"
print "Adchart:",chartobj
if chartobj!=None:
render_height=chartobj['render_height']
mark,=ms
pixelpos=(mark.x,mark.y)
mapsize=adchart.get('mapsize',(210,297))
mapheight_meter=mapsize[1]/1000.0 * proj.scale
mapheight_km=mapheight_meter/1000.0
merc=mapper.latlon2merc((mark.latitude,mark.longitude),17)
pixels=mapper.approx_scale(merc,17,mapheight_km/1.852)
newmerc=(merc[0],merc[1]-pixels)
northrad=proj.north/(180.0/math.pi)
newpixelpos=(pixelpos[0]+render_height*math.sin(northrad),
pixelpos[1]-render_height*math.cos(northrad))
m=AirportMarker()
m.x=newpixelpos[0]
m.y=newpixelpos[1]
latlon=mapper.merc2latlon(newmerc,17)
m.latitude=latlon[0]
m.longitude=latlon[1]
m.weight=1
ms.append(m)
if len(ms)==2 and all(both_lat_lon(x) for x in ms):
print "Have exactly two marks now"
mark1,mark2=ms
lm1,lm2=[mapper.latlon2merc((mark.latitude,mark.longitude),17) for mark in [mark1,mark2]]
ld=(lm2[0]-lm1[0],lm2[1]-lm1[1])
pd=(mark2.x-mark1.x,mark2.y-mark1.y)
lm3=(lm1[0]-ld[1],lm1[1]+ld[0])
pm3=(mark1.x-pd[1],mark1.y+pd[0])
llm3=mapper.merc2latlon(lm3,17)
m=AirportMarker()
m.x=pm3[0]
m.y=pm3[1]
m.latitude,w1=llm3[0],1
m.longitude,w2=llm3[1],1
ms.append(m)
print "delta pixels",pd
print "delta latlon",ld
print "extra end pixels",m.x,m.y
print "extra end latlon",m.latitude,m.longitude
eqns=0
for m in ms:
if both_lat_lon(m): eqns+=2
elif just_lat(m): eqns+=1
elif just_lon(m): eqns+=1
try:
if eqns<4: raise Exception("Unsolvable")
error,A,T=customproj.solve(ms)
matrix=list(A)+list(T)
if proj.matrix:
oldmatrix=list(proj.matrix)
newmatrix=list(A)+list(T)
diff=sum(abs(a-b) for a,b in zip(oldmatrix,newmatrix))
else:
diff=1e30 #enough to trigger update
if diff>1e-12:
proj.matrix=tuple(newmatrix)
proj.updated=datetime.utcnow().replace(microsecond=0)
except Exception,cause:
print "Couldn't solve projection equation %s"%(cause,)
proj.matrix=[1,0,0,1,0,0]
proj.updated=datetime.utcnow().replace(microsecond=0)
meta.Session.add(proj)