Come iniziare con l'analisi spaziale

Creare una tabella con una colonna geometria come colonna VARCHAR. Ad esempio, il comando seguente crea una tabella adatta al caricamento di dati puntuali:

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

Per inserire un punto nella tabella, si devono impartire comandi simili ai seguenti:

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

Il comando seguente visualizza il tipo di geometria dei dati nella colonna_geom:

SELECT inza..ST_GeometryType(the_geom) FROM PointData;

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

Per creare una tabella contenente poligoni, eseguire il seguente comando:

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

Per aggiungere poligoni alla tabella, si impartiscono comandi simili ai seguenti, che specificano i punti che sono i vertici di due poligoni quadrati:

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))'));

Con i poligoni definiti nella tabella Polys, è possibile eseguire interrogazioni di base punto-poligono unendole tramite la funzione ST_Intersects(), come segue:

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

Ad esempio, i comandi seguenti creano una piccola tabella denominata 'geoms e aggiungono un oggetto poligono e un oggetto linea alla tabella:

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)'));

Per visualizzare i dati in 'geoms in forma leggibile, utilizzare una delle funzioni di conversione spaziale:

• Per rendere i dati territoriali come testo, utilizzare 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)

• Per visualizzare i dati spaziali in Keyhole Markup Language (KML), utilizzare 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>

Il tipo di dati delle colonne che memorizzano i dati spaziali dipende dal tipo di cartuccia utilizzata:

Spaziale

I dati spaziali vengono salvati in una colonna VARCHAR.

ESRI spaziale

I dati spaziali vengono salvati in una colonna ST_GEOMETRY.

L'applicazione Feature Manipulation Engine (FME) Workbench, un prodotto di Safe Software Inc, aiuta a preparare i set di dati e a caricarli nel database Netezza Performance Server. Ad esempio, FME Workbench può aiutarvi nei seguenti modi:

• Può rilevare poligoni con più di 3990 vertici e intraprendere un'azione specificata dall'utente, ad esempio tagliare la geometria in geometrie più piccole, saltare e ignorare la geometria, interrompere il processo di caricamento o generalizzare (o smussare) la geometria per renderla abbastanza piccola da rientrare nel campo di destinazione.
• Può rilevare geometrie non semplici, che sono supportate dalla cartuccia ESRI Spatial, ma non dalla cartuccia Spatial. Può filtrarli o tamponarli in modo da poterli trasformare in seguito in geometrie semplici.
• È in grado di caricare dati spaziali binari, rilevando ed eliminando i valori binari noti che non sono supportati dalle funzioni di caricamento.

Per ulteriori informazioni sull'applicazione FME Workbench, consultare la documentazione fornita da Safe Software Inc.

Per registrare le funzioni spaziali e le stored procedure da utilizzare in un database utente, procedere come segue:

1. Creare il nuovo database utente immettendo il seguente comando, dove 'newdb rappresenta il nome del database utente che si desidera creare.

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

2. Registrate nzspatial immettendo il seguente comando:

   % nzcm -d newdb -r nzspatial

3. Eseguite lo script create_inza_db.sh immettendo il seguente comando:

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

Nel database utente, le funzioni spaziali e le stored procedure sono registrate come segue:

• Nello schema INZA, quando è abilitata la funzione Netezza Performance Server di supporto dello schema completo
• Nello schema predefinito, quando lo schema completo supporta la funzionalità di Netezza Performance Server, la funzione è disabilitata

Lo script " convert.sh converte una tabella con una colonna VARCHAR contenente dati geometrici WKB in una nuova tabella con una colonna ST_GEOMETRY contenente dati geometrici binari ESRI. I nomi delle colonne rimangono invariati. Lo script normalizza anche i due SRID non standard nel database well known binary (WKB):

• Il SRID 1111 viene convertito nel SRID specificato (4269 è il valore predefinito).
• Il SRID 1234 viene convertito nel SRID specificato (27700 è il valore predefinito).

La nuova tabella può trovarsi nello stesso database o in un database diverso.

È necessario eseguire lo script 'convert.sh come utente nz. Prima di poter eseguire lo script, devono essere soddisfatte le seguenti condizioni:

• esistono funzioni nzspatial.
• esistono funzioni nzspatial_esri.
• Il database di destinazione è inizializzato INA.

Sintassi