Z/M coordinates are ignored for this calculation. Pass in non-null, geodetic bounding box for function to fill out. LWGEOM must have been built from a GSERIALIZED to provide double aligned point arrays.