Instituto Brasileiro do Meio Ambiente e dos Recursos Naturais Renováveis (IBAMA)¶

In [ ]:

import io

import geopandas as gpd
from owslib.wfs import WebFeatureService

import io import geopandas as gpd from owslib.wfs import WebFeatureService

In [ ]:

wfs = WebFeatureService(
    url='https://siscom.ibama.gov.br/geoserver/ows/',
    version='2.0.0',
)

wfs = WebFeatureService( url='https://siscom.ibama.gov.br/geoserver/ows/', version='2.0.0', )

In [ ]:

list(wfs.contents)

list(wfs.contents)

Sistema Integrado de Gestão Ambiental (SIGA)

In [ ]:

for layer_name, layer in wfs.contents.items():
    # if layer.queryable == 0:
    print(f"Layer: {layer_name}")
    print(f"  Title: {layer.title}")
    print(f"  Abstract: {layer.abstract}")
    print(f"  BoundingBox: {layer.boundingBoxWGS84}")
    print(f"  CRS: {layer.crsOptions}")
    print(f"  Styles: {layer.styles}")
    print(f"  Keywords: {layer.keywords}")
    # print(f"  Queryable: {layer.queryable}")
    # print(f"  Opaque: {layer.opaque}")
    # print(f"  Dimensions: {layer.dimensions}")
    print(f"  MetadataURLs: {layer.metadataUrls}")
    print()

for layer_name, layer in wfs.contents.items(): # if layer.queryable == 0: print(f"Layer: {layer_name}") print(f" Title: {layer.title}") print(f" Abstract: {layer.abstract}") print(f" BoundingBox: {layer.boundingBoxWGS84}") print(f" CRS: {layer.crsOptions}") print(f" Styles: {layer.styles}") print(f" Keywords: {layer.keywords}") # print(f" Queryable: {layer.queryable}") # print(f" Opaque: {layer.opaque}") # print(f" Dimensions: {layer.dimensions}") print(f" MetadataURLs: {layer.metadataUrls}") print()

---

Obter Layers¶

In [ ]:

# Obter os dados no formato GeoJSON (ou outro formato suportado)
response = wfs.getfeature(
    typename='siga:siga_multipolygon',
    # bbox=(173700, 440400, 178700, 441400),
    # srsname='EPSG:28992'
    # srsname='EPSG:4326',
    srsname='EPSG:4674',
    outputFormat='application/json',
)
response

# Obter os dados no formato GeoJSON (ou outro formato suportado) response = wfs.getfeature( typename='siga:siga_multipolygon', # bbox=(173700, 440400, 178700, 441400), # srsname='EPSG:28992' # srsname='EPSG:4326', srsname='EPSG:4674', outputFormat='application/json', ) response

In [ ]:

gdf = gpd.read_file(filename=io.BytesIO(response.read()))
gdf.crs

gdf = gpd.read_file(filename=io.BytesIO(response.read())) gdf.crs

In [ ]:

gdf.info()
gdf.head(2)

gdf.info() gdf.head(2)

In [ ]:

gdf.explore()

gdf.explore()

In [ ]: