GeneratePointGeometry()

int GeneratePointGeometry	(	SHPLOADERSTATE *	state,
		SHPObject *	obj,
		char **	geometry,
		int	force_multi
	)

Generate an allocated geometry string for shapefile object obj using the state parameters if "force_multi" is true, single points will instead be created as multipoints with a single vertice.

Definition at line 227 of file shp2pgsql-core.c.

References shp_loader_state::config, FLAGS_SET_M, FLAGS_SET_Z, shp_loader_state::from_srid, shp_loader_state::has_m, shp_loader_state::has_z, LW_TRUE, lwcollection_as_lwgeom(), lwcollection_construct(), lwfree(), lwgeom_free(), lwgeom_to_hexwkb(), lwgeom_to_wkt(), lwpoint_as_lwgeom(), lwpoint_construct(), lwpoint_construct_empty(), POINT4D::m, malloc(), shp_loader_state::message, MULTIPOINTTYPE, SHPObject::nVertices, SHPObject::padfM, SHPObject::padfX, SHPObject::padfY, SHPObject::padfZ, ptarray_append_point(), ptarray_construct_empty(), SHPLOADERERR, SHPLOADERMSGLEN, SHPLOADEROK, shp_loader_config::use_wkt, WKB_EXTENDED, WKT_EXTENDED, WKT_PRECISION, POINT4D::x, POINT4D::y, and POINT4D::z.

Referenced by ShpLoaderGenerateSQLRowStatement().

{
        LWGEOM **lwmultipoints;
        LWGEOM *lwgeom = NULL;

        POINT4D point4d;

        int dims = 0;
        int u;

        char *mem;
        size_t mem_length;

        FLAGS_SET_Z(dims, state->has_z);
        FLAGS_SET_M(dims, state->has_m);

        /* POINT EMPTY encoded as POINT(NaN NaN) */
        if (obj->nVertices == 1 && isnan(obj->padfX[0]) && isnan(obj->padfY[0]))
        {
                lwgeom = lwpoint_as_lwgeom(lwpoint_construct_empty(state->from_srid, state->has_z, state->has_m));
        }
        /* Not empty */
        else
        {
                /* Allocate memory for our array of LWPOINTs and our dynptarrays */
                lwmultipoints = malloc(sizeof(LWPOINT *) * obj->nVertices);

                /* We need an array of pointers to each of our sub-geometries */
                for (u = 0; u < obj->nVertices; u++)
                {
                        /* Create a ptarray containing a single point */
                        POINTARRAY *pa = ptarray_construct_empty(state->has_z, state->has_m, 1);

                        /* Generate the point */
                        point4d.x = obj->padfX[u];
                        point4d.y = obj->padfY[u];

                        if (state->has_z)
                                point4d.z = obj->padfZ[u];
                        if (state->has_m)
                                point4d.m = obj->padfM[u];

                        /* Add in the point! */
                        ptarray_append_point(pa, &point4d, LW_TRUE);

                        /* Generate the LWPOINT */
                        lwmultipoints[u] = lwpoint_as_lwgeom(lwpoint_construct(state->from_srid, NULL, pa));
                }

                /* If we have more than 1 vertex then we are working on a MULTIPOINT and so generate a MULTIPOINT
                rather than a POINT */
                if ((obj->nVertices > 1) || force_multi)
                {
                        lwgeom = lwcollection_as_lwgeom(lwcollection_construct(MULTIPOINTTYPE, state->from_srid, NULL, obj->nVertices, lwmultipoints));
                }
                else
                {
                        lwgeom = lwmultipoints[0];
                        lwfree(lwmultipoints);
                }
        }

        if (state->config->use_wkt)
        {
                mem = lwgeom_to_wkt(lwgeom, WKT_EXTENDED, WKT_PRECISION, &mem_length);
        }
        else
        {
                mem = lwgeom_to_hexwkb(lwgeom, WKB_EXTENDED, &mem_length);
        }

        if ( !mem )
        {
                snprintf(state->message, SHPLOADERMSGLEN, "unable to write geometry");
                return SHPLOADERERR;
        }

        /* Free all of the allocated items */
        lwgeom_free(lwgeom);

        /* Return the string - everything ok */
        *geometry = mem;

        return SHPLOADEROK;
}

Here is the call graph for this function:

Here is the caller graph for this function: