LWGEOM_angle()

Datum LWGEOM_angle

(

PG_FUNCTION_ARGS

)

Definition at line 2451 of file lwgeom_functions_basic.c.

{
        GSERIALIZED * seri_geoms[4];
        LWGEOM *geom_unser;
        LWPOINT *lwpoint;
        POINT2D points[4];
        double az1,az2 ;
        double result;
        int srids[4];
        int i = 0 ;
        int j = 0;
        int err_code = 0;
        int n_args = PG_NARGS();
 
        /* no deserialize, checking for common error first*/
        for(i=0; i<n_args; i++)
        {
                seri_geoms[i] = PG_GETARG_GSERIALIZED_P(i);
                if (gserialized_is_empty(seri_geoms[i]) )
                {/* empty geom */
                        if (i==3)
                        {
                                n_args = 3 ;
                        }
                        else
                        {
                                err_code = 1 ;
                                break ;
                        }
                } else
                {
                        if(gserialized_get_type(seri_geoms[i]) != POINTTYPE)
                        {/* geom type */
                                err_code = 2 ;
                                break;
                        }
                        else
                        {
                                srids[i] = gserialized_get_srid(seri_geoms[i]) ;
                                if(srids[0] != srids[i])
                                {/* error on srid*/
                                        err_code = 3 ;
                                        break;
                                }
                        }
                }
        }
        if (err_code >0)
                switch (err_code){
                        default: /*always executed*/
                        for (j=0;j<=i;j++)
                                PG_FREE_IF_COPY(seri_geoms[j], j);
                        /*FALLTHROUGH*/
                        case 1:
                        lwpgerror("Empty geometry");
                        PG_RETURN_NULL() ;
                        break;
 
                        case 2:
                        lwpgerror("Argument must be POINT geometries");
                        PG_RETURN_NULL();
                        break;
 
                        case 3:
                        lwpgerror("Operation on mixed SRID geometries");
                        PG_RETURN_NULL();
                        break;
                }
        /* extract points */
        for(i=0; i<n_args; i++)
        {
                geom_unser = lwgeom_from_gserialized(seri_geoms[i]) ;
                lwpoint = lwgeom_as_lwpoint(geom_unser);
                if (!lwpoint)
                {
                        for (j=0;j<n_args;j++)
                                PG_FREE_IF_COPY(seri_geoms[j], j);
                        lwpgerror("Error unserializing geometry");
                        PG_RETURN_NULL() ;
                }
 
                if ( ! getPoint2d_p(lwpoint->point, 0, &points[i]) )
                {
                        /* // can't free serialized geom, it might be needed by lw
                        for (j=0;j<n_args;j++)
                                PG_FREE_IF_COPY(seri_geoms[j], j); */
                        lwpgerror("Error extracting point");
                        PG_RETURN_NULL();
                }
                /* lwfree(geom_unser);don't do, lw may rely on this memory
                lwpoint_free(lwpoint); dont do , this memory is needed ! */
        }
        /* // can't free serialized geom, it might be needed by lw
        for (j=0;j<n_args;j++)
                PG_FREE_IF_COPY(seri_geoms[j], j); */
 
        /* compute azimuth for the 2 pairs of points
         * note that angle is not defined identically for 3 points or 4 points*/
        if (n_args == 3)
        {/* we rely on azimuth to complain if points are identical */
                if ( ! azimuth_pt_pt(&points[0], &points[1], &az1) )
                        PG_RETURN_NULL();
                if ( ! azimuth_pt_pt(&points[2], &points[1], &az2) )
                        PG_RETURN_NULL();
        } else
        {
                if ( ! azimuth_pt_pt(&points[0], &points[1], &az1) )
                        PG_RETURN_NULL();
                if ( ! azimuth_pt_pt(&points[2], &points[3], &az2) )
                        PG_RETURN_NULL();
        }
        result = az2-az1 ;
        result += (result<0) * 2 * M_PI ; /* we dont want negative angle*/
        PG_RETURN_FLOAT8(result);
}

References azimuth_pt_pt(), getPoint2d_p(), gserialized_get_srid(), gserialized_get_type(), gserialized_is_empty(), lwgeom_as_lwpoint(), lwgeom_from_gserialized(), LWPOINT::point, and POINTTYPE.

Here is the call graph for this function: