public static altresponce getAltitude(double lat, double lng, double zoom = 16)
{
short alt = 0;
var answer = new altresponce();
var trytiff = Utilities.GeoTiff.getAltitude(lat, lng);
if (trytiff.currenttype == tiletype.valid)
return trytiff;
//lat += 1 / 1199.0;
//lng -= 1 / 1201f;
// lat -35.115676879882812 double
// lng 117.94178754638671 double
// alt 70 short
int x = (lng < 0) ? (int)(lng - 1) : (int)lng;//(int)Math.Floor(lng);
int y = (lat < 0) ? (int)(lat - 1) : (int)lat; ;//(int)Math.Floor(lat);
string ns;
if (y > 0)
ns = "N";
else
ns = "S";
string ew;
if (x > 0)
ew = "E";
else
ew = "W";
// running tostring at a high rate was costing cpu
if (fnamecache[y] == null)
fnamecache[y] = Math.Abs(y).ToString("00");
if (fnamecache[1000 + x] == null)
fnamecache[1000 + x] = Math.Abs(x).ToString("000");
string filename = ns + fnamecache[y] + ew + fnamecache[1000 + x] + ".hgt";
try
{
if (cache.ContainsKey(filename) || File.Exists(datadirectory + Path.DirectorySeparatorChar + filename))
{ // srtm hgt files
int size = -1;
// add to cache
if (!cache.ContainsKey(filename))
{
FileStream fs = new FileStream(datadirectory + Path.DirectorySeparatorChar + filename, FileMode.Open, FileAccess.Read, FileShare.Read);
if (fs.Length == (1201 * 1201 * 2))
{
size = 1201;
}
else if (fs.Length == (3601 * 3601 * 2))
{
size = 3601;
}
else
return answer;
byte[] altbytes = new byte[2];
short[,] altdata = new short[size, size];
int altlat = 0;
int altlng = 0;
while (fs.Read(altbytes, 0, 2) != 0)
{
altdata[altlat, altlng] = (short)((altbytes[0] << 8) + altbytes[1]);
altlat++;
if (altlat >= size)
{
altlng++;
altlat = 0;
}
}
fs.Close();
cache[filename] = altdata;
}
if (cache[filename].Length == (1201 * 1201))
{
size = 1201;
}
else if (cache[filename].Length == (3601 * 3601))
{
size = 3601;
}
else
return answer;
// remove the base lat long
lat -= y;
lng -= x;
// values should be 0-1199, 1200 is for interpolation
double xf = lng * (size - 2);
double yf = lat * (size - 2);
int x_int = (int)xf;
double x_frac = xf - x_int;
int y_int = (int)yf;
double y_frac = yf - y_int;
y_int = (size - 2) - y_int;
double alt00 = GetAlt(filename, x_int, y_int);
double alt10 = GetAlt(filename, x_int + 1, y_int);
double alt01 = GetAlt(filename, x_int, y_int + 1);
double alt11 = GetAlt(filename, x_int + 1, y_int + 1);
double v1 = avg(alt00, alt10, x_frac);
double v2 = avg(alt01, alt11, x_frac);
double v = avg(v1, v2, -y_frac);
answer.currenttype = tiletype.valid;
answer.alt = v;
return answer;
}
string filename2 = "srtm_" + Math.Round((lng + 2.5 + 180) / 5, 0).ToString("00") + "_" + Math.Round((60 - lat + 2.5) / 5, 0).ToString("00") + ".asc";
if (File.Exists(datadirectory + Path.DirectorySeparatorChar + filename2))
{
StreamReader sr = new StreamReader(readFile(datadirectory + Path.DirectorySeparatorChar + filename2));
int nox = 0;
int noy = 0;
float left = 0;
float top = 0;
int nodata = -9999;
float cellsize = 0;
int rowcounter = 0;
float wantrow = 0;
float wantcol = 0;
while (!sr.EndOfStream)
{
string line = sr.ReadLine();
if (line.StartsWith("ncols"))
{
nox = int.Parse(line.Substring(line.IndexOf(' ')));
//hgtdata = new int[nox * noy];
}
else if (line.StartsWith("nrows"))
{
noy = int.Parse(line.Substring(line.IndexOf(' ')));
//hgtdata = new int[nox * noy];
}
else if (line.StartsWith("xllcorner"))
{
left = float.Parse(line.Substring(line.IndexOf(' ')));
}
else if (line.StartsWith("yllcorner"))
{
top = float.Parse(line.Substring(line.IndexOf(' ')));
}
else if (line.StartsWith("cellsize"))
{
cellsize = float.Parse(line.Substring(line.IndexOf(' ')));
}
else if (line.StartsWith("NODATA_value"))
{
nodata = int.Parse(line.Substring(line.IndexOf(' ')));
}
else
{
string[] data = line.Split(new char[] { ' ' });
if (data.Length == (nox + 1))
{
wantcol = (float)((lng - Math.Round(left, 0)));
wantrow = (float)((lat - Math.Round(top, 0)));
wantrow = (int)(wantrow / cellsize);
wantcol = (int)(wantcol / cellsize);
wantrow = noy - wantrow;
if (rowcounter == wantrow)
{
Console.WriteLine("{0} {1} {2} {3} ans {4} x {5}", lng, lat, left, top, data[(int)wantcol], (nox + wantcol * cellsize));
answer.currenttype = tiletype.valid;
answer.alt = int.Parse(data[(int)wantcol]);
return answer;
}
rowcounter++;
}
}
}
//sr.Close();
answer.currenttype = tiletype.valid;
answer.alt = alt;
return answer;
}
else // get something
{
if (filename.Contains("S00W000") || filename.Contains("S00W001") ||
filename.Contains("S01W000") || filename.Contains("S01W001"))
{
answer.currenttype = tiletype.ocean;
return answer;
}
if (oceantile.Contains(filename))
answer.currenttype = tiletype.ocean;
if (zoom >= 12)
{
if (!Directory.Exists(datadirectory))
Directory.CreateDirectory(datadirectory);
if (requestThread == null)
{
log.Info("Getting " + filename);
lock (objlock)
{
queue.Add(filename);
}
requestThread = new Thread(requestRunner);
requestThread.IsBackground = true;
requestThread.Name = "SRTM request runner";
requestThread.Start();
}
else
{
lock (objlock)
{
if (!queue.Contains(filename))
{
log.Info("Getting " + filename);
queue.Add(filename);
}
}
}
}
}
}
catch { answer.alt = 0; answer.currenttype = tiletype.invalid; }
return answer;
}
public static altresponce getAltitude(double lat, double lng, double zoom = 16)
{
short alt = 0;
var answer = new altresponce();
//lat += 1 / 1199.0;
//lng -= 1 / 1201f;
// lat -35.115676879882812 double
// lng 117.94178754638671 double
// alt 70 short
int x = (lng < 0) ? (int)(lng - 1) : (int)lng; //(int)Math.Floor(lng);
int y = (lat < 0) ? (int)(lat - 1) : (int)lat;; //(int)Math.Floor(lat);
string ns;
if (y > 0)
{
ns = "N";
}
else
{
ns = "S";
}
string ew;
if (x > 0)
{
ew = "E";
}
else
{
ew = "W";
}
// running tostring at a high rate was costing cpu
if (fnamecache[y] == null)
{
fnamecache[y] = Math.Abs(y).ToString("00");
}
if (fnamecache[1000 + x] == null)
{
fnamecache[1000 + x] = Math.Abs(x).ToString("000");
}
string filename = ns + fnamecache[y] + ew + fnamecache[1000 + x] + ".hgt";
try
{
if (cache.ContainsKey(filename) || File.Exists(datadirectory + Path.DirectorySeparatorChar + filename))
{ // srtm hgt files
int size = -1;
// add to cache
if (!cache.ContainsKey(filename))
{
FileStream fs = new FileStream(datadirectory + Path.DirectorySeparatorChar + filename, FileMode.Open, FileAccess.Read, FileShare.Read);
if (fs.Length == (1201 * 1201 * 2))
{
size = 1201;
}
else if (fs.Length == (3601 * 3601 * 2))
{
size = 3601;
}
else
{
return(answer);
}
byte[] altbytes = new byte[2];
short[,] altdata = new short[size, size];
int altlat = 0;
int altlng = 0;
while (fs.Read(altbytes, 0, 2) != 0)
{
altdata[altlat, altlng] = (short)((altbytes[0] << 8) + altbytes[1]);
altlat++;
if (altlat >= size)
{
altlng++;
altlat = 0;
}
}
fs.Close();
cache[filename] = altdata;
}
if (cache[filename].Length == (1201 * 1201))
{
size = 1201;
}
else if (cache[filename].Length == (3601 * 3601))
{
size = 3601;
}
else
{
return(answer);
}
// remove the base lat long
lat -= y;
lng -= x;
// values should be 0-1199, 1200 is for interpolation
double xf = lng * (size - 2);
double yf = lat * (size - 2);
int x_int = (int)xf;
double x_frac = xf - x_int;
int y_int = (int)yf;
double y_frac = yf - y_int;
y_int = (size - 2) - y_int;
double alt00 = GetAlt(filename, x_int, y_int);
double alt10 = GetAlt(filename, x_int + 1, y_int);
double alt01 = GetAlt(filename, x_int, y_int + 1);
double alt11 = GetAlt(filename, x_int + 1, y_int + 1);
double v1 = avg(alt00, alt10, x_frac);
double v2 = avg(alt01, alt11, x_frac);
double v = avg(v1, v2, -y_frac);
answer.currenttype = tiletype.valid;
answer.alt = v;
return(answer);
}
string filename2 = "srtm_" + Math.Round((lng + 2.5 + 180) / 5, 0).ToString("00") + "_" + Math.Round((60 - lat + 2.5) / 5, 0).ToString("00") + ".asc";
if (File.Exists(datadirectory + Path.DirectorySeparatorChar + filename2))
{
StreamReader sr = new StreamReader(readFile(datadirectory + Path.DirectorySeparatorChar + filename2));
int nox = 0;
int noy = 0;
float left = 0;
float top = 0;
int nodata = -9999;
float cellsize = 0;
int rowcounter = 0;
float wantrow = 0;
float wantcol = 0;
while (!sr.EndOfStream)
{
string line = sr.ReadLine();
if (line.StartsWith("ncols"))
{
nox = int.Parse(line.Substring(line.IndexOf(' ')));
//hgtdata = new int[nox * noy];
}
else if (line.StartsWith("nrows"))
{
noy = int.Parse(line.Substring(line.IndexOf(' ')));
//hgtdata = new int[nox * noy];
}
else if (line.StartsWith("xllcorner"))
{
left = float.Parse(line.Substring(line.IndexOf(' ')));
}
else if (line.StartsWith("yllcorner"))
{
top = float.Parse(line.Substring(line.IndexOf(' ')));
}
else if (line.StartsWith("cellsize"))
{
cellsize = float.Parse(line.Substring(line.IndexOf(' ')));
}
else if (line.StartsWith("NODATA_value"))
{
nodata = int.Parse(line.Substring(line.IndexOf(' ')));
}
else
{
string[] data = line.Split(new char[] { ' ' });
if (data.Length == (nox + 1))
{
wantcol = (float)((lng - Math.Round(left, 0)));
wantrow = (float)((lat - Math.Round(top, 0)));
wantrow = (int)(wantrow / cellsize);
wantcol = (int)(wantcol / cellsize);
wantrow = noy - wantrow;
if (rowcounter == wantrow)
{
Console.WriteLine("{0} {1} {2} {3} ans {4} x {5}", lng, lat, left, top, data[(int)wantcol], (nox + wantcol * cellsize));
answer.currenttype = tiletype.valid;
answer.alt = int.Parse(data[(int)wantcol]);
return(answer);
}
rowcounter++;
}
}
}
//sr.Close();
answer.currenttype = tiletype.valid;
answer.alt = alt;
return(answer);
}
else // get something
{
if (oceantile.Contains(filename))
{
answer.currenttype = tiletype.ocean;
}
if (zoom >= 12)
{
if (!Directory.Exists(datadirectory))
{
Directory.CreateDirectory(datadirectory);
}
if (requestThread == null)
{
Console.WriteLine("Getting " + filename);
lock (objlock)
{
queue.Add(filename);
}
requestThread = new Thread(requestRunner);
requestThread.IsBackground = true;
requestThread.Name = "SRTM request runner";
requestThread.Start();
}
else
{
lock (objlock)
{
if (!queue.Contains(filename))
{
Console.WriteLine("Getting " + filename);
queue.Add(filename);
}
}
}
}
}
}
catch { answer.alt = 0; answer.currenttype = tiletype.invalid; }
return(answer);
}
