public Tuple <Second[], Second[], MeterPerSecond[]> CheckPattern(List <Tuple <GeoLocationDecimalDegree, DateTime> > pings)
{
if (pings.Count < 2)
{
return(new Tuple <Second[], Second[], MeterPerSecond[]>(new Second[pings.Count], new Second[pings.Count], new MeterPerSecond[pings.Count]));
}
int ping_count = pings.Count;
GeoLocationDecimalDegree[] ship_ping_locations = new GeoLocationDecimalDegree[ping_count];
for (int ping_index = 0; ping_index < ping_count; ping_index++)
{
ship_ping_locations[ping_index] = pings[ping_index].Item1;
}
Second[] ping_seconds = new Second[ping_count];
MeterPerSecond[] ping_velocities = new MeterPerSecond[ping_count];
//double[] ping_areas = new double[ping_count];
Second[] ping_seconds_in_range = new Second[ping_count];
Meter[,] distance_matrix = ToolsMathFunction.FillArray <GeoLocationDecimalDegree, GeoLocationDecimalDegree, Meter>(new FunctionDistanceGeoLocation(), ship_ping_locations, ship_ping_locations);
Meter ping_distance = new Meter();
ping_seconds[0] = (Second)((pings[1].Item2 - pings[0].Item2).TotalMilliseconds) / 2000.0;
ping_distance = distance_matrix[0, 1] / 2.0;
ping_velocities[0] = ping_distance / ping_seconds[0];
for (int ping_index = 1; ping_index < ping_count - 1; ping_index++)
{
ping_seconds[ping_index] += new Second((pings[ping_index].Item2 - pings[ping_index - 1].Item2).TotalMilliseconds) / 2000.0;
ping_seconds[ping_index] += new Second((pings[ping_index + 1].Item2 - pings[ping_index].Item2).TotalMilliseconds) / 2000.0;
ping_distance = (distance_matrix[ping_index, ping_index - 1] + distance_matrix[ping_index, ping_index + 1]) / 2.0;
ping_velocities[ping_index] = ping_distance / ping_seconds[ping_index];
}
ping_seconds[ping_count - 1] = new Second((pings[ping_count - 1].Item2 - pings[ping_count - 2].Item2).TotalMilliseconds) / 2000.0;
ping_distance = distance_matrix[ping_count - 2, ping_count - 1] / 2.0;
ping_velocities[ping_count - 1] = ping_distance / ping_seconds[ping_count - 1];
for (int index_ping_0 = 0; index_ping_0 < ping_count; index_ping_0++)
{
for (int index_ping_1 = 0; index_ping_1 < ping_count; index_ping_1++)
{
if (distance_matrix[index_ping_0, index_ping_1] < SurveyRangeInMeters)
{
ping_seconds_in_range[index_ping_0] += ping_seconds[index_ping_1];
}
}
}
return(new Tuple <Second[], Second[], MeterPerSecond[]>(ping_seconds, ping_seconds_in_range, ping_velocities));
}
public static Meter ComputeDistanceHaversine(GeoLocationDecimalDegree location_0, GeoLocationDecimalDegree location_1)
{
return(ComputeDistanceHaversine(location_0.Latitude, location_0.Longitude, location_1.Latitude, location_1.Longitude));
}
public static Meter ComputeDistanceEquirectangular(GeoLocationDecimalDegree location_0, GeoLocationDecimalDegree location_1)
{
//Scuks
return(ComputeDistanceEquirectangular(location_0.Latitude, location_0.Longitude, location_1.Latitude, location_1.Longitude));
}
