PostgreSQL计算影响某地的台风

已爬取了十几年来的全国台风数据&＃xff0c;现要获取影响某地的台风&＃xff0c;想法比较粗糙&＃xff0c;就是以某地的中心点为半径和台风路径中心点距离的求解&＃xff0c;一般某地区是以点带面的有个地区半径&＃xff0c;台风一般取七级台风风圈半径。

现在的问题点转换为了&＃xff0c;两点之间求距离。

Postgresql有两种办法可实现&＃xff0c;一种是自定义函数&＃xff0c;一种是利用Postgresql扩展包

第一种自定义函数的方式&＃xff0c;参见代码&＃xff1a;

create or replace function getdistance
(ilngbegin real,ilatbegin real,ilngend real,ilatend real
)
returns float
as
$body$declarevdistance real;vearthradius real;vradlatbegin real;vradlatend real;vradlatdiff real;vradlngdiff real;
begin--地球半径vearthradius:&＃61;6378.137;-经纬度转换vradlatbegin :&＃61; ilatbegin * pi()/180.0;vradlatend :&＃61; ilatend * pi()/180.0;vradlatdiff :&＃61; vradlatbegin - vradlatend;vradlngdiff :&＃61; ilngbegin * pi()/180.0 - ilngend * pi()/180.0;--求距离vdistance :&＃61; 2 * asin(sqrt(power(sin(vradlatdiff / 2), 2) &＃43; cos(vradlatbegin) * cos(vradlatend) * power(sin(vradlngdiff/2),2)));vdistance :&＃61; vdistance * vearthradius*1000;return vdistance;
end;
$body$
language &＃39;plpgsql&＃39; volatile;

Postgresql内置了地理计算的一些功能&＃xff0c;需要在应用时赋权即可。

CREATE EXTENSION cube;
CREATE EXTENSION earthdistance;

这两种方法的调用如下&＃xff1a;

select a.tfid,a.pointlat,pointlng,getdistance(a.pointlat,a.pointlng,19.137154,109.876876) as distance,a.* from typhoon_his_point_info awhere getdistance(a.pointlat,pointlng,19.137154,109.876876)<&＃61;480000  order  by a.tfid,a.pointtime;select a.tfid,a.pointlat,pointlng,earth_distance(ll_to_earth (109.876876, 19.137154),ll_to_earth(a.pointlng, a.pointlat)) as distance,a.*from typhoon_his_point_info awhere getdistance(a.pointlat,pointlng,19.137154,109.876876)<&＃61;480000  order  by a.tfid,a.pointtime;

结论&＃xff1a;两种计算方法的误差在十万分之一&＃xff0c;当然直接使用Postgresql的会更精准一些。

欢迎关注python与大数据分析