I'm trying to use OWSLib to deal with multiple WSF features in python, for now, I'm trying to return all features in sfb WFS that intersects with the specifc feature retrieved in sicar WFS. from owslib.wfs import WebFeatureService, Authentication from owslib.fes2 import PropertyIsEqualTo, Intersects from owslib.etree import etree auth = Authentication() auth.verify = False sicar = WebFeatureService( url="https://geoserver.car.gov.br/geoserver/sicar/ows", auth=auth, version='2.0.0') sfb = WebFeatureService( 'http://sistemas.florestal.gov.br/geoserver/ows', version='2.0.0') car = 'MT-5103700-92A39DC62FDF4D88874F218D41BE1C7F' filter_car = PropertyIsEqualTo('cod_imovel', car).toXML() carInfo = sicar.getfeature(typename='sicar:sicar_imoveis_mt', filter=filter_car, method='post', outputFormat='json') The code above is working fine, the problem starts when i try to get the Intersects geom_property = 'geo_area_imovel'; carInfoStr = carInfo.read().decode('utf-8'); geom = json.loads(carInfoStr)['features'][0]['geometry']['coordinates'] #I'm trying to pass the the coordinates list to Intersects, but i tried with other keys of the object and with the xml response, but doesn't wotk as well. filter_contains = Intersects(geom_property, geom).toXML() intersected = sfb.getfeature(typename='CNFP_orig:APP_Total_simp', filter=filter_contains, outputFormat='json', method='post') The returned errors involve the toXML function, wich depends on the geometry i pass to the intersects, in this case its a string, but need to be the constructor of something of the owslib, but the lib only have Point constructor. I could pass a point if I could retrieve the centroid of geometry, but I don't now how to do this with this lib. BUT, even the example passed in a official PR is not working. the PR with example: https://github.com/geopython/OWSLib/pull/780 The error in this code: Traceback (most recent call last): File "/home/gabrielubuntu/repos/pygeogis/app.py", line 38, in filter_contains = Intersects(geom_property, geom).toXML() File "/home/gabrielubuntu/repos/pygeogis/.venv/lib/python3.10/site-packages/owslib/fes2.py", line 420, in toXML node.append(self.geometry.toXML()) AttributeError: 'dict' object has no attribute 'toXML'

I'm trying to download a lot of datasets from different WFS URLs using ogr2ogr (and to an extent, owslib) in Python. My problem arises when I try to connect to a second WFS. In my test code, I have first tried connecting to one WFS, then downloading the layer. This went fine. Afterwards, I replaced the URL with a second WFS URL. this time I got the error: ConnectionError: ('Connection aborted.', ConnectionResetError(10054 An existing connection was forcibly closed by the remote host)). In the first part of the code I try to extract the version number from the WFS: import subprocess from owslib.wfs import WebFeatureService as wfsmod #connecting to second WFS using owslib url = 'http://wfs2-miljoegis.mim.dk/raastofpaahavet/ows? version=1.1.0&OUTPUTFORMAT=GML2' wfs = wfsmod(url) version = (wfs.identification.version) typenam ="rashavtilladelser20auktion" The connection error happens at line "wfs = wfsmod(url)". I also try to download a specific layer from this url using ogr2ogr/subprocesses: if str(version) == "2.0.0": #layername is the specific title of the layer on the WFS server. layername = "&TYPENAMES="+layername else: layername = "&TYPENAME="+layername cmd = [r"C:\OSGeo4W\bin\ogr2ogr.exe", "-f", "ESRI Shapefile", r"C:\havplan\havplan_data\wfs_downloads\havbrugszoner.shp", "WFS:"+url+"version="+str(version) +"&request=GetFeature"+layername] subprocess.check_call(cmd) as far as I understand, the error happens because I can only connect to one WFS at a time in Python. Is there some way to close the connection to a WFS in owslib and ogr2ogr using a command?

https://www.osgeo.org/foundation-news/tom-kralidis-receives-the-2024-sol-katz-award/ Tom Kralidis was honoured with the 2024 Sol Katz Award, presented on 6 December 2024 at FOSS4G 2024 in Belém, Brazil. This was the 20th year of the award. 2024-12-09 We are honoured to announce that Tom Kralidis is the recipient of the 2024 Sol Katz Award. You may know Tom through his activities within the GeoPython community, specifically through pygeoapi, pycsw, OWSLib, GeoHealthCheck, pygeometa, and th...

Ubuntu OWSLib security update fixes sensitive info exposure issue. Ensure systems are updated to prevent data leaks.

