www.pudn.com > mitab.rar > ogrmultilinestring.cpp
/******************************************************************************
* $Id: ogrmultilinestring.cpp,v 1.13 2004/02/21 15:36:14 warmerda Exp $
*
* Project: OpenGIS Simple Features Reference Implementation
* Purpose: The OGRMultiLineString class.
* Author: Frank Warmerdam, warmerda@home.com
*
******************************************************************************
* Copyright (c) 1999, Frank Warmerdam
*
* Permission is hereby granted, free of charge, to any person obtaining a
* copy of this software and associated documentation files (the "Software"),
* to deal in the Software without restriction, including without limitation
* the rights to use, copy, modify, merge, publish, distribute, sublicense,
* and/or sell copies of the Software, and to permit persons to whom the
* Software is furnished to do so, subject to the following conditions:
*
* The above copyright notice and this permission notice shall be included
* in all copies or substantial portions of the Software.
*
* THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS
* OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
* FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL
* THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
* LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING
* FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER
* DEALINGS IN THE SOFTWARE.
******************************************************************************
*
* $Log: ogrmultilinestring.cpp,v $
* Revision 1.13 2004/02/21 15:36:14 warmerda
* const correctness updates for geometry: bug 289
*
* Revision 1.12 2004/01/16 21:57:16 warmerda
* fixed up EMPTY support
*
* Revision 1.11 2004/01/16 21:20:00 warmerda
* Added EMPTY support
*
* Revision 1.10 2003/09/11 22:47:54 aamici
* add class constructors and destructors where needed in order to
* let the mingw/cygwin binutils produce sensible partially linked objet files
* with 'ld -r'.
*
* Revision 1.9 2003/05/28 19:16:43 warmerda
* fixed up argument names and stuff for docs
*
* Revision 1.8 2002/09/11 13:47:17 warmerda
* preliminary set of fixes for 3D WKB enum
*
* Revision 1.7 2001/07/19 18:25:07 warmerda
* expanded tabs
*
* Revision 1.6 2001/07/18 05:03:05 warmerda
* added CPL_CVSID
*
* Revision 1.5 2001/05/24 18:05:18 warmerda
* substantial fixes to WKT support for MULTIPOINT/LINE/POLYGON
*
* Revision 1.4 1999/12/21 05:45:28 warmerda
* Fixed to check for wkbLineString instead of wkbPoint.
*
* Revision 1.3 1999/11/18 19:02:19 warmerda
* expanded tabs
*
* Revision 1.2 1999/06/25 20:44:43 warmerda
* implemented assignSpatialReference, carry properly
*
* Revision 1.1 1999/05/31 20:42:06 warmerda
* New
*
*/
#include "ogr_geometry.h"
#include "ogr_p.h"
CPL_CVSID("$Id: ogrmultilinestring.cpp,v 1.13 2004/02/21 15:36:14 warmerda Exp $");
/************************************************************************/
/* OGRMultiLineString() */
/************************************************************************/
OGRMultiLineString::OGRMultiLineString()
{
}
/************************************************************************/
/* ~OGRMultiLineString() */
/************************************************************************/
OGRMultiLineString::~OGRMultiLineString()
{
}
/************************************************************************/
/* getGeometryType() */
/************************************************************************/
OGRwkbGeometryType OGRMultiLineString::getGeometryType() const
{
if( getCoordinateDimension() == 3 )
return wkbMultiLineString25D;
else
return wkbMultiLineString;
}
/************************************************************************/
/* getGeometryName() */
/************************************************************************/
const char * OGRMultiLineString::getGeometryName() const
{
return "MULTILINESTRING";
}
/************************************************************************/
/* addGeometryDirectly() */
/************************************************************************/
OGRErr OGRMultiLineString::addGeometryDirectly( OGRGeometry * poNewGeom )
{
if( poNewGeom->getGeometryType() != wkbLineString
&& poNewGeom->getGeometryType() != wkbLineString25D )
return OGRERR_UNSUPPORTED_GEOMETRY_TYPE;
return OGRGeometryCollection::addGeometryDirectly( poNewGeom );
}
/************************************************************************/
/* clone() */
/************************************************************************/
OGRGeometry *OGRMultiLineString::clone() const
{
OGRMultiLineString *poNewGC;
poNewGC = new OGRMultiLineString;
poNewGC->assignSpatialReference( getSpatialReference() );
for( int i = 0; i < getNumGeometries(); i++ )
{
poNewGC->addGeometry( getGeometryRef(i) );
}
return poNewGC;
}
/************************************************************************/
/* importFromWkt() */
/* */
/* Instantiate from well known text format. Currently this is */
/* `MULTILINESTRING ((x y, x y, ...),(x y, ...),...)'. */
/************************************************************************/
OGRErr OGRMultiLineString::importFromWkt( char ** ppszInput )
{
char szToken[OGR_WKT_TOKEN_MAX];
const char *pszInput = *ppszInput;
OGRErr eErr;
/* -------------------------------------------------------------------- */
/* Clear existing rings. */
/* -------------------------------------------------------------------- */
empty();
/* -------------------------------------------------------------------- */
/* Read and verify the ``MULTILINESTRING'' keyword token. */
/* -------------------------------------------------------------------- */
pszInput = OGRWktReadToken( pszInput, szToken );
if( !EQUAL(szToken,getGeometryName()) )
return OGRERR_CORRUPT_DATA;
/* -------------------------------------------------------------------- */
/* The next character should be a ( indicating the start of the */
/* list of linestrings. */
/* -------------------------------------------------------------------- */
pszInput = OGRWktReadToken( pszInput, szToken );
if( szToken[0] != '(' )
return OGRERR_CORRUPT_DATA;
/* -------------------------------------------------------------------- */
/* If the next token is EMPTY, then verify that we have proper */
/* EMPTY format will a trailing closing bracket. */
/* -------------------------------------------------------------------- */
OGRWktReadToken( pszInput, szToken );
if( EQUAL(szToken,"EMPTY") )
{
pszInput = OGRWktReadToken( pszInput, szToken );
pszInput = OGRWktReadToken( pszInput, szToken );
*ppszInput = (char *) pszInput;
if( !EQUAL(szToken,")") )
return OGRERR_CORRUPT_DATA;
else
return OGRERR_NONE;
}
/* ==================================================================== */
/* Read each line in turn. Note that we try to reuse the same */
/* point list buffer from ring to ring to cut down on */
/* allocate/deallocate overhead. */
/* ==================================================================== */
OGRRawPoint *paoPoints = NULL;
int nMaxPoints = 0;
double *padfZ = NULL;
do
{
int nPoints = 0;
/* -------------------------------------------------------------------- */
/* Read points for one line from input. */
/* -------------------------------------------------------------------- */
pszInput = OGRWktReadPoints( pszInput, &paoPoints, &padfZ, &nMaxPoints,
&nPoints );
if( pszInput == NULL )
{
eErr = OGRERR_CORRUPT_DATA;
break;
}
/* -------------------------------------------------------------------- */
/* Create the new line, and add to collection. */
/* -------------------------------------------------------------------- */
OGRLineString *poLine;
poLine = new OGRLineString();
poLine->setPoints( nPoints, paoPoints, padfZ );
eErr = addGeometryDirectly( poLine );
/* -------------------------------------------------------------------- */
/* Read the delimeter following the ring. */
/* -------------------------------------------------------------------- */
pszInput = OGRWktReadToken( pszInput, szToken );
} while( szToken[0] == ',' && eErr == OGRERR_NONE );
/* -------------------------------------------------------------------- */
/* freak if we don't get a closing bracket. */
/* -------------------------------------------------------------------- */
CPLFree( paoPoints );
CPLFree( padfZ );
if( eErr != OGRERR_NONE )
return eErr;
if( szToken[0] != ')' )
return OGRERR_CORRUPT_DATA;
*ppszInput = (char *) pszInput;
return OGRERR_NONE;
}
/************************************************************************/
/* exportToWkt() */
/* */
/* Translate this structure into it's well known text format */
/* equivelent. This could be made alot more CPU efficient! */
/************************************************************************/
OGRErr OGRMultiLineString::exportToWkt( char ** ppszDstText ) const
{
char **papszLines;
int iLine, nCumulativeLength = 0;
OGRErr eErr;
if( getNumGeometries() == 0 )
{
*ppszDstText = CPLStrdup("MULTILINESTRING(EMPTY)");
return OGRERR_NONE;
}
/* -------------------------------------------------------------------- */
/* Build a list of strings containing the stuff for each ring. */
/* -------------------------------------------------------------------- */
papszLines = (char **) CPLCalloc(sizeof(char *),getNumGeometries());
for( iLine = 0; iLine < getNumGeometries(); iLine++ )
{
eErr = getGeometryRef(iLine)->exportToWkt( &(papszLines[iLine]) );
if( eErr != OGRERR_NONE )
return eErr;
CPLAssert( EQUALN(papszLines[iLine],"LINESTRING (", 12) );
nCumulativeLength += strlen(papszLines[iLine] + 11);
}
/* -------------------------------------------------------------------- */
/* Allocate exactly the right amount of space for the */
/* aggregated string. */
/* -------------------------------------------------------------------- */
*ppszDstText = (char *) VSIMalloc(nCumulativeLength+getNumGeometries()+20);
if( *ppszDstText == NULL )
return OGRERR_NOT_ENOUGH_MEMORY;
/* -------------------------------------------------------------------- */
/* Build up the string, freeing temporary strings as we go. */
/* -------------------------------------------------------------------- */
strcpy( *ppszDstText, "MULTILINESTRING (" );
for( iLine = 0; iLine < getNumGeometries(); iLine++ )
{
if( iLine > 0 )
strcat( *ppszDstText, "," );
strcat( *ppszDstText, papszLines[iLine] + 11 );
VSIFree( papszLines[iLine] );
}
strcat( *ppszDstText, ")" );
CPLFree( papszLines );
return OGRERR_NONE;
}