QGIS » QGIS Application

iface = qgis.utils.iface

rasterLayer = iface.addRasterLayer("crs=EPSG:4326&dpiMode=7&featureCount=10&format=image/png&layers=natura2000&styles=&url=http://geodata.nationaalgeoregister.nl/natura2000/ows", "wms master example", "wms" ) 

vectorLayer = iface.addVectorLayer("http://geodata.nationaalgeoregister.nl/bestuurlijkegrenzen/wfs?SERVICE=WFS&VERSION=1.0.0&REQUEST=GetFeature&TYPENAME=bestuurlijkegrenzen:provincies&SRSNAME=EPSG:28992","wfs example","WFS")

#QgsMapLayerRegistry.instance().addMapLayer(vectorLayer)

iter = vectorLayer.getFeatures()

for feature in iter:

    geom=feature.geometry()

    print "Feature ID %d: " % feature.id()

    if geom.type() == QGis.Point:

        x = geom.asPoint()

        print "Point: " + str(x)

    elif geom.type() == QGis.Line:

        x = geom.asPolyline()

        print "Line: %d points" % len(x)

    elif geom.type() == QGis.Polygon:

        x = geom.asPolygon()

        numPts = 0

        for ring in x:

            numPts += len(ring)

        print "Polygon: %d rings with %d points" % (len(x), numPts)

    else:

        print "Unknown"

    #fetch attributes

    #attrs = feature.attributes()

    #print attrs

QgsMapLayerRegistry.instance().removeMapLayer(vectorLayer.id())