oracle - Oracle Spatial 问题：可以使用来自另一个 Oracle 表的查询结果填充 SDO_ORDINATE_ARRAY 吗？

Question

我的业务发展目标是使用 Oracle Spatial 来存储我们的坐标数据。目前，我们从在海洋区域进行测量的科学家那里接收坐标，这些坐标存储在我们的 Oracle 数据库中，作为点的数字对或多边形的长 varchar 数组。但是，我们希望通过使用 Oracle Spatial 来改进对这些数据的管理。

我们从科学家那里获得的坐标信息通常以 CSV 文件的形式出现，附带数据并作为字段条目加载到 Oracle 表中。

我知道我可以手动将顶点输入到 SDO_ORDINATE_ARRAY 中，但是我们经常会在一个 CSV 文件中提供数百个坐标对，这使得手动路径非常低效。

有人可以告诉我是否有办法通过从数据库中已存储信息的其他表中提取信息来填充 SDO_ORDINATE_ARRAY 的内容？

我尝试过的一个例子如下：

名为 GEOMTEST 的测试表由 NAME varchar2(50) COORDS varchar2(4000) COORD_GEOM SDO_GEOMETRY 组成

我用英国威尔士感兴趣的区域“卡迪根湾”填充了名称。COORDS 是我的多边形，作为数组存储在 varchar2 中。这是从 CSV 文件导入的。COORD_GEOM 是我希望将 COORDS 的内容传输到其中的内容。

我试图运行这段代码，但收到一个错误：

插入 geomtest (coord_geom) 值（SDO_GEOMETRY(2003,4326,null,SDO_ELEM_INFO_ARRAY(1,1003,1),SDO_ORDINATE_ARRAY 值（从 geomtest 中选择坐标））；

我使用 Toad 作为我的客户端，错误是“错误：第 18 行。第 120 列，第 18 行结束，第 125 列：找到'值'：保留字不能用作标识符。

我认为这与我在 INSERT 语句的 SDO_ORDINATE_ARRAY 部分中使用 select 语句有关，但不确定如何继续。

我将不胜感激任何建议，

非常感谢

肖恩

标签： oracleoracle19coracle-spatial

Answer 1

不幸的是，您不能直接将包含数字的字符串传递给SDO_ORDINATE_ARRAY构造函数。一种解决方案是编写一个自定义字符串标记器函数，它将坐标字符串解析为单独的数字，并构建一个SDO_ORDINATE_ARRAY对象。这是一个：

create or replace function tokenize (str clob)
return sdo_ordinate_array
is
  s clob := str||',';
  i number;
  j number;
  t sdo_ordinate_array := sdo_ordinate_array();
begin
  i := 1;
  loop
    j := instr(s, ',', i);
    exit when j = 0;
    t.extend();
    t(t.count) := to_number(substr(s,i,j-i));
    i := j+1;
  end loop;
  return t;
end;
/
show errors

这就是它的工作原理。首先让我们用几个例子创建一个简单的表：

drop table geomtest purge;
create table geomtest ( 
  id          number,
  name        varchar2(50 char), 
  coords      clob, 
  coord_geom  sdo_geometry 
);

insert into geomtest (id, name, coords)
values (
  2686,
  'TX/Mitchell', 
  '-101.17416, 32.527592, -101.17417, 32.523998, -101.1836, 32.087082, -100.82121, 32.086479, -100.66497, 32.085278, -100.66024, 32.5252, -101.17416, 32.527592'
);

insert into geomtest (id, name, coords)
values (
  2769,
  'TX/Yoakum',
  '-103.05616, 33.388332, -103.06416, 32.958992, -102.59455, 32.958733, -102.59436, 33.388393, -103.05616, 33.388332'
);
commit;

然后让我们使用 tokenizer 函数来更新几何列：

update geomtest 
set coord_geom = sdo_geometry(2003,4326,null,sdo_elem_info_array(1,1003,1),tokenize(coords));
commit;

检查结果：

SQL> select * from geomtest;

  ID NAME        COORDS                                                                                                                                                        COORD_GEOM(SDO_GTYPE, SDO_SRID, SDO_POINT(X, Y, Z), SDO_ELEM_INFO, SDO_ORDINATES)
---- ----------- ------------------------------------------------------------------------------------------------------------------------------------------------------------- ---------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
2686 TX/Mitchell -101.17416, 32.527592, -101.17417, 32.523998, -101.1836, 32.087082, -100.82121, 32.086479, -100.66497, 32.085278, -100.66024, 32.5252, -101.17416, 32.527592  SDO_GEOMETRY(2003, 4326, NULL, SDO_ELEM_INFO_ARRAY(1, 1003, 1), SDO_ORDINATE_ARRAY(-101.17416, 32.527592, -101.17417, 32.523998, -101.1836, 32.087082, -100.82121, 32.086479, -100.66497, 32.085278, -100.66024, 32.5252, -101.17416, 32.527592))
2769 TX/Yoakum   -103.05616, 33.388332, -103.06416, 32.958992, -102.59455, 32.958733, -102.59436, 33.388393, -103.05616, 33.388332                                             SDO_GEOMETRY(2003, 4326, NULL, SDO_ELEM_INFO_ARRAY(1, 1003, 1), SDO_ORDINATE_ARRAY(-103.05616, 33.388332, -103.06416, 32.958992, -102.59455, 32.958733, -102.59436, 33.388393, -103.05616, 33.388332))

2 rows selected.

笔记

1. 我使用一CLOB列来存储坐标。一个 4000 字节的字符串太小，无法容纳任何严肃的几何图形（除非您的所有形状都非常简单——只有几个点）。

2. 如果进行这种字符串到几何的转换，有更有效的方法，但它们意味着您使用面向几何的字符串表示法：GeoJSON、WKT、GML。这些自然得到 Oracle 的支持。它们还允许更复杂的结构，例如多多边形或带孔的多边形。

编辑：我重写了函数以直接返回一个SDO_GEOMETRY对象。这使它更易于使用：

create or replace function string_to_geom (str clob)
return sdo_geometry
is
  s clob := str||',';
  i number;
  j number;
  t sdo_ordinate_array := sdo_ordinate_array();
begin
  i := 1;
  loop
    j := instr(s, ',', i);
    exit when j = 0;
    t.extend();
    t(t.count) := to_number(substr(s,i,j-i));
    i := j+1;
  end loop;
  return sdo_geometry (2003, 4326, null, sdo_elem_info_array (1,1003,1), t);
end;
/
show errors

像这样使用它：

update geomtest 
set coord_geom = string_to_geom(coords);
commit;