2012年6月7日木曜日

PostGIS:流域界と土地利用の重なり合う部分

流域界を水域系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 件のコメント:

コメントを投稿