カラー出力してみる

奈良湖をずっと真っ黒で描いてましたが、あんまりな見た目なんで、ここは葦原をイメージした緑色＋深いところは青、にしてみましょう。

CREATE TABLE [destテーブル] (rid serial PRIMARY KEY, rast raster);
INSERT INTO [destテーブル](rast) 
 SELECT ST_AddBand(NULL,
        ARRAY[
         ST_Reclass(rast, 1,'[0.0-35.0]:0-0,    (35.0-45.0]:0-0,  (45.0-9999.0]:0', '8BUI', 0), 
         ST_Reclass(rast, 1,'[0.0-35.0]:0-0,    (35.0-45.0]:1-127,(45.0-9999.0]:0', '8BUI', 0), 
         ST_Reclass(rast, 1,'[0.0-35.0]:127-127,(35.0-45.0]:127-1,(45.0-9999.0]:0', '8BUI', 0), 
         ST_Reclass(rast, 1,'[0.0-45.0]:255-255,(45.0-9999.0]:0',                   '8BUI', 0) 
        ]
        ) FROM [srcテーブル];

複数バンドのラスタは、エクスポート時にはバンド1〜4が順にRGBAと判断されますので、これでカラーになります。例によって gdal_translate で出力すると、

アオコ大発生。