市区町村別データの作り方

探し方が悪いのかもだが、ArcGISのデータ以外で市区町村のshpファイルを作ろうとすると結構大変。 しかし、国土数値情報の細かい境界データをダウンロードしてきてQGISでWKTにしてbigqueryにぶち込んで、これで終わり。

select N03_001, N03_004, N03_007, st_astext(st_union_agg( ST_GeogFromText(WKT))) as WKT from boundaries.boundaries group by 1,2,3;```