流域界を水域系CDでディゾルブしたものと土地利用コードメッシュをインターセクトする。
①水系域CDを取得。
SELECT distinct(水系域CD) FROM 水系域ディゾルブテーブル;
②インターセクト
SELECT a2.*,(ST_Intersection(a1.the_geom,a2.the_geom)) As the_geom
FROM 水系域ディゾルブテーブル a1 ,土地利用コード a2
WHERE ST_Intersects(a1.the_geom, a2.the_geom) and a1.W12_002 ='水系域CD';
これで指定した水系域単位ごとに土地利用コードが取得できるハズ!
QGISで見ると流域単位で土地利用コードが集計されてます。

0 件のコメント:
コメントを投稿