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 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 217 218 219 220 221 222 223 224 225 226 227 228 229 230 231 232 233 234 235 236 237 238 239 240 241 242 243 244 245 246 247 248 249 250 251 252 253 254 255 256 257 258 259 260 261 262 263 264 265 266 267 268 269 270 271 272 273 274 275 276 277 278 279 280 281 282 283 284 285 286 287 288 289 290 291 292 293 294 295 296 297 298 299 300 301 302 303 304 305 306 307 308 309 310 311 312 313 314 315 316 317 318 319 320 321 322 323 324 325 326 327 328 329 330 331 332 333 334 335 336 337 338 339 340 341 342 343 344 345 346 347 348 349 350 351 352 353 354 355 356 357 358 359 360 361 362 363 364 365 366 367 368 369 370 371 372 373 374 375 376 377 378 379 380 381 382 383 384 385 386 387 388 389 390 391 392 393 394 395 396 397 398 399 400 401 402 403 404 405 406 407 408 409 410 411 412 413 414 415 416 417 418 419 420 421 422 423 424 425 426 427 428 429 430 431 432 433 434 435 436 437 438 439 440 441 442 443 444 445 446 447
|
<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "http://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd">
<html xmlns="http://www.w3.org/1999/xhtml">
<head>
<meta http-equiv="Content-Type" content="text/xhtml;charset=UTF-8"/>
<meta http-equiv="X-UA-Compatible" content="IE=9"/>
<title>GeographicTranslator: C:/master/geotrans3.7/CCS/src/dtcc/CoordinateSystems/neys/Neys.cpp Source File</title>
<b><p style="text-align:center;color:green;font-size:20px">UNCLASSIFIED</p></b>
<link href="tabs.css" rel="stylesheet" type="text/css"/>
<script type="text/javascript" src="jquery.js"></script>
<script type="text/javascript" src="dynsections.js"></script>
<link href="search/search.css" rel="stylesheet" type="text/css"/>
<script type="text/javascript" src="search/search.js"></script>
<script type="text/javascript">
$(document).ready(function() { searchBox.OnSelectItem(0); });
</script>
<link href="doxygen.css" rel="stylesheet" type="text/css" />
</head>
<body>
<div id="top"><!-- do not remove this div, it is closed by doxygen! -->
<div id="titlearea">
<table cellspacing="0" cellpadding="0">
<tbody>
<tr style="height: 56px;">
<td style="padding-left: 0.5em;">
<div id="projectname">GeographicTranslator
</div>
</td>
</tr>
</tbody>
</table>
</div>
<!-- end header part -->
<!-- Generated by Doxygen 1.8.2 -->
<script type="text/javascript">
var searchBox = new SearchBox("searchBox", "search",false,'Search');
</script>
<div id="navrow1" class="tabs">
<ul class="tablist">
<li><a href="index.html"><span>Main Page</span></a></li>
<li><a href="namespaces.html"><span>Namespaces</span></a></li>
<li><a href="annotated.html"><span>Classes</span></a></li>
<li class="current"><a href="files.html"><span>Files</span></a></li>
<li>
<div id="MSearchBox" class="MSearchBoxInactive">
<span class="left">
<img id="MSearchSelect" src="search/mag_sel.png"
onmouseover="return searchBox.OnSearchSelectShow()"
onmouseout="return searchBox.OnSearchSelectHide()"
alt=""/>
<input type="text" id="MSearchField" value="Search" accesskey="S"
onfocus="searchBox.OnSearchFieldFocus(true)"
onblur="searchBox.OnSearchFieldFocus(false)"
onkeyup="searchBox.OnSearchFieldChange(event)"/>
</span><span class="right">
<a id="MSearchClose" href="javascript:searchBox.CloseResultsWindow()"><img id="MSearchCloseImg" border="0" src="search/close.png" alt=""/></a>
</span>
</div>
</li>
</ul>
</div>
<div id="navrow2" class="tabs2">
<ul class="tablist">
<li><a href="files.html"><span>File List</span></a></li>
<li><a href="globals.html"><span>File Members</span></a></li>
</ul>
</div>
<!-- window showing the filter options -->
<div id="MSearchSelectWindow"
onmouseover="return searchBox.OnSearchSelectShow()"
onmouseout="return searchBox.OnSearchSelectHide()"
onkeydown="return searchBox.OnSearchSelectKey(event)">
<a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(0)"><span class="SelectionMark"> </span>All</a><a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(1)"><span class="SelectionMark"> </span>Classes</a><a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(2)"><span class="SelectionMark"> </span>Namespaces</a><a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(3)"><span class="SelectionMark"> </span>Files</a><a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(4)"><span class="SelectionMark"> </span>Functions</a><a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(5)"><span class="SelectionMark"> </span>Variables</a><a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(6)"><span class="SelectionMark"> </span>Enumerations</a><a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(7)"><span class="SelectionMark"> </span>Enumerator</a><a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(8)"><span class="SelectionMark"> </span>Friends</a><a class="SelectItem" href="javascript:void(0)" onclick="searchBox.OnSelectItem(9)"><span class="SelectionMark"> </span>Macros</a></div>
<!-- iframe showing the search results (closed by default) -->
<div id="MSearchResultsWindow">
<iframe src="javascript:void(0)" frameborder="0"
name="MSearchResults" id="MSearchResults">
</iframe>
</div>
<div id="nav-path" class="navpath">
<ul>
<li class="navelem"><a class="el" href="dir_edf3df7db5b0053f7cba67b155deccfc.html">master</a></li><li class="navelem"><a class="el" href="dir_8b3490cac915c107683d68df3d9c98a6.html">geotrans3.7</a></li><li class="navelem"><a class="el" href="dir_d0951f0e396e7f3ab5db438f6b9f5363.html">CCS</a></li><li class="navelem"><a class="el" href="dir_3b11bca187df89c3daa83d17707827d4.html">src</a></li><li class="navelem"><a class="el" href="dir_ce1728a50ed0f8e5fbefa80fc2874faa.html">dtcc</a></li><li class="navelem"><a class="el" href="dir_5d8cc6cdb4437c7c87ccad9605085cdb.html">CoordinateSystems</a></li><li class="navelem"><a class="el" href="dir_96a567b782bc78a29560926c8b1cd3f9.html">neys</a></li> </ul>
</div>
</div><!-- top -->
<div class="header">
<div class="headertitle">
<div class="title">Neys.cpp</div> </div>
</div><!--header-->
<div class="contents">
<a href="_neys_8cpp.html">Go to the documentation of this file.</a><div class="fragment"><div class="line"><a name="l00001"></a><span class="lineno"> 1</span> <span class="comment">// CLASSIFICATION: UNCLASSIFIED</span></div>
<div class="line"><a name="l00002"></a><span class="lineno"> 2</span> </div>
<div class="line"><a name="l00003"></a><span class="lineno"> 3</span> <span class="comment">/***************************************************************************/</span></div>
<div class="line"><a name="l00004"></a><span class="lineno"> 4</span> <span class="comment">/* RSC IDENTIFIER: NEYS</span></div>
<div class="line"><a name="l00005"></a><span class="lineno"> 5</span> <span class="comment"> *</span></div>
<div class="line"><a name="l00006"></a><span class="lineno"> 6</span> <span class="comment"> * ABSTRACT</span></div>
<div class="line"><a name="l00007"></a><span class="lineno"> 7</span> <span class="comment"> *</span></div>
<div class="line"><a name="l00008"></a><span class="lineno"> 8</span> <span class="comment"> * This component provides conversions between Geodetic coordinates</span></div>
<div class="line"><a name="l00009"></a><span class="lineno"> 9</span> <span class="comment"> * (latitude and longitude in radians) and Ney's (Modified Lambert</span></div>
<div class="line"><a name="l00010"></a><span class="lineno"> 10</span> <span class="comment"> * Conformal Conic) projection coordinates (easting and northing in meters).</span></div>
<div class="line"><a name="l00011"></a><span class="lineno"> 11</span> <span class="comment"> *</span></div>
<div class="line"><a name="l00012"></a><span class="lineno"> 12</span> <span class="comment"> * ERROR HANDLING</span></div>
<div class="line"><a name="l00013"></a><span class="lineno"> 13</span> <span class="comment"> *</span></div>
<div class="line"><a name="l00014"></a><span class="lineno"> 14</span> <span class="comment"> * This component checks parameters for valid values. If an invalid value</span></div>
<div class="line"><a name="l00015"></a><span class="lineno"> 15</span> <span class="comment"> * is found the error code is combined with the current error code using</span></div>
<div class="line"><a name="l00016"></a><span class="lineno"> 16</span> <span class="comment"> * the bitwise or. This combining allows multiple error codes to be</span></div>
<div class="line"><a name="l00017"></a><span class="lineno"> 17</span> <span class="comment"> * returned. The possible error codes are:</span></div>
<div class="line"><a name="l00018"></a><span class="lineno"> 18</span> <span class="comment"> *</span></div>
<div class="line"><a name="l00019"></a><span class="lineno"> 19</span> <span class="comment"> * NEYS_NO_ERROR : No errors occurred in function</span></div>
<div class="line"><a name="l00020"></a><span class="lineno"> 20</span> <span class="comment"> * NEYS_LAT_ERROR : Latitude outside of valid range</span></div>
<div class="line"><a name="l00021"></a><span class="lineno"> 21</span> <span class="comment"> * (-90 to 90 degrees)</span></div>
<div class="line"><a name="l00022"></a><span class="lineno"> 22</span> <span class="comment"> * NEYS_LON_ERROR : Longitude outside of valid range</span></div>
<div class="line"><a name="l00023"></a><span class="lineno"> 23</span> <span class="comment"> * (-180 to 360 degrees)</span></div>
<div class="line"><a name="l00024"></a><span class="lineno"> 24</span> <span class="comment"> * NEYS_EASTING_ERROR : Easting outside of valid range</span></div>
<div class="line"><a name="l00025"></a><span class="lineno"> 25</span> <span class="comment"> * (depends on ellipsoid and projection</span></div>
<div class="line"><a name="l00026"></a><span class="lineno"> 26</span> <span class="comment"> * parameters)</span></div>
<div class="line"><a name="l00027"></a><span class="lineno"> 27</span> <span class="comment"> * NEYS_NORTHING_ERROR : Northing outside of valid range</span></div>
<div class="line"><a name="l00028"></a><span class="lineno"> 28</span> <span class="comment"> * (depends on ellipsoid and projection</span></div>
<div class="line"><a name="l00029"></a><span class="lineno"> 29</span> <span class="comment"> * parameters)</span></div>
<div class="line"><a name="l00030"></a><span class="lineno"> 30</span> <span class="comment"> * NEYS_FIRST_STDP_ERROR : First standard parallel outside of valid</span></div>
<div class="line"><a name="l00031"></a><span class="lineno"> 31</span> <span class="comment"> * range (71 or 74 degrees)</span></div>
<div class="line"><a name="l00032"></a><span class="lineno"> 32</span> <span class="comment"> * NEYS_ORIGIN_LAT_ERROR : Origin latitude outside of valid range</span></div>
<div class="line"><a name="l00033"></a><span class="lineno"> 33</span> <span class="comment"> * (-89 59 59.0 to 89 59 59.0 degrees)</span></div>
<div class="line"><a name="l00034"></a><span class="lineno"> 34</span> <span class="comment"> * NEYS_CENT_MER_ERROR : Central meridian outside of valid range</span></div>
<div class="line"><a name="l00035"></a><span class="lineno"> 35</span> <span class="comment"> * (-180 to 360 degrees)</span></div>
<div class="line"><a name="l00036"></a><span class="lineno"> 36</span> <span class="comment"> * NEYS_A_ERROR : Semi-major axis less than or equal to zero</span></div>
<div class="line"><a name="l00037"></a><span class="lineno"> 37</span> <span class="comment"> * NEYS_INV_F_ERROR : Inverse flattening outside of valid range</span></div>
<div class="line"><a name="l00038"></a><span class="lineno"> 38</span> <span class="comment"> * (250 to 350)</span></div>
<div class="line"><a name="l00039"></a><span class="lineno"> 39</span> <span class="comment"> *</span></div>
<div class="line"><a name="l00040"></a><span class="lineno"> 40</span> <span class="comment"> *</span></div>
<div class="line"><a name="l00041"></a><span class="lineno"> 41</span> <span class="comment"> * REUSE NOTES</span></div>
<div class="line"><a name="l00042"></a><span class="lineno"> 42</span> <span class="comment"> *</span></div>
<div class="line"><a name="l00043"></a><span class="lineno"> 43</span> <span class="comment"> * NEYS is intended for reuse by any application that performs a Ney's (Modified</span></div>
<div class="line"><a name="l00044"></a><span class="lineno"> 44</span> <span class="comment"> * Lambert Conformal Conic) projection or its inverse.</span></div>
<div class="line"><a name="l00045"></a><span class="lineno"> 45</span> <span class="comment"> *</span></div>
<div class="line"><a name="l00046"></a><span class="lineno"> 46</span> <span class="comment"> * REFERENCES</span></div>
<div class="line"><a name="l00047"></a><span class="lineno"> 47</span> <span class="comment"> *</span></div>
<div class="line"><a name="l00048"></a><span class="lineno"> 48</span> <span class="comment"> * Further information on NEYS can be found in the Reuse Manual.</span></div>
<div class="line"><a name="l00049"></a><span class="lineno"> 49</span> <span class="comment"> *</span></div>
<div class="line"><a name="l00050"></a><span class="lineno"> 50</span> <span class="comment"> * NEYS originated from:</span></div>
<div class="line"><a name="l00051"></a><span class="lineno"> 51</span> <span class="comment"> * U.S. Army Topographic Engineering Center</span></div>
<div class="line"><a name="l00052"></a><span class="lineno"> 52</span> <span class="comment"> * Geospatial Information Division</span></div>
<div class="line"><a name="l00053"></a><span class="lineno"> 53</span> <span class="comment"> * 7701 Telegraph Road</span></div>
<div class="line"><a name="l00054"></a><span class="lineno"> 54</span> <span class="comment"> * Alexandria, VA 22310-3864</span></div>
<div class="line"><a name="l00055"></a><span class="lineno"> 55</span> <span class="comment"> *</span></div>
<div class="line"><a name="l00056"></a><span class="lineno"> 56</span> <span class="comment"> * LICENSES</span></div>
<div class="line"><a name="l00057"></a><span class="lineno"> 57</span> <span class="comment"> *</span></div>
<div class="line"><a name="l00058"></a><span class="lineno"> 58</span> <span class="comment"> * None apply to this component.</span></div>
<div class="line"><a name="l00059"></a><span class="lineno"> 59</span> <span class="comment"> *</span></div>
<div class="line"><a name="l00060"></a><span class="lineno"> 60</span> <span class="comment"> * RESTRICTIONS</span></div>
<div class="line"><a name="l00061"></a><span class="lineno"> 61</span> <span class="comment"> *</span></div>
<div class="line"><a name="l00062"></a><span class="lineno"> 62</span> <span class="comment"> * NEYS has no restrictions.</span></div>
<div class="line"><a name="l00063"></a><span class="lineno"> 63</span> <span class="comment"> *</span></div>
<div class="line"><a name="l00064"></a><span class="lineno"> 64</span> <span class="comment"> * ENVIRONMENT</span></div>
<div class="line"><a name="l00065"></a><span class="lineno"> 65</span> <span class="comment"> *</span></div>
<div class="line"><a name="l00066"></a><span class="lineno"> 66</span> <span class="comment"> * NEYS was tested and certified in the following environments:</span></div>
<div class="line"><a name="l00067"></a><span class="lineno"> 67</span> <span class="comment"> *</span></div>
<div class="line"><a name="l00068"></a><span class="lineno"> 68</span> <span class="comment"> * 1. Solaris 2.5 with GCC, version 2.8.1</span></div>
<div class="line"><a name="l00069"></a><span class="lineno"> 69</span> <span class="comment"> * 2. Windows 95 with MS Visual C++, version 6</span></div>
<div class="line"><a name="l00070"></a><span class="lineno"> 70</span> <span class="comment"> *</span></div>
<div class="line"><a name="l00071"></a><span class="lineno"> 71</span> <span class="comment"> * MODIFICATIONS</span></div>
<div class="line"><a name="l00072"></a><span class="lineno"> 72</span> <span class="comment"> *</span></div>
<div class="line"><a name="l00073"></a><span class="lineno"> 73</span> <span class="comment"> * Date Description</span></div>
<div class="line"><a name="l00074"></a><span class="lineno"> 74</span> <span class="comment"> * ---- -----------</span></div>
<div class="line"><a name="l00075"></a><span class="lineno"> 75</span> <span class="comment"> * 8-4-00 Original Code</span></div>
<div class="line"><a name="l00076"></a><span class="lineno"> 76</span> <span class="comment"> * 3-2-07 Original C++ Code</span></div>
<div class="line"><a name="l00077"></a><span class="lineno"> 77</span> <span class="comment"> *</span></div>
<div class="line"><a name="l00078"></a><span class="lineno"> 78</span> <span class="comment"> *</span></div>
<div class="line"><a name="l00079"></a><span class="lineno"> 79</span> <span class="comment"> *</span></div>
<div class="line"><a name="l00080"></a><span class="lineno"> 80</span> <span class="comment"> */</span></div>
<div class="line"><a name="l00081"></a><span class="lineno"> 81</span> </div>
<div class="line"><a name="l00082"></a><span class="lineno"> 82</span> </div>
<div class="line"><a name="l00083"></a><span class="lineno"> 83</span> <span class="comment">/***************************************************************************/</span></div>
<div class="line"><a name="l00084"></a><span class="lineno"> 84</span> <span class="comment">/*</span></div>
<div class="line"><a name="l00085"></a><span class="lineno"> 85</span> <span class="comment"> * INCLUDES</span></div>
<div class="line"><a name="l00086"></a><span class="lineno"> 86</span> <span class="comment"> */</span></div>
<div class="line"><a name="l00087"></a><span class="lineno"> 87</span> </div>
<div class="line"><a name="l00088"></a><span class="lineno"> 88</span> <span class="preprocessor">#include <math.h></span></div>
<div class="line"><a name="l00089"></a><span class="lineno"> 89</span> <span class="preprocessor">#include "<a class="code" href="_lambert_conformal_conic_8h.html">LambertConformalConic.h</a>"</span></div>
<div class="line"><a name="l00090"></a><span class="lineno"> 90</span> <span class="preprocessor">#include "<a class="code" href="_neys_8h.html">Neys.h</a>"</span></div>
<div class="line"><a name="l00091"></a><span class="lineno"> 91</span> <span class="preprocessor">#include "<a class="code" href="_neys_parameters_8h.html">NeysParameters.h</a>"</span></div>
<div class="line"><a name="l00092"></a><span class="lineno"> 92</span> <span class="preprocessor">#include "<a class="code" href="_map_projection_coordinates_8h.html">MapProjectionCoordinates.h</a>"</span></div>
<div class="line"><a name="l00093"></a><span class="lineno"> 93</span> <span class="preprocessor">#include "<a class="code" href="_geodetic_coordinates_8h.html">GeodeticCoordinates.h</a>"</span></div>
<div class="line"><a name="l00094"></a><span class="lineno"> 94</span> <span class="preprocessor">#include "<a class="code" href="_coordinate_conversion_exception_8h.html">CoordinateConversionException.h</a>"</span></div>
<div class="line"><a name="l00095"></a><span class="lineno"> 95</span> <span class="preprocessor">#include "<a class="code" href="_error_messages_8h.html">ErrorMessages.h</a>"</span></div>
<div class="line"><a name="l00096"></a><span class="lineno"> 96</span> </div>
<div class="line"><a name="l00097"></a><span class="lineno"> 97</span> <span class="comment">/*</span></div>
<div class="line"><a name="l00098"></a><span class="lineno"> 98</span> <span class="comment"> * math.h - Standard C math library</span></div>
<div class="line"><a name="l00099"></a><span class="lineno"> 99</span> <span class="comment"> * LambertConformalConic2.h - Is used to convert lambert conformal conic coordinates</span></div>
<div class="line"><a name="l00100"></a><span class="lineno"> 100</span> <span class="comment"> * Neys.h - Is for prototype error checking</span></div>
<div class="line"><a name="l00101"></a><span class="lineno"> 101</span> <span class="comment"> * MapProjectionCoordinates.h - defines map projection coordinates</span></div>
<div class="line"><a name="l00102"></a><span class="lineno"> 102</span> <span class="comment"> * GeodeticCoordinates.h - defines geodetic coordinates</span></div>
<div class="line"><a name="l00103"></a><span class="lineno"> 103</span> <span class="comment"> * CoordinateConversionException.h - Exception handler</span></div>
<div class="line"><a name="l00104"></a><span class="lineno"> 104</span> <span class="comment"> * ErrorMessages.h - Contains exception messages</span></div>
<div class="line"><a name="l00105"></a><span class="lineno"> 105</span> <span class="comment"> */</span></div>
<div class="line"><a name="l00106"></a><span class="lineno"> 106</span> </div>
<div class="line"><a name="l00107"></a><span class="lineno"> 107</span> </div>
<div class="line"><a name="l00108"></a><span class="lineno"> 108</span> <span class="keyword">using namespace </span>MSP::CCS;</div>
<div class="line"><a name="l00109"></a><span class="lineno"> 109</span> </div>
<div class="line"><a name="l00110"></a><span class="lineno"> 110</span> </div>
<div class="line"><a name="l00111"></a><span class="lineno"> 111</span> <span class="comment">/***************************************************************************/</span></div>
<div class="line"><a name="l00112"></a><span class="lineno"> 112</span> <span class="comment">/* DEFINES </span></div>
<div class="line"><a name="l00113"></a><span class="lineno"> 113</span> <span class="comment"> *</span></div>
<div class="line"><a name="l00114"></a><span class="lineno"> 114</span> <span class="comment"> */</span></div>
<div class="line"><a name="l00115"></a><span class="lineno"> 115</span> </div>
<div class="line"><a name="l00116"></a><span class="lineno"><a class="code" href="_neys_8cpp.html#a952eac791b596a61bba0a133a3bb439f"> 116</a></span> <span class="keyword">const</span> <span class="keywordtype">double</span> <a class="code" href="_m_g_r_s_8cpp.html#a598a3330b3c21701223ee0ca14316eca">PI</a> = 3.14159265358979323e0; <span class="comment">/* PI */</span></div>
<div class="line"><a name="l00117"></a><span class="lineno"><a class="code" href="_neys_8cpp.html#aa3c1d8e4403d757bf419703c040db87d"> 117</a></span> <span class="keyword">const</span> <span class="keywordtype">double</span> <a class="code" href="_m_g_r_s_8cpp.html#a27d2bea15134ae7c7f7a2d11379bde53">PI_OVER_2</a> = (<a class="code" href="_m_g_r_s_8cpp.html#a598a3330b3c21701223ee0ca14316eca">PI</a> / 2.0);</div>
<div class="line"><a name="l00118"></a><span class="lineno"><a class="code" href="_neys_8cpp.html#a99168053ae79aafdba6055129c5bb256"> 118</a></span> <span class="keyword">const</span> <span class="keywordtype">double</span> <a class="code" href="_m_g_r_s_8cpp.html#ab7aee7b3221fb08bac0a8360ad6dfeb4">PI_OVER_180</a> = (<a class="code" href="_m_g_r_s_8cpp.html#a598a3330b3c21701223ee0ca14316eca">PI</a> / 180.0);</div>
<div class="line"><a name="l00119"></a><span class="lineno"><a class="code" href="_neys_8cpp.html#a820597b124334bf7fb85214114f4876d"> 119</a></span> <span class="keyword">const</span> <span class="keywordtype">double</span> <a class="code" href="_albers_equal_area_conic_8cpp.html#a820597b124334bf7fb85214114f4876d">TWO_PI</a> = (2.0 * <a class="code" href="_m_g_r_s_8cpp.html#a598a3330b3c21701223ee0ca14316eca">PI</a>);</div>
<div class="line"><a name="l00120"></a><span class="lineno"><a class="code" href="_neys_8cpp.html#ae3709f496dd9694f02743993e8f1707a"> 120</a></span> <span class="keyword">const</span> <span class="keywordtype">double</span> <a class="code" href="_neys_8cpp.html#ae3709f496dd9694f02743993e8f1707a">SEVENTY_ONE</a> = (71.0 * <a class="code" href="_m_g_r_s_8cpp.html#ab7aee7b3221fb08bac0a8360ad6dfeb4">PI_OVER_180</a>); <span class="comment">/* 71 degrees */</span></div>
<div class="line"><a name="l00121"></a><span class="lineno"><a class="code" href="_neys_8cpp.html#a94b9121fe29e16829fab28211e5c2a8e"> 121</a></span> <span class="keyword">const</span> <span class="keywordtype">double</span> <a class="code" href="_neys_8cpp.html#a94b9121fe29e16829fab28211e5c2a8e">SEVENTY_FOUR</a> = (74.0 * <a class="code" href="_m_g_r_s_8cpp.html#ab7aee7b3221fb08bac0a8360ad6dfeb4">PI_OVER_180</a>); <span class="comment">/* 74 degrees */</span></div>
<div class="line"><a name="l00123"></a><span class="lineno"><a class="code" href="_neys_8cpp.html#aa63b6cd7807c756d972c38009904b2a2"> 123</a></span> <span class="keyword">const</span> <span class="keywordtype">double</span> <a class="code" href="_u_t_m_8cpp.html#a7ce12d8a28541150544bc0dfeb789e13">MAX_LAT</a> = (89.999444444444444 * <a class="code" href="_m_g_r_s_8cpp.html#ab7aee7b3221fb08bac0a8360ad6dfeb4">PI_OVER_180</a>); <span class="comment">/* 89 59 58.0 degrees */</span></div>
<div class="line"><a name="l00124"></a><span class="lineno"> 124</span> </div>
<div class="line"><a name="l00125"></a><span class="lineno"> 125</span> </div>
<div class="line"><a name="l00126"></a><span class="lineno"> 126</span> <span class="comment">/************************************************************************/</span></div>
<div class="line"><a name="l00127"></a><span class="lineno"> 127</span> <span class="comment">/* FUNCTIONS </span></div>
<div class="line"><a name="l00128"></a><span class="lineno"> 128</span> <span class="comment"> *</span></div>
<div class="line"><a name="l00129"></a><span class="lineno"> 129</span> <span class="comment"> */</span></div>
<div class="line"><a name="l00130"></a><span class="lineno"> 130</span> </div>
<div class="line"><a name="l00131"></a><span class="lineno"><a class="code" href="class_m_s_p_1_1_c_c_s_1_1_neys.html#a9baf4193928c4ad3422144717a21e0ef"> 131</a></span> <a class="code" href="class_m_s_p_1_1_c_c_s_1_1_neys.html#a9baf4193928c4ad3422144717a21e0ef">Neys::Neys</a>( <span class="keywordtype">double</span> ellipsoidSemiMajorAxis, <span class="keywordtype">double</span> ellipsoidFlattening, <span class="keywordtype">double</span> centralMeridian, <span class="keywordtype">double</span> originLatitude, <span class="keywordtype">double</span> standardParallel, <span class="keywordtype">double</span> falseEasting, <span class="keywordtype">double</span> falseNorthing ) :</div>
<div class="line"><a name="l00132"></a><span class="lineno"> 132</span>  <a class="code" href="class_m_s_p_1_1_c_c_s_1_1_coordinate_system.html">CoordinateSystem</a>(),</div>
<div class="line"><a name="l00133"></a><span class="lineno"> 133</span>  lambertConformalConic2( 0 ),</div>
<div class="line"><a name="l00134"></a><span class="lineno"> 134</span>  Neys_Std_Parallel_1( <a class="code" href="_neys_8cpp.html#ae3709f496dd9694f02743993e8f1707a">SEVENTY_ONE</a> ),</div>
<div class="line"><a name="l00135"></a><span class="lineno"> 135</span>  Neys_Std_Parallel_2( <a class="code" href="_u_t_m_8cpp.html#a7ce12d8a28541150544bc0dfeb789e13">MAX_LAT</a> ),</div>
<div class="line"><a name="l00136"></a><span class="lineno"> 136</span>  Neys_Origin_Long( 0.0 ),</div>
<div class="line"><a name="l00137"></a><span class="lineno"> 137</span>  Neys_Origin_Lat( 80.0 * <a class="code" href="_m_g_r_s_8cpp.html#ab7aee7b3221fb08bac0a8360ad6dfeb4">PI_OVER_180</a> ),</div>
<div class="line"><a name="l00138"></a><span class="lineno"> 138</span>  Neys_False_Easting( 0.0 ),</div>
<div class="line"><a name="l00139"></a><span class="lineno"> 139</span>  Neys_False_Northing( 0.0 ),</div>
<div class="line"><a name="l00140"></a><span class="lineno"> 140</span>  Neys_Delta_Easting( 400000000.0 ),</div>
<div class="line"><a name="l00141"></a><span class="lineno"> 141</span>  Neys_Delta_Northing( 400000000.0 )</div>
<div class="line"><a name="l00142"></a><span class="lineno"> 142</span> {</div>
<div class="line"><a name="l00143"></a><span class="lineno"> 143</span> <span class="comment">/*</span></div>
<div class="line"><a name="l00144"></a><span class="lineno"> 144</span> <span class="comment"> * The constructor receives the ellipsoid parameters and</span></div>
<div class="line"><a name="l00145"></a><span class="lineno"> 145</span> <span class="comment"> * Ney's (Modified Lambert Conformal Conic) projection parameters as inputs, and sets the</span></div>
<div class="line"><a name="l00146"></a><span class="lineno"> 146</span> <span class="comment"> * corresponding state variables. If any errors occur, an exception is thrown </span></div>
<div class="line"><a name="l00147"></a><span class="lineno"> 147</span> <span class="comment"> * with a description of the error.</span></div>
<div class="line"><a name="l00148"></a><span class="lineno"> 148</span> <span class="comment"> *</span></div>
<div class="line"><a name="l00149"></a><span class="lineno"> 149</span> <span class="comment"> * ellipsoidSemiMajorAxis : Semi-major axis of ellipsoid, in meters (input)</span></div>
<div class="line"><a name="l00150"></a><span class="lineno"> 150</span> <span class="comment"> * ellipsoidFlattening : Flattening of ellipsoid (input)</span></div>
<div class="line"><a name="l00151"></a><span class="lineno"> 151</span> <span class="comment"> * centralMeridian : Longitude of origin, in radians (input)</span></div>
<div class="line"><a name="l00152"></a><span class="lineno"> 152</span> <span class="comment"> * originLatitude : Latitude of origin, in radians (input)</span></div>
<div class="line"><a name="l00153"></a><span class="lineno"> 153</span> <span class="comment"> * standardParallel : First standard parallel, in radians (input)</span></div>
<div class="line"><a name="l00154"></a><span class="lineno"> 154</span> <span class="comment"> * falseEasting : False easting, in meters (input)</span></div>
<div class="line"><a name="l00155"></a><span class="lineno"> 155</span> <span class="comment"> * falseNorthing : False northing, in meters (input)</span></div>
<div class="line"><a name="l00156"></a><span class="lineno"> 156</span> <span class="comment"> */</span></div>
<div class="line"><a name="l00157"></a><span class="lineno"> 157</span> </div>
<div class="line"><a name="l00158"></a><span class="lineno"> 158</span>  <span class="keywordtype">double</span> epsilon = 1.0e-10;</div>
<div class="line"><a name="l00159"></a><span class="lineno"> 159</span>  <span class="keywordtype">double</span> inv_f = 1 / ellipsoidFlattening;</div>
<div class="line"><a name="l00160"></a><span class="lineno"> 160</span> </div>
<div class="line"><a name="l00161"></a><span class="lineno"> 161</span>  <span class="keywordflow">if</span> (ellipsoidSemiMajorAxis <= 0.0)</div>
<div class="line"><a name="l00162"></a><span class="lineno"> 162</span>  { <span class="comment">/* Semi-major axis must be greater than zero */</span></div>
<div class="line"><a name="l00163"></a><span class="lineno"> 163</span>  <span class="keywordflow">throw</span> <a class="code" href="class_m_s_p_1_1_c_c_s_1_1_coordinate_conversion_exception.html">CoordinateConversionException</a>( <a class="code" href="class_m_s_p_1_1_c_c_s_1_1_error_messages.html#a75f494a85ff1b96b4ada5511f3bcc914">ErrorMessages::semiMajorAxis</a> );</div>
<div class="line"><a name="l00164"></a><span class="lineno"> 164</span>  }</div>
<div class="line"><a name="l00165"></a><span class="lineno"> 165</span>  <span class="keywordflow">if</span> ((inv_f < 250) || (inv_f > 350))</div>
<div class="line"><a name="l00166"></a><span class="lineno"> 166</span>  { <span class="comment">/* Inverse flattening must be between 250 and 350 */</span></div>
<div class="line"><a name="l00167"></a><span class="lineno"> 167</span>  <span class="keywordflow">throw</span> <a class="code" href="class_m_s_p_1_1_c_c_s_1_1_coordinate_conversion_exception.html">CoordinateConversionException</a>( <a class="code" href="class_m_s_p_1_1_c_c_s_1_1_error_messages.html#ac376fe1588b8d462711455ba1ff73fb0">ErrorMessages::ellipsoidFlattening</a> );</div>
<div class="line"><a name="l00168"></a><span class="lineno"> 168</span> </div>
<div class="line"><a name="l00169"></a><span class="lineno"> 169</span>  }</div>
<div class="line"><a name="l00170"></a><span class="lineno"> 170</span>  <span class="keywordflow">if</span> ((originLatitude < -<a class="code" href="_u_t_m_8cpp.html#a7ce12d8a28541150544bc0dfeb789e13">MAX_LAT</a>) || (originLatitude > <a class="code" href="_u_t_m_8cpp.html#a7ce12d8a28541150544bc0dfeb789e13">MAX_LAT</a>))</div>
<div class="line"><a name="l00171"></a><span class="lineno"> 171</span>  { <span class="comment">/* Origin Latitude out of range */</span></div>
<div class="line"><a name="l00172"></a><span class="lineno"> 172</span>  <span class="keywordflow">throw</span> <a class="code" href="class_m_s_p_1_1_c_c_s_1_1_coordinate_conversion_exception.html">CoordinateConversionException</a>( <a class="code" href="class_m_s_p_1_1_c_c_s_1_1_error_messages.html#a801c71c17f54875371096ac56784f924">ErrorMessages::originLatitude</a> );</div>
<div class="line"><a name="l00173"></a><span class="lineno"> 173</span>  }</div>
<div class="line"><a name="l00174"></a><span class="lineno"> 174</span>  <span class="keywordflow">if</span> ((fabs(standardParallel - <a class="code" href="_neys_8cpp.html#ae3709f496dd9694f02743993e8f1707a">SEVENTY_ONE</a>) > epsilon) && (fabs(standardParallel - <a class="code" href="_neys_8cpp.html#a94b9121fe29e16829fab28211e5c2a8e">SEVENTY_FOUR</a>) > epsilon))</div>
<div class="line"><a name="l00175"></a><span class="lineno"> 175</span>  { <span class="comment">/* First Standard Parallel invalid */</span></div>
<div class="line"><a name="l00176"></a><span class="lineno"> 176</span>  <span class="keywordflow">throw</span> <a class="code" href="class_m_s_p_1_1_c_c_s_1_1_coordinate_conversion_exception.html">CoordinateConversionException</a>( <a class="code" href="class_m_s_p_1_1_c_c_s_1_1_error_messages.html#a42c2dfceebd518fdbc9bd67fdf2dd964">ErrorMessages::standardParallel1</a> );</div>
<div class="line"><a name="l00177"></a><span class="lineno"> 177</span>  }</div>
<div class="line"><a name="l00178"></a><span class="lineno"> 178</span>  <span class="keywordflow">if</span> ((centralMeridian < -<a class="code" href="_m_g_r_s_8cpp.html#a598a3330b3c21701223ee0ca14316eca">PI</a>) || (centralMeridian > <a class="code" href="_albers_equal_area_conic_8cpp.html#a820597b124334bf7fb85214114f4876d">TWO_PI</a>))</div>
<div class="line"><a name="l00179"></a><span class="lineno"> 179</span>  { <span class="comment">/* Origin Longitude out of range */</span></div>
<div class="line"><a name="l00180"></a><span class="lineno"> 180</span>  <span class="keywordflow">throw</span> <a class="code" href="class_m_s_p_1_1_c_c_s_1_1_coordinate_conversion_exception.html">CoordinateConversionException</a>( <a class="code" href="class_m_s_p_1_1_c_c_s_1_1_error_messages.html#a101516ca8b915b90297b9e2abf28bd6f">ErrorMessages:: centralMeridian</a> );</div>
<div class="line"><a name="l00181"></a><span class="lineno"> 181</span>  }</div>
<div class="line"><a name="l00182"></a><span class="lineno"> 182</span> </div>
<div class="line"><a name="l00183"></a><span class="lineno"> 183</span>  <a class="code" href="class_m_s_p_1_1_c_c_s_1_1_coordinate_system.html#ab4b25504df8abc3fb23b5f03eaf6c14e">semiMajorAxis</a> = ellipsoidSemiMajorAxis;</div>
<div class="line"><a name="l00184"></a><span class="lineno"> 184</span>  <a class="code" href="class_m_s_p_1_1_c_c_s_1_1_coordinate_system.html#a103abc45b36f0bc52bd992f2d19bbfc6">flattening</a> = ellipsoidFlattening;</div>
<div class="line"><a name="l00185"></a><span class="lineno"> 185</span> </div>
<div class="line"><a name="l00186"></a><span class="lineno"> 186</span>  Neys_Origin_Lat = originLatitude;</div>
<div class="line"><a name="l00187"></a><span class="lineno"> 187</span>  <span class="keywordflow">if</span> (Neys_Origin_Lat >= 0)</div>
<div class="line"><a name="l00188"></a><span class="lineno"> 188</span>  {</div>
<div class="line"><a name="l00189"></a><span class="lineno"> 189</span>  Neys_Std_Parallel_1 = standardParallel;</div>
<div class="line"><a name="l00190"></a><span class="lineno"> 190</span>  Neys_Std_Parallel_2 = <a class="code" href="_u_t_m_8cpp.html#a7ce12d8a28541150544bc0dfeb789e13">MAX_LAT</a>;</div>
<div class="line"><a name="l00191"></a><span class="lineno"> 191</span>  }</div>
<div class="line"><a name="l00192"></a><span class="lineno"> 192</span>  <span class="keywordflow">else</span></div>
<div class="line"><a name="l00193"></a><span class="lineno"> 193</span>  {</div>
<div class="line"><a name="l00194"></a><span class="lineno"> 194</span>  Neys_Std_Parallel_1 = -standardParallel;</div>
<div class="line"><a name="l00195"></a><span class="lineno"> 195</span>  Neys_Std_Parallel_2 = -<a class="code" href="_u_t_m_8cpp.html#a7ce12d8a28541150544bc0dfeb789e13">MAX_LAT</a>;</div>
<div class="line"><a name="l00196"></a><span class="lineno"> 196</span>  }</div>
<div class="line"><a name="l00197"></a><span class="lineno"> 197</span> </div>
<div class="line"><a name="l00198"></a><span class="lineno"> 198</span>  <span class="keywordflow">if</span> (centralMeridian > <a class="code" href="_m_g_r_s_8cpp.html#a598a3330b3c21701223ee0ca14316eca">PI</a>)</div>
<div class="line"><a name="l00199"></a><span class="lineno"> 199</span>  centralMeridian -= <a class="code" href="_albers_equal_area_conic_8cpp.html#a820597b124334bf7fb85214114f4876d">TWO_PI</a>;</div>
<div class="line"><a name="l00200"></a><span class="lineno"> 200</span>  Neys_Origin_Long = centralMeridian;</div>
<div class="line"><a name="l00201"></a><span class="lineno"> 201</span>  Neys_False_Easting = falseEasting;</div>
<div class="line"><a name="l00202"></a><span class="lineno"> 202</span>  Neys_False_Northing = falseNorthing;</div>
<div class="line"><a name="l00203"></a><span class="lineno"> 203</span> </div>
<div class="line"><a name="l00204"></a><span class="lineno"> 204</span>  lambertConformalConic2 = <span class="keyword">new</span> <a class="code" href="class_m_s_p_1_1_c_c_s_1_1_lambert_conformal_conic.html">LambertConformalConic</a>( <a class="code" href="class_m_s_p_1_1_c_c_s_1_1_coordinate_system.html#ab4b25504df8abc3fb23b5f03eaf6c14e">semiMajorAxis</a>, <a class="code" href="class_m_s_p_1_1_c_c_s_1_1_coordinate_system.html#a103abc45b36f0bc52bd992f2d19bbfc6">flattening</a>, Neys_Origin_Long, Neys_Origin_Lat,</div>
<div class="line"><a name="l00205"></a><span class="lineno"> 205</span>  Neys_Std_Parallel_1, Neys_Std_Parallel_2,</div>
<div class="line"><a name="l00206"></a><span class="lineno"> 206</span>  Neys_False_Easting, Neys_False_Northing );</div>
<div class="line"><a name="l00207"></a><span class="lineno"> 207</span> }</div>
<div class="line"><a name="l00208"></a><span class="lineno"> 208</span> </div>
<div class="line"><a name="l00209"></a><span class="lineno"> 209</span> </div>
<div class="line"><a name="l00210"></a><span class="lineno"><a class="code" href="class_m_s_p_1_1_c_c_s_1_1_neys.html#ac6120b63214e3f4510b86f0537943ca7"> 210</a></span> <a class="code" href="class_m_s_p_1_1_c_c_s_1_1_neys.html#a9baf4193928c4ad3422144717a21e0ef">Neys::Neys</a>( <span class="keyword">const</span> <a class="code" href="class_m_s_p_1_1_c_c_s_1_1_neys.html">Neys</a> &n )</div>
<div class="line"><a name="l00211"></a><span class="lineno"> 211</span> {</div>
<div class="line"><a name="l00212"></a><span class="lineno"> 212</span>  lambertConformalConic2 = <span class="keyword">new</span> <a class="code" href="class_m_s_p_1_1_c_c_s_1_1_lambert_conformal_conic.html">LambertConformalConic</a>( *( n.lambertConformalConic2 ) );</div>
<div class="line"><a name="l00213"></a><span class="lineno"> 213</span>  <a class="code" href="class_m_s_p_1_1_c_c_s_1_1_coordinate_system.html#ab4b25504df8abc3fb23b5f03eaf6c14e">semiMajorAxis</a> = n.<a class="code" href="class_m_s_p_1_1_c_c_s_1_1_coordinate_system.html#ab4b25504df8abc3fb23b5f03eaf6c14e">semiMajorAxis</a>;</div>
<div class="line"><a name="l00214"></a><span class="lineno"> 214</span>  <a class="code" href="class_m_s_p_1_1_c_c_s_1_1_coordinate_system.html#a103abc45b36f0bc52bd992f2d19bbfc6">flattening</a> = n.<a class="code" href="class_m_s_p_1_1_c_c_s_1_1_coordinate_system.html#a103abc45b36f0bc52bd992f2d19bbfc6">flattening</a>;</div>
<div class="line"><a name="l00215"></a><span class="lineno"> 215</span>  Neys_Std_Parallel_1 = n.Neys_Std_Parallel_1; </div>
<div class="line"><a name="l00216"></a><span class="lineno"> 216</span>  Neys_Std_Parallel_2 = n.Neys_Std_Parallel_2; </div>
<div class="line"><a name="l00217"></a><span class="lineno"> 217</span>  Neys_Origin_Long = n.Neys_Origin_Long; </div>
<div class="line"><a name="l00218"></a><span class="lineno"> 218</span>  Neys_Origin_Lat = n.Neys_Origin_Lat; </div>
<div class="line"><a name="l00219"></a><span class="lineno"> 219</span>  Neys_False_Easting = n.Neys_False_Easting; </div>
<div class="line"><a name="l00220"></a><span class="lineno"> 220</span>  Neys_False_Northing = n.Neys_False_Northing; </div>
<div class="line"><a name="l00221"></a><span class="lineno"> 221</span>  Neys_Delta_Easting = n.Neys_Delta_Easting; </div>
<div class="line"><a name="l00222"></a><span class="lineno"> 222</span>  Neys_Delta_Northing = n.Neys_Delta_Northing; </div>
<div class="line"><a name="l00223"></a><span class="lineno"> 223</span> }</div>
<div class="line"><a name="l00224"></a><span class="lineno"> 224</span> </div>
<div class="line"><a name="l00225"></a><span class="lineno"> 225</span> </div>
<div class="line"><a name="l00226"></a><span class="lineno"><a class="code" href="class_m_s_p_1_1_c_c_s_1_1_neys.html#a2d59628bf9e490dc6b8fff627747c8b4"> 226</a></span> <a class="code" href="class_m_s_p_1_1_c_c_s_1_1_neys.html#a2d59628bf9e490dc6b8fff627747c8b4">Neys::~Neys</a>()</div>
<div class="line"><a name="l00227"></a><span class="lineno"> 227</span> {</div>
<div class="line"><a name="l00228"></a><span class="lineno"> 228</span>  <span class="keyword">delete</span> lambertConformalConic2;</div>
<div class="line"><a name="l00229"></a><span class="lineno"> 229</span>  lambertConformalConic2 = 0;</div>
<div class="line"><a name="l00230"></a><span class="lineno"> 230</span> }</div>
<div class="line"><a name="l00231"></a><span class="lineno"> 231</span> </div>
<div class="line"><a name="l00232"></a><span class="lineno"> 232</span> </div>
<div class="line"><a name="l00233"></a><span class="lineno"><a class="code" href="class_m_s_p_1_1_c_c_s_1_1_neys.html#a17cdb82858b4264840131c93574a37b2"> 233</a></span> <a class="code" href="class_m_s_p_1_1_c_c_s_1_1_neys.html">Neys</a>& <a class="code" href="class_m_s_p_1_1_c_c_s_1_1_neys.html#a17cdb82858b4264840131c93574a37b2">Neys::operator=</a>( <span class="keyword">const</span> <a class="code" href="class_m_s_p_1_1_c_c_s_1_1_neys.html">Neys</a> &n )</div>
<div class="line"><a name="l00234"></a><span class="lineno"> 234</span> {</div>
<div class="line"><a name="l00235"></a><span class="lineno"> 235</span>  <span class="keywordflow">if</span>( <span class="keyword">this</span> != &n )</div>
<div class="line"><a name="l00236"></a><span class="lineno"> 236</span>  {</div>
<div class="line"><a name="l00237"></a><span class="lineno"> 237</span>  lambertConformalConic2->operator=( *n.lambertConformalConic2 );</div>
<div class="line"><a name="l00238"></a><span class="lineno"> 238</span>  <a class="code" href="class_m_s_p_1_1_c_c_s_1_1_coordinate_system.html#ab4b25504df8abc3fb23b5f03eaf6c14e">semiMajorAxis</a> = n.<a class="code" href="class_m_s_p_1_1_c_c_s_1_1_coordinate_system.html#ab4b25504df8abc3fb23b5f03eaf6c14e">semiMajorAxis</a>;</div>
<div class="line"><a name="l00239"></a><span class="lineno"> 239</span>  <a class="code" href="class_m_s_p_1_1_c_c_s_1_1_coordinate_system.html#a103abc45b36f0bc52bd992f2d19bbfc6">flattening</a> = n.<a class="code" href="class_m_s_p_1_1_c_c_s_1_1_coordinate_system.html#a103abc45b36f0bc52bd992f2d19bbfc6">flattening</a>;</div>
<div class="line"><a name="l00240"></a><span class="lineno"> 240</span>  Neys_Std_Parallel_1 = n.Neys_Std_Parallel_1; </div>
<div class="line"><a name="l00241"></a><span class="lineno"> 241</span>  Neys_Std_Parallel_2 = n.Neys_Std_Parallel_2; </div>
<div class="line"><a name="l00242"></a><span class="lineno"> 242</span>  Neys_Origin_Long = n.Neys_Origin_Long; </div>
<div class="line"><a name="l00243"></a><span class="lineno"> 243</span>  Neys_Origin_Lat = n.Neys_Origin_Lat; </div>
<div class="line"><a name="l00244"></a><span class="lineno"> 244</span>  Neys_False_Easting = n.Neys_False_Easting; </div>
<div class="line"><a name="l00245"></a><span class="lineno"> 245</span>  Neys_False_Northing = n.Neys_False_Northing; </div>
<div class="line"><a name="l00246"></a><span class="lineno"> 246</span>  Neys_Delta_Easting = n.Neys_Delta_Easting; </div>
<div class="line"><a name="l00247"></a><span class="lineno"> 247</span>  Neys_Delta_Northing = n.Neys_Delta_Northing; </div>
<div class="line"><a name="l00248"></a><span class="lineno"> 248</span>  }</div>
<div class="line"><a name="l00249"></a><span class="lineno"> 249</span> </div>
<div class="line"><a name="l00250"></a><span class="lineno"> 250</span>  <span class="keywordflow">return</span> *<span class="keyword">this</span>;</div>
<div class="line"><a name="l00251"></a><span class="lineno"> 251</span> }</div>
<div class="line"><a name="l00252"></a><span class="lineno"> 252</span> </div>
<div class="line"><a name="l00253"></a><span class="lineno"> 253</span> </div>
<div class="line"><a name="l00254"></a><span class="lineno"><a class="code" href="class_m_s_p_1_1_c_c_s_1_1_neys.html#a466ae200899b9b74f898799b11369aba"> 254</a></span> <a class="code" href="class_m_s_p_1_1_c_c_s_1_1_neys_parameters.html">NeysParameters</a>* <a class="code" href="class_m_s_p_1_1_c_c_s_1_1_neys.html#a466ae200899b9b74f898799b11369aba">Neys::getParameters</a>()<span class="keyword"> const</span></div>
<div class="line"><a name="l00255"></a><span class="lineno"> 255</span> <span class="keyword"></span>{</div>
<div class="line"><a name="l00256"></a><span class="lineno"> 256</span> <span class="comment">/*</span></div>
<div class="line"><a name="l00257"></a><span class="lineno"> 257</span> <span class="comment"> * The function getParameters returns the current ellipsoid</span></div>
<div class="line"><a name="l00258"></a><span class="lineno"> 258</span> <span class="comment"> * parameters and Ney's (Modified Lambert Conformal Conic) projection parameters.</span></div>
<div class="line"><a name="l00259"></a><span class="lineno"> 259</span> <span class="comment"> *</span></div>
<div class="line"><a name="l00260"></a><span class="lineno"> 260</span> <span class="comment"> * ellipsoidSemiMajorAxis : Semi-major axis of ellipsoid, in meters (output)</span></div>
<div class="line"><a name="l00261"></a><span class="lineno"> 261</span> <span class="comment"> * ellipsoidFlattening : Flattening of ellipsoid (output)</span></div>
<div class="line"><a name="l00262"></a><span class="lineno"> 262</span> <span class="comment"> * centralMeridian : Longitude of origin, in radians (output)</span></div>
<div class="line"><a name="l00263"></a><span class="lineno"> 263</span> <span class="comment"> * originLatitude : Latitude of origin, in radians (output)</span></div>
<div class="line"><a name="l00264"></a><span class="lineno"> 264</span> <span class="comment"> * standardParallel : First standard parallel, in radians (output)</span></div>
<div class="line"><a name="l00265"></a><span class="lineno"> 265</span> <span class="comment"> * falseEasting : False easting, in meters (output)</span></div>
<div class="line"><a name="l00266"></a><span class="lineno"> 266</span> <span class="comment"> * falseNorthing : False northing, in meters (output)</span></div>
<div class="line"><a name="l00267"></a><span class="lineno"> 267</span> <span class="comment"> */</span></div>
<div class="line"><a name="l00268"></a><span class="lineno"> 268</span> </div>
<div class="line"><a name="l00269"></a><span class="lineno"> 269</span>  <span class="keywordflow">return</span> <span class="keyword">new</span> <a class="code" href="class_m_s_p_1_1_c_c_s_1_1_neys_parameters.html">NeysParameters</a>( <a class="code" href="class_m_s_p_1_1_c_c_s_1_1_coordinate_type.html#a418b669819e988d0cb43c506beffd8ebadf6d2c4ee9e5612ca18ce33e8c640198">CoordinateType::neys</a>, Neys_Origin_Long, Neys_Origin_Lat, Neys_Std_Parallel_1, Neys_False_Easting, Neys_False_Northing );</div>
<div class="line"><a name="l00270"></a><span class="lineno"> 270</span> }</div>
<div class="line"><a name="l00271"></a><span class="lineno"> 271</span> </div>
<div class="line"><a name="l00272"></a><span class="lineno"> 272</span> </div>
<div class="line"><a name="l00273"></a><span class="lineno"><a class="code" href="class_m_s_p_1_1_c_c_s_1_1_neys.html#ad804f12f75e2436f157bc92dfb263442"> 273</a></span> <a class="code" href="class_m_s_p_1_1_c_c_s_1_1_map_projection_coordinates.html">MSP::CCS::MapProjectionCoordinates</a>* <a class="code" href="class_m_s_p_1_1_c_c_s_1_1_neys.html#ad804f12f75e2436f157bc92dfb263442">Neys::convertFromGeodetic</a>( <a class="code" href="class_m_s_p_1_1_c_c_s_1_1_geodetic_coordinates.html">MSP::CCS::GeodeticCoordinates</a>* geodeticCoordinates )</div>
<div class="line"><a name="l00274"></a><span class="lineno"> 274</span> {</div>
<div class="line"><a name="l00275"></a><span class="lineno"> 275</span> <span class="comment">/*</span></div>
<div class="line"><a name="l00276"></a><span class="lineno"> 276</span> <span class="comment"> * The function convertFromGeodetic converts Geodetic (latitude and</span></div>
<div class="line"><a name="l00277"></a><span class="lineno"> 277</span> <span class="comment"> * longitude) coordinates to Ney's (Modified Lambert Conformal Conic) projection</span></div>
<div class="line"><a name="l00278"></a><span class="lineno"> 278</span> <span class="comment"> * (easting and northing) coordinates, according to the current ellipsoid and</span></div>
<div class="line"><a name="l00279"></a><span class="lineno"> 279</span> <span class="comment"> * Ney's (Modified Lambert Conformal Conic) projection parameters. If any errors </span></div>
<div class="line"><a name="l00280"></a><span class="lineno"> 280</span> <span class="comment"> * occur, an exception is thrown with a description of the error.</span></div>
<div class="line"><a name="l00281"></a><span class="lineno"> 281</span> <span class="comment"> *</span></div>
<div class="line"><a name="l00282"></a><span class="lineno"> 282</span> <span class="comment"> * longitude : Longitude, in radians (input)</span></div>
<div class="line"><a name="l00283"></a><span class="lineno"> 283</span> <span class="comment"> * latitude : Latitude, in radians (input)</span></div>
<div class="line"><a name="l00284"></a><span class="lineno"> 284</span> <span class="comment"> * easting : Easting (X), in meters (output)</span></div>
<div class="line"><a name="l00285"></a><span class="lineno"> 285</span> <span class="comment"> * northing : Northing (Y), in meters (output)</span></div>
<div class="line"><a name="l00286"></a><span class="lineno"> 286</span> <span class="comment"> */</span></div>
<div class="line"><a name="l00287"></a><span class="lineno"> 287</span> </div>
<div class="line"><a name="l00288"></a><span class="lineno"> 288</span>  <span class="keywordtype">double</span> longitude = geodeticCoordinates-><a class="code" href="class_m_s_p_1_1_c_c_s_1_1_geodetic_coordinates.html#a81b352849684e9910b4001be3fec1b0d">longitude</a>();</div>
<div class="line"><a name="l00289"></a><span class="lineno"> 289</span>  <span class="keywordtype">double</span> latitude = geodeticCoordinates-><a class="code" href="class_m_s_p_1_1_c_c_s_1_1_geodetic_coordinates.html#a4214a0369f1e642e1a446f35b7151108">latitude</a>();</div>
<div class="line"><a name="l00290"></a><span class="lineno"> 290</span> </div>
<div class="line"><a name="l00291"></a><span class="lineno"> 291</span>  <span class="keywordflow">if</span> ((latitude < -<a class="code" href="_m_g_r_s_8cpp.html#a27d2bea15134ae7c7f7a2d11379bde53">PI_OVER_2</a>) || (latitude > <a class="code" href="_m_g_r_s_8cpp.html#a27d2bea15134ae7c7f7a2d11379bde53">PI_OVER_2</a>))</div>
<div class="line"><a name="l00292"></a><span class="lineno"> 292</span>  { <span class="comment">/* Latitude out of range */</span></div>
<div class="line"><a name="l00293"></a><span class="lineno"> 293</span>  <span class="keywordflow">throw</span> <a class="code" href="class_m_s_p_1_1_c_c_s_1_1_coordinate_conversion_exception.html">CoordinateConversionException</a>( <a class="code" href="class_m_s_p_1_1_c_c_s_1_1_error_messages.html#a816b9dcb83273fbc1ba1f3c215f26329">ErrorMessages::latitude</a> );</div>
<div class="line"><a name="l00294"></a><span class="lineno"> 294</span>  }</div>
<div class="line"><a name="l00295"></a><span class="lineno"> 295</span>  <span class="keywordflow">if</span> ((longitude < -<a class="code" href="_m_g_r_s_8cpp.html#a598a3330b3c21701223ee0ca14316eca">PI</a>) || (longitude > <a class="code" href="_albers_equal_area_conic_8cpp.html#a820597b124334bf7fb85214114f4876d">TWO_PI</a>))</div>
<div class="line"><a name="l00296"></a><span class="lineno"> 296</span>  { <span class="comment">/* Longitude out of range */</span></div>
<div class="line"><a name="l00297"></a><span class="lineno"> 297</span>  <span class="keywordflow">throw</span> <a class="code" href="class_m_s_p_1_1_c_c_s_1_1_coordinate_conversion_exception.html">CoordinateConversionException</a>( <a class="code" href="class_m_s_p_1_1_c_c_s_1_1_error_messages.html#aca84764ebfd6a968bae4726c33c4c20b">ErrorMessages::longitude</a> );</div>
<div class="line"><a name="l00298"></a><span class="lineno"> 298</span>  }</div>
<div class="line"><a name="l00299"></a><span class="lineno"> 299</span> </div>
<div class="line"><a name="l00300"></a><span class="lineno"> 300</span>  <a class="code" href="class_m_s_p_1_1_c_c_s_1_1_map_projection_coordinates.html">MapProjectionCoordinates</a>* mapProjectionCoordinates =</div>
<div class="line"><a name="l00301"></a><span class="lineno"> 301</span>  lambertConformalConic2-><a class="code" href="class_m_s_p_1_1_c_c_s_1_1_lambert_conformal_conic.html#a2130cbafeb477aba15499fad2d2be60c">convertFromGeodetic</a>( geodeticCoordinates );</div>
<div class="line"><a name="l00302"></a><span class="lineno"> 302</span>  <span class="keywordtype">double</span> easting = mapProjectionCoordinates-><a class="code" href="class_m_s_p_1_1_c_c_s_1_1_map_projection_coordinates.html#a96f2d51e11a3b0b7175e220f9617b16e">easting</a>();</div>
<div class="line"><a name="l00303"></a><span class="lineno"> 303</span>  <span class="keywordtype">double</span> northing = mapProjectionCoordinates-><a class="code" href="class_m_s_p_1_1_c_c_s_1_1_map_projection_coordinates.html#aea94abfb2ad8c80adbbeb7ad1de08720">northing</a>();</div>
<div class="line"><a name="l00304"></a><span class="lineno"> 304</span>  <span class="keyword">delete</span> mapProjectionCoordinates;</div>
<div class="line"><a name="l00305"></a><span class="lineno"> 305</span> </div>
<div class="line"><a name="l00306"></a><span class="lineno"> 306</span>  <span class="keywordflow">return</span> <span class="keyword">new</span> <a class="code" href="class_m_s_p_1_1_c_c_s_1_1_map_projection_coordinates.html">MapProjectionCoordinates</a>( <a class="code" href="class_m_s_p_1_1_c_c_s_1_1_coordinate_type.html#a418b669819e988d0cb43c506beffd8ebadf6d2c4ee9e5612ca18ce33e8c640198">CoordinateType::neys</a>, easting, northing );</div>
<div class="line"><a name="l00307"></a><span class="lineno"> 307</span> }</div>
<div class="line"><a name="l00308"></a><span class="lineno"> 308</span> </div>
<div class="line"><a name="l00309"></a><span class="lineno"> 309</span> </div>
<div class="line"><a name="l00310"></a><span class="lineno"><a class="code" href="class_m_s_p_1_1_c_c_s_1_1_neys.html#aa29f8adb46ab8ee39c5269dabe3e33ea"> 310</a></span> <a class="code" href="class_m_s_p_1_1_c_c_s_1_1_geodetic_coordinates.html">MSP::CCS::GeodeticCoordinates</a>* <a class="code" href="class_m_s_p_1_1_c_c_s_1_1_neys.html#aa29f8adb46ab8ee39c5269dabe3e33ea">Neys::convertToGeodetic</a>(</div>
<div class="line"><a name="l00311"></a><span class="lineno"> 311</span>  <a class="code" href="class_m_s_p_1_1_c_c_s_1_1_map_projection_coordinates.html">MSP::CCS::MapProjectionCoordinates</a>* mapProjectionCoordinates )</div>
<div class="line"><a name="l00312"></a><span class="lineno"> 312</span> {</div>
<div class="line"><a name="l00313"></a><span class="lineno"> 313</span> <span class="comment">/*</span></div>
<div class="line"><a name="l00314"></a><span class="lineno"> 314</span> <span class="comment"> * The function convertToGeodetic converts Ney's (Modified Lambert Conformal</span></div>
<div class="line"><a name="l00315"></a><span class="lineno"> 315</span> <span class="comment"> * Conic) projection (easting and northing) coordinates to Geodetic (latitude)</span></div>
<div class="line"><a name="l00316"></a><span class="lineno"> 316</span> <span class="comment"> * and longitude) coordinates, according to the current ellipsoid and Ney's</span></div>
<div class="line"><a name="l00317"></a><span class="lineno"> 317</span> <span class="comment"> * (Modified Lambert Conformal Conic) projection parameters. If any errors </span></div>
<div class="line"><a name="l00318"></a><span class="lineno"> 318</span> <span class="comment"> * occur, an exception is thrown with a description of the error.</span></div>
<div class="line"><a name="l00319"></a><span class="lineno"> 319</span> <span class="comment"> *</span></div>
<div class="line"><a name="l00320"></a><span class="lineno"> 320</span> <span class="comment"> * easting : Easting (X), in meters (input)</span></div>
<div class="line"><a name="l00321"></a><span class="lineno"> 321</span> <span class="comment"> * northing : Northing (Y), in meters (input)</span></div>
<div class="line"><a name="l00322"></a><span class="lineno"> 322</span> <span class="comment"> * longitude : Longitude, in radians (output)</span></div>
<div class="line"><a name="l00323"></a><span class="lineno"> 323</span> <span class="comment"> * latitude : Latitude, in radians (output)</span></div>
<div class="line"><a name="l00324"></a><span class="lineno"> 324</span> <span class="comment"> */</span></div>
<div class="line"><a name="l00325"></a><span class="lineno"> 325</span> </div>
<div class="line"><a name="l00326"></a><span class="lineno"> 326</span>  <span class="keywordtype">double</span> easting = mapProjectionCoordinates-><a class="code" href="class_m_s_p_1_1_c_c_s_1_1_map_projection_coordinates.html#a96f2d51e11a3b0b7175e220f9617b16e">easting</a>();</div>
<div class="line"><a name="l00327"></a><span class="lineno"> 327</span>  <span class="keywordtype">double</span> northing = mapProjectionCoordinates-><a class="code" href="class_m_s_p_1_1_c_c_s_1_1_map_projection_coordinates.html#aea94abfb2ad8c80adbbeb7ad1de08720">northing</a>();</div>
<div class="line"><a name="l00328"></a><span class="lineno"> 328</span> </div>
<div class="line"><a name="l00329"></a><span class="lineno"> 329</span>  <span class="keywordflow">if</span> ((easting < (Neys_False_Easting - Neys_Delta_Easting))</div>
<div class="line"><a name="l00330"></a><span class="lineno"> 330</span>  ||(easting > (Neys_False_Easting + Neys_Delta_Easting)))</div>
<div class="line"><a name="l00331"></a><span class="lineno"> 331</span>  { <span class="comment">/* Easting out of range */</span></div>
<div class="line"><a name="l00332"></a><span class="lineno"> 332</span>  <span class="keywordflow">throw</span> <a class="code" href="class_m_s_p_1_1_c_c_s_1_1_coordinate_conversion_exception.html">CoordinateConversionException</a>( <a class="code" href="class_m_s_p_1_1_c_c_s_1_1_error_messages.html#a35b927051eb2143a2fda20693abd3290">ErrorMessages::easting</a> );</div>
<div class="line"><a name="l00333"></a><span class="lineno"> 333</span>  }</div>
<div class="line"><a name="l00334"></a><span class="lineno"> 334</span>  <span class="keywordflow">if</span> ((northing < (Neys_False_Northing - Neys_Delta_Northing))</div>
<div class="line"><a name="l00335"></a><span class="lineno"> 335</span>  || (northing > (Neys_False_Northing + Neys_Delta_Northing)))</div>
<div class="line"><a name="l00336"></a><span class="lineno"> 336</span>  { <span class="comment">/* Northing out of range */</span></div>
<div class="line"><a name="l00337"></a><span class="lineno"> 337</span>  <span class="keywordflow">throw</span> <a class="code" href="class_m_s_p_1_1_c_c_s_1_1_coordinate_conversion_exception.html">CoordinateConversionException</a>( <a class="code" href="class_m_s_p_1_1_c_c_s_1_1_error_messages.html#a31cbb504cb4705337394cdda99b6279c">ErrorMessages::northing</a> );</div>
<div class="line"><a name="l00338"></a><span class="lineno"> 338</span>  }</div>
<div class="line"><a name="l00339"></a><span class="lineno"> 339</span> </div>
<div class="line"><a name="l00340"></a><span class="lineno"> 340</span>  <a class="code" href="class_m_s_p_1_1_c_c_s_1_1_geodetic_coordinates.html">GeodeticCoordinates</a>* geodeticCoordinates = </div>
<div class="line"><a name="l00341"></a><span class="lineno"> 341</span>  lambertConformalConic2-><a class="code" href="class_m_s_p_1_1_c_c_s_1_1_lambert_conformal_conic.html#a2e9349aa9c372824052797eb03126c58">convertToGeodetic</a>( mapProjectionCoordinates );</div>
<div class="line"><a name="l00342"></a><span class="lineno"> 342</span>  <span class="keywordtype">double</span> longitude = geodeticCoordinates-><a class="code" href="class_m_s_p_1_1_c_c_s_1_1_geodetic_coordinates.html#a81b352849684e9910b4001be3fec1b0d">longitude</a>();</div>
<div class="line"><a name="l00343"></a><span class="lineno"> 343</span>  <span class="keywordtype">double</span> latitude = geodeticCoordinates-><a class="code" href="class_m_s_p_1_1_c_c_s_1_1_geodetic_coordinates.html#a4214a0369f1e642e1a446f35b7151108">latitude</a>();</div>
<div class="line"><a name="l00344"></a><span class="lineno"> 344</span>  <span class="keyword">delete</span> geodeticCoordinates;</div>
<div class="line"><a name="l00345"></a><span class="lineno"> 345</span> </div>
<div class="line"><a name="l00346"></a><span class="lineno"> 346</span>  <span class="keywordflow">return</span> <span class="keyword">new</span> <a class="code" href="class_m_s_p_1_1_c_c_s_1_1_geodetic_coordinates.html">GeodeticCoordinates</a>( <a class="code" href="class_m_s_p_1_1_c_c_s_1_1_coordinate_type.html#a418b669819e988d0cb43c506beffd8eba49202f41636b7629fc9431e8d1b448d4">CoordinateType::geodetic</a>, longitude, latitude );</div>
<div class="line"><a name="l00347"></a><span class="lineno"> 347</span> }</div>
<div class="line"><a name="l00348"></a><span class="lineno"> 348</span> </div>
<div class="line"><a name="l00349"></a><span class="lineno"> 349</span> </div>
<div class="line"><a name="l00350"></a><span class="lineno"> 350</span> </div>
<div class="line"><a name="l00351"></a><span class="lineno"> 351</span> <span class="comment">// CLASSIFICATION: UNCLASSIFIED</span></div>
</div><!-- fragment --></div><!-- contents -->
<!-- start footer part -->
<hr class="footer"/><address class="footer"><small>
Generated on Tue Feb 16 2016 14:54:02 for GeographicTranslator by doxygen 1.8.2
</small></address>
</body>
</html>
|