def testBoundingSpherePrecision(self):
x = 533
y = 383
z = 9
geodetic = GlobalGeodetic(True)
[minx, miny, maxx, maxy] = geodetic.TileBounds(x, y, z)
ter = TerrainTile(west=minx, south=miny, east=maxx, north=maxy)
ter.fromFile('tests/data/%s_%s_%s.terrain' % (z, x, y))
coords = [LLH2ECEF(*coord) for coord in ter.getVerticesCoordinates()]
sphere = BoundingSphere()
sphere.fromPoints(coords)
for coord in coords:
distance = c3d.distance(sphere.center, coord)
self.assertLessEqual(distance, sphere.radius)
def testBoundingSpherePrecision(self):
x = 533
y = 383
z = 9
geodetic = GlobalGeodetic(True)
[minx, miny, maxx, maxy] = geodetic.TileBounds(x, y, z)
ter = TerrainTile(west=minx, south=miny, east=maxx, north=maxy)
ter.fromFile('tests/data/%s_%s_%s.terrain' % (z, x, y))
llh2ecef = lambda x: LLH2ECEF(x[0], x[1], x[2])
coords = ter.getVerticesCoordinates()
coords = list(map(llh2ecef, coords))
sphere = BoundingSphere()
sphere.fromPoints(coords)
for coord in coords:
distance = c3d.distance(sphere.center, coord)
self.assertLessEqual(distance, sphere.radius)
def testReaderWriter(self):
'''
Circle jerk testing.
We read the file with our reader
We write this data with our writer to a temporary file
We read this temporary file
We compare the results
'''
x = 533
y = 383
z = 9
geodetic = GlobalGeodetic(True)
[minx, miny, maxx, maxy] = geodetic.TileBounds(x, y, z)
ter = TerrainTile(west=minx, south=miny, east=maxx, north=maxy)
ter.fromFile('tests/data/%s_%s_%s.terrain' % (z, x, y))
ter.toFile(self.tmpfile)
self.assertIsInstance(ter.__repr__(), str)
ter2 = TerrainTile(west=minx, south=miny, east=maxx, north=maxy)
ter2.fromFile(self.tmpfile)
self.assertIsInstance(ter2.__repr__(), str)
# check headers
self.assertGreater(len(ter.header), 0)
self.assertEqual(len(ter.header), len(ter2.header))
self.assertEqual(len(ter.header), len(TerrainTile.quantizedMeshHeader))
for k, v in ter.header.items():
self.assertEqual(v, ter2.header[k], 'For k = ' + k)
# check vertices
self.assertGreater(len(ter.u), 0)
self.assertGreater(len(ter.v), 0)
self.assertGreater(len(ter.h), 0)
self.assertEqual(len(ter.u), len(ter2.u))
self.assertEqual(len(ter.v), len(ter2.v))
self.assertEqual(len(ter.h), len(ter2.h))
for i, v in enumerate(ter.u):
self.assertEqual(v, ter2.u[i])
for i, v in enumerate(ter.v):
self.assertEqual(v, ter2.v[i])
for i, v in enumerate(ter.h):
self.assertEqual(v, ter2.h[i])
self.assertEqual(
len(ter.getVerticesCoordinates()),
len(ter2.getVerticesCoordinates())
)
# check indices
self.assertGreater(len(ter.indices), 0)
self.assertEqual(len(ter.indices), len(ter2.indices))
for i, v in enumerate(ter.indices):
self.assertEqual(v, ter2.indices[i], i)
# check edges
self.assertGreater(len(ter.westI), 0)
self.assertEqual(len(ter.westI), len(ter2.westI))
for i, v in enumerate(ter.westI):
self.assertEqual(v, ter2.westI[i], i)
self.assertGreater(len(ter.southI), 0)
self.assertEqual(len(ter.southI), len(ter2.southI))
for i, v in enumerate(ter.southI):
self.assertEqual(v, ter2.southI[i], i)
self.assertGreater(len(ter.eastI), 0)
self.assertEqual(len(ter.eastI), len(ter2.eastI))
for i, v in enumerate(ter.eastI):
self.assertEqual(v, ter2.eastI[i], i)
self.assertGreater(len(ter.northI), 0)
self.assertEqual(len(ter.northI), len(ter2.northI))
for i, v in enumerate(ter.northI):
self.assertEqual(v, ter2.northI[i], i)
self.assertEqual(ter2.getContentType(),
'application/vnd.quantized-mesh')
def testReaderWriter(self):
'''
Circle jerk testing.
We read the file with our reader
We write this data with our writer to a temporary file
We read this temporary file
We compare the results
'''
x = 533
y = 383
z = 9
geodetic = GlobalGeodetic(True)
[minx, miny, maxx, maxy] = geodetic.TileBounds(x, y, z)
ter = TerrainTile(west=minx, south=miny, east=maxx, north=maxy)
ter.fromFile('tests/data/%s_%s_%s.terrain' % (z, x, y))
ter.toFile(self.tmpfile)
self.assertIsInstance(ter.__repr__(), str)
ter2 = TerrainTile(west=minx, south=miny, east=maxx, north=maxy)
ter2.fromFile(self.tmpfile)
self.assertIsInstance(ter2.__repr__(), str)
# check headers
self.assertGreater(len(ter.header), 0)
self.assertEqual(len(ter.header), len(ter2.header))
self.assertEqual(len(ter.header), len(TerrainTile.quantizedMeshHeader))
for k, v in ter.header.iteritems():
self.assertEqual(v, ter2.header[k], 'For k = ' + k)
# check vertices
self.assertGreater(len(ter.u), 0)
self.assertGreater(len(ter.v), 0)
self.assertGreater(len(ter.h), 0)
self.assertEqual(len(ter.u), len(ter2.u))
self.assertEqual(len(ter.v), len(ter2.v))
self.assertEqual(len(ter.h), len(ter2.h))
for i, v in enumerate(ter.u):
self.assertEqual(v, ter2.u[i])
for i, v in enumerate(ter.v):
self.assertEqual(v, ter2.v[i])
for i, v in enumerate(ter.h):
self.assertEqual(v, ter2.h[i])
self.assertEqual(
len(ter.getVerticesCoordinates()),
len(ter2.getVerticesCoordinates())
)
# check indices
self.assertGreater(len(ter.indices), 0)
self.assertEqual(len(ter.indices), len(ter2.indices))
for i, v in enumerate(ter.indices):
self.assertEqual(v, ter2.indices[i], i)
# check edges
self.assertGreater(len(ter.westI), 0)
self.assertEqual(len(ter.westI), len(ter2.westI))
for i, v in enumerate(ter.westI):
self.assertEqual(v, ter2.westI[i], i)
self.assertGreater(len(ter.southI), 0)
self.assertEqual(len(ter.southI), len(ter2.southI))
for i, v in enumerate(ter.southI):
self.assertEqual(v, ter2.southI[i], i)
self.assertGreater(len(ter.eastI), 0)
self.assertEqual(len(ter.eastI), len(ter2.eastI))
for i, v in enumerate(ter.eastI):
self.assertEqual(v, ter2.eastI[i], i)
self.assertGreater(len(ter.northI), 0)
self.assertEqual(len(ter.northI), len(ter2.northI))
for i, v in enumerate(ter.northI):
self.assertEqual(v, ter2.northI[i], i)
self.assertEqual(ter2.getContentType(),
'application/vnd.quantized-mesh')