int GDALJP2Box::DumpReadable( FILE *fpOut, int nIndentLevel )
{
if( fpOut == NULL )
fpOut = stdout;
int i;
for(i=0;i<nIndentLevel;i++)
fprintf( fpOut, " " );
fprintf( fpOut, " Type=%s, Offset=" CPL_FRMT_GIB "/" CPL_FRMT_GIB", Data Size=" CPL_FRMT_GIB,
szBoxType, nBoxOffset, nDataOffset,
GetDataLength() );
if( IsSuperBox() )
{
fprintf( fpOut, " (super)" );
}
fprintf( fpOut, "\n" );
if( IsSuperBox() )
{
GDALJP2Box oSubBox( GetFILE() );
for( oSubBox.ReadFirstChild( this );
strlen(oSubBox.GetType()) > 0;
oSubBox.ReadNextChild( this ) )
{
oSubBox.DumpReadable( fpOut, nIndentLevel + 1 );
}
}
if( EQUAL(GetType(),"uuid") )
{
char *pszHex = CPLBinaryToHex( 16, GetUUID() );
for(i=0;i<nIndentLevel;i++)
fprintf( fpOut, " " );
fprintf( fpOut, " UUID=%s", pszHex );
if( EQUAL(pszHex,"B14BF8BD083D4B43A5AE8CD7D5A6CE03") )
fprintf( fpOut, " (GeoTIFF)" );
if( EQUAL(pszHex,"96A9F1F1DC98402DA7AED68E34451809") )
fprintf( fpOut, " (MSI Worldfile)" );
if( EQUAL(pszHex,"BE7ACFCB97A942E89C71999491E3AFAC") )
fprintf( fpOut, " (XMP)" );
CPLFree( pszHex );
fprintf( fpOut, "\n" );
}
return 0;
}
int GDALJP2Box::DumpReadable( FILE *fpOut )
{
if( fpOut == NULL )
fpOut = stdout;
fprintf( fpOut, " Type=%s, Offset=%d/%d, Data Size=%d",
szBoxType, (int) nBoxOffset, (int) nDataOffset,
(int)(nBoxLength - (nDataOffset - nBoxOffset)) );
if( IsSuperBox() )
{
fprintf( fpOut, " (super)" );
}
fprintf( fpOut, "\n" );
if( IsSuperBox() )
{
GDALJP2Box oSubBox( GetFILE() );
for( oSubBox.ReadFirstChild( this );
strlen(oSubBox.GetType()) > 0;
oSubBox.ReadNextChild( this ) )
{
oSubBox.DumpReadable( fpOut );
}
printf( " (end of %s subboxes)\n", szBoxType );
}
if( EQUAL(GetType(),"uuid") )
{
char *pszHex = CPLBinaryToHex( 16, GetUUID() );
fprintf( fpOut, " UUID=%s", pszHex );
if( EQUAL(pszHex,"B14BF8BD083D4B43A5AE8CD7D5A6CE03") )
fprintf( fpOut, " (GeoTIFF)" );
if( EQUAL(pszHex,"96A9F1F1DC98402DA7AED68E34451809") )
fprintf( fpOut, " (MSI Worldfile)" );
CPLFree( pszHex );
fprintf( fpOut, "\n" );
}
return 0;
}
int GDALJP2Metadata::ReadBoxes( VSILFILE *fpVSIL )
{
GDALJP2Box oBox( fpVSIL );
int iBox = 0;
oBox.ReadFirst();
while( strlen(oBox.GetType()) > 0 )
{
/* -------------------------------------------------------------------- */
/* Collect geotiff box. */
/* -------------------------------------------------------------------- */
if( EQUAL(oBox.GetType(),"uuid")
&& memcmp( oBox.GetUUID(), msi_uuid2, 16 ) == 0 )
{
nGeoTIFFSize = (int) oBox.GetDataLength();
pabyGeoTIFFData = oBox.ReadBoxData();
}
/* -------------------------------------------------------------------- */
/* Collect MSIG box. */
/* -------------------------------------------------------------------- */
if( EQUAL(oBox.GetType(),"uuid")
&& memcmp( oBox.GetUUID(), msig_uuid, 16 ) == 0 )
{
nMSIGSize = (int) oBox.GetDataLength();
pabyMSIGData = oBox.ReadBoxData();
if( nMSIGSize < 70
|| memcmp( pabyMSIGData, "MSIG/", 5 ) != 0 )
{
CPLFree( pabyMSIGData );
pabyMSIGData = NULL;
nMSIGSize = 0;
}
}
/* -------------------------------------------------------------------- */
/* Process asoc box looking for Labelled GML data. */
/* -------------------------------------------------------------------- */
if( EQUAL(oBox.GetType(),"asoc") )
{
GDALJP2Box oSubBox( fpVSIL );
oSubBox.ReadFirstChild( &oBox );
if( EQUAL(oSubBox.GetType(),"lbl ") )
{
char *pszLabel = (char *) oSubBox.ReadBoxData();
if( EQUAL(pszLabel,"gml.data") )
{
CollectGMLData( &oBox );
}
CPLFree( pszLabel );
}
}
/* -------------------------------------------------------------------- */
/* Process simple xml boxes. */
/* -------------------------------------------------------------------- */
if( EQUAL(oBox.GetType(),"xml ") )
{
CPLString osBoxName;
char *pszXML = (char *) oBox.ReadBoxData();
osBoxName.Printf( "BOX_%d", iBox++ );
papszGMLMetadata = CSLSetNameValue( papszGMLMetadata,
osBoxName, pszXML );
CPLFree( pszXML );
}
oBox.ReadNext();
}
return TRUE;
}
int GDALJP2Metadata::ReadBoxes( VSILFILE *fpVSIL )
{
GDALJP2Box oBox( fpVSIL );
int iBox = 0;
oBox.ReadFirst();
while( strlen(oBox.GetType()) > 0 )
{
/* -------------------------------------------------------------------- */
/* Collect geotiff box. */
/* -------------------------------------------------------------------- */
if( EQUAL(oBox.GetType(),"uuid")
&& memcmp( oBox.GetUUID(), msi_uuid2, 16 ) == 0 )
{
nGeoTIFFSize = (int) oBox.GetDataLength();
pabyGeoTIFFData = oBox.ReadBoxData();
}
/* -------------------------------------------------------------------- */
/* Collect MSIG box. */
/* -------------------------------------------------------------------- */
if( EQUAL(oBox.GetType(),"uuid")
&& memcmp( oBox.GetUUID(), msig_uuid, 16 ) == 0 )
{
nMSIGSize = (int) oBox.GetDataLength();
pabyMSIGData = oBox.ReadBoxData();
if( nMSIGSize < 70
|| memcmp( pabyMSIGData, "MSIG/", 5 ) != 0 )
{
CPLFree( pabyMSIGData );
pabyMSIGData = NULL;
nMSIGSize = 0;
}
}
/* -------------------------------------------------------------------- */
/* Collect XMP box. */
/* -------------------------------------------------------------------- */
if( EQUAL(oBox.GetType(),"uuid")
&& memcmp( oBox.GetUUID(), xmp_uuid, 16 ) == 0 &&
pszXMPMetadata == NULL )
{
pszXMPMetadata = (char*) oBox.ReadBoxData();
}
/* -------------------------------------------------------------------- */
/* Process asoc box looking for Labelled GML data. */
/* -------------------------------------------------------------------- */
if( EQUAL(oBox.GetType(),"asoc") )
{
GDALJP2Box oSubBox( fpVSIL );
oSubBox.ReadFirstChild( &oBox );
if( EQUAL(oSubBox.GetType(),"lbl ") )
{
char *pszLabel = (char *) oSubBox.ReadBoxData();
if( EQUAL(pszLabel,"gml.data") )
{
CollectGMLData( &oBox );
}
CPLFree( pszLabel );
}
}
/* -------------------------------------------------------------------- */
/* Process simple xml boxes. */
/* -------------------------------------------------------------------- */
if( EQUAL(oBox.GetType(),"xml ") )
{
CPLString osBoxName;
char *pszXML = (char *) oBox.ReadBoxData();
osBoxName.Printf( "BOX_%d", iBox++ );
papszGMLMetadata = CSLSetNameValue( papszGMLMetadata,
osBoxName, pszXML );
CPLFree( pszXML );
}
/* -------------------------------------------------------------------- */
/* Check for a resd box in jp2h. */
/* -------------------------------------------------------------------- */
if( EQUAL(oBox.GetType(),"jp2h") )
{
GDALJP2Box oSubBox( fpVSIL );
for( oSubBox.ReadFirstChild( &oBox );
strlen(oSubBox.GetType()) > 0;
oSubBox.ReadNextChild( &oBox ) )
{
if( EQUAL(oSubBox.GetType(),"res ") )
{
GDALJP2Box oResBox( fpVSIL );
oResBox.ReadFirstChild( &oSubBox );
// we will use either the resd or resc box, which ever
// happens to be first. Should we prefer resd?
if( oResBox.GetDataLength() == 10 )
{
unsigned char *pabyResData = oResBox.ReadBoxData();
int nVertNum, nVertDen, nVertExp;
int nHorzNum, nHorzDen, nHorzExp;
nVertNum = pabyResData[0] * 256 + pabyResData[1];
nVertDen = pabyResData[2] * 256 + pabyResData[3];
nHorzNum = pabyResData[4] * 256 + pabyResData[5];
nHorzDen = pabyResData[6] * 256 + pabyResData[7];
nVertExp = pabyResData[8];
nHorzExp = pabyResData[9];
// compute in pixels/cm
double dfVertRes =
(nVertNum/(double)nVertDen) * pow(10.0,nVertExp)/100;
double dfHorzRes =
(nHorzNum/(double)nHorzDen) * pow(10.0,nHorzExp)/100;
CPLString osFormatter;
papszMetadata = CSLSetNameValue(
papszMetadata,
"TIFFTAG_XRESOLUTION",
osFormatter.Printf("%g",dfHorzRes) );
papszMetadata = CSLSetNameValue(
papszMetadata,
"TIFFTAG_YRESOLUTION",
osFormatter.Printf("%g",dfVertRes) );
papszMetadata = CSLSetNameValue(
papszMetadata,
"TIFFTAG_RESOLUTIONUNIT",
"3 (pixels/cm)" );
CPLFree( pabyResData );
}
}
}
}
oBox.ReadNext();
}
return TRUE;
}
