要作为列插入表中的栅格记录的像素值

时间:2012-10-30 19:09:55

标签: postgresql pivot postgis crosstab raster

我有一个包含以下列的表:

(ID, row_num, col_num, pix_centroid, pix_val1). 

我有超过1000条记录。我使用以下方法插入数据:

insert into pixelbased (row_num, col_num, pix_centroid, pix_val)
select
    (ST_PixelAsPolygons(rast, 1)).x as X,
    (ST_PixelAsPolygons(rast, 1)).y as Y,
    (ST_Centroid((ST_PixelAsPolygons(rast, 1)).geom)) as geom,
    (ST_PixelAsPolygons(rast, 1)).val as pix_val1
from mytable 
where rid=1`

现在我尝试将所有其他记录作为列插入,_pix_val1_列对我来说很重要。所有其他列将保持不变。换句话说,我希望最终表格包含以下列:

(ID, row_num, col_num, pix_centroid, pix_val1, pix_val2, pix_val3, ....)

有办法吗?

1 个答案:

答案 0 :(得分:2)

如果可能的话,我想将这些数据存储为bytea中的位图。以下是如何获取一系列字节值并将其转换为bytea:

WITH bytes(b) AS (SELECT x % 256 FROM generate_series(1,53000) x)
SELECT ('\x'||string_agg(lpad(to_hex(b),2,'0'),''))::bytea FROM bytes;

您可以使用substr函数访问字节数组的字段或范围。此bytea被组织为线性像素阵列,但您可能会发现将其组织成更传统的位图格式更有用。此外,如果您的像素超过一个字节,您可能需要应对big-endian vs little-endian。您可以在SQL中执行此操作,但在PL / Perl等过程语言中可能会更容易。


如果不这样做,多维数组将是一个合理的选择。

使用generate_series语句替代pix_val字段以方便测试,此查询使用两个聚合传递生成一个二维整数数组:

SELECT ('{'||string_agg(subarray, ',')||'}')::integer[] AS arr
FROM (
   SELECT array_agg(x order by x)::text 
   FROM generate_series(1,53000) x
   GROUP BY width_bucket(x, 1, 53001, 100)
) a(subarray);

array_agg无法聚合数组这一事实使得二维数组的字符串文字形式的不幸使用成为必要。在我看来,这是PostgreSQL中的一个真正的疣;一般来说,它的多维数组很奇怪,与大多数应用程序和语言实现数组的方式不一致。

您可以通过索引来获取数组中的字段。例如:

regress=> SELECT ('{'||string_agg(subarray, ',')||'}')::integer[] AS arr INTO test FROM (SELECT array_agg(x order by x)::text from generate_series(1,53000) x GROUP BY width_bucket(x, 1, 53001, 100)) a(subarray);

regress=> \d test

      Table "public.test"
 Column |   Type    | Modifiers 
--------+-----------+-----------
 arr   | integer[] | 

test包含一个包含两个维度的数组:

regress=> \x
regress=> select array_dims(test.arr), array_ndims(test.arr), array_length(test.arr,1), array_length(test.arr,2) FROM test;
-[ RECORD 1 ]+---------------
array_dims   | [1:100][1:530]
array_ndims  | 2
array_length | 100
array_length | 530

我可以使用两级索引获取元素:

regress=> SELECT test.arr[4][4] FROM test;
 arr  
------
 1594
(1 row)

或带有切片的“列”:

regress=> SELECT test.arr[4:4][1:530] FROM test;

奇怪的是,这仍然是一个二维数组,顶部维度只是一个元素深。如果需要,您可以使用unnestarray_agg将其展平(效率低下)。

PostgreSQL中的二维数组有些奇怪,正如您所看到的那样,但您正在尝试做的事情也是如此。