LWGEOM_azimuth()

Datum LWGEOM_azimuth

(

PG_FUNCTION_ARGS

)

Definition at line 2513 of file lwgeom_functions_basic.c.

{
        GSERIALIZED *geom;
        LWPOINT *lwpoint;
        POINT2D p1, p2;
        double result;
        int32_t srid;
 
        /* Extract first point */
        geom = PG_GETARG_GSERIALIZED_P(0);
        lwpoint = lwgeom_as_lwpoint(lwgeom_from_gserialized(geom));
        if (!lwpoint)
        {
                PG_FREE_IF_COPY(geom, 0);
                lwpgerror("Argument must be POINT geometries");
                PG_RETURN_NULL();
        }
        srid = lwpoint->srid;
        if (!getPoint2d_p(lwpoint->point, 0, &p1))
        {
                PG_FREE_IF_COPY(geom, 0);
                lwpgerror("Error extracting point");
                PG_RETURN_NULL();
        }
        lwpoint_free(lwpoint);
        PG_FREE_IF_COPY(geom, 0);
 
        /* Extract second point */
        geom = PG_GETARG_GSERIALIZED_P(1);
        lwpoint = lwgeom_as_lwpoint(lwgeom_from_gserialized(geom));
        if (!lwpoint)
        {
                PG_FREE_IF_COPY(geom, 1);
                lwpgerror("Argument must be POINT geometries");
                PG_RETURN_NULL();
        }
        if (lwpoint->srid != srid)
        {
                PG_FREE_IF_COPY(geom, 1);
                lwpgerror("Operation on mixed SRID geometries");
                PG_RETURN_NULL();
        }
        if (!getPoint2d_p(lwpoint->point, 0, &p2))
        {
                PG_FREE_IF_COPY(geom, 1);
                lwpgerror("Error extracting point");
                PG_RETURN_NULL();
        }
        lwpoint_free(lwpoint);
        PG_FREE_IF_COPY(geom, 1);
 
        /* Standard return value for equality case */
        if ((p1.x == p2.x) && (p1.y == p2.y))
        {
                PG_RETURN_NULL();
        }
 
        /* Compute azimuth */
        if (!azimuth_pt_pt(&p1, &p2, &result))
        {
                PG_RETURN_NULL();
        }
 
        PG_RETURN_FLOAT8(result);
}

References azimuth_pt_pt(), getPoint2d_p(), lwgeom_as_lwpoint(), lwgeom_from_gserialized(), lwpoint_free(), LWPOINT::point, result, LWPOINT::srid, POINT2D::x, and POINT2D::y.

Here is the call graph for this function: