1.安裝軟件環(huán)境

安裝以下軟件，并安裝GeoServer imagemosaic拓展插件，用來連接并加載PostGIS數(shù)據(jù)庫(kù)中的柵格數(shù)據(jù)。

PostGIS 2.4，PostgreSQL 9.6，GeoServer 2.12，geoserver-2.12-SNAPSHOT-imagemosaic-jdbc-plugin

2.導(dǎo)入柵格數(shù)據(jù)

進(jìn)入PostgreSQL\bin目錄，利用raster2pgsql工具導(dǎo)入柵格數(shù)據(jù)，具體命令如下所示：

<!-- 分塊，切片存儲(chǔ)到PostGIS數(shù)據(jù)庫(kù)中 -->
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

<!-- 以單個(gè)文件方式存儲(chǔ)到PostGIS數(shù)據(jù)庫(kù)中 -->
raster2pgsql -s 4326 -I -C -M E:/personalfile/cs1/alt.tif public.alt | psql -h localhost -p 5432 -U postgres -d rastertest -W

其中各參數(shù)的含義如下：

-c 創(chuàng)建新表，并把柵格數(shù)據(jù)導(dǎo)入到表中，這是默認(rèn)模式
-a 追加?xùn)鸥駭?shù)據(jù)到一個(gè)已經(jīng)存在的表
-d 先drop掉表，然后創(chuàng)建新表，并把柵格數(shù)據(jù)導(dǎo)入到表中
-p 準(zhǔn)備模式，只創(chuàng)建表，不導(dǎo)入數(shù)據(jù).
#柵格處理過程：把適當(dāng)?shù)募s束條件注冊(cè)（寫入）柵格系統(tǒng)表中。
-C 用于柵格約束注冊(cè)——比如柵格的srid, pixelsize 元數(shù)據(jù)信息等等，保證柵格在視圖raster_columns能夠合適地注冊(cè)。
-x 讓設(shè)置最大邊界約束失效，只有在-C參數(shù)同時(shí)使用時(shí)候才能使用這個(gè)參數(shù)。
-r 為規(guī)則的塊設(shè)置約束（空間唯一性約束和覆蓋瓦片）。只有在同時(shí)使用-C 參數(shù)時(shí)，才能使用這個(gè)參數(shù)。
#柵格處理參數(shù)：用于操作輸入柵格數(shù)據(jù)集的可選參數(shù)。
-s 表示<SRID>：指定輸出柵格使用指定的SRID。如果沒有提供或者提供值為0，程序會(huì)檢查柵格的元數(shù)據(jù)（比如柵格
文件）以便決定一個(gè)合適的SRID。
-b 從柵格要抽取的波段位置（下標(biāo)從1開始）。想要抽取不止1個(gè)波段的話，使用逗號(hào)，分隔波段數(shù)字。如果沒有指
定，所有柵格的波段都會(huì)被抽取出來。
-t TILE_SIZE：瓦片大小。 把每一個(gè)表的行表示的柵格切割成瓦片。 TILE_SIZE被表示成WIDTHxHEIGHT的形式，或
者設(shè)置成值為“auto”來讓加載程序使用第一個(gè)柵格參數(shù)計(jì)算合適的瓦片大小，然后把這個(gè)瓦片大小應(yīng)用于所有的柵格
-R, --register 把柵格以文件數(shù)據(jù)形式注冊(cè)（柵格數(shù)據(jù)保存在數(shù)據(jù)庫(kù)之外的文件系統(tǒng)中）。只有柵格的元數(shù)據(jù)和文件路徑
保存在數(shù)據(jù)庫(kù)中（沒有保存像素的數(shù)據(jù)）
-l OVERVIEW_FACTOR ：創(chuàng)建柵格的概覽信息。如果factor值不只一個(gè)，使用逗號(hào)，進(jìn)行分隔。創(chuàng)建的概覽數(shù)據(jù)存
儲(chǔ)在數(shù)據(jù)庫(kù)中，并且不受參數(shù)-R影響。注意你生成的SQL文件會(huì)包含存儲(chǔ)數(shù)據(jù)的主表和概覽表。
-N NODATA ：指定沒有NODATA值的波段的NODATA值。
#用于操作數(shù)據(jù)庫(kù)對(duì)象的可選參數(shù)
-q ：把PostgreSQL的id用引號(hào)包起來。
-f COLUMN :指定目標(biāo)的柵格列名，默認(rèn)是‘rast’。
-F ：創(chuàng)建一個(gè)以文件名為列名的列
-I ：在柵格列上創(chuàng)建一個(gè)GiST索引
-M Vacuum analyze柵格表
-T tablespace ：指定新表的表空間。注意如果索引（包括primary key）依然會(huì)使用默認(rèn)的表空，除非使用-X參數(shù)標(biāo)識(shí)。
-X tablespace ：指定表的索引表空間。如果使用了參數(shù)-I，那么這個(gè)參數(shù)可以應(yīng)用于primary key和空間索引。
-Y ：使用copy語句而不是insert 語句。
-e ：?jiǎn)为?dú)地執(zhí)行每一條語句，不使用一個(gè)大的事務(wù)執(zhí)行所有語句。
-E ENDIAN ：控制柵格輸出時(shí)候生成的二進(jìn)制的字節(jié)順序，如果值為0，表示取XDR，如果值為1，表示取NDR（默認(rèn)取值）。
現(xiàn)在我們只支持NDR了。
-V version： 指定輸出版本的格式，默認(rèn)是0，當(dāng)前值支持0。

創(chuàng)建mosaic數(shù)據(jù)表，并將已導(dǎo)入的柵格數(shù)據(jù)，插入到相應(yīng)字段中，具體SQL語句如下：

/* 新建MOSAIC數(shù)據(jù)表 */
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柵格數(shù)據(jù)表 */
insert into MOSAIC(NAME,TileTable) values ('oek','alt');

導(dǎo)入成功后，如下圖所示：

柵格數(shù)據(jù)表

3.發(fā)布PostGIS柵格數(shù)據(jù)

準(zhǔn)備connect.pgraster.xml.inc文件，配置PostgreSQL連接參數(shù)：

<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>

準(zhǔn)備mapping.pgraster.xml.inc文件，配置mosaic數(shù)據(jù)表參數(shù)：

<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>

準(zhǔn)備oek.pgraster.xml文件，配置柵格數(shù)據(jù)參數(shù)：

<?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>

將準(zhǔn)備好的三個(gè)文件放在同一文件夾中，打開GeoServer，新建柵格數(shù)據(jù)源，如圖所示：

柵格數(shù)據(jù)源

配置連接參數(shù)

發(fā)布柵格數(shù)據(jù)，參數(shù)配置及效果如下圖所示：

發(fā)布參數(shù)配置

柵格服務(wù)效果