Sequence *OGRDODSSequenceLayer::FindSuperSequence( BaseType *poChild )
{
BaseType *poParent;
for( poParent = poChild->get_parent();
poParent != NULL;
poParent = poParent->get_parent() )
{
if( poParent->type() == dods_sequence_c )
{
return dynamic_cast<Sequence *>( poParent );
}
}
return NULL;
}
OGRFeature *OGRDODSSequenceLayer::GetFeature( GIntBig nFeatureId )
{
/* -------------------------------------------------------------------- */
/* Ensure we have the dataset. */
/* -------------------------------------------------------------------- */
if( !ProvideDataDDS() )
return NULL;
Sequence *seq = dynamic_cast<Sequence *>(poTargetVar);
/* -------------------------------------------------------------------- */
/* Figure out what the super and subsequence number this */
/* feature will be, and validate it. If there is not super */
/* sequence the feature id is the subsequence number. */
/* -------------------------------------------------------------------- */
int iSubSeq = -1;
if( nFeatureId < 0 || nFeatureId >= nRecordCount )
return NULL;
if( poSuperSeq == NULL )
iSubSeq = nFeatureId;
else
{
int nSeqOffset = 0, iSuperSeq;
// for now we just scan through till find find out what
// super sequence this in. In the long term we need a better (cached)
// approach that doesn't involve this quadratic cost.
for( iSuperSeq = 0;
iSuperSeq < nSuperSeqCount;
iSuperSeq++ )
{
if( nSeqOffset + panSubSeqSize[iSuperSeq] > nFeatureId )
{
iSubSeq = nFeatureId - nSeqOffset;
break;
}
nSeqOffset += panSubSeqSize[iSuperSeq];
}
CPLAssert( iSubSeq != -1 );
// Make sure we have the right target var ... the one
// corresponding to our current super sequence.
if( iSuperSeq != iLastSuperSeq )
{
iLastSuperSeq = iSuperSeq;
poTargetVar = poSuperSeq->var_value( iSuperSeq, pszSubSeqPath );
seq = dynamic_cast<Sequence *>(poTargetVar);
}
}
/* -------------------------------------------------------------------- */
/* Create the feature being read. */
/* -------------------------------------------------------------------- */
OGRFeature *poFeature;
poFeature = new OGRFeature( poFeatureDefn );
poFeature->SetFID( nFeatureId );
m_nFeaturesRead++;
/* -------------------------------------------------------------------- */
/* Process all the regular data fields. */
/* -------------------------------------------------------------------- */
int iField;
for( iField = 0; iField < poFeatureDefn->GetFieldCount(); iField++ )
{
if( papoFields[iField]->pszPathToSequence )
continue;
BaseType *poFieldVar = GetFieldValue( papoFields[iField], iSubSeq,
NULL );
if( poFieldVar == NULL )
continue;
switch( poFieldVar->type() )
{
case dods_byte_c:
{
signed char byVal;
void *pValPtr = &byVal;
poFieldVar->buf2val( &pValPtr );
poFeature->SetField( iField, byVal );
}
break;
case dods_int16_c:
{
GInt16 nIntVal;
void *pValPtr = &nIntVal;
poFieldVar->buf2val( &pValPtr );
poFeature->SetField( iField, nIntVal );
}
break;
case dods_uint16_c:
{
GUInt16 nIntVal;
void *pValPtr = &nIntVal;
poFieldVar->buf2val( &pValPtr );
poFeature->SetField( iField, nIntVal );
}
break;
case dods_int32_c:
{
GInt32 nIntVal;
void *pValPtr = &nIntVal;
poFieldVar->buf2val( &pValPtr );
poFeature->SetField( iField, nIntVal );
}
break;
case dods_uint32_c:
{
GUInt32 nIntVal;
void *pValPtr = &nIntVal;
poFieldVar->buf2val( &pValPtr );
poFeature->SetField( iField, (int) nIntVal );
}
break;
case dods_float32_c:
poFeature->SetField( iField,
dynamic_cast<Float32 *>(poFieldVar)->value());
break;
case dods_float64_c:
poFeature->SetField( iField,
dynamic_cast<Float64 *>(poFieldVar)->value());
break;
case dods_str_c:
case dods_url_c:
{
string *poStrVal = NULL;
poFieldVar->buf2val( (void **) &poStrVal );
poFeature->SetField( iField, poStrVal->c_str() );
delete poStrVal;
}
break;
default:
break;
}
}
/* -------------------------------------------------------------------- */
/* Handle data nested in sequences. */
/* -------------------------------------------------------------------- */
for( iField = 0; iField < poFeatureDefn->GetFieldCount(); iField++ )
{
OGRDODSFieldDefn *poFD = papoFields[iField];
const char *pszPathFromSubSeq;
if( poFD->pszPathToSequence == NULL )
continue;
CPLAssert( strlen(poFD->pszPathToSequence)
< strlen(poFD->pszFieldName)-1 );
if( strstr(poFD->pszFieldName,poFD->pszPathToSequence) != NULL )
pszPathFromSubSeq =
strstr(poFD->pszFieldName,poFD->pszPathToSequence)
+ strlen(poFD->pszPathToSequence) + 1;
else
continue;
/* -------------------------------------------------------------------- */
/* Get the sequence out of which this variable will be collected. */
/* -------------------------------------------------------------------- */
BaseType *poFieldVar = seq->var_value( iSubSeq,
poFD->pszPathToSequence );
Sequence *poSubSeq;
int nSubSeqCount;
if( poFieldVar == NULL )
continue;
poSubSeq = dynamic_cast<Sequence *>( poFieldVar );
if( poSubSeq == NULL )
continue;
nSubSeqCount = poSubSeq->number_of_rows();
/* -------------------------------------------------------------------- */
/* Allocate array to put values into. */
/* -------------------------------------------------------------------- */
OGRFieldDefn *poOFD = poFeature->GetFieldDefnRef( iField );
int *panIntList = NULL;
double *padfDblList = NULL;
char **papszStrList = NULL;
if( poOFD->GetType() == OFTIntegerList )
{
panIntList = (int *) CPLCalloc(sizeof(int),nSubSeqCount);
}
else if( poOFD->GetType() == OFTRealList )
{
padfDblList = (double *) CPLCalloc(sizeof(double),nSubSeqCount);
}
else if( poOFD->GetType() == OFTStringList )
{
papszStrList = (char **) CPLCalloc(sizeof(char*),nSubSeqCount+1);
}
else
continue;
/* -------------------------------------------------------------------- */
/* Loop, fetching subsequence values. */
/* -------------------------------------------------------------------- */
int iSubIndex;
for( iSubIndex = 0; iSubIndex < nSubSeqCount; iSubIndex++ )
{
poFieldVar = poSubSeq->var_value( iSubIndex, pszPathFromSubSeq );
if( poFieldVar == NULL )
continue;
switch( poFieldVar->type() )
{
case dods_byte_c:
{
signed char byVal;
void *pValPtr = &byVal;
poFieldVar->buf2val( &pValPtr );
panIntList[iSubIndex] = byVal;
}
break;
case dods_int16_c:
{
GInt16 nIntVal;
void *pValPtr = &nIntVal;
poFieldVar->buf2val( &pValPtr );
panIntList[iSubIndex] = nIntVal;
}
break;
case dods_uint16_c:
{
GUInt16 nIntVal;
void *pValPtr = &nIntVal;
poFieldVar->buf2val( &pValPtr );
panIntList[iSubIndex] = nIntVal;
}
break;
case dods_int32_c:
{
GInt32 nIntVal;
void *pValPtr = &nIntVal;
poFieldVar->buf2val( &pValPtr );
panIntList[iSubIndex] = nIntVal;
}
break;
case dods_uint32_c:
{
GUInt32 nIntVal;
void *pValPtr = &nIntVal;
poFieldVar->buf2val( &pValPtr );
panIntList[iSubIndex] = nIntVal;
}
break;
case dods_float32_c:
padfDblList[iSubIndex] =
dynamic_cast<Float32 *>(poFieldVar)->value();
break;
case dods_float64_c:
padfDblList[iSubIndex] =
dynamic_cast<Float64 *>(poFieldVar)->value();
break;
case dods_str_c:
case dods_url_c:
{
string *poStrVal = NULL;
poFieldVar->buf2val( (void **) &poStrVal );
papszStrList[iSubIndex] = CPLStrdup( poStrVal->c_str() );
delete poStrVal;
}
break;
default:
break;
}
}
/* -------------------------------------------------------------------- */
/* Apply back to feature. */
/* -------------------------------------------------------------------- */
if( poOFD->GetType() == OFTIntegerList )
{
poFeature->SetField( iField, nSubSeqCount, panIntList );
CPLFree(panIntList);
}
else if( poOFD->GetType() == OFTRealList )
{
poFeature->SetField( iField, nSubSeqCount, padfDblList );
CPLFree(padfDblList);
}
else if( poOFD->GetType() == OFTStringList )
{
poFeature->SetField( iField, papszStrList );
CSLDestroy( papszStrList );
}
}
/* ==================================================================== */
/* Fetch the geometry. */
/* ==================================================================== */
if( oXField.bValid && oYField.bValid )
{
int iXField = poFeature->GetFieldIndex( oXField.pszFieldName );
int iYField = poFeature->GetFieldIndex( oYField.pszFieldName );
int iZField = -1;
if( oZField.bValid )
iZField = poFeature->GetFieldIndex(oZField.pszFieldName);
/* -------------------------------------------------------------------- */
/* If we can't find the values in attributes then use the more */
/* general mechanism to fetch the value. */
/* -------------------------------------------------------------------- */
if( iXField == -1 || iYField == -1
|| (oZField.bValid && iZField == -1) )
{
poFeature->SetGeometryDirectly(
new OGRPoint( GetFieldValueAsDouble( &oXField, iSubSeq ),
GetFieldValueAsDouble( &oYField, iSubSeq ),
GetFieldValueAsDouble( &oZField, iSubSeq ) ) );
}
/* -------------------------------------------------------------------- */
/* If the fields are list values, then build a linestring. */
/* -------------------------------------------------------------------- */
else if( poFeature->GetFieldDefnRef(iXField)->GetType() == OFTRealList
&& poFeature->GetFieldDefnRef(iYField)->GetType() == OFTRealList )
{
const double *padfX, *padfY, *padfZ = NULL;
int nPointCount, i;
OGRLineString *poLS = new OGRLineString();
padfX = poFeature->GetFieldAsDoubleList( iXField, &nPointCount );
padfY = poFeature->GetFieldAsDoubleList( iYField, &nPointCount );
if( iZField != -1 )
padfZ = poFeature->GetFieldAsDoubleList(iZField,&nPointCount);
poLS->setPoints( nPointCount, (double *) padfX, (double *) padfY,
(double *) padfZ );
// Make a pass clearing out NaN or Inf values.
for( i = 0; i < nPointCount; i++ )
{
double dfX = poLS->getX(i);
double dfY = poLS->getY(i);
double dfZ = poLS->getZ(i);
int bReset = FALSE;
if( OGRDODSIsDoubleInvalid( &dfX ) )
{
dfX = 0.0;
bReset = TRUE;
}
if( OGRDODSIsDoubleInvalid( &dfY ) )
{
dfY = 0.0;
bReset = TRUE;
}
if( OGRDODSIsDoubleInvalid( &dfZ ) )
{
dfZ = 0.0;
bReset = TRUE;
}
if( bReset )
poLS->setPoint( i, dfX, dfY, dfZ );
}
poFeature->SetGeometryDirectly( poLS );
}
/* -------------------------------------------------------------------- */
/* Otherwise build a point. */
/* -------------------------------------------------------------------- */
else
{
poFeature->SetGeometryDirectly(
new OGRPoint(
poFeature->GetFieldAsDouble( iXField ),
poFeature->GetFieldAsDouble( iYField ),
poFeature->GetFieldAsDouble( iZField ) ) );
}
}
return poFeature;
}
Instance MemoryDump::getNextInstance(const QString& component,
const Instance& instance,
KnowledgeSources src) const
{
Instance result;
QString typeString, symbol, offsetString, candidate, arrayIndexString;
bool okay;
// quint32 compatibleCnt = 0;
// A component should have the form (symbol(-offset)?)?symbol(<candidate>)?([index])?
#define SYMBOL "[A-Za-z0-9_]+"
#define NUMBER "\\d+"
QRegExp re(
"^\\s*(?:"
"\\(\\s*"
"(" SYMBOL ")"
"(?:"
"\\s*-\\s*(" SYMBOL ")"
")?"
"\\s*\\)"
")?"
"\\s*(" SYMBOL ")\\s*"
"(?:<\\s*(" NUMBER ")\\s*>\\s*)?"
"((?:\\[\\s*" NUMBER "\\s*\\]\\s*)*)\\s*");
if (!re.exactMatch(component)) {
queryError(QString("Could not parse a part of the query string: %1")
.arg(component));
}
// Set variables according to the matching
typeString = re.cap(1);
offsetString = re.cap(2).trimmed();
symbol = re.cap(3);
candidate = re.cap(4);
arrayIndexString = re.cap(5).trimmed();
int candidateIndex = candidate.isEmpty() ? -1 : candidate.toInt();
// debugmsg(QString("1: %1, 2: %2, 3: %3, 4: %4, 5: %5")
// .arg(re.cap(1))
// .arg(re.cap(2))
// .arg(re.cap(3))
// .arg(re.cap(4))
// .arg(re.cap(5)));
// A candidate index of 0 means to ignore the alternative types
if (candidateIndex == 0)
src = static_cast<KnowledgeSources>(src|ksNoAltTypes);
// If the given instance is Null, we interpret this as the first component
// in the query string and will therefore try to resolve the variable.
if (!instance.isValid()) {
Variable* v = _factory->findVarByName(symbol);
if (!v)
queryError(QString("Variable does not exist: %1").arg(symbol));
if (candidateIndex > 0) {
if (v->altRefTypeCount() < candidateIndex)
queryError(QString("Variable \"%1\" does not have a candidate "
"with index %2")
.arg(symbol)
.arg(candidateIndex));
result = v->altRefTypeInstance(_vmem, candidateIndex - 1);
}
else {
result = v->toInstance(_vmem, BaseType::trLexical, src);
}
}
else {
// Dereference any pointers/arrays first
result = instance.dereference(BaseType::trAnyNonNull);
// Did we get a null instance?
if (!(result.type()->type() & StructOrUnion) &&
(result.isNull() || !result.toPointer()))
queryError(QString("Member \"%1\" is null")
.arg(result.fullName()));
// We have a instance therefore we resolve the member
if (!(result.type()->type() & StructOrUnion))
queryError(QString("Member \"%1\" is not a struct or union")
.arg(result.fullName()));
if (!result.memberExists(symbol))
queryError(QString("Struct \"%1\" has no member named \"%2\"")
.arg(result.typeName())
.arg(symbol));
// Do we have a candidate index?
if (candidateIndex > 0) {
if (result.memberCandidatesCount(symbol) < candidateIndex)
queryError(QString("Member \"%1\" does not have a candidate "
"with index %2")
.arg(symbol)
.arg(candidateIndex));
result = result.memberCandidate(symbol, candidateIndex - 1);
}
else {
result = result.member(symbol, BaseType::trLexical, 0, src);
}
}
if (!result.isValid())
return result;
// Cast the instance if necessary
if (!typeString.isEmpty()) {
quint32 offset = 0;
// Is a offset given?
if (!offsetString.isEmpty()) {
// Is the offset given as string or as int?
offset = offsetString.toUInt(&okay, 10);
if (!okay) {
// String.
BaseType* type = getType(typeString);
if (!type ||
!(type->type() & StructOrUnion))
queryError(QString("The given type \"%1\" is not a struct "
"or union and therefore has no offset")
.arg(typeString));
Structured* structd = dynamic_cast<Structured *>(type);
if (!structd->memberExists(offsetString)) {
queryError(QString("Struct of type \"%1\" has no member "
"named \"%2\"")
.arg(typeString)
.arg(offsetString));
}
else {
StructuredMember* structdMember =
structd->member(offsetString);
offset = structdMember->offset();
}
}
}
// Get address
size_t address;
if (result.type()->type() & (rtPointer))
address = (size_t)result.toPointer() - offset;
else
address = result.address() - offset;
result = getInstanceAt(typeString, address, result.fullNameComponents());
}
// Add array index
if (!arrayIndexString.isEmpty()) {
QRegExp reArrayIndex("\\[\\s*(" NUMBER ")\\s*\\]\\s*");
QStringList matches;
int strpos = 0;
while (strpos < arrayIndexString.size() &&
(strpos = arrayIndexString.indexOf(reArrayIndex, strpos)) >= 0)
{
matches.append(reArrayIndex.cap(1));
strpos += reArrayIndex.cap(0).size();
}
for (int i = 0; i < matches.count(); ++i) {
quint32 arrayIndex = matches[i].toUInt(&okay, 10);
if (okay) {
// Is the result already an instance list?
if (result.isList()) {
InstanceList list(result.toList());
if (arrayIndex < (quint32)list.size())
result = list[arrayIndex];
else
queryError(QString("Given array index %1 is out of bounds.")
.arg(arrayIndex));
}
else {
// Is this a pointer or an array type?
Instance tmp = result.arrayElem(arrayIndex);
if (!tmp.isNull())
result = tmp.dereference(BaseType::trLexical);
// Manually update the address
else {
result.addToAddress(arrayIndex * result.type()->size());
result.setName(QString("%1[%2]").arg(result.name()).arg(arrayIndex));
}
}
}
else {
queryError(QString("Given array index %1 could not be converted "
"to a number.")
.arg(matches[i]));
}
}
}
// Try to dereference this instance as deep as possible
return result.dereference(BaseType::trLexicalAndPointers);
}
bool OGRDODSGridLayer::ProvideDataDDS()
{
if( bDataLoaded )
return poTargetVar != NULL;
const bool bResult = OGRDODSLayer::ProvideDataDDS();
if( !bResult )
return bResult;
for( int iArray=0; iArray < nArrayRefCount; iArray++ )
{
OGRDODSArrayRef *poRef = paoArrayRefs + iArray;
BaseType *poTarget = poDataDDS->var( poRef->pszName );
// Reset ref array pointer to point in DataDDS result.
if( poTarget->type() == dods_grid_c )
{
Grid *poGrid = dynamic_cast<Grid *>( poTarget );
poRef->poArray = dynamic_cast<Array *>(poGrid->array_var());
if( iArray == 0 )
poTargetGrid = poGrid;
}
else if( poTarget->type() == dods_array_c )
{
poRef->poArray = dynamic_cast<Array *>( poTarget );
}
else
{
CPLAssert( false );
return false;
}
if( iArray == 0 )
poTargetArray = poRef->poArray;
// Allocate appropriate raw data array, and pull out data into it.
poRef->pRawData = CPLMalloc( poRef->poArray->width() );
poRef->poArray->buf2val( &(poRef->pRawData) );
}
// Setup pointers to each of the map objects.
if( poTargetGrid != NULL )
{
int iMap = 0;
Grid::Map_iter iterMap;
for( iterMap = poTargetGrid->map_begin();
iterMap != poTargetGrid->map_end();
iterMap++, iMap++ )
{
paoDimensions[iMap].poMap = dynamic_cast<Array *>(*iterMap);
if( paoDimensions[iMap].poMap == NULL )
return false;
paoDimensions[iMap].pRawData =
CPLMalloc( paoDimensions[iMap].poMap->width() );
paoDimensions[iMap].poMap->buf2val( &(paoDimensions[iMap].pRawData) );
}
}
return bResult;
}
OGRDODSGridLayer::OGRDODSGridLayer( OGRDODSDataSource *poDSIn,
const char *pszTargetIn,
AttrTable *poOGRLayerInfoIn ) :
OGRDODSLayer( poDSIn, pszTargetIn, poOGRLayerInfoIn ),
poTargetGrid(NULL),
poTargetArray(NULL),
nArrayRefCount(0),
paoArrayRefs(NULL),
nDimCount(0),
paoDimensions(NULL),
nMaxRawIndex(0)
{
/* -------------------------------------------------------------------- */
/* What is the layer name? */
/* -------------------------------------------------------------------- */
string oLayerName;
const char *pszLayerName = pszTargetIn;
if( poOGRLayerInfo != NULL )
{
oLayerName = poOGRLayerInfo->get_attr( "layer_name" );
if( strlen(oLayerName.c_str()) > 0 )
pszLayerName = oLayerName.c_str();
}
poFeatureDefn = new OGRFeatureDefn( pszLayerName );
poFeatureDefn->Reference();
/* -------------------------------------------------------------------- */
/* Fetch the target variable. */
/* -------------------------------------------------------------------- */
BaseType *poTargVar = poDS->poDDS->var( pszTargetIn );
if( poTargVar->type() == dods_grid_c )
{
poTargetGrid = dynamic_cast<Grid *>( poTargVar );
poTargetArray = dynamic_cast<Array *>(poTargetGrid->array_var());
}
else if( poTargVar->type() == dods_array_c )
{
poTargetGrid = NULL;
poTargetArray = dynamic_cast<Array *>( poTargVar );
}
else
{
CPLAssert( false );
return;
}
/* -------------------------------------------------------------------- */
/* Count arrays in use. */
/* -------------------------------------------------------------------- */
AttrTable *poExtraContainers = NULL;
nArrayRefCount = 1; // primary target.
if( poOGRLayerInfo != NULL )
poExtraContainers = poOGRLayerInfo->find_container("extra_containers");
if( poExtraContainers != NULL )
{
AttrTable::Attr_iter dv_i;
for( dv_i = poExtraContainers->attr_begin();
dv_i != poExtraContainers->attr_end(); dv_i++ )
{
nArrayRefCount++;
}
}
/* -------------------------------------------------------------------- */
/* Collect extra_containers. */
/* -------------------------------------------------------------------- */
paoArrayRefs = new OGRDODSArrayRef[nArrayRefCount];
paoArrayRefs[0].pszName = CPLStrdup( pszTargetIn );
paoArrayRefs[0].poArray = poTargetArray;
nArrayRefCount = 1;
if( poExtraContainers != NULL )
{
AttrTable::Attr_iter dv_i;
for( dv_i = poExtraContainers->attr_begin();
dv_i != poExtraContainers->attr_end(); dv_i++ )
{
const char *pszTargetName=poExtraContainers->get_attr(dv_i).c_str();
BaseType *poExtraTarget = poDS->poDDS->var( pszTargetName );
if( poExtraTarget == NULL )
{
CPLError( CE_Warning, CPLE_AppDefined,
"Unable to find extra_container '%s', skipping.",
pszTargetName );
continue;
}
if( poExtraTarget->type() == dods_array_c )
paoArrayRefs[nArrayRefCount].poArray =
dynamic_cast<Array *>( poExtraTarget );
else if( poExtraTarget->type() == dods_grid_c )
{
Grid *poGrid = dynamic_cast<Grid *>( poExtraTarget );
paoArrayRefs[nArrayRefCount].poArray =
dynamic_cast<Array *>( poGrid->array_var() );
}
else
{
CPLError( CE_Warning, CPLE_AppDefined,
"Target container '%s' is not grid or array, skipping.",
pszTargetName );
continue;
}
paoArrayRefs[nArrayRefCount++].pszName = CPLStrdup(pszTargetName);
}
}
/* -------------------------------------------------------------------- */
/* Collect dimension information. */
/* -------------------------------------------------------------------- */
int iDim;
Array::Dim_iter iterDim;
nDimCount = poTargetArray->dimensions();
paoDimensions = new OGRDODSDim[nDimCount];
nMaxRawIndex = 1;
for( iterDim = poTargetArray->dim_begin(), iDim = 0;
iterDim != poTargetArray->dim_end();
iterDim++, iDim++ )
{
paoDimensions[iDim].pszDimName =
CPLStrdup(poTargetArray->dimension_name(iterDim).c_str());
paoDimensions[iDim].nDimStart =
poTargetArray->dimension_start(iterDim);
paoDimensions[iDim].nDimEnd =
poTargetArray->dimension_stop(iterDim);
paoDimensions[iDim].nDimStride =
poTargetArray->dimension_stride(iterDim);
paoDimensions[iDim].poMap = NULL;
paoDimensions[iDim].nDimEntries =
(paoDimensions[iDim].nDimEnd + 1 - paoDimensions[iDim].nDimStart
+ paoDimensions[iDim].nDimStride - 1)
/ paoDimensions[iDim].nDimStride;
nMaxRawIndex *= paoDimensions[iDim].nDimEntries;
}
/* -------------------------------------------------------------------- */
/* If we are working with a grid, collect the maps. */
/* -------------------------------------------------------------------- */
if( poTargetGrid != NULL )
{
int iMap;
Grid::Map_iter iterMap;
for( iterMap = poTargetGrid->map_begin(), iMap = 0;
iterMap != poTargetGrid->map_end();
iterMap++, iMap++ )
{
paoDimensions[iMap].poMap = dynamic_cast<Array *>(*iterMap);
}
CPLAssert( iMap == nDimCount );
}
/* -------------------------------------------------------------------- */
/* Setup field definitions. The first nDimCount will be the */
/* dimension attributes, and after that comes the actual target */
/* array. */
/* -------------------------------------------------------------------- */
for( iDim = 0; iDim < nDimCount; iDim++ )
{
OGRFieldDefn oField( paoDimensions[iDim].pszDimName, OFTInteger );
if( EQUAL(oField.GetNameRef(), poTargetArray->name().c_str()) )
oField.SetName(CPLSPrintf("%s_i",paoDimensions[iDim].pszDimName));
if( paoDimensions[iDim].poMap != NULL )
{
switch( paoDimensions[iDim].poMap->var()->type() )
{
case dods_byte_c:
case dods_int16_c:
case dods_uint16_c:
case dods_int32_c:
case dods_uint32_c:
oField.SetType( OFTInteger );
break;
case dods_float32_c:
case dods_float64_c:
oField.SetType( OFTReal );
break;
case dods_str_c:
case dods_url_c:
oField.SetType( OFTString );
break;
default:
// Ignore
break;
}
}
poFeatureDefn->AddFieldDefn( &oField );
}
/* -------------------------------------------------------------------- */
/* Setup the array attributes themselves. */
/* -------------------------------------------------------------------- */
int iArray;
for( iArray=0; iArray < nArrayRefCount; iArray++ )
{
OGRDODSArrayRef *poRef = paoArrayRefs + iArray;
OGRFieldDefn oArrayField( poRef->poArray->name().c_str(), OFTInteger );
switch( poRef->poArray->var()->type() )
{
case dods_byte_c:
case dods_int16_c:
case dods_uint16_c:
case dods_int32_c:
case dods_uint32_c:
oArrayField.SetType( OFTInteger );
break;
case dods_float32_c:
case dods_float64_c:
oArrayField.SetType( OFTReal );
break;
case dods_str_c:
case dods_url_c:
oArrayField.SetType( OFTString );
break;
default:
// Ignore
break;
}
poFeatureDefn->AddFieldDefn( &oArrayField );
poRef->iFieldIndex = poFeatureDefn->GetFieldCount() - 1;
}
/* -------------------------------------------------------------------- */
/* X/Y/Z fields. */
/* -------------------------------------------------------------------- */
if( poOGRLayerInfo != NULL )
{
AttrTable *poField = poOGRLayerInfo->find_container("x_field");
if( poField != NULL )
{
oXField.Initialize( poField );
oXField.iFieldIndex =
poFeatureDefn->GetFieldIndex( oXField.pszFieldName );
}
poField = poOGRLayerInfo->find_container("y_field");
if( poField != NULL )
{
oYField.Initialize( poField );
oYField.iFieldIndex =
poFeatureDefn->GetFieldIndex( oYField.pszFieldName );
}
poField = poOGRLayerInfo->find_container("z_field");
if( poField != NULL )
{
oZField.Initialize( poField );
oZField.iFieldIndex =
poFeatureDefn->GetFieldIndex( oZField.pszFieldName );
}
}
/* -------------------------------------------------------------------- */
/* If we have no layerinfo, then check if there are obvious x/y */
/* fields. */
/* -------------------------------------------------------------------- */
else
{
if( poFeatureDefn->GetFieldIndex( "lat" ) != -1
&& poFeatureDefn->GetFieldIndex( "lon" ) != -1 )
{
oXField.Initialize( "lon", "dds" );
oXField.iFieldIndex = poFeatureDefn->GetFieldIndex( "lon" );
oYField.Initialize( "lat", "dds" );
oYField.iFieldIndex = poFeatureDefn->GetFieldIndex( "lat" );
}
else if( poFeatureDefn->GetFieldIndex( "latitude" ) != -1
&& poFeatureDefn->GetFieldIndex( "longitude" ) != -1 )
{
oXField.Initialize( "longitude", "dds" );
oXField.iFieldIndex = poFeatureDefn->GetFieldIndex( "longitude" );
oYField.Initialize( "latitude", "dds" );
oYField.iFieldIndex = poFeatureDefn->GetFieldIndex( "latitude" );
}
}
/* -------------------------------------------------------------------- */
/* Set the layer geometry type if we have point inputs. */
/* -------------------------------------------------------------------- */
if( oZField.iFieldIndex >= 0 )
poFeatureDefn->SetGeomType( wkbPoint25D );
else if( oXField.iFieldIndex >= 0 && oYField.iFieldIndex >= 0 )
poFeatureDefn->SetGeomType( wkbPoint );
else
poFeatureDefn->SetGeomType( wkbNone );
}
int OGRDODSDataSource::Open( const char * pszNewName )
{
CPLAssert( nLayers == 0 );
pszName = CPLStrdup( pszNewName );
/* -------------------------------------------------------------------- */
/* Parse the URL into a base url, projection and constraint */
/* expression. */
/* -------------------------------------------------------------------- */
char *pszWrkURL = CPLStrdup( pszNewName + 5 );
char *pszFound;
pszFound = strstr(pszWrkURL,"&");
if( pszFound )
{
oConstraints = pszFound;
*pszFound = '\0';
}
pszFound = strstr(pszWrkURL,"?");
if( pszFound )
{
oProjection = pszFound+1;
*pszFound = '\0';
}
// Trim common requests.
int nLen = strlen(pszWrkURL);
if( strcmp(pszWrkURL+nLen-4,".das") == 0 )
pszWrkURL[nLen-4] = '\0';
else if( strcmp(pszWrkURL+nLen-4,".dds") == 0 )
pszWrkURL[nLen-4] = '\0';
else if( strcmp(pszWrkURL+nLen-4,".asc") == 0 )
pszWrkURL[nLen-4] = '\0';
else if( strcmp(pszWrkURL+nLen-5,".dods") == 0 )
pszWrkURL[nLen-5] = '\0';
else if( strcmp(pszWrkURL+nLen-5,".html") == 0 )
pszWrkURL[nLen-5] = '\0';
oBaseURL = pszWrkURL;
CPLFree( pszWrkURL );
/* -------------------------------------------------------------------- */
/* Do we want to override the .dodsrc file setting? Only do */
/* the putenv() if there isn't already a DODS_CONF in the */
/* environment. */
/* -------------------------------------------------------------------- */
if( CPLGetConfigOption( "DODS_CONF", NULL ) != NULL
&& getenv("DODS_CONF") == NULL )
{
static char szDODS_CONF[1000];
sprintf( szDODS_CONF, "DODS_CONF=%.980s",
CPLGetConfigOption( "DODS_CONF", "" ) );
putenv( szDODS_CONF );
}
/* -------------------------------------------------------------------- */
/* If we have a overridding AIS file location, apply it now. */
/* -------------------------------------------------------------------- */
if( CPLGetConfigOption( "DODS_AIS_FILE", NULL ) != NULL )
{
string oAISFile = CPLGetConfigOption( "DODS_AIS_FILE", "" );
RCReader::instance()->set_ais_database( oAISFile );
}
/* -------------------------------------------------------------------- */
/* Connect to the server. */
/* -------------------------------------------------------------------- */
string version;
try
{
poConnection = new AISConnect( oBaseURL );
version = poConnection->request_version();
}
catch (Error &e)
{
CPLError(CE_Failure, CPLE_OpenFailed,
"%s", e.get_error_message().c_str() );
return FALSE;
}
/* -------------------------------------------------------------------- */
/* We presume we only work with version 3 servers. */
/* -------------------------------------------------------------------- */
if (version.empty() || version.find("/3.") == string::npos)
{
CPLError( CE_Warning, CPLE_AppDefined,
"I connected to the URL but could not get a DAP 3.x version string\n"
"from the server. I will continue to connect but access may fail.");
}
/* -------------------------------------------------------------------- */
/* Fetch the DAS and DDS info about the server. */
/* -------------------------------------------------------------------- */
try
{
poConnection->request_das( oDAS );
poConnection->request_dds( *poDDS, oProjection + oConstraints );
}
catch (Error &e)
{
CPLError(CE_Failure, CPLE_AppDefined,
"Error fetching DAS or DDS:\n%s",
e.get_error_message().c_str() );
return FALSE;
}
/* -------------------------------------------------------------------- */
/* Do we have any ogr_layer_info attributes in the DAS? If so, */
/* use them to define the layers. */
/* -------------------------------------------------------------------- */
AttrTable::Attr_iter dv_i;
#ifdef LIBDAP_39
AttrTable* poTable = oDAS.container();
if (poTable == NULL)
{
CPLError(CE_Failure, CPLE_AppDefined, "Cannot get container");
return FALSE;
}
#else
AttrTable* poTable = &oDAS;
#endif
for( dv_i = poTable->attr_begin(); dv_i != poTable->attr_end(); dv_i++ )
{
if( EQUALN(poTable->get_name(dv_i).c_str(),"ogr_layer_info",14)
&& poTable->is_container( dv_i ) )
{
AttrTable *poAttr = poTable->get_attr_table( dv_i );
string target_container = poAttr->get_attr( "target_container" );
BaseType *poVar = poDDS->var( target_container.c_str() );
if( poVar == NULL )
{
CPLError( CE_Warning, CPLE_AppDefined,
"Unable to find variable '%s' named in\n"
"ogr_layer_info.target_container, skipping.",
target_container.c_str() );
continue;
}
if( poVar->type() == dods_sequence_c )
AddLayer(
new OGRDODSSequenceLayer(this,
target_container.c_str(),
poAttr) );
else if( poVar->type() == dods_grid_c
|| poVar->type() == dods_array_c )
AddLayer( new OGRDODSGridLayer(this,target_container.c_str(),
poAttr) );
}
}
/* -------------------------------------------------------------------- */
/* Walk through the DODS variables looking for easily targetted */
/* ones. Eventually this will need to be driven by the AIS info. */
/* -------------------------------------------------------------------- */
if( nLayers == 0 )
{
DDS::Vars_iter v_i;
for( v_i = poDDS->var_begin(); v_i != poDDS->var_end(); v_i++ )
{
BaseType *poVar = *v_i;
if( poVar->type() == dods_sequence_c )
AddLayer( new OGRDODSSequenceLayer(this,poVar->name().c_str(),
NULL) );
else if( poVar->type() == dods_grid_c
|| poVar->type() == dods_array_c )
AddLayer( new OGRDODSGridLayer(this,poVar->name().c_str(),
NULL) );
}
}
return TRUE;
}
void KernelSymbolWriter::write()
{
operationStarted();
// Update the time stamp if the symbols have changed
if (_factory->changeClock() != _specs->createdChangeClock) {
_specs->created = QDateTime::currentDateTime();
_specs->createdChangeClock = _factory->changeClock();
}
// Disable compression by default
qint16 flags = 0; // kSym::flagCompressed;
// First, write the header information to the uncompressed device
KernelSymbolStream out(_to);
// out.setKSymVersion(kSym::VERSION_11);
#ifdef WRITE_ASCII_FILE
QFile debugOutFile("/tmp/insight.log");
debugOutFile.open(QIODevice::WriteOnly);
QTextStream dout(&debugOutFile);
#endif
// Write the file header in the following format:
// 1. (qint32) magic number
// 2. (qint16) file version number
// 3. (qint16) flags (currently unused)
// 4. (qint32) Qt's serialization format version (see QDataStream::Version)
out << (qint32) kSym::fileMagic
<< (qint16) out.kSymVersion()
<< (qint16) flags
<< (qint32) out.version();
#ifdef WRITE_ASCII_FILE
dout << QString::fromAscii((char*)(&kSym::fileMagic), sizeof(kSym::fileMagic))
<< " " << kSym::fileVersion << " 0x" << hex << flags
<< dec << " " << out.version() << endl;
#endif
// Write all information from SymFactory in the following format:
// 1. (MemSpecs) data of _specs
// 2.a (qint32) number of compile units
// 2.b (CompileUnit) data of 1st compile unit
// 2.c (CompileUnit) data of 2nd compile unit
// 2.d ...
// 3.a (qint32) number of types
// 3.b (qint32) type (RealType casted to qint32)
// 3.c (subclass of BaseType) data of type
// 3.d (qint32) type (RealType casted to qint32)
// 3.e (subclass of BaseType) data of type
// 3.f ...
// 4.a (qint32) number of id-mappings for types
// 4.b (qint32) 1st source id
// 4.c (qint32) 1st target id
// 4.d (qint32) 2nd source id
// 4.e (qint32) 2nd target id
// 4.f ...
// 5.a (qint32) number of variables
// 5.b (Variable) data of variable
// 5.c (Variable) data of variable
// 5.d ...
// 6.a (qint32) number of ref. types with alternative types
// 6.b (qint32) 1st id of ref. type with alternatives
// 6.c (qint32) number of type alternatives
// 6.d (AltRefType) 1st alternative
// 6.e (AltRefType) 2nd alternative
// 6.f (AltRefType) ...
// 6.g (qint32) 2st id of ref. type with alternatives
// 6.h (qint32) number of type alternatives
// 6.i (AltRefType) 1st alternative
// 6.j (AltRefType) 2nd alternative
// 6.k (AltRefType) ...
// 6.l ...
// 7.a (qint32) number of struct members with alternative types
// 7.b (qint32) 1st id of struct member with alternatives
// 7.c (qint32) id of belonging struct
// 7.d (qint32) number of type alternatives
// 7.e (AltRefType) 1st alternative
// 7.f (AltRefType) 2nd alternative
// 7.g (AltRefType) ...
// 7.h (qint32) 2st id of struct member with alternatives
// 7.i (qint32) id of belonging struct
// 7.j (qint32) number of type alternatives
// 7.k (AltRefType) 1st alternative
// 7.l (AltRefType) 2nd alternative
// 7.m (AltRefType) ...
// 7.l ...
// 8.a (qint32) number of variable with alternative types
// 8.b (qint32) 1st id of variable with alternatives
// 8.c (qint32) number of type alternatives
// 8.d (AltRefType) 1st alternative
// 8.e (AltRefType) 2nd alternative
// 8.f (AltRefType) ...
// 8.g (qint32) 2st id of variable with alternatives
// 8.h (qint32) number of type alternatives
// 8.i (AltRefType) 1st alternative
// 8.j (AltRefType) 2nd alternative
// 8.k (AltRefType) ...
// 8.l ...
try {
QSet<qint32> written_types;
// Write the memory specifications
out << *_specs;
#ifdef WRITE_ASCII_FILE
dout << endl << "# Memory specifications" << endl
<< _specs->toString();
#endif
// Write list of compile units
out << (qint32) _factory->sources().size();
#ifdef WRITE_ASCII_FILE
dout << endl << "# Compile units" << endl
<< _factory->sources().size() << endl;
#endif
CompileUnitIntHash::const_iterator cu_it = _factory->sources().constBegin();
while (cu_it != _factory->sources().constEnd()) {
const CompileUnit* c = cu_it.value();
out << *c;
#ifdef WRITE_ASCII_FILE
dout << "0x" << hex << c->id() << " " << c->name() << endl;
#endif
++cu_it;
checkOperationProgress();
}
// Write list of types
const int types_to_write = _factory->types().size();
out << (qint32) types_to_write;
#ifdef WRITE_ASCII_FILE
dout << endl << "# Types" << endl
<< dec << types_to_write << endl;
#endif
// Make three rounds: first write elementary types, then the
// simple referencing types, finally the structs and unions
for (int round = 0; round < 3; ++round) {
int mask = ElementaryTypes;
switch (round) {
case 1: mask = ReferencingTypes & ~StructOrUnion; break;
case 2: mask = StructOrUnion; break;
}
for (int i = 0; i < _factory->types().size(); i++) {
BaseType* t = _factory->types().at(i);
if (t->type() & mask) {
out << (qint32) t->type();
out << *t;
#ifdef WRITE_ASCII_FILE
dout << "0x" << hex << t->id() << " "
<< realTypeToStr(t->type()) << " "
<< t->name();
RefBaseType* rbt = dynamic_cast<RefBaseType*>(t);
if (rbt)
dout << ", refTypeId = 0x" << rbt->refTypeId();
dout << endl;
#endif
// Remember which types we have written out
written_types.insert(t->id());
}
checkOperationProgress();
}
}
assert(_factory->types().size() == written_types.size());
assert(types_to_write == written_types.size());
// Write list of missing types by ID
const int ids_to_write =
_factory->typesById().size() - _factory->types().size();
out << (qint32)ids_to_write;
#ifdef WRITE_ASCII_FILE
dout << endl << "# Further type relations" << endl
<< dec << ids_to_write
<< endl;
#endif
BaseTypeIntHash::const_iterator bt_id_it = _factory->typesById().constBegin();
int written = 0;
while (bt_id_it != _factory->typesById().constEnd()) {
if (!written_types.contains(bt_id_it.key())) {
out << (qint32) bt_id_it.key() << (qint32) bt_id_it.value()->id();
#ifdef WRITE_ASCII_FILE
dout << hex << "0x" << bt_id_it.key() << " -> 0x"
<< bt_id_it.value()->id() << endl;
#endif
++written;
}
++bt_id_it;
checkOperationProgress();
}
assert(written == ids_to_write);
assert(written_types.size() + written == _factory->typesById().size());
// Write list of variables
out << (qint32) _factory->vars().size();
#ifdef WRITE_ASCII_FILE
dout << endl << "# List of variables" << endl
<< dec << _factory->vars().size() << endl;
#endif
for (int i = 0; i < _factory->vars().size(); i++) {
out << *_factory->vars().at(i);
#ifdef WRITE_ASCII_FILE
dout << hex << "0x" << _factory->vars().at(i)->id() << " "
<< _factory->vars().at(i)->name() << ", refTypeId = 0x"
<< _factory->vars().at(i)->refTypeId() << endl;
#endif
checkOperationProgress();
}
// Find referencing types with alternatives
QList<RefBaseType*> refTypesWithAlt;
MemberList membersWithAlt;
for (int i = 0; i < _factory->types().count(); ++i) {
BaseType* t = _factory->types().at(i);
// Non-structure types
if (t->type() & ReferencingTypes & ~StructOrUnion) {
RefBaseType* rbt = dynamic_cast<RefBaseType*>(t);
if (rbt->altRefTypeCount() > 0)
refTypesWithAlt.append(rbt);
}
// Structure types
else if (t->type() & StructOrUnion) {
Structured* s = dynamic_cast<Structured*>(t);
for (int j = 0; j < s->members().count(); ++j) {
StructuredMember* m = s->members().at(j);
if (m->altRefTypeCount() > 0)
membersWithAlt.append(m);
}
}
checkOperationProgress();
}
// Find variables with type alternatives
VariableList varsWithAlt;
for (int i = 0; i < _factory->vars().size(); i++) {
Variable* v = _factory->vars().at(i);
if (v->altRefTypeCount() > 0)
varsWithAlt.append(v);
checkOperationProgress();
}
// Write list of types with alternative types
out << (qint32) refTypesWithAlt.size();
#ifdef WRITE_ASCII_FILE
dout << endl << "# List of types with alternative types" << endl
<< dec << refTypesWithAlt.size() << endl;
#endif
for (int i = 0; i < refTypesWithAlt.size(); ++i) {
RefBaseType* rbt = refTypesWithAlt.at(i);
out << (qint32) rbt->id();
rbt->writeAltRefTypesTo(out);
#ifdef WRITE_ASCII_FILE
dout << hex << "0x" << rbt->id() << " "
<< dec << rbt->altRefTypeCount()
<< endl;
#endif
checkOperationProgress();
}
// Write list of struct members with alternative types
out << (qint32) (refTypesWithAlt.size() + membersWithAlt.size());
#ifdef WRITE_ASCII_FILE
dout << endl << "# List of struct members with alternative types" << endl
<< dec << membersWithAlt.size() << endl;
#endif
for (int i = 0; i < membersWithAlt.size(); ++i) {
StructuredMember* m = membersWithAlt.at(i);
out << (qint32) m->id()
<< (qint32) m->belongsTo()->id();
m->writeAltRefTypesTo(out);
#ifdef WRITE_ASCII_FILE
dout << hex << "0x" << m->id() << " "
<< hex << "0x" << m->belongsTo()->id() << " "
<< dec << m->altRefTypeCount()
<< endl;
#endif
checkOperationProgress();
}
// Write list of variables with alternative types
out << (qint32) varsWithAlt.size();
#ifdef WRITE_ASCII_FILE
dout << endl << "# List of variables with alternative types" << endl
<< dec << varsWithAlt.size() << endl;
#endif
for (int i = 0; i < varsWithAlt.size(); ++i) {
Variable* v = varsWithAlt.at(i);
out << (qint32) v->id();
v->writeAltRefTypesTo(out);
#ifdef WRITE_ASCII_FILE
dout << hex << "0x" << v->id() << " "
<< dec << v->altRefTypeCount()
<< endl;
#endif
checkOperationProgress();
}
// Since version 17: Write file names containing the orig. symbols
if (out.kSymVersion() >= kSym::VERSION_17)
out <<_factory->origSymFiles();
}
catch (...) {
// Exceptional cleanup
operationStopped();
Console::out() << endl;
throw; // Re-throw exception
}
operationStopped();
QString s("\rReading symbols finished");
if (!_to->isSequential())
s += QString(" (%1 read)").arg(bytesToString(_to->pos()));
s += ".";
shellOut(s, true);
}