int SDTSRawPolygon::Read( DDFRecord * poRecord ) { /* ==================================================================== */ /* Loop over fields in this record, looking for those we */ /* recognise, and need. */ /* ==================================================================== */ for( int iField = 0; iField < poRecord->GetFieldCount(); iField++ ) { DDFField *poField = poRecord->GetField( iField ); if( poField == nullptr ) return FALSE; DDFFieldDefn* poFieldDefn = poField->GetFieldDefn(); if( poFieldDefn == nullptr ) return FALSE; const char *pszFieldName = poFieldDefn->GetName(); if( EQUAL(pszFieldName,"POLY") ) { oModId.Set( poField ); } else if( EQUAL(pszFieldName,"ATID") ) { ApplyATID( poField ); } } return TRUE; }
int SDTSModId::Set( DDFField *poField ) { const char *pachData = poField->GetData(); DDFFieldDefn *poDefn = poField->GetFieldDefn(); if( poDefn->GetSubfieldCount() >= 2 && poDefn->GetSubfield(0)->GetWidth() == 4 ) { if( strlen(pachData) < 5 ) return FALSE; memcpy( szModule, pachData, 4 ); szModule[4] = '\0'; nRecord = atoi( pachData + 4 ); } else { DDFSubfieldDefn *poSF = poField->GetFieldDefn()->FindSubfieldDefn( "MODN" ); if( poSF == nullptr ) return FALSE; int nBytesRemaining; pachData = poField->GetSubfieldData(poSF, &nBytesRemaining); if( pachData == nullptr ) return FALSE; snprintf( szModule, sizeof(szModule), "%s", poSF->ExtractStringData( pachData, nBytesRemaining, nullptr) ); poSF = poField->GetFieldDefn()->FindSubfieldDefn( "RCID" ); if( poSF != nullptr ) { pachData = poField->GetSubfieldData(poSF, &nBytesRemaining); if( pachData != nullptr ) nRecord = poSF->ExtractIntData( pachData, nBytesRemaining, nullptr); } } if( poDefn->GetSubfieldCount() == 3 ) { DDFSubfieldDefn *poSF = poField->GetFieldDefn()->FindSubfieldDefn( "OBRP" ); if( poSF != nullptr ) { int nBytesRemaining; pachData = poField->GetSubfieldData(poSF, &nBytesRemaining); if( pachData != nullptr ) { snprintf( szOBRP, sizeof(szOBRP), "%s", poSF->ExtractStringData( pachData, nBytesRemaining, nullptr) ); } } } return FALSE; }
int SDTSModId::Set( DDFField *poField ) { const char *pachData = poField->GetData(); DDFFieldDefn *poDefn = poField->GetFieldDefn(); if( poDefn->GetSubfieldCount() >= 2 && poDefn->GetSubfield(0)->GetWidth() == 4 ) { memcpy( szModule, pachData, 4 ); szModule[4] = '\0'; nRecord = atoi( pachData + 4 ); } else { DDFSubfieldDefn *poSF; int nBytesRemaining; const char *pachData; poSF = poField->GetFieldDefn()->FindSubfieldDefn( "MODN" ); pachData = poField->GetSubfieldData(poSF, &nBytesRemaining); strncpy( szModule, poSF->ExtractStringData( pachData, nBytesRemaining, NULL), sizeof(szModule) ); szModule[sizeof(szModule)-1] = '\0'; poSF = poField->GetFieldDefn()->FindSubfieldDefn( "RCID" ); pachData = poField->GetSubfieldData(poSF, &nBytesRemaining); nRecord = poSF->ExtractIntData( pachData, nBytesRemaining, NULL); } if( poDefn->GetSubfieldCount() == 3 ) { DDFSubfieldDefn *poSF; poSF = poField->GetFieldDefn()->FindSubfieldDefn( "OBRP" ); if( poSF != NULL ) { int nBytesRemaining; const char *pachData; pachData = poField->GetSubfieldData(poSF, &nBytesRemaining); strncpy( szOBRP, poSF->ExtractStringData( pachData, nBytesRemaining, NULL), sizeof(szOBRP) ); szOBRP[sizeof(szOBRP)-1] = '\0'; } } return FALSE; }
static void ViewRecordField( DDFField * poField ) { int nBytesRemaining; const char *pachFieldData; DDFFieldDefn *poFieldDefn = poField->GetFieldDefn(); // Report general information about the field. printf( " Field %s: %s\n", poFieldDefn->GetName(), poFieldDefn->GetDescription() ); // Get pointer to this fields raw data. We will move through // it consuming data as we report subfield values. pachFieldData = poField->GetData(); nBytesRemaining = poField->GetDataSize(); /* -------------------------------------------------------- */ /* Loop over the repeat count for this fields */ /* subfields. The repeat count will almost */ /* always be one. */ /* -------------------------------------------------------- */ int iRepeat; for( iRepeat = 0; iRepeat < poField->GetRepeatCount(); iRepeat++ ) { /* -------------------------------------------------------- */ /* Loop over all the subfields of this field, advancing */ /* the data pointer as we consume data. */ /* -------------------------------------------------------- */ int iSF; for( iSF = 0; iSF < poFieldDefn->GetSubfieldCount(); iSF++ ) { DDFSubfieldDefn *poSFDefn = poFieldDefn->GetSubfield( iSF ); int nBytesConsumed; nBytesConsumed = ViewSubfield( poSFDefn, pachFieldData, nBytesRemaining ); nBytesRemaining -= nBytesConsumed; pachFieldData += nBytesConsumed; } } }
int SDTSRawPoint::Read( SDTS_IREF * poIREF, DDFRecord * poRecord ) { /* ==================================================================== */ /* Loop over fields in this record, looking for those we */ /* recognise, and need. */ /* ==================================================================== */ for( int iField = 0; iField < poRecord->GetFieldCount(); iField++ ) { DDFField *poField = poRecord->GetField( iField ); if( poField == nullptr ) return FALSE; DDFFieldDefn* poFieldDefn = poField->GetFieldDefn(); if( poFieldDefn == nullptr ) return FALSE; const char *pszFieldName = poFieldDefn->GetName(); if( EQUAL(pszFieldName,"PNTS") ) oModId.Set( poField ); else if( EQUAL(pszFieldName,"ATID") ) ApplyATID( poField ); else if( EQUAL(pszFieldName,"ARID") ) { oAreaId.Set( poField ); } else if( EQUAL(pszFieldName,"SADR") ) { poIREF->GetSADR( poField, 1, &dfX, &dfY, &dfZ ); } } return TRUE; }
int main( int nArgc, char ** papszArgv ) { DDFModule oModule; const char *pszFilename = NULL; int bFSPTHack = FALSE; int bXML = FALSE; /* -------------------------------------------------------------------- */ /* Check arguments. */ /* -------------------------------------------------------------------- */ for( int iArg = 1; iArg < nArgc; iArg++ ) { if( EQUAL(papszArgv[iArg],"-fspt_repeating") ) bFSPTHack = TRUE; else if( EQUAL(papszArgv[iArg],"-xml") ) bXML = TRUE; else pszFilename = papszArgv[iArg]; } if( pszFilename == NULL ) { printf( "Usage: 8211dump [-xml] [-fspt_repeating] filename\n" ); exit( 1 ); } /* -------------------------------------------------------------------- */ /* Open file. */ /* -------------------------------------------------------------------- */ if( !oModule.Open( pszFilename ) ) exit( 1 ); /* -------------------------------------------------------------------- */ /* Apply FSPT hack if required. */ /* -------------------------------------------------------------------- */ if( bFSPTHack ) { DDFFieldDefn *poFSPT = oModule.FindFieldDefn( "FSPT" ); if( poFSPT == NULL ) fprintf( stderr, "unable to find FSPT field to set repeating flag.\n" ); else poFSPT->SetRepeatingFlag( TRUE ); } /* -------------------------------------------------------------------- */ /* Dump header, and all records. */ /* -------------------------------------------------------------------- */ DDFRecord *poRecord; if( bXML ) { printf("<DDFModule>\n"); int nFieldDefnCount = oModule.GetFieldCount(); for( int i = 0; i < nFieldDefnCount; i++ ) { DDFFieldDefn* poFieldDefn = oModule.GetField(i); const char* pszDataStructCode; switch( poFieldDefn->GetDataStructCode() ) { case dsc_elementary: pszDataStructCode = "elementary"; break; case dsc_vector: pszDataStructCode = "vector"; break; case dsc_array: pszDataStructCode = "array"; break; case dsc_concatenated: pszDataStructCode = "concatenated"; break; default: pszDataStructCode = "(unknown)"; break; } const char* pszDataTypeCode; switch( poFieldDefn->GetDataTypeCode() ) { case dtc_char_string: pszDataTypeCode = "char_string"; break; case dtc_implicit_point: pszDataTypeCode = "implicit_point"; break; case dtc_explicit_point: pszDataTypeCode = "explicit_point"; break; case dtc_explicit_point_scaled: pszDataTypeCode = "explicit_point_scaled"; break; case dtc_char_bit_string: pszDataTypeCode = "char_bit_string"; break; case dtc_bit_string: pszDataTypeCode = "bit_string"; break; case dtc_mixed_data_type: pszDataTypeCode = "mixed_data_type"; break; default: pszDataTypeCode = "(unknown)"; break; } printf("<DDFFieldDefn tag=\"%s\" fieldName=\"%s\" arrayDescr=\"%s\" " "formatControls=\"%s\" dataStructCode=\"%s\" dataTypeCode=\"%s\">\n", poFieldDefn->GetName(), poFieldDefn->GetDescription(), poFieldDefn->GetArrayDescr(), poFieldDefn->GetFormatControls(), pszDataStructCode, pszDataTypeCode); int nSubfieldCount = poFieldDefn->GetSubfieldCount(); for( int iSubField = 0; iSubField < nSubfieldCount; iSubField++ ) { DDFSubfieldDefn* poSubFieldDefn = poFieldDefn->GetSubfield(iSubField); printf(" <DDFSubfieldDefn name=\"%s\" format=\"%s\"/>\n", poSubFieldDefn->GetName(), poSubFieldDefn->GetFormat()); } printf("</DDFFieldDefn>\n"); } for( poRecord = oModule.ReadRecord(); poRecord != NULL; poRecord = oModule.ReadRecord() ) { printf("<DDFRecord>\n"); int nFieldCount = poRecord->GetFieldCount(); for( int iField = 0; iField < nFieldCount; iField++ ) { DDFField* poField = poRecord->GetField(iField); DDFFieldDefn* poDefn = poField->GetFieldDefn(); const char* pszFieldName = poDefn->GetName(); printf(" <DDFField name=\"%s\"", pszFieldName); if( poField->GetRepeatCount() > 1 ) printf(" repeatCount=\"%d\"", poField->GetRepeatCount()); int iOffset = 0, nLoopCount; int nRepeatCount = poField->GetRepeatCount(); const char* pachData = poField->GetData(); int nDataSize = poField->GetDataSize(); if( nRepeatCount == 1 && poDefn->GetSubfieldCount() == 0 ) { printf(" value=\"0x"); for( int i = 0; i < nDataSize - 1; i++ ) printf( "%02X", pachData[i] ); printf("\">\n"); } else printf(">\n"); for( nLoopCount = 0; nLoopCount < nRepeatCount; nLoopCount++ ) { for( int iSubField = 0; iSubField < poDefn->GetSubfieldCount(); iSubField++ ) { int nBytesConsumed; DDFSubfieldDefn* poSubFieldDefn = poDefn->GetSubfield(iSubField); const char* pszSubFieldName = poSubFieldDefn->GetName(); printf(" <DDFSubfield name=\"%s\" ", pszSubFieldName); DDFDataType eType = poSubFieldDefn->GetType(); const char* pachSubdata = pachData + iOffset; int nMaxBytes = nDataSize - iOffset; if( eType == DDFFloat ) printf("type=\"float\">%f", poSubFieldDefn->ExtractFloatData( pachSubdata, nMaxBytes, NULL ) ); else if( eType == DDFInt ) printf("type=\"integer\">%d", poSubFieldDefn->ExtractIntData( pachSubdata, nMaxBytes, NULL ) ); else if( eType == DDFBinaryString ) { int nBytes, i; GByte *pabyBString = (GByte *) poSubFieldDefn->ExtractStringData( pachSubdata, nMaxBytes, &nBytes ); printf( "type=\"binary\">0x" ); for( i = 0; i < nBytes; i++ ) printf( "%02X", pabyBString[i] ); } else { GByte* pabyString = (GByte *)poSubFieldDefn->ExtractStringData( pachSubdata, nMaxBytes, NULL ); int bBinary = FALSE; int i; for( i = 0; pabyString[i] != '\0'; i ++ ) { if( pabyString[i] < 32 || pabyString[i] > 127 ) { bBinary = TRUE; break; } } if( bBinary ) { printf( "type=\"binary\">0x" ); for( i = 0; pabyString[i] != '\0'; i ++ ) printf( "%02X", pabyString[i] ); } else { char* pszEscaped = CPLEscapeString((const char*)pabyString, -1, CPLES_XML); printf("type=\"string\">%s", pszEscaped); CPLFree(pszEscaped); } } printf("</DDFSubfield>\n"); poSubFieldDefn->GetDataLength( pachSubdata, nMaxBytes, &nBytesConsumed ); iOffset += nBytesConsumed; } } printf(" </DDFField>\n"); } printf("</DDFRecord>\n"); } printf("</DDFModule>\n"); } else { oModule.Dump( stdout ); long nStartLoc; nStartLoc = VSIFTellL( oModule.GetFP() ); for( poRecord = oModule.ReadRecord(); poRecord != NULL; poRecord = oModule.ReadRecord() ) { printf( "File Offset: %ld\n", nStartLoc ); poRecord->Dump( stdout ); nStartLoc = VSIFTellL( oModule.GetFP() ); } } oModule.Close(); #ifdef DBMALLOC malloc_dump(1); #endif }
int DDFModule::Open( const char * pszFilename, int bFailQuietly ) { static const int nLeaderSize = 24; /* -------------------------------------------------------------------- */ /* Close the existing file if there is one. */ /* -------------------------------------------------------------------- */ if( fpDDF != NULL ) Close(); /* -------------------------------------------------------------------- */ /* Open the file. */ /* -------------------------------------------------------------------- */ fpDDF = VSIFOpenL( pszFilename, "rb" ); if( fpDDF == NULL ) { if( !bFailQuietly ) CPLError( CE_Failure, CPLE_OpenFailed, "Unable to open DDF file `%s'.", pszFilename ); return FALSE; } /* -------------------------------------------------------------------- */ /* Read the 24 byte leader. */ /* -------------------------------------------------------------------- */ char achLeader[nLeaderSize]; if( (int)VSIFReadL( achLeader, 1, nLeaderSize, fpDDF ) != nLeaderSize ) { CPL_IGNORE_RET_VAL(VSIFCloseL( fpDDF )); fpDDF = NULL; if( !bFailQuietly ) CPLError( CE_Failure, CPLE_FileIO, "Leader is short on DDF file `%s'.", pszFilename ); return FALSE; } /* -------------------------------------------------------------------- */ /* Verify that this appears to be a valid DDF file. */ /* -------------------------------------------------------------------- */ int i, bValid = TRUE; for( i = 0; i < nLeaderSize; i++ ) { if( achLeader[i] < 32 || achLeader[i] > 126 ) bValid = FALSE; } if( achLeader[5] != '1' && achLeader[5] != '2' && achLeader[5] != '3' ) bValid = FALSE; if( achLeader[6] != 'L' ) bValid = FALSE; if( achLeader[8] != '1' && achLeader[8] != ' ' ) bValid = FALSE; /* -------------------------------------------------------------------- */ /* Extract information from leader. */ /* -------------------------------------------------------------------- */ if( bValid ) { _recLength = DDFScanInt( achLeader+0, 5 ); _interchangeLevel = achLeader[5]; _leaderIden = achLeader[6]; _inlineCodeExtensionIndicator = achLeader[7]; _versionNumber = achLeader[8]; _appIndicator = achLeader[9]; _fieldControlLength = DDFScanInt(achLeader+10,2); _fieldAreaStart = DDFScanInt(achLeader+12,5); _extendedCharSet[0] = achLeader[17]; _extendedCharSet[1] = achLeader[18]; _extendedCharSet[2] = achLeader[19]; _extendedCharSet[3] = '\0'; _sizeFieldLength = DDFScanInt(achLeader+20,1); _sizeFieldPos = DDFScanInt(achLeader+21,1); _sizeFieldTag = DDFScanInt(achLeader+23,1); if( _recLength < nLeaderSize || _fieldControlLength == 0 || _fieldAreaStart < 24 || _sizeFieldLength == 0 || _sizeFieldPos == 0 || _sizeFieldTag == 0 ) { bValid = FALSE; } } /* -------------------------------------------------------------------- */ /* If the header is invalid, then clean up, report the error */ /* and return. */ /* -------------------------------------------------------------------- */ if( !bValid ) { CPL_IGNORE_RET_VAL(VSIFCloseL( fpDDF )); fpDDF = NULL; if( !bFailQuietly ) CPLError( CE_Failure, CPLE_AppDefined, "File `%s' does not appear to have\n" "a valid ISO 8211 header.\n", pszFilename ); return FALSE; } /* -------------------------------------------------------------------- */ /* Read the whole record info memory. */ /* -------------------------------------------------------------------- */ char *pachRecord = (char *) CPLMalloc(_recLength); memcpy( pachRecord, achLeader, nLeaderSize ); if( (int)VSIFReadL( pachRecord+nLeaderSize, 1, _recLength-nLeaderSize, fpDDF ) != _recLength - nLeaderSize ) { if( !bFailQuietly ) CPLError( CE_Failure, CPLE_FileIO, "Header record is short on DDF file `%s'.", pszFilename ); return FALSE; } /* -------------------------------------------------------------------- */ /* First make a pass counting the directory entries. */ /* -------------------------------------------------------------------- */ int nFieldEntryWidth, nFDCount = 0; nFieldEntryWidth = _sizeFieldLength + _sizeFieldPos + _sizeFieldTag; for( i = nLeaderSize; i < _recLength; i += nFieldEntryWidth ) { if( pachRecord[i] == DDF_FIELD_TERMINATOR ) break; nFDCount++; } /* -------------------------------------------------------------------- */ /* Allocate, and read field definitions. */ /* -------------------------------------------------------------------- */ for( i = 0; i < nFDCount; i++ ) { char szTag[128]; int nEntryOffset = nLeaderSize + i*nFieldEntryWidth; int nFieldLength, nFieldPos; strncpy( szTag, pachRecord+nEntryOffset, _sizeFieldTag ); szTag[_sizeFieldTag] = '\0'; nEntryOffset += _sizeFieldTag; nFieldLength = DDFScanInt( pachRecord+nEntryOffset, _sizeFieldLength ); nEntryOffset += _sizeFieldLength; nFieldPos = DDFScanInt( pachRecord+nEntryOffset, _sizeFieldPos ); if (_fieldAreaStart+nFieldPos < 0 || _recLength - (_fieldAreaStart+nFieldPos) < nFieldLength) { if( !bFailQuietly ) CPLError( CE_Failure, CPLE_FileIO, "Header record invalid on DDF file `%s'.", pszFilename ); CPLFree( pachRecord ); return FALSE; } DDFFieldDefn *poFDefn = new DDFFieldDefn(); if( poFDefn->Initialize( this, szTag, nFieldLength, pachRecord+_fieldAreaStart+nFieldPos ) ) AddField( poFDefn ); else delete poFDefn; } CPLFree( pachRecord ); /* -------------------------------------------------------------------- */ /* Record the current file offset, the beginning of the first */ /* data record. */ /* -------------------------------------------------------------------- */ nFirstRecordOffset = (long)VSIFTellL( fpDDF ); return TRUE; }
void mk_s57() { DDFModule oModule; DDFFieldDefn *poFDefn; oModule.Initialize(); /* -------------------------------------------------------------------- */ /* Create the '0000' definition. */ /* -------------------------------------------------------------------- */ poFDefn = new DDFFieldDefn(); poFDefn->Create( "0000", "", "0001DSIDDSIDDSSI0001DSPM0001VRIDVRIDATTVVRIDVRPCVRIDVRPTVRIDSGCCVRIDSG2DVRIDSG3D0001FRIDFRIDFOIDFRIDATTFFRIDNATFFRIDFFPCFRIDFFPTFRIDFSPCFRIDFSPT", dsc_elementary, dtc_char_string ); oModule.AddField( poFDefn ); /* -------------------------------------------------------------------- */ /* Create the '0001' definition. */ /* -------------------------------------------------------------------- */ poFDefn = new DDFFieldDefn(); poFDefn->Create( "0001", "ISO 8211 Record Identifier", "", dsc_elementary, dtc_bit_string, "(b12)" ); oModule.AddField( poFDefn ); /* -------------------------------------------------------------------- */ /* Create the DSID field. */ /* -------------------------------------------------------------------- */ poFDefn = new DDFFieldDefn(); poFDefn->Create( "DSID", "Data set identification field", "", dsc_vector, dtc_mixed_data_type ); poFDefn->AddSubfield( "RCNM", "b11" ); poFDefn->AddSubfield( "RCID", "b14" ); poFDefn->AddSubfield( "EXPP", "b11" ); poFDefn->AddSubfield( "INTU", "b11" ); poFDefn->AddSubfield( "DSNM", "A" ); poFDefn->AddSubfield( "EDTN", "A" ); poFDefn->AddSubfield( "UPDN", "A" ); poFDefn->AddSubfield( "UADT", "A(8)" ); poFDefn->AddSubfield( "ISDT", "A(8)" ); poFDefn->AddSubfield( "STED", "R(4)" ); poFDefn->AddSubfield( "PRSP", "b11" ); poFDefn->AddSubfield( "PSDN", "A" ); poFDefn->AddSubfield( "PRED", "A" ); poFDefn->AddSubfield( "PROF", "b11" ); poFDefn->AddSubfield( "AGEN", "b12" ); poFDefn->AddSubfield( "COMT", "A" ); oModule.AddField( poFDefn ); /* -------------------------------------------------------------------- */ /* Create the DSSI field. */ /* -------------------------------------------------------------------- */ poFDefn = new DDFFieldDefn(); poFDefn->Create( "DSSI", "Data set structure information field", "", dsc_vector, dtc_mixed_data_type ); poFDefn->AddSubfield( "DSTR", "b11" ); poFDefn->AddSubfield( "AALL", "b11" ); poFDefn->AddSubfield( "NALL", "b11" ); poFDefn->AddSubfield( "NOMR", "b14" ); poFDefn->AddSubfield( "NOCR", "b14" ); poFDefn->AddSubfield( "NOGR", "b14" ); poFDefn->AddSubfield( "NOLR", "b14" ); poFDefn->AddSubfield( "NOIN", "b14" ); poFDefn->AddSubfield( "NOCN", "b14" ); poFDefn->AddSubfield( "NOED", "b14" ); poFDefn->AddSubfield( "NOFA", "b14" ); oModule.AddField( poFDefn ); /* -------------------------------------------------------------------- */ /* Create the DSPM field. */ /* -------------------------------------------------------------------- */ poFDefn = new DDFFieldDefn(); poFDefn->Create( "DSPM", "Data set parameter field", "", dsc_vector, dtc_mixed_data_type ); poFDefn->AddSubfield( "RCNM", "b11" ); poFDefn->AddSubfield( "RCID", "b14" ); poFDefn->AddSubfield( "HDAT", "b11" ); poFDefn->AddSubfield( "VDAT", "b11" ); poFDefn->AddSubfield( "SDAT", "b11" ); poFDefn->AddSubfield( "CSCL", "b14" ); poFDefn->AddSubfield( "DUNI", "b11" ); poFDefn->AddSubfield( "HUNI", "b11" ); poFDefn->AddSubfield( "PUNI", "b11" ); poFDefn->AddSubfield( "COUN", "b11" ); poFDefn->AddSubfield( "COMF", "b14" ); poFDefn->AddSubfield( "SOMF", "b14" ); poFDefn->AddSubfield( "COMT", "A" ); oModule.AddField( poFDefn ); /* -------------------------------------------------------------------- */ /* Create the VRID field. */ /* -------------------------------------------------------------------- */ poFDefn = new DDFFieldDefn(); poFDefn->Create( "VRID", "Vector record identifier field", "", dsc_vector, dtc_mixed_data_type ); poFDefn->AddSubfield( "RCNM", "b11" ); poFDefn->AddSubfield( "RCID", "b14" ); poFDefn->AddSubfield( "RVER", "b12" ); poFDefn->AddSubfield( "RUIN", "b11" ); oModule.AddField( poFDefn ); /* -------------------------------------------------------------------- */ /* Create the ATTV field. */ /* -------------------------------------------------------------------- */ poFDefn = new DDFFieldDefn(); poFDefn->Create( "ATTV", "Vector record attribute field", "", dsc_vector, dtc_mixed_data_type ); /* how do I mark this as repeating? */ poFDefn->AddSubfield( "ATTL", "b12" ); poFDefn->AddSubfield( "ATVL", "A" ); oModule.AddField( poFDefn ); /* -------------------------------------------------------------------- */ /* Create the SG2D field. */ /* -------------------------------------------------------------------- */ poFDefn = new DDFFieldDefn(); poFDefn->Create( "SG2D", "2-D coordinate field", "*", dsc_vector, dtc_mixed_data_type ); /* how do I mark this as repeating? */ poFDefn->AddSubfield( "YCOO", "b24" ); poFDefn->AddSubfield( "XCOO", "b24" ); oModule.AddField( poFDefn ); /* -------------------------------------------------------------------- */ /* Create the SG3D field. */ /* -------------------------------------------------------------------- */ poFDefn = new DDFFieldDefn(); poFDefn->Create( "SG3D", "3-D coordinate (sounding array) field", "*", dsc_vector, dtc_mixed_data_type ); /* how do I mark this as repeating? */ poFDefn->AddSubfield( "YCOO", "b24" ); poFDefn->AddSubfield( "XCOO", "b24" ); poFDefn->AddSubfield( "VE3D", "b24" ); oModule.AddField( poFDefn ); /* -------------------------------------------------------------------- */ // need to add: VRPC, VRPT, SGCC, FRID, FOID, ATTF, NATF, FFPC, // FFPT, FSPC, and FSPT /* -------------------------------------------------------------------- */ /* Create file. */ /* -------------------------------------------------------------------- */ oModule.Dump( stdout ); oModule.Create( "out.000" ); /* -------------------------------------------------------------------- */ /* Create a record. */ /* -------------------------------------------------------------------- */ DDFRecord *poRec = new DDFRecord( &oModule ); DDFField *poField; poField = poRec->AddField( oModule.FindFieldDefn( "0001" ) ); poRec->SetFieldRaw( poField, 0, "\1\0\036", 3 ); poField = poRec->AddField( oModule.FindFieldDefn( "DSID" ) ); poRec->SetIntSubfield ( "DSID", 0, "RCNM", 0, 10 ); poRec->SetIntSubfield ( "DSID", 0, "RCID", 0, 1 ); poRec->SetIntSubfield ( "DSID", 0, "EXPP", 0, 1 ); poRec->SetIntSubfield ( "DSID", 0, "INTU", 0, 4 ); poRec->SetStringSubfield( "DSID", 0, "DSNM", 0, "GB4X0000.000" ); poRec->SetStringSubfield( "DSID", 0, "EDTN", 0, "2" ); poRec->SetStringSubfield( "DSID", 0, "UPDN", 0, "0" ); poRec->SetStringSubfield( "DSID", 0, "UADT", 0, "20010409" ); poRec->SetStringSubfield( "DSID", 0, "ISDT", 0, "20010409" ); poRec->SetFloatSubfield ( "DSID", 0, "STED", 0, 3.1 ); poRec->SetIntSubfield ( "DSID", 0, "PRSP", 0, 1 ); poRec->SetStringSubfield( "DSID", 0, "PSDN", 0, "" ); poRec->SetStringSubfield( "DSID", 0, "PRED", 0, "2.0" ); poRec->SetIntSubfield ( "DSID", 0, "PROF", 0, 1 ); poRec->SetIntSubfield ( "DSID", 0, "AGEN", 0, 540 ); poRec->SetStringSubfield( "DSID", 0, "COMT", 0, "" ); poField = poRec->AddField( oModule.FindFieldDefn( "DSSI" ) ); poRec->SetIntSubfield ( "DSSI", 0, "DSTR", 0, 2 ); poRec->SetIntSubfield ( "DSSI", 0, "AALL", 0, 1 ); poRec->SetIntSubfield ( "DSSI", 0, "NALL", 0, 1 ); poRec->SetIntSubfield ( "DSSI", 0, "NOMR", 0, 22 ); poRec->SetIntSubfield ( "DSSI", 0, "NOCR", 0, 0 ); poRec->SetIntSubfield ( "DSSI", 0, "NOGR", 0, 2141 ); poRec->SetIntSubfield ( "DSSI", 0, "NOLR", 0, 15 ); poRec->SetIntSubfield ( "DSSI", 0, "NOIN", 0, 512 ); poRec->SetIntSubfield ( "DSSI", 0, "NOCN", 0, 2181 ); poRec->SetIntSubfield ( "DSSI", 0, "NOED", 0, 3192 ); poRec->SetIntSubfield ( "DSSI", 0, "NOFA", 0, 0 ); poRec->Write(); delete poRec; /* -------------------------------------------------------------------- */ /* Create a record. */ /* -------------------------------------------------------------------- */ poRec = new DDFRecord( &oModule ); poField = poRec->AddField( oModule.FindFieldDefn( "0001" ) ); poRec->SetFieldRaw( poField, 0, "\2\0\036", 3 ); poField = poRec->AddField( oModule.FindFieldDefn( "DSPM" ) ); poRec->SetIntSubfield ( "DSPM", 0, "RCNM", 0, 20 ); poRec->SetIntSubfield ( "DSPM", 0, "RCID", 0, 1 ); poRec->SetIntSubfield ( "DSPM", 0, "HDAT", 0, 2 ); poRec->SetIntSubfield ( "DSPM", 0, "VDAT", 0, 17 ); poRec->SetIntSubfield ( "DSPM", 0, "SDAT", 0, 23 ); poRec->SetIntSubfield ( "DSPM", 0, "CSCL", 0, 52000 ); poRec->SetIntSubfield ( "DSPM", 0, "DUNI", 0, 1 ); poRec->SetIntSubfield ( "DSPM", 0, "HUNI", 0, 1 ); poRec->SetIntSubfield ( "DSPM", 0, "PUNI", 0, 1 ); poRec->SetIntSubfield ( "DSPM", 0, "COUN", 0, 1 ); poRec->SetIntSubfield ( "DSPM", 0, "COMF", 0, 1000000 ); poRec->SetIntSubfield ( "DSPM", 0, "SOMF", 0, 10 ); poRec->Write(); delete poRec; /* -------------------------------------------------------------------- */ /* Create a record. */ /* -------------------------------------------------------------------- */ poRec = new DDFRecord( &oModule ); poField = poRec->AddField( oModule.FindFieldDefn( "0001" ) ); poRec->SetFieldRaw( poField, 0, "\3\0\036", 3 ); poField = poRec->AddField( oModule.FindFieldDefn( "VRID" ) ); poRec->SetIntSubfield ( "VRID", 0, "RCNM", 0, 110 ); poRec->SetIntSubfield ( "VRID", 0, "RCID", 0, 518 ); poRec->SetIntSubfield ( "VRID", 0, "RVER", 0, 1 ); poRec->SetIntSubfield ( "VRID", 0, "RUIN", 0, 1 ); poField = poRec->AddField( oModule.FindFieldDefn( "SG3D" ) ); poRec->SetIntSubfield ( "SG3D", 0, "YCOO", 0, -325998702 ); poRec->SetIntSubfield ( "SG3D", 0, "XCOO", 0, 612175350 ); poRec->SetIntSubfield ( "SG3D", 0, "VE3D", 0, 174 ); poRec->SetIntSubfield ( "SG3D", 0, "YCOO", 1, -325995189 ); poRec->SetIntSubfield ( "SG3D", 0, "XCOO", 1, 612228812 ); poRec->SetIntSubfield ( "SG3D", 0, "VE3D", 1, 400 ); poRec->Write(); delete poRec; }
void mk_catalog() { DDFModule oModule; DDFFieldDefn *poFDefn; oModule.Initialize(); /* -------------------------------------------------------------------- */ /* Create the '0000' definition. */ /* -------------------------------------------------------------------- */ poFDefn = new DDFFieldDefn(); poFDefn->Create( "0000", "", "0001CATD", dsc_elementary, dtc_char_string ); oModule.AddField( poFDefn ); /* -------------------------------------------------------------------- */ /* Create the '0000' definition. */ /* -------------------------------------------------------------------- */ poFDefn = new DDFFieldDefn(); poFDefn->Create( "0001", "ISO 8211 Record Identifier", "", dsc_elementary, dtc_bit_string, "(b12)" ); oModule.AddField( poFDefn ); /* -------------------------------------------------------------------- */ /* Create the CATD field. */ /* -------------------------------------------------------------------- */ poFDefn = new DDFFieldDefn(); poFDefn->Create( "CATD", "Catalog Directory field", "", dsc_vector, dtc_mixed_data_type ); poFDefn->AddSubfield( "RCNM", "A(2)" ); poFDefn->AddSubfield( "RCID", "I(10)" ); poFDefn->AddSubfield( "FILE", "A" ); poFDefn->AddSubfield( "LFIL", "A" ); poFDefn->AddSubfield( "VOLM", "A" ); poFDefn->AddSubfield( "IMPL", "A(3)" ); poFDefn->AddSubfield( "SLAT", "R" ); poFDefn->AddSubfield( "WLON", "R" ); poFDefn->AddSubfield( "NLAT", "R" ); poFDefn->AddSubfield( "ELON", "R" ); poFDefn->AddSubfield( "CRCS", "A" ); poFDefn->AddSubfield( "COMT", "A" ); oModule.AddField( poFDefn ); oModule.Dump( stdout ); oModule.Create( "out.ddf" ); /* -------------------------------------------------------------------- */ /* Create a record. */ /* -------------------------------------------------------------------- */ DDFRecord *poRec = new DDFRecord( &oModule ); DDFField *poField; poField = poRec->AddField( oModule.FindFieldDefn( "0001" ) ); poRec->SetFieldRaw( poField, 0, "\0\0\036", 3 ); poField = poRec->AddField( oModule.FindFieldDefn( "CATD" ) ); poRec->SetStringSubfield( "CATD", 0, "RCNM", 0, "CD" ); poRec->SetIntSubfield ( "CATD", 0, "RCID", 0, 1 ); poRec->SetStringSubfield( "CATD", 0, "FILE", 0, "CATALOG.030" ); poRec->SetStringSubfield( "CATD", 0, "VOLM", 0, "V01X01" ); poRec->SetStringSubfield( "CATD", 0, "IMPL", 0, "ASC" ); poRec->SetStringSubfield( "CATD", 0, "COMT", 0, "Exchange Set Catalog file" ); poRec->Write(); delete poRec; /* -------------------------------------------------------------------- */ /* Create a record. */ /* -------------------------------------------------------------------- */ poRec = new DDFRecord( &oModule ); poField = poRec->AddField( oModule.FindFieldDefn( "0001" ) ); poRec->SetFieldRaw( poField, 0, "\1\0\036", 3 ); poField = poRec->AddField( oModule.FindFieldDefn( "CATD" ) ); poRec->SetStringSubfield( "CATD", 0, "RCNM", 0, "CD" ); poRec->SetIntSubfield ( "CATD", 0, "RCID", 0, 2 ); poRec->SetStringSubfield( "CATD", 0, "FILE", 0, "No410810.000" ); poRec->SetStringSubfield( "CATD", 0, "VOLM", 0, "V01X01" ); poRec->SetStringSubfield( "CATD", 0, "IMPL", 0, "BIN" ); poRec->SetFloatSubfield ( "CATD", 0, "SLAT", 0, 59.000005 ); poRec->SetFloatSubfield ( "CATD", 0, "WLON", 0, 4.999996 ); poRec->SetFloatSubfield ( "CATD", 0, "NLAT", 0, 59.500003 ); poRec->SetFloatSubfield ( "CATD", 0, "ELON", 0, 5.499997 ); poRec->SetStringSubfield( "CATD", 0, "CRCS", 0, "555C3AD4" ); poRec->Write(); delete poRec; }
static void WriteAttrRecordToDBF( DBFHandle hDBF, int iRecord, SDTSTransfer * poTransfer, DDFField * poSR ) { /* -------------------------------------------------------------------- */ /* Process each subfield in the record. */ /* -------------------------------------------------------------------- */ DDFFieldDefn *poFDefn = poSR->GetFieldDefn(); for( int iSF=0; iSF < poFDefn->GetSubfieldCount(); iSF++ ) { DDFSubfieldDefn *poSFDefn = poFDefn->GetSubfield( iSF ); int iField; int nMaxBytes; const char * pachData = poSR->GetSubfieldData(poSFDefn, &nMaxBytes); /* -------------------------------------------------------------------- */ /* Identify the related DBF field, if any. */ /* -------------------------------------------------------------------- */ for( iField = 0; iField < hDBF->nFields; iField++ ) { if( EQUALN(poSFDefn->GetName(), hDBF->pszHeader+iField*32,10) ) break; } if( iField == hDBF->nFields ) iField = -1; /* -------------------------------------------------------------------- */ /* Handle each of the types. */ /* -------------------------------------------------------------------- */ switch( poSFDefn->GetType() ) { case DDFString: const char *pszValue; pszValue = poSFDefn->ExtractStringData(pachData, nMaxBytes, NULL); if( iField != -1 ) DBFWriteStringAttribute(hDBF, iRecord, iField, pszValue ); break; case DDFFloat: double dfValue; dfValue = poSFDefn->ExtractFloatData(pachData, nMaxBytes, NULL); if( iField != -1 ) DBFWriteDoubleAttribute( hDBF, iRecord, iField, dfValue ); break; case DDFInt: int nValue; nValue = poSFDefn->ExtractIntData(pachData, nMaxBytes, NULL); if( iField != -1 ) DBFWriteIntegerAttribute( hDBF, iRecord, iField, nValue ); break; default: break; } } /* next subfield */ }
static void AddPrimaryAttrToDBFSchema( DBFHandle hDBF, SDTSTransfer *poTransfer, char ** papszModuleList ) { for( int iModule = 0; papszModuleList != NULL && papszModuleList[iModule] != NULL; iModule++ ) { SDTSAttrReader *poAttrReader; /* -------------------------------------------------------------------- */ /* Get a reader on the desired module. */ /* -------------------------------------------------------------------- */ poAttrReader = (SDTSAttrReader *) poTransfer->GetLayerIndexedReader( poTransfer->FindLayer( papszModuleList[iModule] ) ); if( poAttrReader == NULL ) { printf( "Unable to open attribute module %s, skipping.\n" , papszModuleList[iModule] ); continue; } poAttrReader->Rewind(); /* -------------------------------------------------------------------- */ /* Read the first record so we can clone schema information off */ /* of it. */ /* -------------------------------------------------------------------- */ SDTSAttrRecord *poAttrFeature; poAttrFeature = (SDTSAttrRecord *) poAttrReader->GetNextFeature(); if( poAttrFeature == NULL ) { fprintf( stderr, "Didn't find any meaningful attribute records in %s.\n", papszModuleList[iModule] ); continue; } /* -------------------------------------------------------------------- */ /* Clone schema off the first record. Eventually we need to */ /* get the information out of the DDR record, but it isn't */ /* clear to me how to accomplish that with the SDTS++ API. */ /* */ /* The following approach may fail (dramatically) if some */ /* records do not include all subfields. Furthermore, no */ /* effort is made to make DBF field names unique. The SDTS */ /* attributes often have names much beyond the 14 character dbf */ /* limit which may result in non-unique attributes. */ /* -------------------------------------------------------------------- */ DDFFieldDefn *poFDefn = poAttrFeature->poATTR->GetFieldDefn(); int iSF; DDFField *poSR = poAttrFeature->poATTR; for( iSF=0; iSF < poFDefn->GetSubfieldCount(); iSF++ ) { DDFSubfieldDefn *poSFDefn = poFDefn->GetSubfield( iSF ); int nWidth = poSFDefn->GetWidth(); switch( poSFDefn->GetType() ) { case DDFString: if( nWidth == 0 ) { int nMaxBytes; const char * pachData = poSR->GetSubfieldData(poSFDefn, &nMaxBytes); nWidth = strlen(poSFDefn->ExtractStringData(pachData, nMaxBytes, NULL )); } DBFAddField( hDBF, poSFDefn->GetName(), FTString, nWidth, 0 ); break; case DDFInt: if( nWidth == 0 ) nWidth = 9; DBFAddField( hDBF, poSFDefn->GetName(), FTInteger, nWidth, 0 ); break; case DDFFloat: DBFAddField( hDBF, poSFDefn->GetName(), FTDouble, 18, 6 ); break; default: fprintf( stderr, "Dropping attribute `%s' of module `%s'. " "Type unsupported\n", poSFDefn->GetName(), papszModuleList[iModule] ); break; } } if( !poAttrReader->IsIndexed() ) delete poAttrFeature; } /* next module */ }
OGRSDTSLayer::OGRSDTSLayer( SDTSTransfer * poTransferIn, int iLayerIn, OGRSDTSDataSource * poDSIn ) : bPolygonsBuilt(FALSE) { poDS = poDSIn; poTransfer = poTransferIn; iLayer = iLayerIn; poReader = poTransfer->GetLayerIndexedReader( iLayer ); /* -------------------------------------------------------------------- */ /* Define the feature. */ /* -------------------------------------------------------------------- */ int iCATDEntry = poTransfer->GetLayerCATDEntry( iLayer ); poFeatureDefn = new OGRFeatureDefn(poTransfer->GetCATD()->GetEntryModule(iCATDEntry)); SetDescription( poFeatureDefn->GetName() ); poFeatureDefn->Reference(); poFeatureDefn->GetGeomFieldDefn(0)->SetSpatialRef(poDS->GetSpatialRef()); OGRFieldDefn oRecId( "RCID", OFTInteger ); poFeatureDefn->AddFieldDefn( &oRecId ); if( poTransfer->GetLayerType(iLayer) == SLTPoint ) { poFeatureDefn->SetGeomType( wkbPoint ); } else if( poTransfer->GetLayerType(iLayer) == SLTLine ) { poFeatureDefn->SetGeomType( wkbLineString ); oRecId.SetName( "SNID" ); poFeatureDefn->AddFieldDefn( &oRecId ); oRecId.SetName( "ENID" ); poFeatureDefn->AddFieldDefn( &oRecId ); } else if( poTransfer->GetLayerType(iLayer) == SLTPoly ) { poFeatureDefn->SetGeomType( wkbPolygon ); } else if( poTransfer->GetLayerType(iLayer) == SLTAttr ) { poFeatureDefn->SetGeomType( wkbNone ); } /* -------------------------------------------------------------------- */ /* Add schema from referenced attribute records. */ /* -------------------------------------------------------------------- */ char **papszATIDRefs = NULL; if( poTransfer->GetLayerType(iLayer) != SLTAttr ) papszATIDRefs = poReader->ScanModuleReferences(); else papszATIDRefs = CSLAddString( papszATIDRefs, poTransfer->GetCATD()->GetEntryModule(iCATDEntry) ); for( int iTable = 0; papszATIDRefs != NULL && papszATIDRefs[iTable] != NULL; iTable++ ) { SDTSAttrReader *poAttrReader; DDFFieldDefn *poFDefn; /* -------------------------------------------------------------------- */ /* Get the attribute table reader, and the associated user */ /* attribute field. */ /* -------------------------------------------------------------------- */ int nLayerIdx = poTransfer->FindLayer( papszATIDRefs[iTable] ); if( nLayerIdx < 0 ) continue; poAttrReader = (SDTSAttrReader *) poTransfer->GetLayerIndexedReader(nLayerIdx); if( poAttrReader == NULL ) continue; poFDefn = poAttrReader->GetModule()->FindFieldDefn( "ATTP" ); if( poFDefn == NULL ) poFDefn = poAttrReader->GetModule()->FindFieldDefn( "ATTS" ); if( poFDefn == NULL ) continue; /* -------------------------------------------------------------------- */ /* Process each user subfield on the attribute table into an */ /* OGR field definition. */ /* -------------------------------------------------------------------- */ for( int iSF=0; iSF < poFDefn->GetSubfieldCount(); iSF++ ) { DDFSubfieldDefn *poSFDefn = poFDefn->GetSubfield( iSF ); int nWidth = poSFDefn->GetWidth(); char *pszFieldName; if( poFeatureDefn->GetFieldIndex( poSFDefn->GetName() ) != -1 ) pszFieldName = CPLStrdup( CPLSPrintf( "%s_%s", papszATIDRefs[iTable], poSFDefn->GetName() ) ); else pszFieldName = CPLStrdup( poSFDefn->GetName() ); switch( poSFDefn->GetType() ) { case DDFString: { OGRFieldDefn oStrField( pszFieldName, OFTString ); if( nWidth != 0 ) oStrField.SetWidth( nWidth ); poFeatureDefn->AddFieldDefn( &oStrField ); } break; case DDFInt: { OGRFieldDefn oIntField( pszFieldName, OFTInteger ); if( nWidth != 0 ) oIntField.SetWidth( nWidth ); poFeatureDefn->AddFieldDefn( &oIntField ); } break; case DDFFloat: { OGRFieldDefn oRealField( pszFieldName, OFTReal ); // We don't have a precision in DDF files, so we never even // use the width. Otherwise with a precision of zero the // result would look like an integer. poFeatureDefn->AddFieldDefn( &oRealField ); } break; default: break; } CPLFree( pszFieldName ); } /* next iSF (subfield) */ } /* next iTable */ CSLDestroy( papszATIDRefs ); }
static void AssignAttrRecordToFeature( OGRFeature * poFeature, CPL_UNUSED SDTSTransfer * poTransfer, DDFField * poSR ) { /* -------------------------------------------------------------------- */ /* Process each subfield in the record. */ /* -------------------------------------------------------------------- */ DDFFieldDefn *poFDefn = poSR->GetFieldDefn(); for( int iSF=0; iSF < poFDefn->GetSubfieldCount(); iSF++ ) { DDFSubfieldDefn *poSFDefn = poFDefn->GetSubfield( iSF ); int iField; int nMaxBytes; const char * pachData = poSR->GetSubfieldData(poSFDefn, &nMaxBytes); /* -------------------------------------------------------------------- */ /* Indentify this field on the feature. */ /* -------------------------------------------------------------------- */ iField = poFeature->GetFieldIndex( poSFDefn->GetName() ); /* -------------------------------------------------------------------- */ /* Handle each of the types. */ /* -------------------------------------------------------------------- */ switch( poSFDefn->GetType() ) { case DDFString: const char *pszValue; pszValue = poSFDefn->ExtractStringData(pachData, nMaxBytes, NULL); if( iField != -1 ) poFeature->SetField( iField, pszValue ); break; case DDFFloat: double dfValue; dfValue = poSFDefn->ExtractFloatData(pachData, nMaxBytes, NULL); if( iField != -1 ) poFeature->SetField( iField, dfValue ); break; case DDFInt: int nValue; nValue = poSFDefn->ExtractIntData(pachData, nMaxBytes, NULL); if( iField != -1 ) poFeature->SetField( iField, nValue ); break; default: break; } } /* next subfield */ }
int main( int nArgc, char ** papszArgv ) { DDFModule oModule; const char *pszFilename = nullptr; int bFSPTHack = FALSE; for( int iArg = 1; iArg < nArgc; iArg++ ) { if( EQUAL(papszArgv[iArg],"-fspt_repeating") ) bFSPTHack = TRUE; else pszFilename = papszArgv[iArg]; } if( pszFilename == nullptr ) { printf( "Usage: 8211view filename\n" ); exit( 1 ); } /* -------------------------------------------------------------------- */ /* Open the file. Note that by default errors are reported to */ /* stderr, so we don't bother doing it ourselves. */ /* -------------------------------------------------------------------- */ if( !oModule.Open( pszFilename ) ) { exit( 1 ); } if( bFSPTHack ) { DDFFieldDefn *poFSPT = oModule.FindFieldDefn( "FSPT" ); if( poFSPT == nullptr ) fprintf( stderr, "unable to find FSPT field to set repeating flag.\n" ); else poFSPT->SetRepeatingFlag( TRUE ); } /* -------------------------------------------------------------------- */ /* Loop reading records till there are none left. */ /* -------------------------------------------------------------------- */ DDFRecord *poRecord = nullptr; int iRecord = 0; while( (poRecord = oModule.ReadRecord()) != nullptr ) { printf( "Record %d (%d bytes)\n", ++iRecord, poRecord->GetDataSize() ); /* ------------------------------------------------------------ */ /* Loop over each field in this particular record. */ /* ------------------------------------------------------------ */ for( int iField = 0; iField < poRecord->GetFieldCount(); iField++ ) { DDFField *poField = poRecord->GetField( iField ); ViewRecordField( poField ); } } }
int main(int nArgc, char* papszArgv[]) { const char *pszFilename = NULL, *pszOutFilename = NULL; DDFModule oModule; /* -------------------------------------------------------------------- */ /* Check arguments. */ /* -------------------------------------------------------------------- */ for( int iArg = 1; iArg < nArgc; iArg++ ) { if( pszFilename == NULL ) pszFilename = papszArgv[iArg]; else if( pszOutFilename == NULL ) pszOutFilename = papszArgv[iArg]; else { pszFilename = NULL; break; } } if( pszFilename == NULL ) { printf( "Usage: 8211createfromxml filename.xml outfilename\n" ); exit( 1 ); } CPLXMLNode* poRoot = CPLParseXMLFile( pszFilename ); if( poRoot == NULL ) { fprintf(stderr, "Cannot parse XML file '%s'\n", pszFilename); exit( 1 ); } CPLXMLNode* poXMLDDFModule = CPLSearchXMLNode(poRoot, "=DDFModule"); if( poXMLDDFModule == NULL ) { fprintf(stderr, "Cannot find DDFModule node in XML file '%s'\n", pszFilename); exit( 1 ); } /* Compute the size of the DDFField tag */ CPLXMLNode* psIter = poXMLDDFModule->psChild; int nSizeFieldTag = 0; while( psIter != NULL ) { if( psIter->eType == CXT_Element && strcmp(psIter->pszValue, "DDFFieldDefn") == 0 ) { const char* pszTag = CPLGetXMLValue(psIter, "tag", ""); if( nSizeFieldTag == 0 ) nSizeFieldTag = (int)strlen(pszTag); else if( nSizeFieldTag != (int)strlen(pszTag) ) { fprintf(stderr, "All fields have not the same tag size\n"); exit( 1 ); } } psIter = psIter->psNext; } char chInterchangeLevel = '3'; char chLeaderIden = 'L'; char chCodeExtensionIndicator = 'E'; char chVersionNumber = '1'; char chAppIndicator = ' '; const char *pszExtendedCharSet = " ! "; int nSizeFieldLength = 3; int nSizeFieldPos = 4; chInterchangeLevel = CPLGetXMLValue(poXMLDDFModule, "_interchangeLevel", CPLSPrintf("%c", chInterchangeLevel))[0]; chLeaderIden = CPLGetXMLValue(poXMLDDFModule, "_leaderIden", CPLSPrintf("%c", chLeaderIden))[0]; chCodeExtensionIndicator = CPLGetXMLValue(poXMLDDFModule, "_inlineCodeExtensionIndicator", CPLSPrintf("%c", chCodeExtensionIndicator))[0]; chVersionNumber = CPLGetXMLValue(poXMLDDFModule, "_versionNumber", CPLSPrintf("%c", chVersionNumber))[0]; chAppIndicator = CPLGetXMLValue(poXMLDDFModule, "_appIndicator", CPLSPrintf("%c", chAppIndicator))[0]; char szExtendedCharSet[4]; snprintf(szExtendedCharSet, sizeof(szExtendedCharSet), "%s", CPLGetXMLValue(poXMLDDFModule, "_extendedCharSet", pszExtendedCharSet)); pszExtendedCharSet = szExtendedCharSet; nSizeFieldLength = atoi(CPLGetXMLValue(poXMLDDFModule, "_sizeFieldLength", CPLSPrintf("%d", nSizeFieldLength))); nSizeFieldPos = atoi(CPLGetXMLValue(poXMLDDFModule, "_sizeFieldPos", CPLSPrintf("%d", nSizeFieldPos))); nSizeFieldTag = atoi(CPLGetXMLValue(poXMLDDFModule, "_sizeFieldTag", CPLSPrintf("%d", nSizeFieldTag))); oModule.Initialize(chInterchangeLevel, chLeaderIden, chCodeExtensionIndicator, chVersionNumber, chAppIndicator, pszExtendedCharSet, nSizeFieldLength, nSizeFieldPos, nSizeFieldTag); oModule.SetFieldControlLength(atoi(CPLGetXMLValue(poXMLDDFModule, "_fieldControlLength", CPLSPrintf("%d", oModule.GetFieldControlLength())))); int bCreated = FALSE; /* Create DDFFieldDefn and DDFRecord elements */ psIter = poXMLDDFModule->psChild; while( psIter != NULL ) { if( psIter->eType == CXT_Element && strcmp(psIter->pszValue, "DDFFieldDefn") == 0 ) { DDFFieldDefn* poFDefn = new DDFFieldDefn(); DDF_data_struct_code eStructCode = dsc_elementary; const char* pszStructCode = CPLGetXMLValue(psIter, "dataStructCode", ""); if( strcmp(pszStructCode, "elementary") == 0 ) eStructCode = dsc_elementary; else if( strcmp(pszStructCode, "vector") == 0 ) eStructCode = dsc_vector; else if( strcmp(pszStructCode, "array") == 0 ) eStructCode = dsc_array; else if( strcmp(pszStructCode, "concatenated") == 0 ) eStructCode = dsc_concatenated; DDF_data_type_code eTypeCode = dtc_char_string; const char* pszTypeCode = CPLGetXMLValue(psIter, "dataTypeCode", ""); if( strcmp(pszTypeCode, "char_string") == 0 ) eTypeCode = dtc_char_string; else if( strcmp(pszTypeCode, "implicit_point") == 0 ) eTypeCode = dtc_implicit_point; else if( strcmp(pszTypeCode, "explicit_point") == 0 ) eTypeCode = dtc_explicit_point; else if( strcmp(pszTypeCode, "explicit_point_scaled") == 0 ) eTypeCode = dtc_explicit_point_scaled; else if( strcmp(pszTypeCode, "char_bit_string") == 0 ) eTypeCode = dtc_char_bit_string; else if( strcmp(pszTypeCode, "bit_string") == 0 ) eTypeCode = dtc_bit_string; else if( strcmp(pszTypeCode, "mixed_data_type") == 0 ) eTypeCode = dtc_mixed_data_type; const char* pszFormatControls = CPLGetXMLValue(psIter, "formatControls", NULL); if( eStructCode != dsc_elementary ) pszFormatControls = NULL; const char* pszArrayDescr = CPLGetXMLValue(psIter, "arrayDescr", ""); if( eStructCode == dsc_vector ) pszArrayDescr = ""; else if( eStructCode == dsc_array ) pszArrayDescr = "*"; poFDefn->Create( CPLGetXMLValue(psIter, "tag", ""), CPLGetXMLValue(psIter, "fieldName", ""), pszArrayDescr, eStructCode, eTypeCode, pszFormatControls ); CPLXMLNode* psSubIter = psIter->psChild; while( psSubIter != NULL ) { if( psSubIter->eType == CXT_Element && strcmp(psSubIter->pszValue, "DDFSubfieldDefn") == 0 ) { poFDefn->AddSubfield( CPLGetXMLValue(psSubIter, "name", ""), CPLGetXMLValue(psSubIter, "format", "") ); } psSubIter = psSubIter->psNext; } pszFormatControls = CPLGetXMLValue(psIter, "formatControls", NULL); if( pszFormatControls ) poFDefn->SetFormatControls(pszFormatControls); oModule.AddField( poFDefn ); } else if( psIter->eType == CXT_Element && strcmp(psIter->pszValue, "DDFRecord") == 0 ) { //const bool bFirstRecord = !bCreated; if( !bCreated ) { oModule.Create( pszOutFilename ); bCreated = TRUE; } DDFRecord *poRec = new DDFRecord( &oModule ); std::map<std::string, int> oMapField; //if( !bFirstRecord ) // poRec->SetReuseHeader(atoi(CPLGetXMLValue(psIter, "reuseHeader", CPLSPrintf("%d", poRec->GetReuseHeader())))); poRec->SetSizeFieldLength(atoi(CPLGetXMLValue(psIter, "_sizeFieldLength", CPLSPrintf("%d", poRec->GetSizeFieldLength())))); poRec->SetSizeFieldPos(atoi(CPLGetXMLValue(psIter, "_sizeFieldPos", CPLSPrintf("%d", poRec->GetSizeFieldPos())))); poRec->SetSizeFieldTag(atoi(CPLGetXMLValue(psIter, "_sizeFieldTag", CPLSPrintf("%d", poRec->GetSizeFieldTag())))); CPLXMLNode* psSubIter = psIter->psChild; while( psSubIter != NULL ) { if( psSubIter->eType == CXT_Element && strcmp(psSubIter->pszValue, "DDFField") == 0 ) { const char* pszFieldName = CPLGetXMLValue(psSubIter, "name", ""); DDFFieldDefn* poFieldDefn = oModule.FindFieldDefn( pszFieldName ); if( poFieldDefn == NULL ) { fprintf(stderr, "Can't find field '%s'\n", pszFieldName ); exit(1); } int nFieldOcc = oMapField[pszFieldName]; oMapField[pszFieldName] ++ ; DDFField *poField = poRec->AddField( poFieldDefn ); const char* pszValue = CPLGetXMLValue(psSubIter, "value", NULL); if( pszValue != NULL && STARTS_WITH(pszValue, "0x") ) { pszValue += 2; int nDataLen = (int)strlen(pszValue) / 2; char* pabyData = (char*) malloc(nDataLen); for(int i=0;i<nDataLen;i++) { char c; int nHigh, nLow; c = pszValue[2*i]; if( c >= 'A' && c <= 'F' ) nHigh = 10 + c - 'A'; else nHigh = c - '0'; c = pszValue[2*i + 1]; if( c >= 'A' && c <= 'F' ) nLow = 10 + c - 'A'; else nLow = c - '0'; pabyData[i] = (nHigh << 4) + nLow; } poRec->SetFieldRaw( poField, nFieldOcc, (const char *) pabyData, nDataLen ); free(pabyData); } else { CPLXMLNode* psSubfieldIter = psSubIter->psChild; std::map<std::string, int> oMapSubfield; while( psSubfieldIter != NULL ) { if( psSubfieldIter->eType == CXT_Element && strcmp(psSubfieldIter->pszValue, "DDFSubfield") == 0 ) { const char* pszSubfieldName = CPLGetXMLValue(psSubfieldIter, "name", ""); const char* pszSubfieldType = CPLGetXMLValue(psSubfieldIter, "type", ""); const char* pszSubfieldValue = CPLGetXMLValue(psSubfieldIter, NULL, ""); int nOcc = oMapSubfield[pszSubfieldName]; oMapSubfield[pszSubfieldName] ++ ; if( strcmp(pszSubfieldType, "float") == 0 ) { poRec->SetFloatSubfield( pszFieldName, nFieldOcc, pszSubfieldName, nOcc, CPLAtof(pszSubfieldValue) ); } else if( strcmp(pszSubfieldType, "integer") == 0 ) { poRec->SetIntSubfield( pszFieldName, nFieldOcc, pszSubfieldName, nOcc, atoi(pszSubfieldValue) ); } else if( strcmp(pszSubfieldType, "string") == 0 ) { poRec->SetStringSubfield( pszFieldName, nFieldOcc, pszSubfieldName, nOcc, pszSubfieldValue ); } else if( strcmp(pszSubfieldType, "binary") == 0 && STARTS_WITH(pszSubfieldValue, "0x") ) { pszSubfieldValue += 2; int nDataLen = (int)strlen(pszSubfieldValue) / 2; char* pabyData = (char*) malloc(nDataLen); for(int i=0;i<nDataLen;i++) { char c; int nHigh, nLow; c = pszSubfieldValue[2*i]; if( c >= 'A' && c <= 'F' ) nHigh = 10 + c - 'A'; else nHigh = c - '0'; c = pszSubfieldValue[2*i + 1]; if( c >= 'A' && c <= 'F' ) nLow = 10 + c - 'A'; else nLow = c - '0'; pabyData[i] = (nHigh << 4) + nLow; } poRec->SetStringSubfield( pszFieldName, nFieldOcc, pszSubfieldName, nOcc, pabyData, nDataLen ); free(pabyData); } } psSubfieldIter = psSubfieldIter->psNext; } } } psSubIter = psSubIter->psNext; } poRec->Write(); delete poRec; } psIter = psIter->psNext; } CPLDestroyXMLNode(poRoot); oModule.Close(); return 0; }
int S57Writer::CreateS57File( const char *pszFilename ) { DDFModule oModule; DDFFieldDefn *poFDefn; Close(); nNext0001Index = 1; /* -------------------------------------------------------------------- */ /* Create and initialize new module. */ /* -------------------------------------------------------------------- */ poModule = new DDFModule(); poModule->Initialize(); /* -------------------------------------------------------------------- */ /* Create the '0000' definition. */ /* -------------------------------------------------------------------- */ poFDefn = new DDFFieldDefn(); poFDefn->Create( "0000", "", "0001DSIDDSIDDSSI0001DSPM0001VRIDVRIDATTVVRIDVRPCVRIDVRPTVRIDSGCCVRIDSG2DVRIDSG3D0001FRIDFRIDFOIDFRIDATTFFRIDNATFFRIDFFPCFRIDFFPTFRIDFSPCFRIDFSPT", dsc_elementary, dtc_char_string ); poModule->AddField( poFDefn ); /* -------------------------------------------------------------------- */ /* Create the '0001' definition. */ /* -------------------------------------------------------------------- */ poFDefn = new DDFFieldDefn(); poFDefn->Create( "0001", "ISO 8211 Record Identifier", "", dsc_elementary, dtc_bit_string, "(b12)" ); poModule->AddField( poFDefn ); /* -------------------------------------------------------------------- */ /* Create the DSID field. */ /* -------------------------------------------------------------------- */ poFDefn = new DDFFieldDefn(); poFDefn->Create( "DSID", "Data set identification field", "", dsc_vector, dtc_mixed_data_type ); poFDefn->AddSubfield( "RCNM", "b11" ); poFDefn->AddSubfield( "RCID", "b14" ); poFDefn->AddSubfield( "EXPP", "b11" ); poFDefn->AddSubfield( "INTU", "b11" ); poFDefn->AddSubfield( "DSNM", "A" ); poFDefn->AddSubfield( "EDTN", "A" ); poFDefn->AddSubfield( "UPDN", "A" ); poFDefn->AddSubfield( "UADT", "A(8)" ); poFDefn->AddSubfield( "ISDT", "A(8)" ); poFDefn->AddSubfield( "STED", "R(4)" ); poFDefn->AddSubfield( "PRSP", "b11" ); poFDefn->AddSubfield( "PSDN", "A" ); poFDefn->AddSubfield( "PRED", "A" ); poFDefn->AddSubfield( "PROF", "b11" ); poFDefn->AddSubfield( "AGEN", "b12" ); poFDefn->AddSubfield( "COMT", "A" ); poModule->AddField( poFDefn ); /* -------------------------------------------------------------------- */ /* Create the DSSI field. */ /* -------------------------------------------------------------------- */ poFDefn = new DDFFieldDefn(); poFDefn->Create( "DSSI", "Data set structure information field", "", dsc_vector, dtc_mixed_data_type ); poFDefn->AddSubfield( "DSTR", "b11" ); poFDefn->AddSubfield( "AALL", "b11" ); poFDefn->AddSubfield( "NALL", "b11" ); poFDefn->AddSubfield( "NOMR", "b14" ); poFDefn->AddSubfield( "NOCR", "b14" ); poFDefn->AddSubfield( "NOGR", "b14" ); poFDefn->AddSubfield( "NOLR", "b14" ); poFDefn->AddSubfield( "NOIN", "b14" ); poFDefn->AddSubfield( "NOCN", "b14" ); poFDefn->AddSubfield( "NOED", "b14" ); poFDefn->AddSubfield( "NOFA", "b14" ); poModule->AddField( poFDefn ); /* -------------------------------------------------------------------- */ /* Create the DSPM field. */ /* -------------------------------------------------------------------- */ poFDefn = new DDFFieldDefn(); poFDefn->Create( "DSPM", "Data set parameter field", "", dsc_vector, dtc_mixed_data_type ); poFDefn->AddSubfield( "RCNM", "b11" ); poFDefn->AddSubfield( "RCID", "b14" ); poFDefn->AddSubfield( "HDAT", "b11" ); poFDefn->AddSubfield( "VDAT", "b11" ); poFDefn->AddSubfield( "SDAT", "b11" ); poFDefn->AddSubfield( "CSCL", "b14" ); poFDefn->AddSubfield( "DUNI", "b11" ); poFDefn->AddSubfield( "HUNI", "b11" ); poFDefn->AddSubfield( "PUNI", "b11" ); poFDefn->AddSubfield( "COUN", "b11" ); poFDefn->AddSubfield( "COMF", "b14" ); poFDefn->AddSubfield( "SOMF", "b14" ); poFDefn->AddSubfield( "COMT", "A" ); poModule->AddField( poFDefn ); /* -------------------------------------------------------------------- */ /* Create the VRID field. */ /* -------------------------------------------------------------------- */ poFDefn = new DDFFieldDefn(); poFDefn->Create( "VRID", "Vector record identifier field", "", dsc_vector, dtc_mixed_data_type ); poFDefn->AddSubfield( "RCNM", "b11" ); poFDefn->AddSubfield( "RCID", "b14" ); poFDefn->AddSubfield( "RVER", "b12" ); poFDefn->AddSubfield( "RUIN", "b11" ); poModule->AddField( poFDefn ); /* -------------------------------------------------------------------- */ /* Create the VRPC field. */ /* -------------------------------------------------------------------- */ poFDefn = new DDFFieldDefn(); poFDefn->Create( "VRPC", "Vector Record Pointer Control field", "", dsc_vector, dtc_mixed_data_type ); poFDefn->AddSubfield( "VPUI", "b11" ); poFDefn->AddSubfield( "VPIX", "b12" ); poFDefn->AddSubfield( "NVPT", "b12" ); poModule->AddField( poFDefn ); /* -------------------------------------------------------------------- */ /* Create the VRPT field. */ /* -------------------------------------------------------------------- */ poFDefn = new DDFFieldDefn(); poFDefn->Create( "VRPT", "Vector record pointer field", "*", dsc_array, dtc_mixed_data_type ); poFDefn->AddSubfield( "NAME", "B(40)" ); poFDefn->AddSubfield( "ORNT", "b11" ); poFDefn->AddSubfield( "USAG", "b11" ); poFDefn->AddSubfield( "TOPI", "b11" ); poFDefn->AddSubfield( "MASK", "b11" ); poModule->AddField( poFDefn ); /* -------------------------------------------------------------------- */ /* Create the ATTV field. */ /* -------------------------------------------------------------------- */ poFDefn = new DDFFieldDefn(); poFDefn->Create( "ATTV", "Vector record attribute field", "*", dsc_array, dtc_mixed_data_type ); poFDefn->AddSubfield( "ATTL", "b12" ); poFDefn->AddSubfield( "ATVL", "A" ); poModule->AddField( poFDefn ); /* -------------------------------------------------------------------- */ /* Create the SGCC field. */ /* -------------------------------------------------------------------- */ poFDefn = new DDFFieldDefn(); poFDefn->Create( "SGCC", "Coordinate Control Field", "", dsc_vector, dtc_mixed_data_type ); poFDefn->AddSubfield( "CCUI", "b11" ); poFDefn->AddSubfield( "CCIX", "b12" ); poFDefn->AddSubfield( "CCNC", "b12" ); poModule->AddField( poFDefn ); /* -------------------------------------------------------------------- */ /* Create the SG2D field. */ /* -------------------------------------------------------------------- */ poFDefn = new DDFFieldDefn(); poFDefn->Create( "SG2D", "2-D coordinate field", "*", dsc_array, dtc_mixed_data_type ); poFDefn->AddSubfield( "YCOO", "b24" ); poFDefn->AddSubfield( "XCOO", "b24" ); poModule->AddField( poFDefn ); /* -------------------------------------------------------------------- */ /* Create the SG3D field. */ /* -------------------------------------------------------------------- */ poFDefn = new DDFFieldDefn(); poFDefn->Create( "SG3D", "3-D coordinate (sounding array) field", "*", dsc_array, dtc_mixed_data_type ); poFDefn->AddSubfield( "YCOO", "b24" ); poFDefn->AddSubfield( "XCOO", "b24" ); poFDefn->AddSubfield( "VE3D", "b24" ); poModule->AddField( poFDefn ); /* -------------------------------------------------------------------- */ /* Create the FRID field. */ /* -------------------------------------------------------------------- */ poFDefn = new DDFFieldDefn(); poFDefn->Create( "FRID", "Feature record identifier field", "", dsc_vector, dtc_mixed_data_type ); poFDefn->AddSubfield( "RCNM", "b11" ); poFDefn->AddSubfield( "RCID", "b14" ); poFDefn->AddSubfield( "PRIM", "b11" ); poFDefn->AddSubfield( "GRUP", "b11" ); poFDefn->AddSubfield( "OBJL", "b12" ); poFDefn->AddSubfield( "RVER", "b12" ); poFDefn->AddSubfield( "RUIN", "b11" ); poModule->AddField( poFDefn ); /* -------------------------------------------------------------------- */ /* Create the FOID field. */ /* -------------------------------------------------------------------- */ poFDefn = new DDFFieldDefn(); poFDefn->Create( "FOID", "Feature object identifier field", "", dsc_vector, dtc_mixed_data_type ); poFDefn->AddSubfield( "AGEN", "b12" ); poFDefn->AddSubfield( "FIDN", "b14" ); poFDefn->AddSubfield( "FIDS", "b12" ); poModule->AddField( poFDefn ); /* -------------------------------------------------------------------- */ /* Create the ATTF field. */ /* -------------------------------------------------------------------- */ poFDefn = new DDFFieldDefn(); poFDefn->Create( "ATTF", "Feature record attribute field", "*", dsc_array, dtc_mixed_data_type ); poFDefn->AddSubfield( "ATTL", "b12" ); poFDefn->AddSubfield( "ATVL", "A" ); poModule->AddField( poFDefn ); /* -------------------------------------------------------------------- */ /* Create the NATF field. */ /* -------------------------------------------------------------------- */ poFDefn = new DDFFieldDefn(); poFDefn->Create( "NATF", "Feature record national attribute field", "*", dsc_array, dtc_mixed_data_type ); poFDefn->AddSubfield( "ATTL", "b12" ); poFDefn->AddSubfield( "ATVL", "A" ); poModule->AddField( poFDefn ); /* -------------------------------------------------------------------- */ /* Create the FFPC field. */ /* -------------------------------------------------------------------- */ poFDefn = new DDFFieldDefn(); poFDefn->Create( "FFPC", "Feature record to feature object pointer control field", "", dsc_vector, dtc_mixed_data_type ); poFDefn->AddSubfield( "FFUI", "b11" ); poFDefn->AddSubfield( "FFIX", "b12" ); poFDefn->AddSubfield( "NFPT", "b12" ); poModule->AddField( poFDefn ); /* -------------------------------------------------------------------- */ /* Create the FFPT field. */ /* -------------------------------------------------------------------- */ poFDefn = new DDFFieldDefn(); poFDefn->Create( "FFPT", "Feature record to feature object pointer field", "*", dsc_array, dtc_mixed_data_type ); poFDefn->AddSubfield( "LNAM", "B(64)" ); poFDefn->AddSubfield( "RIND", "b11" ); poFDefn->AddSubfield( "COMT", "A" ); poModule->AddField( poFDefn ); /* -------------------------------------------------------------------- */ /* Create the FSPC field. */ /* -------------------------------------------------------------------- */ poFDefn = new DDFFieldDefn(); poFDefn->Create( "FSPC", "Feature record to spatial record pointer control field", "", dsc_vector, dtc_mixed_data_type ); poFDefn->AddSubfield( "FSUI", "b11" ); poFDefn->AddSubfield( "FSIX", "b12" ); poFDefn->AddSubfield( "NSPT", "b12" ); poModule->AddField( poFDefn ); /* -------------------------------------------------------------------- */ /* Create the FSPT field. */ /* -------------------------------------------------------------------- */ poFDefn = new DDFFieldDefn(); poFDefn->Create( "FSPT", "Feature record to spatial record pointer field", "*", dsc_array, dtc_mixed_data_type ); poFDefn->AddSubfield( "NAME", "B(40)" ); poFDefn->AddSubfield( "ORNT", "b11" ); poFDefn->AddSubfield( "USAG", "b11" ); poFDefn->AddSubfield( "MASK", "b11" ); poModule->AddField( poFDefn ); /* -------------------------------------------------------------------- */ /* Create file. */ /* -------------------------------------------------------------------- */ if( !poModule->Create( pszFilename ) ) { delete poModule; poModule = NULL; return FALSE; } return TRUE; }