Przykład kodu w Python zdobywającego URL arkuszy Ortofotomapy z geoportal.gov.pl - wersja bardzo wstępna

requirements.txt

certifi==2020.6.20
chardet==3.0.4
idna==2.10
pyshp==2.1.2
requests==2.24.0
urllib3==1.25.10

script.py

import re

import requests
import shapefile

shapefile_path = 'shapefiles/PL1992_5000_025'
wms_request_base_url = 'https://mapy.geoportal.gov.pl/wss/service/PZGIK/ORTO/WMS/SkorowidzeWgAktualnosci'
wms_request_parameters = {
    'SERVICE': 'WMS',
    'request': 'GetFeatureInfo',
    'version': '1.3.0',
    'layers': 'SkorowidzeOrtofotomapyStarsze,SkorowidzeOrtofotomapy2018,SkorowidzeOrtofotomapy2019,SkorowidzeOrtofotomapy2020',
    'styles': '',
    'crs': 'EPSG:2180',
    'width': 1000,
    'height': 1000,
    'format': 'image/png',
    'transparent': 'true',
    'query_layers': 'SkorowidzeOrtofotomapyStarsze,SkorowidzeOrtofotomapy2018,SkorowidzeOrtofotomapy2019,SkorowidzeOrtofotomapy2020',
    'i': 1,
    'j': 1,
    'INFO_FORMAT': 'text/html'
}
regexp = '{url:"(.+)",godlo:"(.+)", aktualnosc:"(.+)", wielkoscPiksela:"(.+)", ukladWspolrzednych:"(.+)", calyArkuszWyeplnionyTrescia:"(.+)", modulArchiwizacji:"(.+)", zrodloDanych:"(.+)", kolor:"(.+)", numerZgloszeniaPracy:"(.+)", aktualnoscRok:"(.+)"}'
compiled_regexp = re.compile(regexp, re.IGNORECASE)
attribute_names = [
    'url',
    'godlo',
    'aktualnosc',
    'wielkoscPiksela',
    'ukladWspolrzednych',
    'calyArkuszWyeplnionyTrescia',
    'modulArchiwizacji',
    'zrodloDanych',
    'kolor',
    'numerZgloszeniaPracy',
    'aktualnoscRok'
]

if __name__ == '__main__':

    # open shape file
    with shapefile.Reader(shapefile_path) as shp:
        total_rows = shp.numRecords

        # for each record in the shapefile (we are ignoring geometry here)
        for row_num, row in enumerate(shp.iterRecords()):
            # -----------------------------------------------------
            # remove before production use:
            if row_num > 1: break  # for testing, get first n rows
            # -----------------------------------------------------

            # get row as dict so we have attribute names
            r = row.as_dict()

            # create bounding box
            minx = min(r['x1'], r['x2'], r['x3'], r['x4'])
            miny = min(r['y1'], r['y2'], r['y3'], r['y4'])
            maxx = max(r['x1'], r['x2'], r['x3'], r['x4'])
            maxy = max(r['y1'], r['y2'], r['y3'], r['y4'])
            bbox = [minx, miny, maxx, maxy]
            # create comma separated string representation of bbox that we will use for the HTTP request to WMS service
            bbox_str = ','.join([str(el) for el in bbox])
            wms_request_parameters['bbox'] = bbox_str

            print(f'{str(row_num+1).zfill(6)}/{total_rows}', '- godło:', r['godlo'], '- bbox:', bbox)

            # make the HTTP request
            response = requests.get(wms_request_base_url, params=wms_request_parameters)
            if response.ok:
                # construct dictionary with data parsed from http response using regular expressions
                for matches in compiled_regexp.findall(response.text):
                    record = {name: value for name, value in zip(attribute_names, matches)}
                    print(record)

I w jaki sposób teraz to zaimplementować w kodzie?

Nie bardzo wiem co chcesz osiągnąć i z czym masz problem.
Te PDFy to łatwiej chyba wyklikać ze strony, bo nie bardzo widzę scenariusza gdzie możesz chcieć pobierać ich dużo na raz, ale jeżeli tak to może być łatwiej tą wtyczką: https://plugins.qgis.org/plugins/pobieracz_danych_gugik/.
Do wizualizacji samych danych wektorowych BDOT10k w tym stylu GUGIK udostępniał wtyczkę do QGIS: https://plugins.qgis.org/plugins/BDOT10k_GML_SHP_Loader/
A jeżeli chodzi o te identyfikację arkuszy to chyba były pliki z sekcji II: https://www.gov.pl/web/gugik/dane-udostepniane-bez-platnie-do-pobrania-z-serwisu-wwwgeoportalgovpl

Dzięki za wszelkie informacje. Tak, chodzi o pobieranie dużej ilości na raz, chce napisać sobie program python do tego. Qgis też to oczywiście robi, ale potrzebuje robić kilka czynności konkretnych na raz a nie bawić się w qgis :) dzięki za info, ogarnę to wszystko dzisiaj

Mam problem, żeby stworzyć taki program. Udało mi się stworzyć program python który pobiera mapę dla zlokalizowanej działki po numerze id działki przy pomocy ULDK gugik i z wykorzystaniem WMS https://integracja02 gugik.gov.pl/cgi-bin/KrajowaIntegracjaEwidencjiGruntow. Nie potrafię jednak zrobić programu który pobierze jakis podkład pod tą ewidencję działek żeby nałożyć je na siebie i dostać mape na której jest działka oraz podkład typu mapa topograficzna BDOT10k lub G2MOBILE500

@Gecko171 zapomniałem odpowiedzieć 😅 Podkład w sensie chcesz skorzystać z WMS/WMTS to co jest na geoportalu żeby u siebie wyświetlić? Chcesz wygenerować obrazek JPG czy coś z tym kawałkiem mapy plus działką czy coś? Czy coś innego?

Tak. Mam już taki fajny program w python określa współrzędne wierzchołków działek i jej centrum i chciałbym dopisać, żeby pobierało podkład w jpg (mapa topo w stylu BDOT10k lub G2mobile500 mając bbox Tylko nie wiem jak stworzyć zapytanie dla tych wms

Adresy usług są na stronie geoportalu: https://www.geoportal.gov.pl/uslugi/usluga-przegladania-wms
Ogólny opis jak się tworzy zapytania do WMS jest całkiem spoko na wiki: https://pl.wikipedia.org/wiki/Web_Map_Service
Najłatwiej to dodać warstwę w QGIS, odpalić panel z podglądem jakie requesty są wysyłane (f12) i skopiować.

Ok, już widzę. A w jaki sposób sprawdzić i określić jak są ustalane Tilerow i Tilecol? Czy jest sposób na przeliczenie tego z bbox?

Tile row i col to dla WMTS. To trzeba by poczytać jak to się wylicza, ale większość warstw masz też w WMS, w którym wystarczy ci bbox.