public DEMFileDefinition(string name, DEMFileType format, string extension, DEMFileRegistrationMode registration)
{
this.Name = name;
this.Type = format;
this.FileExtension = extension;
this.Registration = registration;
}
public DEMFileDefinition(DEMFileType format, DEMFileRegistrationMode registration)
{
this.Name = null;
this.Type = format;
this.FileExtension = null;
this.Registration = registration;
}
/// <summary>
/// Open specified file
/// </summary>
/// <param name="filePath">If path is rooted (full file name), the specified file will be openened,
/// otherwise the file path will be relative to <see cref="LocalDirectory"/></param>
/// <param name="fileFormat"><see cref="DEMFileType"/> enumeration indicating the file type</param>
/// <returns><see cref="IRasterFile"/> interface for accessing file contents</returns>
public IRasterFile OpenFile(string filePath, DEMFileType fileFormat)
{
if (!Path.IsPathRooted(filePath))
{
filePath = Path.Combine(_localDirectory, filePath);
}
switch (fileFormat)
{
case DEMFileType.GEOTIFF: return(new GeoTiff(filePath));
case DEMFileType.SRTM_HGT: return(new HGTFile(filePath));
case DEMFileType.ASCIIGrid: return(new ASCIIGridFile(filePath));
case DEMFileType.CF_NetCDF: return(new NetCdfFile(filePath));
default:
throw new NotImplementedException($"{fileFormat} file format not implemented.");
}
}
void TestEdges(string dataSetName, double lat, double lon
, string rasterSouthWestName, string rasterSouthEastName
, string rasterNorthWestName, string rasterNorthEastName)
{
DEMDataSet dataSet = DEMDataSet.RegisteredDatasets.FirstOrDefault(d => d.Name == dataSetName);
Assert.NotNull(dataSet);
DEMFileType fileType = dataSet.FileFormat.Type;
int rasterSize = dataSet.PointsPerDegree;
double amountx = (1d / rasterSize) / 4d;
double amounty = (1d / rasterSize) / 4d;
// Regenerates all metadata
//_rasterService.GenerateDirectoryMetadata(dataSet
// , force: true
// , deleteOnError: false
// , maxDegreeOfParallelism: 1);
_elevationService.DownloadMissingFiles(dataSet, lat, lon);
var tiles = _rasterService.GenerateReportForLocation(dataSet, lat, lon);
Assert.True(tiles.Count == 4);
Assert.Single(tiles, t => string.Equals(Path.GetFileName(t.LocalName), rasterSouthWestName, StringComparison.OrdinalIgnoreCase));
Assert.Single(tiles, t => string.Equals(Path.GetFileName(t.LocalName), rasterSouthEastName, StringComparison.OrdinalIgnoreCase));
Assert.Single(tiles, t => string.Equals(Path.GetFileName(t.LocalName), rasterNorthWestName, StringComparison.OrdinalIgnoreCase));
Assert.Single(tiles, t => string.Equals(Path.GetFileName(t.LocalName), rasterNorthEastName, StringComparison.OrdinalIgnoreCase));
if (dataSet.FileFormat.Registration == DEMFileRegistrationMode.Cell)
{
using (var rasterNW = _rasterService.OpenFile(tiles.First(t => string.Equals(rasterNorthWestName, Path.GetFileName(t.LocalName))).LocalName, fileType))
using (var rasterNE = _rasterService.OpenFile(tiles.First(t => string.Equals(rasterNorthEastName, Path.GetFileName(t.LocalName))).LocalName, fileType))
using (var rasterSW = _rasterService.OpenFile(tiles.First(t => string.Equals(rasterSouthWestName, Path.GetFileName(t.LocalName))).LocalName, fileType))
using (var rasterSE = _rasterService.OpenFile(tiles.First(t => string.Equals(rasterSouthEastName, Path.GetFileName(t.LocalName))).LocalName, fileType))
{
var elevNW = rasterNW.GetElevationAtPoint(rasterNW.ParseMetaData(dataSet.FileFormat), rasterSize - 1, rasterSize - 1);
var elevNE = rasterNE.GetElevationAtPoint(rasterNE.ParseMetaData(dataSet.FileFormat), 0, rasterSize - 1);
var elevSW = rasterSW.GetElevationAtPoint(rasterSW.ParseMetaData(dataSet.FileFormat), rasterSize - 1, 0);
var elevSE = rasterSE.GetElevationAtPoint(rasterSE.ParseMetaData(dataSet.FileFormat), 0, 0);
BilinearInterpolator interpolator = new BilinearInterpolator();
var elev0 = interpolator.Interpolate(elevSW, elevSE, elevNW, elevNE, 0.25, 0.25);
var apiElev0 = _elevationService.GetPointElevation(lat + amounty, lon - amountx, dataSet);
Assert.True((elev0 - apiElev0.Elevation.Value) < double.Epsilon);
var elev1 = interpolator.Interpolate(elevSW, elevSE, elevNW, elevNE, 0.75, 0.25);
var apiElev1 = _elevationService.GetPointElevation(lat + amounty, lon + amountx, dataSet);
Assert.True((elev1 - apiElev1.Elevation.Value) < double.Epsilon);
var elev2 = interpolator.Interpolate(elevSW, elevSE, elevNW, elevNE, 0.25, 0.75);
var apiElev2 = _elevationService.GetPointElevation(lat - amounty, lon - amountx, dataSet);
Assert.True((elev2 - apiElev2.Elevation.Value) < double.Epsilon);
var elev3 = interpolator.Interpolate(elevSW, elevSE, elevNW, elevNE, 0.75, 0.75);
var apiElev3 = _elevationService.GetPointElevation(lat - amounty, lon + amountx, dataSet);
Assert.True((elev3 - apiElev3.Elevation.Value) < double.Epsilon);
}
}
else
{
using (var rasterNW = _rasterService.OpenFile(tiles.First(t => string.Equals(rasterNorthWestName, Path.GetFileName(t.LocalName))).LocalName, fileType))
using (var rasterNE = _rasterService.OpenFile(tiles.First(t => string.Equals(rasterNorthEastName, Path.GetFileName(t.LocalName))).LocalName, fileType))
using (var rasterSW = _rasterService.OpenFile(tiles.First(t => string.Equals(rasterSouthWestName, Path.GetFileName(t.LocalName))).LocalName, fileType))
using (var rasterSE = _rasterService.OpenFile(tiles.First(t => string.Equals(rasterSouthEastName, Path.GetFileName(t.LocalName))).LocalName, fileType))
{
// Northen row, west to east
var elevN0 = rasterNW.GetElevationAtPoint(rasterNW.ParseMetaData(dataSet.FileFormat), rasterSize - 1, rasterSize - 1);
var elevN1 = rasterNW.GetElevationAtPoint(rasterNW.ParseMetaData(dataSet.FileFormat), rasterSize, rasterSize - 1);
var elevN2 = rasterNE.GetElevationAtPoint(rasterNE.ParseMetaData(dataSet.FileFormat), 1, rasterSize - 1);
// middle row, west to east
var elevM0 = rasterNW.GetElevationAtPoint(rasterNW.ParseMetaData(dataSet.FileFormat), rasterSize - 1, rasterSize);
var elevM1 = rasterNW.GetElevationAtPoint(rasterNW.ParseMetaData(dataSet.FileFormat), rasterSize, rasterSize);
var elevM2 = rasterNE.GetElevationAtPoint(rasterNE.ParseMetaData(dataSet.FileFormat), 1, rasterSize);
// Sourthen row, west to east
var elevS0 = rasterSW.GetElevationAtPoint(rasterSW.ParseMetaData(dataSet.FileFormat), rasterSize - 1, 1);
var elevS1 = rasterSW.GetElevationAtPoint(rasterSW.ParseMetaData(dataSet.FileFormat), rasterSize, 1);
var elevS2 = rasterSE.GetElevationAtPoint(rasterSE.ParseMetaData(dataSet.FileFormat), 1, 1);
BilinearInterpolator interpolator = new BilinearInterpolator();
var elev0 = interpolator.Interpolate(elevM0, elevM1, elevN0, elevN1, 0.75, 0.75);
var apiElev0 = _elevationService.GetPointElevation(lat + amounty, lon - amountx, dataSet);
Assert.True((elev0 - apiElev0.Elevation.Value) < double.Epsilon);
var elev1 = interpolator.Interpolate(elevM1, elevM2, elevN1, elevN2, 0.25, 0.75);
var apiElev1 = _elevationService.GetPointElevation(lat + amounty, lon + amountx, dataSet);
Assert.True((elev1 - apiElev1.Elevation.Value) < double.Epsilon);
var elev2 = interpolator.Interpolate(elevS0, elevS1, elevM0, elevM1, 0.75, 0.25);
var apiElev2 = _elevationService.GetPointElevation(lat - amounty, lon - amountx, dataSet);
Assert.True((elev2 - apiElev2.Elevation.Value) < double.Epsilon);
var elev3 = interpolator.Interpolate(elevS1, elevS2, elevM1, elevM2, 0.25, 0.25);
var apiElev3 = _elevationService.GetPointElevation(lat - amounty, lon + amountx, dataSet);
Assert.True((elev3 - apiElev3.Elevation.Value) < double.Epsilon);
}
}
}