I need to extract attribute data from all features from a WMS source. So far I've managed to get them "manually" specifying coordinates for each feature using Python and OWSlib like this: from owslib.wms import WebMapService wms = WebMapService() fInfo = wms.getfeatureinfo( layers=['0'], styles=['default'], srs='EPSG:4326', bbox=(-68.68876,-23.09997,-67.82525,-21.12957), # manually set feature coordinates size=(500,500), format='image/jpeg', query_layers=['0'], info_format='application/geojson', xy=(250,250) ) out = open(wms['0'].title.replace(" ", "_")+".json", "wb") out.write(fInfo.read()) out.close() So my problem is exactly that: I need to know beforehand where are all the features and specify a bounding box inside them. Is there a way to list all features or objects from a WMS and download their metadata?

I'm trying to use OWSLib to deal with multiple WSF features in python, for now, I'm trying to return all features in sfb WFS that intersects with the specifc feature retrieved in sicar WFS. from owslib.wfs import WebFeatureService, Authentication from owslib.fes2 import PropertyIsEqualTo, Intersects from owslib.etree import etree auth = Authentication() auth.verify = False sicar = WebFeatureService( url="https://geoserver.car.gov.br/geoserver/sicar/ows", auth=auth, version='2.0.0') sfb = WebFeatureService( 'http://sistemas.florestal.gov.br/geoserver/ows', version='2.0.0') car = 'MT-5103700-92A39DC62FDF4D88874F218D41BE1C7F' filter_car = PropertyIsEqualTo('cod_imovel', car).toXML() carInfo = sicar.getfeature(typename='sicar:sicar_imoveis_mt', filter=filter_car, method='post', outputFormat='json') The code above is working fine, the problem starts when i try to get the Intersects geom_property = 'geo_area_imovel'; carInfoStr = carInfo.read().decode('utf-8'); geom = json.loads(carInfoStr)['features'][0]['geometry']['coordinates'] #I'm trying to pass the the coordinates list to Intersects, but i tried with other keys of the object and with the xml response, but doesn't wotk as well. filter_contains = Intersects(geom_property, geom).toXML() intersected = sfb.getfeature(typename='CNFP_orig:APP_Total_simp', filter=filter_contains, outputFormat='json', method='post') The returned errors involve the toXML function, wich depends on the geometry i pass to the intersects, in this case its a string, but need to be the constructor of something of the owslib, but the lib only have Point constructor. I could pass a point if I could retrieve the centroid of geometry, but I don't now how to do this with this lib. BUT, even the example passed in a official PR is not working. the PR with example: https://github.com/geopython/OWSLib/pull/780 The error in this code: Traceback (most recent call last): File "/home/gabrielubuntu/repos/pygeogis/app.py", line 38, in filter_contains = Intersects(geom_property, geom).toXML() File "/home/gabrielubuntu/repos/pygeogis/.venv/lib/python3.10/site-packages/owslib/fes2.py", line 420, in toXML node.append(self.geometry.toXML()) AttributeError: 'dict' object has no attribute 'toXML'

I'm trying to download a lot of datasets from different WFS URLs using ogr2ogr (and to an extent, owslib) in Python. My problem arises when I try to connect to a second WFS. In my test code, I have first tried connecting to one WFS, then downloading the layer. This went fine. Afterwards, I replaced the URL with a second WFS URL. this time I got the error: ConnectionError: ('Connection aborted.', ConnectionResetError(10054 An existing connection was forcibly closed by the remote host)). In the first part of the code I try to extract the version number from the WFS: import subprocess from owslib.wfs import WebFeatureService as wfsmod #connecting to second WFS using owslib url = 'http://wfs2-miljoegis.mim.dk/raastofpaahavet/ows? version=1.1.0&OUTPUTFORMAT=GML2' wfs = wfsmod(url) version = (wfs.identification.version) typenam ="rashavtilladelser20auktion" The connection error happens at line "wfs = wfsmod(url)". I also try to download a specific layer from this url using ogr2ogr/subprocesses: if str(version) == "2.0.0": #layername is the specific title of the layer on the WFS server. layername = "&TYPENAMES="+layername else: layername = "&TYPENAME="+layername cmd = [r"C:\OSGeo4W\bin\ogr2ogr.exe", "-f", "ESRI Shapefile", r"C:\havplan\havplan_data\wfs_downloads\havbrugszoner.shp", "WFS:"+url+"version="+str(version) +"&request=GetFeature"+layername] subprocess.check_call(cmd) as far as I understand, the error happens because I can only connect to one WFS at a time in Python. Is there some way to close the connection to a WFS in owslib and ogr2ogr using a command?

