创建pgraster的overview表以及overviewconstrants表约束

创建pgraster的overview表以及overviewconstrants表约束

练习postgis官方文档中有关pgraster部分的内容。记录下笔记。

在pgraster导入的时候raster2pgsql这个插件使用的时候，忘记使用-l这个参数来建立overview

这个时候需要手动在写sql添加

1.创建一个overview表

 

CREATE TABLE temp_overviewtable AS SELECT
ST_AddBand(
  ST_MakeEmptyRaster(1024, 1024, 0, 0, 4),
  1, '8BSI'::text, -129, NULL
) r2;

 

2.创建与raster表字段对应关系，使用函数AddOverviewConstraints来实现

 

select addoverviewconstraints('temprature','rast','temp_overviewtable','r2','6')

 

3.创建完成，查询语句查看：

select * from raster_columns
select * from raster_overviews

select st_astiff(st_union(rast)) as rasttiff from public.temprature where filename='wc2.0_2.5m_tmax_01.tif' limit 10

4.关于导出pgraster查看效果：有输入了，当然想看一下输出是很什么样的，而postgis恰恰是去图形化的，压根没有图形可视化的窗口，所以想看效果只能通过写代码看。

怎么写代码，些什么代码？从官网上看有php，aps.net,java的，还有各plpython的方式的这个要求的安装plpython关于这种方式自己去官网试试吧。

5.我选了csdn上的一篇博文：https://blog.csdn.net/theonegis/article/details/55100365

来试，总是报错TypeError: a bytes-like object is required, not 'str'。修改了其中的输出的代码。见下面代码的粉色部分

---

import psycopg2
conn=psycopg2.connect('host=localhost port=5432 user=postgres password=12345 dbname=name')

cur=conn.cursor()

#这里可以把sql拿出来去pg里面试一下看看是不是可用，是不是需要很久的时间，为了效率我这里限制只导出第一幅

cur.execute("select st_astiff(st_union(rast)) as rasttiff from public.temprature where filename='wc2.0_2.5m_tmax_01.tif' limit 1")

rasttiff=cur.fetchone()

#http://initd.org/psycopg/docs/cursor.html可以来psycopg官网看看相关的参数
if rasttiff is  not None:

#输出位置的选择，output是我的输出文件名
    open(r'E:\data\temprature\output.tif','wb').write(rasttiff[0])

cur.close()
conn.close()

6.找到我的e盘的对应位置，看我的图出来了哈哈。

 

具体参考原文如下：

http://postgis.net/docs/RT_AddOverviewConstraints.html

Synopsis

boolean AddOverviewConstraints(name ovschema, name ovtable, name ovcolumn, name refschema, name reftable, name refcolumn, int ovfactor);

boolean AddOverviewConstraints(name ovtable, name ovcolumn, name reftable, name refcolumn, int ovfactor);

Description

Adds constraints on a raster column that are used to display information in the raster_overviews raster catalog.

The ovfactor parameter represents the scale multiplier in the overview column: higher overview factors have lower resolution.

When the ovschema and refschema parameters are omitted, the first table found scanning the search_path will be used.

Availability: 2.0.0

Examples

CREATE TABLE res1 AS SELECT
ST_AddBand(
  ST_MakeEmptyRaster(1000, 1000, 0, 0, 2),
  1, '8BSI'::text, -129, NULL
) r1;

CREATE TABLE res2 AS SELECT
ST_AddBand(
  ST_MakeEmptyRaster(500, 500, 0, 0, 4),
  1, '8BSI'::text, -129, NULL
) r2;

SELECT AddOverviewConstraints('res2', 'r2', 'res1', 'r1', 2);

-- verify if registered correctly in the raster_overviews view --
SELECT o_table_name ot, o_raster_column oc,
       r_table_name rt, r_raster_column rc,
       overview_factor f
FROM raster_overviews WHERE o_table_name = 'res2';
  ot  | oc |  rt  | rc | f
------+----+------+----+---
 res2 | r2 | res1 | r1 | 2
(1 row)
		

See Also

Section 5.2.2, “Raster Overviews”, DropOverviewConstraints, ST_CreateOverview, AddRasterConstraints