protected MapDataElement ReadData(decimal Latitude, decimal Longitude, ReadPointValueDelegate ReadPointValue, decimal DataCoordinateStep)
{
this.readPointValue = ReadPointValue;
//find the lower-left cornerpoint of a quad
decimal LatLow = Latitude.FloorWithSignificance(DataCoordinateStep);
decimal LongLow = Longitude.FloorWithSignificance(DataCoordinateStep);
//SeismicGroundMotionDataPoint pointFound = this.FindPointInDataFileRecursive(LatLow, LongLow, DataStartIndex+1, DataEndIndex);
IMultipleValueDataPoint2D foundPoint = this.FindPointInDataFileRecursive(LatLow, LongLow, DataStartIndex + 1, DataEndIndex);
MapDataElement de = new MapDataElement(Latitude, Longitude);
//add neigboring points for interpolation
#region Neigbouring points for interpolation
#region Special cases
//Corner points
if (foundPoint.DataArrayIndex == DataStartIndex + 1 || foundPoint.DataArrayIndex == NumberOfColumns + 1 || foundPoint.DataArrayIndex == DataEndIndex || foundPoint.DataArrayIndex == DataEndIndex - NumberOfColumns + 1)
{
de.LowerLeftPoint = foundPoint;
de.LowerRightPoint = foundPoint.Clone();
de.LowerRightPoint.Longitude = de.LowerRightPoint.Longitude + CoordinateIncrement;
de.UpperLeftPoint = foundPoint.Clone();
de.UpperLeftPoint.Latitude = de.UpperLeftPoint.Latitude + CoordinateIncrement;
de.UpperRightPoint = foundPoint.Clone();
de.UpperRightPoint.Latitude = de.UpperRightPoint.Latitude + CoordinateIncrement;
de.UpperRightPoint.Longitude = de.UpperRightPoint.Longitude + CoordinateIncrement;
return(de);
}
//bottom row- do nothing
//top row
if (foundPoint.DataArrayIndex < (NumberOfColumns + 1))
{
de.LowerLeftPoint = foundPoint;
IMultipleValueDataPoint2D pNextLongitude = readPointValue(foundPoint.DataArrayIndex + 1);
de.LowerRightPoint = pNextLongitude;
IMultipleValueDataPoint2D NextLatitudePoint = foundPoint.Clone();
NextLatitudePoint.Latitude = foundPoint.Latitude + CoordinateIncrement;
de.UpperLeftPoint = NextLatitudePoint;
IMultipleValueDataPoint2D NextLatitudeAndLongitudePoint = pNextLongitude.Clone();
NextLatitudeAndLongitudePoint.Latitude = foundPoint.Latitude + CoordinateIncrement;
NextLatitudeAndLongitudePoint.Longitude = foundPoint.Longitude + CoordinateIncrement;
de.UpperRightPoint = NextLatitudeAndLongitudePoint;
return(de);
}
//first column -do nothing
//last column
if ((foundPoint.DataArrayIndex - DataStartIndex) % (NumberOfColumns) == 0 && foundPoint.DataArrayIndex > 2)
{
de.LowerLeftPoint = foundPoint;
IMultipleValueDataPoint2D pNextLongitude = foundPoint.Clone();
pNextLongitude.Longitude = foundPoint.Longitude + CoordinateIncrement;
de.LowerRightPoint = pNextLongitude;
IMultipleValueDataPoint2D pNextLatitude = readPointValue(foundPoint.DataArrayIndex - NumberOfColumns);
de.UpperLeftPoint = pNextLatitude;
IMultipleValueDataPoint2D NextLatitudeAndLongitudePoint = pNextLatitude.Clone();
NextLatitudeAndLongitudePoint.Latitude = pNextLatitude.Latitude;
NextLatitudeAndLongitudePoint.Longitude = pNextLongitude.Longitude;
de.UpperRightPoint = NextLatitudeAndLongitudePoint;
return(de);
}
#endregion
#region Typical case
de.LowerLeftPoint = foundPoint;
de.LowerRightPoint = readPointValue(foundPoint.DataArrayIndex + 1); //NextLongitudePoint;
de.UpperLeftPoint = readPointValue(foundPoint.DataArrayIndex - NumberOfColumns); //NextLatitudePoint;
de.UpperRightPoint = readPointValue(foundPoint.DataArrayIndex - NumberOfColumns + 1); //NextLatitudeAndLongitudePoint
#endregion
#endregion
//Interpolation is done by MapDataElement
//if (pointFound!=null)
//{
// decimal SS = GetInterpolatedValue(pointFound, Latitude, Longitude, GroundMotionParameterType.SS);
// decimal S1 = GetInterpolatedValue(pointFound, Latitude, Longitude, GroundMotionParameterType.S1);
// return new SeismicGroundMotionDataPoint(Latitude, Longitude, SS, S1);
//}
//return new SeismicGroundMotionDataPoint();
return(de);
}
protected MapDataElement ReadData(decimal Latitude, decimal Longitude, ReadPointValueDelegate ReadPointValue, decimal DataCoordinateStep)
{
this.readPointValue = ReadPointValue;
//find the lower-left cornerpoint of a quad
decimal LatLow = Latitude.FloorWithSignificance(DataCoordinateStep);
decimal LongLow = Longitude.FloorWithSignificance(DataCoordinateStep);
//SeismicGroundMotionDataPoint pointFound = this.FindPointInDataFileRecursive(LatLow, LongLow, DataStartIndex+1, DataEndIndex);
IMultipleValueDataPoint2D foundPoint = this.FindPointInDataFileRecursive(LatLow, LongLow, DataStartIndex + 1, DataEndIndex);
MapDataElement de = new MapDataElement(Latitude, Longitude);
//add neigboring points for interpolation
#region Neigbouring points for interpolation
#region Special cases
//Corner points
if (foundPoint.DataArrayIndex == DataStartIndex + 1 || foundPoint.DataArrayIndex == NumberOfColumns + 1 || foundPoint.DataArrayIndex == DataEndIndex || foundPoint.DataArrayIndex == DataEndIndex - NumberOfColumns + 1)
{
de.LowerLeftPoint = foundPoint;
de.LowerRightPoint = foundPoint.Clone();
de.LowerRightPoint.Longitude = de.LowerRightPoint.Longitude + CoordinateIncrement;
de.UpperLeftPoint = foundPoint.Clone();
de.UpperLeftPoint.Latitude = de.UpperLeftPoint.Latitude + CoordinateIncrement;
de.UpperRightPoint = foundPoint.Clone();
de.UpperRightPoint.Latitude = de.UpperRightPoint.Latitude + CoordinateIncrement;
de.UpperRightPoint.Longitude = de.UpperRightPoint.Longitude + CoordinateIncrement;
return de;
}
//bottom row- do nothing
//top row
if (foundPoint.DataArrayIndex < (NumberOfColumns + 1))
{
de.LowerLeftPoint = foundPoint;
IMultipleValueDataPoint2D pNextLongitude = readPointValue(foundPoint.DataArrayIndex + 1);
de.LowerRightPoint = pNextLongitude;
IMultipleValueDataPoint2D NextLatitudePoint = foundPoint.Clone();
NextLatitudePoint.Latitude = foundPoint.Latitude + CoordinateIncrement;
de.UpperLeftPoint = NextLatitudePoint;
IMultipleValueDataPoint2D NextLatitudeAndLongitudePoint = pNextLongitude.Clone();
NextLatitudeAndLongitudePoint.Latitude = foundPoint.Latitude + CoordinateIncrement;
NextLatitudeAndLongitudePoint.Longitude = foundPoint.Longitude + CoordinateIncrement;
de.UpperRightPoint = NextLatitudeAndLongitudePoint;
return de;
}
//first column -do nothing
//last column
if ((foundPoint.DataArrayIndex - DataStartIndex) % (NumberOfColumns) == 0 && foundPoint.DataArrayIndex > 2)
{
de.LowerLeftPoint = foundPoint;
IMultipleValueDataPoint2D pNextLongitude = foundPoint.Clone();
pNextLongitude.Longitude = foundPoint.Longitude + CoordinateIncrement;
de.LowerRightPoint = pNextLongitude;
IMultipleValueDataPoint2D pNextLatitude = readPointValue(foundPoint.DataArrayIndex - NumberOfColumns);
de.UpperLeftPoint = pNextLatitude;
IMultipleValueDataPoint2D NextLatitudeAndLongitudePoint = pNextLatitude.Clone();
NextLatitudeAndLongitudePoint.Latitude = pNextLatitude.Latitude;
NextLatitudeAndLongitudePoint.Longitude = pNextLongitude.Longitude;
de.UpperRightPoint = NextLatitudeAndLongitudePoint;
return de;
}
#endregion
#region Typical case
de.LowerLeftPoint = foundPoint;
de.LowerRightPoint = readPointValue(foundPoint.DataArrayIndex + 1); //NextLongitudePoint;
de.UpperLeftPoint = readPointValue(foundPoint.DataArrayIndex - NumberOfColumns); //NextLatitudePoint;
de.UpperRightPoint = readPointValue(foundPoint.DataArrayIndex - NumberOfColumns + 1); //NextLatitudeAndLongitudePoint
#endregion
#endregion
//Interpolation is done by MapDataElement
//if (pointFound!=null)
//{
// decimal SS = GetInterpolatedValue(pointFound, Latitude, Longitude, GroundMotionParameterType.SS);
// decimal S1 = GetInterpolatedValue(pointFound, Latitude, Longitude, GroundMotionParameterType.S1);
// return new SeismicGroundMotionDataPoint(Latitude, Longitude, SS, S1);
//}
//return new SeismicGroundMotionDataPoint();
return de;
}