UTMtest.frink

View or download UTMtest.frink in plain text format


// Torture test for UTM to lat-long conversions.

use CoordinateConversions.frink

biggest = 0
errorsum = 0
samples = 1000

for n = 1 to samples
{
   lat = randomFloat[-80,72] degrees
   long = randomFloat[-180,180] degrees

   [x,y] = LatLongToUTM[lat,long,Datum.WGS84]
   [lat1,long1] = UTMToLatLong[x,y,LatLongToUTMZone[lat,long],Datum.WGS84]
   error = sqrt[(lat1-lat)^2 + (long1-long)^2]
   errorsum = errorsum + error
   if (error > biggest)
      biggest = error
//   println[error earthradius]
}

println["Largest error: " + biggest earthradius]
average = errorsum/samples
println["Average error: " + average earthradius]

   


View or download UTMtest.frink in plain text format


This is a program written in the programming language Frink.
For more information, view the Frink Documentation or see More Sample Frink Programs.

Alan Eliasen was born 17775 days, 11 hours, 24 minutes ago.