1.安装软件环境
安装以下软件，并安装GeoServer imagemosaic拓展插件，用来连接并加载PostGIS数据库中的栅格数据。
PostGIS 2.4，PostgreSQL 9.6，GeoServer 2.12，geoserver-2.12-SNAPSHOT-imagemosaic-jdbc-plugin
2.导入栅格数据
进入PostgreSQLin目录，利用raster2pgsql工具导入栅格数据，具体命令如下所示：
<!-- 分块，切片存储到PostGIS数据库中 --> raster2pgsql -s 4326 -I -C -M E:/personalfile/cs1/alt.tif -F -t 256x256 public.alt | psql -h localhost -p 5432 -U postgres -d rastertest -W <!-- 以单个文件方式存储到PostGIS数据库中 --> raster2pgsql -s 4326 -I -C -M E:/personalfile/cs1/alt.tif public.alt | psql -h localhost -p 5432 -U postgres -d rastertest -W
其中各参数的含义如下：
-c 创建新表，并把栅格数据导入到表中，这是默认模式 -a 追加栅格数据到一个已经存在的表 -d 先drop掉表，然后创建新表，并把栅格数据导入到表中 -p 准备模式，只创建表，不导入数据. #栅格处理过程：把适当的约束条件注册（写入）栅格系统表中。 -C 用于栅格约束注册——比如栅格的srid, pixelsize 元数据信息等等，保证栅格在视图raster_columns能够合适地注册。 -x 让设置最大边界约束失效，只有在-C参数同时使用时候才能使用这个参数。 -r 为规则的块设置约束（空间唯一性约束和覆盖瓦片）。只有在同时使用-C 参数时，才能使用这个参数。 #栅格处理参数：用于操作输入栅格数据集的可选参数。 -s 表示<SRID>：指定输出栅格使用指定的SRID。如果没有提供或者提供值为0，程序会检查栅格的元数据（比如栅格 文件）以便决定一个合适的SRID。 -b 从栅格要抽取的波段位置（下标从1开始）。想要抽取不止1个波段的话，使用逗号，分隔波段数字。如果没有指 定，所有栅格的波段都会被抽取出来。 -t TILE_SIZE：瓦片大小。 把每一个表的行表示的栅格切割成瓦片。 TILE_SIZE被表示成WIDTHxHEIGHT的形式，或 者设置成值为“auto”来让加载程序使用第一个栅格参数计算合适的瓦片大小，然后把这个瓦片大小应用于所有的栅格 -R, --register 把栅格以文件数据形式注册（栅格数据保存在数据库之外的文件系统中）。只有栅格的元数据和文件路径 保存在数据库中（没有保存像素的数据） -l OVERVIEW_FACTOR ：创建栅格的概览信息。如果factor值不只一个，使用逗号，进行分隔。创建的概览数据存 储在数据库中，并且不受参数-R影响。注意你生成的SQL文件会包含存储数据的主表和概览表。 -N NODATA ：指定没有NODATA值的波段的NODATA值。 #用于操作数据库对象的可选参数 -q ：把PostgreSQL的id用引号包起来。 -f COLUMN :指定目标的栅格列名，默认是‘rast’。 -F ：创建一个以文件名为列名的列 -I ：在栅格列上创建一个GiST索引 -M Vacuum analyze栅格表 -T tablespace ：指定新表的表空间。注意如果索引（包括primary key）依然会使用默认的表空，除非使用-X参数标识。 -X tablespace ：指定表的索引表空间。如果使用了参数-I，那么这个参数可以应用于primary key和空间索引。 -Y ：使用copy语句而不是insert 语句。 -e ：单独地执行每一条语句，不使用一个大的事务执行所有语句。 -E ENDIAN ：控制栅格输出时候生成的二进制的字节顺序，如果值为0，表示取XDR，如果值为1，表示取NDR（默认取值）。 现在我们只支持NDR了。 -V version： 指定输出版本的格式，默认是0，当前值支持0。
创建mosaic数据表，并将已导入的栅格数据，插入到相应字段中，具体SQL语句如下：
/* 新建MOSAIC数据表 */ create table MOSAIC (NAME varchar(254) not null, TileTable varchar(254)not null, minX FLOAT8,minY FLOAT8, maxX FLOAT8, maxY FLOAT8,resX FLOAT8, resY FLOAT8, primary key (NAME,TileTable)); /* 插入alt栅格数据表 */ insert into MOSAIC(NAME,TileTable) values ('oek','alt');
导入成功后，如下图所示：
3.发布PostGIS栅格数据
准备connect.pgraster.xml.inc文件，配置PostgreSQL连接参数：
<connect> <dstype value="DBCP"/> <username value="postgres"/> <password value="123456"/> <jdbcUrl value="jdbc:postgresql://localhost:5432/rastertest"/> <driverClassName value="org.postgresql.Driver"/> <maxActive value="10"/> <maxIdle value="0"/> </connect>
准备mapping.pgraster.xml.inc文件，配置mosaic数据表参数：
<spatialExtension name="pgraster"/> <mapping> <masterTable name="mosaic" > <coverageNameAttribute name="name"/> <maxXAttribute name="maxX"/> <maxYAttribute name="maxY"/> <minXAttribute name="minX"/> <minYAttribute name="minY"/> <resXAttribute name="resX"/> <resYAttribute name="resY"/> <tileTableNameAtribute name="tiletable" /> </masterTable> <tileTable> <blobAttributeName name="rast" /> </tileTable> </mapping>
准备oek.pgraster.xml文件，配置栅格数据参数：
<?xml version="1.0" encoding="UTF-8" standalone="no"?> <!DOCTYPE ImageMosaicJDBCConfig [ <!ENTITY mapping PUBLIC "mapping" "mapping.pgraster.xml.inc"> <!ENTITY connect PUBLIC "connect" "connect.pgraster.xml.inc"> ]> <config version="1.0"> <coverageName name="oek"/> <coordsys name="EPSG:4326"/> <!-- interpolation 1 = nearest neighbour, 2 = bipolar, 3 = bicubic --> <scaleop interpolation="1"/> <axisOrder ignore="false"/> &mapping; &connect; </config>
将准备好的三个文件放在同一文件夹中，打开GeoServer，新建栅格数据源，如图所示：
发布栅格数据，参数配置及效果如下图所示：
