Changeset 9ae30cc in tol for tol_pkg/GisTools/test/test_0002/test.tol
- Timestamp:
- Jan 2, 2010, 4:34:37 PM (17 years ago)
- Branches:
- 1941-docker-image-to-test, devel, master
- Children:
- 42020ce
- Parents:
- f9b4992
- File:
-
- 1 edited
-
tol_pkg/GisTools/test/test_0002/test.tol (modified) (1 diff)
Legend:
- Unmodified
- Added
- Removed
-
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.
