Oracle空间问题：可以使用另一个Oracle表的查询结果填充SDO_ORDINATE_ARRAY吗？

Question

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

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

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

有人可以告诉我，如果有一种方法可以通过从已经存储信息的数据库中的其他表中提取信息来填充SDO_ORDINATE_ARRAY的内容吗？

下面是我尝试过的一个例子:

名为GEOMTEST的测试表由NAME varchar 2(50)COORDS varchar 2(4000)COORD _ GEOM SDO _ GEOMETRY组成

我用英国威尔士的“羊毛衫湾”这个有趣的区域来填充这个名字。COORDS是我在varchar2中存储为数组的多边形。它是从CSV文件导入的。COORD_GEOM是我希望将COORDS的内容转换成的内容。

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

将值（SDO_GEOMETRY（20034326，null，SDO_ ELEM_INFO_ARRAY（11003,1））和SDOUOr座标_；

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

我假设这与我在INSERT语句的SDO_。

我很感激任何建议，

多谢

肖恩

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;

然后，让我们使用分词器函数来更新几何列：

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.

笔记

>

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

如果进行字符串到几何的转换，有更有效的方法，但它们意味着您使用面向几何的字符串符号：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;