给定了经纬度的一张my_latlng表，和一个栅格my_grid表，怎么实现my_latlng表回填栅格id？

时间：2017-08-22 01:47:29 阅读：214 评论：0 收藏：0 [点我收藏+]

场景：

假设我们拥有一个拥有了一系列经纬度的表my_latlng(lat string,lng string)表，还有一张给定的栅格表my_grid(gridid bigint,centerlng double,centerlat double,gridx int,gridy int,minlng double,maxlng double,minlat double,maxlat double)并且栅格的为一个边长为５m的正方形，其中：

gridid   ：栅格id
centerlng：栅格中心点经度
centerlat ：栅格中心点纬度
gridx    ：栅格x轴方向的坐标位置
gridy    ：栅格y轴方向的坐标位置

需求：给my_latlng表找它所落的栅格的id，如果my_latlng中的经纬度，在总体栅格以外，就不参与运算。

解决方案一：

由于栅格有最大、最小经纬度，因此可以直接使用栅格的经纬度范围来给my_latlng表回填栅格id：

select t11.gridid,t10.lat,t10.lng
from my_latlng t10 inner join my_grid t11
where t10.lat>=t11.minlat and t10.lat<=t11.maxlat 
and t10.lng>=t11.minlng and t10.lng<t11.minlng;

缺陷：该种方案缺陷inner join 是没有on条件的，如果在ｈｉｖｅ中是没有办法把>=,>,<,<=符号给写到inner join 中 on条件上的，语法问题吧。

因此，上边这条语句是执行的一个cross join，如果my_latlng表有1000wt条记录，而my_grid有10000w条记录时，这样的一个cross join 在加上 where条件，就会导致这个数据在集群中1000spark套餐(假设说1spark套餐：１vcore cpu+12g memory+500g disk。)，5个小时也无法分析出结果。

解决方案二：

我们知道经纬度小数点第５位代表的基本就是米单位，纬度３０°时，经度每变化０.００００１相当于变化１.１m。因此，我们可以粗略的认为５ｍ的栅格在经度、纬度上的变化为０.００００５个单位的变化。

因此，利用上边的这个特性我们可以有以下方案：

第一步、可以找到某些栅格距离自己纬度和经度变化接近５ｍ的周围８＋１个左右的栅格；

(
rpad(t10.lat+0.00005,7,‘0‘)=rpad(t11.centerlat,7,‘0‘)  
or rpad(t10.lat,7,‘0‘)=rpad(t11.centerlat+0.00005,7,‘0‘) 
or rpad(t10.lat,7,‘0‘)=rpad(t11.centerlat,7,‘0‘)
)
and
(
rpad(t10.lng+0.00005,8,‘0‘)=rpad(t11.centerlng,8,‘0‘) 
or rpad(t10.lng,8,‘0‘)=rpad(t11.centerlng+0.00005,8,‘0‘) 
or rpad(t10.lng,8,‘0‘)=rpad(t11.centerlng,8,‘0‘)
)

备注：我们计算范围：经度范围１００.０到１８０.０，纬度范围为：１０.０到９０.０。

第二步、从第一步中的栅格中挑选距离自己最近的一个栅格作为自己归属栅格。

(
(cast(t10.lng as double)-t11.centerlng)*(cast(t10.lng as double)-t11.centerlng)
+(cast(t10.lat as double)-t11.centerlat)*(cast(t10.lat as double)-t11.centerlat)
) distans

但是上边的程序如果放在ｈｉｖｅ中的语句因该是这么写：

select t11.gridid,t10.lat,t10.lng,(
　　(cast(t10.lng as double)-t11.centerlng)*(cast(t10.lng as double)-t11.centerlng)
　　+(cast(t10.lat as double)-t11.centerlat)*(cast(t10.lat as double)-t11.centerlat)) distans
from my_latlng t10 inner join my_grid t11
where (
　　rpad(t10.lat+0.00005,7,‘0‘)=rpad(t11.centerlat,7,‘0‘)  
　　or rpad(t10.lat,7,‘0‘)=rpad(t11.centerlat+0.00005,7,‘0‘) 
　　or rpad(t10.lat,7,‘0‘)=rpad(t11.centerlat,7,‘0‘)
)
and(
　　rpad(t10.lng+0.00005,8,‘0‘)=rpad(t11.centerlng,8,‘0‘) 
　　or rpad(t10.lng,8,‘0‘)=rpad(t11.centerlng+0.00005,8,‘0‘) 
　　or rpad(t10.lng,8,‘0‘)=rpad(t11.centerlng,8,‘0‘)
);

但是上边的程序是有以下两个问题：

问题１）inner join 没有 on 条件，原因是ｗｈｅｒｅ中语句不允许写到ｏｎ中，也是ｈｉｖｅ的语法问题；

问题２）上边这段代码也是执行的cross join，而此自然执行的也很慢。

好的事情是：

１）从这段代码中，我们是可以把多个语句拆分出９个语句，而且这９个语句是可以具有ｏｎ条件的；