2013-04-16 43 views
2

鉴于我有一张表格,其中包含机场列表及其关联的纬度和长度坐标,我将如何着手获取距另一坐标100英里内的机场列表?我的最终目标是按城市搜索用户,并在距城市经纬度100英里范围内找到机场。我拥有的所有数据都不知道从哪里开始。在半径范围内计算100英里

我的数据在SQL Server中。

+0

多少英里是一个纬度和经度之间?这就是你需要找出 –

+0

你正在使用什么版本的SQL服务器?如果您拥有SQL Server 2008,则可以使用GeoCoordinates。 –

+0

http://stackoverflow.com/questions/862856/storing-and-querying-gps-coordinates-effectively也包含一些很好的信息。 –

回答

5
create function dbo.F_GREAT_CIRCLE_DISTANCE 
(
    @Latitude1 float, 
    @Longitude1 float, 
    @Latitude2 float, 
    @Longitude2 float 
) 
returns float as 
begin 
    declare @radius float 

    declare @lon1 float 
    declare @lon2 float 
    declare @lat1 float 
    declare @lat2 float 

    declare @a float 
    declare @distance float 

    -- Sets average radius of Earth in Miles 
    set @radius = 3956 

    -- Convert degrees to radians 
    set @lon1 = radians(@Longitude1) 
    set @lon2 = radians(@Longitude2) 
    set @lat1 = radians(@Latitude1) 
    set @lat2 = radians(@Latitude2) 

    set @a = sqrt(square(sin((@[email protected])/2.0E)) + (cos(@lat1) * cos(@lat2) * square(sin((@[email protected])/2.0E)))) 
    set @distance = @radius * (2.0E *asin(case when 1.0E < @a then 1.0E else @a end)) 

    return @distance 
end 
+0

+1这很好 – Brad

+0

我们可以在MySQL中做到这一点吗? – LoneWOLFs

+0

@Markus是Miles的结果吗? – LoneWOLFs

0

万一你的机场在0,0, 30英里范围内的所有机场应该在一个圆圈内。 测试如果点是一个圆内,则需要测试每个正方形方程: SQRT(X * X + Y * Y)< R, 其中X,Y - 各自增宽的位置,R =100英里

4

您可以使用System.Device.Location.GeoCoordinate

GeoCoordinate city = ....; 
double distance = .......; 
List<GeoCoordinate> airports = .......; //load from db. 

var found = airports.Where(c=>c.GetDistanceTo(city) < distance); 

编辑

对于那些谁可能有兴趣在实际工作代码

var airports1 = Airports.Search(51, 0, 100).ToList(); //~London 
var airports2 = Airports.Search(40.714623, -74.006605,100).ToList(); //~NY 

public class Airports 
{ 
    public class Airport 
    { 
     public string Name; 
     public GeoCoordinate Location; 

     public override string ToString() 
     { 
      return Name; 
     } 
    } 

    static Lazy<List<Airport>> _Airports = new Lazy<List<Airport>>(() => 
     { 
      using (var wc = new WebClient()) 
      { 
       var json = wc.DownloadString("http://www.flightradar24.com/AirportDataService2.php"); 
       var jObj = new JavaScriptSerializer().Deserialize<Dictionary<string,string[]>>(json); 

       return jObj.Values 
          .Select(j => new Airport 
          { 
           Name = (string)j[2], 
           Location = new GeoCoordinate(double.Parse(j[3].ToString(), CultureInfo.InvariantCulture),double.Parse(j[4].ToString(), CultureInfo.InvariantCulture)) 
          }) 
          .ToList(); 
      } 
     }, true); 

    public static IEnumerable<Airport> Search(double lat,double lon, double distInMiles) 
    { 
     var loc = new GeoCoordinate(lat, lon); 
     return _Airports.Value.Where(c => c.Location.GetDistanceTo(loc) < distInMiles * 1609.344); 
    } 
} 
+1

我认为OP正在T-SQL中寻求解决方案。 – Codeman

+0

@ Pheonixblade9可能是这样,但它已被C#标记。 –

1

你可以做这样的事情:

--SELECT Points closer than @radius 
select * from #Points p 
WHERE power(p.x - @locX, 2) + POWER(p.y - @locY, 2) < POWER(@radius,2) 

(这已被无耻地从http://www.sqlservercentral.com/Forums/Topic1228111-8-1.aspx#bm1229043撕开,通过谷歌搜索SQL + Pyhtagoras找到)。

+0

不幸的是,两条经度线之间的距离因纬度而异,所以这种解决方案在实践中根本不起作用。 –

+0

另外,不要忘记,地球表面的距离测量的是地球的曲线*的距离*,而不是两点间的距离*。如果你直接去中国,到中国的速度要快很多,但挖掘这条隧道将会是一些工作。 –

0

下面是“为乌鸦两点之间飞行距离下式: 半正矢

formula:  a = sin²(Δφ/2) + cos(φ1).cos(φ2).sin²(Δλ/2) 
c = 2.atan2(√a, √(1−a)) 
d = R.c 

其中φ是纬度,λ是经度,R是地球的半径(平均半径=6371公里) 音符角度需要以弧度传递以触发函数!

您可以将该代码转换成SQL

JavaScript: 
var R = 6371; // km 
var dLat = (lat2-lat1).toRad(); 
var dLon = (lon2-lon1).toRad(); 
var lat1 = lat1.toRad(); 
var lat2 = lat2.toRad(); 

var a = Math.sin(dLat/2) * Math.sin(dLat/2) + 
     Math.sin(dLon/2) * Math.sin(dLon/2) * Math.cos(lat1) * Math.cos(lat2); 
var c = 2 * Math.atan2(Math.sqrt(a), Math.sqrt(1-a)); 
var d = R * c; 

这里的T-SQL页面的链接数学函数: http://msdn.microsoft.com/en-us/library/ms188919.aspx