def test_multiline_single_31370(self):
"""Test the WithinDistance filter with a GML containing a single
multiline geometry in EPSG:31370.
Test whether the generated XML is correct.
"""
with open('tests/data/util/location/multiline_single_31370.gml',
'r') as gml_file:
gml = gml_file.read()
f = GmlFilter(gml, WithinDistance, {'distance': 100})
f.set_geometry_column('geom')
xml = f.toXML()
assert clean_xml(etree.tostring(xml).decode('utf8')) == clean_xml(
'<ogc:DWithin><ogc:PropertyName>geom</ogc:PropertyName><gml'
':MultiCurve srsName="urn:ogc:def:crs:EPSG::31370"><gml'
':curveMember><gml:LineString><gml:posList>108210.673003586 '
'194850.535439444 108454.928328293 195031.757131969 '
'108746.458877137 194834.777031398</gml:posList></gml:LineString'
'></gml:curveMember><gml:curveMember><gml:LineString><gml'
':posList>109164.056690347 195055.394744037 109211.331914484 '
'194661.434542896 109416.191219077 '
'194440.816830258</gml:posList></gml:LineString></gml'
':curveMember></gml:MultiCurve><gml:Distance '
'units="meter">100.000000</gml:Distance></ogc:DWithin>')
def test_polyon_multiple_31370_stable(self):
"""Test the Within filter with a GML containing multiple
polygon geometries in EPSG:31370.
Test whether the generated XML is correct and stable.
"""
with open('tests/data/util/location/polygon_multiple_31370.gml',
'r') as gml_file:
gml = gml_file.read()
for i in range(10):
f = GmlFilter(gml, Within)
f.set_geometry_column('geom')
xml = f.toXML()
assert clean_xml(etree.tostring(xml).decode('utf8')) == clean_xml(
'<ogc:Or><ogc:Within><ogc:PropertyName>geom</ogc'
':PropertyName><gml:Polygon '
'srsName="urn:ogc:def:crs:EPSG::31370"><gml:exterior'
'><gml:LinearRing><gml:posList>108636.150020818 '
'194960.844295764 108911.922161617 194291.111953824 '
'109195.573506438 195118.42837622 108636.150020818 '
'194960.844295764</gml:posList></gml:LinearRing></gml'
':exterior></gml:Polygon></ogc:Within><ogc:Within><ogc'
':PropertyName>geom</ogc:PropertyName><gml:Polygon '
'srsName="urn:ogc:def:crs:EPSG::31370"><gml:exterior'
'><gml:LinearRing><gml:posList>107485.786233486 '
'196741.544404921 107840.350414513 196339.704999757 '
'108297.344247837 196843.974057217 107485.786233486 '
'196741.544404921</gml:posList></gml:LinearRing></gml'
':exterior></gml:Polygon></ogc:Within></ogc:Or>')
def test_multipolygon_single_31370(self):
"""Test the Within filter with a GML containing a single
multipolygon geometry in EPSG:31370.
Test whether the generated XML is correct.
"""
with open('tests/data/util/location/multipolygon_single_31370.gml',
'r') as gml_file:
gml = gml_file.read()
f = GmlFilter(gml, Within)
f.set_geometry_column('geom')
xml = f.toXML()
assert clean_xml(etree.tostring(xml).decode('utf8')) == clean_xml(
'<ogc:Within><ogc:PropertyName>geom</ogc:PropertyName><gml'
':MultiSurface srsName="urn:ogc:def:crs:EPSG::31370"><gml'
':surfaceMember><gml:Polygon><gml:exterior><gml:LinearRing><gml'
':posList>108588.874796681 195015.998723923 108911.922161617 '
'194251.71593371 109195.573506438 195134.186784266 '
'108588.874796681 195015.998723923</gml:posList></gml:LinearRing'
'></gml:exterior></gml:Polygon></gml:surfaceMember><gml'
':surfaceMember><gml:Polygon><gml:exterior><gml:LinearRing><gml'
':posList>109140.419078278 195748.764698045 109400.432811031 '
'195307.529272768 109597.412911602 195772.402310114 '
'109140.419078278 195748.764698045</gml:posList></gml:LinearRing'
'></gml:exterior></gml:Polygon></gml:surfaceMember></gml'
':MultiSurface></ogc:Within>')
def test_line_multiple_31370_stable(self):
"""Test the WithinDistance filter with a GML containing multiple
line geometries in EPSG:31370.
Test whether the generated XML is correct and stable.
"""
with open('tests/data/util/location/line_multiple_31370.gml',
'r') as gml_file:
gml = gml_file.read()
for i in range(10):
f = GmlFilter(gml, WithinDistance, {'distance': 100})
f.set_geometry_column('geom')
xml = f.toXML()
assert clean_xml(etree.tostring(xml).decode('utf8')) == clean_xml(
'<ogc:Or><ogc:DWithin><ogc:PropertyName>geom</ogc'
':PropertyName><gml :LineString '
'srsName="urn:ogc:def:crs:EPSG::31370"><gml:posList>108344'
'.619471974 195008.119519901 108801.613305297 '
'194842.656235421 109077.385446096 '
'195094.790764152</gml:posList></gml:LineString><gml'
':Distance units="meter">100.000000</gml:Distance></ogc'
':DWithin><ogc:DWithin><ogc:PropertyName>geom</ogc'
':PropertyName><gml :LineString '
'srsName="urn:ogc:def:crs:EPSG::31370"><gml:posList>108194'
'.91459554 196032.416042867 108714.942061046 '
'195528.146985407 108911.922161617 '
'195528.146985407</gml:posList></gml:LineString><gml'
':Distance units="meter">100.000000</gml:Distance>'
'</ogc:DWithin></ogc:Or>')
def test_point_multiple_31370_stable(self):
"""Test the WithinDistance filter with a GML containing multiple
point geometries in EPSG:31370.
Test whether the generated XML is correct and stable.
"""
with open('tests/data/util/location/point_multiple_31370.gml',
'r') as gml_file:
gml = gml_file.read()
for i in range(10):
f = GmlFilter(gml, WithinDistance, {'distance': 100})
f.set_geometry_column('geom')
xml = f.toXML()
assert clean_xml(etree.tostring(xml).decode('utf8')) == clean_xml(
'<ogc:Or><ogc:DWithin><ogc:PropertyName>geom</ogc'
':PropertyName><gml:Point '
'srsName="urn:ogc:def:crs:EPSG::31370"><gml:pos>109124'
'.660670233 194937.206683695</gml:pos></gml:Point><gml'
':Distance units="meter">100.000000</gml:Distance></ogc'
':DWithin><ogc:DWithin><ogc:PropertyName>geom</ogc'
':PropertyName><gml:Point '
'srsName="urn:ogc:def:crs:EPSG::31370"><gml:pos>109234'
'.969526552 195772.402310114</gml:pos></gml:Point><gml'
':Distance units="meter">100.000000</gml:Distance>'
'</ogc:DWithin></ogc:Or>')
def test_multipolygon_multiple_31370_stable(self):
"""Test the Within filter with a GML containing multiple
multipolygon geometries in EPSG:31370.
Test whether the generated XML is correct.
"""
with open('tests/data/util/location/multipolygon_multiple_31370.gml',
'r') as gml_file:
gml = gml_file.read()
for i in range(10):
f = GmlFilter(gml, Within)
f.set_geometry_column('geom')
xml = f.toXML()
assert clean_xml(etree.tostring(xml).decode('utf8')) == clean_xml(
'<ogc:Or><ogc:Within><ogc:PropertyName>geom</ogc'
':PropertyName><gml:MultiSurface '
'srsName="urn:ogc:def:crs:EPSG::31370"><gml:surfaceMember'
'><gml:Polygon><gml:exterior><gml:LinearRing><gml:posList'
'>107564.578273715 196646.993956647 107785.195986354 '
'196386.980223894 107966.417678878 197143.383810084 '
'107564.578273715 '
'196646.993956647</gml:posList></gml:LinearRing'
'></gml:exterior></gml:Polygon></gml:surfaceMember><gml'
':surfaceMember><gml:Polygon><gml:exterior><gml:LinearRing'
'><gml:posList>108384.015492088 197214.29664629 '
'108447.04912427 196489.40987619 108785.854897252 '
'197269.45107445 108384.015492088 '
'197214.29664629</gml:posList></gml:LinearRing'
'></gml:exterior></gml:Polygon></gml:surfaceMember></gml'
':MultiSurface></ogc:Within><ogc:Within><ogc:PropertyName'
'>geom</ogc:PropertyName><gml:MultiSurface '
'srsName="urn:ogc:def:crs:EPSG::31370"><gml:surfaceMember'
'><gml:Polygon><gml:exterior><gml:LinearRing><gml:posList'
'>108588.874796681 195015.998723923 108911.922161617 '
'194251.71593371 109195.573506438 195134.186784266 '
'108588.874796681 '
'195015.998723923</gml:posList></gml:LinearRing'
'></gml:exterior></gml:Polygon></gml:surfaceMember><gml'
':surfaceMember><gml:Polygon><gml:exterior><gml:LinearRing'
'><gml:posList>109140.419078278 195748.764698045 '
'109400.432811031 195307.529272768 109597.412911602 '
'195772.402310114 109140.419078278 '
'195748.764698045</gml:posList></gml:LinearRing'
'></gml:exterior></gml:Polygon></gml:surfaceMember></gml'
':MultiSurface></ogc:Within></ogc:Or>')
def test_polygon_and_box(self):
"""Test the combination of a Within filter with a GML containing a
single polygon geometry in EPSG:31370 and a Box.
Test whether the generated XML is correct.
"""
with open('tests/data/util/location/polygon_single_31370.gml',
'r') as gml_file:
gml = gml_file.read()
polygon = GmlFilter(gml, Within)
box = Box(94720, 186910, 112220, 202870)
location_filter = Or([polygon, Within(box)])
set_geometry_column(location_filter, 'geom')
xml = location_filter.toXML()
assert clean_xml(etree.tostring(xml).decode('utf8')) == clean_xml(
'<ogc:Or><ogc:Within><ogc:PropertyName>geom</ogc:PropertyName'
'><gml:Polygon srsName="urn:ogc:def:crs:EPSG::31370"><gml'
':exterior><gml:LinearRing><gml:posList>108636.150020818 '
'194960.844295764 108911.922161617 194291.111953824 '
'109195.573506438 195118.42837622 108636.150020818 '
'194960.844295764</gml:posList></gml:LinearRing></gml:exterior'
'></gml:Polygon></ogc:Within><ogc:Within><ogc:PropertyName>geom'
'</ogc:PropertyName><gml:Envelope srsDimension="2" '
'srsName="http://www.opengis.net/gml/srs/epsg.xml#31370"><gml'
':lowerCorner>94720.000000 '
'186910.000000</gml:lowerCorner><gml:upperCorner>112220.000000 '
'202870.000000</gml:upperCorner></gml:Envelope></ogc'
':Within></ogc:Or>')
def test_multiline_multiple_31370_stable(self):
"""Test the WithinDistance filter with a GML containing multiple
multiline geometries in EPSG:31370.
Test whether the generated XML is correct.
"""
with open('tests/data/util/location/multiline_multiple_31370.gml',
'r') as gml_file:
gml = gml_file.read()
for i in range(10):
f = GmlFilter(gml, WithinDistance, {'distance': 100})
f.set_geometry_column('geom')
xml = f.toXML()
assert clean_xml(etree.tostring(xml).decode('utf8')) == clean_xml(
'<ogc:Or><ogc:DWithin><ogc:PropertyName>geom</ogc'
':PropertyName><gml:MultiCurve '
'srsName="urn:ogc:def:crs:EPSG::31370"><gml:curveMember><gml'
':LineString><gml:posList>108210.673003586 194850.535439444 '
'108454.928328293 195031.757131969 108746.458877137 '
'194834.777031398</gml:posList></gml:LineString></gml'
':curveMember><gml:curveMember><gml:LineString><gml:posList'
'>109164.056690347 195055.394744037 109211.331914484 '
'194661.434542896 109416.191219077 '
'194440.816830258</gml:posList></gml:LineString></gml'
':curveMember></gml:MultiCurve><gml:Distance '
'units="meter">100.000000</gml:Distance></ogc:DWithin><ogc'
':DWithin><ogc:PropertyName>geom</ogc:PropertyName><gml'
':MultiCurve '
'srsName="urn:ogc:def:crs:EPSG::31370"><gml:curveMember><gml'
':LineString><gml:posList>108226.431411631 196095.44967505 '
'108384.015492088 196276.671367574 108580.995592658 '
'196048.174450913</gml:posList></gml:LineString></gml'
':curveMember><gml:curveMember><gml:LineString><gml:posList'
'>108911.922161617 196379.101019871 109030.110221959 '
'196552.443508373 109282.244750689 '
'196607.597936533</gml:posList></gml:LineString></gml'
':curveMember></gml:MultiCurve><gml:Distance '
'units="meter">100.000000</gml:Distance></ogc:DWithin>'
'</ogc:Or>')
def test_point_single_31370(self):
"""Test the WithinDistance filter with a GML containing a single
point geometry in EPSG:31370.
Test whether the generated XML is correct.
"""
with open('tests/data/util/location/point_single_31370.gml',
'r') as gml_file:
gml = gml_file.read()
f = GmlFilter(gml, WithinDistance, {'distance': 100})
f.set_geometry_column('geom')
xml = f.toXML()
assert clean_xml(etree.tostring(xml).decode('utf8')) == clean_xml(
'<ogc:DWithin><ogc:PropertyName>geom</ogc:PropertyName><gml'
':Point srsName="urn:ogc:def:crs:EPSG::31370"><gml:pos>109124'
'.660670233 194937.206683695</gml:pos></gml:Point><gml:Distance '
'units="meter">100.000000</gml:Distance></ogc:DWithin>')
def test_point_multiple_31370(self):
"""Test the WithinDistance filter with a GML containing multiple
point geometries in EPSG:31370.
Test whether the generated XML is correct.
"""
with open('tests/data/util/location/point_multiple_31370.gml',
'r') as gml_file:
gml = gml_file.read()
f = GmlFilter(gml, WithinDistance, {'distance': 100})
f.set_geometry_column('geom')
xml = f.toXML()
assert xml.tag == '{http://www.opengis.net/ogc}Or'
assert len(list(xml)) == 2
points = [
'109124.660670233 194937.206683695',
'109234.969526552 195772.402310114'
]
for f in xml:
assert f.tag == '{http://www.opengis.net/ogc}DWithin'
point = f.find('./{http://www.opengis.net/gml}Point')
assert point.get('srsName') == 'urn:ogc:def:crs:EPSG::31370'
posList = point.find('./{http://www.opengis.net/gml}pos').text
assert posList in points
points.remove(posList)
assert len(points) == 0
def test_line_single_31370(self):
"""Test the WithinDistance filter with a GML containing a single
line geometry in EPSG:31370.
Test whether the generated XML is correct.
"""
with open('tests/data/util/location/line_single_31370.gml',
'r') as gml_file:
gml = gml_file.read()
f = GmlFilter(gml, WithinDistance, {'distance': 100})
f.set_geometry_column('geom')
xml = f.toXML()
assert clean_xml(etree.tostring(xml).decode('utf8')) == clean_xml(
'<ogc:DWithin><ogc:PropertyName>geom</ogc:PropertyName><gml'
':LineString srsName="urn:ogc:def:crs:EPSG::31370"><gml:posList'
'>108344.619471974 195008.119519901 108801.613305297 '
'194842.656235421 109077.385446096 '
'195094.790764152</gml:posList></gml:LineString><gml:Distance '
'units="meter">100.000000</gml:Distance></ogc:DWithin>')
def test_polygon_single_31370(self):
"""Test the Within filter with a GML containing a single
polygon geometry in EPSG:31370.
Test whether the generated XML is correct.
"""
with open('tests/data/util/location/polygon_single_31370.gml',
'r') as gml_file:
gml = gml_file.read()
f = GmlFilter(gml, Within)
f.set_geometry_column('geom')
xml = f.toXML()
assert clean_xml(etree.tostring(xml).decode('utf8')) == clean_xml(
'<ogc:Within><ogc:PropertyName>geom</ogc:PropertyName><gml'
':Polygon srsName="urn:ogc:def:crs:EPSG::31370"><gml:exterior'
'><gml:LinearRing><gml:posList>108636.150020818 194960.844295764 '
'108911.922161617 194291.111953824 109195.573506438 '
'195118.42837622 108636.150020818 '
'194960.844295764</gml:posList></gml:LinearRing></gml:exterior'
'></gml:Polygon></ogc:Within>')
def test_multipoint_single_31370(self):
"""Test the WithinDistance filter with a GML containing a single
multipoint geometry in EPSG:31370.
Test whether the generated XML is correct.
"""
with open('tests/data/util/location/multipoint_single_31370.gml',
'r') as gml_file:
gml = gml_file.read()
f = GmlFilter(gml, WithinDistance, {'distance': 100})
f.set_geometry_column('geom')
xml = f.toXML()
assert clean_xml(etree.tostring(xml).decode('utf8')) == clean_xml(
'<ogc:DWithin><ogc:PropertyName>geom</ogc:PropertyName><gml'
':MultiPoint srsName="urn:ogc:def:crs:EPSG::31370"><gml'
':pointMember><gml:Point><gml:pos>108770.096489206 '
'194992.361111855</gml:pos></gml:Point></gml:pointMember><gml'
':pointMember><gml:Point><gml:pos>109045.868630005 '
'194929.327479672</gml:pos></gml:Point></gml:pointMember></gml'
':MultiPoint><gml:Distance '
'units="meter">100.000000</gml:Distance></ogc:DWithin>')
def test_multipoint_multiple_31370_stable(self):
"""Test the WithinDistance filter with a GML containing multiple
multipoint geometries in EPSG:31370.
Test whether the generated XML is correct and stable.
"""
with open('tests/data/util/location/multipoint_multiple_31370.gml',
'r') as gml_file:
gml = gml_file.read()
for i in range(10):
f = GmlFilter(gml, WithinDistance, {'distance': 100})
f.set_geometry_column('geom')
xml = f.toXML()
assert clean_xml(etree.tostring(xml).decode('utf8')) == clean_xml(
'<ogc:Or><ogc:DWithin><ogc:PropertyName>geom</ogc'
':PropertyName><gml:MultiPoint '
'srsName="urn:ogc:def:crs:EPSG::31370"><gml:pointMember><gml'
':Point><gml:pos>108770.096489206 '
'194992.361111855</gml:pos></gml:Point></gml:pointMember'
'><gml:pointMember><gml:Point><gml:pos>109045.868630005 '
'194929.327479672</gml:pos></gml:Point></gml:pointMember'
'></gml:MultiPoint><gml:Distance '
'units="meter">100.000000</gml:Distance></ogc:DWithin><ogc'
':DWithin><ogc:PropertyName>geom</ogc:PropertyName><gml'
':MultiPoint '
'srsName="urn:ogc:def:crs:EPSG::31370"><gml:pointMember><gml'
':Point><gml:pos>108825.250917366 '
'195433.596537133</gml:pos></gml:Point></gml:pointMember'
'><gml:pointMember><gml:Point><gml:pos>108738.579673115 '
'195614.818229658</gml:pos></gml:Point></gml:pointMember'
'></gml:MultiPoint><gml:Distance '
'units="meter">100.000000</gml:Distance></ogc:DWithin>'
'</ogc:Or>')
def test_multipolygon_multiple_31370(self):
"""Test the Within filter with a GML containing multiple
multipolygon geometries in EPSG:31370.
Test whether the generated XML is correct.
"""
with open('tests/data/util/location/multipolygon_multiple_31370.gml',
'r') as gml_file:
gml = gml_file.read()
f = GmlFilter(gml, Within)
f.set_geometry_column('geom')
xml = f.toXML()
assert xml.tag == '{http://www.opengis.net/ogc}Or'
assert len(list(xml)) == 2
multipolygons = [
[
'107564.578273715 196646.993956647 107785.195986354 '
'196386.980223894 107966.417678878 197143.383810084 '
'107564.578273715 196646.993956647',
'108384.015492088 197214.29664629 108447.04912427 '
'196489.40987619 108785.854897252 197269.45107445 '
'108384.015492088 197214.29664629'
],
[
'108588.874796681 195015.998723923 108911.922161617 '
'194251.71593371 109195.573506438 195134.186784266 '
'108588.874796681 195015.998723923',
'109140.419078278 195748.764698045 109400.432811031 '
'195307.529272768 109597.412911602 195772.402310114 '
'109140.419078278 195748.764698045'
]
]
for f in xml:
assert f.tag == '{http://www.opengis.net/ogc}Within'
multicurve = f.find('./{http://www.opengis.net/gml}MultiSurface')
assert multicurve.get('srsName') == 'urn:ogc:def:crs:EPSG::31370'
posList = [
p.text for p in multicurve.findall(
'.//{http://www.opengis.net/gml}posList')
]
multipolygons.remove(posList)
assert len(multipolygons) == 0
def test_multiline_multiple_31370(self):
"""Test the WithinDistance filter with a GML containing multiple
multiline geometries in EPSG:31370.
Test whether the generated XML is correct.
"""
with open('tests/data/util/location/multiline_multiple_31370.gml',
'r') as gml_file:
gml = gml_file.read()
f = GmlFilter(gml, WithinDistance, {'distance': 100})
f.set_geometry_column('geom')
xml = f.toXML()
assert xml.tag == '{http://www.opengis.net/ogc}Or'
assert len(list(xml)) == 2
multilines = [[
'108210.673003586 194850.535439444 108454.928328293 '
'195031.757131969 108746.458877137 194834.777031398',
'109164.056690347 195055.394744037 109211.331914484 '
'194661.434542896 109416.191219077 194440.816830258'
],
[
'108226.431411631 196095.44967505 108384.015492088 '
'196276.671367574 108580.995592658 196048.174450913',
'108911.922161617 196379.101019871 109030.110221959 '
'196552.443508373 109282.244750689 196607.597936533'
]]
for f in xml:
assert f.tag == '{http://www.opengis.net/ogc}DWithin'
multicurve = f.find('./{http://www.opengis.net/gml}MultiCurve')
assert multicurve.get('srsName') == 'urn:ogc:def:crs:EPSG::31370'
posList = [
p.text for p in multicurve.findall(
'.//{http://www.opengis.net/gml}posList')
]
multilines.remove(posList)
assert len(multilines) == 0
def test_multipoint_multiple_31370(self):
"""Test the WithinDistance filter with a GML containing multiple
multipoint geometries in EPSG:31370.
Test whether the generated XML is correct.
"""
with open('tests/data/util/location/multipoint_multiple_31370.gml',
'r') as gml_file:
gml = gml_file.read()
f = GmlFilter(gml, WithinDistance, {'distance': 100})
f.set_geometry_column('geom')
xml = f.toXML()
assert xml.tag == '{http://www.opengis.net/ogc}Or'
assert len(list(xml)) == 2
multipoints = [[
'108770.096489206 194992.361111855',
'109045.868630005 194929.327479672'
],
[
'108825.250917366 195433.596537133',
'108738.579673115 195614.818229658'
]]
for f in xml:
assert f.tag == '{http://www.opengis.net/ogc}DWithin'
multipoint = f.find('./{http://www.opengis.net/gml}MultiPoint')
assert multipoint.get('srsName') == 'urn:ogc:def:crs:EPSG::31370'
posList = [
p.text for p in multipoint.findall(
'.//{http://www.opengis.net/gml}pos')
]
multipoints.remove(posList)
assert len(multipoints) == 0
def test_polyon_multiple_disjoint_31370(self):
"""Test the Disjoint filter with the And combinator with a GML
containing multiple polygon geometries in EPSG:31370.
Test whether the generated XML is correct.
"""
with open('tests/data/util/location/polygon_multiple_31370.gml',
'r') as gml_file:
gml = gml_file.read()
f = GmlFilter(gml, Disjoint, combinator=And)
f.set_geometry_column('geom')
xml = f.toXML()
assert xml.tag == '{http://www.opengis.net/ogc}And'
assert len(list(xml)) == 2
polygons = [
'108636.150020818 194960.844295764 108911.922161617 '
'194291.111953824 109195.573506438 195118.42837622 '
'108636.150020818 194960.844295764',
'107485.786233486 196741.544404921 107840.350414513 '
'196339.704999757 108297.344247837 196843.974057217 '
'107485.786233486 196741.544404921'
]
for f in xml:
assert f.tag == '{http://www.opengis.net/ogc}Disjoint'
point = f.find('./{http://www.opengis.net/gml}Polygon')
assert point.get('srsName') == 'urn:ogc:def:crs:EPSG::31370'
posList = point.find('.//{http://www.opengis.net/gml}posList').text
assert posList in polygons
polygons.remove(posList)
assert len(polygons) == 0
def test_line_multiple_31370(self):
"""Test the WithinDistance filter with a GML containing multiple
line geometries in EPSG:31370.
Test whether the generated XML is correct.
"""
with open('tests/data/util/location/line_multiple_31370.gml',
'r') as gml_file:
gml = gml_file.read()
f = GmlFilter(gml, WithinDistance, {'distance': 100})
f.set_geometry_column('geom')
xml = f.toXML()
assert xml.tag == '{http://www.opengis.net/ogc}Or'
assert len(list(xml)) == 2
lines = [
'108344.619471974 195008.119519901 108801.613305297 '
'194842.656235421 109077.385446096 195094.790764152',
'108194.91459554 196032.416042867 108714.942061046 '
'195528.146985407 108911.922161617 195528.146985407'
]
for f in xml:
assert f.tag == '{http://www.opengis.net/ogc}DWithin'
point = f.find('./{http://www.opengis.net/gml}LineString')
assert point.get('srsName') == 'urn:ogc:def:crs:EPSG::31370'
posList = point.find('./{http://www.opengis.net/gml}posList').text
assert posList in lines
lines.remove(posList)
assert len(lines) == 0