GeoServer釋出PostGIS資料庫中的柵格資料

 

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，新建柵格資料來源，如圖所示：

 

柵格資料來源

 

配置連線引數

 

釋出柵格資料，引數配置及效果如下圖所示：

 

釋出引數配置