capi_read.c

/*
* # GEOS C API example 1
*
* Reads two WKT representations and calculates the
* intersection, prints it out, and cleans up.
*/

/* To print to stdout */
#include <stdio.h>
#include <stdarg.h>

/* Only the CAPI header is required */
#include <geos_c.h>

/*
* GEOS requires two message handlers to return
* error and notice message to the calling program.
*
*   typedef void(* GEOSMessageHandler) (const char *fmt,...)
*
* Here we stub out an example that just prints the
* messages to stdout.
*/
static void
geos_message_handler(const char* fmt, ...)
{
    va_list ap;
    va_start(ap, fmt);
    vprintf (fmt, ap);
    va_end(ap);
}


int main()
{
    /* Send notice and error messages to our stdout handler */
    initGEOS(geos_message_handler, geos_message_handler);

    /* Two squares that overlap */
    const char* wkt_a = "POLYGON((0 0, 10 0, 10 10, 0 10, 0 0))";
    const char* wkt_b = "POLYGON((5 5, 15 5, 15 15, 5 15, 5 5))";

    /* Read the WKT into geometry objects */
    GEOSWKTReader* reader = GEOSWKTReader_create();
    GEOSGeometry* geom_a = GEOSWKTReader_read(reader, wkt_a);
    GEOSGeometry* geom_b = GEOSWKTReader_read(reader, wkt_b);

    /* Calculate the intersection */
    GEOSGeometry* inter = GEOSIntersection(geom_a, geom_b);

    /* Convert result to WKT */
    GEOSWKTWriter* writer = GEOSWKTWriter_create();
    /* Trim trailing zeros off output */
    GEOSWKTWriter_setTrim(writer, 1);
    char* wkt_inter = GEOSWKTWriter_write(writer, inter);

    /* Print answer */
    printf("Geometry A:         %s\n", wkt_a);
    printf("Geometry B:         %s\n", wkt_b);
    printf("Intersection(A, B): %s\n", wkt_inter);

    /* Clean up everything we allocated */
    GEOSWKTReader_destroy(reader);
    GEOSWKTWriter_destroy(writer);
    GEOSGeom_destroy(geom_a);
    GEOSGeom_destroy(geom_b);
    GEOSGeom_destroy(inter);
    GEOSFree(wkt_inter);

    /* Clean up the global context */
    finishGEOS();

    /* Done */
    return 0;
}