What is valid geometry?
Well, this checks for basic spatial integrity depending on the data type, such as polygons not having inner and outer rings crossing, and multilines not retracing back over the same points.
Part of the documentation from the OGC standards for simple feature access describe validity of polygons being:
- Polygons are topologically closed
- The boundary of a Polygon consists of a set of LinearRings that make up its exterior and interior boundaries
- No two Rings in the boundary cross, and the Rings in the boundary of a Polygon may intersect at a Point but only as a tangent
- A Polygon may not have cut lines, spikes or punctures
- The interior of every Polygon is a connected point set
- The exterior of a Polygon with 1 or more holes is not connected. Each hole defines a connected component of the exterior
STIsValid Method
The STIsValid method can be used to detect if you have shape features that are invalid. The method will return 0 if the feature is invalid
e.g. The following TSQL command will count the number of records in the table with invalid shape features
select count(id) from cadastre where geom.STIsValid() = 0
MakeValid Method
Luckily for us there is a method in SQL Server to fix invalid shape data - MakeValid.
Be aware that MakeValid could alter your shape in a way that you consider incorrect, depending on the original problem with the shape. So you should really investigate any shapes that are invalid using STIsValid to find them, and then if using MakeValid to resolve them, check the resulting shapes produced by MakeValid before performing an update.
My sample dataset only had a handful of errors caused by overshoots etc, so I was able to resolve my data with the following command
update cadastre set geom = geom.MakeValid() where geom.STIsValid() = 0
If you have a lot of invalid shapes you may want to perform commits in some kind of unit of work, such as this script that commits every 1000 updated rows
declare @rowCount int set @rowCount = 9999 while @rowCount > 0 begin begin tran update top(1000) cadastre set geom = geom.MakeValid() where geom.STIsValid() = 0 set @rowCount = @@rowcount commit end
Thank you for your post, it was a great help to solve my problems regarding geometries.
ReplyDeleteThanks - exactly what I needed. Stan
ReplyDelete