1
0
mirror of https://github.com/tumic0/GPXSee.git synced 2024-11-28 05:34:47 +01:00

Compare commits

..

No commits in common. "ff50ffa0032fb151ccdf717c53136d72e6a051a3" and "501fd03511bb9270867cdc1ff989a3f81ec9d82b" have entirely different histories.

4 changed files with 83 additions and 80 deletions

View File

@ -7,9 +7,6 @@
#define MARKER_SIZE_SMALL 8
#define MARKER_SIZE_LARGE 16
#define PROJ(object, parent) \
((object).isNull() ? (parent) : (object))
static int markerSize(const QString &str)
{
if (str == "small")
@ -113,7 +110,7 @@ static bool isWS(char c)
return (c == 0x20 || c == 0x09 || c == 0x0A || c == 0x0D) ? true : false;
}
static bool possiblyJSONObject(QFile *file)
static bool isJSONObject(QFile *file)
{
char c;
@ -129,21 +126,12 @@ static bool possiblyJSONObject(QFile *file)
return false;
}
bool GeoJSONParser::a2c(const QJsonArray &data, const Projection &proj,
Coordinates &c)
static Coordinates coordinates(const QJsonArray &data, const Projection &proj)
{
c = (data.count() >= 2 && data.at(0).isDouble() && data.at(1).isDouble())
? proj.xy2ll(PointD(data.at(0).toDouble(), data.at(1).toDouble()))
: Coordinates();
if (c.isValid())
return true;
else {
QJsonDocument doc(QJsonDocument::fromVariant(data.toVariantList()));
_errorString = QString("%1: invalid coordinates")
.arg(QString::fromUtf8(doc.toJson(QJsonDocument::Compact)));
return false;
}
if (data.count() >= 2 && data.at(0).isDouble() && data.at(1).isDouble())
return proj.xy2ll(PointD(data.at(0).toDouble(), data.at(1).toDouble()));
else
return Coordinates();
}
bool GeoJSONParser::crs(const QJsonObject &object, Projection &proj)
@ -209,9 +197,11 @@ bool GeoJSONParser::point(const QJsonObject &object, const Projection &parent,
if (!crs(object, proj))
return false;
Coordinates c;
if (!a2c(coordinates, PROJ(proj, parent), c))
Coordinates c(::coordinates(coordinates, proj.isNull() ? parent : proj));
if (!c.isValid()) {
_errorString = "Invalid Point coordinates";
return false;
}
Waypoint waypoint(c);
if (coordinates.count() == 3 && coordinates.at(2).isDouble())
@ -236,7 +226,6 @@ bool GeoJSONParser::multiPoint(const QJsonObject &object,
Projection proj;
if (!crs(object, proj))
return false;
Coordinates c;
for (int i = 0; i < coordinates.size(); i++) {
if (!coordinates.at(i).isArray()) {
@ -244,8 +233,11 @@ bool GeoJSONParser::multiPoint(const QJsonObject &object,
return false;
} else {
QJsonArray data(coordinates.at(i).toArray());
if (!a2c(data, PROJ(proj, parent), c))
Coordinates c(::coordinates(data, proj.isNull() ? parent : proj));
if (!c.isValid()) {
_errorString = "Invalid MultiPoint coordinates";
return false;
}
Waypoint waypoint(c);
if (data.count() == 3 && data.at(2).isDouble())
@ -274,8 +266,7 @@ bool GeoJSONParser::lineString(const QJsonObject &object,
Projection proj;
if (!crs(object, proj))
return false;
SegmentData segment;
Coordinates c;
SegmentData sd;
for (int i = 0; i < coordinates.size(); i++) {
if (!coordinates.at(i).isArray()) {
@ -284,16 +275,19 @@ bool GeoJSONParser::lineString(const QJsonObject &object,
}
QJsonArray data(coordinates.at(i).toArray());
if (!a2c(data, PROJ(proj, parent), c))
Coordinates c(::coordinates(data, proj.isNull() ? parent : proj));
if (!c.isValid()) {
_errorString = "Invalid LineString coordinates";
return false;
}
Trackpoint t(c);
if (data.count() == 3 && data.at(2).isDouble())
t.setElevation(data.at(2).toDouble());
segment.append(t);
sd.append(t);
}
TrackData track(segment);
TrackData track(sd);
setTrackProperties(track, properties);
tracks.append(track);
@ -317,7 +311,6 @@ bool GeoJSONParser::multiLineString(const QJsonObject &object,
if (!crs(object, proj))
return false;
TrackData track;
Coordinates c;
for (int i = 0; i < coordinates.size(); i++) {
if (!coordinates.at(i).isArray()) {
@ -334,8 +327,11 @@ bool GeoJSONParser::multiLineString(const QJsonObject &object,
}
QJsonArray data(ls.at(j).toArray());
if (!a2c(data, PROJ(proj, parent), c))
Coordinates c(::coordinates(data, proj.isNull() ? parent : proj));
if (!c.isValid()) {
_errorString = "Invalid MultiLineString coordinates";
return false;
}
Trackpoint t(c);
if (data.count() == 3 && data.at(2).isDouble())
@ -368,8 +364,7 @@ bool GeoJSONParser::polygon(const QJsonObject &object, const Projection &parent,
Projection proj;
if (!crs(object, proj))
return false;
::Polygon poly;
Coordinates c;
::Polygon pg;
for (int i = 0; i < coordinates.size(); i++) {
if (!coordinates.at(i).isArray()) {
@ -387,16 +382,18 @@ bool GeoJSONParser::polygon(const QJsonObject &object, const Projection &parent,
}
QJsonArray point(lr.at(j).toArray());
if (!a2c(point, PROJ(proj, parent), c))
Coordinates c(::coordinates(point, proj.isNull() ? parent : proj));
if (!c.isValid()) {
_errorString = "Invalid Polygon linear ring coordinates";
return false;
}
data.append(c);
}
poly.append(data);
pg.append(data);
}
Area area(poly);
Area area(pg);
setAreaProperties(area, properties);
areas.append(area);
@ -419,14 +416,14 @@ bool GeoJSONParser::multiPolygon(const QJsonObject &object,
if (!crs(object, proj))
return false;
Area area;
Coordinates c;
for (int i = 0; i < coordinates.size(); i++) {
if (!coordinates.at(i).isArray()) {
_errorString = "Invalid MultiPolygon data";
return false;
} else {
::Polygon poly;
::Polygon pg;
QJsonArray polygon(coordinates.at(i).toArray());
for (int j = 0; j < polygon.size(); j++) {
@ -445,16 +442,18 @@ bool GeoJSONParser::multiPolygon(const QJsonObject &object,
}
QJsonArray point(lr.at(k).toArray());
if (!a2c(point, PROJ(proj, parent), c))
Coordinates c(::coordinates(point, proj.isNull() ? parent : proj));
if (!c.isValid()) {
_errorString = "Invalid MultiPolygon linear ring coordinates";
return false;
}
data.append(c);
}
poly.append(data);
pg.append(data);
}
area.append(poly);
area.append(pg);
}
}
@ -483,36 +482,37 @@ bool GeoJSONParser::geometryCollection(const QJsonObject &object,
switch (type(geometry)) {
case Point:
if (!point(geometry, PROJ(proj, parent), properties,
if (!point(geometry, proj.isNull() ? parent : proj, properties,
waypoints))
return false;
break;
case MultiPoint:
if (!multiPoint(geometry, PROJ(proj, parent), properties,
waypoints))
if (!multiPoint(geometry, proj.isNull() ? parent : proj,
properties, waypoints))
return false;
break;
case LineString:
if (!lineString(geometry, PROJ(proj, parent), properties,
tracks))
if (!lineString(geometry, proj.isNull() ? parent : proj,
properties, tracks))
return false;
break;
case MultiLineString:
if (!multiLineString(geometry, PROJ(proj, parent), properties,
tracks))
if (!multiLineString(geometry, proj.isNull() ? parent : proj,
properties, tracks))
return false;
break;
case Polygon:
if (!polygon(geometry, PROJ(proj, parent), properties, areas))
return false;
break;
case MultiPolygon:
if (!multiPolygon(geometry, PROJ(proj, parent), properties,
if (!polygon(geometry, proj.isNull() ? parent : proj, properties,
areas))
return false;
break;
case MultiPolygon:
if (!multiPolygon(geometry, proj.isNull() ? parent : proj,
properties, areas))
return false;
break;
case GeometryCollection:
if (!geometryCollection(geometry, PROJ(proj, parent),
if (!geometryCollection(geometry, proj.isNull() ? parent : proj,
properties, tracks, areas, waypoints))
return false;
break;
@ -542,22 +542,26 @@ bool GeoJSONParser::feature(const QJsonObject &object, const Projection &parent,
switch (type(geometry)) {
case Point:
return point(geometry, PROJ(proj, parent), properties, waypoints);
case MultiPoint:
return multiPoint(geometry, PROJ(proj, parent), properties,
return point(geometry, proj.isNull() ? parent : proj, properties,
waypoints);
case MultiPoint:
return multiPoint(geometry, proj.isNull() ? parent : proj,
properties, waypoints);
case LineString:
return lineString(geometry, PROJ(proj, parent), properties, tracks);
return lineString(geometry, proj.isNull() ? parent : proj,
properties, tracks);
case MultiLineString:
return multiLineString(geometry, PROJ(proj, parent), properties,
tracks);
return multiLineString(geometry, proj.isNull() ? parent : proj,
properties, tracks);
case GeometryCollection:
return geometryCollection(geometry, PROJ(proj, parent), properties,
tracks, areas, waypoints);
return geometryCollection(geometry, proj.isNull() ? parent : proj,
properties, tracks, areas, waypoints);
case Polygon:
return polygon(geometry, PROJ(proj, parent), properties, areas);
return polygon(geometry, proj.isNull() ? parent : proj, properties,
areas);
case MultiPolygon:
return multiPolygon(geometry, PROJ(proj, parent), properties, areas);
return multiPolygon(geometry, proj.isNull() ? parent : proj,
properties, areas);
default:
_errorString = geometry["type"].toString()
+ ": invalid/missing Feature geometry";
@ -580,8 +584,8 @@ bool GeoJSONParser::featureCollection(const QJsonObject &object,
return false;
for (int i = 0; i < features.size(); i++)
if (!feature(features.at(i).toObject(), PROJ(proj, parent), tracks,
areas, waypoints))
if (!feature(features.at(i).toObject(), proj.isNull() ? parent : proj,
tracks, areas, waypoints))
return false;
return true;
@ -593,8 +597,8 @@ bool GeoJSONParser::parse(QFile *file, QList<TrackData> &tracks,
{
Q_UNUSED(routes);
if (!possiblyJSONObject(file)) {
_errorString = "Not a GeoJSON file";
if (!isJSONObject(file)) {
_errorString = "Not a JSON file";
return false;
} else
file->reset();
@ -603,8 +607,8 @@ bool GeoJSONParser::parse(QFile *file, QList<TrackData> &tracks,
QJsonDocument doc(QJsonDocument::fromJson(file->readAll(), &error));
if (doc.isNull()) {
_errorString = QString("JSON parse error on offset %1: %2")
.arg(QString::number(error.offset), error.errorString());
_errorString = "JSON parse error: " + error.errorString() + " ["
+ QString::number(error.offset) + "]";
return false;
}

View File

@ -30,7 +30,6 @@ private:
FeatureCollection
};
bool a2c(const QJsonArray &data, const Projection &proj, Coordinates &c);
Type type(const QJsonObject &json);
bool crs(const QJsonObject &object, Projection &proj);
bool point(const QJsonObject &object, const Projection &parent,

View File

@ -8,6 +8,7 @@
#include "proj/polarstereographic.h"
#include "proj/obliquestereographic.h"
#include "proj/polyconic.h"
#include "proj/latlon.h"
#include "datum.h"
#include "gcs.h"
#include "pcs.h"
@ -39,7 +40,7 @@ Projection::Method::Method(int id)
Projection::Projection(const PCS &pcs)
: _gcs(pcs.gcs()), _ct(0), _units(pcs.conversion().units()),
_cs(pcs.conversion().cs())
_cs(pcs.conversion().cs()), _geographic(false)
{
const Ellipsoid &ellipsoid = _gcs.datum().ellipsoid();
const Projection::Setup &setup = pcs.conversion().setup();
@ -114,7 +115,7 @@ Projection::Projection(const PCS &pcs)
}
Projection::Projection(const GCS &gcs, const CoordinateSystem &cs)
: _gcs(gcs), _units(LinearUnits(9001)), _cs(cs)
: _gcs(gcs), _units(LinearUnits(9001)), _cs(cs), _geographic(true)
{
_ct = new LatLon(gcs.angularUnits());
}
@ -124,6 +125,7 @@ Projection::Projection(const Projection &p)
_gcs = p._gcs;
_units = p._units;
_ct = p._ct ? p._ct->clone() : 0;
_geographic = p._geographic;
_cs = p._cs;
}
@ -140,6 +142,7 @@ Projection &Projection::operator=(const Projection &p)
_gcs = p._gcs;
_units = p._units;
_ct = p._ct ? p._ct->clone() : 0;
_geographic = p._geographic;
_cs = p._cs;
}
@ -152,7 +155,7 @@ bool Projection::operator==(const Projection &p) const
return false;
return (*_ct == *p._ct && _gcs == p._gcs && _units == p._units
&& _cs == p._cs);
&& _cs == p._cs && _geographic == p._geographic);
}
PointD Projection::ll2xy(const Coordinates &c) const

View File

@ -7,7 +7,6 @@
#include "linearunits.h"
#include "coordinatesystem.h"
#include "gcs.h"
#include "proj/latlon.h"
class PCS;
class CT;
@ -71,7 +70,7 @@ public:
int _id;
};
Projection() : _ct(0) {}
Projection() : _ct(0), _geographic(false) {}
Projection(const Projection &p);
Projection(const PCS &pcs);
Projection(const GCS &gcs, const CoordinateSystem &cs
@ -91,10 +90,7 @@ public:
// and except of WMTS/WMS it is not needed.
return (_gcs.isValid() && _ct != 0 && _units.isValid());
}
bool isGeographic() const
{
return (dynamic_cast<const LatLon*>(_ct) != 0);
}
bool isGeographic() const {return _geographic;}
PointD ll2xy(const Coordinates &c) const;
Coordinates xy2ll(const PointD &p) const;
@ -107,6 +103,7 @@ private:
const CT *_ct;
LinearUnits _units;
CoordinateSystem _cs;
bool _geographic;
};
#ifndef QT_NO_DEBUG