mirror of
https://git.osgeo.org/gitea/postgis/postgis
synced 2024-10-24 00:52:40 +00:00
93087d39dc
git-svn-id: http://svn.osgeo.org/postgis/trunk@8053 b70326c6-7e19-0410-871a-916f4a2858ee
187 lines
5.6 KiB
C
187 lines
5.6 KiB
C
/**********************************************************************
|
|
* $Id$
|
|
*
|
|
* PostGIS - Spatial Types for PostgreSQL
|
|
* http://postgis.refractions.net
|
|
*
|
|
* Copyright 2006 Corporacion Autonoma Regional de Santander
|
|
* Eduin Carrillo <yecarrillo@cas.gov.co>
|
|
* Copyright 2010 Paul Ramsey <pramsey@cleverelephant.ca>
|
|
*
|
|
* This is free software; you can redistribute and/or modify it under
|
|
* the terms of hte GNU General Public Licence. See the COPYING file.
|
|
*
|
|
**********************************************************************/
|
|
|
|
#include "liblwgeom_internal.h"
|
|
#include "stringbuffer.h"
|
|
|
|
static int lwgeom_to_kml2_sb(const LWGEOM *geom, int precision, const char *prefix, stringbuffer_t *sb);
|
|
static int lwpoint_to_kml2_sb(const LWPOINT *point, int precision, const char *prefix, stringbuffer_t *sb);
|
|
static int lwline_to_kml2_sb(const LWLINE *line, int precision, const char *prefix, stringbuffer_t *sb);
|
|
static int lwpoly_to_kml2_sb(const LWPOLY *poly, int precision, const char *prefix, stringbuffer_t *sb);
|
|
static int lwcollection_to_kml2_sb(const LWCOLLECTION *col, int precision, const char *prefix, stringbuffer_t *sb);
|
|
static int ptarray_to_kml2_sb(const POINTARRAY *pa, int precision, stringbuffer_t *sb);
|
|
|
|
/*
|
|
* KML 2.2.0
|
|
*/
|
|
|
|
/* takes a GEOMETRY and returns a KML representation */
|
|
char*
|
|
lwgeom_to_kml2(const LWGEOM *geom, int precision, const char *prefix)
|
|
{
|
|
stringbuffer_t *sb;
|
|
int rv;
|
|
char *kml;
|
|
|
|
/* Can't do anything with empty */
|
|
if( lwgeom_is_empty(geom) )
|
|
return NULL;
|
|
|
|
sb = stringbuffer_create();
|
|
rv = lwgeom_to_kml2_sb(geom, precision, prefix, sb);
|
|
|
|
if ( rv == LW_FAILURE )
|
|
{
|
|
stringbuffer_destroy(sb);
|
|
return NULL;
|
|
}
|
|
|
|
kml = stringbuffer_getstringcopy(sb);
|
|
stringbuffer_destroy(sb);
|
|
|
|
return kml;
|
|
}
|
|
|
|
static int
|
|
lwgeom_to_kml2_sb(const LWGEOM *geom, int precision, const char *prefix, stringbuffer_t *sb)
|
|
{
|
|
switch (geom->type)
|
|
{
|
|
case POINTTYPE:
|
|
return lwpoint_to_kml2_sb((LWPOINT*)geom, precision, prefix, sb);
|
|
|
|
case LINETYPE:
|
|
return lwline_to_kml2_sb((LWLINE*)geom, precision, prefix, sb);
|
|
|
|
case POLYGONTYPE:
|
|
return lwpoly_to_kml2_sb((LWPOLY*)geom, precision, prefix, sb);
|
|
|
|
case MULTIPOINTTYPE:
|
|
case MULTILINETYPE:
|
|
case MULTIPOLYGONTYPE:
|
|
return lwcollection_to_kml2_sb((LWCOLLECTION*)geom, precision, prefix, sb);
|
|
|
|
default:
|
|
lwerror("lwgeom_to_kml2: '%s' geometry type not supported", lwtype_name(geom->type));
|
|
return LW_FAILURE;
|
|
}
|
|
}
|
|
|
|
static int
|
|
ptarray_to_kml2_sb(const POINTARRAY *pa, int precision, stringbuffer_t *sb)
|
|
{
|
|
int i, j;
|
|
int dims = FLAGS_GET_Z(pa->flags) ? 3 : 2;
|
|
POINT4D pt;
|
|
double *d;
|
|
|
|
for ( i = 0; i < pa->npoints; i++ )
|
|
{
|
|
getPoint4d_p(pa, i, &pt);
|
|
d = (double*)(&pt);
|
|
if ( i ) stringbuffer_append(sb," ");
|
|
for (j = 0; j < dims; j++)
|
|
{
|
|
if ( j ) stringbuffer_append(sb,",");
|
|
if( fabs(d[j]) < OUT_MAX_DOUBLE )
|
|
{
|
|
if ( stringbuffer_aprintf(sb, "%.*f", precision, d[j]) < 0 ) return LW_FAILURE;
|
|
}
|
|
else
|
|
{
|
|
if ( stringbuffer_aprintf(sb, "%g", d[j]) < 0 ) return LW_FAILURE;
|
|
}
|
|
stringbuffer_trim_trailing_zeroes(sb);
|
|
}
|
|
}
|
|
return LW_SUCCESS;
|
|
}
|
|
|
|
|
|
static int
|
|
lwpoint_to_kml2_sb(const LWPOINT *point, int precision, const char *prefix, stringbuffer_t *sb)
|
|
{
|
|
/* Open point */
|
|
if ( stringbuffer_aprintf(sb, "<%sPoint><%scoordinates>", prefix, prefix) < 0 ) return LW_FAILURE;
|
|
/* Coordinate array */
|
|
if ( ptarray_to_kml2_sb(point->point, precision, sb) == LW_FAILURE ) return LW_FAILURE;
|
|
/* Close point */
|
|
if ( stringbuffer_aprintf(sb, "</%scoordinates></%sPoint>", prefix, prefix) < 0 ) return LW_FAILURE;
|
|
return LW_SUCCESS;
|
|
}
|
|
|
|
static int
|
|
lwline_to_kml2_sb(const LWLINE *line, int precision, const char *prefix, stringbuffer_t *sb)
|
|
{
|
|
/* Open linestring */
|
|
if ( stringbuffer_aprintf(sb, "<%sLineString><%scoordinates>", prefix, prefix) < 0 ) return LW_FAILURE;
|
|
/* Coordinate array */
|
|
if ( ptarray_to_kml2_sb(line->points, precision, sb) == LW_FAILURE ) return LW_FAILURE;
|
|
/* Close linestring */
|
|
if ( stringbuffer_aprintf(sb, "</%scoordinates></%sLineString>", prefix, prefix) < 0 ) return LW_FAILURE;
|
|
|
|
return LW_SUCCESS;
|
|
}
|
|
|
|
static int
|
|
lwpoly_to_kml2_sb(const LWPOLY *poly, int precision, const char *prefix, stringbuffer_t *sb)
|
|
{
|
|
int i, rv;
|
|
|
|
/* Open polygon */
|
|
if ( stringbuffer_aprintf(sb, "<%sPolygon>", prefix) < 0 ) return LW_FAILURE;
|
|
for ( i = 0; i < poly->nrings; i++ )
|
|
{
|
|
/* Inner or outer ring opening tags */
|
|
if( i )
|
|
rv = stringbuffer_aprintf(sb, "<%sinnerBoundaryIs><%sLinearRing><%scoordinates>", prefix, prefix, prefix);
|
|
else
|
|
rv = stringbuffer_aprintf(sb, "<%souterBoundaryIs><%sLinearRing><%scoordinates>", prefix, prefix, prefix);
|
|
if ( rv < 0 ) return LW_FAILURE;
|
|
|
|
/* Coordinate array */
|
|
if ( ptarray_to_kml2_sb(poly->rings[i], precision, sb) == LW_FAILURE ) return LW_FAILURE;
|
|
|
|
/* Inner or outer ring closing tags */
|
|
if( i )
|
|
rv = stringbuffer_aprintf(sb, "</%scoordinates></%sLinearRing></%sinnerBoundaryIs>", prefix, prefix, prefix);
|
|
else
|
|
rv = stringbuffer_aprintf(sb, "</%scoordinates></%sLinearRing></%souterBoundaryIs>", prefix, prefix, prefix);
|
|
if ( rv < 0 ) return LW_FAILURE;
|
|
}
|
|
/* Close polygon */
|
|
if ( stringbuffer_aprintf(sb, "</%sPolygon>", prefix) < 0 ) return LW_FAILURE;
|
|
|
|
return LW_SUCCESS;
|
|
}
|
|
|
|
static int
|
|
lwcollection_to_kml2_sb(const LWCOLLECTION *col, int precision, const char *prefix, stringbuffer_t *sb)
|
|
{
|
|
int i, rv;
|
|
|
|
/* Open geometry */
|
|
if ( stringbuffer_aprintf(sb, "<%sMultiGeometry>", prefix) < 0 ) return LW_FAILURE;
|
|
for ( i = 0; i < col->ngeoms; i++ )
|
|
{
|
|
rv = lwgeom_to_kml2_sb(col->geoms[i], precision, prefix, sb);
|
|
if ( rv == LW_FAILURE ) return LW_FAILURE;
|
|
}
|
|
/* Close geometry */
|
|
if ( stringbuffer_aprintf(sb, "</%sMultiGeometry>", prefix) < 0 ) return LW_FAILURE;
|
|
|
|
return LW_SUCCESS;
|
|
}
|