Updated Load location data API

backend/api-gateway/Models/LocationData.cs

@@ -1,9 +1,13 @@
-namespace APIGateway.Models
+using System.Text.Json.Serialization;
+
+namespace APIGateway.Models
 {
     public class LocationDataRequest
     {
-        public string DatasetId { get; set; } = string.Empty;
-        public List<Coordinate> Location { get; set; } = new List<Coordinate>();
+        [JsonPropertyName("datasetId")]
+        public string DatasetId { get; set; }
+        [JsonPropertyName("location")]
+        public List<List<List<double>>> Location { get; set; }
     }
 
     public class Coordinate
@@ -14,7 +18,7 @@ public class Coordinate
 
     public class LocationDataResponse
     {
-        public List<DatasetItem> CurrentDatasetData { get; set; } = new List<DatasetItem>();
+        public List<DatasetItem> IndividualData { get; set; } = new List<DatasetItem>();
         public List<DatasetItem> GeneralData { get; set; } = new List<DatasetItem>();
     }
 

backend/src/BIE.Core/BIE.Core.API/ApiHelper.cs

@@ -5,12 +5,83 @@
 using Accord.Math;
 using BIE.Core.API.Controllers;
 using BieMetadata;
+using NetTopologySuite.Geometries;
+using ProjNet.CoordinateSystems;
+using ProjNet.CoordinateSystems.Transformations;
+using ProjNet.IO.CoordinateSystems;
 
 namespace BIE.Core.API;
 
 public static class ApiHelper
 {
     private static CultureInfo sCultureInfo = new CultureInfo("en-US");
+    private static string epsg27700 = @"PROJCS[""North_America_Albers_Equal_Area_Conic"",
+                        GEOGCS[""NAD83"",
+                            DATUM[""North_American_Datum_1983"",
+                                SPHEROID[""GRS 1980"",6378137,298.257222101,
+                                    AUTHORITY[""EPSG"",""7019""]],
+                                AUTHORITY[""EPSG"",""6269""]],
+                            PRIMEM[""Greenwich"",0,
+                                AUTHORITY[""EPSG"",""8901""]],
+                            UNIT[""degree"",0.0174532925199433,
+                                AUTHORITY[""EPSG"",""9122""]],
+                            AUTHORITY[""EPSG"",""4269""]],
+                        PROJECTION[""Albers_Conic_Equal_Area""],
+                        PARAMETER[""latitude_of_center"",40],
+                        PARAMETER[""longitude_of_center"",-96],
+                        PARAMETER[""standard_parallel_1"",20],
+                        PARAMETER[""standard_parallel_2"",60],
+                        PARAMETER[""false_easting"",0],
+                        PARAMETER[""false_northing"",0],
+                        UNIT[""metre"",1,
+                            AUTHORITY[""EPSG"",""9001""]],
+                        AXIS[""Easting"",EAST],
+                        AXIS[""Northing"",NORTH],
+                        AUTHORITY[""ESRI"",""102008""]
+                        "; // see http://epsg.io/27700
+    private static ICoordinateTransformation mTransformation = new CoordinateTransformationFactory().
+        CreateFromCoordinateSystems(GeographicCoordinateSystem.WGS84, 
+        (CoordinateSystem)CoordinateSystemWktReader.Parse(epsg27700));
+
+    public static double getDistance(double longitude, double latitude, Point p)
+    {
+        double[] sourceInWGS84_arr = { longitude, latitude };
+        double[] sourceInEPSG27700_arr = mTransformation.MathTransform.Transform(sourceInWGS84_arr);
+        Coordinate sourceInEPSG27700 = new Coordinate(sourceInEPSG27700_arr[0], sourceInEPSG27700_arr[1]);
+
+        double[] targetPointInWGS84 = { p.Y, p.X };
+        double[] targetPointInEpsg27700 = mTransformation.MathTransform.Transform(targetPointInWGS84);
+        Coordinate targetInEPSG27700 = new Coordinate(targetPointInEpsg27700[0], targetPointInEpsg27700[1]);
+        var distance = sourceInEPSG27700.Distance(targetInEPSG27700);
+        return distance;
+    }
+
+    public static double getArea(Polygon p)
+    {
+        // convert each point of the polygon into a new coordinate (revert x and y
+        // convert it into EPSG27700
+        // return area
+        var coordinates = p.Coordinates;
+        Console.WriteLine($"coordinates {coordinates}");
+        var transformedCoordinates = new Coordinate[coordinates.Length];
+        // Transform each coordinate
+        for (int i = 0; i < coordinates.Length; i++)
+        {
+            double[] pointWGS84 = { coordinates[i].Y, coordinates[i].X }; // Switch X and Y
+            double[] pointEPSG27700 = mTransformation.MathTransform.Transform(pointWGS84);
+            transformedCoordinates[i] = new Coordinate(pointEPSG27700[0], pointEPSG27700[1]);
+            Console.WriteLine($"transformedCoordinates[i] {transformedCoordinates[i]}");
+        }
+        var geometryFactory = new GeometryFactory(new PrecisionModel(), 27700);
+        var transformedPolygon = new Polygon(new LinearRing(transformedCoordinates), geometryFactory);
+        Console.WriteLine($"area {transformedPolygon.Area}");
+        Console.WriteLine($"transformed coords {transformedPolygon.Coordinates}");
+
+
+
+        // Return the area of the transformed polygon
+        return transformedPolygon.Area;
+    }
 
     /// <summary>
     /// Get the Bounding Box from the query parameters.
@@ -53,6 +124,38 @@ public static string GetPolygonFromBoundingBox(BoundingBox boundingBox)
             $" {bottomLong} {bottomLat}))";
     }
 
+    public static List<string> ConvertToWktPolygons(List<List<List<double>>> locations)
+    {
+        var culture = new CultureInfo("en-US");
+        var wktPolygons = new List<string>();
+
+        foreach (var polygon in locations)
+        {
+            var wktString = "POLYGON((";
+            foreach (var point in polygon)
+            {
+                if (point.Count != 2)
+                {
+                    throw new ArgumentException("Each point should have exactly two coordinates.");
+                }
+
+                var longitude = point[0].ToString(culture);
+                var latitude = point[1].ToString(culture);
+                wktString += $"{longitude} {latitude}, ";
+            }
+
+            // Close the polygon by repeating the first point
+            var firstPoint = polygon[0];
+            var firstLongitude = firstPoint[0].ToString(culture);
+            var firstLatitude = firstPoint[1].ToString(culture);
+            wktString += $"{firstLongitude} {firstLatitude}))";
+
+            wktPolygons.Add(wktString);
+        }
+
+        return wktPolygons;
+    }
+
     /// <summary>
     /// Get the polygon of the bounding box given in the queryparameters
     /// </summary>
@@ -84,7 +187,7 @@ public static bool BoxIntersection(BoundingBox box1, BoundingBox box2)
         var right2 = box2.maxX;
         var top2 = box2.maxY;
 
-        Console.WriteLine($"left1: {left1}, left2: {left2}");
+        //Console.WriteLine($"left1: {left1}, left2: {left2}");
 
         return !(right1 < left2 || right2 < left1 || top1 < bottom2 || top2 < bottom1);
     }

backend/src/BIE.Core/BIE.Core.API/BIE.Core.API.csproj

@@ -16,7 +16,9 @@
     <PackageReference Include="Accord.Math" Version="3.8.0" />
     <PackageReference Include="Microsoft.VisualStudio.Azure.Containers.Tools.Targets" Version="1.11.1" />
     <PackageReference Include="MongoDB.Driver" Version="2.26.0" />
+    <PackageReference Include="NetTopologySuite" Version="2.5.0" />
     <PackageReference Include="Newtonsoft.Json" Version="13.0.3" />
+    <PackageReference Include="ProjNet" Version="2.0.0" />
     <PackageReference Include="Swashbuckle.AspNetCore" Version="5.6.3" />
   </ItemGroup>