Od danych punktowych przez interpolację do danych dla obszarów

Dzień dobry! Mam dosyć złożony problem i niezbyt mogę znaleźć rozwiązanko w Internecie (zwykle szukam w j. angielskim przez Google i się udaje).

Przykładowo mam dane o zanieczyszczeniach powietrza mierzone w stacjach pomiarowych (punkty). Mam również podział administracyjny (powiaty albo podregiony). Chciałbym wykonać interpolację na podstawie danych ze stacji, a następnie do każdej z jednostek przestrzennych (załóżmy do każdego powiatu) przypisać wartość średnią na podstawie przeprowadzonej interpolacji. Zaznaczyłem na czerwono to, czego kompletnie nie umiem zrobić.

Dysponuję QGIS-em i Surferem, ale jeżeli ktoś zna sposób w innym softwarze, mogę spróbować i w innym. Może znacie jakiś tutorialik?