使用 linq 或亚音速计算距离

Question

我从搜索页面得到了这个 MySQL 语句，用户输入邮政编码，它会在输入的邮政编码 15 英里内找到最近的 stiocklist。

 SELECT *  , (
(
ACOS( SIN( "+SENTLNG +" * PI( ) /180 ) * SIN( s_lat * PI( ) /180 ) + COS( " + SENTLNG +" * PI( ) /180 ) * COS( s_lat * PI( ) /180 ) *  COS( (
" + SENTLANG + " - s_lng
) * PI( ) /180 ) ) *180 / PI( )
) *60 * 1.1515
) AS distance_miles
FROM new_stockists
WHERE s_lat IS NOT NULL
HAVING distance_miles <15
ORDER BY distance_miles ASC
LIMIT 0 , 15  

但现在我正在使用 linq 和 subsonic，但不知道如何在 linq 或 subsonic 中做到这一点 我们将非常感谢您的帮助，也请注意，我必须发送动态的发件人地址，即页面顶部提到的邮政编码，我会致电谷歌，从他们那里获取给定邮政编码的 lng 和 lat 。

Answer 1

我建议（如果可能）获取邮政编码的中心点，然后应用半正弦公式来查找特定半径内的所有商店。

这里的公式以公里为单位。
您必须更改相关数字，它将适用于英里。
例如：将 6371.392896 转换为英里。

DECLARE @radiusInKm AS FLOAT
将 @lat2Compare 声明为 FLOAT
将@long2Compare声明为FLOAT
设置@radiusInKm = 5.000
SET @lat2Compare = insert_your_lat_to_compare_here
SET @long2Compare = insert_you_long_to_compare_here

SELECT * FROM insert_your_table_here WITH(NOLOCK)
WHERE (6371.392896*2*ATN2(SQRT((sin((弧度(GeoLatitude - @lat2Compare)) / 2) * sin((弧度(GeoLatitude - @lat2Compare)) / 2)) + (cos(弧度(GeoLatitude)) * cos(弧度(@lat2Compare)) * sin(弧度(GeoLongitude - @long2Compare)/2) * sin(弧度(GeoLongitude - @long2Compare)/2)))
, SQRT(1-((sin((弧度(GeoLatitude - @lat2Compare)) / 2) * sin((弧度(GeoLatitude - @lat2Compare)) / 2)) + (cos(弧度(GeoLatitude)) * cos(弧度(@lat2Compare)) * sin(弧度(GeoLongitude - @long2Compare)/2) * sin(弧度(GeoLongitude - @long2Compare)/2)))
))) <= @radiusInKm

如果您想在 C# 中执行半正矢公式，则

double resultDistance = 0.0;
双avgRadiusOfEarth = 6371.392896； //地球半径不同，我取平均值。

//半正矢公式
//距离= R * 2 * aTan2（A的平方根，1 - A的平方根）
// 其中 A = 正弦平方（纬度差 / 2）+（纬度 1 的余弦 * 纬度 2 的余弦 * 正弦平方（经度差 / 2））
// 且 R = 地球周长

double DifferenceInLat = DegreeToRadian(currentLatitude - latitudeToCompare);
双差InLong = DegreeToRadian(currentLongitude - longtitudeToCompare);
double aInnerFormula = Math.Cos(DegreeToRadian(currentLatitude)) * Math.Cos(DegreeToRadian(latitudeToCompare)) * Math.Sin(differenceInLong / 2) * Math.Sin(differenceInLong / 2);
double aFormula = (Math.Sin((differenceInLat) / 2) * Math.Sin((differenceInLat) / 2)) + (aInnerFormula);
resultDistance = avgRadiusOfEarth * 2 * Math.Atan2(Math.Sqrt(aFormula), Math.Sqrt(1 - aFormula));

DegreesToRadian 是我自定义创建的函数。
这是一个简单的 1 行“Math.PI * 角度 / 180.0”

对于 LINQ，您可以使用 C# 数学函数以及所有 C# 检查。例如：!= 等于不等于等
请参阅以下示例。
它并不完整，因此请根据您的喜好进行调整。

var linqQuery = 来自 linqCollection 中
在此处插入您的收藏
其中 s_lat != 什么都没有
选择 Math.ACos(Math.Sin(DegreesToRadian(sentlng))
* Math.Pi / 180))

请查找下面的 MSDN 链接以获取所有简单的 LINQ 示例。尝试一下，希望这有帮助

我的博客条目 - SQL Haversine

MSDN - 101 LINQ 示例

Answer 2

使用 SQL 的这一部分在 SQL Server 中创建一个新视图：

SELECT * , (your equation here) as distance
FROM new_stocklists
WHERE s_lat is not NULL

然后您可以为您的视图创建一个 Linq 对象（SQLMetal 将执行此操作或 Visual Studio 中的 Linq to SQL）。然后您可以使用 Linq 来查询该视图。假设您的对象是 StockDistance：

var list = db.StockDistance.Where(x=>x.distance<15)
    .OrderBy(x=>x.distance)
    .Take(15);

Answer 3

您可以在 MS SQL 中创建一个存储过程来执行您的查询的操作，然后从您的应用程序调用该查询。 Linq 确实支持存储过程 - 有点像这样

partial class StockistsDataContext
{
    [Function(Name = "dbo.NewStockistsByDistance")]
    public ISingleResult<NewStockist> NewStockistsByDistance(
        [Parameter(DbType = "Int", Name = "s_lat")] int lat,
        [Parameter(DbType = "Int", Name = "s_lng")] int lng)
    {
        var result = this.ExecuteMethodCall(this, ((MethodInfo)(MethodInfo.GetCurrentMethod())), lat, lng);
        return ((ISingleResult<NewStockist>)(result.ReturnValue));
    }
}

回退到存储过程有点令人恼火，但我不认为 Linq 解析 Queryables 表达式树的部分可以检测相关数学函数并将它们映射到 SQL。