27 #include "../postgis_config.h"
35 #include "access/gist.h"
36 #include "access/itup.h"
39 #include "utils/elog.h"
40 #include "mb/pg_wchar.h"
41 #include "lib/stringinfo.h"
42 #include "utils/array.h"
43 #include "utils/builtins.h"
44 #include "utils/lsyscache.h"
48 #include "lwgeom_pg.h"
50 #include "lwgeom_transform.h"
53 #include "access/htup_details.h"
86 char *input = PG_GETARG_CSTRING(0);
87 int32 geom_typmod = -1;
94 if ( (PG_NARGS()>2) && (!PG_ARGISNULL(2)) ) {
95 geom_typmod = PG_GETARG_INT32(2);
101 if (
str[0] ==
'\0' ) {
102 ereport(ERROR,(errmsg(
"parse error - invalid geometry")));
107 if( strncasecmp(
str,
"SRID=",5) == 0 )
111 while ( tmp && *tmp !=
';' )
115 if ( tmp && *(tmp+1) ==
'0' )
131 size_t hexsize = strlen(
str);
148 PG_PARSER_ERROR(lwg_parser_result);
151 lwgeom = lwg_parser_result.
geom;
158 if ( geom_typmod >= 0 )
161 POSTGIS_DEBUG(3,
"typmod and geometry were consistent");
165 POSTGIS_DEBUG(3,
"typmod was -1");
171 PG_RETURN_POINTER(ret);
201 GSERIALIZED *pg_lwgeom = PG_GETARG_GSERIALIZED_P(0);
202 text *format_text = PG_GETARG_TEXT_P(1);
205 char *format_str = NULL;
207 char * formatted_str;
208 text * formatted_text;
215 lwpgerror(
"Only points are supported, you tried type %s.",
lwtype_name(geom_type));
220 if (format_text == NULL) {
221 lwpgerror(
"ST_AsLatLonText: invalid format string (null");
226 lwpgerror(
"ST_AsLatLonText: invalid coordinate");
231 assert(format_str != NULL);
235 tmp = (
char *)pg_do_encoding_conversion(
236 (uint8_t *)format_str, strlen(format_str), GetDatabaseEncoding(), PG_UTF8);
238 if ( tmp != format_str ) {
245 assert(formatted_str != NULL);
249 tmp = (
char *)pg_do_encoding_conversion(
250 (uint8_t *)formatted_str, strlen(formatted_str),
251 PG_UTF8, GetDatabaseEncoding());
253 if ( tmp != formatted_str) {
254 pfree(formatted_str);
259 formatted_text = cstring_to_text(formatted_str);
260 pfree(formatted_str);
262 PG_RETURN_POINTER(formatted_text);
284 PG_RETURN_CSTRING(hexwkb);
303 if ( (PG_NARGS()>1) && (!PG_ARGISNULL(1)) )
305 type = PG_GETARG_TEXT_P(1);
307 if ( ! strncmp(VARDATA(
type),
"xdr", 3) ||
308 ! strncmp(VARDATA(
type),
"XDR", 3) )
324 text_size = hexwkb_size - 1 + VARHDRSZ;
325 result = palloc(text_size);
326 memcpy(VARDATA(result), hexwkb, hexwkb_size - 1);
327 SET_VARSIZE(result, text_size);
331 PG_FREE_IF_COPY(geom, 0);
332 PG_RETURN_TEXT_P(result);
358 result = cstring_to_text(hexwkb);
362 PG_FREE_IF_COPY(geom, 0);
363 PG_RETURN_TEXT_P(result);
376 bytea *bytea_wkb = PG_GETARG_BYTEA_P(0);
379 uint8_t *wkb = (uint8_t*)VARDATA(bytea_wkb);
383 lwpgerror(
"Unable to parse WKB");
385 if ((PG_NARGS() > 1) && (!PG_ARGISNULL(1)))
387 int32 srid = PG_GETARG_INT32(1);
396 PG_FREE_IF_COPY(bytea_wkb, 0);
397 PG_RETURN_POINTER(geom);
407 bytea *bytea_twkb = PG_GETARG_BYTEA_P(0);
410 uint8_t *twkb = (uint8_t*)VARDATA(bytea_twkb);
419 PG_FREE_IF_COPY(bytea_twkb, 0);
420 PG_RETURN_POINTER(geom);
438 if ( (PG_NARGS()>1) && (!PG_ARGISNULL(1)) )
440 type = PG_GETARG_TEXT_P(1);
442 if ( ! strncmp(VARDATA(
type),
"xdr", 3) ||
443 ! strncmp(VARDATA(
type),
"XDR", 3) )
452 wkb_size= VARSIZE_ANY_EXHDR(geom);
460 result = palloc(wkb_size + VARHDRSZ);
461 memcpy(VARDATA(result), wkb, wkb_size);
462 SET_VARSIZE(result, wkb_size+VARHDRSZ);
466 PG_FREE_IF_COPY(geom, 0);
467 PG_RETURN_BYTEA_P(result);
483 if ( PG_ARGISNULL(0) ) PG_RETURN_NULL();
485 geom = PG_GETARG_GSERIALIZED_P(0);
491 if ( PG_NARGS() > 1 && ! PG_ARGISNULL(1) )
492 sp.precision_xy = PG_GETARG_INT32(1);
495 if ( PG_NARGS() > 2 && ! PG_ARGISNULL(2) )
496 sp.precision_z = PG_GETARG_INT32(2);
499 if ( PG_NARGS() > 3 && ! PG_ARGISNULL(3) )
500 sp.precision_m = PG_GETARG_INT32(3);
506 if ( PG_NARGS() > 4 && ! PG_ARGISNULL(4) && PG_GETARG_BOOL(4) )
510 if ( PG_NARGS() > 5 && ! PG_ARGISNULL(5) && PG_GETARG_BOOL(5) )
518 result = palloc(twkb_size + VARHDRSZ);
519 memcpy(VARDATA(result), twkb, twkb_size);
520 SET_VARSIZE(result, twkb_size + VARHDRSZ);
522 PG_RETURN_BYTEA_P(result);
529 ArrayType *arr_geoms = NULL;
530 ArrayType *arr_ids = NULL;
531 int num_geoms, num_ids, i = 0;
533 ArrayIterator iter_geoms, iter_ids;
534 bool null_geom, null_id;
535 Datum val_geom, val_id;
537 int is_homogeneous =
true;
538 uint32_t subtype = 0;
542 int64_t *idlist = NULL;
551 if ( PG_NARGS() < 2 || PG_ARGISNULL(0) || PG_ARGISNULL(1) )
554 arr_geoms = PG_GETARG_ARRAYTYPE_P(0);
555 arr_ids = PG_GETARG_ARRAYTYPE_P(1);
557 num_geoms = ArrayGetNItems(ARR_NDIM(arr_geoms), ARR_DIMS(arr_geoms));
558 num_ids = ArrayGetNItems(ARR_NDIM(arr_ids), ARR_DIMS(arr_ids));
560 if ( num_geoms != num_ids )
562 elog(ERROR,
"size of geometry[] and integer[] arrays must match");
569 iter_geoms = array_create_iterator(arr_geoms, 0, NULL);
570 iter_ids = array_create_iterator(arr_ids, 0, NULL);
572 while( array_iterate(iter_geoms, &val_geom, &null_geom) &&
573 array_iterate(iter_ids, &val_id, &null_id) )
578 if ( null_geom || null_id )
580 elog(NOTICE,
"ST_AsTWKB skipping NULL entry at position %d", i);
585 uid = DatumGetInt64(val_id);
595 idlist = palloc0(num_geoms *
sizeof(int64_t));
601 elog(ERROR,
"Geometries have different dimensionality");
602 PG_FREE_IF_COPY(arr_geoms, 0);
603 PG_FREE_IF_COPY(arr_ids, 1);
614 is_homogeneous =
false;
622 array_free_iterator(iter_geoms);
623 array_free_iterator(iter_ids);
627 elog(NOTICE,
"No valid geometry - id pairs found");
628 PG_FREE_IF_COPY(arr_geoms, 0);
629 PG_FREE_IF_COPY(arr_ids, 1);
632 if ( is_homogeneous )
641 if ( PG_NARGS() > 2 && ! PG_ARGISNULL(2) )
642 sp.precision_xy = PG_GETARG_INT32(2);
645 if ( PG_NARGS() > 3 && ! PG_ARGISNULL(3) )
646 sp.precision_z = PG_GETARG_INT32(3);
649 if ( PG_NARGS() > 4 && ! PG_ARGISNULL(4) )
650 sp.precision_m = PG_GETARG_INT32(4);
656 if ( PG_NARGS() > 5 && ! PG_ARGISNULL(5) && PG_GETARG_BOOL(5) )
660 if ( PG_NARGS() > 6 && ! PG_ARGISNULL(6) && PG_GETARG_BOOL(6) )
666 sp.precision_xy, sp.precision_z, sp.precision_m,
670 result = palloc(twkb_size + VARHDRSZ);
671 memcpy(VARDATA(result), twkb, twkb_size);
672 SET_VARSIZE(result, twkb_size + VARHDRSZ);
678 PG_FREE_IF_COPY(arr_geoms, 0);
679 PG_FREE_IF_COPY(arr_ids, 1);
681 PG_RETURN_BYTEA_P(result);
697 PG_FREE_IF_COPY(geom, 0);
698 PG_RETURN_POINTER(result);
709 PG_RETURN_POINTER(geom);
718 elog(ERROR,
"%s",
string);
730 text *wkt_text = PG_GETARG_TEXT_P(0);
738 result = DirectFunctionCall1(
LWGEOM_in, CStringGetDatum(wkt));
744 PG_RETURN_DATUM(result);
758 StringInfo buf = (StringInfo) PG_GETARG_POINTER(0);
759 int32 geom_typmod = -1;
763 if ( (PG_NARGS()>2) && (!PG_ARGISNULL(2)) ) {
764 geom_typmod = PG_GETARG_INT32(2);
773 buf->cursor = buf->len;
778 if ( geom_typmod >= 0 )
781 POSTGIS_DEBUG(3,
"typmod and geometry were consistent");
785 POSTGIS_DEBUG(3,
"typmod was -1");
789 PG_RETURN_POINTER(geom);
797 POSTGIS_DEBUG(2,
"LWGEOM_send called");
810 POSTGIS_DEBUG(2,
"LWGEOM_to_bytea called");
825 POSTGIS_DEBUG(2,
"LWGEOM_from_bytea start");
827 result = (
GSERIALIZED *)DatumGetPointer(DirectFunctionCall1(
830 PG_RETURN_POINTER(result);
GSERIALIZED * postgis_valid_typmod(GSERIALIZED *gser, int32_t typmod)
Check the consistency of the metadata we want to enforce in the typmod: srid, type and dimensionality...
int gserialized_has_bbox(const GSERIALIZED *g)
Check if a GSERIALIZED has a bounding box without deserializing first.
int32_t gserialized_get_srid(const GSERIALIZED *g)
Extract the SRID from the serialized form (it is packed into three bytes so this is a handy function)...
LWGEOM * lwgeom_from_gserialized(const GSERIALIZED *g)
Allocate a new LWGEOM from a GSERIALIZED.
uint32_t gserialized_get_type(const GSERIALIZED *g)
Extract the geometry type from the serialized form (it hides in the anonymous data area,...
GSERIALIZED * gserialized_drop_gbox(GSERIALIZED *g)
Remove the bounding box from a GSERIALIZED.
#define LW_PARSER_CHECK_ALL
LWGEOM * lwcollection_as_lwgeom(const LWCOLLECTION *obj)
uint32_t lwtype_get_collectiontype(uint8_t type)
Given an lwtype number, what homogeneous collection can hold it?
int32_t lwgeom_get_srid(const LWGEOM *geom)
Return SRID number.
uint8_t * lwgeom_to_twkb_with_idlist(const LWGEOM *geom, int64_t *idlist, uint8_t variant, int8_t precision_xy, int8_t precision_z, int8_t precision_m, size_t *twkb_size)
Convert LWGEOM to a char* in TWKB format.
void lwgeom_set_srid(LWGEOM *geom, int32_t srid)
Set the SRID on an LWGEOM For collections, only the parent gets an SRID, all the children get SRID_UN...
void lwgeom_free(LWGEOM *geom)
#define LW_PARSER_CHECK_NONE
LWGEOM * lwgeom_from_twkb(const uint8_t *twkb, size_t twkb_size, char check)
WKB inputs must have a declared size, to prevent malformed WKB from reading off the end of the memory...
void lwgeom_parser_result_init(LWGEOM_PARSER_RESULT *parser_result)
int lwgeom_parse_wkt(LWGEOM_PARSER_RESULT *parser_result, char *wktstr, int parse_flags)
Parse a WKT geometry string into an LWGEOM structure.
uint8_t * lwgeom_to_wkb(const LWGEOM *geom, uint8_t variant, size_t *size_out)
Convert LWGEOM to a char* in WKB format.
uint8_t * bytes_from_hexbytes(const char *hexbuf, size_t hexsize)
int lwgeom_needs_bbox(const LWGEOM *geom)
Check whether or not a lwgeom is big enough to warrant a bounding box.
int lwgeom_has_z(const LWGEOM *geom)
Return LW_TRUE if geometry has Z ordinates.
#define POINTTYPE
LWTYPE numbers, used internally by PostGIS.
int lwgeom_isfinite(const LWGEOM *lwgeom)
Check if a LWGEOM has any non-finite (NaN or Inf) coordinates.
char * lwgeom_to_hexwkb(const LWGEOM *geom, uint8_t variant, size_t *size_out)
LWCOLLECTION * lwcollection_construct_empty(uint8_t type, int32_t srid, char hasz, char hasm)
char * lwpoint_to_latlon(const LWPOINT *p, const char *format)
const char * lwtype_name(uint8_t type)
Return the type name string associated with a type number (e.g.
void lwcollection_free(LWCOLLECTION *col)
#define TWKB_DEFAULT_PRECISION
LWCOLLECTION * lwcollection_add_lwgeom(LWCOLLECTION *col, const LWGEOM *geom)
Appends geom to the collection managed by col.
LWGEOM * lwgeom_from_wkb(const uint8_t *wkb, const size_t wkb_size, const char check)
WKB inputs must have a declared size, to prevent malformed WKB from reading off the end of the memory...
int lwgeom_has_m(const LWGEOM *geom)
Return LW_TRUE if geometry has M ordinates.
void lwgeom_add_bbox(LWGEOM *lwgeom)
Compute a bbox if not already computed.
void lwgeom_parser_result_free(LWGEOM_PARSER_RESULT *parser_result)
uint8_t * lwgeom_to_twkb(const LWGEOM *geom, uint8_t variant, int8_t precision_xy, int8_t precision_z, int8_t precision_m, size_t *twkb_size)
This library is the generic geometry handling section of PostGIS.
Datum WKBFromLWGEOM(PG_FUNCTION_ARGS)
Datum parse_WKT_lwgeom(PG_FUNCTION_ARGS)
Datum TWKBFromLWGEOM(PG_FUNCTION_ARGS)
PG_FUNCTION_INFO_V1(LWGEOM_in)
Datum LWGEOM_in(PG_FUNCTION_ARGS)
Datum LWGEOM_to_latlon(PG_FUNCTION_ARGS)
Datum LWGEOM_out(PG_FUNCTION_ARGS)
Datum LWGEOM_to_bytea(PG_FUNCTION_ARGS)
Datum LWGEOMFromTWKB(PG_FUNCTION_ARGS)
Datum LWGEOM_from_bytea(PG_FUNCTION_ARGS)
Datum TWKBFromLWGEOMArray(PG_FUNCTION_ARGS)
Datum LWGEOM_send(PG_FUNCTION_ARGS)
Datum LWGEOM_dropBBOX(PG_FUNCTION_ARGS)
Datum LWGEOM_to_text(PG_FUNCTION_ARGS)
Datum LWGEOMFromEWKB(PG_FUNCTION_ARGS)
Datum LWGEOM_addBBOX(PG_FUNCTION_ARGS)
void elog_ERROR(const char *string)
Datum LWGEOM_asHEXEWKB(PG_FUNCTION_ARGS)
Datum LWGEOM_recv(PG_FUNCTION_ARGS)
static uint32_t lwgeom_get_type(const LWGEOM *geom)
Return LWTYPE number.
char * text_to_cstring(const text *textptr)
GSERIALIZED * geometry_serialize(LWGEOM *lwgeom)
Parser result structure: returns the result of attempting to convert (E)WKT/(E)WKB to LWGEOM.