/* +---------------------------------------------------------------------------+
| The Mobile Robot Programming Toolkit (MRPT) C++ library |
| |
| http://www.mrpt.org/ |
| |
| Copyright (C) 2005-2011 University of Malaga |
| |
| This software was written by the Machine Perception and Intelligent |
| Robotics Lab, University of Malaga (Spain). |
| Contact: Jose-Luis Blanco <jlblanco@ctima.uma.es> |
| |
| This file is part of the MRPT project. |
| |
| MRPT is free software: you can redistribute it and/or modify |
| it under the terms of the GNU General Public License as published by |
| the Free Software Foundation, either version 3 of the License, or |
| (at your option) any later version. |
| |
| MRPT is distributed in the hope that it will be useful, |
| but WITHOUT ANY WARRANTY; without even the implied warranty of |
| MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the |
| GNU General Public License for more details. |
| |
| You should have received a copy of the GNU General Public License |
| along with MRPT. If not, see <http://www.gnu.org/licenses/>. |
| |
+---------------------------------------------------------------------------+ */
#include <mrpt/base.h>
#include <mrpt/topography.h>
using namespace mrpt::utils;
using namespace mrpt::math;
using namespace mrpt::topography;
using namespace std;
#include <mrpt/examples_config.h>
std::vector<string> names;
std::vector<CMatrixD> results;
const double TH = 0.02;
void exampleResults()
{
names.resize( 6 ); names[0] = "MAND"; names[1] = "OROZ"; names[2] = "SANC"; names[3] = "ARJO"; names[4] = "CANI"; names[5] = "HALC";
results.resize( 25 );
// Example 01
results[0] = CMatrixD(6,3);
results[0].set_unsafe(0,0,435699.5503); results[0].set_unsafe(0,1,4194142.0571); results[0].set_unsafe(0,2,515.5790);
results[0].set_unsafe(1,0,443948.2295); results[0].set_unsafe(1,1,4160754.7114); results[0].set_unsafe(1,2,1489.1880);
results[0].set_unsafe(2,0,436542.0710); results[0].set_unsafe(2,1,4177367.6694); results[0].set_unsafe(2,2,1054.1890);
results[0].set_unsafe(3,0,406284.9288); results[0].set_unsafe(3,1,4176492.3033); results[0].set_unsafe(3,2,569.1820);
results[0].set_unsafe(4,0,406546.3268); results[0].set_unsafe(4,1,4156628.3641); results[0].set_unsafe(4,2,897.6260);
results[0].set_unsafe(5,0,408051.6599); results[0].set_unsafe(5,1,4193733.3066); results[0].set_unsafe(5,2,530.4460);
// Example 02
results[1] = CMatrixD(6,3);
results[1].set_unsafe(0,0,435810.797); results[1].set_unsafe(0,1,4194348.753); results[1].set_unsafe(0,2,437.516);
results[1].set_unsafe(1,0,444059.747); results[1].set_unsafe(1,1,4160961.358); results[1].set_unsafe(1,2,1410.518);
results[1].set_unsafe(2,0,436653.444); results[1].set_unsafe(2,1,4177574.314); results[1].set_unsafe(2,2,975.764);
results[1].set_unsafe(3,0,406396.235); results[1].set_unsafe(3,1,4176698.718); results[1].set_unsafe(3,2,490.219);
results[1].set_unsafe(4,0,406657.783); results[1].set_unsafe(4,1,4156834.717); results[1].set_unsafe(4,2,818.224);
results[1].set_unsafe(5,0,408162.836); results[1].set_unsafe(5,1,4193939.783); results[1].set_unsafe(5,2,451.896);
// Example 03
results[2] = CMatrixD(6,3);
results[2].set_unsafe(0,0,435811.252); results[2].set_unsafe(0,1,4194347.983); results[2].set_unsafe(0,2,515.579);
results[2].set_unsafe(1,0,444059.929); results[2].set_unsafe(1,1,4160960.606); results[2].set_unsafe(1,2,1489.188);
results[2].set_unsafe(2,0,436653.768); results[2].set_unsafe(2,1,4177573.581); results[2].set_unsafe(2,2,1054.189);
results[2].set_unsafe(3,0,406396.600); results[2].set_unsafe(3,1,4176698.222); results[2].set_unsafe(3,2,569.182);
results[2].set_unsafe(4,0,406657.992); results[2].set_unsafe(4,1,4156834.266); results[2].set_unsafe(4,2,897.626);
results[2].set_unsafe(5,0,408163.337); results[2].set_unsafe(5,1,4193939.241); results[2].set_unsafe(5,2,530.446);
// Example 04
results[3] = CMatrixD(6,3);
results[3].set_unsafe(0,0,435811.073); results[3].set_unsafe(0,1,4194347.990); results[3].set_unsafe(0,2,466.175);
results[3].set_unsafe(1,0,444059.962); results[3].set_unsafe(1,1,4160960.412); results[3].set_unsafe(1,2,1439.784);
results[3].set_unsafe(2,0,436653.708); results[3].set_unsafe(2,1,4177573.506); results[3].set_unsafe(2,2,1004.785);
results[3].set_unsafe(3,0,406396.665); results[3].set_unsafe(3,1,4176698.318); results[3].set_unsafe(3,2,519.778);
results[3].set_unsafe(4,0,406658.202); results[3].set_unsafe(4,1,4156834.269); results[3].set_unsafe(4,2,848.222);
results[3].set_unsafe(5,0,408163.269); results[3].set_unsafe(5,1,4193939.405); results[3].set_unsafe(5,2,481.042);
// Example 05
results[4] = CMatrixD(6,3);
results[4].set_unsafe(0,0,435811.251); results[4].set_unsafe(0,1,4194347.974); results[4].set_unsafe(0,2,465.655);
results[4].set_unsafe(1,0,444059.934); results[4].set_unsafe(1,1,4160960.612); results[4].set_unsafe(1,2,1439.740);
results[4].set_unsafe(2,0,436653.771); results[4].set_unsafe(2,1,4177573.582); results[4].set_unsafe(2,2,1004.556);
results[4].set_unsafe(3,0,406396.596); results[4].set_unsafe(3,1,4176698.222); results[4].set_unsafe(3,2,520.037);
results[4].set_unsafe(4,0,406657.988); results[4].set_unsafe(4,1,4156834.272); results[4].set_unsafe(4,2,848.836);
results[4].set_unsafe(5,0,408163.337); results[4].set_unsafe(5,1,4193939.239); results[4].set_unsafe(5,2,480.961);
// Example 06
results[5] = CMatrixD(6,3);
results[5].set_unsafe(0,0,435811.250); results[5].set_unsafe(0,1,4194347.994); results[5].set_unsafe(0,2,437.516);
results[5].set_unsafe(1,0,444059.921); results[5].set_unsafe(1,1,4160960.599); results[5].set_unsafe(1,2,1410.518);
results[5].set_unsafe(2,0,436653.764); results[5].set_unsafe(2,1,4177573.581); results[5].set_unsafe(2,2,975.764);
results[5].set_unsafe(3,0,406396.606); results[5].set_unsafe(3,1,4176698.225); results[5].set_unsafe(3,2,490.219);
results[5].set_unsafe(4,0,406657.997); results[5].set_unsafe(4,1,4156834.260); results[5].set_unsafe(4,2,818.224);
results[5].set_unsafe(5,0,408163.340); results[5].set_unsafe(5,1,4193939.242); results[5].set_unsafe(5,2,451.896);
// Example 07
results[6] = CMatrixD(6,3);
results[6].set_unsafe(0,0,435811.074); results[6].set_unsafe(0,1,4194347.994); results[6].set_unsafe(0,2,466.791);
results[6].set_unsafe(1,0,444059.961); results[6].set_unsafe(1,1,4160960.410); results[6].set_unsafe(1,2,1439.793);
results[6].set_unsafe(2,0,436653.708); results[6].set_unsafe(2,1,4177573.505); results[6].set_unsafe(2,2,1005.039);
results[6].set_unsafe(3,0,406396.667); results[6].set_unsafe(3,1,4176698.320); results[6].set_unsafe(3,2,519.494);
results[6].set_unsafe(4,0,406658.203); results[6].set_unsafe(4,1,4156834.270); results[6].set_unsafe(4,2,847.499);
results[6].set_unsafe(5,0,408163.267); results[6].set_unsafe(5,1,4193939.401); results[6].set_unsafe(5,2,481.171);
// Example 08
results[7] = CMatrixD(6,3);
results[7].set_unsafe(0,0,435811.250); results[7].set_unsafe(0,1,4194347.974); results[7].set_unsafe(0,2,465.657);
results[7].set_unsafe(1,0,444059.935); results[7].set_unsafe(1,1,4160960.612); results[7].set_unsafe(1,2,1439.739);
results[7].set_unsafe(2,0,436653.772); results[7].set_unsafe(2,1,4177573.582); results[7].set_unsafe(2,2,1004.556);
results[7].set_unsafe(3,0,406396.597); results[7].set_unsafe(3,1,4176698.223); results[7].set_unsafe(3,2,520.037);
results[7].set_unsafe(4,0,406657.988); results[7].set_unsafe(4,1,4156834.271); results[7].set_unsafe(4,2,848.837);
results[7].set_unsafe(5,0,408163.339); results[7].set_unsafe(5,1,4193939.239); results[7].set_unsafe(5,2,480.959);
// Example 10
results[9] = CMatrixD(6,3);
results[9].set_unsafe(0,0,435811.2512); results[9].set_unsafe(0,1,4194347.9756); results[9].set_unsafe(0,2,465.6558);
results[9].set_unsafe(1,0,444059.9337); results[9].set_unsafe(1,1,4160960.6109); results[9].set_unsafe(1,2,1439.7390);
results[9].set_unsafe(2,0,436653.7721); results[9].set_unsafe(2,1,4177573.5812); results[9].set_unsafe(2,2,1004.5560);
results[9].set_unsafe(3,0,406396.5970); results[9].set_unsafe(3,1,4176698.2216); results[9].set_unsafe(3,2,520.0376);
results[9].set_unsafe(4,0,406657.9894); results[9].set_unsafe(4,1,4156834.2734); results[9].set_unsafe(4,2,848.8373);
results[9].set_unsafe(5,0,408163.3366); results[9].set_unsafe(5,1,4193939.2373); results[9].set_unsafe(5,2,480.9593);
// Example 11
results[10] = CMatrixD(6,3);
results[10].set_unsafe(0,0,435811.2512); results[10].set_unsafe(0,1,4194347.9756); results[10].set_unsafe(0,2,465.6558);
results[10].set_unsafe(1,0,444059.9337); results[10].set_unsafe(1,1,4160960.6109); results[10].set_unsafe(1,2,1439.7390);
results[10].set_unsafe(2,0,436653.7721); results[10].set_unsafe(2,1,4177573.5812); results[10].set_unsafe(2,2,1004.5560);
results[10].set_unsafe(3,0,406396.5970); results[10].set_unsafe(3,1,4176698.2216); results[10].set_unsafe(3,2,520.0376);
results[10].set_unsafe(4,0,406657.9894); results[10].set_unsafe(4,1,4156834.2734); results[10].set_unsafe(4,2,848.8373);
results[10].set_unsafe(5,0,408163.3366); results[10].set_unsafe(5,1,4193939.2373); results[10].set_unsafe(5,2,480.9593);
// Example 12
results[11] = CMatrixD(6,3);
results[11].set_unsafe(0,0,435811.156); results[11].set_unsafe(0,1,4194347.892); results[11].set_unsafe(0,2,466.100);
results[11].set_unsafe(1,0,444059.719); results[11].set_unsafe(1,1,4160960.489); results[11].set_unsafe(1,2,1439.389);
results[11].set_unsafe(2,0,436653.606); results[11].set_unsafe(2,1,4177573.398); results[11].set_unsafe(2,2,1004.594);
results[11].set_unsafe(3,0,406396.045); results[11].set_unsafe(3,1,4176697.926); results[11].set_unsafe(3,2,520.238);
results[11].set_unsafe(4,0,406658.812); results[11].set_unsafe(4,1,4156834.183); results[11].set_unsafe(4,2,849.050);
results[11].set_unsafe(5,0,408163.542); results[11].set_unsafe(5,1,4193940.012); results[11].set_unsafe(5,2,480.414);
// Example 13
results[12] = CMatrixD(6,3);
results[12].set_unsafe(0,0,435811.095); results[12].set_unsafe(0,1,4194348.146); results[12].set_unsafe(0,2,465.884);
results[12].set_unsafe(1,0,444059.906); results[12].set_unsafe(1,1,4160960.276); results[12].set_unsafe(1,2,1439.564);
results[12].set_unsafe(2,0,436653.748); results[12].set_unsafe(2,1,4177573.521); results[12].set_unsafe(2,2,1004.572);
results[12].set_unsafe(3,0,406396.375); results[12].set_unsafe(3,1,4176698.033); results[12].set_unsafe(3,2,520.147);
results[12].set_unsafe(4,0,406658.370); results[12].set_unsafe(4,1,4156834.517); results[12].set_unsafe(4,2,848.945);
results[12].set_unsafe(5,0,408163.386); results[12].set_unsafe(5,1,4193939.406); results[12].set_unsafe(5,2,480.674);
}
void displayResults( const vector_double &thisResults, const unsigned int &example, const bool color )
{
cout << "Example " << example << ". Results" << setprecision(16) << endl;
cout << "----------------------------------------------------------------------------" << endl;
cout << "N\t" << "X\t\t\t" << "Y\t\t\t" << "Z" << endl;
cout << "----------------------------------------------------------------------------" << endl;
if( color )
{
const size_t N = thisResults.size()/3;
for( unsigned int i = 0; i < N; ++i )
{
cout << names[i] << "\t";
if( fabs( thisResults[3*i+0]-results[example-1].get_unsafe(i,0)) < TH )
mrpt::system::setConsoleColor(CONCOL_GREEN);
else
mrpt::system::setConsoleColor(CONCOL_RED);
cout << thisResults[3*i+0] << "\t";
if( fabs( thisResults[3*i+1]-results[example-1].get_unsafe(i,1)) < TH )
mrpt::system::setConsoleColor(CONCOL_GREEN);
else
mrpt::system::setConsoleColor(CONCOL_RED);
cout << thisResults[3*i+1] << "\t";
if( fabs( thisResults[3*i+2]-results[example-1].get_unsafe(i,2)) < TH )
mrpt::system::setConsoleColor(CONCOL_GREEN);
else
mrpt::system::setConsoleColor(CONCOL_RED);
cout << thisResults[3*i+2] << endl;
mrpt::system::setConsoleColor(CONCOL_NORMAL);
}
cout << endl;
}
else
{
const size_t N = thisResults.size()/3;
for( unsigned int i = 0; i < N; ++i )
cout << names[i] << "\t" << thisResults[3*i+0] << "\t" << thisResults[3*i+1] << "\t" << thisResults[3*i+2] << endl;
cout << endl;
}
}
void TestCoordinatesConversions()
{
// Initial UTM coordinates (X,Y,Z)
mrpt::math::TPoint3D UTMCoords;
UTMCoords.x = 435500;
UTMCoords.y = 4194142.057;
UTMCoords.z = 515.579;
const int huso = 30;
const char hem = 'n';
TGeodeticCoords GeodeticCoords;
TEllipsoid ellip = TEllipsoid::Ellipsoid_WGS84();
// Dump to console of input values
cout << "INPUT DATA: " << endl << "----------------------------------" << endl;
cout << "UTM Coordinates: " << endl;
cout << "X = " << setprecision(15) << UTMCoords.x << endl;
cout << "Y = " << UTMCoords.y << endl;
cout << "Z = " << UTMCoords.z << endl;
cout << "Time zone: " << huso << endl;
cout << "Hemisphery: " << hem << endl;
cout << "Ellipsoid: " << ellip.name << endl;
cout << "----------------------------------" << endl << endl;
// UTM to Geodetic ...
UTMToGeodetic( UTMCoords, huso, hem, GeodeticCoords, ellip );
cout << "UTM to Geodetic" << endl;
cout << "Geodetic Coordinates:" << endl;
cout << "Lon = " << TCoords( GeodeticCoords.lon ) << " [" << GeodeticCoords.lon << "]" << endl;
cout << "Lat = " << TCoords( GeodeticCoords.lat ) << " [" << GeodeticCoords.lon << "]" << endl;
cout << "H = " << GeodeticCoords.height << endl;
cout << "----------------------------------" << endl;
// Geodetic to Geocentric ...
TGeocentricCoords GeocentricCoords;
geodeticToGeocentric( GeodeticCoords, GeocentricCoords, ellip );
cout << "Geodetic to Geocentric" << endl;
cout << "Geocentric Coordinates:" << endl;
cout << "X = " << GeocentricCoords.x << endl;
cout << "Y = " << GeocentricCoords.y << endl;
cout << "Z = " << GeocentricCoords.z << endl;
cout << "----------------------------------" << endl;
// 7 parameter transformation (X,Y,Z) -> (X',Y',Z')
// TDatum7Params datum( dx, dy, dz, rx(sec), ry(sec), rz(sec), ds(ppm) );
TDatum7Params datum( 109.4714, 106.7407, 141.7916, 4.7834, 7.9668, -5.3771, -6.690000 );
TGeocentricCoords GeocentricCoords2;
transform7params( GeocentricCoords, datum, GeocentricCoords2 );
//transform7params( TGeocentricCoords(5029475.945,-328201.0396,3896351.728), datum, GeocentricCoords2 );
cout << "transform7params" << endl;
cout << "Transformed Geocentric Coordinates:" << endl;
cout << "X = " << GeocentricCoords2.x << endl;
cout << "Y = " << GeocentricCoords2.y << endl;
cout << "Z = " << GeocentricCoords2.z << endl;
cout << "----------------------------------" << endl;
// ... and back to Geodetic
TGeodeticCoords GeodeticCoords2;
ellip = TEllipsoid::Ellipsoid_Hough_1960();
geocentricToGeodetic( GeocentricCoords2, GeodeticCoords2, ellip );
cout << "Geocentric to Geodetic" << endl;
cout << "Geodetic Coordinates:" << endl;
cout << "Lon = " << TCoords( GeodeticCoords2.lon ) << " [" << GeodeticCoords2.lon << "]" << endl;
cout << "Lat = " << TCoords( GeodeticCoords2.lat ) << " [" << GeodeticCoords2.lat << "]" << endl;
cout << "H = " << GeodeticCoords2.height << endl;
cout << "----------------------------------" << endl;
// ... and back to UTM
TUTMCoords UTMCoords2;
int time_zone2;
char latitude_band2;
//geodeticToUTM( GeodeticCoords2, UTMCoords2, time_zone2, latitude_band2, ellip );
geodeticToUTM( TGeodeticCoords(37.89604181,-3.72987289,542.8624741), UTMCoords2, time_zone2, latitude_band2, TEllipsoid(6378270,6356794.343,"USER") );
cout << "Geodetic to UTM" << endl;
cout << "UTM Coordinates:" << endl;
cout << "X = " << UTMCoords2.x << endl;
cout << "Y = " << UTMCoords2.y << endl;
cout << "Z = " << UTMCoords2.z << endl;
cout << "Time zone: " << time_zone2 << endl;
cout << "----------------------------------" << endl;
// 10 parameter transformation (X,Y,Z) -> (X',Y',Z')
// Helmert 2D transformation (X,Y) -> (X',Y')
// Helmert 3D transformation (X,Y,Z) -> (X',Y',Z')
// 1D transformation (X,Y,Z) -> (X')
// Interpolation (X,Y) -> (X',Y')
}
void Examples_01()
{
TGeodeticCoords MAND_92633( TCoords(37,53,33.072573).getDecimalValue(), TCoords(-3,43,52.68965).getDecimalValue(), 515.579 );
TGeodeticCoords OROZ_96965( TCoords(37,35,31.75505).getDecimalValue(), TCoords(-3,38,5.70754).getDecimalValue(), 1489.188 );
TGeodeticCoords SANC_94744( TCoords(37,44,29.04253).getDecimalValue(), TCoords(-3,43,12.90285).getDecimalValue(), 1054.189 );
TGeodeticCoords ARJO_94633( TCoords(37,43,51.28918).getDecimalValue(), TCoords(-4,3,48.64503).getDecimalValue(), 569.182 );
TGeodeticCoords CANI_96833( TCoords(37,33,6.94237).getDecimalValue(), TCoords(-4,3,28.81611).getDecimalValue(), 897.626 );
TGeodeticCoords HALC_92543( TCoords(37,53,11.26838).getDecimalValue(), TCoords(-4,2,44.35794).getDecimalValue(), 530.446 );
// A2
int UTMZone;
char UTMBand;
TUTMCoords MAND_UTM, OROZ_UTM, SANC_UTM, ARJO_UTM, CANI_UTM, HALC_UTM;
geodeticToUTM( MAND_92633, MAND_UTM, UTMZone, UTMBand, TEllipsoid::Ellipsoid_WGS84() );
geodeticToUTM( OROZ_96965, OROZ_UTM, UTMZone, UTMBand, TEllipsoid::Ellipsoid_WGS84() );
geodeticToUTM( SANC_94744, SANC_UTM, UTMZone, UTMBand, TEllipsoid::Ellipsoid_WGS84() );
geodeticToUTM( ARJO_94633, ARJO_UTM, UTMZone, UTMBand, TEllipsoid::Ellipsoid_WGS84() );
geodeticToUTM( CANI_96833, CANI_UTM, UTMZone, UTMBand, TEllipsoid::Ellipsoid_WGS84() );
geodeticToUTM( HALC_92543, HALC_UTM, UTMZone, UTMBand, TEllipsoid::Ellipsoid_WGS84() );
vector_double thisResults;
thisResults.resize(18);
thisResults[0] = MAND_UTM.x; thisResults[1] = MAND_UTM.y; thisResults[2] = MAND_UTM.z;
thisResults[3] = OROZ_UTM.x; thisResults[4] = OROZ_UTM.y; thisResults[5] = OROZ_UTM.z;
thisResults[6] = SANC_UTM.x; thisResults[7] = SANC_UTM.y; thisResults[8] = SANC_UTM.z;
thisResults[9] = ARJO_UTM.x; thisResults[10] = ARJO_UTM.y; thisResults[11] = ARJO_UTM.z;
thisResults[12] = CANI_UTM.x; thisResults[13] = CANI_UTM.y; thisResults[14] = CANI_UTM.z;
thisResults[15] = HALC_UTM.x; thisResults[16] = HALC_UTM.y; thisResults[17] = HALC_UTM.z;
displayResults( thisResults, 1, true );
}
void Examples_02()
{
TGeodeticCoords MAND_92633( TCoords(37,53,33.072573).getDecimalValue(), TCoords(-3,43,52.68965).getDecimalValue(), 515.579 );
TGeodeticCoords OROZ_96965( TCoords(37,35,31.75505).getDecimalValue(), TCoords(-3,38,5.70754).getDecimalValue(), 1489.188 );
TGeodeticCoords SANC_94744( TCoords(37,44,29.04253).getDecimalValue(), TCoords(-3,43,12.90285).getDecimalValue(), 1054.189 );
TGeodeticCoords ARJO_94633( TCoords(37,43,51.28918).getDecimalValue(), TCoords(-4,3,48.64503).getDecimalValue(), 569.182 );
TGeodeticCoords CANI_96833( TCoords(37,33,6.94237).getDecimalValue(), TCoords(-4,3,28.81611).getDecimalValue(), 897.626 );
TGeodeticCoords HALC_92543( TCoords(37,53,11.26838).getDecimalValue(), TCoords(-4,2,44.35794).getDecimalValue(), 530.446 );
const TEllipsoid ellip = TEllipsoid::Ellipsoid_Internacional_1924();
// A3
TGeocentricCoords MAND_Geoc, OROZ_Geoc, SANC_Geoc, ARJO_Geoc, CANI_Geoc, HALC_Geoc;
geodeticToGeocentric( MAND_92633, MAND_Geoc, TEllipsoid::Ellipsoid_WGS84() );
geodeticToGeocentric( OROZ_96965, OROZ_Geoc, TEllipsoid::Ellipsoid_WGS84() );
geodeticToGeocentric( SANC_94744, SANC_Geoc, TEllipsoid::Ellipsoid_WGS84() );
geodeticToGeocentric( ARJO_94633, ARJO_Geoc, TEllipsoid::Ellipsoid_WGS84() );
geodeticToGeocentric( CANI_96833, CANI_Geoc, TEllipsoid::Ellipsoid_WGS84() );
geodeticToGeocentric( HALC_92543, HALC_Geoc, TEllipsoid::Ellipsoid_WGS84() );
//A5
TGeocentricCoords MAND_Geoc2, OROZ_Geoc2, SANC_Geoc2, ARJO_Geoc2, CANI_Geoc2, HALC_Geoc2;
TDatum7Params datum( 131.0320, 100.2510, 163.3540, -1.2438, -0.0195, -1.1436, -9.3900 );
transform7params( MAND_Geoc, datum, MAND_Geoc2 );
transform7params( OROZ_Geoc, datum, OROZ_Geoc2 );
transform7params( SANC_Geoc, datum, SANC_Geoc2 );
transform7params( ARJO_Geoc, datum, ARJO_Geoc2 );
transform7params( CANI_Geoc, datum, CANI_Geoc2 );
transform7params( HALC_Geoc, datum, HALC_Geoc2 );
//A4
TGeodeticCoords MAND_Geod, OROZ_Geod, SANC_Geod, ARJO_Geod, CANI_Geod, HALC_Geod;
geocentricToGeodetic( MAND_Geoc2, MAND_Geod, ellip );
geocentricToGeodetic( OROZ_Geoc2, OROZ_Geod, ellip );
geocentricToGeodetic( SANC_Geoc2, SANC_Geod, ellip );
geocentricToGeodetic( ARJO_Geoc2, ARJO_Geod, ellip );
geocentricToGeodetic( CANI_Geoc2, CANI_Geod, ellip );
geocentricToGeodetic( HALC_Geoc2, HALC_Geod, ellip );
//A2
int UTMZone;
char UTMBand;
TUTMCoords MAND_UTM, OROZ_UTM, SANC_UTM, ARJO_UTM, CANI_UTM, HALC_UTM;
geodeticToUTM( MAND_Geod, MAND_UTM, UTMZone, UTMBand, ellip );
geodeticToUTM( OROZ_Geod, OROZ_UTM, UTMZone, UTMBand, ellip );
geodeticToUTM( SANC_Geod, SANC_UTM, UTMZone, UTMBand, ellip );
geodeticToUTM( ARJO_Geod, ARJO_UTM, UTMZone, UTMBand, ellip );
geodeticToUTM( CANI_Geod, CANI_UTM, UTMZone, UTMBand, ellip );
geodeticToUTM( HALC_Geod, HALC_UTM, UTMZone, UTMBand, ellip );
vector_double thisResults;
thisResults.resize(18);
thisResults[0] = MAND_UTM.x; thisResults[1] = MAND_UTM.y; thisResults[2] = MAND_UTM.z;
thisResults[3] = OROZ_UTM.x; thisResults[4] = OROZ_UTM.y; thisResults[5] = OROZ_UTM.z;
thisResults[6] = SANC_UTM.x; thisResults[7] = SANC_UTM.y; thisResults[8] = SANC_UTM.z;
thisResults[9] = ARJO_UTM.x; thisResults[10] = ARJO_UTM.y; thisResults[11] = ARJO_UTM.z;
thisResults[12] = CANI_UTM.x; thisResults[13] = CANI_UTM.y; thisResults[14] = CANI_UTM.z;
thisResults[15] = HALC_UTM.x; thisResults[16] = HALC_UTM.y; thisResults[17] = HALC_UTM.z;
displayResults( thisResults, 2, true );
}
void Examples_03()
{
// TOPCON
TGeodeticCoords MAND_92633( TCoords(37,53,33.072573).getDecimalValue(), TCoords(-3,43,52.68965).getDecimalValue(), 515.579 );
TGeodeticCoords OROZ_96965( TCoords(37,35,31.75505).getDecimalValue(), TCoords(-3,38,5.70754).getDecimalValue(), 1489.188 );
TGeodeticCoords SANC_94744( TCoords(37,44,29.04253).getDecimalValue(), TCoords(-3,43,12.90285).getDecimalValue(), 1054.189 );
TGeodeticCoords ARJO_94633( TCoords(37,43,51.28918).getDecimalValue(), TCoords(-4,3,48.64503).getDecimalValue(), 569.182 );
TGeodeticCoords CANI_96833( TCoords(37,33,6.94237).getDecimalValue(), TCoords(-4,3,28.81611).getDecimalValue(), 897.626 );
TGeodeticCoords HALC_92543( TCoords(37,53,11.26838).getDecimalValue(), TCoords(-4,2,44.35794).getDecimalValue(), 530.446 );
const TEllipsoid ellip = TEllipsoid::Ellipsoid_WGS84();
//A2
int UTMZone;
char UTMBand;
TUTMCoords MAND_UTM, OROZ_UTM, SANC_UTM, ARJO_UTM, CANI_UTM, HALC_UTM;
geodeticToUTM( MAND_92633, MAND_UTM, UTMZone, UTMBand, ellip );
geodeticToUTM( OROZ_96965, OROZ_UTM, UTMZone, UTMBand, ellip );
geodeticToUTM( SANC_94744, SANC_UTM, UTMZone, UTMBand, ellip );
geodeticToUTM( ARJO_94633, ARJO_UTM, UTMZone, UTMBand, ellip );
geodeticToUTM( CANI_96833, CANI_UTM, UTMZone, UTMBand, ellip );
geodeticToUTM( HALC_92543, HALC_UTM, UTMZone, UTMBand, ellip );
//A7
TPoint2D MAND_point, OROZ_point, SANC_point, ARJO_point, CANI_point, HALC_point;
TDatumHelmert2D_TOPCON datum( 1.0000008551, 0.0000002900, 110.1121641064, 202.4657352579 );
transformHelmert2D_TOPCON( TPoint2D( MAND_UTM.x, MAND_UTM.y ), datum, MAND_point );
transformHelmert2D_TOPCON( TPoint2D( OROZ_UTM.x, OROZ_UTM.y ), datum, OROZ_point );
transformHelmert2D_TOPCON( TPoint2D( SANC_UTM.x, SANC_UTM.y ), datum, SANC_point );
transformHelmert2D_TOPCON( TPoint2D( ARJO_UTM.x, ARJO_UTM.y ), datum, ARJO_point );
transformHelmert2D_TOPCON( TPoint2D( CANI_UTM.x, CANI_UTM.y ), datum, CANI_point );
transformHelmert2D_TOPCON( TPoint2D( HALC_UTM.x, HALC_UTM.y ), datum, HALC_point );
vector_double thisResults;
thisResults.resize(18);
thisResults[0] = MAND_point.x; thisResults[1] = MAND_point.y; thisResults[2] = MAND_UTM.z;
thisResults[3] = OROZ_point.x; thisResults[4] = OROZ_point.y; thisResults[5] = OROZ_UTM.z;
thisResults[6] = SANC_point.x; thisResults[7] = SANC_point.y; thisResults[8] = SANC_UTM.z;
thisResults[9] = ARJO_point.x; thisResults[10] = ARJO_point.y; thisResults[11] = ARJO_UTM.z;
thisResults[12] = CANI_point.x; thisResults[13] = CANI_point.y; thisResults[14] = CANI_UTM.z;
thisResults[15] = HALC_point.x; thisResults[16] = HALC_point.y; thisResults[17] = HALC_UTM.z;
displayResults( thisResults, 3, true );
}
void Examples_04()
{
// TOPCON
TGeodeticCoords MAND_92633( TCoords(37,53,33.072573).getDecimalValue(), TCoords(-3,43,52.68965).getDecimalValue(), 515.579 );
TGeodeticCoords OROZ_96965( TCoords(37,35,31.75505).getDecimalValue(), TCoords(-3,38,5.70754).getDecimalValue(), 1489.188 );
TGeodeticCoords SANC_94744( TCoords(37,44,29.04253).getDecimalValue(), TCoords(-3,43,12.90285).getDecimalValue(), 1054.189 );
TGeodeticCoords ARJO_94633( TCoords(37,43,51.28918).getDecimalValue(), TCoords(-4,3,48.64503).getDecimalValue(), 569.182 );
TGeodeticCoords CANI_96833( TCoords(37,33,6.94237).getDecimalValue(), TCoords(-4,3,28.81611).getDecimalValue(), 897.626 );
TGeodeticCoords HALC_92543( TCoords(37,53,11.26838).getDecimalValue(), TCoords(-4,2,44.35794).getDecimalValue(), 530.446 );
const TEllipsoid ellip = TEllipsoid::Ellipsoid_WGS84();
//A2
int UTMZone;
char UTMBand;
TUTMCoords MAND_UTM, OROZ_UTM, SANC_UTM, ARJO_UTM, CANI_UTM, HALC_UTM;
geodeticToUTM( MAND_92633, MAND_UTM, UTMZone, UTMBand, ellip );
geodeticToUTM( OROZ_96965, OROZ_UTM, UTMZone, UTMBand, ellip );
geodeticToUTM( SANC_94744, SANC_UTM, UTMZone, UTMBand, ellip );
geodeticToUTM( ARJO_94633, ARJO_UTM, UTMZone, UTMBand, ellip );
geodeticToUTM( CANI_96833, CANI_UTM, UTMZone, UTMBand, ellip );
geodeticToUTM( HALC_92543, HALC_UTM, UTMZone, UTMBand, ellip );
//A7
TUTMCoords MAND_point, OROZ_point, SANC_point, ARJO_point, CANI_point, HALC_point;
TDatumHelmert3D_TOPCON datum( 0.9999969509, -0.0000070551, 142.4402790576, -0.0000060616, 1.0000054540, 185.6991268722, -49.4041666667 );
transformHelmert3D_TOPCON( MAND_UTM, datum, MAND_point );
transformHelmert3D_TOPCON( OROZ_UTM, datum, OROZ_point );
transformHelmert3D_TOPCON( SANC_UTM, datum, SANC_point );
transformHelmert3D_TOPCON( ARJO_UTM, datum, ARJO_point );
transformHelmert3D_TOPCON( CANI_UTM, datum, CANI_point );
transformHelmert3D_TOPCON( HALC_UTM, datum, HALC_point );
vector_double thisResults;
thisResults.resize(18);
thisResults[0] = MAND_point.x; thisResults[1] = MAND_point.y; thisResults[2] = MAND_point.z;
thisResults[3] = OROZ_point.x; thisResults[4] = OROZ_point.y; thisResults[5] = OROZ_point.z;
thisResults[6] = SANC_point.x; thisResults[7] = SANC_point.y; thisResults[8] = SANC_point.z;
thisResults[9] = ARJO_point.x; thisResults[10] = ARJO_point.y; thisResults[11] = ARJO_point.z;
thisResults[12] = CANI_point.x; thisResults[13] = CANI_point.y; thisResults[14] = CANI_point.z;
thisResults[15] = HALC_point.x; thisResults[16] = HALC_point.y; thisResults[17] = HALC_point.z;
displayResults( thisResults, 4, true );
}
void Examples_05()
{
// TOPCON
TGeodeticCoords MAND_92633( TCoords(37,53,33.072573).getDecimalValue(), TCoords(-3,43,52.68965).getDecimalValue(), 515.579 );
TGeodeticCoords OROZ_96965( TCoords(37,35,31.75505).getDecimalValue(), TCoords(-3,38,5.70754).getDecimalValue(), 1489.188 );
TGeodeticCoords SANC_94744( TCoords(37,44,29.04253).getDecimalValue(), TCoords(-3,43,12.90285).getDecimalValue(), 1054.189 );
TGeodeticCoords ARJO_94633( TCoords(37,43,51.28918).getDecimalValue(), TCoords(-4,3,48.64503).getDecimalValue(), 569.182 );
TGeodeticCoords CANI_96833( TCoords(37,33,6.94237).getDecimalValue(), TCoords(-4,3,28.81611).getDecimalValue(), 897.626 );
TGeodeticCoords HALC_92543( TCoords(37,53,11.26838).getDecimalValue(), TCoords(-4,2,44.35794).getDecimalValue(), 530.446 );
const TEllipsoid ellip = TEllipsoid( 6378137, 6356752.31424518, "USER" );
//A2
int UTMZone;
char UTMBand;
TUTMCoords MAND_UTM, OROZ_UTM, SANC_UTM, ARJO_UTM, CANI_UTM, HALC_UTM;
geodeticToUTM( MAND_92633, MAND_UTM, UTMZone, UTMBand, ellip );
geodeticToUTM( OROZ_96965, OROZ_UTM, UTMZone, UTMBand, ellip );
geodeticToUTM( SANC_94744, SANC_UTM, UTMZone, UTMBand, ellip );
geodeticToUTM( ARJO_94633, ARJO_UTM, UTMZone, UTMBand, ellip );
geodeticToUTM( CANI_96833, CANI_UTM, UTMZone, UTMBand, ellip );
geodeticToUTM( HALC_92543, HALC_UTM, UTMZone, UTMBand, ellip );
//A5
TUTMCoords MAND_point, OROZ_point, SANC_point, ARJO_point, CANI_point, HALC_point;
TDatum7Params_TOPCON datum(
109.035107801496, 202.610276063935, 32.727057571807,
0.9999999999, 0.0000005460, 0.0000156068,
-0.0000005463, 0.9999999998, 0.0000180851,
-0.0000156068, -0.0000180851, 0.9999999997,
0.842993545636 );
transform7params_TOPCON( MAND_UTM, datum, MAND_point );
transform7params_TOPCON( OROZ_UTM, datum, OROZ_point );
transform7params_TOPCON( SANC_UTM, datum, SANC_point );
transform7params_TOPCON( ARJO_UTM, datum, ARJO_point );
transform7params_TOPCON( CANI_UTM, datum, CANI_point );
transform7params_TOPCON( HALC_UTM, datum, HALC_point );
vector_double thisResults;
thisResults.resize(18);
thisResults[0] = MAND_point.x; thisResults[1] = MAND_point.y; thisResults[2] = MAND_point.z;
thisResults[3] = OROZ_point.x; thisResults[4] = OROZ_point.y; thisResults[5] = OROZ_point.z;
thisResults[6] = SANC_point.x; thisResults[7] = SANC_point.y; thisResults[8] = SANC_point.z;
thisResults[9] = ARJO_point.x; thisResults[10] = ARJO_point.y; thisResults[11] = ARJO_point.z;
thisResults[12] = CANI_point.x; thisResults[13] = CANI_point.y; thisResults[14] = CANI_point.z;
thisResults[15] = HALC_point.x; thisResults[16] = HALC_point.y; thisResults[17] = HALC_point.z;
displayResults( thisResults, 5, true );
}
void Examples_06()
{
TGeodeticCoords MAND_92633( TCoords(37,53,33.072573).getDecimalValue(), TCoords(-3,43,52.68965).getDecimalValue(), 515.579 );
TGeodeticCoords OROZ_96965( TCoords(37,35,31.75505).getDecimalValue(), TCoords(-3,38,5.70754).getDecimalValue(), 1489.188 );
TGeodeticCoords SANC_94744( TCoords(37,44,29.04253).getDecimalValue(), TCoords(-3,43,12.90285).getDecimalValue(), 1054.189 );
TGeodeticCoords ARJO_94633( TCoords(37,43,51.28918).getDecimalValue(), TCoords(-4,3,48.64503).getDecimalValue(), 569.182 );
TGeodeticCoords CANI_96833( TCoords(37,33,6.94237).getDecimalValue(), TCoords(-4,3,28.81611).getDecimalValue(), 897.626 );
TGeodeticCoords HALC_92543( TCoords(37,53,11.26838).getDecimalValue(), TCoords(-4,2,44.35794).getDecimalValue(), 530.446 );
const TEllipsoid ellip = TEllipsoid::Ellipsoid_Internacional_1924();
// A3
TGeocentricCoords MAND_Geoc, OROZ_Geoc, SANC_Geoc, ARJO_Geoc, CANI_Geoc, HALC_Geoc;
geodeticToGeocentric( MAND_92633, MAND_Geoc, ellip );
geodeticToGeocentric( OROZ_96965, OROZ_Geoc, ellip );
geodeticToGeocentric( SANC_94744, SANC_Geoc, ellip );
geodeticToGeocentric( ARJO_94633, ARJO_Geoc, ellip );
geodeticToGeocentric( CANI_96833, CANI_Geoc, ellip );
geodeticToGeocentric( HALC_92543, HALC_Geoc, ellip );
//A5
TGeocentricCoords MAND_Geoc2, OROZ_Geoc2, SANC_Geoc2, ARJO_Geoc2, CANI_Geoc2, HALC_Geoc2;
TDatum7Params datum( 131.0320, 100.2510, 163.3540, -1.2438, -0.0195, -1.1436, -9.3900 );
transform7params( MAND_Geoc, datum, MAND_Geoc2 );
transform7params( OROZ_Geoc, datum, OROZ_Geoc2 );
transform7params( SANC_Geoc, datum, SANC_Geoc2 );
transform7params( ARJO_Geoc, datum, ARJO_Geoc2 );
transform7params( CANI_Geoc, datum, CANI_Geoc2 );
transform7params( HALC_Geoc, datum, HALC_Geoc2 );
//A4
TGeodeticCoords MAND_Geod, OROZ_Geod, SANC_Geod, ARJO_Geod, CANI_Geod, HALC_Geod;
geocentricToGeodetic( MAND_Geoc2, MAND_Geod, ellip );
geocentricToGeodetic( OROZ_Geoc2, OROZ_Geod, ellip );
geocentricToGeodetic( SANC_Geoc2, SANC_Geod, ellip );
geocentricToGeodetic( ARJO_Geoc2, ARJO_Geod, ellip );
geocentricToGeodetic( CANI_Geoc2, CANI_Geod, ellip );
geocentricToGeodetic( HALC_Geoc2, HALC_Geod, ellip );
//A2
int UTMZone;
char UTMBand;
TUTMCoords MAND_UTM, OROZ_UTM, SANC_UTM, ARJO_UTM, CANI_UTM, HALC_UTM;
geodeticToUTM( MAND_Geod, MAND_UTM, UTMZone, UTMBand, ellip );
geodeticToUTM( OROZ_Geod, OROZ_UTM, UTMZone, UTMBand, ellip );
geodeticToUTM( SANC_Geod, SANC_UTM, UTMZone, UTMBand, ellip );
geodeticToUTM( ARJO_Geod, ARJO_UTM, UTMZone, UTMBand, ellip );
geodeticToUTM( CANI_Geod, CANI_UTM, UTMZone, UTMBand, ellip );
geodeticToUTM( HALC_Geod, HALC_UTM, UTMZone, UTMBand, ellip );
vector_double thisResults;
thisResults.resize(18);
thisResults[0] = MAND_UTM.x; thisResults[1] = MAND_UTM.y; thisResults[2] = MAND_UTM.z;
thisResults[3] = OROZ_UTM.x; thisResults[4] = OROZ_UTM.y; thisResults[5] = OROZ_UTM.z;
thisResults[6] = SANC_UTM.x; thisResults[7] = SANC_UTM.y; thisResults[8] = SANC_UTM.z;
thisResults[9] = ARJO_UTM.x; thisResults[10] = ARJO_UTM.y; thisResults[11] = ARJO_UTM.z;
thisResults[12] = CANI_UTM.x; thisResults[13] = CANI_UTM.y; thisResults[14] = CANI_UTM.z;
thisResults[15] = HALC_UTM.x; thisResults[16] = HALC_UTM.y; thisResults[17] = HALC_UTM.z;
displayResults( thisResults, 6, false);
//A7
TPoint2D MAND_point, OROZ_point, SANC_point, ARJO_point, CANI_point, HALC_point;
TDatumHelmert2D_TOPCON dat( 0.9999980621, 0.0000078756, -31.7350304827, 10.8014994871 );
transformHelmert2D_TOPCON( TPoint2D( MAND_UTM.x, MAND_UTM.y ), dat, MAND_point );
transformHelmert2D_TOPCON( TPoint2D( OROZ_UTM.x, OROZ_UTM.y ), dat, OROZ_point );
transformHelmert2D_TOPCON( TPoint2D( SANC_UTM.x, SANC_UTM.y ), dat, SANC_point );
transformHelmert2D_TOPCON( TPoint2D( ARJO_UTM.x, ARJO_UTM.y ), dat, ARJO_point );
transformHelmert2D_TOPCON( TPoint2D( CANI_UTM.x, CANI_UTM.y ), dat, CANI_point );
transformHelmert2D_TOPCON( TPoint2D( HALC_UTM.x, HALC_UTM.y ), dat, HALC_point );
thisResults[0] = MAND_point.x; thisResults[1] = MAND_point.y; thisResults[2] = MAND_UTM.z;
thisResults[3] = OROZ_point.x; thisResults[4] = OROZ_point.y; thisResults[5] = OROZ_UTM.z;
thisResults[6] = SANC_point.x; thisResults[7] = SANC_point.y; thisResults[8] = SANC_UTM.z;
thisResults[9] = ARJO_point.x; thisResults[10] = ARJO_point.y; thisResults[11] = ARJO_UTM.z;
thisResults[12] = CANI_point.x; thisResults[13] = CANI_point.y; thisResults[14] = CANI_UTM.z;
thisResults[15] = HALC_point.x; thisResults[16] = HALC_point.y; thisResults[17] = HALC_UTM.z;
displayResults( thisResults, 6, true );
}
void Examples_07()
{
TGeodeticCoords MAND_92633( TCoords(37,53,33.072573).getDecimalValue(), TCoords(-3,43,52.68965).getDecimalValue(), 515.579 );
TGeodeticCoords OROZ_96965( TCoords(37,35,31.75505).getDecimalValue(), TCoords(-3,38,5.70754).getDecimalValue(), 1489.188 );
TGeodeticCoords SANC_94744( TCoords(37,44,29.04253).getDecimalValue(), TCoords(-3,43,12.90285).getDecimalValue(), 1054.189 );
TGeodeticCoords ARJO_94633( TCoords(37,43,51.28918).getDecimalValue(), TCoords(-4,3,48.64503).getDecimalValue(), 569.182 );
TGeodeticCoords CANI_96833( TCoords(37,33,6.94237).getDecimalValue(), TCoords(-4,3,28.81611).getDecimalValue(), 897.626 );
TGeodeticCoords HALC_92543( TCoords(37,53,11.26838).getDecimalValue(), TCoords(-4,2,44.35794).getDecimalValue(), 530.446 );
const TEllipsoid ellip = TEllipsoid::Ellipsoid_Internacional_1924();
// A3
TGeocentricCoords MAND_Geoc, OROZ_Geoc, SANC_Geoc, ARJO_Geoc, CANI_Geoc, HALC_Geoc;
geodeticToGeocentric( MAND_92633, MAND_Geoc, ellip );
geodeticToGeocentric( OROZ_96965, OROZ_Geoc, ellip );
geodeticToGeocentric( SANC_94744, SANC_Geoc, ellip );
geodeticToGeocentric( ARJO_94633, ARJO_Geoc, ellip );
geodeticToGeocentric( CANI_96833, CANI_Geoc, ellip );
geodeticToGeocentric( HALC_92543, HALC_Geoc, ellip );
//A5
TGeocentricCoords MAND_Geoc2, OROZ_Geoc2, SANC_Geoc2, ARJO_Geoc2, CANI_Geoc2, HALC_Geoc2;
TDatum7Params datum( 131.0320, 100.2510, 163.3540, -1.2438, -0.0195, -1.1436, -9.3900 );
transform7params( MAND_Geoc, datum, MAND_Geoc2 );
transform7params( OROZ_Geoc, datum, OROZ_Geoc2 );
transform7params( SANC_Geoc, datum, SANC_Geoc2 );
transform7params( ARJO_Geoc, datum, ARJO_Geoc2 );
transform7params( CANI_Geoc, datum, CANI_Geoc2 );
transform7params( HALC_Geoc, datum, HALC_Geoc2 );
//A4
TGeodeticCoords MAND_Geod, OROZ_Geod, SANC_Geod, ARJO_Geod, CANI_Geod, HALC_Geod;
geocentricToGeodetic( MAND_Geoc2, MAND_Geod, ellip );
geocentricToGeodetic( OROZ_Geoc2, OROZ_Geod, ellip );
geocentricToGeodetic( SANC_Geoc2, SANC_Geod, ellip );
geocentricToGeodetic( ARJO_Geoc2, ARJO_Geod, ellip );
geocentricToGeodetic( CANI_Geoc2, CANI_Geod, ellip );
geocentricToGeodetic( HALC_Geoc2, HALC_Geod, ellip );
//A2
int UTMZone;
char UTMBand;
TUTMCoords MAND_UTM, OROZ_UTM, SANC_UTM, ARJO_UTM, CANI_UTM, HALC_UTM;
geodeticToUTM( MAND_Geod, MAND_UTM, UTMZone, UTMBand, ellip );
geodeticToUTM( OROZ_Geod, OROZ_UTM, UTMZone, UTMBand, ellip );
geodeticToUTM( SANC_Geod, SANC_UTM, UTMZone, UTMBand, ellip );
geodeticToUTM( ARJO_Geod, ARJO_UTM, UTMZone, UTMBand, ellip );
geodeticToUTM( CANI_Geod, CANI_UTM, UTMZone, UTMBand, ellip );
geodeticToUTM( HALC_Geod, HALC_UTM, UTMZone, UTMBand, ellip );
vector_double thisResults;
thisResults.resize(18);
thisResults[0] = MAND_UTM.x; thisResults[1] = MAND_UTM.y; thisResults[2] = MAND_UTM.z;
thisResults[3] = OROZ_UTM.x; thisResults[4] = OROZ_UTM.y; thisResults[5] = OROZ_UTM.z;
thisResults[6] = SANC_UTM.x; thisResults[7] = SANC_UTM.y; thisResults[8] = SANC_UTM.z;
thisResults[9] = ARJO_UTM.x; thisResults[10] = ARJO_UTM.y; thisResults[11] = ARJO_UTM.z;
thisResults[12] = CANI_UTM.x; thisResults[13] = CANI_UTM.y; thisResults[14] = CANI_UTM.z;
thisResults[15] = HALC_UTM.x; thisResults[16] = HALC_UTM.y; thisResults[17] = HALC_UTM.z;
displayResults( thisResults, 7, false);
//A8
TPoint3D MAND_point, OROZ_point, SANC_point, ARJO_point, CANI_point, HALC_point;
TDatumHelmert3D_TOPCON dat( 0.9999944446, 0.0000005236, 0.5022460081, -0.0000136628, 1.0000023039, -4.4675776973, 29.2747378118 );
transformHelmert3D_TOPCON( MAND_UTM, dat, MAND_point );
transformHelmert3D_TOPCON( OROZ_UTM, dat, OROZ_point );
transformHelmert3D_TOPCON( SANC_UTM, dat, SANC_point );
transformHelmert3D_TOPCON( ARJO_UTM, dat, ARJO_point );
transformHelmert3D_TOPCON( CANI_UTM, dat, CANI_point );
transformHelmert3D_TOPCON( HALC_UTM, dat, HALC_point );
thisResults[0] = MAND_point.x; thisResults[1] = MAND_point.y; thisResults[2] = MAND_point.z;
thisResults[3] = OROZ_point.x; thisResults[4] = OROZ_point.y; thisResults[5] = OROZ_point.z;
thisResults[6] = SANC_point.x; thisResults[7] = SANC_point.y; thisResults[8] = SANC_point.z;
thisResults[9] = ARJO_point.x; thisResults[10] = ARJO_point.y; thisResults[11] = ARJO_point.z;
thisResults[12] = CANI_point.x; thisResults[13] = CANI_point.y; thisResults[14] = CANI_point.z;
thisResults[15] = HALC_point.x; thisResults[16] = HALC_point.y; thisResults[17] = HALC_point.z;
displayResults( thisResults, 7, true );
}
void Examples_08()
{
TGeodeticCoords MAND_92633( TCoords(37,53,33.072573).getDecimalValue(), TCoords(-3,43,52.68965).getDecimalValue(), 515.579 );
TGeodeticCoords OROZ_96965( TCoords(37,35,31.75505).getDecimalValue(), TCoords(-3,38,5.70754).getDecimalValue(), 1489.188 );
TGeodeticCoords SANC_94744( TCoords(37,44,29.04253).getDecimalValue(), TCoords(-3,43,12.90285).getDecimalValue(), 1054.189 );
TGeodeticCoords ARJO_94633( TCoords(37,43,51.28918).getDecimalValue(), TCoords(-4,3,48.64503).getDecimalValue(), 569.182 );
TGeodeticCoords CANI_96833( TCoords(37,33,6.94237).getDecimalValue(), TCoords(-4,3,28.81611).getDecimalValue(), 897.626 );
TGeodeticCoords HALC_92543( TCoords(37,53,11.26838).getDecimalValue(), TCoords(-4,2,44.35794).getDecimalValue(), 530.446 );
const TEllipsoid ellip = TEllipsoid::Ellipsoid_Internacional_1924();
// A3
TGeocentricCoords MAND_Geoc, OROZ_Geoc, SANC_Geoc, ARJO_Geoc, CANI_Geoc, HALC_Geoc;
geodeticToGeocentric( MAND_92633, MAND_Geoc, ellip );
geodeticToGeocentric( OROZ_96965, OROZ_Geoc, ellip );
geodeticToGeocentric( SANC_94744, SANC_Geoc, ellip );
geodeticToGeocentric( ARJO_94633, ARJO_Geoc, ellip );
geodeticToGeocentric( CANI_96833, CANI_Geoc, ellip );
geodeticToGeocentric( HALC_92543, HALC_Geoc, ellip );
//A5
TGeocentricCoords MAND_Geoc2, OROZ_Geoc2, SANC_Geoc2, ARJO_Geoc2, CANI_Geoc2, HALC_Geoc2;
TDatum7Params datum( 131.0320, 100.2510, 163.3540, -1.2438, -0.0195, -1.1436, -9.3900 );
transform7params( MAND_Geoc, datum, MAND_Geoc2 );
transform7params( OROZ_Geoc, datum, OROZ_Geoc2 );
transform7params( SANC_Geoc, datum, SANC_Geoc2 );
transform7params( ARJO_Geoc, datum, ARJO_Geoc2 );
transform7params( CANI_Geoc, datum, CANI_Geoc2 );
transform7params( HALC_Geoc, datum, HALC_Geoc2 );
//A4
TGeodeticCoords MAND_Geod, OROZ_Geod, SANC_Geod, ARJO_Geod, CANI_Geod, HALC_Geod;
geocentricToGeodetic( MAND_Geoc2, MAND_Geod, ellip );
geocentricToGeodetic( OROZ_Geoc2, OROZ_Geod, ellip );
geocentricToGeodetic( SANC_Geoc2, SANC_Geod, ellip );
geocentricToGeodetic( ARJO_Geoc2, ARJO_Geod, ellip );
geocentricToGeodetic( CANI_Geoc2, CANI_Geod, ellip );
geocentricToGeodetic( HALC_Geoc2, HALC_Geod, ellip );
//A2
int UTMZone;
char UTMBand;
TUTMCoords MAND_UTM, OROZ_UTM, SANC_UTM, ARJO_UTM, CANI_UTM, HALC_UTM;
geodeticToUTM( MAND_Geod, MAND_UTM, UTMZone, UTMBand, ellip );
geodeticToUTM( OROZ_Geod, OROZ_UTM, UTMZone, UTMBand, ellip );
geodeticToUTM( SANC_Geod, SANC_UTM, UTMZone, UTMBand, ellip );
geodeticToUTM( ARJO_Geod, ARJO_UTM, UTMZone, UTMBand, ellip );
geodeticToUTM( CANI_Geod, CANI_UTM, UTMZone, UTMBand, ellip );
geodeticToUTM( HALC_Geod, HALC_UTM, UTMZone, UTMBand, ellip );
vector_double thisResults;
thisResults.resize(18);
thisResults[0] = MAND_UTM.x; thisResults[1] = MAND_UTM.y; thisResults[2] = MAND_UTM.z;
thisResults[3] = OROZ_UTM.x; thisResults[4] = OROZ_UTM.y; thisResults[5] = OROZ_UTM.z;
thisResults[6] = SANC_UTM.x; thisResults[7] = SANC_UTM.y; thisResults[8] = SANC_UTM.z;
thisResults[9] = ARJO_UTM.x; thisResults[10] = ARJO_UTM.y; thisResults[11] = ARJO_UTM.z;
thisResults[12] = CANI_UTM.x; thisResults[13] = CANI_UTM.y; thisResults[14] = CANI_UTM.z;
thisResults[15] = HALC_UTM.x; thisResults[16] = HALC_UTM.y; thisResults[17] = HALC_UTM.z;
displayResults( thisResults, 8, false);
//A5
TGeocentricCoords MAND_Geoc3, OROZ_Geoc3, SANC_Geoc3, ARJO_Geoc3, CANI_Geoc3, HALC_Geoc3;
//TDatum7Params dat( -34.095438458048, 10.843234125336, 212.169692142148, -8.348308019638, 6.751893142065, -1.738975581025, -1.898114919019 );
TDatum7Params_TOPCON dat( -34.095438458048, 10.843234125336, 212.169692142148,
0.9999999994, 0.0000084308, 0.0000327341,
-0.0000084321, 0.9999999991, 0.0000404737,
-0.0000327338, -0.0000404737, 0.9999999986,
-1.898114919019 );
transform7params_TOPCON( MAND_UTM, dat, MAND_Geoc3 );
transform7params_TOPCON( OROZ_UTM, dat, OROZ_Geoc3 );
transform7params_TOPCON( SANC_UTM, dat, SANC_Geoc3 );
transform7params_TOPCON( ARJO_UTM, dat, ARJO_Geoc3 );
transform7params_TOPCON( CANI_UTM, dat, CANI_Geoc3 );
transform7params_TOPCON( HALC_UTM, dat, HALC_Geoc3 );
thisResults[0] = MAND_Geoc3.x; thisResults[1] = MAND_Geoc3.y; thisResults[2] = MAND_Geoc3.z;
thisResults[3] = OROZ_Geoc3.x; thisResults[4] = OROZ_Geoc3.y; thisResults[5] = OROZ_Geoc3.z;
thisResults[6] = SANC_Geoc3.x; thisResults[7] = SANC_Geoc3.y; thisResults[8] = SANC_Geoc3.z;
thisResults[9] = ARJO_Geoc3.x; thisResults[10] = ARJO_Geoc3.y; thisResults[11] = ARJO_Geoc3.z;
thisResults[12] = CANI_Geoc3.x; thisResults[13] = CANI_Geoc3.y; thisResults[14] = CANI_Geoc3.z;
thisResults[15] = HALC_Geoc3.x; thisResults[16] = HALC_Geoc3.y; thisResults[17] = HALC_Geoc3.z;
displayResults( thisResults, 8, true );
}
void Examples_10()
{
TGeodeticCoords MAND_92633( TCoords(37,53,33.072573).getDecimalValue(), TCoords(-3,43,52.68965).getDecimalValue(), 515.579 );
TGeodeticCoords OROZ_96965( TCoords(37,35,31.75505).getDecimalValue(), TCoords(-3,38,5.70754).getDecimalValue(), 1489.188 );
TGeodeticCoords SANC_94744( TCoords(37,44,29.04253).getDecimalValue(), TCoords(-3,43,12.90285).getDecimalValue(), 1054.189 );
TGeodeticCoords ARJO_94633( TCoords(37,43,51.28918).getDecimalValue(), TCoords(-4,3,48.64503).getDecimalValue(), 569.182 );
TGeodeticCoords CANI_96833( TCoords(37,33,6.94237).getDecimalValue(), TCoords(-4,3,28.81611).getDecimalValue(), 897.626 );
TGeodeticCoords HALC_92543( TCoords(37,53,11.26838).getDecimalValue(), TCoords(-4,2,44.35794).getDecimalValue(), 530.446 );
const TEllipsoid ellip = TEllipsoid::Ellipsoid_Hayford_1909();
// A3
TGeocentricCoords MAND_Geoc, OROZ_Geoc, SANC_Geoc, ARJO_Geoc, CANI_Geoc, HALC_Geoc;
geodeticToGeocentric( MAND_92633, MAND_Geoc, ellip );
geodeticToGeocentric( OROZ_96965, OROZ_Geoc, ellip );
geodeticToGeocentric( SANC_94744, SANC_Geoc, ellip );
geodeticToGeocentric( ARJO_94633, ARJO_Geoc, ellip );
geodeticToGeocentric( CANI_96833, CANI_Geoc, ellip );
geodeticToGeocentric( HALC_92543, HALC_Geoc, ellip );
//A5
TGeocentricCoords MAND_Geoc2, OROZ_Geoc2, SANC_Geoc2, ARJO_Geoc2, CANI_Geoc2, HALC_Geoc2;
TDatum7Params datum( 284.2535, -116.9549, -34.8027, 4.78337, 7.96684, -5.37707, -6.6900 );
transform7params( MAND_Geoc, datum, MAND_Geoc2 );
transform7params( OROZ_Geoc, datum, OROZ_Geoc2 );
transform7params( SANC_Geoc, datum, SANC_Geoc2 );
transform7params( ARJO_Geoc, datum, ARJO_Geoc2 );
transform7params( CANI_Geoc, datum, CANI_Geoc2 );
transform7params( HALC_Geoc, datum, HALC_Geoc2 );
//A4
TGeodeticCoords MAND_Geod, OROZ_Geod, SANC_Geod, ARJO_Geod, CANI_Geod, HALC_Geod;
geocentricToGeodetic( MAND_Geoc2, MAND_Geod, ellip );
geocentricToGeodetic( OROZ_Geoc2, OROZ_Geod, ellip );
geocentricToGeodetic( SANC_Geoc2, SANC_Geod, ellip );
geocentricToGeodetic( ARJO_Geoc2, ARJO_Geod, ellip );
geocentricToGeodetic( CANI_Geoc2, CANI_Geod, ellip );
geocentricToGeodetic( HALC_Geoc2, HALC_Geod, ellip );
//A2
int UTMZone;
char UTMBand;
TUTMCoords MAND_UTM, OROZ_UTM, SANC_UTM, ARJO_UTM, CANI_UTM, HALC_UTM;
geodeticToUTM( MAND_Geod, MAND_UTM, UTMZone, UTMBand, ellip );
geodeticToUTM( OROZ_Geod, OROZ_UTM, UTMZone, UTMBand, ellip );
geodeticToUTM( SANC_Geod, SANC_UTM, UTMZone, UTMBand, ellip );
geodeticToUTM( ARJO_Geod, ARJO_UTM, UTMZone, UTMBand, ellip );
geodeticToUTM( CANI_Geod, CANI_UTM, UTMZone, UTMBand, ellip );
geodeticToUTM( HALC_Geod, HALC_UTM, UTMZone, UTMBand, ellip );
vector_double thisResults;
thisResults.resize(18);
thisResults[0] = MAND_UTM.x; thisResults[1] = MAND_UTM.y; thisResults[2] = MAND_UTM.z;
thisResults[3] = OROZ_UTM.x; thisResults[4] = OROZ_UTM.y; thisResults[5] = OROZ_UTM.z;
thisResults[6] = SANC_UTM.x; thisResults[7] = SANC_UTM.y; thisResults[8] = SANC_UTM.z;
thisResults[9] = ARJO_UTM.x; thisResults[10] = ARJO_UTM.y; thisResults[11] = ARJO_UTM.z;
thisResults[12] = CANI_UTM.x; thisResults[13] = CANI_UTM.y; thisResults[14] = CANI_UTM.z;
thisResults[15] = HALC_UTM.x; thisResults[16] = HALC_UTM.y; thisResults[17] = HALC_UTM.z;
displayResults( thisResults, 10, true);
}
void Examples_11()
{
TGeodeticCoords MAND_92633( TCoords(37,53,33.072573).getDecimalValue(), TCoords(-3,43,52.68965).getDecimalValue(), 515.579 );
TGeodeticCoords OROZ_96965( TCoords(37,35,31.75505).getDecimalValue(), TCoords(-3,38,5.70754).getDecimalValue(), 1489.188 );
TGeodeticCoords SANC_94744( TCoords(37,44,29.04253).getDecimalValue(), TCoords(-3,43,12.90285).getDecimalValue(), 1054.189 );
TGeodeticCoords ARJO_94633( TCoords(37,43,51.28918).getDecimalValue(), TCoords(-4,3,48.64503).getDecimalValue(), 569.182 );
TGeodeticCoords CANI_96833( TCoords(37,33,6.94237).getDecimalValue(), TCoords(-4,3,28.81611).getDecimalValue(), 897.626 );
TGeodeticCoords HALC_92543( TCoords(37,53,11.26838).getDecimalValue(), TCoords(-4,2,44.35794).getDecimalValue(), 530.446 );
const TEllipsoid ellip = TEllipsoid::Ellipsoid_Hayford_1909();
// A3
TGeocentricCoords MAND_Geoc, OROZ_Geoc, SANC_Geoc, ARJO_Geoc, CANI_Geoc, HALC_Geoc;
geodeticToGeocentric( MAND_92633, MAND_Geoc, ellip );
geodeticToGeocentric( OROZ_96965, OROZ_Geoc, ellip );
geodeticToGeocentric( SANC_94744, SANC_Geoc, ellip );
geodeticToGeocentric( ARJO_94633, ARJO_Geoc, ellip );
geodeticToGeocentric( CANI_96833, CANI_Geoc, ellip );
geodeticToGeocentric( HALC_92543, HALC_Geoc, ellip );
//A6
TGeocentricCoords MAND_Geoc2, OROZ_Geoc2, SANC_Geoc2, ARJO_Geoc2, CANI_Geoc2, HALC_Geoc2;
TDatum10Params datum( 109.4714, 106.7407, 141.7916, 5039726.4242, -341417.6071, 3882515.5524, 4.78337, 7.96684, -5.37707, -6.6900 );
transform10params( MAND_Geoc, datum, MAND_Geoc2 );
transform10params( OROZ_Geoc, datum, OROZ_Geoc2 );
transform10params( SANC_Geoc, datum, SANC_Geoc2 );
transform10params( ARJO_Geoc, datum, ARJO_Geoc2 );
transform10params( CANI_Geoc, datum, CANI_Geoc2 );
transform10params( HALC_Geoc, datum, HALC_Geoc2 );
//A4
TGeodeticCoords MAND_Geod, OROZ_Geod, SANC_Geod, ARJO_Geod, CANI_Geod, HALC_Geod;
geocentricToGeodetic( MAND_Geoc2, MAND_Geod, ellip );
geocentricToGeodetic( OROZ_Geoc2, OROZ_Geod, ellip );
geocentricToGeodetic( SANC_Geoc2, SANC_Geod, ellip );
geocentricToGeodetic( ARJO_Geoc2, ARJO_Geod, ellip );
geocentricToGeodetic( CANI_Geoc2, CANI_Geod, ellip );
geocentricToGeodetic( HALC_Geoc2, HALC_Geod, ellip );
//A2
int UTMZone;
char UTMBand;
TUTMCoords MAND_UTM, OROZ_UTM, SANC_UTM, ARJO_UTM, CANI_UTM, HALC_UTM;
geodeticToUTM( MAND_Geod, MAND_UTM, UTMZone, UTMBand, ellip );
geodeticToUTM( OROZ_Geod, OROZ_UTM, UTMZone, UTMBand, ellip );
geodeticToUTM( SANC_Geod, SANC_UTM, UTMZone, UTMBand, ellip );
geodeticToUTM( ARJO_Geod, ARJO_UTM, UTMZone, UTMBand, ellip );
geodeticToUTM( CANI_Geod, CANI_UTM, UTMZone, UTMBand, ellip );
geodeticToUTM( HALC_Geod, HALC_UTM, UTMZone, UTMBand, ellip );
vector_double thisResults;
thisResults.resize(18);
thisResults[0] = MAND_UTM.x; thisResults[1] = MAND_UTM.y; thisResults[2] = MAND_UTM.z;
thisResults[3] = OROZ_UTM.x; thisResults[4] = OROZ_UTM.y; thisResults[5] = OROZ_UTM.z;
thisResults[6] = SANC_UTM.x; thisResults[7] = SANC_UTM.y; thisResults[8] = SANC_UTM.z;
thisResults[9] = ARJO_UTM.x; thisResults[10] = ARJO_UTM.y; thisResults[11] = ARJO_UTM.z;
thisResults[12] = CANI_UTM.x; thisResults[13] = CANI_UTM.y; thisResults[14] = CANI_UTM.z;
thisResults[15] = HALC_UTM.x; thisResults[16] = HALC_UTM.y; thisResults[17] = HALC_UTM.z;
displayResults( thisResults, 11, true);
}
void Examples_12()
{
TGeodeticCoords MAND_92633( TCoords(37,53,33.072573).getDecimalValue(), TCoords(-3,43,52.68965).getDecimalValue(), 515.579 );
TGeodeticCoords OROZ_96965( TCoords(37,35,31.75505).getDecimalValue(), TCoords(-3,38,5.70754).getDecimalValue(), 1489.188 );
TGeodeticCoords SANC_94744( TCoords(37,44,29.04253).getDecimalValue(), TCoords(-3,43,12.90285).getDecimalValue(), 1054.189 );
TGeodeticCoords ARJO_94633( TCoords(37,43,51.28918).getDecimalValue(), TCoords(-4,3,48.64503).getDecimalValue(), 569.182 );
TGeodeticCoords CANI_96833( TCoords(37,33,6.94237).getDecimalValue(), TCoords(-4,3,28.81611).getDecimalValue(), 897.626 );
TGeodeticCoords HALC_92543( TCoords(37,53,11.26838).getDecimalValue(), TCoords(-4,2,44.35794).getDecimalValue(), 530.446 );
const TEllipsoid ellip = TEllipsoid::Ellipsoid_Hayford_1909();
//A2
int UTMZone;
char UTMBand;
TUTMCoords MAND_UTM, OROZ_UTM, SANC_UTM, ARJO_UTM, CANI_UTM, HALC_UTM;
geodeticToUTM( MAND_92633, MAND_UTM, UTMZone, UTMBand, ellip );
geodeticToUTM( OROZ_96965, OROZ_UTM, UTMZone, UTMBand, ellip );
geodeticToUTM( SANC_94744, SANC_UTM, UTMZone, UTMBand, ellip );
geodeticToUTM( ARJO_94633, ARJO_UTM, UTMZone, UTMBand, ellip );
geodeticToUTM( CANI_96833, CANI_UTM, UTMZone, UTMBand, ellip );
geodeticToUTM( HALC_92543, HALC_UTM, UTMZone, UTMBand, ellip );
//A7
TPoint2D MAND_point, OROZ_point, SANC_point, ARJO_point, CANI_point, HALC_point;
TDatumHelmert2D dat( 15.8778, 0.4491, TCoords(0,32,8.75881).getDecimalValue() /*degrees*/, -325.9604, 422956.6976, 4176709.7722 );
transformHelmert2D( TPoint2D( MAND_UTM.x, MAND_UTM.y ), dat, MAND_point );
transformHelmert2D( TPoint2D( OROZ_UTM.x, OROZ_UTM.y ), dat, OROZ_point );
transformHelmert2D( TPoint2D( SANC_UTM.x, SANC_UTM.y ), dat, SANC_point );
transformHelmert2D( TPoint2D( ARJO_UTM.x, ARJO_UTM.y ), dat, ARJO_point );
transformHelmert2D( TPoint2D( CANI_UTM.x, CANI_UTM.y ), dat, CANI_point );
transformHelmert2D( TPoint2D( HALC_UTM.x, HALC_UTM.y ), dat, HALC_point );
//A9
TPoint3D MAND_point2, OROZ_point2, SANC_point2, ARJO_point2, CANI_point2, HALC_point2;
TDatum1DTransf datum( -0.00001847 /*rad*/, -0.00001536/*rad*/, -49.4041, 0 );
transform1D( TPoint3D( MAND_point.x, MAND_point.y, MAND_UTM.z ), datum, MAND_point2 );
transform1D( TPoint3D( OROZ_point.x, OROZ_point.y, OROZ_UTM.z ), datum, OROZ_point2 );
transform1D( TPoint3D( SANC_point.x, SANC_point.y, SANC_UTM.z ), datum, SANC_point2 );
transform1D( TPoint3D( ARJO_point.x, ARJO_point.y, ARJO_UTM.z ), datum, ARJO_point2 );
transform1D( TPoint3D( CANI_point.x, CANI_point.y, CANI_UTM.z ), datum, CANI_point2 );
transform1D( TPoint3D( HALC_point.x, HALC_point.y, HALC_UTM.z ), datum, HALC_point2 );
vector_double thisResults;
thisResults.resize(18);
thisResults[0] = MAND_point2.x; thisResults[1] = MAND_point2.y; thisResults[2] = MAND_point2.z;
thisResults[3] = OROZ_point2.x; thisResults[4] = OROZ_point2.y; thisResults[5] = OROZ_point2.z;
thisResults[6] = SANC_point2.x; thisResults[7] = SANC_point2.y; thisResults[8] = SANC_point2.z;
thisResults[9] = ARJO_point2.x; thisResults[10] = ARJO_point2.y; thisResults[11] = ARJO_point2.z;
thisResults[12] = CANI_point2.x; thisResults[13] = CANI_point2.y; thisResults[14] = CANI_point2.z;
thisResults[15] = HALC_point2.x; thisResults[16] = HALC_point2.y; thisResults[17] = HALC_point2.z;
displayResults( thisResults, 12, true);
}
void Examples_13()
{
TGeodeticCoords MAND_92633( TCoords(37,53,33.072573).getDecimalValue(), TCoords(-3,43,52.68965).getDecimalValue(), 515.579 );
TGeodeticCoords OROZ_96965( TCoords(37,35,31.75505).getDecimalValue(), TCoords(-3,38,5.70754).getDecimalValue(), 1489.188 );
TGeodeticCoords SANC_94744( TCoords(37,44,29.04253).getDecimalValue(), TCoords(-3,43,12.90285).getDecimalValue(), 1054.189 );
TGeodeticCoords ARJO_94633( TCoords(37,43,51.28918).getDecimalValue(), TCoords(-4,3,48.64503).getDecimalValue(), 569.182 );
TGeodeticCoords CANI_96833( TCoords(37,33,6.94237).getDecimalValue(), TCoords(-4,3,28.81611).getDecimalValue(), 897.626 );
TGeodeticCoords HALC_92543( TCoords(37,53,11.26838).getDecimalValue(), TCoords(-4,2,44.35794).getDecimalValue(), 530.446 );
const TEllipsoid ellip = TEllipsoid::Ellipsoid_Hayford_1909();
//A2
int UTMZone;
char UTMBand;
TUTMCoords MAND_UTM, OROZ_UTM, SANC_UTM, ARJO_UTM, CANI_UTM, HALC_UTM;
geodeticToUTM( MAND_92633, MAND_UTM, UTMZone, UTMBand, ellip );
geodeticToUTM( OROZ_96965, OROZ_UTM, UTMZone, UTMBand, ellip );
geodeticToUTM( SANC_94744, SANC_UTM, UTMZone, UTMBand, ellip );
geodeticToUTM( ARJO_94633, ARJO_UTM, UTMZone, UTMBand, ellip );
geodeticToUTM( CANI_96833, CANI_UTM, UTMZone, UTMBand, ellip );
geodeticToUTM( HALC_92543, HALC_UTM, UTMZone, UTMBand, ellip );
//A10
TPoint3D MAND_point, OROZ_point, SANC_point, ARJO_point, CANI_point, HALC_point;
TDatumTransfInterpolation dat( 423150.0514, 4197393.9437, -334.2393 /*ppm*/, -314.2785 /*ppm*/, 1.33958 /*secs*/ );
transfInterpolation( MAND_UTM, dat, MAND_point );
transfInterpolation( OROZ_UTM, dat, OROZ_point );
transfInterpolation( SANC_UTM, dat, SANC_point );
transfInterpolation( ARJO_UTM, dat, ARJO_point );
transfInterpolation( CANI_UTM, dat, CANI_point );
transfInterpolation( HALC_UTM, dat, HALC_point );
//A9
TPoint3D MAND_point2, OROZ_point2, SANC_point2, ARJO_point2, CANI_point2, HALC_point2;
TDatum1DTransf datum( -0.00001822 /*rad*/, -0.00001304/*rad*/, -49.8094, 0 );
transform1D( MAND_point, datum, MAND_point2 );
transform1D( OROZ_point, datum, OROZ_point2 );
transform1D( SANC_point, datum, SANC_point2 );
transform1D( ARJO_point, datum, ARJO_point2 );
transform1D( CANI_point, datum, CANI_point2 );
transform1D( HALC_point, datum, HALC_point2 );
vector_double thisResults;
thisResults.resize(18);
thisResults[0] = MAND_point2.x; thisResults[1] = MAND_point2.y; thisResults[2] = MAND_point2.z;
thisResults[3] = OROZ_point2.x; thisResults[4] = OROZ_point2.y; thisResults[5] = OROZ_point2.z;
thisResults[6] = SANC_point2.x; thisResults[7] = SANC_point2.y; thisResults[8] = SANC_point2.z;
thisResults[9] = ARJO_point2.x; thisResults[10] = ARJO_point2.y; thisResults[11] = ARJO_point2.z;
thisResults[12] = CANI_point2.x; thisResults[13] = CANI_point2.y; thisResults[14] = CANI_point2.z;
thisResults[15] = HALC_point2.x; thisResults[16] = HALC_point2.y; thisResults[17] = HALC_point2.z;
displayResults( thisResults, 13, true);
}
int main(int argc, char **argv)
{
try
{
exampleResults();
TestCoordinatesConversions();
mrpt::system::setConsoleColor(CONCOL_BLUE);
cout << "**************************************" << endl;
cout << "\tTOPCON with th = " << TH << endl;
cout << "**************************************" << endl;
mrpt::system::setConsoleColor(CONCOL_NORMAL);
Examples_01();
Examples_02();
Examples_03();
Examples_04();
Examples_05();
Examples_06();
Examples_07();
Examples_08();
mrpt::system::setConsoleColor(CONCOL_BLUE);
cout << "**************************************" << endl;
cout << "\tLEICA with th = " << TH << endl;
cout << "**************************************" << endl;
mrpt::system::setConsoleColor(CONCOL_NORMAL);
cout << "Example 09 missing in pdf" << endl;
Examples_10();
Examples_11();
Examples_12();
Examples_13();
return 0;
} catch (std::exception &e)
{
std::cout << "MRPT exception caught: " << e.what() << std::endl;
return -1;
}
catch (...)
{
printf("Another exception!!");
return -1;
}
}
