EVOLUTION-MANAGER

Edit File: georaster_dataset.cpp

/****************************************************************************** * * Name: georaster_dataset.cpp * Project: Oracle Spatial GeoRaster Driver * Purpose: Implement GeoRasterDataset Methods * Author: Ivan Lucena [ivan.lucena at oracle.com] * ****************************************************************************** * Copyright (c) 2008, Ivan Lucena <ivan dot lucena at oracle dot com> * Copyright (c) 2013, Even Rouault <even dot rouault at mines-paris dot org> * * 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 "cpl_error.h" #include "cpl_vsi_virtual.h" #include "gdaljp2metadata.h" #include "cpl_list.h" #include "gdal.h" #include "gdal_frmts.h" #include "gdal_priv.h" #include "ogr_spatialref.h" #include "georaster_priv.h" CPL_CVSID("$Id: georaster_dataset.cpp 40616 2017-11-02 12:40:51Z ilucena $") // --------------------------------------------------------------------------- // GeoRasterDataset() // --------------------------------------------------------------------------- GeoRasterDataset::GeoRasterDataset() { bGeoTransform = false; bForcedSRID = false; poGeoRaster = NULL; papszSubdatasets = NULL; adfGeoTransform[0] = 0.0; adfGeoTransform[1] = 1.0; adfGeoTransform[2] = 0.0; adfGeoTransform[3] = 0.0; adfGeoTransform[4] = 0.0; adfGeoTransform[5] = 1.0; pszProjection = NULL; poMaskBand = NULL; bApplyNoDataArray = false; poJP2Dataset = NULL; } // --------------------------------------------------------------------------- // ~GeoRasterDataset() // --------------------------------------------------------------------------- GeoRasterDataset::~GeoRasterDataset() { FlushCache(); poGeoRaster->FlushMetadata(); delete poGeoRaster; if( poMaskBand ) { delete poMaskBand; } if( poJP2Dataset ) { delete poJP2Dataset; } CPLFree( pszProjection ); CSLDestroy( papszSubdatasets ); } // --------------------------------------------------------------------------- // Identify() // --------------------------------------------------------------------------- int GeoRasterDataset::Identify( GDALOpenInfo* poOpenInfo ) { // ------------------------------------------------------------------- // Verify georaster prefix // ------------------------------------------------------------------- char* pszFilename = poOpenInfo->pszFilename; if( STARTS_WITH_CI(pszFilename, "georaster:") == false && STARTS_WITH_CI(pszFilename, "geor:") == false ) { return false; } return true; } // --------------------------------------------------------------------------- // Open() // --------------------------------------------------------------------------- GDALDataset* GeoRasterDataset::Open( GDALOpenInfo* poOpenInfo ) { // ------------------------------------------------------------------- // It should not have an open file pointer. // ------------------------------------------------------------------- if( poOpenInfo->fpL != NULL ) { return NULL; } // ------------------------------------------------------------------- // Check identification string and usage // ------------------------------------------------------------------- if( ! Identify( poOpenInfo ) ) { return NULL; } // ------------------------------------------------------------------- // Create a GeoRaster wrapper object // ------------------------------------------------------------------- GeoRasterWrapper* poGRW = GeoRasterWrapper::Open( poOpenInfo->pszFilename, poOpenInfo->eAccess == GA_Update ); if( ! poGRW ) { return NULL; } // ------------------------------------------------------------------- // Create a corresponding GDALDataset // ------------------------------------------------------------------- GeoRasterDataset *poGRD = new GeoRasterDataset(); if( ! poGRD ) { return NULL; } poGRD->eAccess = poOpenInfo->eAccess; poGRD->poGeoRaster = poGRW; // ------------------------------------------------------------------- // List Subdatasets // ------------------------------------------------------------------- if( ! poGRW->bUniqueFound ) { if( poGRD->eAccess == GA_ReadOnly ) { poGRD->SetSubdatasets( poGRW ); if( CSLCount( poGRD->papszSubdatasets ) == 0 ) { delete poGRD; poGRD = NULL; } } return (GDALDataset*) poGRD; } // ------------------------------------------------------------------- // Assign GeoRaster information // ------------------------------------------------------------------- poGRD->poGeoRaster = poGRW; poGRD->nRasterXSize = poGRW->nRasterColumns; poGRD->nRasterYSize = poGRW->nRasterRows; poGRD->nBands = poGRW->nRasterBands; if( poGRW->bIsReferenced ) { poGRD->adfGeoTransform[1] = poGRW->dfXCoefficient[0]; poGRD->adfGeoTransform[2] = poGRW->dfXCoefficient[1]; poGRD->adfGeoTransform[0] = poGRW->dfXCoefficient[2]; poGRD->adfGeoTransform[4] = poGRW->dfYCoefficient[0]; poGRD->adfGeoTransform[5] = poGRW->dfYCoefficient[1]; poGRD->adfGeoTransform[3] = poGRW->dfYCoefficient[2]; } // ------------------------------------------------------------------- // Copy RPC values to RPC metadata domain // ------------------------------------------------------------------- if( poGRW->phRPC ) { char **papszRPC_MD = RPCInfoToMD( poGRW->phRPC ); char **papszSanitazed = NULL; int i = 0; int n = CSLCount( papszRPC_MD ); for( i = 0; i < n; i++ ) { if ( STARTS_WITH_CI(papszRPC_MD[i], "MIN_LAT") || STARTS_WITH_CI(papszRPC_MD[i], "MIN_LONG") || STARTS_WITH_CI(papszRPC_MD[i], "MAX_LAT") || STARTS_WITH_CI(papszRPC_MD[i], "MAX_LONG") ) { continue; } papszSanitazed = CSLAddString( papszSanitazed, papszRPC_MD[i] ); } poGRD->SetMetadata( papszSanitazed, "RPC" ); CSLDestroy( papszRPC_MD ); CSLDestroy( papszSanitazed ); } // ------------------------------------------------------------------- // Open for JPEG 2000 compression for reading // ------------------------------------------------------------------- if( EQUAL( poGRW->sCompressionType.c_str(), "JP2-F" ) && poGRD->eAccess == GA_ReadOnly ) { poGRD->JP2_Open( poOpenInfo->eAccess ); if( ! poGRD->poJP2Dataset ) { delete poGRD; return NULL; } } // ------------------------------------------------------------------- // Load mask band // ------------------------------------------------------------------- poGRW->bHasBitmapMask = EQUAL( "TRUE", CPLGetXMLValue( poGRW->phMetadata, "layerInfo.objectLayer.bitmapMask", "FALSE" ) ); if( poGRW->bHasBitmapMask ) { poGRD->poMaskBand = new GeoRasterRasterBand( poGRD, 0, DEFAULT_BMP_MASK ); } // ------------------------------------------------------------------- // Check for filter Nodata environment variable, default is YES // ------------------------------------------------------------------- const char *pszGEOR_FILTER_NODATA = CPLGetConfigOption( "GEOR_FILTER_NODATA_VALUES", "NO" ); if( ! EQUAL(pszGEOR_FILTER_NODATA, "NO") ) { poGRD->bApplyNoDataArray = true; } // ------------------------------------------------------------------- // Create bands // ------------------------------------------------------------------- int i = 0; for( i = 1; i <= poGRD->nBands; i++ ) { poGRD->SetBand( i, new GeoRasterRasterBand( poGRD, i, 0 , poGRD->poJP2Dataset) ); } // ------------------------------------------------------------------- // Set IMAGE_STRUCTURE metadata information // ------------------------------------------------------------------- if( poGRW->nBandBlockSize == 1 ) { poGRD->SetMetadataItem( "INTERLEAVE", "BSQ", "IMAGE_STRUCTURE" ); } else { if( EQUAL( poGRW->sInterleaving.c_str(), "BSQ" ) ) { poGRD->SetMetadataItem( "INTERLEAVE", "BSQ", "IMAGE_STRUCTURE" ); } else if( EQUAL( poGRW->sInterleaving.c_str(), "BIP" ) ) { poGRD->SetMetadataItem( "INTERLEAVE", "PIB", "IMAGE_STRUCTURE" ); } else if( EQUAL( poGRW->sInterleaving.c_str(), "BIL" ) ) { poGRD->SetMetadataItem( "INTERLEAVE", "BIL", "IMAGE_STRUCTURE" ); } } poGRD->SetMetadataItem( "COMPRESSION", CPLGetXMLValue( poGRW->phMetadata, "rasterInfo.compression.type", "NONE" ), "IMAGE_STRUCTURE" ); if( STARTS_WITH_CI(poGRW->sCompressionType.c_str(), "JPEG") ) { poGRD->SetMetadataItem( "COMPRESSION_QUALITY", CPLGetXMLValue( poGRW->phMetadata, "rasterInfo.compression.quality", "undefined" ), "IMAGE_STRUCTURE" ); } if( EQUAL( poGRW->sCellDepth.c_str(), "1BIT" ) ) { poGRD->SetMetadataItem( "NBITS", "1", "IMAGE_STRUCTURE" ); } if( EQUAL( poGRW->sCellDepth.c_str(), "2BIT" ) ) { poGRD->SetMetadataItem( "NBITS", "2", "IMAGE_STRUCTURE" ); } if( EQUAL( poGRW->sCellDepth.c_str(), "4BIT" ) ) { poGRD->SetMetadataItem( "NBITS", "4", "IMAGE_STRUCTURE" ); } // ------------------------------------------------------------------- // Set Metadata on "ORACLE" domain // ------------------------------------------------------------------- char* pszDoc = CPLSerializeXMLTree( poGRW->phMetadata ); poGRD->SetMetadataItem( "TABLE_NAME", CPLSPrintf( "%s%s", poGRW->sSchema.c_str(), poGRW->sTable.c_str()), "ORACLE" ); poGRD->SetMetadataItem( "COLUMN_NAME", poGRW->sColumn.c_str(), "ORACLE" ); poGRD->SetMetadataItem( "RDT_TABLE_NAME", poGRW->sDataTable.c_str(), "ORACLE" ); poGRD->SetMetadataItem( "RASTER_ID", CPLSPrintf( "%lld", poGRW->nRasterId ), "ORACLE" ); poGRD->SetMetadataItem( "SRID", CPLSPrintf( "%lld", poGRW->nSRID ), "ORACLE" ); poGRD->SetMetadataItem( "WKT", poGRW->sWKText.c_str(), "ORACLE" ); poGRD->SetMetadataItem( "COMPRESSION", poGRW->sCompressionType.c_str(), "ORACLE" ); poGRD->SetMetadataItem( "METADATA", pszDoc, "ORACLE" ); CPLFree( pszDoc ); // ------------------------------------------------------------------- // Return a GDALDataset // ------------------------------------------------------------------- return (GDALDataset*) poGRD; } // --------------------------------------------------------------------------- // JP2Open() // --------------------------------------------------------------------------- void GeoRasterDataset::JP2_Open( GDALAccess /* eAccess */ ) { GDALDriver* poJP2Driver = NULL; static const char * const apszDrivers[] = { "JP2OPENJPEG", "JP2ECW", "JP2MRSID", "JPEG2000", "JP2KAK", NULL }; // Find at least one available JP2 driver for( int iDriver = 0; apszDrivers[iDriver] != NULL; iDriver++ ) { poJP2Driver = (GDALDriver*) GDALGetDriverByName(apszDrivers[iDriver]); if( poJP2Driver ) { break; } } // If JP2 driver is installed, try to open the LOB via VSIOCILOB handler poJP2Dataset = NULL; if( poJP2Driver ) { CPLString osDSName; osDSName.Printf( "/vsiocilob/%s,%s,%s,%s,%lld,noext", poGeoRaster->poConnection->GetUser(), poGeoRaster->poConnection->GetPassword(), poGeoRaster->poConnection->GetServer(), poGeoRaster->sDataTable.c_str(), poGeoRaster->nRasterId ); CPLPushErrorHandler( CPLQuietErrorHandler ); poJP2Dataset = (GDALDataset*) GDALOpenEx( osDSName.c_str(), GDAL_OF_RASTER, apszDrivers, NULL, NULL ); CPLPopErrorHandler(); if( ! poJP2Dataset ) { CPLString osLastErrorMsg(CPLGetLastErrorMsg()); CPLError( CE_Failure, CPLE_AppDefined, "Unable to open JPEG2000 image within GeoRaster dataset.\n%s", osLastErrorMsg.c_str() ); } } else { CPLError( CE_Failure, CPLE_AppDefined, "Unable to open JPEG2000 image within GeoRaster dataset.\n%s", "No JPEG2000 capable driver (JP2OPENJPEG, " "JP2ECW, JP2MRSID, etc...) is available." ); } } // --------------------------------------------------------------------------- // JP2CreateCopy() // --------------------------------------------------------------------------- void GeoRasterDataset::JP2_CreateCopy( GDALDataset* poJP2DS, char** papszOptions, int* pnResolutions, GDALProgressFunc pfnProgress, void* pProgressData ) { GDALDriver* poJP2Driver = NULL; static const char * const apszDrivers[] = { "JP2OPENJPEG", "JP2ECW", "JP2MRSID", "JPEG2000", "JP2KAK", NULL }; // Find at least one available JP2 driver for( int iDriver = 0; apszDrivers[iDriver] != NULL; iDriver++ ) { poJP2Driver = (GDALDriver*) GDALGetDriverByName(apszDrivers[iDriver]); if( poJP2Driver ) { break; } } // If a JP2 driver is installed calls driver's CreateCopy poJP2Dataset = NULL; if( poJP2Driver ) { char** papszOpt = NULL; const char* pszFetched = CSLFetchNameValue( papszOptions, "JP2_BLOCKXSIZE" ); if( pszFetched ) { papszOpt = CSLAddNameValue( papszOpt, "BLOCKXSIZE", pszFetched ); papszOpt = CSLAddNameValue( papszOpt, "TILE_HEIGHT", pszFetched ); } CPLDebug("GEOR","JP2_BLOCKXSIZE %s", pszFetched ); pszFetched = CSLFetchNameValue( papszOptions, "JP2_BLOCKYSIZE" ); if( pszFetched ) { papszOpt = CSLAddNameValue( papszOpt, "BLOCKYSIZE", pszFetched ); papszOpt = CSLAddNameValue( papszOpt, "TILE_WIDTH", pszFetched ); } pszFetched = CSLFetchNameValue( papszOptions, "JP2_QUALITY" ); if( pszFetched ) { papszOpt = CSLAddNameValue( papszOpt, "QUALITY", pszFetched ); if( STARTS_WITH_CI( pszFetched, "100" ) ) { papszOpt = CSLAddNameValue( papszOpt, "REVERSIBLE", "TRUE" ); } poGeoRaster->nCompressQuality = atoi( pszFetched ); } else { poGeoRaster->nCompressQuality = 25; // JP2OpenJPEG default... } pszFetched = CSLFetchNameValue( papszOptions, "JP2_REVERSIBLE" ); if( pszFetched ) { papszOpt = CSLAddNameValue( papszOpt, "REVERSIBLE", pszFetched ); } pszFetched = CSLFetchNameValue( papszOptions, "JP2_RESOLUTIONS" ); if( pszFetched ) { papszOpt = CSLAddNameValue( papszOpt, "RESOLUTIONS", pszFetched ); papszOpt = CSLAddNameValue( papszOpt, "RESOLUTIONS_LEVELS", pszFetched ); papszOpt = CSLAddNameValue( papszOpt, "LAYERS", pszFetched ); } pszFetched = CSLFetchNameValue( papszOptions, "JP2_PROGRESSION" ); if( pszFetched ) { papszOpt = CSLAddNameValue( papszOpt, "PROGRESSION", pszFetched ); } papszOpt = CSLAddNameValue( papszOpt, "CODEC", "JP2" ); papszOpt = CSLAddNameValue( papszOpt, "GeoJP2", "NO" ); papszOpt = CSLAddNameValue( papszOpt, "GMLJP2", "NO" ); papszOpt = CSLAddNameValue( papszOpt, "YCBCR420", "NO" ); papszOpt = CSLAddNameValue( papszOpt, "TARGET", "0" ); CPLPushErrorHandler( CPLQuietErrorHandler ); CPLString osDSName; osDSName.Printf( "/vsiocilob/%s,%s,%s,%s,%lld,noext", poGeoRaster->poConnection->GetUser(), poGeoRaster->poConnection->GetPassword(), poGeoRaster->poConnection->GetServer(), poGeoRaster->sDataTable.c_str(), poGeoRaster->nRasterId ); poJP2Dataset = (GDALDataset*) GDALCreateCopy( poJP2Driver, osDSName.c_str(), poJP2DS, false, (char**) papszOpt, pfnProgress, pProgressData ); CPLPopErrorHandler(); CSLDestroy( papszOpt ); if( ! poJP2Dataset ) { CPLString osLastErrorMsg(CPLGetLastErrorMsg()); CPLError( CE_Failure, CPLE_AppDefined, "Unable to copy JPEG2000 image within GeoRaster dataset.\n%s", osLastErrorMsg.c_str() ); return; } } else { CPLError( CE_Failure, CPLE_AppDefined, "Unable to copy JPEG2000 image within GeoRaster dataset.\n%s", "No JPEG2000 capable driver (JP2OPENJPEG, " "JP2ECW, JP2MRSID, etc...) is available." ); return; } // Retrieve the number of resolutions based on the number of overviews CPLPushErrorHandler( CPLQuietErrorHandler ); *pnResolutions = poJP2Dataset->GetRasterBand(1)->GetOverviewCount() + 1; delete poJP2Dataset; CPLPopErrorHandler(); // Avoid showing warning regards writing aux.xml file poJP2Dataset = NULL; } // --------------------------------------------------------------------------- // JP2_CopyDirect() // --------------------------------------------------------------------------- boolean GeoRasterDataset::JP2_CopyDirect( const char* pszJP2Filename, int* pnResolutions, GDALProgressFunc pfnProgress, void* pProgressData ) { char** papszFileList = GetFileList(); if( CSLCount(papszFileList) == 0 ) { return false; } VSILFILE *fpInput = VSIFOpenL( pszJP2Filename, "r" ); VSILFILE *fpOutput = VSIFOpenL( papszFileList[0], "wb" ); size_t nCache = (size_t) ( GDALGetCacheMax() * 0.25 ); void *pBuffer = (GByte*) VSIMalloc( sizeof(GByte) * nCache ); GDALJP2Box oBox( fpInput ); (void) oBox.ReadFirst(); GUInt32 nLBox; GUInt32 nTBox; int nBoxCount = 0; while( strlen(oBox.GetType()) > 0 ) { nBoxCount++; if( EQUAL( oBox.GetType(), "jp " ) || EQUAL( oBox.GetType(), "ftyp" ) || EQUAL( oBox.GetType(), "jp2h" ) ) { size_t nDataLength = (size_t) oBox.GetDataLength(); size_t nSize = VSIFReadL( pBuffer, 1, nDataLength, fpInput); if ( nSize != nDataLength ) { CPLError( CE_Warning, CPLE_AppDefined, "amount read differs from JP2 Box data length" ); } nLBox = CPL_MSBWORD32( (int) nDataLength + 8 ); memcpy( &nTBox, oBox.GetType(), 4 ); VSIFWriteL( &nLBox, 4, 1, fpOutput ); VSIFWriteL( &nTBox, 4, 1, fpOutput ); VSIFWriteL( pBuffer, 1, nSize, fpOutput ); } if( EQUAL( oBox.GetType(), "jp2c" ) ) { size_t nCount = 0; size_t nSize = 0; size_t nDataLength = oBox.GetDataLength(); nLBox = CPL_MSBWORD32( (int) nDataLength + 8 ); memcpy( &nTBox, oBox.GetType(), 4 ); VSIFWriteL( &nLBox, 4, 1, fpOutput ); VSIFWriteL( &nTBox, 4, 1, fpOutput ); while( nCount < nDataLength ) { size_t nChunk = (size_t) MIN( nCache, nDataLength - nCount ); nSize = VSIFReadL( pBuffer, 1, nChunk, fpInput ); if ( nSize != nChunk ) { CPLError( CE_Warning, CPLE_AppDefined, "amount read differs from JP2 data length" ); } VSIFWriteL( pBuffer, 1, nSize, fpOutput ); nCount += nSize; pfnProgress( (float) nCount / (float) nDataLength, NULL, pProgressData ); } } if( ! oBox.ReadNext() ) { break; } } VSIFCloseL( fpInput ); VSIFCloseL( fpOutput ); CSLDestroy( papszFileList ); CPLFree( pBuffer ); // Retrieve the number of resolutions based on the number of overviews JP2_Open( GA_ReadOnly ); if( poJP2Dataset ) { *pnResolutions = poJP2Dataset->GetRasterBand(1)->GetOverviewCount() + 1; delete poJP2Dataset; poJP2Dataset = NULL; } return (nBoxCount > 0); } // --------------------------------------------------------------------------- // JPG_CopyDirect() // --------------------------------------------------------------------------- boolean GeoRasterDataset::JPEG_CopyDirect( const char* pszJPGFilename, GDALProgressFunc pfnProgress, void* pProgressData ) { OWConnection* poConnection = poGeoRaster->poConnection; OWStatement* poStmt = NULL; OCILobLocator* poLocator; poStmt = poConnection->CreateStatement( CPLSPrintf( "select rasterblock from %s where rasterid = %lld " "and rownum = 1 for update", poGeoRaster->sDataTable.c_str(), poGeoRaster->nRasterId ) ); poStmt->Define( &poLocator ); if( poStmt->Execute() ) { VSILFILE *fpInput = VSIFOpenL( pszJPGFilename, "r" ); size_t nCache = (size_t) ( GDALGetCacheMax() * 0.25 ); void *pBuffer = (GByte*) VSIMalloc( sizeof(GByte) * nCache ); VSIFSeekL( fpInput, 0L, SEEK_END); size_t nCount = 0; size_t nSize = 0; size_t nDataLength = VSIFTellL( fpInput ); VSIFSeekL( fpInput, 0L, SEEK_SET ); GUIntBig nWrite = (GUIntBig) 0; GUIntBig nCurOff = (GUIntBig) 0; while( nCount < nDataLength ) { size_t nChunk = (size_t) MIN( nCache, nDataLength - nCount ); nSize = VSIFReadL( pBuffer, 1, nChunk, fpInput ); if ( nSize != nChunk ) { CPLError( CE_Warning, CPLE_AppDefined, "amount read differs from JPG length" ); } nWrite = poStmt->WriteBlob( poLocator, (void*) pBuffer, (nCurOff + 1), nSize ); nCurOff += nWrite; nCount += nSize; pfnProgress( (float) nCount / (float) nDataLength, NULL, pProgressData ); } VSIFCloseL( fpInput ); CPLFree( pBuffer ); return true; } if( poLocator ) { OWStatement::Free( &poLocator, 1 ); } if( poStmt ) { delete poStmt; } return false; } // --------------------------------------------------------------------------- // GetFileList() // --------------------------------------------------------------------------- char** GeoRasterDataset::GetFileList() { char** papszFileList = NULL; if( EQUAL( poGeoRaster->sCompressionType.c_str(), "JP2-F" ) ) { CPLString osDSName; osDSName.Printf( "/vsiocilob/%s,%s,%s,%s,%lld,noext", this->poGeoRaster->poConnection->GetUser(), this->poGeoRaster->poConnection->GetPassword(), this->poGeoRaster->poConnection->GetServer(), this->poGeoRaster->sDataTable.c_str(), this->poGeoRaster->nRasterId ); papszFileList = CSLAddString( papszFileList, osDSName.c_str() ); } return papszFileList; } // --------------------------------------------------------------------------- // Create() // --------------------------------------------------------------------------- GDALDataset *GeoRasterDataset::Create( const char *pszFilename, int nXSize, int nYSize, int nBands, GDALDataType eType, char **papszOptions ) { // ------------------------------------------------------------------- // Check for supported Data types // ------------------------------------------------------------------- const char* pszCellDepth = OWSetDataType( eType ); if( EQUAL( pszCellDepth, "Unknown" ) ) { CPLError( CE_Failure, CPLE_AppDefined, "Attempt to create GeoRaster with unsupported data type (%s)", GDALGetDataTypeName( eType ) ); return NULL; } // ------------------------------------------------------------------- // Open the Dataset // ------------------------------------------------------------------- GeoRasterDataset* poGRD = NULL; poGRD = (GeoRasterDataset*) GDALOpen( pszFilename, GA_Update ); if( ! poGRD ) { return NULL; } // ------------------------------------------------------------------- // Get the GeoRaster // ------------------------------------------------------------------- GeoRasterWrapper* poGRW = poGRD->poGeoRaster; if( ! poGRW ) { delete poGRD; return NULL; } // ------------------------------------------------------------------- // Set basic information and default values // ------------------------------------------------------------------- poGRW->nRasterColumns = nXSize; poGRW->nRasterRows = nYSize; poGRW->nRasterBands = nBands; poGRW->sCellDepth = pszCellDepth; poGRW->nRowBlockSize = DEFAULT_BLOCK_ROWS; poGRW->nColumnBlockSize = DEFAULT_BLOCK_COLUMNS; poGRW->nBandBlockSize = 1; if( poGRW->bUniqueFound ) { poGRW->PrepareToOverwrite(); } // ------------------------------------------------------------------- // Check the create options to use in initialization // ------------------------------------------------------------------- const char* pszFetched = ""; char* pszDescription = NULL; char* pszInsert = NULL; int nQuality = -1; if( ! poGRW->sTable.empty() ) { pszFetched = CSLFetchNameValue( papszOptions, "DESCRIPTION" ); if( pszFetched ) { pszDescription = CPLStrdup( pszFetched ); } } if( poGRW->sTable.empty() ) { poGRW->sTable = "GDAL_IMPORT"; poGRW->sDataTable = "GDAL_RDT"; } if( poGRW->sColumn.empty() ) { poGRW->sColumn = "RASTER"; } pszFetched = CSLFetchNameValue( papszOptions, "INSERT" ); if( pszFetched ) { pszInsert = CPLStrdup( pszFetched ); } pszFetched = CSLFetchNameValue( papszOptions, "BLOCKXSIZE" ); if( pszFetched ) { poGRW->nColumnBlockSize = atoi( pszFetched ); } pszFetched = CSLFetchNameValue( papszOptions, "BLOCKYSIZE" ); if( pszFetched ) { poGRW->nRowBlockSize = atoi( pszFetched ); } pszFetched = CSLFetchNameValue( papszOptions, "NBITS" ); if( pszFetched != NULL ) { poGRW->sCellDepth = CPLSPrintf( "%dBIT", atoi( pszFetched ) ); } pszFetched = CSLFetchNameValue( papszOptions, "COMPRESS" ); if( pszFetched != NULL && ( EQUAL( pszFetched, "JPEG-F" ) || EQUAL( pszFetched, "JP2-F" ) || EQUAL( pszFetched, "DEFLATE" ) ) ) { poGRW->sCompressionType = pszFetched; } else { poGRW->sCompressionType = "NONE"; } pszFetched = CSLFetchNameValue( papszOptions, "QUALITY" ); if( pszFetched ) { poGRW->nCompressQuality = atoi( pszFetched ); nQuality = poGRW->nCompressQuality; } pszFetched = CSLFetchNameValue( papszOptions, "INTERLEAVE" ); bool bInterleve_ind = false; if( pszFetched ) { bInterleve_ind = true; if( EQUAL( pszFetched, "BAND" ) || EQUAL( pszFetched, "BSQ" ) ) { poGRW->sInterleaving = "BSQ"; } if( EQUAL( pszFetched, "LINE" ) || EQUAL( pszFetched, "BIL" ) ) { poGRW->sInterleaving = "BIL"; } if( EQUAL( pszFetched, "PIXEL" ) || EQUAL( pszFetched, "BIP" ) ) { poGRW->sInterleaving = "BIP"; } } else { if( EQUAL( poGRW->sCompressionType.c_str(), "NONE" ) == false ) { poGRW->sInterleaving = "BIP"; } } pszFetched = CSLFetchNameValue( papszOptions, "BLOCKBSIZE" ); if( pszFetched ) { poGRW->nBandBlockSize = atoi( pszFetched ); } else { if( nBands == 3 || nBands == 4 ) { poGRW->nBandBlockSize = nBands; } } if( bInterleve_ind == false && ( poGRW->nBandBlockSize == 3 || poGRW->nBandBlockSize == 4 ) ) { poGRW->sInterleaving = "BIP"; } if( STARTS_WITH_CI(poGRW->sCompressionType.c_str(), "JPEG") ) { if( ! EQUAL( poGRW->sInterleaving.c_str(), "BIP" ) ) { CPLError( CE_Warning, CPLE_IllegalArg, "compress=JPEG assumes interleave=BIP" ); poGRW->sInterleaving = "BIP"; } } pszFetched = CSLFetchNameValue( papszOptions, "BLOCKING" ); if( pszFetched ) { if( EQUAL( pszFetched, "NO" ) ) { poGRW->bBlocking = false; } if( EQUAL( pszFetched, "OPTIMALPADDING" ) ) { if( poGRW->poConnection->GetVersion() < 11 ) { CPLError( CE_Warning, CPLE_IllegalArg, "BLOCKING=OPTIMALPADDING not supported on Oracle older than 11g" ); } else { poGRW->bAutoBlocking = true; poGRW->bBlocking = true; } } } // ------------------------------------------------------------------- // Validate options // ------------------------------------------------------------------- if( pszDescription && poGRW->bUniqueFound ) { CPLError( CE_Failure, CPLE_IllegalArg, "Option (DESCRIPTION) cannot be used on a existing GeoRaster." ); delete poGRD; return NULL; } if( pszInsert && poGRW->bUniqueFound ) { CPLError( CE_Failure, CPLE_IllegalArg, "Option (INSERT) cannot be used on a existing GeoRaster." ); delete poGRD; return NULL; } /* Compression JPEG-B is deprecated. It should be able to read but not * to create new GeoRaster on databases with that compression option. * * TODO: Remove that options on next release. */ if( EQUAL( poGRW->sCompressionType.c_str(), "JPEG-B" ) ) { CPLError( CE_Failure, CPLE_IllegalArg, "Option (COMPRESS=%s) is deprecated and cannot be used.", poGRW->sCompressionType.c_str() ); delete poGRD; return NULL; } if( EQUAL( poGRW->sCompressionType.c_str(), "JPEG-F" ) ) { /* JPEG-F can only compress byte data type */ if( eType != GDT_Byte ) { CPLError( CE_Failure, CPLE_IllegalArg, "Option (COMPRESS=%s) can only be used with Byte data type.", poGRW->sCompressionType.c_str() ); delete poGRD; return NULL; } /* JPEG-F can compress one band per block or 3 for RGB * or 4 for RGBA. */ if( ( poGRW->nBandBlockSize != 1 && poGRW->nBandBlockSize != 3 && poGRW->nBandBlockSize != 4 ) || ( ( poGRW->nBandBlockSize != 1 && ( poGRW->nBandBlockSize != poGRW->nRasterBands ) ) ) ) { CPLError( CE_Failure, CPLE_IllegalArg, "Option (COMPRESS=%s) requires BLOCKBSIZE to be 1 (for any " "number of bands), 3 (for 3 bands RGB) and 4 (for 4 bands RGBA).", poGRW->sCompressionType.c_str() ); delete poGRD; return NULL; } // There is a limit on how big a compressed block can be. if( ( poGRW->nColumnBlockSize * poGRW->nRowBlockSize * poGRW->nBandBlockSize * ( GDALGetDataTypeSize( eType ) / 8 ) ) > ( 50 * 1024 * 1024 ) ) { CPLError( CE_Failure, CPLE_IllegalArg, "Option (COMPRESS=%s) each data block must not exceed 50Mb. " "Consider reducing BLOCK{X,Y,B}XSIZE.", poGRW->sCompressionType.c_str() ); delete poGRD; return NULL; } } if( EQUAL( poGRW->sCompressionType.c_str(), "DEFLATE" ) ) { if( ( poGRW->nColumnBlockSize * poGRW->nRowBlockSize * poGRW->nBandBlockSize * ( GDALGetDataTypeSize( eType ) / 8 ) ) > ( 1024 * 1024 * 1024 ) ) { CPLError( CE_Failure, CPLE_IllegalArg, "For (COMPRESS=%s) each data block must not exceed 1Gb. " "Consider reducing BLOCK{X,Y,B}XSIZE.", poGRW->sCompressionType.c_str() ); delete poGRD; return NULL; } } // When the compression is JP2-F it should be just one block if( EQUAL( poGRW->sCompressionType.c_str(), "JP2-F" ) ) { poGRW->nRowBlockSize = poGRW->nRasterRows; poGRW->nColumnBlockSize = poGRW->nRasterColumns; poGRW->nBandBlockSize = poGRW->nRasterBands; poGRW->bBlocking = false; } pszFetched = CSLFetchNameValue( papszOptions, "OBJECTTABLE" ); if( pszFetched ) { int nVersion = poGRW->poConnection->GetVersion(); if( nVersion <= 11 ) { CPLError( CE_Failure, CPLE_IllegalArg, "Driver create-option OBJECTTABLE not " "supported on Oracle %d", nVersion ); delete poGRD; return NULL; } } poGRD->poGeoRaster->bCreateObjectTable = CPLFetchBool( papszOptions, "OBJECTTABLE", false ); // ------------------------------------------------------------------- // Create a SDO_GEORASTER object on the server // ------------------------------------------------------------------- const bool bSuccess = poGRW->Create( pszDescription, pszInsert, poGRW->bUniqueFound ); CPLFree( pszInsert ); CPLFree( pszDescription ); if( ! bSuccess ) { delete poGRD; return NULL; } // ------------------------------------------------------------------- // Prepare an identification string // ------------------------------------------------------------------- char szStringId[OWTEXT]; strcpy( szStringId, CPLSPrintf( "georaster:%s,%s,%s,%s,%lld", poGRW->poConnection->GetUser(), poGRW->poConnection->GetPassword(), poGRW->poConnection->GetServer(), poGRW->sDataTable.c_str(), poGRW->nRasterId ) ); delete poGRD; poGRD = (GeoRasterDataset*) GDALOpen( szStringId, GA_Update ); if( ! poGRD ) { return NULL; } // ------------------------------------------------------------------- // Load additional options // ------------------------------------------------------------------- pszFetched = CSLFetchNameValue( papszOptions, "VATNAME" ); if( pszFetched ) { poGRW->sValueAttributeTab = pszFetched; } pszFetched = CSLFetchNameValue( papszOptions, "SRID" ); if( pszFetched ) { poGRD->bForcedSRID = true; poGRD->poGeoRaster->SetGeoReference( atoi( pszFetched ) ); } poGRD->poGeoRaster->bGenSpatialExtent = CPLFetchBool( papszOptions, "SPATIALEXTENT", TRUE ); pszFetched = CSLFetchNameValue( papszOptions, "EXTENTSRID" ); if( pszFetched ) { poGRD->poGeoRaster->nExtentSRID = atoi( pszFetched ); } pszFetched = CSLFetchNameValue( papszOptions, "COORDLOCATION" ); if( pszFetched ) { if( EQUAL( pszFetched, "CENTER" ) ) { poGRD->poGeoRaster->eModelCoordLocation = MCL_CENTER; } else if( EQUAL( pszFetched, "UPPERLEFT" ) ) { poGRD->poGeoRaster->eModelCoordLocation = MCL_UPPERLEFT; } else { CPLError( CE_Warning, CPLE_IllegalArg, "Incorrect COORDLOCATION (%s)", pszFetched ); } } if ( nQuality > 0 ) { poGRD->poGeoRaster->nCompressQuality = nQuality; } pszFetched = CSLFetchNameValue( papszOptions, "GENPYRAMID" ); if( pszFetched != NULL ) { if (!(EQUAL(pszFetched, "NN") || EQUAL(pszFetched, "BILINEAR") || EQUAL(pszFetched, "BIQUADRATIC") || EQUAL(pszFetched, "CUBIC") || EQUAL(pszFetched, "AVERAGE4") || EQUAL(pszFetched, "AVERAGE16"))) { CPLError( CE_Warning, CPLE_IllegalArg, "Wrong resample method for pyramid (%s)", pszFetched); } poGRD->poGeoRaster->bGenPyramid = true; poGRD->poGeoRaster->sPyramidResampling = pszFetched; } pszFetched = CSLFetchNameValue( papszOptions, "GENPYRLEVELS" ); if( pszFetched != NULL ) { poGRD->poGeoRaster->bGenPyramid = true; poGRD->poGeoRaster->nPyramidLevels = atoi(pszFetched); } // ------------------------------------------------------------------- // Return a new Dataset // ------------------------------------------------------------------- return (GDALDataset*) poGRD; } // --------------------------------------------------------------------------- // CreateCopy() // --------------------------------------------------------------------------- GDALDataset *GeoRasterDataset::CreateCopy( const char* pszFilename, GDALDataset* poSrcDS, int bStrict, char** papszOptions, GDALProgressFunc pfnProgress, void* pProgressData ) { (void) bStrict; int nBands = poSrcDS->GetRasterCount(); if (nBands == 0) { CPLError( CE_Failure, CPLE_NotSupported, "GeoRaster driver does not support source dataset with zero band.\n"); return NULL; } GDALRasterBand* poBand = poSrcDS->GetRasterBand( 1 ); GDALDataType eType = poBand->GetRasterDataType(); // ----------------------------------------------------------- // Create a GeoRaster on the server or select one to overwrite // ----------------------------------------------------------- GeoRasterDataset *poDstDS = (GeoRasterDataset *) GeoRasterDataset::Create( pszFilename, poSrcDS->GetRasterXSize(), poSrcDS->GetRasterYSize(), poSrcDS->GetRasterCount(), eType, papszOptions ); if( poDstDS == NULL ) { return NULL; } // ----------------------------------------------------------- // Copy information to the dataset // ----------------------------------------------------------- double adfTransform[6]; if ( poSrcDS->GetGeoTransform( adfTransform ) == CE_None ) { if ( ! ( adfTransform[0] == 0.0 && adfTransform[1] == 1.0 && adfTransform[2] == 0.0 && adfTransform[3] == 0.0 && adfTransform[4] == 0.0 && adfTransform[5] == 1.0 ) ) { poDstDS->SetGeoTransform( adfTransform ); if( ! poDstDS->bForcedSRID ) /* forced by create option SRID */ { poDstDS->SetProjection( poSrcDS->GetProjectionRef() ); } } } // -------------------------------------------------------------------- // Copy GCPs // -------------------------------------------------------------------- if( poSrcDS->GetGCPCount() > 0 ) { poDstDS->SetGCPs( poSrcDS->GetGCPCount(), poSrcDS->GetGCPs(), poSrcDS->GetGCPProjection() ); } // -------------------------------------------------------------------- // Copy RPC // -------------------------------------------------------------------- char **papszRPCMetadata = GDALGetMetadata( poSrcDS, "RPC" ); if ( papszRPCMetadata != NULL ) { poDstDS->poGeoRaster->phRPC = (GDALRPCInfo*) VSIMalloc( sizeof(GDALRPCInfo) ); GDALExtractRPCInfo( papszRPCMetadata, poDstDS->poGeoRaster->phRPC ); } // -------------------------------------------------------------------- // Copy information to the raster bands // -------------------------------------------------------------------- int bHasNoDataValue = FALSE; double dfNoDataValue = 0.0; double dfMin = 0.0; double dfMax = 0.0; double dfStdDev = 0.0; double dfMean = 0.0; double dfMedian = 0.0; double dfMode = 0.0; int iBand = 0; for( iBand = 1; iBand <= poSrcDS->GetRasterCount(); iBand++ ) { GDALRasterBand* poSrcBand = poSrcDS->GetRasterBand( iBand ); GeoRasterRasterBand* poDstBand = (GeoRasterRasterBand*) poDstDS->GetRasterBand( iBand ); // ---------------------------------------------------------------- // Copy Color Table // ---------------------------------------------------------------- GDALColorTable* poColorTable = poSrcBand->GetColorTable(); if( poColorTable ) { poDstBand->SetColorTable( poColorTable ); } // ---------------------------------------------------------------- // Copy statistics information, without median and mode. // ---------------------------------------------------------------- if( poSrcBand->GetStatistics( false, false, &dfMin, &dfMax, &dfMean, &dfStdDev ) == CE_None ) { poDstBand->SetStatistics( dfMin, dfMax, dfMean, dfStdDev ); /* That will not be recorded in the GeoRaster metadata since it * doesn't have median and mode, so those values are only useful * at runtime. */ } // ---------------------------------------------------------------- // Copy statistics metadata information, including median and mode. // ---------------------------------------------------------------- const char *pszMin = poSrcBand->GetMetadataItem( "STATISTICS_MINIMUM" ); const char *pszMax = poSrcBand->GetMetadataItem( "STATISTICS_MAXIMUM" ); const char *pszMean = poSrcBand->GetMetadataItem( "STATISTICS_MEAN" ); const char *pszMedian = poSrcBand->GetMetadataItem( "STATISTICS_MEDIAN" ); const char *pszMode = poSrcBand->GetMetadataItem( "STATISTICS_MODE" ); const char *pszStdDev = poSrcBand->GetMetadataItem( "STATISTICS_STDDEV" ); const char *pszSkipFX = poSrcBand->GetMetadataItem( "STATISTICS_SKIPFACTORX" ); const char *pszSkipFY = poSrcBand->GetMetadataItem( "STATISTICS_SKIPFACTORY" ); if ( pszMin != NULL && pszMax != NULL && pszMean != NULL && pszMedian != NULL && pszMode != NULL && pszStdDev != NULL ) { dfMin = CPLScanDouble( pszMin, MAX_DOUBLE_STR_REP ); dfMax = CPLScanDouble( pszMax, MAX_DOUBLE_STR_REP ); dfMean = CPLScanDouble( pszMean, MAX_DOUBLE_STR_REP ); dfMedian = CPLScanDouble( pszMedian, MAX_DOUBLE_STR_REP ); dfMode = CPLScanDouble( pszMode, MAX_DOUBLE_STR_REP ); if ( ! ( ( dfMin > dfMax ) || ( dfMean > dfMax ) || ( dfMean < dfMin ) || ( dfMedian > dfMax ) || ( dfMedian < dfMin ) || ( dfMode > dfMax ) || ( dfMode < dfMin ) ) ) { if ( ! pszSkipFX ) { pszSkipFX = pszSkipFY != NULL ? pszSkipFY : "1"; } poDstBand->poGeoRaster->SetStatistics( iBand, pszMin, pszMax, pszMean, pszMedian, pszMode, pszStdDev, pszSkipFX ); } } // ---------------------------------------------------------------- // Copy Raster Attribute Table (RAT) // ---------------------------------------------------------------- GDALRasterAttributeTableH poRAT = GDALGetDefaultRAT( poSrcBand ); if( poRAT != NULL ) { poDstBand->SetDefaultRAT( (GDALRasterAttributeTable*) poRAT ); } // ---------------------------------------------------------------- // Copy NoData Value // ---------------------------------------------------------------- dfNoDataValue = poSrcBand->GetNoDataValue( &bHasNoDataValue ); if( bHasNoDataValue ) { poDstBand->SetNoDataValue( dfNoDataValue ); } } // -------------------------------------------------------------------- // Copy actual imagery. // -------------------------------------------------------------------- int nXSize = poDstDS->GetRasterXSize(); int nYSize = poDstDS->GetRasterYSize(); int nBlockXSize = 0; int nBlockYSize = 0; poDstDS->GetRasterBand( 1 )->GetBlockSize( &nBlockXSize, &nBlockYSize ); // -------------------------------------------------------------------- // JP2-F has one block with full image size. Use tile size instead // -------------------------------------------------------------------- const char* pszFetched = CSLFetchNameValue( papszOptions, "COMPRESS" ); if( pszFetched != NULL && EQUAL( pszFetched, "JP2-F" ) ) { nBlockXSize = DEFAULT_JP2_TILE_COLUMNS; nBlockYSize = DEFAULT_JP2_TILE_ROWS; pszFetched = CSLFetchNameValue( papszOptions, "JP2_BLOCKXSIZE" ); if( pszFetched != NULL ) { nBlockXSize = atoi( pszFetched ); } pszFetched = CSLFetchNameValue( papszOptions, "JP2_BLOCKYSIZE" ); if( pszFetched != NULL ) { nBlockYSize = atoi( pszFetched ); } } // -------------------------------------------------------------------- // Allocate memory buffer to read one block from one band // -------------------------------------------------------------------- void *pData = VSI_MALLOC3_VERBOSE( nBlockXSize, nBlockYSize, GDALGetDataTypeSizeBytes(eType) ); if( pData == NULL ) { delete poDstDS; return NULL; } int iYOffset = 0; int iXOffset = 0; int iXBlock = 0; int iYBlock = 0; int nBlockCols = 0; int nBlockRows = 0; CPLErr eErr = CE_None; int nPixelSize = GDALGetDataTypeSize( poSrcDS->GetRasterBand(1)->GetRasterDataType() ) / 8; if( EQUAL( poDstDS->poGeoRaster->sCompressionType.c_str(), "JPEG-F" ) && nBlockXSize == nXSize && nBlockYSize == nYSize ) { // -------------------------------------------------------------------- // Load JPEG avoiding decompression/compression - direct copy // -------------------------------------------------------------------- const char* pszDriverName = poSrcDS->GetDriverName(); if ( EQUAL( pszDriverName, "JPEG" ) ) { char** papszFileList = poSrcDS->GetFileList(); if ( poDstDS->JPEG_CopyDirect( papszFileList[0], pfnProgress, pProgressData ) ) { CPLDebug("GEOR","JPEG Direct copy succeed"); } } } else if( EQUAL( poDstDS->poGeoRaster->sCompressionType.c_str(), "JP2-F" ) ) { // -------------------------------------------------------------------- // Load JP2K avoiding decompression/compression - direct copy // -------------------------------------------------------------------- boolean bJP2CopyDirectSucceed = false; const char* pszDriverName = poSrcDS->GetDriverName(); int nJP2Resolution = -1; if ( EQUAL( pszDriverName, "JP2OpenJPEG" ) && poSrcDS->GetRasterBand(1)->GetColorTable() == NULL ) { // --------------------------------------------------------------- // Try to load the JP2 file directly // --------------------------------------------------------------- char** papszFileList = poSrcDS->GetFileList(); bJP2CopyDirectSucceed = poDstDS->JP2_CopyDirect( papszFileList[0], &nJP2Resolution, pfnProgress, pProgressData ); } if( ! bJP2CopyDirectSucceed ) { // --------------------------------------------------------------- // Use VSIOCILOB to load using a resident JP2 driver // --------------------------------------------------------------- poDstDS->JP2_CreateCopy( poSrcDS, /* JP2 dataset */ papszOptions, /* options list */ &nJP2Resolution, /* returned resolution */ pfnProgress, /* progress function */ pProgressData ); /* progress data */ } // Number of pyramid levels is the number of resolutions - 1 poDstDS->poGeoRaster->SetMaxLevel( MAX( 1, nJP2Resolution - 1 ) ); } else if( poDstDS->poGeoRaster->nBandBlockSize == 1) { // ---------------------------------------------------------------- // Band order // ---------------------------------------------------------------- int nBandCount = poSrcDS->GetRasterCount(); for( iBand = 1; iBand <= nBandCount; iBand++ ) { GDALRasterBand *poSrcBand = poSrcDS->GetRasterBand( iBand ); GDALRasterBand *poDstBand = poDstDS->GetRasterBand( iBand ); for( iYOffset = 0, iYBlock = 0; iYOffset < nYSize; iYOffset += nBlockYSize, iYBlock++ ) { for( iXOffset = 0, iXBlock = 0; iXOffset < nXSize; iXOffset += nBlockXSize, iXBlock++ ) { nBlockCols = MIN( nBlockXSize, nXSize - iXOffset ); nBlockRows = MIN( nBlockYSize, nYSize - iYOffset ); eErr = poSrcBand->RasterIO( GF_Read, iXOffset, iYOffset, nBlockCols, nBlockRows, pData, nBlockCols, nBlockRows, eType, nPixelSize, nPixelSize * nBlockXSize, NULL ); if( eErr != CE_None ) { return NULL; } eErr = poDstBand->WriteBlock( iXBlock, iYBlock, pData ); if( eErr != CE_None ) { return NULL; } } if( ( eErr == CE_None ) && ( ! pfnProgress( ( ( iBand - 1) / (float) nBandCount ) + ( iYOffset + nBlockRows ) / (float) (nYSize * nBandCount), NULL, pProgressData ) ) ) { eErr = CE_Failure; CPLError( CE_Failure, CPLE_UserInterrupt, "User terminated CreateCopy()" ); } } } } else { // ---------------------------------------------------------------- // Block order // ---------------------------------------------------------------- poDstDS->poGeoRaster->SetWriteOnly( true ); for( iYOffset = 0, iYBlock = 0; iYOffset < nYSize; iYOffset += nBlockYSize, iYBlock++ ) { for( iXOffset = 0, iXBlock = 0; iXOffset < nXSize; iXOffset += nBlockXSize, iXBlock++ ) { nBlockCols = MIN( nBlockXSize, nXSize - iXOffset ); nBlockRows = MIN( nBlockYSize, nYSize - iYOffset ); for( iBand = 1; iBand <= poSrcDS->GetRasterCount(); iBand++ ) { GDALRasterBand *poSrcBand = poSrcDS->GetRasterBand( iBand ); GDALRasterBand *poDstBand = poDstDS->GetRasterBand( iBand ); eErr = poSrcBand->RasterIO( GF_Read, iXOffset, iYOffset, nBlockCols, nBlockRows, pData, nBlockCols, nBlockRows, eType, nPixelSize, nPixelSize * nBlockXSize, NULL ); if( eErr != CE_None ) { return NULL; } eErr = poDstBand->WriteBlock( iXBlock, iYBlock, pData ); if( eErr != CE_None ) { return NULL; } } } if( ( eErr == CE_None ) && ( ! pfnProgress( ( iYOffset + nBlockRows ) / (double) nYSize, NULL, pProgressData ) ) ) { eErr = CE_Failure; CPLError( CE_Failure, CPLE_UserInterrupt, "User terminated CreateCopy()" ); } } } CPLFree( pData ); // -------------------------------------------------------------------- // Finalize // -------------------------------------------------------------------- poDstDS->FlushCache(); if( pfnProgress ) { CPLDebug("GEOR", "Output dataset: (georaster:%s/%s@%s,%s,%lld) on %s%s,%s", poDstDS->poGeoRaster->poConnection->GetUser(), poDstDS->poGeoRaster->poConnection->GetPassword(), poDstDS->poGeoRaster->poConnection->GetServer(), poDstDS->poGeoRaster->sDataTable.c_str(), poDstDS->poGeoRaster->nRasterId, poDstDS->poGeoRaster->sSchema.c_str(), poDstDS->poGeoRaster->sTable.c_str(), poDstDS->poGeoRaster->sColumn.c_str() ); } return poDstDS; } // --------------------------------------------------------------------------- // IRasterIO() // --------------------------------------------------------------------------- CPLErr GeoRasterDataset::IRasterIO( GDALRWFlag eRWFlag, int nXOff, int nYOff, int nXSize, int nYSize, void *pData, int nBufXSize, int nBufYSize, GDALDataType eBufType, int nBandCount, int *panBandMap, GSpacing nPixelSpace, GSpacing nLineSpace, GSpacing nBandSpace, GDALRasterIOExtraArg* psExtraArg ) { if( EQUAL( poGeoRaster->sCompressionType.c_str(), "JP2-F" ) ) { if( poJP2Dataset ) { return poJP2Dataset->RasterIO( eRWFlag, nXOff, nYOff, nXSize, nYSize, pData, nBufXSize, nBufYSize, eBufType, nBandCount, panBandMap, nPixelSpace, nLineSpace, nBandSpace, psExtraArg ); } else { return CE_Failure; } } else { if( poGeoRaster->nBandBlockSize > 1 ) { return GDALDataset::BlockBasedRasterIO( eRWFlag, nXOff, nYOff, nXSize, nYSize, pData, nBufXSize, nBufYSize, eBufType, nBandCount, panBandMap, nPixelSpace, nLineSpace, nBandSpace, psExtraArg ); } else { return GDALDataset::IRasterIO( eRWFlag, nXOff, nYOff, nXSize, nYSize, pData, nBufXSize, nBufYSize, eBufType, nBandCount, panBandMap, nPixelSpace, nLineSpace, nBandSpace, psExtraArg ); } } } // --------------------------------------------------------------------------- // FlushCache() // --------------------------------------------------------------------------- void GeoRasterDataset::FlushCache() { GDALDataset::FlushCache(); } // --------------------------------------------------------------------------- // GetGeoTransform() // --------------------------------------------------------------------------- CPLErr GeoRasterDataset::GetGeoTransform( double *padfTransform ) { if( poGeoRaster->phRPC ) { return CE_Failure; } if( poGeoRaster->nSRID == 0 ) { return CE_Failure; } memcpy( padfTransform, adfGeoTransform, sizeof(double) * 6 ); bGeoTransform = true; return CE_None; } // --------------------------------------------------------------------------- // GetProjectionRef() // --------------------------------------------------------------------------- const char* GeoRasterDataset::GetProjectionRef( void ) { if( poGeoRaster->phRPC ) { return ""; } if( ! poGeoRaster->bIsReferenced ) { return ""; } if( poGeoRaster->nSRID == UNKNOWN_CRS || poGeoRaster->nSRID == 0 ) { return ""; } if( pszProjection ) { return pszProjection; } OGRSpatialReference oSRS; // -------------------------------------------------------------------- // Check if the SRID is a valid EPSG code // -------------------------------------------------------------------- CPLPushErrorHandler( CPLQuietErrorHandler ); if( oSRS.importFromEPSG( poGeoRaster->nSRID ) == OGRERR_NONE ) { /* * Ignores the WKT from Oracle and use the one from GDAL's * EPSG tables. That would ensure that other drivers/software * will recognize the parameters. */ if( oSRS.exportToWkt( &pszProjection ) == OGRERR_NONE ) { CPLPopErrorHandler(); return pszProjection; } } CPLPopErrorHandler(); // -------------------------------------------------------------------- // Try to interpreter the WKT text // -------------------------------------------------------------------- poGeoRaster->QueryWKText(); char* pszWKText = CPLStrdup( poGeoRaster->sWKText ); if( ! ( oSRS.importFromWkt( &pszWKText ) == OGRERR_NONE && oSRS.GetRoot() ) ) { return pszWKText; } // ---------------------------------------------------------------- // Decorate with Authority name // ---------------------------------------------------------------- if( strlen(poGeoRaster->sAuthority) > 0 ) { oSRS.SetAuthority(oSRS.GetRoot()->GetValue(), poGeoRaster->sAuthority.c_str(), poGeoRaster->nSRID); } int nSpher = OWParseEPSG( oSRS.GetAttrValue("GEOGCS|DATUM|SPHEROID") ); if( nSpher > 0 ) { oSRS.SetAuthority( "GEOGCS|DATUM|SPHEROID", "EPSG", nSpher ); } int nDatum = OWParseEPSG( oSRS.GetAttrValue("GEOGCS|DATUM") ); if( nDatum > 0 ) { oSRS.SetAuthority( "GEOGCS|DATUM", "EPSG", nDatum ); } // ---------------------------------------------------------------- // Checks for Projection info // ---------------------------------------------------------------- const char *pszProjName = oSRS.GetAttrValue( "PROJECTION" ); if( pszProjName ) { int nProj = OWParseEPSG( pszProjName ); // ---------------------------------------------------------------- // Decorate with EPSG Authority // ---------------------------------------------------------------- if( nProj > 0 ) { oSRS.SetAuthority( "PROJECTION", "EPSG", nProj ); } // ---------------------------------------------------------------- // Translate projection names to GDAL's standards // ---------------------------------------------------------------- if ( EQUAL( pszProjName, "Transverse Mercator" ) ) { oSRS.SetProjection( SRS_PT_TRANSVERSE_MERCATOR ); } else if ( EQUAL( pszProjName, "Albers Conical Equal Area" ) ) { oSRS.SetProjection( SRS_PT_ALBERS_CONIC_EQUAL_AREA ); } else if ( EQUAL( pszProjName, "Azimuthal Equidistant" ) ) { oSRS.SetProjection( SRS_PT_AZIMUTHAL_EQUIDISTANT ); } else if ( EQUAL( pszProjName, "Miller Cylindrical" ) ) { oSRS.SetProjection( SRS_PT_MILLER_CYLINDRICAL ); } else if ( EQUAL( pszProjName, "Hotine Oblique Mercator" ) ) { oSRS.SetProjection( SRS_PT_HOTINE_OBLIQUE_MERCATOR ); } else if ( EQUAL( pszProjName, "Wagner IV" ) ) { oSRS.SetProjection( SRS_PT_WAGNER_IV ); } else if ( EQUAL( pszProjName, "Wagner VII" ) ) { oSRS.SetProjection( SRS_PT_WAGNER_VII ); } else if ( EQUAL( pszProjName, "Eckert IV" ) ) { oSRS.SetProjection( SRS_PT_ECKERT_IV ); } else if ( EQUAL( pszProjName, "Eckert VI" ) ) { oSRS.SetProjection( SRS_PT_ECKERT_VI ); } else if ( EQUAL( pszProjName, "New Zealand Map Grid" ) ) { oSRS.SetProjection( SRS_PT_NEW_ZEALAND_MAP_GRID ); } else if ( EQUAL( pszProjName, "Lambert Conformal Conic" ) ) { oSRS.SetProjection( SRS_PT_LAMBERT_CONFORMAL_CONIC_1SP ); } else if ( EQUAL( pszProjName, "Lambert Azimuthal Equal Area" ) ) { oSRS.SetProjection( SRS_PT_LAMBERT_AZIMUTHAL_EQUAL_AREA ); } else if ( EQUAL( pszProjName, "Van der Grinten" ) ) { oSRS.SetProjection( SRS_PT_VANDERGRINTEN ); } else if ( EQUAL( pszProjName, "Lambert Conformal Conic (Belgium 1972)" ) ) { oSRS.SetProjection( SRS_PT_LAMBERT_CONFORMAL_CONIC_2SP_BELGIUM ); } else if ( EQUAL( pszProjName, "Cylindrical Equal Area" ) ) { oSRS.SetProjection( SRS_PT_CYLINDRICAL_EQUAL_AREA ); } else if ( EQUAL( pszProjName, "Interrupted Goode Homolosine" ) ) { oSRS.SetProjection( SRS_PT_GOODE_HOMOLOSINE ); } } oSRS.exportToWkt( &pszProjection ); return pszProjection; } // --------------------------------------------------------------------------- // SetGeoTransform() // --------------------------------------------------------------------------- CPLErr GeoRasterDataset::SetGeoTransform( double *padfTransform ) { memcpy( adfGeoTransform, padfTransform, sizeof( double ) * 6 ); poGeoRaster->dfXCoefficient[0] = adfGeoTransform[1]; poGeoRaster->dfXCoefficient[1] = adfGeoTransform[2]; poGeoRaster->dfXCoefficient[2] = adfGeoTransform[0]; poGeoRaster->dfYCoefficient[0] = adfGeoTransform[4]; poGeoRaster->dfYCoefficient[1] = adfGeoTransform[5]; poGeoRaster->dfYCoefficient[2] = adfGeoTransform[3]; bGeoTransform = true; return CE_None; } // --------------------------------------------------------------------------- // SetProjection() // --------------------------------------------------------------------------- CPLErr GeoRasterDataset::SetProjection( const char *pszProjString ) { OGRSpatialReference oSRS; char* pszWKT = CPLStrdup( pszProjString ); OGRErr eOGRErr = oSRS.importFromWkt( &pszWKT ); if( eOGRErr != OGRERR_NONE ) { poGeoRaster->SetGeoReference( UNKNOWN_CRS ); return CE_Failure; } // -------------------------------------------------------------------- // Try to extract EPGS authority code // -------------------------------------------------------------------- const char *pszAuthName = NULL; const char *pszAuthCode = NULL; if( oSRS.IsGeographic() ) { pszAuthName = oSRS.GetAuthorityName( "GEOGCS" ); pszAuthCode = oSRS.GetAuthorityCode( "GEOGCS" ); } else if( oSRS.IsProjected() ) { pszAuthName = oSRS.GetAuthorityName( "PROJCS" ); pszAuthCode = oSRS.GetAuthorityCode( "PROJCS" ); } if( pszAuthName != NULL && pszAuthCode != NULL ) { if( EQUAL( pszAuthName, "ORACLE" ) || EQUAL( pszAuthName, "EPSG" ) ) { poGeoRaster->SetGeoReference( atoi( pszAuthCode ) ); return CE_None; } } // ---------------------------------------------------------------- // Convert SRS into old style format (SF-SQL 1.0) // ---------------------------------------------------------------- OGRSpatialReference *poSRS2 = oSRS.Clone(); poSRS2->StripCTParms(); double dfAngularUnits = poSRS2->GetAngularUnits( NULL ); if( fabs(dfAngularUnits - 0.0174532925199433) < 0.0000000000000010 ) { /* match the precision used on Oracle for that particular value */ poSRS2->SetAngularUnits( "Decimal Degree", 0.0174532925199433 ); } char* pszCloneWKT = NULL; if( poSRS2->exportToWkt( &pszCloneWKT ) != OGRERR_NONE ) { delete poSRS2; return CE_Failure; } const char *pszProjName = poSRS2->GetAttrValue( "PROJECTION" ); if( pszProjName ) { // ---------------------------------------------------------------- // Translate projection names to Oracle's standards // ---------------------------------------------------------------- if ( EQUAL( pszProjName, SRS_PT_TRANSVERSE_MERCATOR ) ) { poSRS2->SetProjection( "Transverse Mercator" ); } else if ( EQUAL( pszProjName, SRS_PT_ALBERS_CONIC_EQUAL_AREA ) ) { poSRS2->SetProjection( "Albers Conical Equal Area" ); } else if ( EQUAL( pszProjName, SRS_PT_AZIMUTHAL_EQUIDISTANT ) ) { poSRS2->SetProjection( "Azimuthal Equidistant" ); } else if ( EQUAL( pszProjName, SRS_PT_MILLER_CYLINDRICAL ) ) { poSRS2->SetProjection( "Miller Cylindrical" ); } else if ( EQUAL( pszProjName, SRS_PT_HOTINE_OBLIQUE_MERCATOR ) ) { poSRS2->SetProjection( "Hotine Oblique Mercator" ); } else if ( EQUAL( pszProjName, SRS_PT_WAGNER_IV ) ) { poSRS2->SetProjection( "Wagner IV" ); } else if ( EQUAL( pszProjName, SRS_PT_WAGNER_VII ) ) { poSRS2->SetProjection( "Wagner VII" ); } else if ( EQUAL( pszProjName, SRS_PT_ECKERT_IV ) ) { poSRS2->SetProjection( "Eckert IV" ); } else if ( EQUAL( pszProjName, SRS_PT_ECKERT_VI ) ) { poSRS2->SetProjection( "Eckert VI" ); } else if ( EQUAL( pszProjName, SRS_PT_NEW_ZEALAND_MAP_GRID ) ) { poSRS2->SetProjection( "New Zealand Map Grid" ); } else if ( EQUAL( pszProjName, SRS_PT_LAMBERT_CONFORMAL_CONIC_1SP ) ) { poSRS2->SetProjection( "Lambert Conformal Conic" ); } else if ( EQUAL( pszProjName, SRS_PT_LAMBERT_AZIMUTHAL_EQUAL_AREA ) ) { poSRS2->SetProjection( "Lambert Azimuthal Equal Area" ); } else if ( EQUAL( pszProjName, SRS_PT_VANDERGRINTEN ) ) { poSRS2->SetProjection( "Van der Grinten" ); } else if ( EQUAL( pszProjName, SRS_PT_LAMBERT_CONFORMAL_CONIC_2SP_BELGIUM ) ) { poSRS2->SetProjection( "Lambert Conformal Conic (Belgium 1972)" ); } else if ( EQUAL( pszProjName, SRS_PT_CYLINDRICAL_EQUAL_AREA ) ) { poSRS2->SetProjection( "Cylindrical Equal Area" ); } else if ( EQUAL( pszProjName, SRS_PT_GOODE_HOMOLOSINE ) ) { poSRS2->SetProjection( "Interrupted Goode Homolosine" ); } // ---------------------------------------------------------------- // Translate projection's parameters to Oracle's standards // ---------------------------------------------------------------- char* pszStart = NULL; CPLFree( pszCloneWKT ); if( poSRS2->exportToWkt( &pszCloneWKT ) != OGRERR_NONE ) { delete poSRS2; return CE_Failure; } if( ( pszStart = strstr(pszCloneWKT, SRS_PP_AZIMUTH) ) != NULL ) { strncpy( pszStart, "Azimuth", strlen(SRS_PP_AZIMUTH) ); } if( ( pszStart = strstr(pszCloneWKT, SRS_PP_CENTRAL_MERIDIAN) ) != NULL ) { strncpy( pszStart, "Central_Meridian", strlen(SRS_PP_CENTRAL_MERIDIAN) ); } if( ( pszStart = strstr(pszCloneWKT, SRS_PP_FALSE_EASTING) ) != NULL ) { strncpy( pszStart, "False_Easting", strlen(SRS_PP_FALSE_EASTING) ); } if( ( pszStart = strstr(pszCloneWKT, SRS_PP_FALSE_NORTHING) ) != NULL ) { strncpy( pszStart, "False_Northing", strlen(SRS_PP_FALSE_NORTHING) ); } if( ( pszStart = strstr(pszCloneWKT, SRS_PP_LATITUDE_OF_CENTER) ) != NULL ) { strncpy( pszStart, "Latitude_Of_Center", strlen(SRS_PP_LATITUDE_OF_CENTER) ); } if( ( pszStart = strstr(pszCloneWKT, SRS_PP_LATITUDE_OF_ORIGIN) ) != NULL ) { strncpy( pszStart, "Latitude_Of_Origin", strlen(SRS_PP_LATITUDE_OF_ORIGIN) ); } if( ( pszStart = strstr(pszCloneWKT, SRS_PP_LONGITUDE_OF_CENTER) ) != NULL ) { strncpy( pszStart, "Longitude_Of_Center", strlen(SRS_PP_LONGITUDE_OF_CENTER) ); } if( ( pszStart = strstr(pszCloneWKT, SRS_PP_PSEUDO_STD_PARALLEL_1) ) != NULL ) { strncpy( pszStart, "Pseudo_Standard_Parallel_1", strlen(SRS_PP_PSEUDO_STD_PARALLEL_1) ); } if( ( pszStart = strstr(pszCloneWKT, SRS_PP_SCALE_FACTOR) ) != NULL ) { strncpy( pszStart, "Scale_Factor", strlen(SRS_PP_SCALE_FACTOR) ); } if( ( pszStart = strstr(pszCloneWKT, SRS_PP_STANDARD_PARALLEL_1) ) != NULL ) { strncpy( pszStart, "Standard_Parallel_1", strlen(SRS_PP_STANDARD_PARALLEL_1) ); } if( ( pszStart = strstr(pszCloneWKT, SRS_PP_STANDARD_PARALLEL_2) ) != NULL ) { strncpy( pszStart, "Standard_Parallel_2", strlen(SRS_PP_STANDARD_PARALLEL_2) ); } if( ( pszStart = strstr(pszCloneWKT, SRS_PP_STANDARD_PARALLEL_2) ) != NULL ) { strncpy( pszStart, "Standard_Parallel_2", strlen(SRS_PP_STANDARD_PARALLEL_2) ); } // ---------------------------------------------------------------- // Fix Unit name // ---------------------------------------------------------------- if( ( pszStart = strstr(pszCloneWKT, "metre") ) != NULL ) { strncpy( pszStart, SRS_UL_METER, strlen(SRS_UL_METER) ); } } // -------------------------------------------------------------------- // Tries to find a SRID compatible with the WKT // -------------------------------------------------------------------- OWConnection* poConnection = poGeoRaster->poConnection; OWStatement* poStmt = NULL; int nNewSRID = 0; const char *pszFuncName = "FIND_GEOG_CRS"; if( poSRS2->IsProjected() ) { pszFuncName = "FIND_PROJ_CRS"; } poStmt = poConnection->CreateStatement( CPLSPrintf( "DECLARE\n" " LIST SDO_SRID_LIST;" "BEGIN\n" " SELECT SDO_CS.%s('%s', null) into LIST FROM DUAL;\n" " IF LIST.COUNT() > 0 then\n" " SELECT LIST(1) into :out from dual;\n" " ELSE\n" " SELECT 0 into :out from dual;\n" " END IF;\n" "END;", pszFuncName, pszCloneWKT ) ); poStmt->BindName( ":out", &nNewSRID ); CPLPushErrorHandler( CPLQuietErrorHandler ); if( poStmt->Execute() ) { CPLPopErrorHandler(); if ( nNewSRID > 0 ) { poGeoRaster->SetGeoReference( nNewSRID ); CPLFree( pszCloneWKT ); return CE_None; } } // -------------------------------------------------------------------- // Search by simplified WKT or insert it as a user defined SRS // -------------------------------------------------------------------- int nCounter = 0; poStmt = poConnection->CreateStatement( CPLSPrintf( "SELECT COUNT(*) FROM MDSYS.CS_SRS WHERE WKTEXT = '%s'", pszCloneWKT)); poStmt->Define( &nCounter ); CPLPushErrorHandler( CPLQuietErrorHandler ); if( poStmt->Execute() && nCounter > 0 ) { poStmt = poConnection->CreateStatement( CPLSPrintf( "SELECT SRID FROM MDSYS.CS_SRS WHERE WKTEXT = '%s'", pszCloneWKT)); poStmt->Define( &nNewSRID ); if( poStmt->Execute() ) { CPLPopErrorHandler(); poGeoRaster->SetGeoReference( nNewSRID ); CPLFree( pszCloneWKT ); return CE_None; } } CPLPopErrorHandler(); poStmt = poConnection->CreateStatement( CPLSPrintf( "DECLARE\n" " MAX_SRID NUMBER := 0;\n" "BEGIN\n" " SELECT MAX(SRID) INTO MAX_SRID FROM MDSYS.CS_SRS;\n" " MAX_SRID := MAX_SRID + 1;\n" " INSERT INTO MDSYS.CS_SRS (SRID, WKTEXT, CS_NAME)\n" " VALUES (MAX_SRID, '%s', '%s');\n" " SELECT MAX_SRID INTO :out FROM DUAL;\n" "END;", pszCloneWKT, oSRS.GetRoot()->GetChild(0)->GetValue() ) ); poStmt->BindName( ":out", &nNewSRID ); CPLErr eError = CE_None; CPLPushErrorHandler( CPLQuietErrorHandler ); if( poStmt->Execute() ) { CPLPopErrorHandler(); poGeoRaster->SetGeoReference( nNewSRID ); } else { CPLPopErrorHandler(); poGeoRaster->SetGeoReference( UNKNOWN_CRS ); CPLError( CE_Warning, CPLE_UserInterrupt, "Insufficient privileges to insert reference system to " "table MDSYS.CS_SRS." ); eError = CE_Warning; } CPLFree( pszCloneWKT ); return eError; } /************************************************************************/ /* GetMetadataDomainList() */ /************************************************************************/ char **GeoRasterDataset::GetMetadataDomainList() { return BuildMetadataDomainList(GDALDataset::GetMetadataDomainList(), TRUE, "SUBDATASETS", NULL); } // --------------------------------------------------------------------------- // GetMetadata() // --------------------------------------------------------------------------- char **GeoRasterDataset::GetMetadata( const char *pszDomain ) { if( pszDomain != NULL && STARTS_WITH_CI(pszDomain, "SUBDATASETS") ) return papszSubdatasets; else return GDALDataset::GetMetadata( pszDomain ); } // --------------------------------------------------------------------------- // Delete() // --------------------------------------------------------------------------- CPLErr GeoRasterDataset::Delete( const char* pszFilename ) { (void) pszFilename; /*** GeoRasterDataset* poGRD = NULL; poGRD = (GeoRasterDataset*) GDALOpen( pszFilename, GA_Update ); if( ! poGRD ) { return CE_Failure; } if( ! poGRD->poGeoRaster->Delete() ) { return CE_Failure; } ***/ return CE_None; } // --------------------------------------------------------------------------- // SetSubdatasets() // --------------------------------------------------------------------------- void GeoRasterDataset::SetSubdatasets( GeoRasterWrapper* poGRW ) { OWConnection* poConnection = poGRW->poConnection; OWStatement* poStmt = NULL; // ----------------------------------------------------------- // List all the GeoRaster Tables of that User/Database // ----------------------------------------------------------- if( poGRW->sTable.empty() && poGRW->sColumn.empty() ) { poStmt = poConnection->CreateStatement( "SELECT DISTINCT TABLE_NAME, OWNER FROM ALL_SDO_GEOR_SYSDATA\n" " ORDER BY TABLE_NAME ASC" ); char szTable[OWNAME]; char szOwner[OWNAME]; poStmt->Define( szTable ); poStmt->Define( szOwner ); if( poStmt->Execute() ) { int nCount = 1; do { papszSubdatasets = CSLSetNameValue( papszSubdatasets, CPLSPrintf( "SUBDATASET_%d_NAME", nCount ), CPLSPrintf( "geor:%s/%s@%s,%s.%s", poConnection->GetUser(), poConnection->GetPassword(), poConnection->GetServer(), szOwner, szTable ) ); papszSubdatasets = CSLSetNameValue( papszSubdatasets, CPLSPrintf( "SUBDATASET_%d_DESC", nCount ), CPLSPrintf( "%s.Table=%s", szOwner, szTable ) ); nCount++; } while( poStmt->Fetch() ); } return; } // ----------------------------------------------------------- // List all the GeoRaster Columns of that Table // ----------------------------------------------------------- if( ! poGRW->sTable.empty() && poGRW->sColumn.empty() ) { poStmt = poConnection->CreateStatement( CPLSPrintf( "SELECT DISTINCT COLUMN_NAME, OWNER FROM ALL_SDO_GEOR_SYSDATA\n" " WHERE TABLE_NAME = UPPER('%s')\n" " ORDER BY COLUMN_NAME ASC", poGRW->sTable.c_str() ) ); char szColumn[OWNAME]; char szOwner[OWNAME]; poStmt->Define( szColumn ); poStmt->Define( szOwner ); if( poStmt->Execute() ) { int nCount = 1; do { papszSubdatasets = CSLSetNameValue( papszSubdatasets, CPLSPrintf( "SUBDATASET_%d_NAME", nCount ), CPLSPrintf( "geor:%s/%s@%s,%s.%s,%s", poConnection->GetUser(), poConnection->GetPassword(), poConnection->GetServer(), szOwner, poGRW->sTable.c_str(), szColumn ) ); papszSubdatasets = CSLSetNameValue( papszSubdatasets, CPLSPrintf( "SUBDATASET_%d_DESC", nCount ), CPLSPrintf( "Table=%s.%s Column=%s", szOwner, poGRW->sTable.c_str(), szColumn ) ); nCount++; } while( poStmt->Fetch() ); } return; } // ----------------------------------------------------------- // List all the rows that contains GeoRaster on Table/Column/Where // ----------------------------------------------------------- CPLString osAndWhere = ""; if( ! poGRW->sWhere.empty() ) { osAndWhere = CPLSPrintf( "AND %s", poGRW->sWhere.c_str() ); } poStmt = poConnection->CreateStatement( CPLSPrintf( "SELECT T.%s.RASTERDATATABLE, T.%s.RASTERID, \n" " extractValue(t.%s.metadata, " "'/georasterMetadata/rasterInfo/dimensionSize[@type=\"ROW\"]/size','%s'),\n" " extractValue(t.%s.metadata, " "'/georasterMetadata/rasterInfo/dimensionSize[@type=\"COLUMN\"]/size','%s'),\n" " extractValue(t.%s.metadata, " "'/georasterMetadata/rasterInfo/dimensionSize[@type=\"BAND\"]/size','%s'),\n" " extractValue(t.%s.metadata, " "'/georasterMetadata/rasterInfo/cellDepth','%s'),\n" " extractValue(t.%s.metadata, " "'/georasterMetadata/spatialReferenceInfo/SRID','%s')\n" " FROM %s%s T\n" " WHERE %s IS NOT NULL %s\n" " ORDER BY T.%s.RASTERDATATABLE ASC,\n" " T.%s.RASTERID ASC", poGRW->sColumn.c_str(), poGRW->sColumn.c_str(), poGRW->sColumn.c_str(), OW_XMLNS, poGRW->sColumn.c_str(), OW_XMLNS, poGRW->sColumn.c_str(), OW_XMLNS, poGRW->sColumn.c_str(), OW_XMLNS, poGRW->sColumn.c_str(), OW_XMLNS, poGRW->sSchema.c_str(), poGRW->sTable.c_str(), poGRW->sColumn.c_str(), osAndWhere.c_str(), poGRW->sColumn.c_str(), poGRW->sColumn.c_str() ) ); char szDataTable[OWNAME]; char szRasterId[OWNAME]; char szRows[OWNAME]; char szColumns[OWNAME]; char szBands[OWNAME]; char szCellDepth[OWNAME]; char szSRID[OWNAME]; poStmt->Define( szDataTable ); poStmt->Define( szRasterId ); poStmt->Define( szRows ); poStmt->Define( szColumns ); poStmt->Define( szBands ); poStmt->Define( szCellDepth ); poStmt->Define( szSRID ); if( poStmt->Execute() ) { int nCount = 1; do { papszSubdatasets = CSLSetNameValue( papszSubdatasets, CPLSPrintf( "SUBDATASET_%d_NAME", nCount ), CPLSPrintf( "geor:%s/%s@%s,%s,%s", poConnection->GetUser(), poConnection->GetPassword(), poConnection->GetServer(), szDataTable, szRasterId ) ); const char* pszXBands = ""; if( ! EQUAL( szBands, "" ) ) { pszXBands = CPLSPrintf( "x%s", szBands ); } papszSubdatasets = CSLSetNameValue( papszSubdatasets, CPLSPrintf( "SUBDATASET_%d_DESC", nCount ), CPLSPrintf( "[%sx%s%s] CellDepth=%s SRID=%s", szRows, szColumns, pszXBands, szCellDepth, szSRID ) ); nCount++; } while( poStmt->Fetch() ); } } int GeoRasterDataset::GetGCPCount() { if ( poGeoRaster ) { return poGeoRaster->nGCPCount; } return 0; } // --------------------------------------------------------------------------- // SetGCPs() // --------------------------------------------------------------------------- CPLErr GeoRasterDataset::SetGCPs( int nGCPCountIn, const GDAL_GCP *pasGCPListIn, const char *pszGCPProjection ) { if( GetAccess() == GA_Update ) { poGeoRaster->SetGCP( nGCPCountIn, pasGCPListIn ); SetProjection( pszGCPProjection ); } else { CPLError(CE_Failure, CPLE_NotSupported, "SetGCPs() is only supported on GeoRaster insert or update."); return CE_Failure; } return CE_None; } const GDAL_GCP* GeoRasterDataset::GetGCPs() { if( poGeoRaster->nGCPCount > 0 && poGeoRaster->pasGCPList ) { return poGeoRaster->pasGCPList; } return NULL; } // --------------------------------------------------------------------------- // GetGCPProjection() // --------------------------------------------------------------------------- const char* GeoRasterDataset::GetGCPProjection() { if( poGeoRaster && poGeoRaster->nGCPCount > 0 ) return pszProjection; else return ""; } // --------------------------------------------------------------------------- // IBuildOverviews() // --------------------------------------------------------------------------- CPLErr GeoRasterDataset::IBuildOverviews( const char* pszResampling, int nOverviews, int* panOverviewList, int nListBands, int* panBandList, GDALProgressFunc pfnProgress, void* pProgressData ) { (void) panBandList; (void) nListBands; if( EQUAL( poGeoRaster->sCompressionType.c_str(), "JP2-F" ) ) { return CE_None; // Ignore it, JP2 automatically has overviews } // --------------------------------------------------------------- // Can't update on read-only access mode // --------------------------------------------------------------- if( GetAccess() != GA_Update ) { CPLError( CE_Failure, CPLE_AppDefined, "Can't build overviews/pyramids on read-only access." ); return CE_Failure; } // --------------------------------------------------------------- // Uses internal sdo_generatePyramid at PL/SQL? // --------------------------------------------------------------- bool bInternal = true; const char *pszGEOR_INTERNAL_PYR = CPLGetConfigOption( "GEOR_INTERNAL_PYR", "YES" ); if( EQUAL(pszGEOR_INTERNAL_PYR, "NO") ) { bInternal = false; } // ----------------------------------------------------------- // Pyramids applies to the whole dataset not to a specific band // ----------------------------------------------------------- if( nBands < GetRasterCount()) { CPLError( CE_Failure, CPLE_AppDefined, "Invalid GeoRaster Pyramids band selection" ); return CE_Failure; } // --------------------------------------------------------------- // Initialize progress reporting // --------------------------------------------------------------- if( ! pfnProgress( 0.1, NULL, pProgressData ) ) { CPLError( CE_Failure, CPLE_UserInterrupt, "User terminated" ); return CE_Failure; } // --------------------------------------------------------------- // Clear existing overviews // --------------------------------------------------------------- if( nOverviews == 0 ) { poGeoRaster->DeletePyramid(); return CE_None; } // ----------------------------------------------------------- // Pyramids levels can not be treated individually // ----------------------------------------------------------- if( nOverviews > 0 ) { int i; for( i = 1; i < nOverviews; i++ ) { // ----------------------------------------------------------- // Power of 2, starting on 2, e.g. 2, 4, 8, 16, 32, 64, 128 // ----------------------------------------------------------- if( panOverviewList[0] != 2 || ( panOverviewList[i] != panOverviewList[i-1] * 2 ) ) { CPLError( CE_Failure, CPLE_AppDefined, "Invalid GeoRaster Pyramids levels." ); return CE_Failure; } } } // ----------------------------------------------------------- // Re-sampling method: // NN, BILINEAR, AVERAGE4, AVERAGE16 and CUBIC // ----------------------------------------------------------- char szMethod[OWNAME]; if( EQUAL( pszResampling, "NEAREST" ) ) { strcpy( szMethod, "NN" ); } else if( STARTS_WITH_CI(pszResampling, "AVERAGE") ) { strcpy( szMethod, "AVERAGE4" ); } else { CPLError( CE_Failure, CPLE_AppDefined, "Invalid resampling method" ); return CE_Failure; } // ----------------------------------------------------------- // Generate pyramids on poGeoRaster // ----------------------------------------------------------- if( ! poGeoRaster->GeneratePyramid( nOverviews, szMethod, bInternal ) ) { CPLError( CE_Failure, CPLE_AppDefined, "Error generating pyramid" ); return CE_Failure; } // ----------------------------------------------------------- // If Pyramid was done internally on the server exit here // ----------------------------------------------------------- if( bInternal ) { pfnProgress( 1 , NULL, pProgressData ); return CE_None; } // ----------------------------------------------------------- // Load the pyramids data using GDAL methods // ----------------------------------------------------------- CPLErr eErr = CE_None; int i = 0; for( i = 0; i < nBands; i++ ) { GeoRasterRasterBand* poBand = (GeoRasterRasterBand*) papoBands[i]; // ------------------------------------------------------- // Clean up previous overviews // ------------------------------------------------------- int j = 0; if( poBand->nOverviewCount && poBand->papoOverviews ) { for( j = 0; j < poBand->nOverviewCount; j++ ) { delete poBand->papoOverviews[j]; } CPLFree( poBand->papoOverviews ); } // ------------------------------------------------------- // Create new band's overviews list // ------------------------------------------------------- poBand->nOverviewCount = poGeoRaster->nPyramidMaxLevel; poBand->papoOverviews = (GeoRasterRasterBand**) VSIMalloc( sizeof(GeoRasterRasterBand*) * poBand->nOverviewCount ); for( j = 0; j < poBand->nOverviewCount; j++ ) { poBand->papoOverviews[j] = new GeoRasterRasterBand( (GeoRasterDataset*) this, ( i + 1 ), ( j + 1 ) ); } } // ----------------------------------------------------------- // Load band's overviews // ----------------------------------------------------------- for( i = 0; i < nBands; i++ ) { GeoRasterRasterBand* poBand = (GeoRasterRasterBand*) papoBands[i]; void *pScaledProgressData = GDALCreateScaledProgress( i / (double) nBands, ( i + 1) / (double) nBands, pfnProgress, pProgressData ); eErr = GDALRegenerateOverviews( (GDALRasterBandH) poBand, poBand->nOverviewCount, (GDALRasterBandH*) poBand->papoOverviews, pszResampling, GDALScaledProgress, pScaledProgressData ); GDALDestroyScaledProgress( pScaledProgressData ); } return eErr; } // --------------------------------------------------------------------------- // CreateMaskBand() // --------------------------------------------------------------------------- CPLErr GeoRasterDataset::CreateMaskBand( int /*nFlags*/ ) { if( ! poGeoRaster->InitializeMask( DEFAULT_BMP_MASK, poGeoRaster->nRowBlockSize, poGeoRaster->nColumnBlockSize, poGeoRaster->nTotalRowBlocks, poGeoRaster->nTotalColumnBlocks, poGeoRaster->nTotalBandBlocks ) ) { return CE_Failure; } poGeoRaster->bHasBitmapMask = true; return CE_None; } /*****************************************************************************/ /* GDALRegister_GEOR */ /*****************************************************************************/ void CPL_DLL GDALRegister_GEOR() { if( !GDAL_CHECK_VERSION( "GeoRaster driver" ) ) return; if( GDALGetDriverByName( "GeoRaster" ) != NULL ) return; GDALDriver *poDriver = new GDALDriver(); poDriver->SetDescription( "GeoRaster" ); poDriver->SetMetadataItem( GDAL_DCAP_RASTER, "YES" ); poDriver->SetMetadataItem( GDAL_DMD_LONGNAME, "Oracle Spatial GeoRaster" ); poDriver->SetMetadataItem( GDAL_DMD_HELPTOPIC, "frmt_georaster.html" ); poDriver->SetMetadataItem( GDAL_DMD_SUBDATASETS, "YES" ); poDriver->SetMetadataItem( GDAL_DMD_CREATIONDATATYPES, "Byte UInt16 Int16 UInt32 Int32 Float32 " "Float64 CFloat32 CFloat64" ); poDriver->SetMetadataItem( GDAL_DMD_CREATIONOPTIONLIST, "<CreationOptionList>" " <Option name='DESCRIPTION' type='string' description='Table Description'/>" " <Option name='INSERT' type='string' description='Column Values'/>" " <Option name='BLOCKXSIZE' type='int' description='Column Block Size' " "default='512'/>" " <Option name='BLOCKYSIZE' type='int' description='Row Block Size' " "default='512'/>" " <Option name='BLOCKBSIZE' type='int' description='Band Block Size'/>" " <Option name='BLOCKING' type='string-select' default='YES'>" " <Value>YES</Value>" " <Value>NO</Value>" " <Value>OPTIMALPADDING</Value>" " </Option>" " <Option name='SRID' type='int' description='Overwrite EPSG code'/>" " <Option name='GENPYRAMID' type='string-select' " " description='Generate Pyramid, inform resampling method'>" " <Value>NN</Value>" " <Value>BILINEAR</Value>" " <Value>BIQUADRATIC</Value>" " <Value>CUBIC</Value>" " <Value>AVERAGE4</Value>" " <Value>AVERAGE16</Value>" " </Option>" " <Option name='GENPYRLEVELS' type='int' description='Number of pyramid level to generate'/>" " <Option name='OBJECTTABLE' type='boolean' " "description='Create RDT as object table'/>" " <Option name='SPATIALEXTENT' type='boolean' " "description='Generate Spatial Extent' " "default='TRUE'/>" " <Option name='EXTENTSRID' type='int' description='Spatial ExtentSRID code'/>" " <Option name='COORDLOCATION' type='string-select' default='CENTER'>" " <Value>CENTER</Value>" " <Value>UPPERLEFT</Value>" " </Option>" " <Option name='VATNAME' type='string' description='Value Attribute Table Name'/>" " <Option name='NBITS' type='int' description='BITS for sub-byte " "data types (1,2,4) bits'/>" " <Option name='INTERLEAVE' type='string-select'>" " <Value>BSQ</Value>" " <Value>BIP</Value>" " <Value>BIL</Value>" " </Option>" " <Option name='COMPRESS' type='string-select'>" " <Value>NONE</Value>" " <Value>JPEG-F</Value>" " <Value>JP2-F</Value>" " <Value>DEFLATE</Value>" " </Option>" " <Option name='QUALITY' type='int' description='JPEG quality 0..100' " "default='75'/>" " <Option name='JP2_QUALITY' type='string' description='For JP2-F compression, single quality value or comma separated list " "of increasing quality values for several layers, each in the 0-100 range' default='25'/>" " <Option name='JP2_BLOCKXSIZE' type='int' description='For JP2 compression, tile Width' default='1024'/>" " <Option name='JP2_BLOCKYSIZE' type='int' description='For JP2 compression, tile Height' default='1024'/>" " <Option name='JP2_REVERSIBLE' type='boolean' description='For JP2-F compression, True if the compression is reversible' default='false'/>" " <Option name='JP2_RESOLUTIONS' type='int' description='For JP2-F compression, Number of resolutions.' min='1' max='30'/>" " <Option name='JP2_PROGRESSION' type='string-select' description='For JP2-F compression, progression order' default='LRCP'>" " <Value>LRCP</Value>" " <Value>RLCP</Value>" " <Value>RPCL</Value>" " <Value>PCRL</Value>" " <Value>CPRL</Value>" " </Option>" "</CreationOptionList>" ); poDriver->pfnOpen = GeoRasterDataset::Open; poDriver->pfnCreate = GeoRasterDataset::Create; poDriver->pfnCreateCopy = GeoRasterDataset::CreateCopy; poDriver->pfnIdentify = GeoRasterDataset::Identify; poDriver->pfnDelete = GeoRasterDataset::Delete; GetGDALDriverManager()->RegisterDriver( poDriver ); VSIInstallOCILobHandler(); }