/****************************************************************************** * $Id: ogrdodslayer.cpp 10646 2007-01-18 02:38:10Z warmerdam $ * * Project:  OGR/DODS Interface * Purpose:  Implements OGRDODSLayer class. * Author:   Frank Warmerdam, warmerdam@pobox.com * ****************************************************************************** * Copyright (c) 2004, Frank Warmerdam <warmerdam@pobox.com> * * 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. ****************************************************************************/#include "ogr_dods.h"#include "cpl_conv.h"CPL_CVSID("$Id: ogrdodslayer.cpp 10646 2007-01-18 02:38:10Z warmerdam $");/************************************************************************//*                            OGRDODSLayer()                            *//************************************************************************/OGRDODSLayer::OGRDODSLayer( OGRDODSDataSource *poDSIn,                             const char *pszTargetIn,                            AttrTable *poOGRLayerInfoIn ){    poDS = poDSIn;    poFeatureDefn = NULL;    pszQuery = NULL;    pszFIDColumn = NULL;    poSRS = NULL;    iNextShapeId = 0;    pszTarget = CPLStrdup( pszTargetIn );    papoFields = NULL;    bDataLoaded = FALSE;    poConnection = NULL;    poTargetVar = NULL;    poOGRLayerInfo = poOGRLayerInfoIn;    bKnowExtent = FALSE;    poDataDDS = new DataDDS( poDSIn->poBTF );/* ==================================================================== *//*      Harvest some metadata if available.                             *//* ==================================================================== */    if( poOGRLayerInfo != NULL )    {        string oMValue;/* -------------------------------------------------------------------- *//*      spatial_ref                                                     *//* -------------------------------------------------------------------- */        oMValue = poOGRLayerInfo->get_attr( "spatial_ref" );        if( oMValue.length() > 0 )        {            poSRS = new OGRSpatialReference();            if( poSRS->SetFromUserInput( oMValue.c_str() ) != OGRERR_NONE )            {                CPLError( CE_Warning, CPLE_AppDefined,                           "Ignoring unreconised SRS '%s'",                           oMValue.c_str() );                delete poSRS;                poSRS = NULL;            }        }/* -------------------------------------------------------------------- *//*      Layer extents.                                                  *//* -------------------------------------------------------------------- */        AttrTable *poLayerExt=poOGRLayerInfo->find_container("layer_extents");        if( poLayerExt != NULL )        {            bKnowExtent = TRUE;            sExtent.MinX = atof(poLayerExt->get_attr("x_min").c_str());            sExtent.MaxX = atof(poLayerExt->get_attr("x_max").c_str());            sExtent.MinY = atof(poLayerExt->get_attr("y_min").c_str());            sExtent.MaxY = atof(poLayerExt->get_attr("y_max").c_str());        }    }/* ==================================================================== *//*      Are we actually referencing a nested subsequence?  If so, we    *//*      need to find the super sequence so we can do the layered        *//*      stepping.                                                       *//* ==================================================================== */}/************************************************************************//*                           ~OGRDODSLayer()                            *//************************************************************************/OGRDODSLayer::~OGRDODSLayer(){    if( m_nFeaturesRead > 0 && poFeatureDefn != NULL )    {        CPLDebug( "DODS", "%d features read on layer '%s'.",                  (int) m_nFeaturesRead,                   poFeatureDefn->GetName() );    }    if( papoFields != NULL )    {        int iField;        for( iField = 0; iField < poFeatureDefn->GetFieldCount(); iField++ )            delete papoFields[iField];        CPLFree( papoFields );    }    if( poSRS != NULL )        poSRS->Release();    if( poFeatureDefn != NULL )        poFeatureDefn->Release();    CPLFree( pszFIDColumn );    pszFIDColumn = NULL;    CPLFree( pszTarget );    pszTarget = NULL;    if( poConnection != NULL )        delete poConnection;    delete poDataDDS;}/************************************************************************//*                            ResetReading()                            *//************************************************************************/void OGRDODSLayer::ResetReading(){    iNextShapeId = 0;}/************************************************************************//*                           GetNextFeature()                           *//************************************************************************/OGRFeature *OGRDODSLayer::GetNextFeature(){    OGRFeature *poFeature;    for( poFeature = GetFeature( iNextShapeId++ );          poFeature != NULL;         poFeature = GetFeature( iNextShapeId++ ) )    {        if( FilterGeometry( poFeature->GetGeometryRef() )            && (m_poAttrQuery == NULL                || m_poAttrQuery->Evaluate( poFeature )) )            return poFeature;        delete poFeature;    }    return NULL;}/************************************************************************//*                           TestCapability()                           *//************************************************************************/int OGRDODSLayer::TestCapability( const char * pszCap ){    return FALSE;}/************************************************************************//*                           GetSpatialRef()                            *//************************************************************************/OGRSpatialReference *OGRDODSLayer::GetSpatialRef(){    return poSRS;}/************************************************************************//*                           ProvideDataDDS()                           *//************************************************************************/int OGRDODSLayer::ProvideDataDDS(){    if( bDataLoaded )        return poTargetVar != NULL;    bDataLoaded = TRUE;    try    {        poConnection = new AISConnect( poDS->oBaseURL );        CPLDebug( "DODS", "request_data(%s,%s)",                  poDS->oBaseURL.c_str(),                  (poDS->oProjection + poDS->oConstraints).c_str() );        // We may need to use custom constraints here.         poConnection->request_data( *poDataDDS,                                     poDS->oProjection + poDS->oConstraints );    }    catch (Error &e)    {        CPLError( CE_Failure, CPLE_AppDefined,                  "DataDDS request failed:\n%s",                   e.get_error_message().c_str() );        return FALSE;    }    poTargetVar = poDataDDS->var( pszTarget );    return poTargetVar != NULL;}/************************************************************************//*                             GetExtent()                              *//************************************************************************/OGRErr OGRDODSLayer::GetExtent(OGREnvelope *psExtent, int bForce){    OGRErr eErr;    if( bKnowExtent )    {        *psExtent = this->sExtent;        return OGRERR_NONE;    }    if( !bForce )        return OGRERR_FAILURE;    eErr = OGRLayer::GetExtent( &sExtent, bForce );    if( eErr == OGRERR_NONE )    {        bKnowExtent = TRUE;        *psExtent = sExtent;    }    return eErr;}