Function msGMLWriteWFSQuery

mapgml.c:1347–1561

** msGMLWriteWFSQuery() ** ** Similar to msGMLWriteQuery() but tuned for use with WFS 1.0.0 */

** Similar to msGMLWriteQuery() but tuned for use with WFS 1.0.0
*/
int msGMLWriteWFSQuery(mapObj *map, FILE *stream, int startindex, int maxfeatures, 
                       char *default_namespace_prefix, int outputformat)
{
#ifdef USE_WFS_SVR
  int status;
  int i,j,k;
  layerObj *lp=NULL;
  shapeObj shape;
  rectObj  resultBounds = {-1.0,-1.0,-1.0,-1.0};
  int features = 0;
  int currentfeature =0;

  gmlGroupListObj *groupList=NULL;
  gmlItemListObj *itemList=NULL;
  gmlConstantListObj *constantList=NULL;
  gmlGeometryListObj *geometryList=NULL;
  gmlItemObj *item=NULL;
  gmlConstantObj *constant=NULL;

  char *namespace_prefix=NULL;
  const char *axis = NULL;
  int bSwapAxis = 0;
  double tmp;
  const char *srsMap =  NULL;

  msInitShape(&shape);

  /*add a check to see if the map projection is set to be north-east*/
  for( i = 0; i < map->projection.numargs; i++ )
  {
      if( strstr(map->projection.args[i],"epsgaxis=") != NULL )
      {
          axis = strstr(map->projection.args[i],"=") + 1;
          break;
      }
  }

  if (axis && strcasecmp(axis,"ne") == 0 )
    bSwapAxis = 1;

  
  /* Need to start with BBOX of the whole resultset */
  if (msGetQueryResultBounds(map, &resultBounds) > 0)
  {
      if (bSwapAxis)
      {
          tmp = resultBounds.minx;
          resultBounds.minx =  resultBounds.miny;
          resultBounds.miny = tmp;

          tmp = resultBounds.maxx;
          resultBounds.maxx =  resultBounds.maxy;
          resultBounds.maxy = tmp;

      }
      srsMap = msOWSGetEPSGProj(&(map->projection), NULL, "FGO", MS_TRUE);
      if (!srsMap)
        msOWSGetEPSGProj(&(map->projection), &(map->web.metadata), "FGO", MS_TRUE);