在两个位置之间的线上的不同位置处减去Lat \ Lng坐标

时间:2012-10-27 11:43:32

标签: c# sql geolocation sql-server-2008-r2 geometry

我有两组lat / lng坐标,在视觉上我画了一条线。 有没有办法在你走路时推断出纬度/经度坐标?沿着这条线

所以沿着线10%,lat / lng将是

沿着线的20%,lat / lng将是

沿着线的30%,纬度/经度

等..

我希望sql中的地理内容可能有一个简单的解决方案..

理想情况下,在SQL(SQL Server 2008 R2数据库)中,或者 - 如果太棘手 - 可能在C#

任何帮助将不胜感激

由于

2 个答案:

答案 0 :(得分:1)

C#中的这个库提供了许多可能有用的工具:

http://www.gavaghan.org/blog/free-source-code/geodesy-library-vincentys-formula/

EDIT 这使用地图投影,因此它不使用直线。您需要注意纬度之间的距离或多或少是恒定的,但经度之间的距离取决于纬度。

除了上面的库之外,还有一些更简单的近似值,用于测量两点之间的距离,以及给定距离处的点的坐标,并从起点开始。

这是一些测量两点之间距离的代码:

double lat1 = FSConvert.DegreesToRadians(start.Latitude.Decimal);
                        double lon1 = FSConvert.DegreesToRadians(start.Longitude.Decimal) * -1;
                        double lat2 = FSConvert.DegreesToRadians(end.Latitude.Decimal);
                        double lon2 = FSConvert.DegreesToRadians(end.Longitude.Decimal) * -1;
                        double result = Math.Acos(Math.Sin(lat1) * Math.Sin(lat2) + Math.Cos(lat1) * Math.Cos(lat2) * Math.Cos(lon1 - lon2));
                        distance = FSConvert.DradtoMeters(result);

这是一些可以解决问题的代码:

if (start == null || end == null)
                        {
                            return 0.0f;
                        }

                        double lat1 = FSConvert.DegreesToRadians(start.Latitude.Decimal);
                        double lon1 = FSConvert.DegreesToRadians(start.Longitude.Decimal) * -1;
                        double lat2 = FSConvert.DegreesToRadians(end.Latitude.Decimal);
                        double lon2 = FSConvert.DegreesToRadians(end.Longitude.Decimal) * -1;
                        double y = Math.Atan2(Math.Sin(lon1 - lon2) * Math.Cos(lat2), Math.Cos(lat1) * Math.Sin(lat2) - Math.Sin(lat1) * Math.Cos(lat2) * Math.Cos(lon1 - lon2));
                        const double x = 2 * Math.PI;
                        double result = y - x * Math.Floor(y / x);
                        bearing = FSConvert.RadiansToDegrees(result);

这是一些代码,用于根据距离和来自原点的方位计算终端坐标

double lat1 = FSConvert.DegreesToRadians(start.Latitude.Decimal);
                        double lon1 = FSConvert.DegreesToRadians(start.Longitude.Decimal) * -1;
                        double d = FSConvert.MetersToDrad(distance);
                        double tc = FSConvert.DegreesToRadians(bearing);

                        double lat = Math.Asin(Math.Sin(lat1) * Math.Cos(d) + Math.Cos(lat1) * Math.Sin(d) * Math.Cos(tc));
                        double lon = ((lon1 - Math.Asin(Math.Sin(tc) * Math.Sin(d) / Math.Cos(lat)) + Math.PI) % (2 * Math.PI)) - Math.PI;

                        var returnPoint = new FSPoint {
                            Latitude = {
                                Decimal = FSConvert.RadiansToDegrees(lat)
                            },
                            Longitude = {
                                Decimal = FSConvert.RadiansToDegrees(lon) * -1
                            }
                        };
                        point = returnPoint;

对FSConvert的引用只是改变RadiansToDegrees等等 - 这是微不足道的。 FSPoint只是一个lat / long结构。

所以过程是:

  1. 计算两点之间的距离和方位
  2. 除了
  3. 之外,将距离除以10 0r
  4. 使用增量距离和方位计算终端坐标。
  5. 坚持走线路#39;直到你到达另一端
  6. 如果您需要更多帮助,请告诉我

答案 1 :(得分:0)

在JAVA中,您可以使用此代码并且效果非常好(从上面的ScruffyDuck的回复转换)

double x2 = next.latitude; 
                double y2 = next.longitude;
                double x1 = current.latitude;
                double y1 = current.longitude;

                double lat1 = Math.toRadians(x1);
                double lon1 = Math.toRadians(y1) * -1;
                double lat2 = Math.toRadians(x2);
                double lon2 = Math.toRadians(y2) * -1;
                double distance = Math.acos(Math.sin(lat1) * Math.sin(lat2) + Math.cos(lat1) * Math.cos(lat2) * Math.cos(lon1 - lon2));

                double var1 = Math.atan2(Math.sin(lon1 - lon2) * Math.cos(lat2), Math.cos(lat1) * Math.sin(lat2) - Math.sin(lat1) * Math.cos(lat2) * Math.cos(lon1 - lon2));
                double var2 = 2 * Math.PI;
                double bearing = var1 - var2 * Math.floor(var1 / var2);

                double lat = Math.asin(Math.sin(lat1) * Math.cos(distance) + Math.cos(lat1) * Math.sin(distance) * Math.cos(bearing));
                double lon = ((lon1 - Math.asin(Math.sin(bearing) * Math.sin(distance) / Math.cos(lat)) + Math.PI) % (2 * Math.PI)) - Math.PI;

                double x = Math.toDegrees(lat);
                double y = Math.toDegrees(lon) * -1;
                //****************

                Location l = new Location("temp");
                l.setLatitude(x);
                l.setLongitude(y);