对接搜救模型

Models/SimulationCommon/SeaSJ.cs

@@ -24,6 +24,19 @@ public class Current
     public double Direction { get; set; }
 }
 
+public class NCread
+{
+        
+    //各数组来源文件'Text_readNC.cs'；时间范围：2024-06-04T00:00:00 ... 2024-06-04T23:00:00；
+    //地点范围:落水点为中心100公里
+    public float[] longitudeArray = new float[10];//经度一维数组来源自文件'Text_readNC.cs'
+    public float[] latitudeArray = new float[7];//纬度一维数组来源自文件'Text_readNC.cs'
+    public float[][][] u10Array ;//风的10米U（向东）分量三维数组来源自文件'Text_readNC.cs'
+    public float[][][] v10Array;//风的10米V（向北）分量三维数组来源自文件'Text_readNC.cs'
+    public float[][][] p140208Array;//海洋上空的自由对流速度三维数组来源自文件'Text_readNC.cs'
+    public float[][][] mwdArray;//平均波向（单位：度；0度表示北方，90度表示东方）三维数组来源自文件'Text_readNC.cs'       
+
+}
 
 public class SeaSJ
 {
@@ -38,15 +51,15 @@ public class SeaSJ
     }
 
     // 交付类方法
-    public static List<double[]> GetDrift(double[] initialPosition, double dt, double totalTime)
+    public static List<double[]> GetDrift(NCread nCread,double[] initialPosition, double dt, double totalTime)
     {
         // run 获取轨迹的函数
         // initialPosition --> 初始位置; dt --> 时间步长; totalTime --> 总时长
-        List<double[]> trajectory = CalculateDriftTrajectory(initialPosition, dt, totalTime);
+        List<double[]> trajectory = CalculateDriftTrajectory(nCread,initialPosition, dt, totalTime);
         return trajectory;
     }
 
-    public static List<double[]> CalculateDriftTrajectory(double[] initialPosition, double dt, double totalTime)
+    public static List<double[]> CalculateDriftTrajectory(NCread nCread,double[] initialPosition, double dt, double totalTime)
     {
         int timeSteps = (int)(totalTime / dt);
         List<double[]> trajectory = new List<double[]>();
@@ -57,8 +70,8 @@ public class SeaSJ
             double[] currentPos = trajectory[t - 1];
 
             // 动态获取当前位置的风力和洋流数据
-            double[] windVelocity = GetWindVelocityFromAPI(currentPos[0], currentPos[1], t);
-            double[] currentVelocity = GetCurrentVelocityFromAPI(currentPos[0], currentPos[1], t);
+            double[] windVelocity = GetWindVelocityFromAPI(nCread,currentPos[0], currentPos[1], t);
+            double[] currentVelocity = GetCurrentVelocityFromAPI(nCread,currentPos[0], currentPos[1], t);
 
             // 计算漂移速度
             double[] driftVelocity =
@@ -80,49 +93,89 @@ public class SeaSJ
         return trajectory;
     }
 
