最近有一批分幅的影像数据，每幅大概100-200MB，之前都是以散文件的形式直接放在磁盘里面。然后用QGIS加载，保存一个工程文件。同事说需要那个区域的影像的时候，我就打开QGIS，拉框选择，然后查看文件名，再根据文件名找到数据，然后发给同事。
时间久了，觉得这种方法比较low，也比较低效，人工操作的内容比较多，不好写到脚本里面去。正好这段时间手头没什么事情，就想想看有什么好办法，将影像管理起来，操作更加自动化。

最开始想到的是，自己编写一个脚本，遍历文件夹中所有影像，使用GDAL读影像的范围，然后把影像的名字，路径，范围等信息写在一张数据表里面。以后需要哪个范围的数据，直接一个SQL语句，返回区域内影像的路径。然后按照路径打包数据，发给同事。
但是后来发现一些问题：1. GDAL不太熟，而且不同类型的影像数据投影还不一样，空间范围要重投影，2. 既然都用到数据库来存储代表了影像空间范围的几何体，为什么不直接用数据库来存影像啊。3. 将数据发给同事的过程不好自动化，如果同事要的不是整张影像而是一部分的话，还需要进行裁剪。想到这些，就开始在网上查找可以存储影像的数据库。

作为一个GISer，对于PostgreSQL（与PostGIS）还是不陌生的。就看看PG的raster模块应该可以一试。PostGIS中支持使用raster类型来表示栅格数据。可以使用raster2pgsql来进行导入，支持的数据格式为GDAL支持的数据格式。

oracle 导入数据暂时还没找着现成的工具，就使用sql plus和sql进行导入了。网上搜了下，代码如下所示。

--创建一个自增的ID序列
CREATE SEQUENCE land_image_id_seq
    increment by 1
    start with 1
    nomaxvalue
    minvalue 1
    NOCYCLE;

--创建GeoRaster表（GeoRaster表就是带有SDO_GEORASTER类型列的表）
CREATE TABLE land(
    id NUMBER PRIMARY KEY,
    name VARCHAR(30),
    image SDO_GEORASTER
);

--创建一个自增的触发器
CREATE OR REPLACE TRIGGER land_image_id_trig
    BEFORE INSERT ON land
    FOR EACH ROW
    BEGIN
        SELECT land_image_id_seq.nextval INTO :new.id FROM dual;
    END;
/

--创建GeoRaster Data Table（RDT）表用来存储GeoRaste影像数据
CREATE TABLE rdt_land OF SDO_RASTER
(
    PRIMARY KEY(RASTERID, PYRAMIDLEVEL, BANDBLOCKNUMBER, ROWBLOCKNUMBER, COLUMNBLOCKNUMBER)
)
LOB(RASTERBLOCK) STORE AS (NOCACHE NOLOGGING);

--插入GeoRaster表一条记录（初始化GeoRaster对象）
INSERT INTO land (name, image) VALUES ('0', SDO_GEOR.INIT('rdt_land'));
call dbms_java.grant_permission('MDSYS','SYS:java.io.FilePermission', '/opt/test_data/land0.tif', 'read');
-- 导入Raster Data Table表中
DECLARE
    img SDO_GEORASTER;
BEGIN
    --Select the georater column
    SELECT image INTO img FROM land WHERE name='0' FOR UPDATE;
    --Import into the georater object
    SDO_GEOR.IMPORTFROM(img, 'blocksize=(256, 256)', 'TIFF', 'file', '/opt/test_data/land0.tif');
    --update the column
    UPDATE land SET image=img WHERE name='0';

    COMMIT;
END;
/

虽然网上资料比较多，但是自己操作起来，过程还是比较复杂的，得先创建栅格表和块表，然后获取数据的读权限，再进行导入。调试了很久才能成功导入数据，对于接触oracle影像模块的我来说，几乎每一步操作都会卡壳，还好从网上都能找到解决方案，总体来说，使用起来感觉还是不直观，需要书写大段的SQL语句，查询较多的参数配置。而且对于数据量大于2G的影像，总是导入失败。搜了下，有可能使用sql语句不能直接导入2G以上的影像。。鉴于如此多不便之处，暂时不采用这个方案了吧。

翻阅网上的资料可以使用 raster2pgsql来进行导入，大概用法如下所示：

The raster2pgsql is a raster loader executable that loads GDAL supported raster formats into sql suitable for loading into a PostGIS raster table. It is capable of loading folders of raster files as well as creating overviews of rasters.

后面还有10多个参数说明，太长了，就不列在这一篇博客里了。感觉功能很强大，但是有点复杂，不利于上手。

EXAMPLE USAGE:
raster2pgsql raster_options_go_here raster_file someschema.sometable > out.sql

Example：
raster2pgsql -s 4326 -I -C -M *.tif -F -t 100x100 public.demelevation > elev.sql
psql -d gisdb -f elev.sql

导入多个图像
raster2pgsql -s 4326 -C -l 2,4 -I -F -t 2700x2700 gray_50m_partial*.tif data_import.gray_50m_partial | psql -h localhost -p 5434 -U postgres -d mastering_postgis

运行raster2pgsql，可以生成sql语句，SQL语句大概如下所示：

1 BEGIN;
2 CREATE TABLE "public"."demelevation" ("rid" serial PRIMARY KEY,"rast" raster,"filename" text);
3 INSERT INTO "public"."demelevation" ("rast","filename") VALUES ('0100000300F..省略若干..223E4C32224056'::raster,'land.tif');
...中间省略若干行...
10587 CREATE INDEX ON "public"."demelevation" USING gist (st_convexhull("rast"));
10588 ANALYZE "public"."demelevation";
10589 SELECT AddRasterConstraints('public','demelevation','rast',TRUE,TRUE,TRUE,TRUE,TRUE,TRUE,FALSE,TRUE,TRUE,TRUE,TRUE,TRUE);
10590 END;
10591 VACUUM ANALYZE "public"."demelevation";

一个300MB的数据，生成的导入sql语句大概600MB，10000多行。看来要合理设置分片大小。

目前看来oracle各种操作有点复杂，需要书写大段的SQL语句。使用postgis需要使用raster2sql语句来生成sql，感觉不是很直观，如果不进行切分的话，存储不了过大的数据。Ganos整体来说，使用起来比较简便，一条SQL语句就可以导入影像，文档是中文的，也比较通俗易懂，所有操作都在云上进行，维护应该比本地简单，也符合目前云计算的发展趋势。
后续准备对Ganos for PG进行进一步的评估，从功能和性价比方面看是否合适目前的需求。

数据是导入进去了，但是能不能直接从数据库查看这些导入的数据呢？从论坛里面看到好像可以使用WMS等影像服务来查看，具体不知道怎么做，以后等有这个需求再去深入问问。