bool CoordinateSystemConnector::loadMetaData(IlwisObject* data)
{
Ilwis3Connector::loadMetaData(data);
CoordinateSystem *csy = static_cast<CoordinateSystem *>(data);
QString ellipsoideName;
IEllipsoid ell = getEllipsoid();
GeodeticDatum *datum = getDatum(ellipsoideName);
if ( !ell.isValid() && ellipsoideName != sUNDEF){
QString ellres = QString("ilwis://tables/ellipsoid?code=%1").arg(ellipsoideName);
if (!ell.prepare(ellres)) {
return ERROR1("No ellipsoid for this code %1",ellipsoideName);
}
}
if ( type() == itCONVENTIONALCOORDSYSTEM ) {
ConventionalCoordinateSystem *csycc = static_cast<ConventionalCoordinateSystem *>(csy);
IProjection proj = getProjection();
if ( !proj.isValid()) {
return ERROR1(ERR_NO_INITIALIZED_1, "projection");
}
csycc->setDatum(datum);
csycc->setEllipsoid(ell);
csycc->setProjection(proj);
proj->setCoordinateSystem(csycc);
proj->setParameter(Projection::pvELLCODE, ell->toProj4());
csycc->prepare();
} else if ( type() == itUNKNOWN){
//TODO other types of csy
}
return true;
}
void GeodeticDatum::fromCode(const QString &gcode)
{
InternalDatabaseConnection stmt;
QString query = QString("Select * from datum where code='%1' or code='%2'").arg(gcode).arg('D' + gcode); // temporary workaround (only working for wgs84) because file datums.csv does not have datumnames that match the ones in epsg.pcs
if (stmt.exec(query)) {
if ( stmt.next()) {
QString area = stmt.value(stmt.record().indexOf("area")).toString();
QString geocode = stmt.value(stmt.record().indexOf("code")).toString();
QString ellipsoid = stmt.value(stmt.record().indexOf("ellipsoid")).toString();
double dx = stmt.value(stmt.record().indexOf("dx")).toDouble();
double dy = stmt.value(stmt.record().indexOf("dy")).toDouble();
double dz = stmt.value(stmt.record().indexOf("dz")).toDouble();
setArea(area);
code(geocode);
setWktName(stmt.value(stmt.record().indexOf("wkt")).toString());
IEllipsoid ellips;
ellips.prepare("code=ellipsoid:" + ellipsoid);
set3TransformationParameters(dx, dy, dz, ellips);
} else {
kernel()->issues()->log(TR("No datum for this code %1").arg(gcode));
}
} else {
kernel()->issues()->logSql(stmt.lastError());
}
}
bool CoordinateSystemConnector::loadMetaData(IlwisObject* data, const IOOptions& options)
{
Ilwis3Connector::loadMetaData(data, options);
CoordinateSystem *csy = static_cast<CoordinateSystem *>(data);
QString ellipsoideName;
IEllipsoid ell = getEllipsoid();
GeodeticDatum *datum = getDatum(ellipsoideName);
if ( !ell.isValid() && ellipsoideName != sUNDEF){
QString ellres = QString("ilwis://tables/ellipsoid?code=%1").arg(ellipsoideName);
if (!ell.prepare(ellres)) {
return ERROR1("No ellipsoid for this code %1",ellipsoideName);
}
}
QString cb = _odf->value("CoordSystem", "CoordBounds");
QStringList cbparts = cb.split(" ");
if ( cbparts.size() == 4 && cbparts[0] != "-1e+308") {
bool ok1, ok2, ok3, ok4;
Envelope box( Coordinate(
cbparts[0].toDouble(&ok1),
cbparts[1].toDouble(&ok2)),
Coordinate(
cbparts[2].toDouble(&ok3),
cbparts[3].toDouble(&ok4)));
if ( !( ok1 && ok2 && ok3 && ok4)) {
return ERROR2(ERR_NO_INITIALIZED_2, TR("envelop"), csy->name());
}
csy->envelope(box);
} else {
QString type = _odf->value("CoordSystem", "Type");
if ( type == "LatLon") {
Envelope box(Coordinate(-180,-90), Coordinate(180,90));
csy->envelope(box);
}
}
if ( type() == itCONVENTIONALCOORDSYSTEM ) {
ConventionalCoordinateSystem *csycc = static_cast<ConventionalCoordinateSystem *>(csy);
IProjection proj = getProjection(csycc);
if ( !proj.isValid()) {
return ERROR1(ERR_NO_INITIALIZED_1, "projection");
}
csycc->setDatum(datum);
csycc->setEllipsoid(ell);
csycc->setProjection(proj);
proj->setCoordinateSystem(csycc);
proj->setParameter(Projection::pvELLCODE, ell->toProj4());
csycc->prepare();
} else if ( type() == itUNKNOWN){
//TODO: other types of csy
}
return true;
}
GeodeticDatum *InternalIlwisObjectFactory::createDatum(const Resource& resource, const IOOptions &options) const {
QString query;
if ( resource.code() != sUNDEF) {
QString code = resource.code();
if ( code != "") {
query = QString("Select * from datum where code = '%1'").arg(code);
}
}
if ( resource["area"] != sUNDEF) {
QString name = resource.name();
QString area = resource["area"].toString();
query = QString("Select * from datum where name='%1' and area='%1'").arg(name, area);
}
if ( query == "")
return 0;
InternalDatabaseConnection db;
if (db.exec(query) && db.next()) {
GeodeticDatum *datum = new GeodeticDatum();
QSqlRecord rec = db.record();
datum->name(rec.field("name").value().toString());
datum->setDescription(rec.field("description").value().toString());
datum->setAuthority(rec.field("authority").value().toString());
datum->setArea(rec.field("area").value().toString());
datum->code(rec.field("code").value().toString());
QString ellips = rec.field("code").value().toString();
IEllipsoid ell;
QString ellres = QString("code=ellipsoid:%1").arg(ellips);
ell.prepare(ellres);
datum->set3TransformationParameters(rec.field("dx").value().toDouble(),
rec.field("dy").value().toDouble(),
rec.field("dz").value().toDouble(),
ell);
return datum;
}
return 0;
}
IEllipsoid CoordinateSystemConnector::getEllipsoid() {
QString ell = _odf->value("CoordSystem","Ellipsoid");
if ( ell == "?")
return IEllipsoid();
if ( ell == "User Defined") {
double invf = _odf->value("Ellipsoid","1/f").toDouble();
double majoraxis = _odf->value("Ellipsoid","a").toDouble();
Ellipsoid *ellips = new Ellipsoid();
QString newName = ellips->setEllipsoid(majoraxis,invf);;
ellips->name(newName);
IEllipsoid ellipsoid(ellips);
return ellipsoid;
}
IEllipsoid ellipsoid;
QString code = name2Code(ell, "ellipsoid");
if ( code == sUNDEF)
return IEllipsoid();
QString resource = QString("ilwis://tables/ellipsoid?code=%1").arg(code);
ellipsoid.prepare(resource);
return ellipsoid;
}