-    public static double[] GetWindVelocityFromAPI(double latitude, double longitude, int time)
+    public static double[] GetWindVelocityFromAPI(NCread windNCread,double latitude, double longitude, int time)
     {
-        // 风力数据来自于 double latitude, double longitude, string time
-        WeatherResponse weatherResponse = new WeatherResponse() { };
+        float[] longitudeArray = windNCread.longitudeArray;//经度一维数组来源自文件'Text_readNC.cs'
+        float[] latitudeArray = windNCread.latitudeArray;//纬度一维数组来源自文件'Text_readNC.cs'
+        float[][][] u10Array = windNCread.u10Array;//风的10米U（向东）分量三维数组来源自文件'Text_readNC.cs'
+        float[][][] v10Array = windNCread.v10Array;//风的10米V（向北）分量三维数组来源自文件'Text_readNC.cs'
+
+        int longitudeNum = 0;
+        int latitudeNum = 0;
+
+        //定义NC文件中读取不到的坐标点海洋数据信息（后续可采用插值法等）
+        //经度连续化
+        for (int i = 0; i < 10; i++)
+        {
+            if (longitude >= longitudeArray[i] && longitude < longitudeArray[i + 1])
+            {
+                longitude = longitudeArray[i];
+                longitudeNum = i;
+            }
+        }
 
-        double windSpeed = weatherResponse.Wind.Speed;
-        double windDirection = weatherResponse.Wind.Deg;
+        //纬度连续化
+        for (int i = 0; i < 7; i++)
+        {
+            if (latitude >= latitudeArray[i] && latitude < latitudeArray[i + 1])
+            {
+                latitude = latitudeArray[i];
+                latitudeNum = i;
+            }
+        }
 
-        // 将风速和风向转换为X和Y分量（单位：海里每小时）
-        double windSpeedInKnots = windSpeed / 1.852;
-        double windDirectionInRadians = windDirection * (Math.PI / 180);
-        double windX = windSpeedInKnots * Math.Cos(windDirectionInRadians);
-        double windY = windSpeedInKnots * Math.Sin(windDirectionInRadians);
+        double windX = (double)u10Array[time][latitudeNum][longitudeNum];
+        double windY = (double)v10Array[time][latitudeNum][ longitudeNum];
 
         return new double[] { windX, windY };
     }
 
-    public static double[] GetCurrentVelocityFromAPI(double latitude, double longitude, int time)
+    public static double[] GetCurrentVelocityFromAPI(NCread windNCread,double latitude, double longitude, int time)
     {
-        // 洋流数据来自于 double latitude, double longitude, string time
-        CurrentResponse currentResponse = new CurrentResponse() { };
+        float[] longitudeArray = windNCread.longitudeArray;//经度一维数组来源自文件'Text_readNC.cs'
+        float[] latitudeArray = windNCread.latitudeArray;//纬度一维数组来源自文件'Text_readNC.cs'
+        float[][][] u10Array = windNCread.u10Array;//风的10米U（向东）分量三维数组来源自文件'Text_readNC.cs'
+        float[][][] v10Array = windNCread.v10Array;//风的10米V（向北）分量三维数组来源自文件'Text_readNC.cs'
+
+        int longitudeNum = 0;
+        int latitudeNum = 0;
+
+        //定义NC文件中读取不到的坐标点海洋数据信息（后续可采用插值法等）
+        //经度连续化
+        for (int i = 0; i < 10; i++)
+        {
+            if (longitude >= longitudeArray[i] && longitude < longitudeArray[i + 1])
+            {
+                longitude = longitudeArray[i];
+                longitudeNum = i;
+            }
+        }
 
-        double currentSpeed = currentResponse.Current.Speed;
-        double currentDirection = currentResponse.Current.Direction;
+        //纬度连续化
+        for (int i = 0; i < 7; i++)
+        {
+            if (latitude >= latitudeArray[i] && latitude < latitudeArray[i + 1])
+            {
+                latitude = latitudeArray[i];
+                latitudeNum = i;
+            }
+        }
 
-        // 将洋流速度和方向转换为X和Y分量（单位：海里每小时）
-        double currentSpeedInKnots = currentSpeed / 1.852;
-        double currentDirectionInRadians = currentDirection * (Math.PI / 180);
-        double currentX = currentSpeedInKnots * Math.Cos(currentDirectionInRadians);
-        double currentY = currentSpeedInKnots * Math.Sin(currentDirectionInRadians);
+        double windX = (double)u10Array[time][ latitudeNum][longitudeNum];
+        double windY = (double)v10Array[time][latitudeNum][ longitudeNum];
 
-        return new double[] { currentX, currentY };
+        return new double[] { windX, windY };
     }
     
     // 交付部分方法代码
