[已解决]代码提速问题！麻烦下(过程已详细补充清楚)

suxiong · 发表于 2014-1-7 16:08

本帖最后由 suxiong 于 2014-1-8 09:35 编辑

QQ截图20140108093341.jpg

复件 ne.zip (37.68 KB, 下载次数: 5)

香川群子 · 发表于 2014-1-7 16:37

呵呵，不知道你的代码有什么用处，能否简单介绍一下……

计算的话，规则什么的也应该说一下。

香川群子 · 发表于 2014-1-7 16:38

看数据像是经纬度距离计算。

suxiong · 发表于 2014-1-8 09:35

香川群子发表于 2014-1-7 16:37
呵呵，不知道你的代码有什么用处，能否简单介绍一下……

已补充了！麻烦看下！

suxiong · 发表于 2014-1-8 09:36

香川群子发表于 2014-1-7 16:38
看数据像是经纬度距离计算。

嗯，是的！

suxiong · 发表于 2014-1-8 09:39

香川群子发表于 2014-1-7 16:38
看数据像是经纬度距离计算。

嗯，就是哈！只是计算最有点大。

计算量是第一列*第三列；后面还要找最小值。

香川群子 · 发表于 2014-1-8 13:43

如果只需要最终计算结果，可以忽略中间计算过程的话，算法可以改进：

1. 简化计算得到近似值 (可以大大提高计算速度)
2. 求得对于A点的最短距离（近似值）
3. 对于和近似值接近的小部分B点值进行精度值计算，并从中得到最短距离点

另外，对B列点事先排序，也可大大提高计算效率。
（用二分法逼近中点）

香川群子 · 发表于 2014-1-8 14:07

提高效率的具体算法改进：

1. 把经纬度值预处理转为长整型数值以便快速计算（比双精度计算快很多）
2. 仅需计算长整形数值的经度差、维度差的平方的和（近似直角斜边）
3. 比较该值的大小，取其中最小的100个或者20个值的位置
4. 对这些筛选过的位置进行精确计算即可。

香川群子 · 发表于 2014-1-8 14:15

你再上一个12000行的数据附件吧。

现在的附件只有50多行用不了多少时间没有意义。

香川群子 · 发表于 2014-1-8 14:54

已知经纬度两点距离的公式
两点距离=6378137*2*ASIN(SQRT(SUMSQ(SIN((RADIANS(起点纬度) - RADIANS(终点纬度))/2))+ COS(RADIANS(起点纬度))*COS( RADIANS(终点纬度))*SUMSQ(SIN((RADIANS(起点经度)-RADIANS(终点经度))/2))))

或者正确的弧度公式是：
C = sin(lat1) * sin(lat2) + cos(lat1) * cos(lat2) * cos(lon2-lon1)

		自动登录	找回密码
密码			注册