/****************************************************************************** * $Id: ogrntfdatasource.cpp 10646 2007-01-18 02:38:10Z warmerdam $ * * Project:  UK NTF Reader * Purpose:  Implements OGRNTFDataSource class * Author:   Frank Warmerdam, warmerdam@pobox.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. ****************************************************************************/#include "ntf.h"#include "cpl_conv.h"#include "cpl_string.h"CPL_CVSID("$Id: ogrntfdatasource.cpp 10646 2007-01-18 02:38:10Z warmerdam $");/************************************************************************//*                          OGRNTFDataSource()                          *//************************************************************************/OGRNTFDataSource::OGRNTFDataSource(){    nLayers = 0;    papoLayers = NULL;    nNTFFileCount = 0;    papoNTFFileReader = NULL;    pszName = NULL;    iCurrentReader = -1;    iCurrentFC = 0;    nFCCount = 0;    papszFCNum = NULL;    papszFCName = NULL;    poFCLayer = NULL;    papszOptions = NULL;    poSpatialRef = new OGRSpatialReference( "PROJCS[\"OSGB 1936 / British National Grid\",GEOGCS[\"OSGB 1936\",DATUM[\"OSGB_1936\",SPHEROID[\"Airy 1830\",6377563.396,299.3249646,AUTHORITY[\"EPSG\",\"7001\"]],AUTHORITY[\"EPSG\",\"6277\"]],PRIMEM[\"Greenwich\",0,AUTHORITY[\"EPSG\",\"8901\"]],UNIT[\"degree\",0.0174532925199433],AUTHORITY[\"EPSG\",\"4277\"]],PROJECTION[\"Transverse_Mercator\"],PARAMETER[\"latitude_of_origin\",49],PARAMETER[\"central_meridian\",-2],PARAMETER[\"scale_factor\",0.999601272],PARAMETER[\"false_easting\",400000],PARAMETER[\"false_northing\",-100000],UNIT[\"metre\",1,AUTHORITY[\"EPSG\",\"9001\"]],AUTHORITY[\"EPSG\",\"27700\"]]" );/* -------------------------------------------------------------------- *//*      Allow initialization of options from the environment.           *//* -------------------------------------------------------------------- */    if( getenv("OGR_NTF_OPTIONS") != NULL )    {        papszOptions =             CSLTokenizeStringComplex( getenv("OGR_NTF_OPTIONS"), ",",                                      FALSE, FALSE );    }}/************************************************************************//*                         ~OGRNTFDataSource()                          *//************************************************************************/OGRNTFDataSource::~OGRNTFDataSource(){    int         i;    for( i = 0; i < nNTFFileCount; i++ )        delete papoNTFFileReader[i];    CPLFree( papoNTFFileReader );    for( i = 0; i < nLayers; i++ )        delete papoLayers[i];    if( poFCLayer != NULL )        delete poFCLayer;        CPLFree( papoLayers );    CPLFree( pszName );    CSLDestroy( papszOptions );    CSLDestroy( papszFCNum );    CSLDestroy( papszFCName );    if( poSpatialRef )        poSpatialRef->Release();}/************************************************************************//*                           TestCapability()                           *//************************************************************************/int OGRNTFDataSource::TestCapability( const char * ){    return FALSE;}/************************************************************************//*                           GetNamedLayer()                            *//************************************************************************/OGRNTFLayer * OGRNTFDataSource::GetNamedLayer( const char * pszName ){    for( int i = 0; i < nLayers; i++ )    {        if( EQUAL(papoLayers[i]->GetLayerDefn()->GetName(),pszName) )            return (OGRNTFLayer *) papoLayers[i];    }    return NULL;}/************************************************************************//*                              AddLayer()                              *//************************************************************************/void OGRNTFDataSource::AddLayer( OGRLayer * poNewLayer ){    papoLayers = (OGRLayer **)        CPLRealloc( papoLayers, sizeof(void*) * ++nLayers );        papoLayers[nLayers-1] = poNewLayer;}/************************************************************************//*                              GetLayer()                              *//************************************************************************/OGRLayer *OGRNTFDataSource::GetLayer( int iLayer ){    if( iLayer < 0 || iLayer > nLayers )        return NULL;    else if( iLayer == nLayers )        return poFCLayer;    else        return papoLayers[iLayer];}/************************************************************************//*                           GetLayerCount()                            *//************************************************************************/int OGRNTFDataSource::GetLayerCount(){    if( poFCLayer == NULL )        return nLayers;    else        return nLayers + 1;}/************************************************************************//*                                Open()                                *//************************************************************************/int OGRNTFDataSource::Open( const char * pszFilename, int bTestOpen,                            char ** papszLimitedFileList ){    VSIStatBuf      stat;    char            **papszFileList = NULL;    pszName = CPLStrdup( pszFilename );/* -------------------------------------------------------------------- *//*      Is the given path a directory or a regular file?                *//* -------------------------------------------------------------------- */    if( CPLStat( pszFilename, &stat ) != 0         || (!VSI_ISDIR(stat.st_mode) && !VSI_ISREG(stat.st_mode)) )    {        if( !bTestOpen )            CPLError( CE_Failure, CPLE_AppDefined,                   "%s is neither a file or directory, NTF access failed.\n",                      pszFilename );        return FALSE;    }    /* -------------------------------------------------------------------- *//*      Build a list of filenames we figure are NTF files.              *//* -------------------------------------------------------------------- */    if( VSI_ISREG(stat.st_mode) )    {        papszFileList = CSLAddString( NULL, pszFilename );    }    else    {        char      **candidateFileList = CPLReadDir( pszFilename );        int         i;        for( i = 0;              candidateFileList != NULL && candidateFileList[i] != NULL;              i++ )         {            if( papszLimitedFileList != NULL                && CSLFindString(papszLimitedFileList,                                 candidateFileList[i]) == -1 )            {                continue;            }                        if( strlen(candidateFileList[i]) > 4              && EQUALN(candidateFileList[i] + strlen(candidateFileList[i])-4,                       ".ntf",4) )            {                char       fullFilename[2048];                sprintf( fullFilename, "%s%c%s",                          pszFilename,#ifdef WIN32                         '\\',#else                         '/',#endif                         candidateFileList[i] );                papszFileList = CSLAddString( papszFileList, fullFilename );            }        }        CSLDestroy( candidateFileList );        if( CSLCount(papszFileList) == 0 )        {            if( !bTestOpen )                CPLError( CE_Failure, CPLE_OpenFailed,                          "No candidate NTF files (.ntf) found in\n"                          "directory: %s",                          pszFilename );            return FALSE;        }    }/* -------------------------------------------------------------------- *//*      Loop over all these files trying to open them.  In testopen     *//*      mode we first read the first 80 characters, to verify that      *//*      it looks like an NTF file.  Note that we don't keep the file    *//*      open ... we don't want to occupy alot of file handles when      *//*      handling a whole directory.                                     *//* -------------------------------------------------------------------- */    int         i;    papoNTFFileReader = (NTFFileReader **)        CPLCalloc(sizeof(void*), CSLCount(papszFileList));        for( i = 0; papszFileList[i] != NULL; i++ )    {        if( bTestOpen )        {            char        szHeader[80];            FILE        *fp;            int         j;