-    public static Tuple<Point2f[], double> getminEnclosingRect(List<Tuple<double, double>> latLonList)
+    public static Tuple<Point2f[], double> getminEnclosingRect(List<double[]> latLonList)
     {
         // 转换经纬度为墨卡托投影坐标，并添加到点集合
         List<Point2f> pointList = new List<Point2f>();
         foreach (var latLon in latLonList)
         {
-            double x = Rectangular_Area_Search_Function.MokatuoLat(latLon.Item1);
-            double y = Rectangular_Area_Search_Function.MokatuoLon(latLon.Item2);
+            double x = Rectangular_Area_Search_Function.MokatuoLat(latLon[0]);
+            double y = Rectangular_Area_Search_Function.MokatuoLon(latLon[1]);
             pointList.Add(new Point2f((float)x, (float)y));
         }
 

Models/SimulationCommon/Text_readNC.cs

@@ -0,0 +1,239 @@
+using System.Net;
+using Newtonsoft.Json.Linq;
+
+namespace SimulationCommon;
+
+public class Text_readNC
+{
+    public NCread windNCread = new NCread();
+    public void GetNCData()
+    {
+        string url1 = "http://10.130.100.5:7785//rescue-platform-service/api/v1/dem/getNcInfos?ncPath=D:/data/haiyang/DATA.nc&varName=v10&centerLon=116.41992593821296&centerLat=40.18801994965735";
+
+        // Create a GET request to the specified URL
+        HttpWebRequest request1 = (HttpWebRequest)WebRequest.Create(url1);
+        request1.Method = "GET";
+
+        // Send request and get response
+        using (HttpWebResponse response1 = (HttpWebResponse)request1.GetResponse())
+        {
+            using (StreamReader reader1 = new StreamReader(response1.GetResponseStream()))
+            {
+                string result1 = reader1.ReadToEnd();
+
+                // Parse the JSON string into a JObject
+                JObject resultObject1 = JObject.Parse(result1);
+
+                // Extract the 'data' field from the parsed JObject
+                JArray dataArray1 = (JArray)resultObject1["data"];
+
+                // Determine the dimensions of the data array
+                int timeSteps = dataArray1.Count;
+                int latitudes = ((JArray)dataArray1[0]).Count;
+                int longitudes = ((JArray)dataArray1[0][0]).Count;
+
+                // Initialize float[][][] array
+                float[][][] v10Array = new float[timeSteps][][];
+
+                // Convert each element to float and populate v10Array
+                for (int t = 0; t < timeSteps; t++)
+                {
+                    v10Array[t] = new float[latitudes][];
+
+                    for (int lat = 0; lat < latitudes; lat++)
+                    {
+                        v10Array[t][lat] = new float[longitudes];
+
+                        for (int lon = 0; lon < longitudes; lon++)
+                        {
+                            v10Array[t][lat][lon] = (float)dataArray1[t][lat][lon];
+                        }
+                    }
+                }
+
+                JArray lonArray = (JArray)resultObject1["data_lon"];
+                JArray lanArray = (JArray)resultObject1["data_lat"];
+                //JArray tArray = (JArray)resultObject["data_time"];
+
+                int lons = lonArray.Count;
+                int lans = lanArray.Count;
+                //int times = tArray.Count;
+
+                float[] longitudeArray = new float[lons];
+                float[] latitudeArray = new float[lans];
+                //float[] timeArray = new float[times];
+
+                for (int i = 0; i < lons; i++)
+                {
+                    longitudeArray[i] = (float)lonArray[i];
+                }
+
+                for (int j = 0; j < lans; j++)
+                {
+                    latitudeArray[j] = (float)lanArray[j];
+                }
+
+                /* for (int k = 0; k < lons; k++)
+                   {
+                       timeArray[k] = (float)tArray[k];
+                   }*/
+
+                windNCread.latitudeArray = latitudeArray;
+                windNCread.longitudeArray = longitudeArray;
+                windNCread.v10Array = v10Array;
+                // Example: Print the first value of the array
+                Console.WriteLine("a value of v10Array: " + v10Array[0][6][9]);
+                Console.WriteLine("a value of longitudeArray: " + longitudeArray[0]);
+                Console.WriteLine("a value of latitudeArray: " + latitudeArray[0]);
+                //Console.WriteLine("a value of timeArray: " + timeArray[0]);
+            }
+        }
+
+        string url2 = "http://10.130.100.5:7785//rescue-platform-service/api/v1/dem/getNcInfos?ncPath=D:/data/haiyang/DATA.nc&varName=u10&centerLon=116.41992593821296&centerLat=40.18801994965735";
+
+        // Create a GET request to the specified URL
+        HttpWebRequest request2 = (HttpWebRequest)WebRequest.Create(url2);
+        request1.Method = "GET";
+
+        // Send request and get response
+        using (HttpWebResponse response2 = (HttpWebResponse)request2.GetResponse())
+        {
+            using (StreamReader reader2 = new StreamReader(response2.GetResponseStream()))
+            {
+                string result2 = reader2.ReadToEnd();
+
+                // Parse the JSON string into a JObject
+                JObject resultObject2 = JObject.Parse(result2);
+
+                // Extract the 'data' field from the parsed JObject
+                JArray dataArray2 = (JArray)resultObject2["data"];
+
+                // Determine the dimensions of the data array
+                int timeSteps2 = dataArray2.Count;
+                int latitudes2 = ((JArray)dataArray2[0]).Count;
+                int longitudes2 = ((JArray)dataArray2[0][0]).Count;
+
+                // Initialize float[][][] array
+                float[][][] u10Array = new float[timeSteps2][][];
+
+                // Convert each element to float and populate v10Array
+                for (int t = 0; t < timeSteps2; t++)
+                {
+                    u10Array[t] = new float[latitudes2][];
+
+                    for (int lat = 0; lat < latitudes2; lat++)
+                    {
+                        u10Array[t][lat] = new float[longitudes2];
+
+                        for (int lon = 0; lon < longitudes2; lon++)
+                        {
+                            u10Array[t][lat][lon] = (float)dataArray2[t][lat][lon];
+                        }
+                    }
+                }
+                
+             
+                windNCread.u10Array = u10Array;
+                Console.WriteLine("a value of u10Array: " + u10Array[0][6][9]);
+
+            }
+        }
+
+        string url3 = "http://10.130.100.5:7785//rescue-platform-service/api/v1/dem/getNcInfos?ncPath=D:/data/haiyang/DATA.nc&varName=p140208&centerLon=116.41992593821296&centerLat=40.18801994965735";
+
+        // Create a GET request to the specified URL
+        HttpWebRequest request3 = (HttpWebRequest)WebRequest.Create(url3);
+        request1.Method = "GET";
+
+        // Send request and get response
+        using (HttpWebResponse response3 = (HttpWebResponse)request3.GetResponse())
+        {
+            using (StreamReader reader3 = new StreamReader(response3.GetResponseStream()))
+            {
+                string result3 = reader3.ReadToEnd();
+
+                // Parse the JSON string into a JObject
+                JObject resultObject3 = JObject.Parse(result3);
+
+                // Extract the 'data' field from the parsed JObject
+                JArray dataArray3 = (JArray)resultObject3["data"];
+
+                // Determine the dimensions of the data array
+                int timeSteps3 = dataArray3.Count;
+                int latitudes3 = ((JArray)dataArray3[0]).Count;
+                int longitudes3 = ((JArray)dataArray3[0][0]).Count;
+
+                // Initialize float[][][] array
+                float[][][] p140208Array = new float[timeSteps3][][];
+
+                // Convert each element to float and populate v10Array
+                for (int t = 0; t < timeSteps3; t++)
+                {
+                    p140208Array[t] = new float[latitudes3][];
+
+                    for (int lat = 0; lat < latitudes3; lat++)
+                    {
+                        p140208Array[t][lat] = new float[longitudes3];
+
+                        for (int lon = 0; lon < longitudes3; lon++)
+                        {
+                            p140208Array[t][lat][lon] = (float)dataArray3[t][lat][lon];
+                        }
+                    }
+                }
+                
+                windNCread.p140208Array = p140208Array;
+                Console.WriteLine("a value of p140208Array: " + p140208Array[2][3][9]);
+            }
+        }
+
+        string url4 = "http://10.130.100.5:7785//rescue-platform-service/api/v1/dem/getNcInfos?ncPath=D:/data/haiyang/DATA.nc&varName=mwd&centerLon=116.41992593821296&centerLat=40.18801994965735";
+
+        // Create a GET request to the specified URL
+        HttpWebRequest request4 = (HttpWebRequest)WebRequest.Create(url4);
+        request1.Method = "GET";
+
+        // Send request and get response
+        using (HttpWebResponse response4 = (HttpWebResponse)request4.GetResponse())
+        {
+            using (StreamReader reader4 = new StreamReader(response4.GetResponseStream()))
+            {
+                string result4 = reader4.ReadToEnd();
+
+                // Parse the JSON string into a JObject
+                JObject resultObject4 = JObject.Parse(result4);
+
+                // Extract the 'data' field from the parsed JObject
+                JArray dataArray4 = (JArray)resultObject4["data"];
+
+                // Determine the dimensions of the data array
+                int timeSteps4 = dataArray4.Count;
+                int latitudes4 = ((JArray)dataArray4[0]).Count;
+                int longitudes4 = ((JArray)dataArray4[0][0]).Count;
+
+                // Initialize float[][][] array
+                float[][][] mwdArray = new float[timeSteps4][][];
+
+                // Convert each element to float and populate v10Array
+                for (int t = 0; t < timeSteps4; t++)
+                {
+                    mwdArray[t] = new float[latitudes4][];
+
+                    for (int lat = 0; lat < latitudes4; lat++)
+                    {
+                        mwdArray[t][lat] = new float[longitudes4];
+
+                        for (int lon = 0; lon < longitudes4; lon++)
+                        {
+                            mwdArray[t][lat][lon] = (float)dataArray4[t][lat][lon];
+                        }
+                    }
+                }
+                
+                windNCread.mwdArray = mwdArray;
+                Console.WriteLine("a value of mwdArray: " + mwdArray[2][3][9]);
+                Console.ReadKey();
+            }
+        }
+    }
+}

SimulationServer/Entity/AircraftSJ.cs

@@ -1,6 +1,7 @@
 using KYFramework;
 using Model;
 using MuShiApp;
+using OpenCvSharp;
 using SimulationCommon;
 using SimulationSingleServer.Utils;
 using Unity.Mathematics;
@@ -26,21 +27,19 @@ public class AircraftSJ : AircraftEntity
     public override void Start()
     {
         //TODO 计算 AirRoute[]
-
-        float[] longitudeArray = null;//：经度一维数组；
-        float[] latitudeArray = null;//：纬度一维数组；
-        float[][][] u10Arrayu = null;// 10Array：风的10米U分量三维数组_24（时间）*7（纬度）*10（经度）；
-        
-        float[][][] v10Array = null;// v10Array：风的10米V分量三维数组_24（时间）*7（纬度）*10（经度）；
-        
-        float[][][] p140208Array = null;//：海洋上空的自由对流速度_24（时间）*7（纬度）*10（经度）；
-        float[][][] mwdArray = null;//：平均波向_24（时间）*7（纬度）*10（经度）；
+        double[] initialPosition = { 40.18801994965735, 116.41992593821296 };
+        double dt = 1;
+        double totalTime = 24.0;
         
-        WeatherResponse weatherResponse = new WeatherResponse() { };
+        Text_readNC text_ReadNC = new Text_readNC();
+        text_ReadNC.GetNCData();
+        var nCread = text_ReadNC.windNCread;
+        //漂移轨迹
+        List<double[]> trajectory = SeaSJ.CalculateDriftTrajectory(nCread, initialPosition,dt, totalTime);
+        Tuple<Point2f[], double> temp = SeaSJ.getminEnclosingRect(trajectory);
         
         
         
-        //