Scribd
Upload
26 views1 page

Point Test

The document contains code for testing transformations between coordinate systems. It defines functions for calculating distance and creating point objects in different coordinate systems (XYZ, geocentric, geodetic). The tests validate that distances are equal to 0 when transforming between the different coordinate systems for sample data, and when adding navigational offsets.

Uploaded by

kpalamartchouk
Copyright
© Attribution Non-Commercial (BY-NC)
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
26 views1 page

Point Test

The document contains code for testing transformations between coordinate systems. It defines functions for calculating distance and creating point objects in different coordinate systems (XYZ, geocentric, geodetic). The tests validate that distances are equal to 0 when transforming between the different coordinate systems for sample data, and when adding navigational offsets.

Uploaded by

kpalamartchouk
Copyright
© Attribution Non-Commercial (BY-NC)
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 1

#!/usr/bin/python us=0).add_neu(matrix(((0, 0, 1))).T) − matrix((( 0, 0, 1))).

T), 0)
# self.assertAlmostEqual(dist(point(latitude = 0, longitude = 270, earthRadi
import unittest us=0).add_neu(matrix(((1, 0, 0))).T) − matrix((( 0, 0, 1))).T), 0)
from point import point # self.assertAlmostEqual(dist(point(latitude = 0, longitude = 270, earthRadi
from math import radians us=0).add_neu(matrix(((0, 1, 0))).T) − matrix((( 1, 0, 0))).T), 0)
# self.assertAlmostEqual(dist(point(latitude = 0, longitude = 270, earthRadi
def dist(xx, yy): us=0).add_neu(matrix(((0, 0, 1))).T) − matrix((( 0,−1, 0))).T), 0)
return sum((x − y)**2 for (x, y) in zip(xx, yy)) # self.assertAlmostEqual(dist(point(latitude = 90, longitude = 270, earthRadi
us=0).add_neu(matrix(((1, 0, 0))).T) − matrix((( 0, 1, 0))).T), 0)
class transformationsTest(unittest.TestCase): # self.assertAlmostEqual(dist(point(latitude = 90, longitude = 270, earthRadi
def test_xyz_to_geo_real_data(self): us=0).add_neu(matrix(((0, 1, 0))).T) − matrix((( 1, 0, 0))).T), 0)
WSRT_XYZ__TRUE = (3828735.77084, 443304.96442, 5064884.71882) # self.assertAlmostEqual(dist(point(latitude = 90, longitude = 270, earthRadi
WSRT_GEOC_TRUE = (radians(52.729318000), radians(6.604505677), 6364651.9237) us=0).add_neu(matrix(((0, 0, 1))).T) − matrix((( 0, 0, 1))).T), 0)
WSRT_GEOD_TRUE = (radians(52.91461174), radians(6.60450568), 82.2314)

WSRT_XYZ = point(WSRT_XYZ__TRUE, kind=’xyz’)


WSRT_GEOC = point(WSRT_GEOC_TRUE, kind=’geoc’)
WSRT_GEOD = point(WSRT_GEOD_TRUE, kind=’geod’)
self.assertAlmostEqual(dist(WSRT_XYZ, WSRT_GEOC), 0.0, places=6)
self.assertAlmostEqual(dist(WSRT_XYZ, WSRT_GEOD), 0.0, places=6)

self.assertAlmostEqual(dist(WSRT_XYZ.geoc, WSRT_GEOC_TRUE), 0.0, places=6)


self.assertAlmostEqual(dist(WSRT_XYZ.geod, WSRT_GEOD_TRUE), 0.0, places=6)
if __name__ == "__main__":
unittest.main()
#
#class transformationsTest(unittest.TestCase):
# def test_add_neu(self):
# self.assertAlmostEqual(dist(point(latitude = 0, longitude = 0, earthRadi
us=0).add_neu(matrix(((1, 0, 0))).T) − matrix((( 0, 0, 1))).T), 0)
# self.assertAlmostEqual(dist(point(latitude = 0, longitude = 0, earthRadi
us=0).add_neu(matrix(((0, 1, 0))).T) − matrix((( 0, 1, 0))).T), 0)
# self.assertAlmostEqual(dist(point(latitude = 0, longitude = 0, earthRadi
us=0).add_neu(matrix(((0, 0, 1))).T) − matrix((( 1, 0, 0))).T), 0)
# self.assertAlmostEqual(dist(point(latitude = 90, longitude = 0, earthRadi
us=0).add_neu(matrix(((1, 0, 0))).T) − matrix(((−1, 0, 0))).T), 0)
# self.assertAlmostEqual(dist(point(latitude = 90, longitude = 0, earthRadi
us=0).add_neu(matrix(((0, 1, 0))).T) − matrix((( 0, 1, 0))).T), 0)
# self.assertAlmostEqual(dist(point(latitude = 90, longitude = 0, earthRadi
us=0).add_neu(matrix(((0, 0, 1))).T) − matrix((( 0, 0, 1))).T), 0)
# self.assertAlmostEqual(dist(point(latitude = 0, longitude = 90, earthRadi
us=0).add_neu(matrix(((1, 0, 0))).T) − matrix((( 0, 0, 1))).T), 0)
# self.assertAlmostEqual(dist(point(latitude = 0, longitude = 90, earthRadi
us=0).add_neu(matrix(((0, 1, 0))).T) − matrix(((−1, 0, 0))).T), 0)
# self.assertAlmostEqual(dist(point(latitude = 0, longitude = 90, earthRadi
us=0).add_neu(matrix(((0, 0, 1))).T) − matrix((( 0, 1, 0))).T), 0)
# self.assertAlmostEqual(dist(point(latitude = 90, longitude = 90, earthRadi
us=0).add_neu(matrix(((1, 0, 0))).T) − matrix((( 0,−1, 0))).T), 0)
# self.assertAlmostEqual(dist(point(latitude = 90, longitude = 90, earthRadi
us=0).add_neu(matrix(((0, 1, 0))).T) − matrix(((−1, 0, 0))).T), 0)
# self.assertAlmostEqual(dist(point(latitude = 90, longitude = 90, earthRadi
us=0).add_neu(matrix(((0, 0, 1))).T) − matrix((( 0, 0, 1))).T), 0)
# self.assertAlmostEqual(dist(point(latitude = 0, longitude = 180, earthRadi
us=0).add_neu(matrix(((1, 0, 0))).T) − matrix((( 0, 0, 1))).T), 0)
# self.assertAlmostEqual(dist(point(latitude = 0, longitude = 180, earthRadi
us=0).add_neu(matrix(((0, 1, 0))).T) − matrix((( 0,−1, 0))).T), 0)
# self.assertAlmostEqual(dist(point(latitude = 0, longitude = 180, earthRadi
us=0).add_neu(matrix(((0, 0, 1))).T) − matrix(((−1, 0, 0))).T), 0)
# self.assertAlmostEqual(dist(point(latitude = 90, longitude = 180, earthRadi
us=0).add_neu(matrix(((1, 0, 0))).T) − matrix((( 1, 0, 0))).T), 0)
# self.assertAlmostEqual(dist(point(latitude = 90, longitude = 180, earthRadi
us=0).add_neu(matrix(((0, 1, 0))).T) − matrix((( 0,−1, 0))).T), 0)
# self.assertAlmostEqual(dist(point(latitude = 90, longitude = 180, earthRadi

You might also like