I need to extract attribute data from all features from a WMS source. So far I've managed to get them "manually" specifying coordinates for each feature using Python and OWSlib like this: from owslib.wms import WebMapService wms = WebMapService() fInfo = wms.getfeatureinfo( layers=['0'], styles=['default'], srs='EPSG:4326', bbox=(-68.68876,-23.09997,-67.82525,-21.12957), # manually set feature coordinates size=(500,500), format='image/jpeg', query_layers=['0'], info_format='application/geojson', xy=(250,250) ) out = open(wms['0'].title.replace(" ", "_")+".json", "wb") out.write(fInfo.read()) out.close() So my problem is exactly that: I need to know beforehand where are all the features and specify a bounding box inside them. Is there a way to list all features or objects from a WMS and download their metadata?

I'm trying to use OWSLib to deal with multiple WSF features in python, for now, I'm trying to return all features in sfb WFS that intersects with the specifc feature retrieved in sicar WFS. from owslib.wfs import WebFeatureService, Authentication from owslib.fes2 import PropertyIsEqualTo, Intersects from owslib.etree import etree auth = Authentication() auth.verify = False sicar = WebFeatureService( url="https://geoserver.car.gov.br/geoserver/sicar/ows", auth=auth, version='2.0.0') sfb = WebFeatureService( 'http://sistemas.florestal.gov.br/geoserver/ows', version='2.0.0') car = 'MT-5103700-92A39DC62FDF4D88874F218D41BE1C7F' filter_car = PropertyIsEqualTo('cod_imovel', car).toXML() carInfo = sicar.getfeature(typename='sicar:sicar_imoveis_mt', filter=filter_car, method='post', outputFormat='json') The code above is working fine, the problem starts when i try to get the Intersects geom_property = 'geo_area_imovel'; carInfoStr = carInfo.read().decode('utf-8'); geom = json.loads(carInfoStr)['features'][0]['geometry']['coordinates'] #I'm trying to pass the the coordinates list to Intersects, but i tried with other keys of the object and with the xml response, but doesn't wotk as well. filter_contains = Intersects(geom_property, geom).toXML() intersected = sfb.getfeature(typename='CNFP_orig:APP_Total_simp', filter=filter_contains, outputFormat='json', method='post') The returned errors involve the toXML function, wich depends on the geometry i pass to the intersects, in this case its a string, but need to be the constructor of something of the owslib, but the lib only have Point constructor. I could pass a point if I could retrieve the centroid of geometry, but I don't now how to do this with this lib. BUT, even the example passed in a official PR is not working. the PR with example: https://github.com/geopython/OWSLib/pull/780 The error in this code: Traceback (most recent call last): File "/home/gabrielubuntu/repos/pygeogis/app.py", line 38, in filter_contains = Intersects(geom_property, geom).toXML() File "/home/gabrielubuntu/repos/pygeogis/.venv/lib/python3.10/site-packages/owslib/fes2.py", line 420, in toXML node.append(self.geometry.toXML()) AttributeError: 'dict' object has no attribute 'toXML'

I'm trying to download a lot of datasets from different WFS URLs using ogr2ogr (and to an extent, owslib) in Python. My problem arises when I try to connect to a second WFS. In my test code, I have first tried connecting to one WFS, then downloading the layer. This went fine. Afterwards, I replaced the URL with a second WFS URL. this time I got the error: ConnectionError: ('Connection aborted.', ConnectionResetError(10054 An existing connection was forcibly closed by the remote host)). In the first part of the code I try to extract the version number from the WFS: import subprocess from owslib.wfs import WebFeatureService as wfsmod #connecting to second WFS using owslib url = 'http://wfs2-miljoegis.mim.dk/raastofpaahavet/ows? version=1.1.0&OUTPUTFORMAT=GML2' wfs = wfsmod(url) version = (wfs.identification.version) typenam ="rashavtilladelser20auktion" The connection error happens at line "wfs = wfsmod(url)". I also try to download a specific layer from this url using ogr2ogr/subprocesses: if str(version) == "2.0.0": #layername is the specific title of the layer on the WFS server. layername = "&TYPENAMES="+layername else: layername = "&TYPENAME="+layername cmd = [r"C:\OSGeo4W\bin\ogr2ogr.exe", "-f", "ESRI Shapefile", r"C:\havplan\havplan_data\wfs_downloads\havbrugszoner.shp", "WFS:"+url+"version="+str(version) +"&request=GetFeature"+layername] subprocess.check_call(cmd) as far as I understand, the error happens because I can only connect to one WFS at a time in Python. Is there some way to close the connection to a WFS in owslib and ogr2ogr using a command?

