access_WFS_with_python

# -*- coding: utf-8 -*-
"""
Created on Mon Mar 22 18:13:38 2021

@author: mlc
"""

#import fiona; # help(fiona.open)
#import shapely as st
#import matplotlib.pyplot as plt
#import numpy as np
#from matplotlib import interactive, style

import pandas as pd
import geopandas as geopd
import os
from time import process_time 
from owslib.wfs import WebFeatureService # <-------------------------------------------------
# a instalação no conda tem erro: fazer pip install --force-reinstall pyproj
# também tive 
import pyproj
from requests import Request


# test pyproj  # https://stackoverflow.com/questions/55390492/runtimeerror-bno-arguments-in-initialization-list
pyproj.Proj("+init=epsg:3857")

#C:\Users\mlc\OneDrive - Universidade de Lisboa\Documents\geomatica-sigdr-2020-2021\Trabalho_avaliacao\InvFlorNacional
os.chdir("C:\\Users\\mlc\\OneDrive - Universidade de Lisboa\\Documents\\geomatica-sigdr-2020-2021\\Trabalho_avaliacao\\InvFlorNacional")
os.getcwd()

# URL for WFS backend
url = "http://si.icnf.pt/wfs/rnap"
url = "http://si.icnf.pt/wfs/ifn_2015"

# See details about this particular WFS
# -------------------------------------

# Initialize <--------------------------------------------  create connection
wfs = WebFeatureService(url=url)

# Service provider 
print(wfs.identification.title)

# Get WFS version
print(wfs.version)

# Available methods
print([operation.name for operation in wfs.operations])

# Available data layers
print(list(wfs.contents))

# Print all metadata of all layers
for layer, meta in wfs.items():
    print(meta.__dict__)

# Get data from WFS
# -----------------

# Fetch the last available layer (as an example) --> 'vaestoruutu:vaki2017_5km'
layer = list(wfs.contents)[-1]

# Specify the parameters for fetching the data
params = dict(service='WFS', version="1.0.0", request='GetFeature', typeName=layer, outputFormat='json')

# Parse the URL with parameters # <----------------------------- create request
q = Request('GET', url, params=params).prepare().url

# Read data from URL  # <--------------------------------------  read data with geopandas
data = geopd.read_file(q)

data.columns