鈍足ランナーのIT日記

走るのが好きな5流のITエンジニアのブログ。

趣味の範囲は広いけど、どれも中途半端なクソブロガー楽しめるWebアプリを作ってあっと言わせたい。サーバーサイドPerl(Mojolicious)、クライアントサイドVue.js。Arduinoにも触手を伸ばす予定。

三島市のコインランドリーのボロネイ図作り、人口をカウントする

コインランドリー用のテーブルを用意しておく

前に述べた通り。QGISへ取り込む。

ボロネイ図を描写する

入力にコインランドリーのベクターを指定する。
そうするとシェープファイルが吐き出される。(volonay.shp)
そのシェープファイルを入力に新たなテーブルを作る

shp2pgsql  -s 4326 volonoy.shp volonoy > volonoy.sql
psql -d test -f volonoy.sql

人口カウント用のSQLを発行する

select  v.id, (select name from volonoy b where b.id = v.id), sum(tblt000609)  from volonoy v
 inner join  population   p  
on st_intersects( p.geom, 
v.geom ) group by v.id order by sum desc

f:id:kechiya:20160928001929p:plain
f:id:kechiya:20160928001944p:plain