queryDemo.py

import opals
from osgeo import ogr
from opals import Query
 
qry = Query.Query()
qry.commons.screenLogLevel = opals.Types.LogLevel.error
qry.inFile = ["strip21_ne.tif","strip21_nw.tif","strip21_se.tif","strip21_sw.tif"]
qry.inGeometry = ["LINESTRING(529635 5338727, 529635 5338626, 529763 5338626, 529826 5338711)"]
qry.run()
outGeom = qry.outGeometry
line3d = ogr.CreateGeometryFromWkt(outGeom[0])
print("\nResulting 3D line geometry: {}".format(outGeom[0]))
print("\nEnvelope 2D: {}".format(line3d.GetEnvelope()))
print("\nEnvelope 3D: {}".format(line3d.GetEnvelope3D()))
print("\nLine length: {}".format(line3d.Length()))