Implement a persistent cache for navigation data.
[fg:flightgear.git] / src / GUI / MapWidget.cxx
1 #ifdef HAVE_CONFIG_H
2 #  include "config.h"
3 #endif
4
5 #include "MapWidget.hxx"
6
7 #include <sstream>
8 #include <algorithm> // for std::sort
9 #include <plib/puAux.h>
10
11 #include <simgear/sg_inlines.h>
12 #include <simgear/misc/strutils.hxx>
13 #include <simgear/magvar/magvar.hxx>
14 #include <simgear/timing/sg_time.hxx> // for magVar julianDate
15 #include <simgear/structure/exception.hxx>
16
17 #include <Main/globals.hxx>
18 #include <Main/fg_props.hxx>
19 #include <Autopilot/route_mgr.hxx>
20 #include <Navaids/positioned.hxx>
21 #include <Navaids/navrecord.hxx>
22 #include <Navaids/navlist.hxx>
23 #include <Navaids/fix.hxx>
24 #include <Airports/simple.hxx>
25 #include <Airports/runways.hxx>
26 #include <Main/fg_os.hxx>      // fgGetKeyModifiers()
27 #include <Navaids/routePath.hxx>
28
29 const char* RULER_LEGEND_KEY = "ruler-legend";
30
31 /* equatorial and polar earth radius */
32 const float rec  = 6378137;          // earth radius, equator (?)
33 const float rpol = 6356752.314f;      // earth radius, polar   (?)
34
35 /************************************************************************
36   some trigonometric helper functions
37   (translated more or less directly from Alexei Novikovs perl original)
38 *************************************************************************/
39
40 //Returns Earth radius at a given latitude (Ellipsoide equation with two equal axis)
41 static float earth_radius_lat( float lat )
42 {
43   double a = cos(lat)/rec;
44   double b = sin(lat)/rpol;
45   return 1.0f / sqrt( a * a + b * b );
46 }
47
48 ///////////////////////////////////////////////////////////////////////////
49
50 static puBox makePuBox(int x, int y, int w, int h)
51 {
52   puBox r;
53   r.min[0] = x;
54   r.min[1] = y;
55   r.max[0] =  x + w;
56   r.max[1] = y + h;
57   return r;
58 }
59
60 static bool puBoxIntersect(const puBox& a, const puBox& b)
61 {
62   int x0 = SG_MAX2(a.min[0], b.min[0]);
63   int y0 = SG_MAX2(a.min[1], b.min[1]);
64   int x1 = SG_MIN2(a.max[0], b.max[0]);
65   int y1 = SG_MIN2(a.max[1], b.max[1]);
66
67   return (x0 <= x1) && (y0 <= y1);
68 }
69
70 class MapData;
71 typedef std::vector<MapData*> MapDataVec;
72
73 class MapData
74 {
75 public:
76   static const int HALIGN_LEFT = 1;
77   static const int HALIGN_CENTER = 2;
78   static const int HALIGN_RIGHT = 3;
79
80   static const int VALIGN_TOP = 1 << 4;
81   static const int VALIGN_CENTER = 2 << 4;
82   static const int VALIGN_BOTTOM = 3 << 4;
83
84   MapData(int priority) :
85     _dirtyText(true),
86     _age(0),
87     _priority(priority),
88     _width(0),
89     _height(0),
90     _offsetDir(HALIGN_LEFT | VALIGN_CENTER),
91     _offsetPx(10),
92     _dataVisible(false)
93   {
94   }
95
96   void setLabel(const std::string& label)
97   {
98     if (label == _label) {
99       return; // common case, and saves invalidation
100     }
101
102     _label = label;
103     _dirtyText = true;
104   }
105
106   void setText(const std::string &text)
107   {
108     if (_rawText == text) {
109       return; // common case, and saves invalidation
110     }
111
112     _rawText = text;
113     _dirtyText = true;
114   }
115
116   void setDataVisible(bool vis) {
117     if (vis == _dataVisible) {
118       return;
119     }
120
121     if (_rawText.empty()) {
122       vis = false;
123     }
124
125     _dataVisible = vis;
126     _dirtyText = true;
127   }
128
129   static void setFont(puFont f)
130   {
131     _font = f;
132     _fontHeight = f.getStringHeight();
133     _fontDescender = f.getStringDescender();
134   }
135
136   static void setPalette(puColor* pal)
137   {
138     _palette = pal;
139   }
140
141   void setPriority(int pri)
142   {
143     _priority = pri;
144   }
145
146   int priority() const
147   { return _priority; }
148
149   void setAnchor(const SGVec2d& anchor)
150   {
151     _anchor = anchor;
152   }
153
154   void setOffset(int direction, int px)
155   {
156     if ((_offsetPx == px) && (_offsetDir == direction)) {
157       return;
158     }
159
160     _dirtyOffset = true;
161     _offsetDir = direction;
162     _offsetPx = px;
163   }
164
165   bool isClipped(const puBox& vis) const
166   {
167     validate();
168     if ((_width < 1) || (_height < 1)) {
169       return true;
170     }
171
172     return !puBoxIntersect(vis, box());
173   }
174
175   bool overlaps(const MapDataVec& l) const
176   {
177     validate();
178     puBox b(box());
179
180     MapDataVec::const_iterator it;
181     for (it = l.begin(); it != l.end(); ++it) {
182       if (puBoxIntersect(b, (*it)->box())) {
183         return true;
184       }
185     } // of list iteration
186
187     return false;
188   }
189
190   puBox box() const
191   {
192     validate();
193     return makePuBox(
194       _anchor.x() + _offset.x(),
195       _anchor.y() + _offset.y(),
196       _width, _height);
197   }
198
199   void draw()
200   {
201     validate();
202
203     int xx = _anchor.x() + _offset.x();
204     int yy = _anchor.y() + _offset.y();
205
206     if (_dataVisible) {
207       puBox box(makePuBox(0,0,_width, _height));
208       int border = 1;
209       box.draw(xx, yy, PUSTYLE_DROPSHADOW, _palette, FALSE, border);
210
211       // draw lines
212       int lineHeight = _fontHeight;
213       int xPos = xx + MARGIN;
214       int yPos = yy + _height - (lineHeight + MARGIN);
215       glColor3f(0.8, 0.8, 0.8);
216
217       for (unsigned int ln=0; ln<_lines.size(); ++ln) {
218         _font.drawString(_lines[ln].c_str(), xPos, yPos);
219         yPos -= lineHeight + LINE_LEADING;
220       }
221     } else {
222       glColor3f(0.8, 0.8, 0.8);
223       _font.drawString(_label.c_str(), xx, yy + _fontDescender);
224     }
225   }
226
227   void age()
228   {
229     ++_age;
230   }
231
232   void resetAge()
233   {
234     _age = 0;
235   }
236
237   bool isExpired() const
238   { return (_age > 100); }
239
240   static bool order(MapData* a, MapData* b)
241   {
242     return a->_priority > b->_priority;
243   }
244 private:
245   void validate() const
246   {
247     if (!_dirtyText) {
248       if (_dirtyOffset) {
249         computeOffset();
250       }
251
252       return;
253     }
254
255     if (_dataVisible) {
256       measureData();
257     } else {
258       measureLabel();
259     }
260
261     computeOffset();
262     _dirtyText = false;
263   }
264
265   void measureData() const
266   {
267     _lines = simgear::strutils::split(_rawText, "\n");
268   // measure text to find width and height
269     _width = -1;
270     _height = 0;
271
272     for (unsigned int ln=0; ln<_lines.size(); ++ln) {
273       _height += _fontHeight;
274       if (ln > 0) {
275         _height += LINE_LEADING;
276       }
277
278       int lw = _font.getStringWidth(_lines[ln].c_str());
279       _width = std::max(_width, lw);
280     } // of line measurement
281
282     if ((_width < 1) || (_height < 1)) {
283       // will be clipped
284       return;
285     }
286
287     _height += MARGIN * 2;
288     _width += MARGIN * 2;
289   }
290
291   void measureLabel() const
292   {
293     if (_label.empty()) {
294       _width = _height = -1;
295       return;
296     }
297
298     _height = _fontHeight;
299     _width = _font.getStringWidth(_label.c_str());
300   }
301
302   void computeOffset() const
303   {
304     _dirtyOffset = false;
305     if ((_width <= 0) || (_height <= 0)) {
306       return;
307     }
308
309     int hOffset = 0;
310     int vOffset = 0;
311
312     switch (_offsetDir & 0x0f) {
313     default:
314     case HALIGN_LEFT:
315       hOffset = _offsetPx;
316       break;
317
318     case HALIGN_CENTER:
319       hOffset = -(_width>>1);
320       break;
321
322     case HALIGN_RIGHT:
323       hOffset = -(_offsetPx + _width);
324       break;
325     }
326
327     switch (_offsetDir & 0xf0) {
328     default:
329     case VALIGN_TOP:
330       vOffset = -(_offsetPx + _height);
331       break;
332
333     case VALIGN_CENTER:
334       vOffset = -(_height>>1);
335       break;
336
337     case VALIGN_BOTTOM:
338       vOffset = _offsetPx;
339       break;
340     }
341
342     _offset = SGVec2d(hOffset, vOffset);
343   }
344
345   static const int LINE_LEADING = 3;
346         static const int MARGIN = 3;
347
348   mutable bool _dirtyText;
349   mutable bool _dirtyOffset;
350   int _age;
351   std::string _rawText;
352   std::string _label;
353   mutable std::vector<std::string> _lines;
354   int _priority;
355   mutable int _width, _height;
356   SGVec2d _anchor;
357   int _offsetDir;
358   int _offsetPx;
359   mutable SGVec2d _offset;
360   bool _dataVisible;
361
362   static puFont _font;
363   static puColor* _palette;
364   static int _fontHeight;
365   static int _fontDescender;
366 };
367
368 puFont MapData::_font;
369 puColor* MapData::_palette;
370 int MapData::_fontHeight = 0;
371 int MapData::_fontDescender = 0;
372
373 ///////////////////////////////////////////////////////////////////////////
374
375 const int MAX_ZOOM = 12;
376 const int SHOW_DETAIL_ZOOM = 8;
377 const int CURSOR_PAN_STEP = 32;
378
379 MapWidget::MapWidget(int x, int y, int maxX, int maxY) :
380   puObject(x,y,maxX, maxY)
381 {
382   _route = static_cast<FGRouteMgr*>(globals->get_subsystem("route-manager"));
383   _gps = fgGetNode("/instrumentation/gps");
384
385   _width = maxX - x;
386   _height = maxY - y;
387   _hasPanned = false;
388   _orthoAzimuthProject = false;
389   
390   MapData::setFont(legendFont);
391   MapData::setPalette(colour);
392
393   _magVar = new SGMagVar();
394 }
395
396 MapWidget::~MapWidget()
397 {
398   delete _magVar;
399   clearData();
400 }
401
402 void MapWidget::setProperty(SGPropertyNode_ptr prop)
403 {
404   _root = prop;
405   int zoom = _root->getIntValue("zoom", -1);
406   if (zoom < 0) {
407     _root->setIntValue("zoom", 6); // default zoom
408   }
409   
410 // expose MAX_ZOOM to the UI
411   _root->setIntValue("max-zoom", MAX_ZOOM);
412   _root->setBoolValue("centre-on-aircraft", true);
413   _root->setBoolValue("draw-data", false);
414   _root->setBoolValue("magnetic-headings", true);
415 }
416
417 void MapWidget::setSize(int w, int h)
418 {
419   puObject::setSize(w, h);
420
421   _width = w;
422   _height = h;
423
424 }
425
426 void MapWidget::doHit( int button, int updown, int x, int y )
427 {
428   puObject::doHit(button, updown, x, y);
429   if (updown == PU_DRAG) {
430     handlePan(x, y);
431     return;
432   }
433
434   if (button == 3) { // mouse-wheel up
435     zoomIn();
436   } else if (button == 4) { // mouse-wheel down
437     zoomOut();
438   }
439
440   if (button != active_mouse_button) {
441     return;
442   }
443
444   _hitLocation = SGVec2d(x - abox.min[0], y - abox.min[1]);
445
446   if (updown == PU_UP) {
447     puDeactivateWidget();
448   } else if (updown == PU_DOWN) {
449     puSetActiveWidget(this, x, y);
450
451     if (fgGetKeyModifiers() & KEYMOD_CTRL) {
452       _clickGeod = unproject(_hitLocation - SGVec2d(_width>>1, _height>>1));
453     }
454   }
455 }
456
457 void MapWidget::handlePan(int x, int y)
458 {
459   SGVec2d delta = SGVec2d(x, y) - _hitLocation;
460   pan(delta);
461   _hitLocation = SGVec2d(x,y);
462 }
463
464 int MapWidget::checkKey (int key, int updown )
465 {
466   if ((updown == PU_UP) || !isVisible () || !isActive () || (window != puGetWindow())) {
467     return FALSE ;
468   }
469
470   switch (key)
471   {
472
473   case PU_KEY_UP:
474     pan(SGVec2d(0, -CURSOR_PAN_STEP));
475     break;
476
477   case PU_KEY_DOWN:
478     pan(SGVec2d(0, CURSOR_PAN_STEP));
479     break ;
480
481   case PU_KEY_LEFT:
482     pan(SGVec2d(CURSOR_PAN_STEP, 0));
483     break;
484
485   case PU_KEY_RIGHT:
486     pan(SGVec2d(-CURSOR_PAN_STEP, 0));
487     break;
488
489   case '-':
490     zoomOut();
491
492     break;
493
494   case '=':
495     zoomIn();
496     break;
497
498   default :
499     return FALSE;
500   }
501
502   return TRUE ;
503 }
504
505 void MapWidget::pan(const SGVec2d& delta)
506 {
507   _hasPanned = true; 
508   _projectionCenter = unproject(-delta);
509 }
510
511 int MapWidget::zoom() const
512 {
513   int z = _root->getIntValue("zoom");
514   SG_CLAMP_RANGE(z, 0, MAX_ZOOM);
515   return z;
516 }
517
518 void MapWidget::zoomIn()
519 {
520   if (zoom() >= MAX_ZOOM) {
521     return;
522   }
523
524   _root->setIntValue("zoom", zoom() + 1);
525 }
526
527 void MapWidget::zoomOut()
528 {
529   if (zoom() <= 0) {
530     return;
531   }
532
533   _root->setIntValue("zoom", zoom() - 1);
534 }
535
536 void MapWidget::draw(int dx, int dy)
537 {
538   _aircraft = SGGeod::fromDeg(fgGetDouble("/position/longitude-deg"),
539     fgGetDouble("/position/latitude-deg"));
540     
541   bool mag = _root->getBoolValue("magnetic-headings");
542   if (mag != _magneticHeadings) {
543     clearData(); // flush cached data text, since it often includes heading
544     _magneticHeadings =  mag;
545   }
546   
547   if (_hasPanned) {
548       _root->setBoolValue("centre-on-aircraft", false);
549       _hasPanned = false;
550   }
551   else if (_root->getBoolValue("centre-on-aircraft")) {
552     _projectionCenter = _aircraft;
553   }
554
555   double julianDate = globals->get_time_params()->getJD();
556   _magVar->update(_projectionCenter, julianDate);
557
558   bool aircraftUp = _root->getBoolValue("aircraft-heading-up");
559   if (aircraftUp) {
560     _upHeading = fgGetDouble("/orientation/heading-deg");
561   } else {
562     _upHeading = 0.0;
563   }
564
565   _cachedZoom = MAX_ZOOM - zoom();
566   SGGeod topLeft = unproject(SGVec2d(_width/2, _height/2));
567   // compute draw range, including a fudge factor for ILSs and other 'long'
568   // symbols
569   _drawRangeNm = SGGeodesy::distanceNm(_projectionCenter, topLeft) + 10.0;
570
571 // drawing operations
572   GLint sx = (int) abox.min[0],
573     sy = (int) abox.min[1];
574   glScissor(dx + sx, dy + sy, _width, _height);
575   glEnable(GL_SCISSOR_TEST);
576
577   glMatrixMode(GL_MODELVIEW);
578   glPushMatrix();
579   // cetere drawing about the widget center (which is also the
580   // projection centre)
581   glTranslated(dx + sx + (_width/2), dy + sy + (_height/2), 0.0);
582
583   drawLatLonGrid();
584
585   if (aircraftUp) {
586     int textHeight = legendFont.getStringHeight() + 5;
587
588     // draw heading line
589     SGVec2d loc = project(_aircraft);
590     glColor3f(1.0, 1.0, 1.0);
591     drawLine(loc, SGVec2d(loc.x(), (_height / 2) - textHeight));
592
593     int displayHdg;
594     if (_magneticHeadings) {
595       displayHdg = (int) fgGetDouble("/orientation/heading-magnetic-deg");
596     } else {
597       displayHdg = (int) _upHeading;
598     }
599
600     double y = (_height / 2) - textHeight;
601     char buf[16];
602     ::snprintf(buf, 16, "%d", displayHdg);
603     int sw = legendFont.getStringWidth(buf);
604     legendFont.drawString(buf, loc.x() - sw/2, y);
605   }
606
607   drawAirports();
608   drawNavaids();
609   drawTraffic();
610   drawGPSData();
611   drawNavRadio(fgGetNode("/instrumentation/nav[0]", false));
612   drawNavRadio(fgGetNode("/instrumentation/nav[1]", false));
613   paintAircraftLocation(_aircraft);
614   paintRoute();
615   paintRuler();
616
617   drawData();
618
619   glPopMatrix();
620   glDisable(GL_SCISSOR_TEST);
621 }
622
623 void MapWidget::paintRuler()
624 {
625   if (_clickGeod == SGGeod()) {
626     return;
627   }
628
629   SGVec2d acftPos = project(_aircraft);
630   SGVec2d clickPos = project(_clickGeod);
631
632   glColor4f(0.0, 1.0, 1.0, 0.6);
633   drawLine(acftPos, clickPos);
634
635   circleAtAlt(clickPos, 8, 10, 5);
636
637   double dist, az, az2;
638   SGGeodesy::inverse(_aircraft, _clickGeod, az, az2, dist);
639   char buffer[1024];
640         ::snprintf(buffer, 1024, "%03d/%.1fnm",
641                 displayHeading(az), dist * SG_METER_TO_NM);
642
643   MapData* d = getOrCreateDataForKey((void*) RULER_LEGEND_KEY);
644   d->setLabel(buffer);
645   d->setAnchor(clickPos);
646   d->setOffset(MapData::VALIGN_TOP | MapData::HALIGN_CENTER, 15);
647   d->setPriority(20000);
648
649
650 }
651
652 void MapWidget::paintAircraftLocation(const SGGeod& aircraftPos)
653 {
654   SGVec2d loc = project(aircraftPos);
655
656   double hdg = fgGetDouble("/orientation/heading-deg");
657
658   glLineWidth(2.0);
659   glColor4f(1.0, 1.0, 0.0, 1.0);
660   glPushMatrix();
661   glTranslated(loc.x(), loc.y(), 0.0);
662   glRotatef(hdg - _upHeading, 0.0, 0.0, -1.0);
663
664   const SGVec2d wingspan(12, 0);
665   const SGVec2d nose(0, 8);
666   const SGVec2d tail(0, -14);
667   const SGVec2d tailspan(4,0);
668
669   drawLine(-wingspan, wingspan);
670   drawLine(nose, tail);
671   drawLine(tail - tailspan, tail + tailspan);
672
673   glPopMatrix();
674   glLineWidth(1.0);
675 }
676
677 void MapWidget::paintRoute()
678 {
679   if (_route->numWaypts() < 2) {
680     return;
681   }
682
683   RoutePath path(_route->flightPlan());
684
685 // first pass, draw the actual lines
686   glLineWidth(2.0);
687
688   for (int w=0; w<_route->numWaypts(); ++w) {
689     SGGeodVec gv(path.pathForIndex(w));
690     if (gv.empty()) {
691       continue;
692     }
693
694     if (w < _route->currentIndex()) {
695       glColor4f(0.5, 0.5, 0.5, 0.7);
696     } else {
697       glColor4f(1.0, 0.0, 1.0, 1.0);
698     }
699
700     flightgear::WayptRef wpt(_route->wayptAtIndex(w));
701     if (wpt->flag(flightgear::WPT_MISS)) {
702       glEnable(GL_LINE_STIPPLE);
703       glLineStipple(1, 0x00FF);
704     }
705
706     glBegin(GL_LINE_STRIP);
707     for (unsigned int i=0; i<gv.size(); ++i) {
708       SGVec2d p = project(gv[i]);
709       glVertex2d(p.x(), p.y());
710     }
711
712     glEnd();
713     glDisable(GL_LINE_STIPPLE);
714   }
715
716   glLineWidth(1.0);
717 // second pass, draw waypoint symbols and data
718   for (int w=0; w < _route->numWaypts(); ++w) {
719     flightgear::WayptRef wpt(_route->wayptAtIndex(w));
720     SGGeod g = path.positionForIndex(w);
721     if (g == SGGeod()) {
722       continue; // Vectors or similar
723     }
724
725     SGVec2d p = project(g);
726     glColor4f(1.0, 0.0, 1.0, 1.0);
727     circleAtAlt(p, 8, 12, 5);
728
729     std::ostringstream legend;
730     legend << wpt->ident();
731     if (wpt->altitudeRestriction() != flightgear::RESTRICT_NONE) {
732       legend << '\n' << SGMiscd::roundToInt(wpt->altitudeFt()) << '\'';
733     }
734
735     if (wpt->speedRestriction() == flightgear::SPEED_RESTRICT_MACH) {
736       legend << '\n' << wpt->speedMach() << "M";
737     } else if (wpt->speedRestriction() != flightgear::RESTRICT_NONE) {
738       legend << '\n' << SGMiscd::roundToInt(wpt->speedKts()) << "Kts";
739     }
740
741     MapData* d = getOrCreateDataForKey(reinterpret_cast<void*>(w * 2));
742     d->setText(legend.str());
743     d->setLabel(wpt->ident());
744     d->setAnchor(p);
745     d->setOffset(MapData::VALIGN_TOP | MapData::HALIGN_CENTER, 15);
746     d->setPriority(w < _route->currentIndex() ? 9000 : 12000);
747
748   } // of second waypoint iteration
749 }
750
751 /**
752  * Round a SGGeod to an arbitrary precision.
753  * For example, passing precision of 0.5 will round to the nearest 0.5 of
754  * a degree in both lat and lon - passing in 3.0 rounds to the nearest 3 degree
755  * multiple, and so on.
756  */
757 static SGGeod roundGeod(double precision, const SGGeod& g)
758 {
759   double lon = SGMiscd::round(g.getLongitudeDeg() / precision);
760   double lat = SGMiscd::round(g.getLatitudeDeg() / precision);
761
762   return SGGeod::fromDeg(lon * precision, lat * precision);
763 }
764
765 bool MapWidget::drawLineClipped(const SGVec2d& a, const SGVec2d& b)
766 {
767   double minX = SGMiscd::min(a.x(), b.x()),
768     minY = SGMiscd::min(a.y(), b.y()),
769     maxX = SGMiscd::max(a.x(), b.x()),
770     maxY = SGMiscd::max(a.y(), b.y());
771
772   int hh = _height >> 1, hw = _width >> 1;
773
774   if ((maxX < -hw) || (minX > hw) || (minY > hh) || (maxY < -hh)) {
775     return false;
776   }
777
778   glVertex2dv(a.data());
779   glVertex2dv(b.data());
780   return true;
781 }
782
783 SGVec2d MapWidget::gridPoint(int ix, int iy)
784 {
785         int key = (ix + 0x7fff) | ((iy + 0x7fff) << 16);
786         GridPointCache::iterator it = _gridCache.find(key);
787         if (it != _gridCache.end()) {
788                 return it->second;
789         }
790
791         SGGeod gp = SGGeod::fromDeg(
792     _gridCenter.getLongitudeDeg() + ix * _gridSpacing,
793                 _gridCenter.getLatitudeDeg() + iy * _gridSpacing);
794
795         SGVec2d proj = project(gp);
796         _gridCache[key] = proj;
797         return proj;
798 }
799
800 void MapWidget::drawLatLonGrid()
801 {
802   _gridSpacing = 1.0;
803   _gridCenter = roundGeod(_gridSpacing, _projectionCenter);
804   _gridCache.clear();
805
806   int ix = 0;
807   int iy = 0;
808
809   glColor4f(0.8, 0.8, 0.8, 0.4);
810   glBegin(GL_LINES);
811   bool didDraw;
812   do {
813     didDraw = false;
814     ++ix;
815     ++iy;
816
817     for (int x = -ix; x < ix; ++x) {
818       didDraw |= drawLineClipped(gridPoint(x, -iy), gridPoint(x+1, -iy));
819       didDraw |= drawLineClipped(gridPoint(x, iy), gridPoint(x+1, iy));
820       didDraw |= drawLineClipped(gridPoint(x, -iy), gridPoint(x, -iy + 1));
821       didDraw |= drawLineClipped(gridPoint(x, iy), gridPoint(x, iy - 1));
822
823     }
824
825     for (int y = -iy; y < iy; ++y) {
826       didDraw |= drawLineClipped(gridPoint(-ix, y), gridPoint(-ix, y+1));
827       didDraw |= drawLineClipped(gridPoint(-ix, y), gridPoint(-ix + 1, y));
828       didDraw |= drawLineClipped(gridPoint(ix, y), gridPoint(ix, y+1));
829       didDraw |= drawLineClipped(gridPoint(ix, y), gridPoint(ix - 1, y));
830     }
831
832     if (ix > 30) {
833       break;
834     }
835   } while (didDraw);
836
837   glEnd();
838 }
839
840 void MapWidget::drawGPSData()
841 {
842   std::string gpsMode = _gps->getStringValue("mode");
843
844   SGGeod wp0Geod = SGGeod::fromDeg(
845         _gps->getDoubleValue("wp/wp[0]/longitude-deg"),
846         _gps->getDoubleValue("wp/wp[0]/latitude-deg"));
847
848   SGGeod wp1Geod = SGGeod::fromDeg(
849         _gps->getDoubleValue("wp/wp[1]/longitude-deg"),
850         _gps->getDoubleValue("wp/wp[1]/latitude-deg"));
851
852 // draw track line
853   double gpsTrackDeg = _gps->getDoubleValue("indicated-track-true-deg");
854   double gpsSpeed = _gps->getDoubleValue("indicated-ground-speed-kt");
855   double az2;
856
857   if (gpsSpeed > 3.0) { // only draw track line if valid
858     SGGeod trackRadial;
859     SGGeodesy::direct(_aircraft, gpsTrackDeg, _drawRangeNm * SG_NM_TO_METER, trackRadial, az2);
860
861     glColor4f(1.0, 1.0, 0.0, 1.0);
862     glEnable(GL_LINE_STIPPLE);
863     glLineStipple(1, 0x00FF);
864     drawLine(project(_aircraft), project(trackRadial));
865     glDisable(GL_LINE_STIPPLE);
866   }
867
868   if (gpsMode == "dto") {
869     SGVec2d wp0Pos = project(wp0Geod);
870     SGVec2d wp1Pos = project(wp1Geod);
871
872     glColor4f(1.0, 0.0, 1.0, 1.0);
873     drawLine(wp0Pos, wp1Pos);
874
875   }
876
877   if (_gps->getBoolValue("scratch/valid")) {
878     // draw scratch data
879
880   }
881 }
882
883 class MapAirportFilter : public FGAirport::AirportFilter
884 {
885 public:
886   MapAirportFilter(SGPropertyNode_ptr nd)
887   {
888     _heliports = nd->getBoolValue("show-heliports", false);
889     _hardRunwaysOnly = nd->getBoolValue("hard-surfaced-airports", true);
890     _minLengthFt = fgGetDouble("/sim/navdb/min-runway-length-ft", 2000);
891   }
892
893   virtual FGPositioned::Type maxType() const {
894     return _heliports ? FGPositioned::HELIPORT : FGPositioned::AIRPORT;
895   }
896
897   virtual bool passAirport(FGAirport* aApt) const {
898     if (_hardRunwaysOnly) {
899       return aApt->hasHardRunwayOfLengthFt(_minLengthFt);
900     }
901
902     return true;
903   }
904
905 private:
906   bool _heliports;
907   bool _hardRunwaysOnly;
908   double _minLengthFt;
909 };
910
911 void MapWidget::drawAirports()
912 {
913   MapAirportFilter af(_root);
914   FGPositioned::List apts = FGPositioned::findWithinRange(_projectionCenter, _drawRangeNm, &af);
915   for (unsigned int i=0; i<apts.size(); ++i) {
916     drawAirport((FGAirport*) apts[i].get());
917   }
918 }
919
920 class NavaidFilter : public FGPositioned::Filter
921 {
922 public:
923   NavaidFilter(bool fixesEnabled, bool navaidsEnabled) :
924     _fixes(fixesEnabled),
925     _navaids(navaidsEnabled)
926   {}
927
928   virtual bool pass(FGPositioned* aPos) const {
929     if (_fixes && (aPos->type() == FGPositioned::FIX)) {
930       // ignore fixes which end in digits - expirmental
931       if (aPos->ident().length() > 4 && isdigit(aPos->ident()[3]) && isdigit(aPos->ident()[4])) {
932         return false;
933       }
934     }
935
936     return true;
937   }
938
939   virtual FGPositioned::Type minType() const {
940     return _fixes ? FGPositioned::FIX : FGPositioned::NDB;
941   }
942
943   virtual FGPositioned::Type maxType() const {
944     return _navaids ? FGPositioned::VOR : FGPositioned::FIX;
945   }
946
947 private:
948   bool _fixes, _navaids;
949 };
950
951 void MapWidget::drawNavaids()
952 {
953   bool fixes = _root->getBoolValue("draw-fixes");
954   NavaidFilter f(fixes, _root->getBoolValue("draw-navaids"));
955
956   if (f.minType() <= f.maxType()) {
957     FGPositioned::List navs = FGPositioned::findWithinRange(_projectionCenter, _drawRangeNm, &f);
958
959     glLineWidth(1.0);
960     for (unsigned int i=0; i<navs.size(); ++i) {
961       FGPositioned::Type ty = navs[i]->type();
962       if (ty == FGPositioned::NDB) {
963         drawNDB(false, (FGNavRecord*) navs[i].get());
964       } else if (ty == FGPositioned::VOR) {
965         drawVOR(false, (FGNavRecord*) navs[i].get());
966       } else if (ty == FGPositioned::FIX) {
967         drawFix((FGFix*) navs[i].get());
968       }
969     } // of navaid iteration
970   } // of navaids || fixes are drawn test
971 }
972
973 void MapWidget::drawNDB(bool tuned, FGNavRecord* ndb)
974 {
975   SGVec2d pos = project(ndb->geod());
976
977   if (tuned) {
978     glColor3f(0.0, 1.0, 1.0);
979   } else {
980     glColor3f(0.0, 0.0, 0.0);
981   }
982
983   glEnable(GL_LINE_STIPPLE);
984   glLineStipple(1, 0x00FF);
985   circleAt(pos, 20, 6);
986   circleAt(pos, 20, 10);
987   glDisable(GL_LINE_STIPPLE);
988
989   if (validDataForKey(ndb)) {
990     setAnchorForKey(ndb, pos);
991     return;
992   }
993
994   char buffer[1024];
995         ::snprintf(buffer, 1024, "%s\n%s %3.0fKhz",
996                 ndb->name().c_str(), ndb->ident().c_str(),ndb->get_freq()/100.0);
997
998   MapData* d = createDataForKey(ndb);
999   d->setPriority(40);
1000   d->setLabel(ndb->ident());
1001   d->setText(buffer);
1002   d->setOffset(MapData::HALIGN_CENTER | MapData::VALIGN_BOTTOM, 10);
1003   d->setAnchor(pos);
1004
1005 }
1006
1007 void MapWidget::drawVOR(bool tuned, FGNavRecord* vor)
1008 {
1009   SGVec2d pos = project(vor->geod());
1010   if (tuned) {
1011     glColor3f(0.0, 1.0, 1.0);
1012   } else {
1013     glColor3f(0.0, 0.0, 1.0);
1014   }
1015
1016   circleAt(pos, 6, 8);
1017
1018   if (validDataForKey(vor)) {
1019     setAnchorForKey(vor, pos);
1020     return;
1021   }
1022
1023   char buffer[1024];
1024         ::snprintf(buffer, 1024, "%s\n%s %6.3fMhz",
1025                 vor->name().c_str(), vor->ident().c_str(),
1026     vor->get_freq() / 100.0);
1027
1028   MapData* d = createDataForKey(vor);
1029   d->setText(buffer);
1030   d->setLabel(vor->ident());
1031   d->setPriority(tuned ? 10000 : 100);
1032   d->setOffset(MapData::HALIGN_CENTER | MapData::VALIGN_BOTTOM, 12);
1033   d->setAnchor(pos);
1034 }
1035
1036 void MapWidget::drawFix(FGFix* fix)
1037 {
1038   SGVec2d pos = project(fix->geod());
1039   glColor3f(0.0, 0.0, 0.0);
1040   circleAt(pos, 3, 6);
1041
1042   if (_cachedZoom > SHOW_DETAIL_ZOOM) {
1043     return; // hide fix labels beyond a certain zoom level
1044   }
1045
1046   if (validDataForKey(fix)) {
1047     setAnchorForKey(fix, pos);
1048     return;
1049   }
1050
1051   MapData* d = createDataForKey(fix);
1052   d->setLabel(fix->ident());
1053   d->setPriority(20);
1054   d->setOffset(MapData::VALIGN_CENTER | MapData::HALIGN_LEFT, 10);
1055   d->setAnchor(pos);
1056 }
1057
1058 void MapWidget::drawNavRadio(SGPropertyNode_ptr radio)
1059 {
1060   if (!radio || radio->getBoolValue("slaved-to-gps", false)
1061         || !radio->getBoolValue("in-range", false)) {
1062     return;
1063   }
1064
1065   if (radio->getBoolValue("nav-loc", false)) {
1066     drawTunedLocalizer(radio);
1067   }
1068
1069   // identify the tuned station - unfortunately we don't get lat/lon directly,
1070   // need to do the frequency search again
1071   double mhz = radio->getDoubleValue("frequencies/selected-mhz", 0.0);
1072
1073   FGNavRecord* nav = FGNavList::findByFreq(mhz, _aircraft,
1074                                            FGNavList::navFilter());
1075   if (!nav || (nav->ident() != radio->getStringValue("nav-id"))) {
1076     // mismatch between navradio selection logic and ours!
1077     return;
1078   }
1079
1080   glLineWidth(1.0);
1081   drawVOR(true, nav);
1082
1083   SGVec2d pos = project(nav->geod());
1084   SGGeod range;
1085   double az2;
1086   double trueRadial = radio->getDoubleValue("radials/target-radial-deg");
1087   SGGeodesy::direct(nav->geod(), trueRadial, nav->get_range() * SG_NM_TO_METER, range, az2);
1088   SGVec2d prange = project(range);
1089
1090   SGVec2d norm = normalize(prange - pos);
1091   SGVec2d perp(norm.y(), -norm.x());
1092
1093   circleAt(pos, 64, length(prange - pos));
1094   drawLine(pos, prange);
1095
1096 // draw to/from arrows
1097   SGVec2d midPoint = (pos + prange) * 0.5;
1098   if (radio->getBoolValue("from-flag")) {
1099     norm = -norm;
1100     perp = -perp;
1101   }
1102
1103   int sz = 10;
1104   SGVec2d arrowB = midPoint - (norm * sz) + (perp * sz);
1105   SGVec2d arrowC = midPoint - (norm * sz) - (perp * sz);
1106   drawLine(midPoint, arrowB);
1107   drawLine(arrowB, arrowC);
1108   drawLine(arrowC, midPoint);
1109
1110   drawLine(pos, (2 * pos) - prange); // reciprocal radial
1111 }
1112
1113 void MapWidget::drawTunedLocalizer(SGPropertyNode_ptr radio)
1114 {
1115   double mhz = radio->getDoubleValue("frequencies/selected-mhz", 0.0);
1116   FGNavRecord* loc = FGNavList::findByFreq(mhz, _aircraft, FGNavList::locFilter());
1117   if (!loc || (loc->ident() != radio->getStringValue("nav-id"))) {
1118     // mismatch between navradio selection logic and ours!
1119     return;
1120   }
1121
1122   if (loc->runway()) {
1123     drawILS(true, loc->runway());
1124   }
1125 }
1126
1127 /*
1128 void MapWidget::drawObstacle(FGPositioned* obs)
1129 {
1130   SGVec2d pos = project(obs->geod());
1131   glColor3f(0.0, 0.0, 0.0);
1132   glLineWidth(2.0);
1133   drawLine(pos, pos + SGVec2d());
1134 }
1135 */
1136
1137 void MapWidget::drawAirport(FGAirport* apt)
1138 {
1139         // draw tower location
1140         SGVec2d towerPos = project(apt->getTowerLocation());
1141
1142   if (_cachedZoom <= SHOW_DETAIL_ZOOM) {
1143     glColor3f(1.0, 1.0, 1.0);
1144     glLineWidth(1.0);
1145
1146     drawLine(towerPos + SGVec2d(3, 0), towerPos + SGVec2d(3, 10));
1147     drawLine(towerPos + SGVec2d(-3, 0), towerPos + SGVec2d(-3, 10));
1148     drawLine(towerPos + SGVec2d(-6, 20), towerPos + SGVec2d(-3, 10));
1149     drawLine(towerPos + SGVec2d(6, 20), towerPos + SGVec2d(3, 10));
1150     drawLine(towerPos + SGVec2d(-6, 20), towerPos + SGVec2d(6, 20));
1151   }
1152
1153   if (validDataForKey(apt)) {
1154     setAnchorForKey(apt, towerPos);
1155   } else {
1156     char buffer[1024];
1157     ::snprintf(buffer, 1024, "%s\n%s",
1158       apt->ident().c_str(), apt->name().c_str());
1159
1160     MapData* d = createDataForKey(apt);
1161     d->setText(buffer);
1162     d->setLabel(apt->ident());
1163     d->setPriority(100 + scoreAirportRunways(apt));
1164     d->setOffset(MapData::VALIGN_TOP | MapData::HALIGN_CENTER, 6);
1165     d->setAnchor(towerPos);
1166   }
1167
1168   if (_cachedZoom > SHOW_DETAIL_ZOOM) {
1169     return;
1170   }
1171
1172   for (unsigned int r=0; r<apt->numRunways(); ++r) {
1173     FGRunway* rwy = apt->getRunwayByIndex(r);
1174                 if (!rwy->isReciprocal()) {
1175                         drawRunwayPre(rwy);
1176                 }
1177   }
1178
1179         for (unsigned int r=0; r<apt->numRunways(); ++r) {
1180                 FGRunway* rwy = apt->getRunwayByIndex(r);
1181                 if (!rwy->isReciprocal()) {
1182                         drawRunway(rwy);
1183                 }
1184
1185                 if (rwy->ILS()) {
1186                         drawILS(false, rwy);
1187                 }
1188         } // of runway iteration
1189
1190 }
1191
1192 int MapWidget::scoreAirportRunways(FGAirport* apt)
1193 {
1194   bool needHardSurface = _root->getBoolValue("hard-surfaced-airports", true);
1195   double minLength = _root->getDoubleValue("min-runway-length-ft", 2000.0);
1196
1197   int score = 0;
1198   unsigned int numRunways(apt->numRunways());
1199   for (unsigned int r=0; r<numRunways; ++r) {
1200     FGRunway* rwy = apt->getRunwayByIndex(r);
1201     if (rwy->isReciprocal()) {
1202       continue;
1203     }
1204
1205     if (needHardSurface && !rwy->isHardSurface()) {
1206       continue;
1207     }
1208
1209     if (rwy->lengthFt() < minLength) {
1210       continue;
1211     }
1212
1213     int scoreLength = SGMiscd::roundToInt(rwy->lengthFt() / 200.0);
1214     score += scoreLength;
1215   } // of runways iteration
1216
1217   return score;
1218 }
1219
1220 void MapWidget::drawRunwayPre(FGRunway* rwy)
1221 {
1222   SGVec2d p1 = project(rwy->begin());
1223         SGVec2d p2 = project(rwy->end());
1224
1225   glLineWidth(4.0);
1226   glColor3f(1.0, 0.0, 1.0);
1227         drawLine(p1, p2);
1228 }
1229
1230 void MapWidget::drawRunway(FGRunway* rwy)
1231 {
1232         // line for runway
1233         // optionally show active, stopway, etc
1234         // in legend, show published heading and length
1235         // and threshold elevation
1236
1237   SGVec2d p1 = project(rwy->begin());
1238         SGVec2d p2 = project(rwy->end());
1239   glLineWidth(2.0);
1240   glColor3f(1.0, 1.0, 1.0);
1241   SGVec2d inset = normalize(p2 - p1) * 2;
1242
1243         drawLine(p1 + inset, p2 - inset);
1244
1245   if (validDataForKey(rwy)) {
1246     setAnchorForKey(rwy, (p1 + p2) * 0.5);
1247     return;
1248   }
1249   
1250         char buffer[1024];
1251         ::snprintf(buffer, 1024, "%s/%s\n%03d/%03d\n%.0f'",
1252                 rwy->ident().c_str(),
1253                 rwy->reciprocalRunway()->ident().c_str(),
1254                 displayHeading(rwy->headingDeg()),
1255                 displayHeading(rwy->reciprocalRunway()->headingDeg()),
1256                 rwy->lengthFt());
1257
1258   MapData* d = createDataForKey(rwy);
1259   d->setText(buffer);
1260   d->setLabel(rwy->ident() + "/" + rwy->reciprocalRunway()->ident());
1261   d->setPriority(50);
1262   d->setOffset(MapData::HALIGN_CENTER | MapData::VALIGN_BOTTOM, 12);
1263   d->setAnchor((p1 + p2) * 0.5);
1264 }
1265
1266 void MapWidget::drawILS(bool tuned, FGRunway* rwy)
1267 {
1268         // arrow, tip centered on the landing threshold
1269   // using LOC transmitter position would be more accurate, but
1270   // is visually cluttered
1271         // arrow width is based upon the computed localizer width
1272
1273         FGNavRecord* loc = rwy->ILS();
1274         double halfBeamWidth = loc->localizerWidth() * 0.5;
1275         SGVec2d t = project(rwy->threshold());
1276         SGGeod locEnd;
1277         double rangeM = loc->get_range() * SG_NM_TO_METER;
1278         double radial = loc->get_multiuse();
1279   SG_NORMALIZE_RANGE(radial, 0.0, 360.0);
1280         double az2;
1281
1282 // compute the three end points at the widge end of the arrow
1283         SGGeodesy::direct(loc->geod(), radial, -rangeM, locEnd, az2);
1284         SGVec2d endCentre = project(locEnd);
1285
1286         SGGeodesy::direct(loc->geod(), radial + halfBeamWidth, -rangeM * 1.1, locEnd, az2);
1287         SGVec2d endR = project(locEnd);
1288
1289         SGGeodesy::direct(loc->geod(), radial - halfBeamWidth, -rangeM * 1.1, locEnd, az2);
1290         SGVec2d endL = project(locEnd);
1291
1292 // outline two triangles
1293   glLineWidth(1.0);
1294   if (tuned) {
1295     glColor3f(0.0, 1.0, 1.0);
1296   } else {
1297     glColor3f(0.0, 0.0, 1.0);
1298         }
1299
1300   glBegin(GL_LINE_LOOP);
1301                 glVertex2dv(t.data());
1302                 glVertex2dv(endCentre.data());
1303                 glVertex2dv(endL.data());
1304         glEnd();
1305         glBegin(GL_LINE_LOOP);
1306                 glVertex2dv(t.data());
1307                 glVertex2dv(endCentre.data());
1308                 glVertex2dv(endR.data());
1309         glEnd();
1310
1311         if (validDataForKey(loc)) {
1312     setAnchorForKey(loc, endR);
1313     return;
1314   }
1315
1316         char buffer[1024];
1317         ::snprintf(buffer, 1024, "%s\n%s\n%03d - %3.2fMHz",
1318                 loc->ident().c_str(), loc->name().c_str(),
1319     displayHeading(radial),
1320     loc->get_freq()/100.0);
1321
1322   MapData* d = createDataForKey(loc);
1323   d->setPriority(40);
1324   d->setLabel(loc->ident());
1325   d->setText(buffer);
1326   d->setOffset(MapData::HALIGN_CENTER | MapData::VALIGN_BOTTOM, 10);
1327   d->setAnchor(endR);
1328 }
1329
1330 void MapWidget::drawTraffic()
1331 {
1332   if (!_root->getBoolValue("draw-traffic")) {
1333     return;
1334   }
1335
1336   if (_cachedZoom > SHOW_DETAIL_ZOOM) {
1337     return;
1338   }
1339
1340   const SGPropertyNode* ai = fgGetNode("/ai/models", true);
1341
1342   for (int i = 0; i < ai->nChildren(); ++i) {
1343     const SGPropertyNode *model = ai->getChild(i);
1344     // skip bad or dead entries
1345     if (!model || model->getIntValue("id", -1) == -1) {
1346       continue;
1347     }
1348
1349     const std::string& name(model->getName());
1350     SGGeod pos = SGGeod::fromDegFt(
1351       model->getDoubleValue("position/longitude-deg"),
1352       model->getDoubleValue("position/latitude-deg"),
1353       model->getDoubleValue("position/altitude-ft"));
1354
1355     double dist = SGGeodesy::distanceNm(_projectionCenter, pos);
1356     if (dist > _drawRangeNm) {
1357       continue;
1358     }
1359
1360     double heading = model->getDoubleValue("orientation/true-heading-deg");
1361     if ((name == "aircraft") || (name == "multiplayer") ||
1362         (name == "wingman") || (name == "tanker")) {
1363       drawAIAircraft(model, pos, heading);
1364     } else if ((name == "ship") || (name == "carrier") || (name == "escort")) {
1365       drawAIShip(model, pos, heading);
1366     }
1367   } // of ai/models iteration
1368 }
1369
1370 void MapWidget::drawAIAircraft(const SGPropertyNode* model, const SGGeod& pos, double hdg)
1371 {
1372
1373   SGVec2d p = project(pos);
1374
1375   glColor3f(0.0, 0.0, 0.0);
1376   glLineWidth(2.0);
1377   circleAt(p, 4, 6.0); // black diamond
1378
1379 // draw heading vector
1380   int speedKts = static_cast<int>(model->getDoubleValue("velocities/true-airspeed-kt"));
1381   if (speedKts > 1) {
1382     glLineWidth(1.0);
1383
1384     const double dt = 15.0 / (3600.0); // 15 seconds look-ahead
1385     double distanceM = speedKts * SG_NM_TO_METER * dt;
1386
1387     SGGeod advance;
1388     double az2;
1389     SGGeodesy::direct(pos, hdg, distanceM, advance, az2);
1390
1391     drawLine(p, project(advance));
1392   }
1393
1394
1395   // draw callsign / altitude / speed
1396   char buffer[1024];
1397         ::snprintf(buffer, 1024, "%s\n%d'\n%dkts",
1398                 model->getStringValue("callsign", "<>"),
1399                 static_cast<int>(pos.getElevationFt() / 50.0) * 50,
1400     speedKts);
1401
1402   MapData* d = getOrCreateDataForKey((void*) model);
1403   d->setText(buffer);
1404   d->setLabel(model->getStringValue("callsign", "<>"));
1405   d->setPriority(speedKts > 5 ? 60 : 10); // low priority for parked aircraft
1406   d->setOffset(MapData::VALIGN_CENTER | MapData::HALIGN_LEFT, 10);
1407   d->setAnchor(p);
1408
1409 }
1410
1411 void MapWidget::drawAIShip(const SGPropertyNode* model, const SGGeod& pos, double hdg)
1412 {
1413   SGVec2d p = project(pos);
1414
1415   glColor3f(0.0, 0.0, 0.5);
1416   glLineWidth(2.0);
1417   circleAt(p, 4, 6.0); // blue diamond (to differentiate from aircraft.
1418
1419 // draw heading vector
1420   int speedKts = static_cast<int>(model->getDoubleValue("velocities/speed-kts"));
1421   if (speedKts > 1) {
1422     glLineWidth(1.0);
1423
1424     const double dt = 15.0 / (3600.0); // 15 seconds look-ahead
1425     double distanceM = speedKts * SG_NM_TO_METER * dt;
1426
1427     SGGeod advance;
1428     double az2;
1429     SGGeodesy::direct(pos, hdg, distanceM, advance, az2);
1430
1431     drawLine(p, project(advance));
1432   }
1433
1434   // draw callsign / speed
1435   char buffer[1024];
1436         ::snprintf(buffer, 1024, "%s\n%dkts",
1437                 model->getStringValue("name", "<>"),
1438     speedKts);
1439
1440   MapData* d = getOrCreateDataForKey((void*) model);
1441   d->setText(buffer);
1442   d->setLabel(model->getStringValue("name", "<>"));
1443   d->setPriority(speedKts > 2 ? 30 : 10); // low priority for slow moving ships
1444   d->setOffset(MapData::VALIGN_CENTER | MapData::HALIGN_LEFT, 10);
1445   d->setAnchor(p);
1446 }
1447
1448 SGVec2d MapWidget::project(const SGGeod& geod) const
1449 {
1450   SGVec2d p;
1451   double r = earth_radius_lat(geod.getLatitudeRad());
1452   
1453   if (_orthoAzimuthProject) {
1454     // http://mathworld.wolfram.com/OrthographicProjection.html
1455     double cosTheta = cos(geod.getLatitudeRad());
1456     double sinDLambda = sin(geod.getLongitudeRad() - _projectionCenter.getLongitudeRad());
1457     double cosDLambda = cos(geod.getLongitudeRad() - _projectionCenter.getLongitudeRad());
1458     double sinTheta1 = sin(_projectionCenter.getLatitudeRad());
1459     double sinTheta = sin(geod.getLatitudeRad());
1460     double cosTheta1 = cos(_projectionCenter.getLatitudeRad());
1461     
1462     p = SGVec2d(cosTheta * sinDLambda,
1463                 (cosTheta1 * sinTheta) - (sinTheta1 * cosTheta * cosDLambda)) * r * currentScale();
1464     
1465   } else {
1466     // Sanson-Flamsteed projection, relative to the projection center
1467     double lonDiff = geod.getLongitudeRad() - _projectionCenter.getLongitudeRad(),
1468       latDiff = geod.getLatitudeRad() - _projectionCenter.getLatitudeRad();
1469
1470     p = SGVec2d(cos(geod.getLatitudeRad()) * lonDiff, latDiff) * r * currentScale();
1471       
1472   }
1473   
1474 // rotate as necessary
1475   double cost = cos(_upHeading * SG_DEGREES_TO_RADIANS),
1476     sint = sin(_upHeading * SG_DEGREES_TO_RADIANS);
1477   double rx = cost * p.x() - sint * p.y();
1478   double ry = sint * p.x() + cost * p.y();
1479   return SGVec2d(rx, ry);
1480 }
1481
1482 SGGeod MapWidget::unproject(const SGVec2d& p) const
1483 {
1484   // unrotate, if necessary
1485   double cost = cos(-_upHeading * SG_DEGREES_TO_RADIANS),
1486     sint = sin(-_upHeading * SG_DEGREES_TO_RADIANS);
1487   SGVec2d ur(cost * p.x() - sint * p.y(),
1488              sint * p.x() + cost * p.y());
1489
1490   double r = earth_radius_lat(_projectionCenter.getLatitudeRad());
1491   SGVec2d unscaled = ur * (1.0 / (currentScale() * r));
1492
1493   if (_orthoAzimuthProject) {
1494       double phi = length(p);
1495       double c = asin(phi);
1496       double sinTheta1 = sin(_projectionCenter.getLatitudeRad());
1497       double cosTheta1 = cos(_projectionCenter.getLatitudeRad());
1498       
1499       double lat = asin(cos(c) * sinTheta1 + ((unscaled.y() * sin(c) * cosTheta1) / phi));
1500       double lon = _projectionCenter.getLongitudeRad() + 
1501         atan((unscaled.x()* sin(c)) / (phi  * cosTheta1 * cos(c) - unscaled.y() * sinTheta1 * sin(c)));
1502       return SGGeod::fromRad(lon, lat);
1503   } else {
1504       double lat = unscaled.y() + _projectionCenter.getLatitudeRad();
1505       double lon = (unscaled.x() / cos(lat)) + _projectionCenter.getLongitudeRad();
1506       return SGGeod::fromRad(lon, lat);
1507   }
1508 }
1509
1510 double MapWidget::currentScale() const
1511 {
1512   return 1.0 / pow(2.0, _cachedZoom);
1513 }
1514
1515 void MapWidget::circleAt(const SGVec2d& center, int nSides, double r)
1516 {
1517   glBegin(GL_LINE_LOOP);
1518   double advance = (SGD_PI * 2) / nSides;
1519   glVertex2d(center.x(), center.y() + r);
1520   double t=advance;
1521   for (int i=1; i<nSides; ++i) {
1522     glVertex2d(center.x() + (sin(t) * r), center.y() + (cos(t) * r));
1523     t += advance;
1524   }
1525   glEnd();
1526 }
1527
1528 void MapWidget::circleAtAlt(const SGVec2d& center, int nSides, double r, double r2)
1529 {
1530   glBegin(GL_LINE_LOOP);
1531   double advance = (SGD_PI * 2) / nSides;
1532   glVertex2d(center.x(), center.y() + r);
1533   double t=advance;
1534   for (int i=1; i<nSides; ++i) {
1535     double rr = (i%2 == 0) ? r : r2;
1536     glVertex2d(center.x() + (sin(t) * rr), center.y() + (cos(t) * rr));
1537     t += advance;
1538   }
1539   glEnd();
1540 }
1541
1542 void MapWidget::drawLine(const SGVec2d& p1, const SGVec2d& p2)
1543 {
1544   glBegin(GL_LINES);
1545     glVertex2dv(p1.data());
1546     glVertex2dv(p2.data());
1547   glEnd();
1548 }
1549
1550 void MapWidget::drawLegendBox(const SGVec2d& pos, const std::string& t)
1551 {
1552         std::vector<std::string> lines(simgear::strutils::split(t, "\n"));
1553         const int LINE_LEADING = 4;
1554         const int MARGIN = 4;
1555
1556 // measure
1557         int maxWidth = -1, totalHeight = 0;
1558         int lineHeight = legendFont.getStringHeight();
1559
1560         for (unsigned int ln=0; ln<lines.size(); ++ln) {
1561                 totalHeight += lineHeight;
1562                 if (ln > 0) {
1563                         totalHeight += LINE_LEADING;
1564                 }
1565
1566                 int lw = legendFont.getStringWidth(lines[ln].c_str());
1567                 maxWidth = std::max(maxWidth, lw);
1568         } // of line measurement
1569
1570         if (maxWidth < 0) {
1571                 return; // all lines are empty, don't draw
1572         }
1573
1574         totalHeight += MARGIN * 2;
1575
1576 // draw box
1577         puBox box;
1578         box.min[0] = 0;
1579         box.min[1] = -totalHeight;
1580         box.max[0] = maxWidth + (MARGIN * 2);
1581         box.max[1] = 0;
1582         int border = 1;
1583         box.draw (pos.x(), pos.y(), PUSTYLE_DROPSHADOW, colour, FALSE, border);
1584
1585 // draw lines
1586         int xPos = pos.x() + MARGIN;
1587         int yPos = pos.y() - (lineHeight + MARGIN);
1588         glColor3f(0.8, 0.8, 0.8);
1589
1590         for (unsigned int ln=0; ln<lines.size(); ++ln) {
1591                 legendFont.drawString(lines[ln].c_str(), xPos, yPos);
1592                 yPos -= lineHeight + LINE_LEADING;
1593         }
1594 }
1595
1596 void MapWidget::drawData()
1597 {
1598   std::sort(_dataQueue.begin(), _dataQueue.end(), MapData::order);
1599
1600   int hw = _width >> 1,
1601     hh = _height >> 1;
1602   puBox visBox(makePuBox(-hw, -hh, _width, _height));
1603
1604   unsigned int d = 0;
1605   int drawn = 0;
1606   std::vector<MapData*> drawQueue;
1607
1608   bool drawData = _root->getBoolValue("draw-data");
1609   const int MAX_DRAW_DATA = 25;
1610   const int MAX_DRAW = 50;
1611
1612   for (; (d < _dataQueue.size()) && (drawn < MAX_DRAW); ++d) {
1613     MapData* md = _dataQueue[d];
1614     md->setDataVisible(drawData);
1615
1616     if (md->isClipped(visBox)) {
1617       continue;
1618     }
1619
1620     if (md->overlaps(drawQueue)) {
1621       if (drawData) { // overlapped with data, let's try just the label
1622         md->setDataVisible(false);
1623         if (md->overlaps(drawQueue)) {
1624           continue;
1625         }
1626       } else {
1627         continue;
1628       }
1629     } // of overlaps case
1630
1631     drawQueue.push_back(md);
1632     ++drawn;
1633     if (drawData && (drawn >= MAX_DRAW_DATA)) {
1634       drawData = false;
1635     }
1636   }
1637
1638   // draw lowest-priority first, so higher-priorty items appear on top
1639   std::vector<MapData*>::reverse_iterator r;
1640   for (r = drawQueue.rbegin(); r!= drawQueue.rend(); ++r) {
1641     (*r)->draw();
1642   }
1643
1644   _dataQueue.clear();
1645   KeyDataMap::iterator it = _mapData.begin();
1646   for (; it != _mapData.end(); ) {
1647     it->second->age();
1648     if (it->second->isExpired()) {
1649       delete it->second;
1650       KeyDataMap::iterator cur = it++;
1651       _mapData.erase(cur);
1652     } else {
1653       ++it;
1654     }
1655   } // of expiry iteration
1656 }
1657
1658 bool MapWidget::validDataForKey(void* key)
1659 {
1660   KeyDataMap::iterator it = _mapData.find(key);
1661   if (it == _mapData.end()) {
1662     return false; // no valid data for the key!
1663   }
1664
1665   it->second->resetAge(); // mark data as valid this frame
1666   _dataQueue.push_back(it->second);
1667   return true;
1668 }
1669
1670 void MapWidget::setAnchorForKey(void* key, const SGVec2d& anchor)
1671 {
1672   KeyDataMap::iterator it = _mapData.find(key);
1673   if (it == _mapData.end()) {
1674     throw sg_exception("no valid data for key!");
1675   }
1676
1677   it->second->setAnchor(anchor);
1678 }
1679
1680 MapData* MapWidget::getOrCreateDataForKey(void* key)
1681 {
1682   KeyDataMap::iterator it = _mapData.find(key);
1683   if (it == _mapData.end()) {
1684     return createDataForKey(key);
1685   }
1686
1687   it->second->resetAge(); // mark data as valid this frame
1688   _dataQueue.push_back(it->second);
1689   return it->second;
1690 }
1691
1692 MapData* MapWidget::createDataForKey(void* key)
1693 {
1694   KeyDataMap::iterator it = _mapData.find(key);
1695   if (it != _mapData.end()) {
1696     throw sg_exception("duplicate data requested for key!");
1697   }
1698
1699   MapData* d =  new MapData(0);
1700   _mapData[key] = d;
1701   _dataQueue.push_back(d);
1702   d->resetAge();
1703   return d;
1704 }
1705
1706 void MapWidget::clearData()
1707 {
1708   KeyDataMap::iterator it = _mapData.begin();
1709   for (; it != _mapData.end(); ++it) {
1710     delete it->second;
1711   }
1712   
1713   _mapData.clear();
1714 }
1715
1716 int MapWidget::displayHeading(double h) const
1717 {
1718   if (_magneticHeadings) {
1719     h -= _magVar->get_magvar() * SG_RADIANS_TO_DEGREES;
1720   }
1721   
1722   SG_NORMALIZE_RANGE(h, 0.0, 360.0);
1723   return SGMiscd::roundToInt(h);
1724 }