PostGIS批量导入栅格数据

PostGIS批量导入栅格数据,第1张

概述如果我们有一批以文件存储的影像数据如何利用PostGIS批量的导出到PostgreSQL数据库中进行管理呢? (单个数据的导入参见我上篇博文:PostGIS导入导出栅格数据) 我的实验环境如下: OS: Ubuntu 16.04 LTS PostgreSQL:9.5.5 PostGIS: 2.2 我的实验数据是全球气象数据(不过,我只下载了中国部分19区的平均气温数据(TIF格式)做示例演示) 下

如果我们有一批以文件存储的影像数据如何利用PostGIS批量的导出到Postgresql数据库中进行管理呢?
(单个数据的导入参见我上篇博文:PostGIS导入导出栅格数据)

我的实验环境如下:
OS: Ubuntu 16.04 LTS
Postgresql:9.5.5
PostGIS: 2.2

我的实验数据是全球气象数据(不过,我只下载了中国部分19区的平均气温数据(TIF格式)做示例演示)
下载网址:http://www.worldclim.org/tiles.php
下载好的数据进行解压可以看到包含了12个月的平均气温数据:

使用GDAL的命令行工具gdalinfo查看其元数据,可以看到其采用的是wgs84坐标系,大小是3600×3600。

我们可以在raster2pgsql命令中使用通配符批量的导入数据到Postgresql数据库,命令如下:
(具体参数可自行查看命令的帮助文档)

raster2pgsql -s 4326 -I -C -M ./*.tif -F -t 256x256 staging.tmean_19 | psql -h localhost -p 5432 -U postgres -d postgis_in_action -W

需要说明的是:
-t参数指定了数据的分块大小,原始的3600×3600以256×256的Tile进行划分存储。
-F参数在关系表中增加名为filename的列用于保存原始文件的名称,这在我们进行Tile的合并的时候很有用。

进入数据库可以查看我们导入的数据表

通过sql进行查看,发现每个文件被分割成了225块。(但是(3600*3600)/(256*256)应该是197,不是225,我也不知道为什么被分割成了225块???)

最后,说说怎么导出数据到文件。和上一篇博文一样,我们还是使用Python脚本进行导出。这里主要用到了两个函数ST_AsTIFF和ST_Union。ST_AsTIFF上一篇中我们已经使用过,用于将Postgresql中的raster类型导出为TIFF格式对应的bytea类型,而ST_Union对于栅格数据来说,用于将分块的影像重新拼接为完整的一副影像。
详细代码如下:

# -*- Coding: utf-8 -*-import psycopg2# Connect to an existing databaseconn = psycopg2.connect('host=localhost port=5432 user=postgres password=password dbname=postgis_in_action')# Open a cursor to perform database operationscur = conn.cursor()# Execute sql query# cur.execute("SELECT ST_AsTIFF(rast,'LZW') AS rasttiff FROM staging.wsIEarth WHERE rID=1;")cur.execute("SELECT ST_AsTIFF(ST_Union(rast),'LZW') AS rasttiff FROM staging.tmean_19 WHERE filename='tmean1_19.tif';")# Fetch data as Python objectsrasttiff = cur.fetchone()# Write data to fileif rasttiff is not None:    open('/home/theone/Desktop/tmean1.tif','wb').write(str(rasttiff[0]))# Close communication with the databasecur.close()conn.close()

导出的结果可以在QGIS中查看和原来的影像一模一样:

总结

以上是内存溢出为你收集整理的PostGIS批量导入栅格数据全部内容,希望文章能够帮你解决PostGIS批量导入栅格数据所遇到的程序开发问题。

如果觉得内存溢出网站内容还不错,欢迎将内存溢出网站推荐给程序员好友。

欢迎分享,转载请注明来源:内存溢出

原文地址: http://outofmemory.cn/sjk/1182025.html

(0)
打赏 微信扫一扫 微信扫一扫 支付宝扫一扫 支付宝扫一扫
上一篇 2022-06-02
下一篇 2022-06-02

发表评论

登录后才能评论

评论列表(0条)

保存