Added support for Orux maps

gpxsee.pro

@@ -117,6 +117,7 @@ HEADERS += src/common/config.h \
     src/data/twonavparser.h \
     src/map/IMG/section.h \
     src/map/gemfmap.h \
+    src/map/oruxmap.h \
     src/map/osmdroidmap.h \
     src/map/proj/polyconic.h \
     src/map/proj/webmercator.h \
@@ -321,6 +322,7 @@ SOURCES += src/main.cpp \
     src/GUI/passwordedit.cpp \
     src/data/twonavparser.cpp \
     src/map/gemfmap.cpp \
+    src/map/oruxmap.cpp \
     src/map/osmdroidmap.cpp \
     src/map/proj/polyconic.cpp \
     src/map/proj/webmercator.cpp \

src/map/maplist.cpp

@@ -18,6 +18,7 @@
 #include "qctmap.h"
 #include "osmdroidmap.h"
 #include "gemfmap.h"
+#include "oruxmap.h"
 #include "invalidmap.h"
 #include "maplist.h"
 
@@ -52,6 +53,7 @@ MapList::ParserMap MapList::parsers()
 	map.insert("qct", &QCTMap::create);
 	map.insert("sqlite", &OsmdroidMap::create);
 	map.insert("gemf", &GEMFMap::create);
+	map.insert("otrk2.xml", &OruxMap::create);
 
 	return map;
 }
