Added support for sqlitedb maps (Locus, OsmAnd, RMaps)

gpxsee.pro

@@ -109,6 +109,7 @@ HEADERS += src/common/config.h \
     src/map/projection.h \
     src/map/ellipsoid.h \
     src/map/datum.h \
+    src/map/sqlitemap.h \
     src/map/webmercator.h \
     src/map/transversemercator.h \
     src/map/latlon.h \
@@ -289,6 +290,7 @@ SOURCES += src/main.cpp \
     src/map/emptymap.cpp \
     src/map/ozimap.cpp \
     src/map/polyconic.cpp \
+    src/map/sqlitemap.cpp \
     src/map/tar.cpp \
     src/map/atlas.cpp \
     src/map/ozf.cpp \

src/map/maplist.cpp

@@ -13,6 +13,7 @@
 #include "bsbmap.h"
 #include "kmzmap.h"
 #include "aqmmap.h"
+#include "sqlitemap.h"
 #include "invalidmap.h"
 #include "maplist.h"
 
@@ -53,6 +54,8 @@ Map *MapList::loadFile(const QString &path, bool *terminate)
 		map = new KMZMap(path);
 	else if (suffix == "aqm")
 		map = new AQMMap(path);
+	else if (suffix == "sqlitedb")
+		map = new SqliteMap(path);
 
 	return map ? map : new InvalidMap(path, "Unknown file format");
 }
@@ -106,6 +109,8 @@ QString MapList::formats()
 	  + qApp->translate("MapList", "KMZ maps") + " (*.kmz);;"
 	  + qApp->translate("MapList", "OziExplorer maps") + " (*.map);;"
 	  + qApp->translate("MapList", "MBTiles maps") + " (*.mbtiles);;"
+	  + qApp->translate("MapList", "Sqlite maps (Locus, OsmAnd, RMaps)")
+		+ " (*.sqlitedb);;"
 	  + qApp->translate("MapList", "TrekBuddy maps/atlases") + " (*.tar *.tba);;"
 	  + qApp->translate("MapList", "GeoTIFF images") + " (*.tif *.tiff);;"
 	  + qApp->translate("MapList", "TwoNav maps") + " (*.rmap *.rtmap);;"
@@ -116,7 +121,7 @@ QStringList MapList::filter()
 {
 	QStringList filter;
 	filter << "*.aqm" << "*.gmap" << "*.gmapi" << "*.img" << "*.jnx" << "*.kap"
-	  << "*.kmz" << "*.map" << "*.mbtiles" << "*.rmap" << "*.rtmap" << "*.tar"
-	  << "*.tba" << "*.tif" << "*.tiff"  << "*.xml";
+	  << "*.kmz" << "*.map" << "*.mbtiles" << "*.rmap" << "*.rtmap"
+	  << "*.sqlitedb" << "*.tar" << "*.tba" << "*.tif" << "*.tiff"  << "*.xml";
 	return filter;
 }

src/map/sqlitemap.cpp

