空间分析入门

创建一个表格，几何图形列为 VARCHAR 列。 例如，下面的命令创建了一个适合加载点数据的表格：

CREATE TABLE PointData (PointID integer, the_geom varchar(200));

要在表格中插入点，请执行与以下命令类似的命令：

INSERT INTO PointData VALUES (1, inza..ST_Point(3423, 4356));
INSERT INTO PointData VALUES (1, inza..ST_WKTToSQL('Point (3423 4356)'));

以下命令将显示_geom 列中数据的几何类型：

SELECT inza..ST_GeometryType(the_geom) FROM PointData;

ST_GEOMETRYTYPE
---------------------
ST_Point
ST_Point

要创建包含多边形的表格，请执行以下命令：

CREATE TABLE Polys (PolyID integer, the_geom varchar(200));

要在表格中添加多边形，可发出与以下命令类似的命令，指定点作为两个正多边形的顶点：

INSERT INTO Polys VALUES (1, inza..ST_WKTToSQL('Polygon ((1000 1000, 1000 5000, 4000 5000, 4000 1000, 1000 1000))'));
INSERT INTO Polys VALUES (2, inza..ST_WKTToSQL('Polygon ((100 100, 100 500, 400 500, 400 100, 100 100))'));

利用 Polys 表中定义的多边形，可以通过 ST_Intersects()函数进行连接，从而发出基本的多边形内点查询，如下所示：

SELECT PointID, PolyID FROM PointData AS a, Polys AS b WHERE inza..ST_Intersects(a.the_geom, b.the_geom);

 POINTID | POLYID
---------+--------
       1 |      1
       1 |      1

例如，下面的命令将创建一个名为 "geoms的小表，并在表中添加多边形和线字符串对象：

CREATE TABLE geoms (PolyID INTEGER, the_geom VARCHAR(64000));
INSERT INTO geoms VALUES (1, inza..ST_WKTtoSQL('Polygon ((1 1, 1 4, 3.5 2.5, 6 4, 6 1, 1 1))'));
INSERT INTO geoms VALUES (2, inza..ST_WKTtoSQL('Linestring (1 1, 1 4, 6 4, 6 1)'));

要以可读形式显示 "geoms中的数据，请使用空间转换函数之一：

• 要将空间数据呈现为文本，请使用ST_AsText:

  SELECT inza..ST_AsText(the_geom) FROM geoms;
  
  ST_ASTEXT
  ---------------------------------------------------------------------
  POLYGON ((1 1, 1 4, 3.5 2.5, 6 4, 6 1, 1 1))
  LINESTRING (1 1, 1 4, 6 4, 6 1)

• 要以 Keyhole 标记语言（KML）呈现空间数据，请使用ST_AsKML:

  SELECT inza..ST_AsKML(the_geom) FROM geoms;
  ST_ASKML 
  ---------------------------------------------------------------------
  <LineString><coordinates>1,1 1,4 6,4 6,1</coordinates></LineString>
  <Polygon><outerBoundaryIs><LinearRing><coordinates>1,1 1,4 3.5,2.5 6,4 6,1 1,1</coordinates></LinearRing></outerBoundaryIs></Polygon>

存储空间数据的列的数据类型取决于正在使用的盒的类型：

空间

空间数据保存在 VARCHAR 列中。

空间 ESRI

空间数据保存在 ST_GEOMETRY 列中。

功能处理引擎（FME）工作台应用程序是 Safe Software Inc、 帮助准备数据集并将其加载到Netezza Performance Server数据库中。 例如，FME 工作台可以在以下方面为您提供帮助：

• 它可以检测到顶点超过 3990 个的多边形，并采取用户指定的操作，如将几何体切成更小的几何体、跳过并忽略几何体、中止加载过程，或对几何体进行泛化（或平滑），使其足够小以适应目标区域。
• 它可以检测非简单几何图形，ESRI 空间插件支持这些几何图形，但空间插件不支持。 它可以将其过滤掉，或者进行缓冲，以便随后将其转化为简单的几何图形。
• 它可以加载二进制空间数据，同时检测和转义加载功能不支持的已知二进制值。

有关 FME Workbench 应用程序的更多信息，请参阅 Safe Software Inc.

要在用户数据库中注册空间函数和存储过程，请执行以下步骤：

1. 输入以下命令创建新用户数据库，其中 "newdb代表要创建的用户数据库名称。

   % nzsql -c "create database newdb;"
    CREATE DATABASE

2. 输入以下命令注册 nzspatial：

   % nzcm -d newdb -r nzspatial

3. 输入以下命令运行create_inza_db.sh脚本：

   % /nz/export/ae/utilities/bin/create_inza_db.sh newdb

在用户数据库中，空间函数和存储过程的注册方式如下：

• 在 INZA 模式中，当启用全模式支持Netezza Performance Server功能时
• 在默认模式中，当全模式支持Netezza Performance Server功能禁用时

convert.sh脚本将包含 WKB 几何数据 VARCHAR 列的表转换为包含 ESRI 二进制几何数据 ST_GEOMETRY 列的新表。 列名保持不变。 该脚本还对已知二进制 (WKB) 数据库中的两个非标准 SRID 进行了规范化处理：

• SRID 1111 将转换为您指定的 SRID（默认为 4269）。
• SRID 1234 将转换为您指定的 SRID（默认为 27700）。

新表可以在同一个数据库中，也可以在不同的数据库中。

您必须以用户 nz 的身份运行 "convert.sh脚本。 在运行脚本之前，必须满足以下条件：

• nzspatial 函数。
• nzspatial_esri 函数。
• 目标数据库已 INA 初始化。

语法