using System;

using System.Collections.Generic;

using static System.Math;

using NetTopologySuite.Geometries;

using static NetTopologySuite.Algorithm.AngleUtility;

public class Program

{

    public static void Main()

    {

        //////////FOR RESULTS LOG//////////////////////

        var questionOne = CheckArea();

        var questionTwo = isColumInTerrain();

        var questionThree = spaceBetweenColumns();

        var questionFour = CheckIfColumnIsRectangular();

        for(int i=0;i<questionOne.Count;i++) //print question 1 results

        {

            Console.WriteLine("Question 1 : Column " + i + " comply with maximum area? " + questionOne[i]);

        }

        for(int i=0;i<questionTwo.Count;i++) //print question 2 results

        {

            Console.WriteLine("Question 2 : Column " + i + " is inside terrain? " + questionTwo[i]);

        }

        for(int i=0;i<questionThree.Count;i++) //print question 3 results

        {

            Console.WriteLine("Question 3 : Column " + i + " have 7m distance to all columns? " + questionThree[i]);

        }

        for(int i=0;i<questionFour.Count;i++) //print question 4 results

        {

            Console.WriteLine("Question 4 : Column " + i + " in rectangular? " + questionFour[i]);

        }

    }

​

    //////////BASE DATA - POLYGON CREATION//////////////////////

    public static Polygon CreatePolygon(double[] xValues, double[] yValues) 

    { 

        var coordinates = new List<Coordinate>(); 

        for (int i = 0; i < xValues.Length; i++) 

        { 

            coordinates.Add(new Coordinate(xValues[i], yValues[i])); 

        } 

        var ring = new LinearRing(coordinates.ToArray()); 

        return new Polygon(ring); 

    } 

​

    //////////BASE DATA - COLUMNS//////////////////////

    public static List<Polygon> LoadColumns() 

    { 

        var columns = new List<Polygon>(); 

​

        columns.Add(CreatePolygon(new[] { 2.0, 6.0, 6.0, 2.0, 2.0 }, new[] { 3.0, 3.0, 5.0, 5.0, 3.0 })); 

        columns.Add(CreatePolygon(new[] { 9.0, 12.0, 10.0, 7.0, 9.0 }, new[] { 5.0, 8.0, 10.0, 7.0, 5.0 })); 

        columns.Add(CreatePolygon(new[] { 17.0, 17.0, 19.0, 19.0, 17.0 }, new[] { 1.0, -3.0, -3.0, 1.0, 1.0 })); 

        columns.Add(CreatePolygon(new[] { 19.0, 19.0, 24.0, 24.0, 19.0 }, new[] { 8.0, 5.0, 5.0, 6.0, 8.0 })); 

​

        return columns; 

    } 

​

​

    //////////BASE DATA - TERRAIN//////////////////////

    public static Polygon LoadTerrain() 

    { 

        return CreatePolygon(new[] { 1.0, 28.0, 23.0, -1.0, 1.0 }, new[] { 9.0, 9.0, -6.0, -2.0, 9.0 }); 

    }

​

    //////////QUESTION 1//////////////////////

    public static List<bool> CheckArea()

    {

        List<Polygon> columns = LoadColumns(); // get column 2D cross sections 

        double maxArea = 10.0; //maximum area

        List<bool> result = new List<bool>(); // result will be stored in this list 

        foreach(var column in columns) 

            { 

                var valid = true; 

                if(column.Area > maxArea)  //area property

                { 

                     valid = false; 

                } 

​

                result.Add(valid); //return list of results

            } 

        return result; 

    }

    //////////QUESTION 2//////////////////////

    public static List<bool> isColumInTerrain()

    {

        List<Polygon> columns = LoadColumns(); // get column 2D cross sections 

        Polygon terrain = LoadTerrain(); //terrain polygon

        List<bool> result = new List<bool>(); // result will be stored in this list     

        foreach(var column in columns) //loop columns

            { 

                var contains = true; 

                terrain.Contains(column); //Contains check if one geometry contains another (from docs-> class polygon)             

                result.Add(contains); //store result

            } 

        return result; //return list of results

    }

    //////////QUESTION 3//////////////////////

    public static List<bool> spaceBetweenColumns()

    {

        List<Polygon> columns = LoadColumns(); // get column 2D cross sections 

        double minDistance = 7.0;

        List<bool> result = new List<bool>(); // result will be stored in this list     

        for (int i = 0; i < columns.Count; i++) //loop columns with access to index

        { 

            bool comparisonResult = true;

            for (int j = 0; j < columns.Count; j++) //loop columns again with access to index

            {   

                if (i != j ) // dont compare to itself

                {

                    double distance = columns[i].Distance(columns[j]);

                    if (distance < minDistance)

                    {

                        comparisonResult = false;

                    }

                }

            }

            result.Add(comparisonResult); //return if loops ended

        }

        return result; //return list of results

    }

    //////////QUESTION 4//////////////////////

    public static List<bool> CheckIfColumnIsRectangular ()

    {

        List<Polygon> columns = LoadColumns(); // get column 2D cross sections 

        List<bool> result = new List<bool>(); // result will be stored in this list     

        double straightAngle = 90;

        foreach(var column in columns) //loop columns

            { 

                var coordinates = column.Coordinates; //get coordinates of this column

                bool isRectangular = true; // initial value

                if (coordinates.Length == 5) //4 vertices + closure

                {

                    for (int i = 0; i < coordinates.Length; i++) //loop columns with access to index

                    { 

                        Coordinate coord = coordinates[i];

                        Coordinate nextCoord = coordinates[(i+1)%4]; 

                        Coordinate beforeCoord = coordinates[(i+3)%4];

​

                        //compute angle between coordinates using lib helper. Could be done usinng dot product and arccos in more steps

                        double angle = AngleBetween(beforeCoord, coord, nextCoord)*(180/Math.PI); //in radians transformed in degrees

​

                        if (angle != straightAngle)

                        {

                            isRectangular = false;

                            result.Add(isRectangular); //return result to list

                            break;

                        }               

                    }

                }

                result.Add(isRectangular); // return result to list if loop didnt break

            } 

        return result; //return list of results

    }

}

​

​

​