DM_spatial_query.py

## @package python.demo.DM_spatial_query
# "
# example script showing how to load an OPALS datamanager 
# and perform a simple spatial query
# (script works in python 2 and 3)
# 
 
from __future__ import print_function  # print function syntax as in python 3
import sys
 
from opals import pyDM
 
def DM_spatial_query(filename, maxOutput=1000):
  # pyDM.Datamanager.load parameter: filename(string), readOnly(bool) threadSafety(bool)
  dm = pyDM.Datamanager.load(filename, True, False)
  if not dm:
    print("Unable to open ODM '" + filename + "'")
    sys.exit(1)
 
  limit = dm.getLimit()
 
  # output some pyDM statistics
  print("ODM contains", dm.sizePoint(), "points")
  print("2D-limit (%.3f," % limit.xmin, "%.3f) -" % limit.ymin, "(%.3f," % limit.xmax, "%.3f)" % limit.ymax)
 
  # spatial query test
  print("Perform spatial query")
  pts = dm.searchPoint(limit, True)  # search all points inside manager limit
  assert len(pts) == dm.sizePoint()  # check if points found is equal to points within manager
  print(len(pts), "points found")
 
  print("Output points...")
  i = 0
  for pt in pts:
    i = i + 1
    print("Point", "%5d" % i, "%.3f" % pt.x, "%.3f" % pt.y, "%.3f" % pt.z)
    if i == maxOutput:
      break
 
  print(i, "of", len(pts), "points listed")
 
if len(sys.argv) == 1:
    print("ODM parameter missing")
    sys.exit(-1)
 
DM_spatial_query(sys.argv[1])