def test(self):
"""
Test GenSimpleDBApp with 1-D data in 2-D space.
"""
from spatialdata.spatialdb.generator.GenSimpleDBApp import GenSimpleDBApp
from spatialdata.geocoords.CSCart import CSCart
from spatialdata.spatialdb.SimpleDB import SimpleDB
app = GenSimpleDBApp()
app.run()
# Test write using query
db = SimpleDB()
db.inventory.iohandler.inventory.filename = \
"data/gen1Din2D.spatialdb"
db.inventory.iohandler._configure()
db.inventory.label = "test"
db.inventory.queryType = "nearest"
db._configure()
qlocs = numpy.array( [[-2.0, 2.0],
[ 3.0, -4.0],
[ 0.0, 0.7]],
numpy.float64)
dataE = numpy.array( [[-0.15, 3.45],
[2.4, 6.4],
[-0.6, 3.45]], numpy.float64)
errE = [0, 0, 0]
from spatialdata.geocoords.CSCart import CSCart
cs = CSCart()
cs.inventory.spaceDim = 2
cs._configure()
db.open()
db.queryVals(["two", "one"])
data = numpy.zeros(dataE.shape, dtype=numpy.float64)
err = []
nlocs = qlocs.shape[0]
for i in xrange(nlocs):
e = db.query(data[i,:], qlocs[i,:], cs)
err.append(e)
db.close()
self.assertEqual(len(errE), len(err))
for vE, v in zip(errE, err):
self.assertEqual(vE, v)
self.assertEqual(len(dataE.shape), len(data.shape))
for dE, d in zip(dataE.shape, data.shape):
self.assertEqual(dE, d)
for vE, v in zip(numpy.reshape(dataE, -1), numpy.reshape(data, -1)):
self.assertAlmostEqual(vE, v, 6)
return
def test(self):
"""
Test GenSimpleDBApp with 1-D data in 2-D space.
"""
from spatialdata.spatialdb.generator.GenSimpleDBApp import GenSimpleDBApp
from spatialdata.geocoords.CSCart import CSCart
from spatialdata.spatialdb.SimpleDB import SimpleDB
app = GenSimpleDBApp()
app.run()
# Test write using query
db = SimpleDB()
db.inventory.iohandler.inventory.filename = "data/gen1Din2D.spatialdb"
db.inventory.iohandler._configure()
db.inventory.label = "test"
db.inventory.queryType = "nearest"
db._configure()
qlocs = numpy.array([[-2.0, 2.0], [3.0, -4.0], [0.0, 0.7]],
numpy.float64)
dataE = numpy.array([[-0.15, 3.45], [2.4, 6.4], [-0.6, 3.45]],
numpy.float64)
errE = [0, 0, 0]
from spatialdata.geocoords.CSCart import CSCart
cs = CSCart()
cs.inventory.spaceDim = 2
cs._configure()
db.open()
db.setQueryValues(["two", "one"])
data = numpy.zeros(dataE.shape, dtype=numpy.float64)
err = []
nlocs = qlocs.shape[0]
for i in range(nlocs):
e = db.query(data[i, :], qlocs[i, :], cs)
err.append(e)
db.close()
self.assertEqual(len(errE), len(err))
for vE, v in zip(errE, err):
self.assertEqual(vE, v)
self.assertEqual(len(dataE.shape), len(data.shape))
for dE, d in zip(dataE.shape, data.shape):
self.assertEqual(dE, d)
for vE, v in zip(numpy.reshape(dataE, -1), numpy.reshape(data, -1)):
self.assertAlmostEqual(vE, v, 6)
#!/usr/bin/env nemesis
#
# ======================================================================
#
# Brad T. Aagaard, U.S. Geological Survey
#
# This code was developed as part of the Computational Infrastructure
# for Geodynamics (http://geodynamics.org).
#
# Copyright (c) 2010-2022 University of California, Davis
#
# See LICENSE.md for license information.
#
# ======================================================================
#
# @file spatialdata/applications/gensimpledb.py
#
# @brief Generic use of GenSimpleDBApp.
# ----------------------------------------------------------------------
if __name__ == '__main__':
from spatialdata.spatialdb.generator.GenSimpleDBApp import GenSimpleDBApp
app = GenSimpleDBApp()
app.run()
# End of file
#!/usr/bin/env python # # ====================================================================== # # Brad T. Aagaard, U.S. Geological Survey # # This code was developed as part of the Computational Infrastructure # for Geodynamics (http://geodynamics.org). # # Copyright (c) 2010-2017 University of California, Davis # # See COPYING for license information. # # ====================================================================== # ## @file spatialdata/applications/gensimpledb.py ## @brief Generic use of GenSimpleDBApp. # ---------------------------------------------------------------------- if __name__ == '__main__': from spatialdata.spatialdb.generator.GenSimpleDBApp import GenSimpleDBApp app = GenSimpleDBApp() app.run() # End of file
