• 周一. 8月 15th, 2022

5G编程聚合网

5G时代下一个聚合的编程学习网

热门标签

GeoServer发布PostGIS数据库中的栅格数据

admin

11月 28, 2021

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,新建栅格数据源,如图所示:

 发布栅格数据,参数配置及效果如下图所示:

发表回复

您的电子邮箱地址不会被公开。