Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- # Konwertujemy m^2 do km^2
- def convert_m2_to_km2(liczba):
- return liczba/10**6
- # Szukamy listę województw oraz listę powierzchni
- def get_data(vlayer, nr=6):
- names = []
- areas = []
- d = qgis.core.QgsDistanceArea()
- d.setEllipsoid('GRS80')
- d.setEllipsoidalMode(True)
- # Dostawca danych:
- provider = vlayer.dataProvider()
- # Pobranie wszystkich obiektów:
- feature_list = provider.getFeatures()
- # Potrzebujemy wartości True, aby pętla się wykonała:
- feat = qgis.core.QgsFeature()
- # Iteracyjne pobieranie informacji o poszczególnych elementach (według konfiguracji dostarczonej przez funkcję select)
- while feature_list.nextFeature(feat):
- # Pobranie listy atrybutów danego elementu:
- attrs = feat.attributes()
- names.append(attrs[nr].encode('utf-8'))
- # Otrzymujemy powierzchnie w m^2, żeby było w km^2 konwertujemy
- area = convert_m2_to_km2(d.measure(feat.geometry()))
- areas.append(area)
- return names, areas
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement