#!/usr/bin/env python

import Numeric
import MLab
import string

f = open('NED_test.out', 'r')

count_missing_information = 0
count_missing_coverage = 0
count_total = 0

differences_list = []
for line in f.readlines():
	count_total += 1
	if 'missing coverage' in line:
		count_missing_coverage += 1
		continue
	if 'missing information' in line:
		count_missing_information += 1
		continue

	try:
		(ident, lon, lat, FAA, NED, difference_s) = string.split(line, ' ')
	except:
		continue

	difference = float(difference_s)

	differences_list.append(difference)

	if abs(difference) > 100:
		print line,

differences = Numeric.array(differences_list)
print 'count_total=%d' % (count_total)
print 'count_missing_coverage=%d' % (count_missing_coverage)
print 'count_missing_information=%d' % (count_missing_information)
print 'differences min=%f' % (MLab.min(differences))
print 'differences max=%f' % (MLab.max(differences))
print 'differences mean=%f' % (MLab.mean(differences))
print 'differences standard deviation=%f' % (MLab.std(differences))