Changeset 9ae30cc in tol
- Timestamp:
- Jan 2, 2010, 4:34:37 PM (17 years ago)
- Branches:
- 1941-docker-image-to-test, devel, master
- Children:
- 42020ce
- Parents:
- f9b4992
- Location:
- tol_pkg/GisTools
- Files:
-
- 2 edited
-
GisTools.tol (modified) (3 diffs)
-
test/test_0002/test.tol (modified) (1 diff)
Legend:
- Unmodified
- Added
- Removed
-
tol_pkg/GisTools/GisTools.tol
rf9b4992 r9ae30cc 124 124 125 125 //////////////////////////////////////////////////////////////////////////// 126 Static Set datum = Copy(Empty); 127 Static Real BuildDatum(Real void) 128 //////////////////////////////////////////////////////////////////////////// 129 { 130 Set datum := [[ 126 Static Set _datum = Copy(Empty); 127 Static Set Datum(Real void) 128 //////////////////////////////////////////////////////////////////////////// 129 { 130 If(Card(_datum),_datum, { 131 Set _datum := [[ 131 132 @UTMFormula::New("NAD83", 6378137, 6356752.3142, "Global"); 132 133 @UTMFormula::New("WGS84", 6378137, 6356752.3142, "Global"); … … 142 143 @UTMFormula::New("Bessel.1841", 6377397.2, 6356079.0, "Central Europe, Chile, Indonesia"); 143 144 @UTMFormula::New("Everest.1830", 6377276.3, 6356075.4, "South Asia") ]]; 144 Real SetIndexByName( datum);145 True145 Real SetIndexByName(_datum); 146 _datum}) 146 147 }; 147 148 … … 340 341 }; 341 342 342 343 Real _UTMFormula.BuildDatum = @UTMFormula::BuildDatum(?);344 345 343 ////////////////////////////////////////////////////////////////////////////// 346 344 //Abstract class to handle with a bidimensional point. -
tol_pkg/GisTools/test/test_0002/test.tol
rf9b4992 r9ae30cc 21 21 //#Embed "../../../_embed_all.tol"; 22 22 23 Set GisTools::@UTMFormula::datum;23 Set datum = GisTools::@UTMFormula::Datum(0); 24 24 /* */ 25 25 26 Matrix data = Include("sample.bmt")[1];27 VMatrix lon = Mat2VMat(SubCol(data,[[1]]));28 VMatrix lat = Mat2VMat(SubCol(data,[[2]]));29 VMatrix deg = lon|lat;30 VMatrix lonZone = Mat2VMat(SubCol(data,[[3]]));31 VMatrix latZone = Floor(lat/8)*8;32 VMatrix x = Mat2VMat(SubCol(data,[[4]]));33 VMatrix y = Mat2VMat(SubCol(data,[[5]]));34 VMatrix utm = x|y;26 Matrix data = Include("sample.bmt")[1]; 27 VMatrix lon = Mat2VMat(SubCol(data,[[1]])); 28 VMatrix lat = Mat2VMat(SubCol(data,[[2]])); 29 VMatrix deg = lon|lat; 30 VMatrix lonZone = Mat2VMat(SubCol(data,[[3]])); 31 VMatrix latZone = Floor(lat/8)*8; 32 VMatrix x = Mat2VMat(SubCol(data,[[4]])); 33 VMatrix y = Mat2VMat(SubCol(data,[[5]])); 34 VMatrix utm = x|y; 35 35 36 /* * / 37 VMatrix deg = Mat2VMat(Row(80.05111111,57.38666667)); 38 VMatrix utm = Mat2VMat(Row(442957.26,6360827.14)); 39 /* * / 40 VMatrix deg = { 41 Matrix aux= ( 42 (80.05111111,-57.38666667), 43 (80.05111111,10.38666667), 44 (80.05111111,-10.38666667), 45 (80.05111111,57.38666667)); 46 Mat2VMat(aux)}; 47 VMatrix utm = { 48 Matrix aux= ( 49 (442957.26,3639172.86), 50 (396126.87,1148319.33), 51 (396126.87,8851680.67), 52 (442957.26,6360827.14)); 53 Mat2VMat(aux)}; 54 /* */ 55 Real zone = 44; 56 VMatrix utm_ = GisTools::@UTMFormula::datum::WGS84::deg2utm(deg,lonZone); 57 VMatrix deg_ = GisTools::@UTMFormula::datum::WGS84::utm2deg(utm,lonZone,latZone); 36 VMatrix utm_ = datum::WGS84::deg2utm(deg,lonZone); 37 VMatrix deg_ = datum::WGS84::utm2deg(utm,lonZone,latZone); 58 38 59 39 VMatrix cmp = deg | deg_ | utm | utm_ | lonZone | latZone; 60 40 Real deg.quality = Max(0,1-VMatMax(Abs(deg-deg_))); 61 Real utm.quality = Max(0,1-VMatMax(Round(Abs(utm-utm_)))) 41 Real utm.quality = Max(0,1-VMatMax(Round(Abs(utm-utm_)))); 62 42 63 /* 43 64 44 Real numErr1 = Copy(NError); 65 45 Real numWar1 = Copy(NWarning); 66 46 67 Set partialResults_ = [[numErr0, numErr1 ]];47 Set partialResults_ = [[numErr0, numErr1,deg.quality,utm.quality]]; 68 48 69 49 //This is a messure of the success of the test 70 Real quality_ = And(numErr1 == numErr0) ;50 Real quality_ = And(numErr1 == numErr0)*deg.quality*utm.quality; 71 51 72 52 //Return the results
Note: See TracChangeset
for help on using the changeset viewer.