@@ -0,0 +1,274 @@
+#include <QSqlQuery>
+#include <QSqlRecord>
+#include <QSqlField>
+#include <QPainter>
+#include <QPixmapCache>
+#include <QImageReader>
+#include <QBuffer>
+#include <QtConcurrent>
+#include "osm.h"
+#include "sqlitemap.h"
+
+
+#define META_TYPE(type) static_cast<QMetaType::Type>(type)
+
+class SqliteTile
+{
+public:
+	SqliteTile(const QPoint &xy, const QByteArray &data, const QString &key)
+	  : _xy(xy), _data(data), _key(key) {}
+
+	const QPoint &xy() const {return _xy;}
+	const QString &key() const {return _key;}
+	const QPixmap &pixmap() const {return _pixmap;}
+
+	void load() {_pixmap.loadFromData(_data);}
+
+private:
+	QPoint _xy;
+	QByteArray _data;
+	QString _key;
+	QPixmap _pixmap;
+};
+
+
+SqliteMap::SqliteMap(const QString &fileName, QObject *parent)
+  : Map(fileName, parent), _mapRatio(1.0), _valid(false)
+{
+	_db = QSqlDatabase::addDatabase("QSQLITE", fileName);
+	_db.setDatabaseName(fileName);
+	_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(4).name() != "image"
+	  || META_TYPE(r.field(4).type()) != QMetaType::QByteArray) {
+		_errorString = "Invalid table format";
+		return;
+	}
+
+	{
+		QSqlQuery query("SELECT min(z), max(z) FROM tiles", _db);
+		if (!query.first()) {
+			_errorString = "Empty tile set";
+			return;
+		}
+		_zooms = Range(17 - query.value(1).toInt(), 17 - query.value(0).toInt());
+		if (_zooms.min() < 0 || !_zooms.isValid()) {
+			_errorString = "Invalid zoom levels";
+			return;
+		}
+	}
+	_zoom = _zooms.max();
+
+	{
+		QString sql = QString("SELECT min(x), min(y), max(x), max(y) FROM tiles"
+		  " WHERE z = %1").arg(17 - _zooms.min());
+		QSqlQuery query(sql, _db);
+		query.first();
+
+		double minX = OSM::index2mercator(qMin((1<<_zooms.min()) - 1,
+		  qMax(0, query.value(0).toInt())), _zooms.min());
+		double minY = OSM::index2mercator(qMin((1<<_zooms.min()) - 1,
+		  qMax(0, query.value(1).toInt())), _zooms.min());
+		double maxX = OSM::index2mercator(qMin((1<<_zooms.min()) - 1,
+		  qMax(0, query.value(2).toInt())) + 1, _zooms.min());
+		double maxY = OSM::index2mercator(qMin((1<<_zooms.min()) - 1,
+		  qMax(0, query.value(3).toInt())) + 1, _zooms.min());
+		Coordinates tl(OSM::m2ll(QPointF(minX, -minY)));
+		Coordinates br(OSM::m2ll(QPointF(maxX, -maxY)));
+		// Workaround of broken zoom levels 0 and 1 due to numerical instability
+		tl.rlat() = qMin(tl.lat(), OSM::BOUNDS.top());
+		br.rlat() = qMax(br.lat(), OSM::BOUNDS.bottom());
+		_bounds = RectC(tl, br);
+	}
+
+	{
+		QString sql = QString("SELECT image FROM tiles LIMIT 1");
+		QSqlQuery query(sql, _db);
+		query.first();
+
+		QByteArray data = query.value(0).toByteArray();
+		QBuffer buffer(&data);
+		QImageReader reader(&buffer);
+		QSize tileSize(reader.size());
+
+		if (!tileSize.isValid() || tileSize.width() != tileSize.height()) {
+			_errorString = "Unsupported/invalid tile images";
+			return;
+		}
+		_tileSize = tileSize.width();
+	}
+
+	_db.close();
+
+	_valid = true;
+}
+
+void SqliteMap::load()
+{
+	_db.open();
+}
+
+void SqliteMap::unload()
+{
+	_db.close();
+}
+
+QRectF SqliteMap::bounds()
+{
+	return QRectF(ll2xy(_bounds.topLeft()), ll2xy(_bounds.bottomRight()));
+}
+
+int SqliteMap::limitZoom(int zoom) const
+{
+	if (zoom < _zooms.min())
+		return _zooms.min();
+	if (zoom > _zooms.max())
+		return _zooms.max();
+
+	return zoom;
+}
+
+int SqliteMap::zoomFit(const QSize &size, const RectC &rect)
+{
+	if (!rect.isValid())
+		_zoom = _zooms.max();
+	else {
+		QRectF tbr(OSM::ll2m(rect.topLeft()), OSM::ll2m(rect.bottomRight()));
+		QPointF sc(tbr.width() / size.width(), tbr.height() / size.height());
+		_zoom = limitZoom(OSM::scale2zoom(qMax(sc.x(), -sc.y()) / _mapRatio,
+		  _tileSize));
+	}
+
+	return _zoom;
+}
+
+qreal SqliteMap::resolution(const QRectF &rect)
+{
+	return OSM::resolution(rect.center(), _zoom, _tileSize);
+}
+
+int SqliteMap::zoomIn()
+{
+	_zoom = qMin(_zoom + 1, _zooms.max());
+	return _zoom;
+}
+
+int SqliteMap::zoomOut()
+{
+	_zoom = qMax(_zoom - 1, _zooms.min());
+	return _zoom;
+}
+
+void SqliteMap::setDevicePixelRatio(qreal deviceRatio, qreal mapRatio)
+{
+	Q_UNUSED(deviceRatio);
+	_mapRatio = mapRatio;
+}
+
+qreal SqliteMap::tileSize() const
+{
+	return (_tileSize / _mapRatio);
+}
+
+QByteArray SqliteMap::tileData(int zoom, const QPoint &tile) const
+{
+	QSqlQuery query(_db);
+	query.prepare("SELECT image FROM tiles WHERE z=:zoom AND x=:x AND y=:y");
+	query.bindValue(":zoom", 17 - zoom);
+	query.bindValue(":x", tile.x());
+	query.bindValue(":y", tile.y());
+	query.exec();
+
+	if (query.first())
+		return query.value(0).toByteArray();
+
+	return QByteArray();
+}
+
+void SqliteMap::draw(QPainter *painter, const QRectF &rect, Flags flags)
+{
+	Q_UNUSED(flags);
+	qreal scale = OSM::zoom2scale(_zoom, _tileSize);
+	QRectF b(bounds());
+
+
+	QPoint tile = OSM::mercator2tile(QPointF(rect.topLeft().x() * scale,
+	  -rect.topLeft().y() * scale) * _mapRatio, _zoom);
+	QPointF tl(floor(rect.left() / tileSize())
+	  * tileSize(), floor(rect.top() / tileSize()) * tileSize());
+
+	QSizeF s(qMin(rect.right() - tl.x(), b.width()),
+	  qMin(rect.bottom() - tl.y(), b.height()));
+	int width = ceil(s.width() / tileSize());
+	int height = ceil(s.height() / tileSize());
+
+
+	QList<SqliteTile> tiles;
+
+	for (int i = 0; i < width; i++) {
+		for (int j = 0; j < height; j++) {
+			QPixmap pm;
+			QPoint t(tile.x() + i, tile.y() + j);
+			QString key = path() + "-" + QString::number(_zoom) + "_"
+			  + QString::number(t.x()) + "_" + QString::number(t.y());
+
+			if (QPixmapCache::find(key, &pm)) {
+				QPointF tp(qMax(tl.x(), b.left()) + (t.x() - tile.x())
+				  * tileSize(), qMax(tl.y(), b.top()) + (t.y() - tile.y())
+				  * tileSize());
+				drawTile(painter, pm, tp);
+			} else {
+				tiles.append(SqliteTile(t, tileData(_zoom, t), key));
+			}
+		}
+	}
+
+	QFuture<void> future = QtConcurrent::map(tiles, &SqliteTile::load);
+	future.waitForFinished();
+
+	for (int i = 0; i < tiles.size(); i++) {
+		const SqliteTile &mt = tiles.at(i);
+		QPixmap pm(mt.pixmap());
+		if (pm.isNull())
+			continue;
+
+		QPixmapCache::insert(mt.key(), pm);
+
+		QPointF tp(qMax(tl.x(), b.left()) + (mt.xy().x() - tile.x())
+		  * tileSize(), qMax(tl.y(), b.top()) + (mt.xy().y() - tile.y())
+		  * tileSize());
+		drawTile(painter, pm, tp);
+	}
+}
+
+void SqliteMap::drawTile(QPainter *painter, QPixmap &pixmap, QPointF &tp)
+{
+	pixmap.setDevicePixelRatio(_mapRatio);
+	painter->drawPixmap(tp, pixmap);
+}
+
+QPointF SqliteMap::ll2xy(const Coordinates &c)
+{
+	qreal scale = OSM::zoom2scale(_zoom, _tileSize);
+	QPointF m = OSM::ll2m(c);
+	return QPointF(m.x() / scale, m.y() / -scale) / _mapRatio;
+}
+
+Coordinates SqliteMap::xy2ll(const QPointF &p)
+{
+	qreal scale = OSM::zoom2scale(_zoom, _tileSize);
+	return OSM::m2ll(QPointF(p.x() * scale, -p.y() * scale) * _mapRatio);
+}

src/map/sqlitemap.h

@@ -0,0 +1,53 @@
+#ifndef SQLITEMAP_H
+#define SQLITEMAP_H
+
+#include <QSqlDatabase>
+#include "common/range.h"
+#include "map.h"
+
+class SqliteMap : public Map
+{
+public:
+	SqliteMap(const QString &fileName, QObject *parent = 0);
+
+	QRectF bounds();
+	RectC llBounds() {return _bounds;}
+	qreal resolution(const QRectF &rect);
+
+	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 load();
+	void unload();
+	void setDevicePixelRatio(qreal deviceRatio, qreal mapRatio);
+
+	bool isValid() const {return _valid;}
+	QString errorString() const {return _errorString;}
+
+private:
+	int limitZoom(int zoom) const;
+	qreal tileSize() const;
+	QByteArray tileData(int zoom, const QPoint &tile) const;
+	void drawTile(QPainter *painter, QPixmap &pixmap, QPointF &tp);
+
+	QSqlDatabase _db;
+
+	RectC _bounds;
+	Range _zooms;
+	int _zoom;
+	int _tileSize;
+	qreal _mapRatio;
+
+	bool _valid;
+	QString _errorString;
+};
+
+#endif // SQLITEMAP_H