forked from alpinechough/Alpinechough.Srtm
-
Notifications
You must be signed in to change notification settings - Fork 0
/
SrtmDataCell.cs
166 lines (142 loc) · 5.16 KB
/
SrtmDataCell.cs
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
#region MIT License
// MIT License
// Copyright (c) 2012 Alpine Chough Software.
//
// Permission is hereby granted, free of charge, to any person obtaining a copy
// of this software and associated documentation files (the "Software"), to deal
// in the Software without restriction, including without limitation the rights
// to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
// copies of the Software, and to permit persons to whom the Software is
// furnished to do so, subject to the following conditions:
//
// The above copyright notice and this permission notice shall be included in all
// copies or substantial portions of the Software.
//
// THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
// IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
// FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
// AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
// LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
// OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
// SOFTWARE.
#endregion
using System;
using System.IO;
namespace Alpinechough.Srtm
{
/// <summary>
/// SRTM data cell.
/// </summary>
/// <exception cref='FileNotFoundException'>
/// Is thrown when a file path argument specifies a file that does not exist.
/// </exception>
/// <exception cref='ArgumentException'>
/// Is thrown when an argument passed to a method is invalid.
/// </exception>
/// <exception cref='ArgumentOutOfRangeException'>
/// Is thrown when an argument passed to a method is invalid because it is outside the allowable range of values as
/// specified by the method.
/// </exception>
public class SrtmDataCell
{
#region Lifecycle
/// <summary>
/// Initializes a new instance of the <see cref="Alpinechough.Srtm.SrtmDataCell"/> class.
/// </summary>
/// <param name='filepath'>
/// Filepath.
/// </param>
/// <exception cref='FileNotFoundException'>
/// Is thrown when a file path argument specifies a file that does not exist.
/// </exception>
/// <exception cref='ArgumentException'>
/// Is thrown when an argument passed to a method is invalid.
/// </exception>
public SrtmDataCell (string filepath)
{
if (!File.Exists (filepath))
throw new FileNotFoundException ("File not found.", filepath);
if (string.Compare (".hgt", Path.GetExtension (filepath), StringComparison.CurrentCultureIgnoreCase) != 0)
throw new ArgumentException ("Invalid extension.", filepath);
string filename = Path.GetFileNameWithoutExtension (filepath).ToLower ();
string[] fileCoordinate = filename.Split (new [] { 'e', 'w' });
if (fileCoordinate.Length != 2)
throw new ArgumentException ("Invalid filename.", filepath);
fileCoordinate [0] = fileCoordinate [0].TrimStart (new [] { 'n', 's' });
Latitude = int.Parse (fileCoordinate [0]);
if (filename.Contains ("s"))
Latitude *= -1;
Longitude = int.Parse (fileCoordinate [1]);
if (filename.Contains ("w"))
Longitude *= -1;
HgtData = File.ReadAllBytes (filepath);
switch (HgtData.Length) {
case 1201*1201*2: // SRTM-3
PointsPerCell = 1201;
break;
case 3601*3601*2: // SRTM-1
PointsPerCell = 3601;
break;
default:
throw new ArgumentException ("Invalid file size.", filepath);
}
}
#endregion
#region Properties
/// <summary>
/// Gets or sets the hgt data.
/// </summary>
/// <value>
/// The hgt data.
/// </value>
private byte[] HgtData { get; set; }
/// <summary>
/// Gets or sets the points per cell.
/// </summary>
/// <value>
/// The points per cell.
/// </value>
private int PointsPerCell { get; set; }
/// <summary>
/// Gets or sets the latitude of the srtm data file.
/// </summary>
/// <value>
/// The latitude.
/// </value>
public int Latitude { get; private set; }
/// <summary>
/// Gets or sets the longitude of the srtm data file.
/// </summary>
/// <value>
/// The longitude.
/// </value>
public int Longitude { get; private set; }
#endregion
#region Public Methods
/// <summary>
/// Gets the height.
/// </summary>
/// <returns>
/// The height.
/// </returns>
/// <param name='coordinates'>
/// Coordinates.
/// </param>
/// <exception cref='ArgumentOutOfRangeException'>
/// Is thrown when an argument passed to a method is invalid because it is outside the allowable range of values as
/// specified by the method.
/// </exception>
public int GetHeight (IGeographicalCoordinates coordinates)
{
int localLat = (int)((coordinates.Latitude - Latitude) * PointsPerCell);
int localLon = (int)(((coordinates.Longitude - Longitude)) * PointsPerCell);
int bytesPos = ((PointsPerCell - localLat - 1) * PointsPerCell * 2) + localLon * 2;
Console.WriteLine (bytesPos);
if (bytesPos < 0 || bytesPos > PointsPerCell * PointsPerCell * 2)
throw new ArgumentOutOfRangeException ("Coordinates out of range.", "coordinates");
// Motorola "big-endian" order with the most significant byte first
return (HgtData [bytesPos]) << 8 | HgtData [bytesPos + 1];
}
#endregion
}
}