/****************************************************************************** * * Project: MSSQL Spatial driver * Purpose: Implements OGRMSSQLGeometryValidator class to create valid *SqlGeometries. Author: Tamas Szekeres, szekerest at gmail.com * ****************************************************************************** * Copyright (c) 2010, Tamas Szekeres * * SPDX-License-Identifier: MIT ****************************************************************************/ #include "cpl_conv.h" #include "ogr_mssqlspatial.h" /************************************************************************/ /* OGRMSSQLGeometryValidator() */ /************************************************************************/ OGRMSSQLGeometryValidator::OGRMSSQLGeometryValidator(const OGRGeometry *poGeom, int geomColumnType) : poOriginalGeometry(poGeom), nGeomColumnType(geomColumnType), bIsValid(IsValid(poGeom)) { } /************************************************************************/ /* IsValidLatLon() */ /************************************************************************/ static double MakeValidLatitude(double latitude) { if (latitude < -90) return -90; if (latitude > 90.0) return 90.0; return latitude; } static double MakeValidLongitude(double longitude) { if (longitude < -15069.0) return -15069.0; if (longitude > 15069.0) return 15069.0; return longitude; } bool OGRMSSQLGeometryValidator::IsValidLatLon(double longitude, double latitude) { if (MakeValidLatitude(latitude) != latitude) { if (poValidGeometry == nullptr) CPLError(CE_Warning, CPLE_NotSupported, "Latitude values must be between -90 and 90 degrees"); return false; } if (MakeValidLongitude(longitude) != longitude) { if (poValidGeometry == nullptr) CPLError( CE_Warning, CPLE_NotSupported, "Longitude values must be between -15069 and 15069 degrees"); return false; } return true; } /************************************************************************/ /* IsValidCircularZ() */ /************************************************************************/ bool OGRMSSQLGeometryValidator::IsValidCircularZ(double z1, double z2) { if (z1 != z2) { if (poValidGeometry == nullptr) CPLError(CE_Warning, CPLE_NotSupported, "Circular arc segments with Z values must have equal Z " "value for all 3 points"); return false; } return true; } /************************************************************************/ /* IsValidPolygonRingCount() */ /************************************************************************/ bool OGRMSSQLGeometryValidator::IsValidPolygonRingCount(const OGRCurve *poGeom) { if (poGeom->getNumPoints() < 4) { if (poValidGeometry == nullptr) CPLError( CE_Warning, CPLE_NotSupported, "Each ring of a polygon must contain at least four points"); return false; } return true; } /************************************************************************/ /* IsValidPolygonRingClosed() */ /************************************************************************/ bool OGRMSSQLGeometryValidator::IsValidPolygonRingClosed(const OGRCurve *poGeom) { if (poGeom->get_IsClosed() == FALSE) { if (poValidGeometry == nullptr) CPLError(CE_Warning, CPLE_NotSupported, "Each ring of a polygon must have the same start and end " "points."); return false; } return true; } /************************************************************************/ /* ValidatePoint() */ /************************************************************************/ bool OGRMSSQLGeometryValidator::IsValid(const OGRPoint *poGeom) { if (poGeom->IsEmpty()) return true; if (nGeomColumnType == MSSQLCOLTYPE_GEOGRAPHY) return IsValidLatLon(poGeom->getX(), poGeom->getY()); return true; } void OGRMSSQLGeometryValidator::MakeValid(OGRPoint *poGeom) { if (poGeom->IsEmpty()) return; if (nGeomColumnType == MSSQLCOLTYPE_GEOGRAPHY) { poGeom->setX(MakeValidLongitude(poGeom->getX())); poGeom->setY(MakeValidLatitude(poGeom->getY())); } } /************************************************************************/ /* ValidateMultiPoint() */ /************************************************************************/ bool OGRMSSQLGeometryValidator::IsValid(const OGRMultiPoint *poGeom) { if (nGeomColumnType == MSSQLCOLTYPE_GEOGRAPHY) { for (const auto point : *poGeom) { if (!IsValid(point)) return false; } } return true; } void OGRMSSQLGeometryValidator::MakeValid(OGRMultiPoint *poGeom) { if (nGeomColumnType == MSSQLCOLTYPE_GEOGRAPHY) { for (auto point : *poGeom) { MakeValid(point); } } } /************************************************************************/ /* ValidateSimpleCurve() */ /************************************************************************/ bool OGRMSSQLGeometryValidator::IsValid(const OGRSimpleCurve *poGeom) { if (nGeomColumnType == MSSQLCOLTYPE_GEOGRAPHY) { const int numPoints = poGeom->getNumPoints(); for (int i = 0; i < numPoints; i++) { if (!IsValidLatLon(poGeom->getX(i), poGeom->getY(i))) return false; } } return true; } void OGRMSSQLGeometryValidator::MakeValid(OGRSimpleCurve *poGeom) { if (nGeomColumnType == MSSQLCOLTYPE_GEOGRAPHY) { const int numPoints = poGeom->getNumPoints(); for (int i = 0; i < numPoints; i++) { poGeom->setPoint(i, MakeValidLongitude(poGeom->getX(i)), MakeValidLatitude(poGeom->getY(i))); } } } /************************************************************************/ /* ValidateCircularString() */ /************************************************************************/ bool OGRMSSQLGeometryValidator::IsValid(const OGRCircularString *poGeom) { if (!IsValid(poGeom->toSimpleCurve())) return false; if (poGeom->Is3D()) { const int numPoints = poGeom->getNumPoints(); for (int i = 1; i < numPoints; i++) { if (!IsValidCircularZ(poGeom->getZ(i), poGeom->getZ(0))) { return false; } } } return true; } void OGRMSSQLGeometryValidator::MakeValid(OGRCircularString *poGeom) { MakeValid(poGeom->toSimpleCurve()); if (poGeom->Is3D()) { const int numPoints = poGeom->getNumPoints(); for (int i = 1; i < numPoints; i++) { poGeom->setZ(i, poGeom->getZ(0)); } } } /************************************************************************/ /* ValidateCompoundCurve() */ /************************************************************************/ bool OGRMSSQLGeometryValidator::IsValid(const OGRCompoundCurve *poGeom) { for (const auto poCurve : *poGeom) { switch (wkbFlatten(poCurve->getGeometryType())) { case wkbLineString: if (!IsValid(poCurve->toLineString())) return false; break; case wkbCircularString: if (!IsValid(poCurve->toCircularString())) return false; break; default: break; } } return true; } void OGRMSSQLGeometryValidator::MakeValid(OGRCompoundCurve *poGeom) { for (auto poCurve : *poGeom) { switch (wkbFlatten(poCurve->getGeometryType())) { case wkbLineString: MakeValid(poCurve->toLineString()); break; case wkbCircularString: MakeValid(poCurve->toCircularString()); break; default: break; } } } /************************************************************************/ /* ValidateMultiLineString() */ /************************************************************************/ bool OGRMSSQLGeometryValidator::IsValid(const OGRMultiLineString *poGeom) { if (nGeomColumnType == MSSQLCOLTYPE_GEOGRAPHY) { for (const auto part : *poGeom) { if (!IsValid(part)) return false; } } return true; } void OGRMSSQLGeometryValidator::MakeValid(OGRMultiLineString *poGeom) { if (nGeomColumnType == MSSQLCOLTYPE_GEOGRAPHY) { for (auto part : *poGeom) { MakeValid(part); } } } /************************************************************************/ /* ValidatePolygon() */ /************************************************************************/ void OGRMSSQLGeometryValidator::MakeValid(OGRPolygon *poGeom) { OGRMSSQLGeometryValidator::MakeValid(poGeom->toCurvePolygon()); poGeom->closeRings(); } /************************************************************************/ /* ValidateCurvePolygon() */ /************************************************************************/ bool OGRMSSQLGeometryValidator::IsValid(const OGRCurvePolygon *poGeom) { if (poGeom->IsEmpty()) return true; for (const auto part : *poGeom) { if (!IsValid(part)) return false; if (!IsValidPolygonRingCount(part)) return false; if (!IsValidPolygonRingClosed(part)) return false; } return true; } void OGRMSSQLGeometryValidator::MakeValid(OGRCurvePolygon *poGeom) { if (poGeom->IsEmpty()) return; for (auto part : *poGeom) { MakeValid(part); } } /************************************************************************/ /* ValidateMultiPolygon() */ /************************************************************************/ bool OGRMSSQLGeometryValidator::IsValid(const OGRMultiPolygon *poGeom) { for (const auto part : *poGeom) { if (!IsValid(part)) return false; } return true; } void OGRMSSQLGeometryValidator::MakeValid(OGRMultiPolygon *poGeom) { for (auto part : *poGeom) { MakeValid(part); } } /************************************************************************/ /* ValidateGeometryCollection() */ /************************************************************************/ bool OGRMSSQLGeometryValidator::IsValid(const OGRGeometryCollection *poGeom) { for (const auto part : *poGeom) { if (!IsValid(part)) return false; } return true; } void OGRMSSQLGeometryValidator::MakeValid(OGRGeometryCollection *poGeom) { for (auto part : *poGeom) { MakeValid(part); } } /************************************************************************/ /* ValidateGeometry() */ /************************************************************************/ bool OGRMSSQLGeometryValidator::IsValid(const OGRGeometry *poGeom) { if (!poGeom) return false; switch (wkbFlatten(poGeom->getGeometryType())) { case wkbPoint: return IsValid(poGeom->toPoint()); case wkbLineString: return IsValid(poGeom->toSimpleCurve()); case wkbPolygon: return IsValid(poGeom->toPolygon()); case wkbCurvePolygon: return IsValid(poGeom->toCurvePolygon()); case wkbMultiPoint: return IsValid(poGeom->toMultiPoint()); case wkbMultiLineString: return IsValid(poGeom->toMultiLineString()); case wkbCircularString: return IsValid(poGeom->toCircularString()); case wkbCompoundCurve: return IsValid(poGeom->toCompoundCurve()); case wkbMultiPolygon: return IsValid(poGeom->toMultiPolygon()); case wkbGeometryCollection: return IsValid(poGeom->toGeometryCollection()); default: break; } return false; } void OGRMSSQLGeometryValidator::MakeValid(OGRGeometry *poGeom) { if (!poGeom) return; switch (wkbFlatten(poGeom->getGeometryType())) { case wkbPoint: MakeValid(poGeom->toPoint()); break; case wkbLineString: MakeValid(poGeom->toSimpleCurve()); break; case wkbPolygon: MakeValid(poGeom->toPolygon()); break; case wkbCurvePolygon: MakeValid(poGeom->toCurvePolygon()); break; case wkbMultiPoint: MakeValid(poGeom->toMultiPoint()); break; case wkbMultiLineString: MakeValid(poGeom->toMultiLineString()); break; case wkbCircularString: MakeValid(poGeom->toCircularString()); break; case wkbCompoundCurve: MakeValid(poGeom->toCompoundCurve()); break; case wkbMultiPolygon: MakeValid(poGeom->toMultiPolygon()); break; case wkbGeometryCollection: MakeValid(poGeom->toGeometryCollection()); break; default: break; } } bool OGRMSSQLGeometryValidator::ValidateGeometry(OGRGeometry *poGeom) { poValidGeometry.reset(); if (!IsValid(poGeom)) { poValidGeometry.reset(poGeom->clone()); MakeValid(poValidGeometry.get()); return false; } return true; } /************************************************************************/ /* GetValidGeometryRef() */ /************************************************************************/ const OGRGeometry *OGRMSSQLGeometryValidator::GetValidGeometryRef() const { if (bIsValid || poOriginalGeometry == nullptr) return poOriginalGeometry; if (poValidGeometry) { CPLError(CE_Warning, CPLE_NotSupported, "Invalid geometry has been converted from %s to %s.", poOriginalGeometry->getGeometryName(), poValidGeometry->getGeometryName()); } else { CPLError(CE_Warning, CPLE_NotSupported, "Invalid geometry has been converted from %s to null.", poOriginalGeometry->getGeometryName()); } return poValidGeometry.get(); }