postgis实现点数据热力图

参考

https://blog.csdn.net/honor_zhang/article/details/95066077

插入测试数据

select * from gisdata.tbl_pos;



insert into gisdata.tbl_pos values ( random()*100000, md5(random()::text), random()*1000, 
point((random()*10000::int), (random()*10000::int)) );  
   
  
  select   
  width_bucket(pos[0], 0, 10001, 100) p1,  -- x轴落在哪列bucket  
  width_bucket(pos[1], 0, 10001, 100) p2,  -- y轴落在哪列bucket  
  avg(pos[0]),  avg(pos[1]) ,
  min(pos[0]),  
  max(pos[0]),  
  count(*)  
from gisdata.tbl_pos  
  group by p1,p2; 


cluster  kmeans;;