private static readonly short VALUE_SIZE = 2;//Bytes
public static void ReadTiff(string filename, GpsLocation filterMin, GpsLocation filterMax, GpsLocation myLocation, int skipFactor, List <GpsLocation> eleData)
{
try
{
using (Tiff tiff = Tiff.Open(filename, "r"))
{
if (tiff == null)
{
throw new SystemException("File open error.");
}
var(width, height) = GetImageSize(tiff);
var(pixelSizeX, pixelSizeY) = GetPixelSize(tiff);
var startLocation = GetStartLocation(tiff);
var scanline = new byte[tiff.ScanlineSize()];
for (int i = 0; i < height; i++)
{
tiff.ReadScanline(scanline, i);
var latitude = startLocation.Latitude + (pixelSizeY * i);
if (scanline.Length / 2 != width)
{
throw new SystemException("Invalid GeoTiff data");
}
if (latitude >= filterMin.Latitude && latitude <= filterMax.Latitude)
{
for (var j = 0; j < scanline.Length / VALUE_SIZE; j++)
{
var longitude = startLocation.Longitude + (pixelSizeX * j);
byte loByte = scanline[j * VALUE_SIZE + 0];
byte hiByte = scanline[j * VALUE_SIZE + 1];
var alt = hiByte * 256 + loByte;
if (longitude >= filterMin.Longitude && longitude <= filterMax.Longitude)
{
if (i % skipFactor == 0 && j % skipFactor == 0)
{
var ep = new GpsLocation {
Altitude = alt, Latitude = latitude, Longitude = longitude
};
ep.QuickDistance(myLocation);
ep.QuickBearing(myLocation);
ep.GetVerticalViewAngle(myLocation);
eleData.Add(ep);
}
}
}
}
}
}
}
catch (Exception ex)
{
throw new SystemException("GeoTiff parsing error", ex);
}
}
private static readonly double MAX_LON = Dg2Rad(180d); // PI
public GeoPoint(double latitude, double longitude, double alt, GpsLocation myLocation = null)
{
Latitude = latitude;
Longitude = longitude;
Altitude = (UInt32)alt;
if (myLocation != null)
{
Distance = myLocation.QuickDistance(myLocation);
Bearing = myLocation.QuickBearing(myLocation);
VerticalAngle = GetVerticalAngle(Distance, alt - myLocation.Altitude);
}
}