I need to extract attribute data from all features from a WMS source. So far I've managed to get them "manually" specifying coordinates for each feature using Python and OWSlib like this: from owslib.wms import WebMapService wms = WebMapService() fInfo = wms.getfeatureinfo( layers=['0'], styles=['default'], srs='EPSG:4326', bbox=(-68.68876,-23.09997,-67.82525,-21.12957), # manually set feature coordinates size=(500,500), format='image/jpeg', query_layers=['0'], info_format='application/geojson', xy=(250,250) ) out = open(wms['0'].title.replace(" ", "_")+".json", "wb") out.write(fInfo.read()) out.close() So my problem is exactly that: I need to know beforehand where are all the features and specify a bounding box inside them. Is there a way to list all features or objects from a WMS and download their metadata?

I'm trying to use OWSLib to deal with multiple WSF features in python, for now, I'm trying to return all features in sfb WFS that intersects with the specifc feature retrieved in sicar WFS. from owslib.wfs import WebFeatureService, Authentication from owslib.fes2 import PropertyIsEqualTo, Intersects from owslib.etree import etree auth = Authentication() auth.verify = False sicar = WebFeatureService( url="https://geoserver.car.gov.br/geoserver/sicar/ows", auth=auth, version='2.0.0') sfb = WebFeatureService( 'http://sistemas.florestal.gov.br/geoserver/ows', version='2.0.0') car = 'MT-5103700-92A39DC62FDF4D88874F218D41BE1C7F' filter_car = PropertyIsEqualTo('cod_imovel', car).toXML() carInfo = sicar.getfeature(typename='sicar:sicar_imoveis_mt', filter=filter_car, method='post', outputFormat='json') The code above is working fine, the problem starts when i try to get the Intersects geom_property = 'geo_area_imovel'; carInfoStr = carInfo.read().decode('utf-8'); geom = json.loads(carInfoStr)['features'][0]['geometry']['coordinates'] #I'm trying to pass the the coordinates list to Intersects, but i tried with other keys of the object and with the xml response, but doesn't wotk as well. filter_contains = Intersects(geom_property, geom).toXML() intersected = sfb.getfeature(typename='CNFP_orig:APP_Total_simp', filter=filter_contains, outputFormat='json', method='post') The returned errors involve the toXML function, wich depends on the geometry i pass to the intersects, in this case its a string, but need to be the constructor of something of the owslib, but the lib only have Point constructor. I could pass a point if I could retrieve the centroid of geometry, but I don't now how to do this with this lib. BUT, even the example passed in a official PR is not working. the PR with example: https://github.com/geopython/OWSLib/pull/780 The error in this code: Traceback (most recent call last): File "/home/gabrielubuntu/repos/pygeogis/app.py", line 38, in filter_contains = Intersects(geom_property, geom).toXML() File "/home/gabrielubuntu/repos/pygeogis/.venv/lib/python3.10/site-packages/owslib/fes2.py", line 420, in toXML node.append(self.geometry.toXML()) AttributeError: 'dict' object has no attribute 'toXML'

I'm trying to download a lot of datasets from different WFS URLs using ogr2ogr (and to an extent, owslib) in Python. My problem arises when I try to connect to a second WFS. In my test code, I have first tried connecting to one WFS, then downloading the layer. This went fine. Afterwards, I replaced the URL with a second WFS URL. this time I got the error: ConnectionError: ('Connection aborted.', ConnectionResetError(10054 An existing connection was forcibly closed by the remote host)). In the first part of the code I try to extract the version number from the WFS: import subprocess from owslib.wfs import WebFeatureService as wfsmod #connecting to second WFS using owslib url = 'http://wfs2-miljoegis.mim.dk/raastofpaahavet/ows? version=1.1.0&OUTPUTFORMAT=GML2' wfs = wfsmod(url) version = (wfs.identification.version) typenam ="rashavtilladelser20auktion" The connection error happens at line "wfs = wfsmod(url)". I also try to download a specific layer from this url using ogr2ogr/subprocesses: if str(version) == "2.0.0": #layername is the specific title of the layer on the WFS server. layername = "&TYPENAMES="+layername else: layername = "&TYPENAME="+layername cmd = [r"C:\OSGeo4W\bin\ogr2ogr.exe", "-f", "ESRI Shapefile", r"C:\havplan\havplan_data\wfs_downloads\havbrugszoner.shp", "WFS:"+url+"version="+str(version) +"&request=GetFeature"+layername] subprocess.check_call(cmd) as far as I understand, the error happens because I can only connect to one WFS at a time in Python. Is there some way to close the connection to a WFS in owslib and ogr2ogr using a command?

