余計な部分をすてる

大きいラスタの操作はけっこう時間かかります。ラスタをインポートするとき、普通に入れると元ファイルの図郭が丸々入りますが、インポート後にいろいろ調べ物をするには必要最小限のラスタタイルだけに絞って扱ったほうが何かと便利です。邪魔な生駒山西側の大阪湾部分とか木津川周辺の低地はすてててしまいましょう。てきとうな矩形以外はさくっと削除。

DELETE FROM [テーブル]
    WHERE NOT ST_Intersects(ST_GeomFromText('POLYGON*1', 4326), rast);

前回紹介したピクセル単位の操作は若干ややこしかったですが、単にラスタ単位で行選択するだけだったらこんなふうに割と簡単に書けますし、時間もかかりません。ただこれは、当然行単位（タイル単位）のDELETE操作ですので、タイル単位でしか消せません。

*1: 135.6767 34.68, 135.6767 34.51, 135.8400 34.51, 135.8400 34.68, 135.6767 34.68