属性で ベクタデータ をサブセットする

昨日夜10時
「明日までにxxフォルダ内の等高線ベクタデータから100m間隔のデータを作っといて」
とさっくり依頼された。
数百ファイルにわかれて30GBくらいあるんだけど・・・

finds.jpの WMSじゃダメ? と一応聞いてみたが
ローカルに持っておきたいとのこと。

しょうがないので以下の内容(作者さまに感謝)をちょっと変えて
3の倍数でアホになるアレで以下のような処理を行った。
https://code.google.com/p/geospatialpython/source/browse/trunk/subset.py

あ、でも OGR を使えばもっと簡単だったのかも。