@@ -62,7 +64,7 @@ Map *MapList::loadFile(const QString &path, const Projection &proj, bool *isDir)
 {
 	ParserMap::iterator it;
 	QFileInfo fi(path);
-	QString suffix(fi.suffix().toLower());
+	QString suffix(fi.completeSuffix().toLower());
 	Map *map = 0;
 	QStringList errors;
 
@@ -161,6 +163,7 @@ QString MapList::formats()
 	  + qApp->translate("MapList", "Mapsforge maps") + " (*.map);;"
 	  + qApp->translate("MapList", "OziExplorer maps") + " (*.map);;"
 	  + qApp->translate("MapList", "MBTiles maps") + " (*.mbtiles);;"
+	  + qApp->translate("MapList", "Orux maps") + " (*.otrk2.xml);;"
 	  + qApp->translate("MapList", "QuickChart maps") + " (*.qct);;"
 	  + qApp->translate("MapList", "TwoNav maps") + " (*.rmap *.rtmap);;"
 	  + qApp->translate("MapList", "Osmdroid SQLite maps") + " (*.sqlite);;"

src/map/oruxmap.cpp

@@ -0,0 +1,483 @@
+#include <QPainter>
+#include <QPixmapCache>
+#include <QSqlQuery>
+#include <QSqlRecord>
+#include <QSqlField>
+#include <QFile>
+#include <QFileInfo>
+#include <QDir>
+#include <QXmlStreamReader>
+#include "pcs.h"
+#include "utm.h"
+#include "oruxmap.h"
+
+
+#define META_TYPE(type) static_cast<QMetaType::Type>(type)
+
+static bool intAttr(QXmlStreamReader &reader, const QXmlStreamAttributes &attr,
+  const QString &name, int &val)
+{
+	bool ok;
+	if (!attr.hasAttribute(name)) {
+		reader.raiseError(QString("Missing %1 attribute").arg(name));
+		return false;
+	}
+	val = attr.value(name).toInt(&ok);
+	if (!ok) {
+		reader.raiseError(QString("Invalid %1 attribute").arg(name));
+		return false;
+	}
+
+	return true;
+}
+
+static bool dblAttr(QXmlStreamReader &reader, const QXmlStreamAttributes &attr,
+  const QString &name, double &val)
+{
+	bool ok;
+	if (!attr.hasAttribute(name)) {
+		reader.raiseError(QString("Missing %1 attribute").arg(name));
+		return false;
+	}
+	val = attr.value(name).toDouble(&ok);
+	if (!ok) {
+		reader.raiseError(QString("Invalid %1 attribute").arg(name));
+		return false;
+	}
+
+	return true;
+}
+
+static bool strAttr(QXmlStreamReader &reader, const QXmlStreamAttributes &attr,
+  const QString &name, QString &val)
+{
+	if (!attr.hasAttribute(name)) {
+		reader.raiseError(QString("Missing %1 attribute").arg(name));
+		return false;
+	}
+
+	val = attr.value(name).toString();
+
+	return true;
+}
+
+static PointD corner2point(const QString &name, const QSize &size)
+{
+	if (name == "TL")
+		return PointD(0, 0);
+	else if (name == "BR")
+		return PointD(size.width(), size.height());
+	else if (name == "TR")
+		return PointD(size.width(), 0);
+	else if (name == "BL")
+		return PointD(0, size.height());
+	else
+		return PointD();
+}
+
+static Projection::Setup params2setup(const QStringList &list)
+{
+	double params[7] = {0, 0, 1, 0, 0, 0, 0};
+	bool ok;
+
+	for (int i = 1; i < list.size(); i++) {
+		if (!list.at(i).isEmpty()) {
+			params[i - 1] = list.at(i).toDouble(&ok);
+			if (!ok)
+				return Projection::Setup();
+		}
+	}
+
+	return Projection::Setup(params[0], params[1], params[2], params[3],
+	  params[4], params[5], params[6]);
+}
+
+static Projection::Setup utm2setup(const QStringList &list)
+{
+	bool ok;
+	if (list.size() < 2)
+		return Projection::Setup();
+	int zone = list.at(1).toInt(&ok);
+	return ok ? UTM::setup(zone) : Projection::Setup();
+}
+
+static GCS createGCS(const QString &datum)
+{
+	QStringList dl(datum.split(':'));
+	return (GCS::gcs(dl.first()));
+}
+
+static Projection createProjection(const GCS &gcs, const QString &name)
+{
+	PCS pcs;
+	QStringList pl(name.split(','));
+
+	if (pl.first() == "Latitude/Longitude")
+		return Projection(gcs);
+	else if (pl.first() == "UTM")
+		pcs = PCS(gcs, 9807, utm2setup(pl), 9001);
+	else if (pl.first() == "Mercator")
+		pcs = PCS(gcs, 1024, Projection::Setup(), 9001);
+	else if (pl.first() == "Transverse Mercator")
+		pcs = PCS(gcs, 9807, params2setup(pl), 9001);
+	else if (pl.first() == "Lambert Conformal Conic")
+		pcs = PCS(gcs, 9802, params2setup(pl), 9001);
+	else if (pl.first() == "Albers Equal Area")
+		pcs = PCS(gcs, 9822, params2setup(pl), 9001);
+	else if (pl.first() == "(A)Lambert Azimuthual Equal Area")
+		pcs = PCS(gcs, 9820, params2setup(pl), 9001);
+	else if (pl.first() == "Polyconic (American)")
+		pcs = PCS(gcs, 9818, params2setup(pl), 9001);
+	else if (pl.first() == "(NZTM2) New Zealand TM 2000")
+		pcs = PCS(gcs, 9807, Projection::Setup(0, 173.0, 0.9996, 1600000,
+		  10000000, NAN, NAN), 9001);
+	else if (pl.first() == "(BNG) British National Grid")
+		pcs = PCS(gcs, 9807, Projection::Setup(49, -2, 0.999601, 400000,
+		  -100000, NAN, NAN), 9001);
+	else if (pl.first() == "(IG) Irish Grid")
+		pcs = PCS(gcs, 9807, Projection::Setup(53.5, -8, 1.000035, 200000,
+		  250000, NAN, NAN), 9001);
+	else if (pl.first() == "(SG) Swedish Grid")
+		pcs = PCS(gcs, 9807, Projection::Setup(0, 15.808278, 1, 1500000, 0, NAN,
+		  NAN), 9001);
+	else if (pl.first() == "(I) France Zone I")
+		pcs = PCS(gcs, 9802, Projection::Setup(49.5, 2.337229, NAN, 600000,
+		  1200000, 48.598523, 50.395912), 9001);
+	else if (pl.first() == "(II) France Zone II")
+		pcs = PCS(gcs, 9802, Projection::Setup(46.8, 2.337229, NAN, 600000,
+		  2200000, 45.898919, 47.696014), 9001);
+	else if (pl.first() == "(III) France Zone III")
+		pcs = PCS(gcs, 9802, Projection::Setup(44.1, 2.337229, NAN, 600000,
+		  3200000, 43.199291, 44.996094), 9001);
+	else if (pl.first() == "(IV) France Zone IV")
+		pcs = PCS(gcs, 9802, Projection::Setup(42.165, 2.337229, NAN, 234.358,
+		  4185861.369, 41.560388, 42.767663), 9001);
+	else if (pl.first() == "(VICGRID) Victoria Australia")
+		pcs = PCS(gcs, 9802, Projection::Setup(-37, 145, NAN, 2500000, 4500000,
+		  -36, -38), 9001);
+	else if (pl.first() == "(VG94) VICGRID94 Victoria Australia")
+		pcs = PCS(gcs, 9802, Projection::Setup(-37, 145, NAN, 2500000, 2500000,
+		  -36, -38), 9001);
+	else if (pl.first() == "(SUI) Swiss Grid")
+		pcs = PCS(gcs, 9815, Projection::Setup(46.570866, 7.26225, 1.0, 600000,
+		  200000, 90.0, 90.0), 9001);
+	else
+		return Projection();
+
+	return Projection(pcs);
+}
+
+static Transform computeTransformation(const Projection &proj,
+  const QList<CalibrationPoint> &points)
+{
+	QList<ReferencePoint> rp;
+
+	for (int i = 0; i < points.size(); i++)
+		rp.append(points.at(i).rp(proj));
+
+	return Transform(rp);
+}
+
+void OruxMap::calibrationPoints(QXmlStreamReader &reader, const QSize &size,
+  QList<CalibrationPoint> &points)
+{
+	while (reader.readNextStartElement()) {
+		if (reader.name() == QLatin1String("CalibrationPoint")) {
+			double lon, lat;
+			QString corner;
+
+			QXmlStreamAttributes attr = reader.attributes();
+			if (!dblAttr(reader, attr, "lat", lat))
+				return;
+			if (!dblAttr(reader, attr, "lon", lon))
+				return;
+			if (!strAttr(reader, attr, "corner", corner))
+				return;
+
+			CalibrationPoint p(corner2point(corner, size), Coordinates(lon, lat));
+			points.append(p);
+
+			reader.readElementText();
+		} else
+			reader.skipCurrentElement();
+	}
+}
+
+void OruxMap::mapCalibration(QXmlStreamReader &reader, int level)
+{
+	int zoom;
+	QSize tileSize, size;
+	QString datum, projection;
+	QList<CalibrationPoint> points;
+	Projection proj;
+
+	QXmlStreamAttributes attr = reader.attributes();
+	if (!intAttr(reader, attr, "layerLevel", zoom))
+		return;
+
+	while (reader.readNextStartElement()) {
+		if (reader.name() == QLatin1String("OruxTracker"))
+			oruxTracker(reader, level + 1);
+		else if (reader.name() == QLatin1String("MapName")) {
+			QString name(reader.readElementText());
+			if (!level)
+				_name = name;
+		} else if (reader.name() == QLatin1String("MapChunks")) {
+			int xMax, yMax, width, height;
+
+			QXmlStreamAttributes attr = reader.attributes();
+			if (!intAttr(reader, attr, "xMax", xMax))
+				return;
+			if (!intAttr(reader, attr, "yMax", yMax))
+				return;
+			if (!intAttr(reader, attr, "img_width", width))
+				return;
+			if (!intAttr(reader, attr, "img_height", height))
+				return;
+			if (!strAttr(reader, attr, "datum", datum))
+				return;
+			if (!strAttr(reader, attr, "projection", projection))
+				return;
+
+			tileSize = QSize(width, height);
+			size = QSize(xMax * width, yMax * height);
+			GCS gcs(createGCS(datum));
+			if (!gcs.isValid()) {
+				reader.raiseError(QString("%1: invalid/unknown datum")
+				  .arg(datum));
+				return;
+			}
+			proj = createProjection(gcs, projection);
+			if (!proj.isValid()) {
+				reader.raiseError(QString("%1: invalid/unknown projection")
+				  .arg(projection));
+				return;
+			}
+
+			reader.readElementText();
+		} else if (reader.name() == QLatin1String("CalibrationPoints"))
+			calibrationPoints(reader, size, points);
+		else
+			reader.skipCurrentElement();
+	}
+
+	if (tileSize.isValid()) {
+		Transform t(computeTransformation(proj, points));
+		_zooms.append(Zoom(zoom, tileSize, size, proj, t));
+	}
+}
+
+void OruxMap::oruxTracker(QXmlStreamReader &reader, int level)
+{
+	while (reader.readNextStartElement()) {
+		if (reader.name() == QLatin1String("MapCalibration"))
+			mapCalibration(reader, level);
+		else
+			reader.skipCurrentElement();
+	}
+}
+
+bool OruxMap::readXML(const QString &path)
+{
+	QFile file(path);
+
+	if (!file.open(QFile::ReadOnly | QFile::Text))
+		return false;
+
+	QXmlStreamReader reader(&file);
+	if (reader.readNextStartElement()) {
+		if (reader.name() == QLatin1String("OruxTracker"))
+			oruxTracker(reader, 0);
+		else
+			reader.raiseError("Not a Orux map calibration file");
+	}
+	if (reader.error()) {
+		_errorString = QString("%1: %2").arg(reader.lineNumber())
+		  .arg(reader.errorString());
+		return false;
+	}
+
+	return true;
+}
+
+
+OruxMap::OruxMap(const QString &fileName, QObject *parent)
+  : Map(fileName, parent), _zoom(0), _mapRatio(1.0), _valid(false)
+{
+	if (!readXML(fileName))
+		return;
+
+	if (_zooms.isEmpty()) {
+		_errorString = "No usable zoom level found";
+		return;
+	}
+
+	QFileInfo fi(fileName);
+	QDir dir(fi.absoluteDir());
+	QString dbFile(dir.absoluteFilePath("OruxMapsImages.db"));
+	if (!QFileInfo::exists(dbFile)) {
+		_errorString = "Image DB file not found";
+		return;
+	}
+
+	_db = QSqlDatabase::addDatabase("QSQLITE", dbFile);
+	_db.setDatabaseName(dbFile);
+	_db.setConnectOptions("QSQLITE_OPEN_READONLY");
+
+	if (!_db.open()) {
+		_errorString = "Error opening database file";
+		return;
+	}
+
+	QSqlRecord r = _db.record("tiles");
+	if (r.isEmpty()
+	  || r.field(0).name() != "x"
+	  || META_TYPE(r.field(0).type()) != QMetaType::Int
+	  || r.field(1).name() != "y"
+	  || META_TYPE(r.field(1).type()) != QMetaType::Int
+	  || r.field(2).name() != "z"
+	  || META_TYPE(r.field(2).type()) != QMetaType::Int
+	  || r.field(3).name() != "image"
+	  || META_TYPE(r.field(3).type()) != QMetaType::QByteArray) {
+		_errorString = "Invalid table format";
+		return;
+	}
+
+	_db.close();
+	_valid = true;
+}
+
+int OruxMap::zoomFit(const QSize &size, const RectC &rect)
+{
+	if (!rect.isValid())
+		_zoom = _zooms.size() - 1;
+	else {
+		for (int i = 1; i < _zooms.size(); i++) {
+			_zoom = i;
+			QRect sbr(QPoint(ll2xy(rect.topLeft()).toPoint()),
+			  QPoint(ll2xy(rect.bottomRight()).toPoint()));
+			if (sbr.size().width() >= size.width() || sbr.size().height()
+			  >= size.height()) {
+				_zoom--;
+				break;
+			}
+		}
+	}
+
+	return _zoom;
+}
+
+int OruxMap::zoomIn()
+{
+	_zoom = qMin(_zoom + 1, _zooms.size() - 1);
+	return _zoom;
+}
+
+int OruxMap::zoomOut()
+{
+	_zoom = qMax(_zoom - 1, 0);
+	return _zoom;
+}
+
+void OruxMap::load()
+{
+	_db.open();
+}
+
+void OruxMap::unload()
+{
+	_db.close();
+}
+
+void OruxMap::setDevicePixelRatio(qreal deviceRatio, qreal mapRatio)
+{
+	Q_UNUSED(deviceRatio);
+	_mapRatio = mapRatio;
+}
+
+QPixmap OruxMap::tile(int zoom, int x, int y) const
+{
+	QSqlQuery query(_db);
+	query.prepare("SELECT image FROM tiles WHERE z=:z AND x=:x AND y=:y");
+	query.bindValue(":z", zoom);
+	query.bindValue(":x", x);
+	query.bindValue(":y", y);
+	query.exec();
+	if (!query.first())
+		return QPixmap();
+
+	QImage img(QImage::fromData(query.value(0).toByteArray()));
+	return QPixmap::fromImage(img);
+}
+
+void OruxMap::draw(QPainter *painter, const QRectF &rect, Flags flags)
+{
+	Q_UNUSED(flags);
+	const Zoom &z = _zooms.at(_zoom);
+	QSizeF ts(z.tileSize.width() / _mapRatio, z.tileSize.height() / _mapRatio);
+	QPointF tl(floor(rect.left() / ts.width()) * ts.width(),
+	  floor(rect.top() / ts.height()) * ts.height());
+
+	QSizeF s(rect.right() - tl.x(), rect.bottom() - tl.y());
+	for (int i = 0; i < ceil(s.width() / ts.width()); i++) {
+		for (int j = 0; j < ceil(s.height() / ts.height()); j++) {
+			int x = round(tl.x() * _mapRatio + i * z.tileSize.width());
+			int y = round(tl.y() * _mapRatio + j * z.tileSize.height());
+
+			QPixmap pixmap;
+			QString key = path() + "/" + QString::number(z.zoom)
+			  + "_" + QString::number(x/z.tileSize.width())
+			  + "_" + QString::number(y/z.tileSize.height());
+			if (!QPixmapCache::find(key, &pixmap)) {
+				pixmap = tile(z.zoom, x/z.tileSize.width(), y/z.tileSize.height());
+				if (!pixmap.isNull())
+					QPixmapCache::insert(key, pixmap);
+			}
+
+			if (pixmap.isNull())
+				qWarning("%s: error loading tile image", qPrintable(key));
+			else {
+				pixmap.setDevicePixelRatio(_mapRatio);
+				QPointF tp(tl.x() + i * ts.width(), tl.y() + j * ts.height());
+				painter->drawPixmap(tp, pixmap);
+			}
+		}
+	}
+}
+
+QRectF OruxMap::bounds()
+{
+	const Zoom &z = _zooms.at(_zoom);
+	return QRectF(QPointF(0, 0), z.size / _mapRatio);
+}
+
+QPointF OruxMap::ll2xy(const Coordinates &c)
+{
+	const Zoom &z = _zooms.at(_zoom);
+	QPointF p(z.transform.proj2img(z.projection.ll2xy(c)));
+	return (p / _mapRatio);
+}
+
+Coordinates OruxMap::xy2ll(const QPointF &p)
+{
+	const Zoom &z = _zooms.at(_zoom);
+	return z.projection.xy2ll(z.transform.img2proj(p * _mapRatio));
+}
+
+Map *OruxMap::create(const QString &path, const Projection &, bool *isDir)
+{
+	if (isDir)
+		*isDir = false;
+
+	return new OruxMap(path);
+}
+
+#ifndef QT_NO_DEBUG
+QDebug operator<<(QDebug dbg, const OruxMap::Zoom &zoom)
+{
+	dbg.nospace() << "Zoom(" << zoom.zoom << ", " << zoom.tileSize << ", "
+	  << zoom.size << ")";
+
+	return dbg.space();
+}
+#endif // QT_NO_DEBUG

src/map/oruxmap.h

@@ -0,0 +1,77 @@
+#ifndef ORUXMAP_H
+#define ORUXMAP_H
+
+#include <QDebug>
+#include <QSqlDatabase>
+#include "map.h"
+#include "projection.h"
+#include "transform.h"
+#include "calibrationpoint.h"
+
+class QXmlStreamReader;
+
+class OruxMap : public Map
+{
+	Q_OBJECT
+
+public:
+	OruxMap(const QString &fileName, QObject *parent = 0);
+
+	QString name() const {return _name;}
+
+	QRectF bounds();
+
+	int zoom() const {return _zoom;}
+	void setZoom(int zoom) {_zoom = zoom;}
+	int zoomFit(const QSize &size, const RectC &rect);
+	int zoomIn();
+	int zoomOut();
+
+	QPointF ll2xy(const Coordinates &c);
+	Coordinates xy2ll(const QPointF &p);
+
+	void draw(QPainter *painter, const QRectF &rect, Flags flags);
+
+	void setDevicePixelRatio(qreal deviceRatio, qreal mapRatio);
+	void load();
+	void unload();
+
+	bool isValid() const {return _valid;}
+	QString errorString() const {return _errorString;}
+
+	static Map *create(const QString &path, const Projection &, bool *isDir);
+
+private:
+	struct Zoom {
+		Zoom(int zoom, const QSize &tileSize, const QSize &size,
+		  const Projection &proj, const Transform &transform)
+		  : zoom(zoom), tileSize(tileSize), size(size), projection(proj),
+		  transform(transform) {}
+
+		int zoom;
+		QSize tileSize;
+		QSize size;
+		Projection projection;
+		Transform transform;
+	};
+
+	bool readXML(const QString &path);
+	void oruxTracker(QXmlStreamReader &reader, int level);
+	void mapCalibration(QXmlStreamReader &reader, int level);
+	void calibrationPoints(QXmlStreamReader &reader, const QSize &size,
+	  QList<CalibrationPoint> &points);
+	QPixmap tile(int zoom, int x, int y) const;
+
+	friend QDebug operator<<(QDebug dbg, const Zoom &zoom);
+
+	QString _name;
+	QList<Zoom> _zooms;
+	QSqlDatabase _db;
+	int _zoom;
+	qreal _mapRatio;
+
+	bool _valid;
+	QString _errorString;
+};
+
+#endif // ORUXMAP_H