Rasterテーブルを手動でつくってみる

インポートしたままだとピクセル値が32bit floatで扱いづらい。そんな精度もいらないので 8bit uint（メートル単位）のテーブル用意します。

CREATE TABLE [8BUIテーブル名] (rid serial PRIMARY KEY, rast raster);

簡単すぎて説明することがないｗ。PostGISの既存ユーザーさんは、AddRasterColumn()的な関数を使ってないことに気づくかも。
ここに値を入れるわけですが、今回は湖をつくるのが目的で、山のほうのデータはいらないので、標高50mより上は「データなし」にしてみます*1。

INSERT INTO [8BUIテーブル名](rast)
 SELECT ST_Reclass(rast, 1,'[0.0-50.0]:0-50, (50.0-9999.0]:255', '8BUI', 255) rast
   FROM [元テーブル名];

every valley〜 shall be exalted ♪
不等号の書き方が微妙なのは仕様です。

*1:ST_AddBandでバンド追加してupdateでもいいんだが、他のツールでの扱いやすさを考えて＆説明のストーリ上、別テーブルでやってみた