static void do_geom_test(char * in, char * out)
{
LWGEOM *g, *h;
char *tmp;
g = lwgeom_from_wkt(in, LW_PARSER_CHECK_NONE);
h = lwgeom_homogenize(g);
tmp = lwgeom_to_ewkt(h);
if (strcmp(tmp, out))
fprintf(stderr, "\nIn: %s\nOut: %s\nExp: %s\n",
in, tmp, out);
CU_ASSERT_STRING_EQUAL(tmp, out);
lwfree(tmp);
lwgeom_free(g);
/* See http://trac.osgeo.org/postgis/ticket/1104 */
lwgeom_free(h);
}
Datum geom_from_kml(PG_FUNCTION_ARGS)
{
GSERIALIZED *geom;
LWGEOM *lwgeom, *hlwgeom;
xmlDocPtr xmldoc;
text *xml_input;
int xml_size;
char *xml;
bool hasz=true;
xmlNodePtr xmlroot=NULL;
/* Get the KML stream */
if (PG_ARGISNULL(0)) PG_RETURN_NULL();
xml_input = PG_GETARG_TEXT_P(0);
xml = text2cstring(xml_input);
xml_size = VARSIZE(xml_input) - VARHDRSZ;
/* Begin to Parse XML doc */
xmlInitParser();
xmldoc = xmlParseMemory(xml, xml_size);
if (!xmldoc || (xmlroot = xmlDocGetRootElement(xmldoc)) == NULL)
{
xmlFreeDoc(xmldoc);
xmlCleanupParser();
lwerror("invalid KML representation");
}
lwgeom = parse_kml(xmlroot, &hasz);
/* Homogenize geometry result if needed */
if (lwgeom->type == COLLECTIONTYPE)
{
hlwgeom = lwgeom_homogenize(lwgeom);
lwgeom_release(lwgeom);
lwgeom = hlwgeom;
}
lwgeom_add_bbox(lwgeom);
/* KML geometries could be either 2 or 3D
*
* So we deal with 3D in all structures allocation, and flag hasz
* to false if we met once a missing Z dimension
* In this case, we force recursive 2D.
*/
if (!hasz)
{
LWGEOM *tmp = lwgeom_force_2d(lwgeom);
lwgeom_free(lwgeom);
lwgeom = tmp;
}
geom = geometry_serialize(lwgeom);
lwgeom_free(lwgeom);
xmlFreeDoc(xmldoc);
xmlCleanupParser();
PG_RETURN_POINTER(geom);
}