123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183 |
- /**************************************************************************
- *
- * mpolyfromshape.c
- *
- * This program demonstrates how to insert spatial data into a database
- * in ESRI shapefile format using the Informix Spatial DataBlade's
- * SE_MpolyFromShape 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 mpolyfromshape mpolyfromshape.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;
- Geometry geom = {0};
- Point pt[30];
- int offsets[10];
- int max_alloced = 0;
- int lot_id;
- char *lot_shape_buf;
- int lot_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");
- rc = SQLExecDirect (hstmt, (UCHAR *) "DROP TABLE lots", SQL_NTS);
- rc = SQLExecDirect (hstmt, (UCHAR *) "CREATE TABLE lots (lot_id integer, lot ST_Geometry)", SQL_NTS);
- returncode_check (NULL, hstmt, rc, "SQLExecDirect");
- /* Initialize entries for non-spatial columns in the table */
- lot_id = 203;
- /* The first polygon's exterior ring. */
- offsets[0] = 0;
- pt[0].x = 10; pt[0].y = 10;
- pt[1].x = 10; pt[1].y = 20;
- pt[2].x = 20; pt[2].y = 20;
- pt[3].x = 20; pt[3].y = 10;
- pt[4].x = 10; pt[4].y = 10;
- /* The first polygon's interior ring. */
- offsets[1] = 5;
- pt[5].x = 14; pt[5].y = 14;
- pt[6].x = 14; pt[6].y = 16;
- pt[7].x = 16; pt[7].y = 16;
- pt[8].x = 16; pt[8].y = 14;
- pt[9].x = 14; pt[9].y = 14;
- /* The second polygon's exterior ring. */
- offsets[2] = 10;
- pt[10].x = 100; pt[10].y = 100;
- pt[11].x = 100; pt[11].y = 200;
- pt[12].x = 200; pt[12].y = 200;
- pt[13].x = 200; pt[13].y = 100;
- pt[14].x = 100; pt[14].y = 100;
- /* The third polygons exterior ring. */
- offsets[3] = 15;
- pt[15].x = 50; pt[15].y = 50;
- pt[16].x = 50; pt[16].y = 60;
- pt[17].x = 60; pt[17].y = 60;
- pt[18].x = 60; pt[18].y = 50;
- pt[19].x = 50; pt[19].y = 50;
-
- /* The third polygon's interior ring. */
- offsets[4] = 20;
- pt[20].x = 52; pt[20].y = 52;
- pt[21].x = 52; pt[21].y = 54;
- pt[22].x = 54; pt[22].y = 54;
- pt[23].x = 54; pt[23].y = 52;
- pt[24].x = 52; pt[24].y = 52;
- /* The third polygon's interior ring. */
- offsets[5] = 25;
- pt[25].x = 56; pt[25].y = 56;
- pt[26].x = 56; pt[26].y = 58;
- pt[27].x = 58; pt[27].y = 58;
- pt[28].x = 58; pt[28].y = 56;
- pt[29].x = 56; pt[29].y = 56;
-
- geom.type = geomMultiPolygon;
- geom.num_points = 30;
- geom.num_parts = 6;
- geom.offsets = offsets;
- geom.pt = pt;
- /* Convert the polygon to an ESRI shape. */
- geom_to_shape (&geom, &max_alloced, &lot_shape_len, &lot_shape_buf);
-
- /* INFORMIX_EXTEST_BEGIN mpolyfromshape */
- /* Create the SQL insert statement to populate the lots table.
- * The question marks are parameter markers that indicate the
- * column values that will be inserted at run time. */
- sprintf(sql_stmt,
- "INSERT INTO lots (lot_id,lot)"
- "VALUES(?, SE_MpolyFromShape(?, %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 lot_id to the first parameter. */
- pcbvalue1 = 0;
- rc = SQLBindParameter (hstmt, 1, SQL_PARAM_INPUT, SQL_C_SLONG,
- SQL_INTEGER, 0, 0,
- &lot_id, 0, &pcbvalue1);
- returncode_check (NULL, hstmt, rc, "SQLBindParameter");
- /* Bind the lot geometry to the second parameter. */
- pcbvalue2 = lot_shape_len;
- rc = SQLBindParameter (hstmt, 2, SQL_PARAM_INPUT, SQL_C_BINARY,
- SQL_INFX_UDT_LVARCHAR, lot_shape_len, 0,
- lot_shape_buf, lot_shape_len, &pcbvalue2);
- returncode_check (NULL, hstmt, rc, "SQLBindParameter");
- /* Execute the insert statement. */
- rc = SQLExecute (hstmt);
- returncode_check (NULL, hstmt, rc, "SQLExecute");
- /* INFORMIX_EXTEST_END mpolyfromshape */
- 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");
- }
|