I'm trying to get proj strings from cartopy crs', such as cartopy.crs.PlateCarree(). I can manually infer the proj string from the init() function in the source code but I need to query it, so I can apply it to a tiff file with gdal. I'm trying to avoid adding modules outside from scipy, numpy, gdal, cartopy, OWSLib, MatPlotLib, and Pillow. Is there a solution?

I need to extract attribute data from all features from a WMS source. So far I've managed to get them "manually" specifying coordinates for each feature using Python and OWSlib like this: from owslib.wms import WebMapService wms = WebMapService() fInfo = wms.getfeatureinfo( layers=['0'], styles=['default'], srs='EPSG:4326', bbox=(-68.68876,-23.09997,-67.82525,-21.12957), # manually set feature coordinates size=(500,500), format='image/jpeg', query_layers=['0'], info_format='application/geojson', xy=(250,250) ) out = open(wms['0'].title.replace(" ", "_")+".json", "wb") out.write(fInfo.read()) out.close() So my problem is exactly that: I need to know beforehand where are all the features and specify a bounding box inside them. Is there a way to list all features or objects from a WMS and download their metadata?

I'm trying to use OWSLib to deal with multiple WSF features in python, for now, I'm trying to return all features in sfb WFS that intersects with the specifc feature retrieved in sicar WFS. from owslib.wfs import WebFeatureService, Authentication from owslib.fes2 import PropertyIsEqualTo, Intersects from owslib.etree import etree auth = Authentication() auth.verify = False sicar = WebFeatureService( url="https://geoserver.car.gov.br/geoserver/sicar/ows", auth=auth, version='2.0.0') sfb = WebFeatureService( 'http://sistemas.florestal.gov.br/geoserver/ows', version='2.0.0') car = 'MT-5103700-92A39DC62FDF4D88874F218D41BE1C7F' filter_car = PropertyIsEqualTo('cod_imovel', car).toXML() carInfo = sicar.getfeature(typename='sicar:sicar_imoveis_mt', filter=filter_car, method='post', outputFormat='json') The code above is working fine, the problem starts when i try to get the Intersects geom_property = 'geo_area_imovel'; carInfoStr = carInfo.read().decode('utf-8'); geom = json.loads(carInfoStr)['features'][0]['geometry']['coordinates'] #I'm trying to pass the the coordinates list to Intersects, but i tried with other keys of the object and with the xml response, but doesn't wotk as well. filter_contains = Intersects(geom_property, geom).toXML() intersected = sfb.getfeature(typename='CNFP_orig:APP_Total_simp', filter=filter_contains, outputFormat='json', method='post') The returned errors involve the toXML function, wich depends on the geometry i pass to the intersects, in this case its a string, but need to be the constructor of something of the owslib, but the lib only have Point constructor. I could pass a point if I could retrieve the centroid of geometry, but I don't now how to do this with this lib. BUT, even the example passed in a official PR is not working. the PR with example: https://github.com/geopython/OWSLib/pull/780 The error in this code: Traceback (most recent call last): File "/home/gabrielubuntu/repos/pygeogis/app.py", line 38, in filter_contains = Intersects(geom_property, geom).toXML() File "/home/gabrielubuntu/repos/pygeogis/.venv/lib/python3.10/site-packages/owslib/fes2.py", line 420, in toXML node.append(self.geometry.toXML()) AttributeError: 'dict' object has no attribute 'toXML'

I'm trying to download a lot of datasets from different WFS URLs using ogr2ogr (and to an extent, owslib) in Python. My problem arises when I try to connect to a second WFS. In my test code, I have first tried connecting to one WFS, then downloading…

I need to extract attribute data from all features from a WMS source. So far I've managed to get them "manually" specifying coordinates for each feature using Python and OWSlib like this: from owslib.wms import WebMapService wms = WebMapService() fInfo…

I'm trying to use OWSLib to deal with multiple WSF features in python, for now, I'm trying to return all features in sfb WFS that intersects with the specifc feature retrieved in sicar WFS. from owslib.wfs import WebFeatureService, Authentication from…