OGRErr OGRILI1Layer::CreateFeature( OGRFeature *poFeature ) {
static long tid = -1; //system generated TID (must be unique within table)
VSIFPrintf( poDS->GetTransferFile(), "OBJE" );
if ( !EQUAL(poFeatureDefn->GetFieldDefn(0)->GetNameRef(), "TID") )
{
//Input is not generated from an Interlis 1 source
if (poFeature->GetFID() != OGRNullFID)
tid = poFeature->GetFID();
else
++tid;
VSIFPrintf( poDS->GetTransferFile(), " %ld", tid );
//Embedded geometry
if( poFeature->GetGeometryRef() != NULL )
{
OGRGeometry *poGeometry = poFeature->GetGeometryRef();
// 2D Point
if( poGeometry->getGeometryType() == wkbPoint )
{
OGRPoint *poPoint = (OGRPoint *) poGeometry;
VSIFPrintf( poDS->GetTransferFile(), " %s", d2str(poPoint->getX()) );
VSIFPrintf( poDS->GetTransferFile(), " %s", d2str(poPoint->getY()) );
}
// 3D Point
else if( poGeometry->getGeometryType() == wkbPoint25D )
{
OGRPoint *poPoint = (OGRPoint *) poGeometry;
VSIFPrintf( poDS->GetTransferFile(), " %s", d2str(poPoint->getX()) );
VSIFPrintf( poDS->GetTransferFile(), " %s", d2str(poPoint->getY()) );
VSIFPrintf( poDS->GetTransferFile(), " %s", d2str(poPoint->getZ()) );
}
}
}
// Write all fields.
for(int iField = 0; iField < poFeatureDefn->GetFieldCount(); iField++ )
{
if ( !EQUAL(poFeatureDefn->GetFieldDefn(iField)->GetNameRef(), "ILI_Geometry") )
{
if ( poFeature->IsFieldSet( iField ) )
{
const char *pszRaw = poFeature->GetFieldAsString( iField );
VSIFPrintf( poDS->GetTransferFile(), " %s", pszRaw );
}
else
{
VSIFPrintf( poDS->GetTransferFile(), " @" );
}
}
}
VSIFPrintf( poDS->GetTransferFile(), "\n" );
// Write out Geometry
if( poFeature->GetGeometryRef() != NULL )
{
if (EQUAL(poFeatureDefn->GetFieldDefn(poFeatureDefn->GetFieldCount()-1)->GetNameRef(), "ILI_Geometry"))
{
//Write original ILI geometry
VSIFPrintf( poDS->GetTransferFile(), "%s", poFeature->GetFieldAsString( poFeatureDefn->GetFieldCount()-1 ) );
}
else
{
//Convert to ILI geometry
GeometryAppend(poFeature->GetGeometryRef());
}
}
return OGRERR_NONE;
}
int OGRILI1Layer::GeometryAppend( OGRGeometry *poGeometry )
{
/* -------------------------------------------------------------------- */
/* 2D Point */
/* -------------------------------------------------------------------- */
if( poGeometry->getGeometryType() == wkbPoint )
{
/* embedded in from non-geometry fields */
}
/* -------------------------------------------------------------------- */
/* 3D Point */
/* -------------------------------------------------------------------- */
else if( poGeometry->getGeometryType() == wkbPoint25D )
{
/* embedded in from non-geometry fields */
}
/* -------------------------------------------------------------------- */
/* LineString and LinearRing */
/* -------------------------------------------------------------------- */
else if( poGeometry->getGeometryType() == wkbLineString
|| poGeometry->getGeometryType() == wkbLineString25D )
{
AppendCoordinateList( (OGRLineString *) poGeometry, poDS );
}
/* -------------------------------------------------------------------- */
/* Polygon */
/* -------------------------------------------------------------------- */
else if( poGeometry->getGeometryType() == wkbPolygon
|| poGeometry->getGeometryType() == wkbPolygon25D )
{
OGRPolygon *poPolygon = (OGRPolygon *) poGeometry;
if( poPolygon->getExteriorRing() != NULL )
{
if( !GeometryAppend( poPolygon->getExteriorRing() ) )
return FALSE;
}
for( int iRing = 0; iRing < poPolygon->getNumInteriorRings(); iRing++ )
{
OGRLinearRing *poRing = poPolygon->getInteriorRing(iRing);
if( !GeometryAppend( poRing ) )
return FALSE;
}
}
/* -------------------------------------------------------------------- */
/* MultiPolygon */
/* -------------------------------------------------------------------- */
else if( wkbFlatten(poGeometry->getGeometryType()) == wkbMultiPolygon
|| wkbFlatten(poGeometry->getGeometryType()) == wkbMultiLineString
|| wkbFlatten(poGeometry->getGeometryType()) == wkbMultiPoint
|| wkbFlatten(poGeometry->getGeometryType()) == wkbGeometryCollection )
{
OGRGeometryCollection *poGC = (OGRGeometryCollection *) poGeometry;
int iMember;
if( wkbFlatten(poGeometry->getGeometryType()) == wkbMultiPolygon )
{
}
else if( wkbFlatten(poGeometry->getGeometryType()) == wkbMultiLineString )
{
}
else if( wkbFlatten(poGeometry->getGeometryType()) == wkbMultiPoint )
{
}
else
{
}
for( iMember = 0; iMember < poGC->getNumGeometries(); iMember++)
{
OGRGeometry *poMember = poGC->getGeometryRef( iMember );
if( !GeometryAppend( poMember ) )
return FALSE;
}
}
else
return FALSE;
return TRUE;
}
int OGRILI1Layer::GeometryAppend( OGRGeometry *poGeometry )
{
//CPLDebug( "OGR_ILI", "OGRILI1Layer::GeometryAppend OGRGeometryType: %s", OGRGeometryTypeToName(poGeometry->getGeometryType()));
/* -------------------------------------------------------------------- */
/* 2D Point */
/* -------------------------------------------------------------------- */
if( poGeometry->getGeometryType() == wkbPoint )
{
/* embedded in from non-geometry fields */
}
/* -------------------------------------------------------------------- */
/* 3D Point */
/* -------------------------------------------------------------------- */
else if( poGeometry->getGeometryType() == wkbPoint25D )
{
/* embedded in from non-geometry fields */
}
/* -------------------------------------------------------------------- */
/* LineString and LinearRing */
/* -------------------------------------------------------------------- */
else if( poGeometry->getGeometryType() == wkbLineString
|| poGeometry->getGeometryType() == wkbLineString25D )
{
AppendCoordinateList( (OGRLineString *) poGeometry, poDS );
}
/* -------------------------------------------------------------------- */
/* Polygon */
/* -------------------------------------------------------------------- */
else if( poGeometry->getGeometryType() == wkbPolygon
|| poGeometry->getGeometryType() == wkbPolygon25D )
{
OGRPolygon *poPolygon = (OGRPolygon *) poGeometry;
if( poPolygon->getExteriorRing() != NULL )
{
if( !GeometryAppend( poPolygon->getExteriorRing() ) )
return FALSE;
}
for( int iRing = 0; iRing < poPolygon->getNumInteriorRings(); iRing++ )
{
OGRLinearRing *poRing = poPolygon->getInteriorRing(iRing);
if( !GeometryAppend( poRing ) )
return FALSE;
}
}
/* -------------------------------------------------------------------- */
/* MultiPolygon */
/* -------------------------------------------------------------------- */
else if( wkbFlatten(poGeometry->getGeometryType()) == wkbMultiPolygon
|| wkbFlatten(poGeometry->getGeometryType()) == wkbMultiLineString
|| wkbFlatten(poGeometry->getGeometryType()) == wkbMultiPoint
|| wkbFlatten(poGeometry->getGeometryType()) == wkbGeometryCollection
|| wkbFlatten(poGeometry->getGeometryType()) == wkbMultiCurve
|| wkbFlatten(poGeometry->getGeometryType()) == wkbMultiCurveZ )
{
OGRGeometryCollection *poGC = (OGRGeometryCollection *) poGeometry;
int iMember;
if( wkbFlatten(poGeometry->getGeometryType()) == wkbMultiPolygon )
{
}
else if( wkbFlatten(poGeometry->getGeometryType()) == wkbMultiLineString )
{
}
else if( wkbFlatten(poGeometry->getGeometryType()) == wkbMultiPoint )
{
}
else
{
}
for( iMember = 0; iMember < poGC->getNumGeometries(); iMember++)
{
OGRGeometry *poMember = poGC->getGeometryRef( iMember );
if( !GeometryAppend( poMember ) )
return FALSE;
}
}
else if( poGeometry->getGeometryType() == wkbCompoundCurve
|| poGeometry->getGeometryType() == wkbCompoundCurveZ )
{
AppendCoumpoundCurve( ( OGRCompoundCurve *) poGeometry, poDS );
} else {
CPLError(CE_Warning, CPLE_AppDefined, "Skipping unknown geometry type '%s'", OGRGeometryTypeToName(poGeometry->getGeometryType()));
return FALSE;
}
return TRUE;
}
OGRErr OGRILI1Layer::ICreateFeature( OGRFeature *poFeature ) {
static long tid = -1; //system generated TID (must be unique within table)
VSIFPrintf( poDS->GetTransferFile(), "OBJE" );
if ( poFeatureDefn->GetFieldCount() && !EQUAL(poFeatureDefn->GetFieldDefn(0)->GetNameRef(), "TID") )
{
//Input is not generated from an Interlis 1 source
if (poFeature->GetFID() != OGRNullFID)
tid = (int)poFeature->GetFID();
else
++tid;
VSIFPrintf( poDS->GetTransferFile(), " %ld", tid );
//Embedded geometry
if( poFeature->GetGeometryRef() != NULL )
{
OGRGeometry *poGeometry = poFeature->GetGeometryRef();
// 2D Point
if( poGeometry->getGeometryType() == wkbPoint )
{
OGRPoint *poPoint = (OGRPoint *) poGeometry;
VSIFPrintf( poDS->GetTransferFile(), " %s", d2str(poPoint->getX()) );
VSIFPrintf( poDS->GetTransferFile(), " %s", d2str(poPoint->getY()) );
}
// 3D Point
else if( poGeometry->getGeometryType() == wkbPoint25D )
{
OGRPoint *poPoint = (OGRPoint *) poGeometry;
VSIFPrintf( poDS->GetTransferFile(), " %s", d2str(poPoint->getX()) );
VSIFPrintf( poDS->GetTransferFile(), " %s", d2str(poPoint->getY()) );
VSIFPrintf( poDS->GetTransferFile(), " %s", d2str(poPoint->getZ()) );
}
}
}
// Write all fields.
for(int iField = 0; iField < poFeatureDefn->GetFieldCount(); iField++ )
{
if ( poFeature->IsFieldSet( iField ) )
{
const char *pszRaw = poFeature->GetFieldAsString( iField );
if (poFeatureDefn->GetFieldDefn( iField )->GetType() == OFTString) {
//Interlis 1 encoding is ISO 8859-1 (Latin1) -> Recode from UTF-8
char* pszString = CPLRecode(pszRaw, CPL_ENC_UTF8, CPL_ENC_ISO8859_1);
//Replace spaces
for(size_t i=0; i<strlen(pszString); i++ ) {
if (pszString[i] == ' ') pszString[i] = '_';
}
VSIFPrintf( poDS->GetTransferFile(), " %s", pszString );
CPLFree( pszString );
} else {
VSIFPrintf( poDS->GetTransferFile(), " %s", pszRaw );
}
}
else
{
VSIFPrintf( poDS->GetTransferFile(), " @" );
}
}
VSIFPrintf( poDS->GetTransferFile(), "\n" );
// Write out Geometry
if( poFeature->GetGeometryRef() != NULL )
{
GeometryAppend(poFeature->GetGeometryRef());
}
return OGRERR_NONE;
}
int OGRILI1Layer::GeometryAppend( OGRGeometry *poGeometry )
{
#ifdef DEBUG_VERBOSE
CPLDebug( "OGR_ILI", "OGRILI1Layer::GeometryAppend OGRGeometryType: %s",
OGRGeometryTypeToName(poGeometry->getGeometryType()) );
#endif
/* -------------------------------------------------------------------- */
/* 2D Point */
/* -------------------------------------------------------------------- */
if( poGeometry->getGeometryType() == wkbPoint )
{
/* embedded in from non-geometry fields */
}
/* -------------------------------------------------------------------- */
/* 3D Point */
/* -------------------------------------------------------------------- */
else if( poGeometry->getGeometryType() == wkbPoint25D )
{
/* embedded in from non-geometry fields */
}
/* -------------------------------------------------------------------- */
/* LineString and LinearRing */
/* -------------------------------------------------------------------- */
else if( poGeometry->getGeometryType() == wkbLineString
|| poGeometry->getGeometryType() == wkbLineString25D )
{
AppendCoordinateList( poGeometry->toLineString(),
poDS );
}
/* -------------------------------------------------------------------- */
/* Polygon */
/* -------------------------------------------------------------------- */
else if( poGeometry->getGeometryType() == wkbPolygon
|| poGeometry->getGeometryType() == wkbPolygon25D )
{
OGRPolygon *poPolygon = poGeometry->toPolygon();
for( auto&& poRing: poPolygon )
{
if( !GeometryAppend( poRing ) )
return FALSE;
}
}
/* -------------------------------------------------------------------- */
/* MultiPolygon */
/* -------------------------------------------------------------------- */
else if( wkbFlatten(poGeometry->getGeometryType()) == wkbMultiPolygon
|| wkbFlatten(poGeometry->getGeometryType()) == wkbMultiLineString
|| wkbFlatten(poGeometry->getGeometryType()) == wkbMultiPoint
|| wkbFlatten(poGeometry->getGeometryType()) == wkbGeometryCollection
|| wkbFlatten(poGeometry->getGeometryType()) == wkbMultiCurve
|| wkbFlatten(poGeometry->getGeometryType()) == wkbMultiCurveZ )
{
OGRGeometryCollection *poGC = poGeometry->toGeometryCollection();
#if 0
// TODO: Why this large NOP block?
if( wkbFlatten(poGeometry->getGeometryType()) == wkbMultiPolygon )
{
}
else if( wkbFlatten(poGeometry->getGeometryType()) == wkbMultiLineString )
{
}
else if( wkbFlatten(poGeometry->getGeometryType()) == wkbMultiPoint )
{
}
else
{
}
#endif
for( auto&& poMember: poGC )
{
if( !GeometryAppend( poMember ) )
return FALSE;
}
}
else if( poGeometry->getGeometryType() == wkbCompoundCurve
|| poGeometry->getGeometryType() == wkbCompoundCurveZ )
{
AppendCoumpoundCurve( poGeometry->toCompoundCurve(),
poDS );
} else {
CPLError( CE_Warning, CPLE_AppDefined,
"Skipping unknown geometry type '%s'",
OGRGeometryTypeToName(poGeometry->getGeometryType()) );
return FALSE;
}
return TRUE;
}
