analytics/informix/extend/spatial.8.22.FC2/examples/odbc/polyfromshape.c

/**************************************************************************
 *
 * polyfromshape.c
 *
 * This program demonstrates how to insert spatial data into a database
 * in ESRI shapefile format using the Informix Spatial DataBlade's
 * SE_PolyFromShape function.
 *
 * This program creates a table and inserts several polygons into it.
 * For an example program which shows how to insert all spatial types,
 * please see load_shapes.c
 *
 * The following command will likely compile this source:
 *
 *   cc -o polyfromshape polyfromshape.c \
 *        -I$INFORMIXDIR/incl/cli -L$INFORMIXDIR/lib/cli -lifcli -lifdmr
 *
 **************************************************************************/

#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <float.h>

#ifdef WINNT
#include <io.h>
#include <windows.h>
#include <conio.h>
#endif

#include <infxcli.h>     /* ODBC typedefs and data structures */
#include "odbcutils.c"   /* ODBC utility functions */

/*
 * Constants, macros, typedefs, data structures, & functions for 
 * processing ESRI shapefile format data.
 */
#include "commfuncs.h"
#include "commfuncs.c"
#include "shapefuncs.c"


void main (int argc, char **argv)
{
    SQLHDBC       hdbc;
    SQLHENV       henv;
    SQLHSTMT      hstmt;
    char          sql_stmt[200];
    int           rc;
    SDWORD        pcbvalue1, pcbvalue2, pcbvalue3;
    SDWORD        pcbvalue4, pcbvalue5;
    Geometry      geom = {0};
    Point         pt[30];
    int           offsets[5];
    int           max_alloced = 0;
    int           id;
    int           size;
    char         *name, *type;
    int           name_len, type_len;
    char         *zone_shape_buf;
    int           zone_shape_len;
    int           srid;
    
    /* Check for the correct number of arguments entered. */
    if (argc < 2)
    {
        printf ("Usage: %s <datasource> [<srid>]\n", argv[0]);
        exit (1);
    }

    /* Connect to the database. */
    server_connect ((UCHAR *) argv[1], &henv, &hdbc);

    /* Obtain srid to use for inserts */
    srid = (argc > 2) ? atoi(argv[2]) : 0;

    /* Allocate memory for the SQL statement handle and
     * associate the statement handle with the connection handle. */
    rc = SQLAllocHandle (SQL_HANDLE_STMT, hdbc, &hstmt);
	returncode_check (hdbc, (SQLHSTMT)NULL, rc, "SQLAllocHandle");

    sprintf(sql_stmt,
            "DROP TABLE sensitive_areas");
    rc = SQLExecDirect (hstmt, (UCHAR *) sql_stmt, SQL_NTS);

    sprintf(sql_stmt,
            "CREATE TABLE sensitive_areas (id integer, name varchar(128), "
            "size float, type varchar(32), zone ST_Polygon)");
    rc = SQLExecDirect (hstmt, (UCHAR *) sql_stmt, SQL_NTS);
    returncode_check (NULL, hstmt, rc, "SQLExecDirect");

    /* Initialize entries for non-spatial columns in the table */
    id = 289;
    name = "Johnson County Hospital";
    name_len = strlen(name);
    size = 102281.91;
    type = "hospital";
    type_len = strlen(type);

    /* Populate the polygon point array */
    offsets[0] = 0;
    pt[0].x = 22000; pt[0].y = 45000;
    pt[1].x = 22000; pt[1].y = 58000;
    pt[2].x = 28000; pt[2].y = 58000;
    pt[3].x = 28000; pt[3].y = 42000;
    pt[4].x = 25000; pt[4].y = 42000;
    pt[5].x = 25000; pt[5].y = 45000;
    pt[6].x = 22000; pt[6].y = 45000;

    geom.type = geomPolygon;
    geom.num_points = 7;
    geom.num_parts  = 1;
    geom.offsets = offsets;
    geom.pt = pt;

    /* Convert the polygon to an ESRI shape. */
    geom_to_shape (&geom, &max_alloced, &zone_shape_len, &zone_shape_buf);
               

    /* INFORMIX_EXTEST_BEGIN polyfromshape */

    /* Create the SQL insert statement to populate the sensitive_areas
     * table. The question marks are parameter markers that indicate
     * the column values that will be inserted at run time. */
    sprintf(sql_stmt,
            "INSERT INTO sensitive_areas (id, name, size, type, zone) "
            "VALUES(?, ?, ?, ?, SE_PolyFromShape(?, %d))", srid);

    /* Prepare the SQL statement for execution. */
    rc = SQLPrepare (hstmt, (unsigned char *)sql_stmt, SQL_NTS);
    returncode_check (NULL, hstmt, rc, "SQLPrepare");

    /* Bind the id to the first parameter. */
    pcbvalue1 = 0;
    rc = SQLBindParameter (hstmt, 1, SQL_PARAM_INPUT, SQL_C_SLONG,
                           SQL_INTEGER, 0, 0,
                           &id, 0, &pcbvalue1);
    returncode_check (NULL, hstmt, rc, "SQLBindParameter");

    /* Bind the name to the second parameter. */
    pcbvalue2 = name_len;
    rc = SQLBindParameter (hstmt, 2, SQL_PARAM_INPUT, SQL_C_CHAR,
                           SQL_CHAR, name_len, 0,
                           name, 0, &pcbvalue2);
    returncode_check (NULL, hstmt, rc, "SQLBindParameter");

    /* Bind the size to the third parameter. */
    pcbvalue3 = 0;
    rc = SQLBindParameter (hstmt, 3, SQL_PARAM_INPUT, SQL_C_FLOAT,
                           SQL_REAL, 0, 0,
                           &size, 0, &pcbvalue3);
    returncode_check (NULL, hstmt, rc, "SQLBindParameter");

    /* Bind the type to the fourth parameter. */
    pcbvalue4 = type_len;
    rc = SQLBindParameter (hstmt, 4, SQL_PARAM_INPUT, SQL_C_CHAR,
                           SQL_VARCHAR, type_len, 0,
                           type, type_len, &pcbvalue4);
    returncode_check (NULL, hstmt, rc, "SQLBindParameter");

    /* Bind the zone geometry to the fifth parameter. */
    pcbvalue5 = zone_shape_len;
    rc = SQLBindParameter (hstmt, 5, SQL_PARAM_INPUT, SQL_C_BINARY,
                           SQL_INFX_UDT_LVARCHAR, zone_shape_len, 0,
                           zone_shape_buf, zone_shape_len, &pcbvalue5);
    returncode_check (NULL, hstmt, rc, "SQLBindParameter");

    /* Execute the insert statement. */
    rc = SQLExecute (hstmt);
    returncode_check (NULL, hstmt, rc, "SQLExecute");

    /* INFORMIX_EXTEST_END polyfromshape */

    SQLFreeStmt (hstmt, SQL_CLOSE);          /* Close the statement handle */
    SQLFreeHandle (SQL_HANDLE_STMT, hstmt);  /* Free the statement handle */
    SQLDisconnect (hdbc);                    /* Close the connection */
    SQLFreeHandle (SQL_HANDLE_DBC, hdbc);    /* Free the database handle */
    SQLFreeHandle (SQL_HANDLE_ENV, henv);    /* Free the ODBC environment */
    printf( "\nTest Complete\n");
}