ptarray_filterm()

static POINTARRAY* ptarray_filterm ( POINTARRAY *  pa, double  min, double  max, int  returnm  )

static

Definition at line 31 of file lwmval.c.

{
        LWDEBUGF(2, "Entered %s", __func__);
 
        /*Check if M exists
        * This should already be tested earlier so we don't handle it properly.
        * If this happens it is because the function is used in another context than filterM
        and we throw an error*/
        if(!FLAGS_GET_M(pa->flags))
                lwerror("missing m-value in function %s\n",__func__);
 
        /*Dimensions in input geometry*/
        int ndims = FLAGS_NDIMS(pa->flags);
 
        /*Dimensions in result*/
        int res_ndims;
        if(returnm)
                res_ndims = ndims;
        else
                res_ndims = ndims-1;
 
        int pointsize = res_ndims * sizeof(double);
 
        double m_val;
 
        //M-value will always be the last dimension
        int m_pos = ndims-1;
 
        uint32_t i, counter=0;
        for(i=0;i<pa->npoints;i++)
        {
                m_val = *((double*)pa->serialized_pointlist + i*ndims + m_pos);
                if(m_val >= min && m_val <= max)
                        counter++;
        }
 
        POINTARRAY *pa_res = ptarray_construct(FLAGS_GET_Z(pa->flags),returnm * FLAGS_GET_M(pa->flags), counter);
 
        double *res_cursor = (double*) pa_res->serialized_pointlist;
        for(i=0;i<pa->npoints;i++)
        {
                m_val = *((double*)pa->serialized_pointlist + i*ndims + m_pos);
                if(m_val >= min && m_val <= max)
                {
                        memcpy(res_cursor, (double*) pa->serialized_pointlist + i*ndims, pointsize);
                        res_cursor+=res_ndims;
                }
        }
 
        return pa_res;
 
}

References POINTARRAY::flags, FLAGS_GET_M, FLAGS_GET_Z, FLAGS_NDIMS, LWDEBUGF, lwerror(), POINTARRAY::npoints, ptarray_construct(), and POINTARRAY::serialized_pointlist.

Referenced by lwline_filterm(), lwpoint_filterm(), and lwpoly_filterm().

Here is the call graph for this function:

Here is the caller graph for this function: