STRBaseLaneStandard :: STRBaseLaneStandard (const RGString& xmlString) : BaseLaneStandard (xmlString),
mCharacteristicArray (NULL), mDifferenceArray (NULL), mUnnormalizedDifferenceArray(0),
mNormsLeft (NULL), mNormsRight (NULL), mRelativeSizes (NULL),
mSize (0), mDifferenceSize (0), mOmissionSize (0), mOmissionArray (NULL), mNumberOfLargePeaks (0),
mLargeCharacteristicArray (NULL), mLargeDifferenceArray (NULL), mNumberOfLargeDifferences (0),
mFirstIntervalFraction (-1.0), mSmallestIntervalFraction (-1.0), mMaxRelativeHeight (-1),
mCorrelationAcceptanceThreshold (0.993), mCorrelationAutoAcceptanceThreshold (0.9975)
//
//, mLastHalfSize (0), mLastHalfCharacteristicArray (NULL), mLastHalfUnnormalizedDifferenceArray (NULL),
//mLastHalfRelativeSizes (NULL), mLastHalfDifferenceSize (0), mLastHalfNormsLeft (NULL), mLastHalfNormsRight (NULL)
{
RGString TextInput (xmlString);
RGBracketStringSearch CharacteristicsToken ("<Characteristics>", "</Characteristics>", TextInput);
RGString CharacteristicsString;
RGString HeightString;
RGBracketStringSearch HeightToken ("<RelativeHeights>", "</RelativeHeights>", TextInput);
RGBracketStringSearch OmissionSearch ("<Omissions>", "</Omissions>", TextInput);
RGBracketStringSearch AcceptanceThresholdSearch ("<CorrelationAcceptanceThreshold>", "</CorrelationAcceptanceThreshold>", TextInput);
RGBracketStringSearch AutoAcceptanceThresholdSearch ("<CorrelationAutoAcceptanceThreshold>", "</CorrelationAutoAcceptanceThreshold>", TextInput);
RGString OmissionString;
size_t endPosition;
RGString token;
RGString delim;
RGString numberString;
double minChar;
double maxChar;
RGDList tokens;
double value;
RGPDouble* valuePtr;
double previous;
double next;
double cumulative;
int i = 0;
double temp;
int k;
int j;
if (!CharacteristicsToken.FindNextBracketedString (0, endPosition, CharacteristicsString)) {
Msg << "Internal Lane Standard named " << Name << " could not find Characteristics token from string\n";
Msg << TextInput << "\n";
Valid = FALSE;
}
else {
RGStringTokenizer data (CharacteristicsString);
data.AddDelimiter (' ');
data.AddDelimiter (" ");
data.AddDelimiter (" ");
data.AddDelimiter (" ");
data.AddDelimiter (" ");
data.AddRemoveItem (',');
data.AddDelimiter ('\t');
data.AddDelimiter ('\n');
while (data.NextToken (token, delim)) {
if (token.Length () > 0) {
value = token.ConvertToDouble ();
valuePtr = new RGPDouble (value);
if (tokens.Entries () == 0) {
minChar = maxChar = value;
tokens.Append (valuePtr);
}
else if (value < minChar) {
minChar = value;
tokens.Prepend (valuePtr);
}
else if (value > maxChar) {
maxChar = value;
tokens.Append (valuePtr);
}
else
tokens.Insert (valuePtr);
}
}
mSize = (int)tokens.Entries ();
NCharacteristics = mSize;
mCharacteristicArray = new double [mSize];
mDifferenceSize = mSize - 1;
mDifferenceArray = new double [mDifferenceSize];
mUnnormalizedDifferenceArray = new double [mDifferenceSize];
mNormsLeft = new double [mDifferenceSize];
mNormsRight = new double [mDifferenceSize];
//mLastHalfSize = mSize / 2;
//if (2 * mLastHalfSize < mSize)
// mLastHalfSize += 1;
//mLastHalfDifferenceSize = mLastHalfSize - 1;
//mLastHalfCharacteristicArray = new double [mLastHalfSize];
//mLastHalfUnnormalizedDifferenceArray = new double [mLastHalfDifferenceSize];
//mLastHalfNormsLeft = new double [mLastHalfDifferenceSize];
//mLastHalfNormsRight = new double [mLastHalfDifferenceSize];
while (valuePtr = (RGPDouble*)tokens.GetFirst ()) {
mCharacteristicArray [i] = valuePtr->GetDouble ();
delete valuePtr;
i++;
}
//int displacement = mSize - mLastHalfSize; // +0, +1, -1?????
//for (i=0; i<mLastHalfSize; i++)
// mLastHalfCharacteristicArray [i] = mCharacteristicArray [i + displacement];
previous = mCharacteristicArray [0];
for (i=0; i<mDifferenceSize; i++) {
next = mCharacteristicArray [i+1];
mUnnormalizedDifferenceArray [i] = mDifferenceArray [i] = next - previous;
previous = next;
}
//for (i=0; i<mLastHalfSize; i++)
// mLastHalfUnnormalizedDifferenceArray [i] = mUnnormalizedDifferenceArray [i + displacement];
double norm = 0.0;
for (i=0; i<mDifferenceSize; i++)
norm += mDifferenceArray [i] * mDifferenceArray [i];
norm = sqrt (norm);
if (norm > 0.0) {
for (i=0; i<mDifferenceSize; i++)
mDifferenceArray [i] /= norm;
}
cumulative = 0.0;
for (i=0; i<mDifferenceSize; i++) {
temp = mUnnormalizedDifferenceArray [i];
cumulative += temp * temp;
mNormsLeft [i] = cumulative;
}
cumulative = 0.0;
for (i=mDifferenceSize-1; i>=0; i--) {
temp = mUnnormalizedDifferenceArray [i];
cumulative += temp * temp;
mNormsRight [i] = cumulative;
}
//cumulative = 0.0;
//for (i=0; i<mLastHalfDifferenceSize; i++) {
// temp = mLastHalfUnnormalizedDifferenceArray [i];
// cumulative += temp * temp;
// mLastHalfNormsLeft [i] = cumulative;
//}
//cumulative = 0.0;
//for (i=mLastHalfDifferenceSize-1; i>=0; i--) {
// temp = mLastHalfUnnormalizedDifferenceArray [i];
// cumulative += temp * temp;
// mLastHalfNormsRight [i] = cumulative;
//}
MaxCharacteristic = maxChar;
MinCharacteristic = minChar;
double width = mCharacteristicArray [mSize - 1] - mCharacteristicArray [0];
if (width > 0.0)
mFirstIntervalFraction = mUnnormalizedDifferenceArray [0] / width;
else
mFirstIntervalFraction = 1.0;
double temp = mUnnormalizedDifferenceArray [0];
double value;
for (i=0; i<mDifferenceSize; i++) {
value = mUnnormalizedDifferenceArray [i];
if (value < temp)
temp = value;
}
if (width > 0.0)
mSmallestIntervalFraction = temp / width;
else
mSmallestIntervalFraction = 1.0;
if (!HeightToken.FindNextBracketedString (0, endPosition, HeightString)) {
mRelativeSizes = NULL;
// mLastHalfRelativeSizes = NULL;
}
else {
RGStringTokenizer data2 (HeightString);
data2.AddDelimiter (' ');
data2.AddDelimiter (" ");
data2.AddDelimiter (" ");
data2.AddDelimiter (" ");
data2.AddDelimiter (" ");
data2.AddRemoveItem (',');
data2.AddDelimiter ('\t');
data2.AddDelimiter ('\n');
mRelativeSizes = new int [mSize];
for (k=0; k<mSize; k++)
mRelativeSizes [k] = -1;
k = 0;
while (data2.NextToken (token, delim)) {
if (token.Length () > 0) {
token.ToUpper ();
char first = token.GetFirstCharacter ();
switch (first) {
case 'H':
mRelativeSizes [k] = 5;
break;
case 'L':
mRelativeSizes [k] = 1;
break;
case 'M':
if (token.Length () == 1)
mRelativeSizes [k] = 3;
else if (token.GetLastCharacter () == 'H')
mRelativeSizes [k] = 4;
else
mRelativeSizes [k] = 2;
break;
default:
mRelativeSizes [k] = -1;
}
}
k++;
if (k >= mSize)
break;
}
//mLastHalfRelativeSizes = new int [mLastHalfSize];
//for (i=0; i<mLastHalfSize; i++)
// mLastHalfRelativeSizes [i] = mRelativeSizes [i + displacement];
int maxRelativeSize = 0;
int relativeSize;
for (j=0; j<mSize; j++) {
relativeSize = mRelativeSizes [j];
if (relativeSize > maxRelativeSize)
maxRelativeSize = relativeSize;
}
int numberOfHighsInILS = 0;
mMaxRelativeHeight = maxRelativeSize;
for (j=0; j<mSize; j++) {
if (mRelativeSizes [j] == maxRelativeSize)
numberOfHighsInILS++;
}
mNumberOfLargePeaks = numberOfHighsInILS;
mNumberOfLargeDifferences = numberOfHighsInILS - 1;
mLargeCharacteristicArray = new double [numberOfHighsInILS];
if (mNumberOfLargeDifferences > 0)
mLargeDifferenceArray = new double [mNumberOfLargeDifferences];
k = 0;
for (j=0; j<mSize; j++) {
if (mRelativeSizes [j] == maxRelativeSize) {
mLargeCharacteristicArray [k] = mCharacteristicArray [j];
k++;
}
}
if (mLargeDifferenceArray != NULL) {
for (j=0; j<mNumberOfLargeDifferences; j++)
mLargeDifferenceArray [j] = mLargeCharacteristicArray [j + 1] - mLargeCharacteristicArray [j];
}
}
if (OmissionSearch.FindNextBracketedString (0, endPosition, OmissionString)) {
RGStringTokenizer omissionData (OmissionString);
omissionData.AddDelimiter (' ');
omissionData.AddDelimiter (" ");
omissionData.AddDelimiter (" ");
omissionData.AddDelimiter (" ");
omissionData.AddDelimiter (" ");
omissionData.AddRemoveItem (',');
omissionData.AddDelimiter ('\t');
omissionData.AddDelimiter ('\n');
while (omissionData.NextToken (token, delim)) {
if (token.Length () > 0) {
value = token.ConvertToDouble ();
valuePtr = new RGPDouble (value);
if (tokens.Entries () == 0) {
minChar = maxChar = value;
tokens.Append (valuePtr);
}
else if (value < minChar) {
minChar = value;
tokens.Prepend (valuePtr);
}
else if (value > maxChar) {
maxChar = value;
tokens.Append (valuePtr);
}
else
tokens.Insert (valuePtr);
}
}
mOmissionSize = (int)tokens.Entries ();
if (mOmissionSize == 0)
return;
mOmissionArray = new double [mOmissionSize];
i = 0;
while (valuePtr = (RGPDouble*)tokens.GetFirst ()) {
mOmissionArray [i] = valuePtr->GetDouble ();
delete valuePtr;
i++;
}
}
}
cout << "ILS Name = " << Name.GetData () << endl;
if (!AcceptanceThresholdSearch.FindNextBracketedString (0, endPosition, numberString))
cout << "\tCould not find ILS acceptance threshold...using " << mCorrelationAcceptanceThreshold << endl;
else {
mCorrelationAcceptanceThreshold = numberString.ConvertToDouble ();
cout << "\tFrom ILS settings...ILS acceptance threshold = " << mCorrelationAcceptanceThreshold << endl;
}
if (!AutoAcceptanceThresholdSearch.FindNextBracketedString (0, endPosition, numberString))
cout << "\tCould not find ILS auto acceptance threshold...using default value = " << mCorrelationAutoAcceptanceThreshold << endl << endl;
else {
mCorrelationAutoAcceptanceThreshold = numberString.ConvertToDouble ();
cout << "\tFrom ILS settings...ILS auto acceptance threshold = " << mCorrelationAutoAcceptanceThreshold << endl << endl;
}
}
RGPersistent* RGHashTable :: Insert (RGPersistent* p) {
RGDList* list = SelectList (p);
NumberOfEntries++;
return list->Insert (p);
}
