Suchen Sie mit Open Source GIS oder ArcGIS for Desktop nach Polygonen ohne rechten Winkel?

Jetzt digitalisieren wir einige Gebäude in einem bestimmten Bereich.

Die obligatorische Regel für diese Arbeit - in den meisten Fällen sollten Gebäude rechtwinklig sein.

Wir verwenden QGIS mit CAD-Werkzeugen für diese Arbeit, aber manchmal machen wir Fehler und erstellen Polygone mit einer unregelmäßigen Form.

Weiß jemand, wie wir mit Open Source GIS oder ArcGIS solche Polygone ohne rechten Winkel finden können?

Ich kenne kein vorhandenes Tool, aber Sie können eines in ArcPy schreiben oder GDAL / OGR wie folgt verwenden:

• Für jedes Polygon ...
  • Holen Sie sich die Geometrie
  • Folgen Sie der Wicklung des Polygons und berechnen Sie den Innenwinkel an jedem Scheitelpunkt
  • Wenn ein Winkel nicht 90 Grad beträgt, fügen Sie die FID / OID (oder ein anderes Attribut) zu einer Liste von Ablehnungen hinzu
• Drucken Sie die Liste der Ablehnungen

Nachfolgend finden Sie einen möglichen Ansatz. Die Funktion gibt true oder false zurück, je nachdem, ob das Polygon Winkel unterhalb einer bestimmten Größe aufweist oder innerhalb eines Bereichs um einen Zielwinkel liegt. Denken Sie daran, dass dies ein sehr einfacher Ansatz ist und eine geradlinige Digitalisierung voraussetzt. Ich teste auf einen Kreis, aber nicht auf Kurven oder andere Möglichkeiten, die die Funktion beeinträchtigen könnten.

angleTarget = gewünschter Winkel (z. B. 90).

edgeVariance = zulässige Waffel der geraden Linie (z. B. Richtungsänderung um 0,5 Grad zulässig).

angleVariance = zulässige Abweichung des gewünschten Winkels (z. B. 1, wenn 91 Grad in Ordnung sind).

Brian

private static bool AngleWithinTolerance(IPolygon pPoly, double angletarget, double edgeVariance, double angleVariance)
    {
        GeometryEnvironment geometryEnvironment = new GeometryEnvironment();
        IConstructAngle constructAngle = geometryEnvironment as IConstructAngle;
        IPointCollection ptcol = (IPointCollection)pPoly;
        double angle;

        //No circles!
        if (ptcol.PointCount < 3) return false;

        //Check angle made by last point first point and second point in the collection.
        angle = Math.Abs(constructAngle.ConstructThreePoint(ptcol.get_Point(ptcol.PointCount - 2), ptcol.get_Point(0), ptcol.get_Point(1)) * (180/3.14159250439667));
        if (angle < edgeVariance || (angle < angletarget + angleVariance & angle > angletarget - angleVariance))
        {
            //Angle at index 0 is OK - check all other points in collection.
            for (int x = 0; x != ptcol.PointCount - 2; x++)
            {
                angle = Math.Abs(constructAngle.ConstructThreePoint(ptcol.get_Point(x), ptcol.get_Point(x + 1), ptcol.get_Point(x + 2)) * (180 / 3.14159250439667));
                if (angle > edgeVariance & (angle > angletarget + angleVariance || angle < angletarget - angleVariance))
                {
                    return false;
                }
            }
        }
        else
        {
            return false;
        }
        //never failed.
        return true;
    }