PostGIS  3.4.0dev-r@@SVN_REVISION@@

◆ RASTER_getPolygon()

Datum RASTER_getPolygon ( PG_FUNCTION_ARGS  )

Definition at line 980 of file rtpg_geometry.c.

981 {
982  rt_pgraster *pgraster = NULL;
983  rt_raster raster = NULL;
984  int num_bands = 0;
985  int nband = 1;
986  int err;
987  LWMPOLY *surface = NULL;
988  GSERIALIZED *rtn = NULL;
989 
990  /* raster */
991  if (PG_ARGISNULL(0))
992  PG_RETURN_NULL();
993  pgraster = (rt_pgraster *) PG_DETOAST_DATUM(PG_GETARG_DATUM(0));
994 
995  raster = rt_raster_deserialize(pgraster, FALSE);
996  if (!raster) {
997  PG_FREE_IF_COPY(pgraster, 0);
998  elog(ERROR, "RASTER_getPolygon: Could not deserialize raster");
999  PG_RETURN_NULL();
1000  }
1001 
1002  /* num_bands */
1003  num_bands = rt_raster_get_num_bands(raster);
1004  if (num_bands < 1) {
1005  elog(NOTICE, "Raster provided has no bands");
1007  PG_FREE_IF_COPY(pgraster, 0);
1008  PG_RETURN_NULL();
1009  }
1010 
1011  /* band index is 1-based */
1012  if (!PG_ARGISNULL(1))
1013  nband = PG_GETARG_INT32(1);
1014  if (nband < 1 || nband > num_bands) {
1015  elog(NOTICE, "Invalid band index (must use 1-based). Returning NULL");
1017  PG_FREE_IF_COPY(pgraster, 0);
1018  PG_RETURN_NULL();
1019  }
1020 
1021  /* get band surface */
1022  err = rt_raster_surface(raster, nband - 1, &surface);
1024  PG_FREE_IF_COPY(pgraster, 0);
1025 
1026  if (err != ES_NONE) {
1027  elog(ERROR, "RASTER_getPolygon: Could not get raster band's surface");
1028  PG_RETURN_NULL();
1029  }
1030  else if (surface == NULL) {
1031  elog(NOTICE, "Raster is empty or all pixels of band are NODATA. Returning NULL");
1032  PG_RETURN_NULL();
1033  }
1034 
1035  rtn = geometry_serialize(lwmpoly_as_lwgeom(surface));
1036  lwmpoly_free(surface);
1037 
1038  PG_RETURN_POINTER(rtn);
1039 }
#define FALSE
Definition: dbfopen.c:72
void lwmpoly_free(LWMPOLY *mpoly)
Definition: lwmpoly.c:53
LWGEOM * lwmpoly_as_lwgeom(const LWMPOLY *obj)
Definition: lwgeom.c:294
void rt_raster_destroy(rt_raster raster)
Release memory associated to a raster.
Definition: rt_raster.c:86
rt_errorstate rt_raster_surface(rt_raster raster, int nband, LWMPOLY **surface)
Get a raster as a surface (multipolygon).
Definition: rt_geometry.c:355
@ ES_NONE
Definition: librtcore.h:182
uint16_t rt_raster_get_num_bands(rt_raster raster)
Definition: rt_raster.c:376
rt_raster rt_raster_deserialize(void *serialized, int header_only)
Return a raster from a serialized form.
Definition: rt_serialize.c:725
nband
Definition: pixval.py:53
raster
Be careful!! Zeros function's input parameter can be a (height x width) array, not (width x height): ...
Definition: rtrowdump.py:121
Struct definitions.
Definition: librtcore.h:2403

References ES_NONE, FALSE, lwmpoly_as_lwgeom(), lwmpoly_free(), pixval::nband, rtrowdump::raster, rt_raster_deserialize(), rt_raster_destroy(), rt_raster_get_num_bands(), and rt_raster_surface().

Here is the call graph for this function: