When opening QGIS and add connection to the ArcGIS Feature/Map Service like :https://gis.adem.alabama.gov/arcgis/rest/services it giving me the list of all services and the layers within each service.

Is there a way to do it using PyQGIS ?

enter image description here

I know to load layer with:

QgsVectorLayer(uri, layer_name, "arcgisfeatureserver")

I want to get all the layers available on a source using PyQGIS.

1 Answer 1


On the webpage of the REST services URL there is a json link at the top left. Clicking this will show a json of the service metadata, including a list of names. Likewise, a similar link will show a json of all the layers in each service.

enter image description here

The service names and corresponding layers can be extracted like so:

import urllib.request, json 

json_url = 'https://gis.adem.alabama.gov/arcgis/rest/services?f=pjson'

# download the services json
with urllib.request.urlopen(json_url) as url:
    data = json.loads(url.read().decode())
# get services from the json
services = data['services']

# get names and types of services
for i,x in enumerate(services):
    print(i, ' | ',  x['name'], ' | ', x['type'])

# number of service
num = 5

# get name and type from list 
service_name = services[num]['name']
service_type = services[num]['type']

# get layers from service
layers_url = f"https://gis.adem.alabama.gov/arcgis/rest/services/{service_name}/{service_type}/layers?f=pjson"

# download the layers json
with urllib.request.urlopen(layers_url) as url:
    data = json.loads(url.read().decode())

layers = data['layers']

available_layers = []
for layer in layers:
    available_layers.append([layer[k] for k in ['id', 'name', 'type', 'geometryType']])

[print(x) for x in available_layers]

## pick a layer 
layer_index = 0

# check that the chosen index is available
if layer_index in (x[0] for x in available_layers):
    # build the REST service url
    uri = f"crs='EPSG:4326' url='https://gis.adem.alabama.gov/arcgis/rest/services/{service_name}/{service_type}/{layer_index}"
    # make QGIS layer
    lyr = QgsVectorLayer(uri, layer['name'], "arcgisfeatureserver")
    # add to map
    print(f"chosen layer index ({layer_index}) is not in available layers")

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge that you have read and understand our privacy policy and code of conduct.

Not the answer you're looking for? Browse other questions tagged or ask your own question.