1.安裝軟件環境
安裝以下軟件，並安裝GeoServer imagemosaic拓展插件，用來連接並加載PostGIS數據庫中的柵格數據。
PostGIS 2.4，PostgreSQL 9.6，GeoServer 2.12，geoserver-2.12-SNAPSHOT-imagemosaic-jdbc-plugin
2.導入柵格數據
進入PostgreSQL\bin目錄，利用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，新建柵格數據源，如圖所示：
發布柵格數據，參數配置及效果如下圖所示：
