Mam plik rastrowy dotyczący użytkowania gruntów w północnych Włoszech (uprawy / obszary miejskie itp.)
Zdefiniowałem niektóre strefy buforowe wokół niektórych punktów warstwy wektorowej i chciałbym obliczyć% tych obszarów zajmowanych przez każdą wartość pliku rastrowego, aby oszacować średnie wykorzystanie wybranych obszarów (np. 70% upraw, 30% miejskich itp.).
Wiem, że w ArcGIS jest funkcja Analityka przestrzenna, która wykonuje tak zwany „ Obszar tabelaryczny ”, ale nie mogę znaleźć takiej funkcji w QGIS.
Jeśli ktoś mógłby mi pomóc, byłbym bardzo wdzięczny.
Odpowiedzi:
spróbuj użyć „Cross-Classification and Tabulation” SAGA GIS z wtyczki Sextante.
Mam nadzieję, że to pomoże
źródło
Kodowałem wtyczkę QGIS o nazwie LecoS (spójrz w narzędziu do pobierania wtyczek, wtyczkach eksperymentalnych, najpierw zainstaluj Scipy!). Może osiągnąć to, co chcesz (procent pokrycia terenu dla danego rastra), jednak do tej pory nie znalazłem czasu na kodowanie nakładki wielokątów w pythonie, dlatego musisz wcześniej przyciąć raster indywidualnie dla każdego bufora i funkcji wektorowej. Zobacz przypadek użycia tutaj .
źródło
W tym przypadku można zastosować klasyfikację krzyżową i tabulację, jak wspomniano powyżej. Aby go użyć, najpierw musisz mieć warstwę rastrową z wartościami int. W zależności od twojej warstwy użytkowania, będziesz musiał wykonać pewne przeklasyfikowanie lub przekodować, aby przygotować swój raster.
Innym założeniem jest to, że twoje wielokąty (strefy buforowe) muszą mieć pewną wartość int. Jeśli nie, musisz go utworzyć (z formułą kalkulatora $ id + 1).
Rasteryzuj warstwę wektorową, używając kodu int jako wartości rastrowej. Tutaj musisz podać pewne warunki: użyj tej samej wartości w pikselach i tego samego rozszerzenia, co raster, z którego chcesz wyodrębnić informacje.
Teraz skorzystaj z narzędzi Cross-Classification and Tabulation z narzędzi SAGA. Musisz wziąć pod uwagę, że utworzy to tabelę xy z liczbą pikseli każdej kategorii użycia w każdym wielokącie. Tutaj musisz wstawić „maksymalną liczbę klas”, co oznacza, że na przykład, jeśli masz raster z 10 kategoriami i 25 klasami wielokątów, twoje maksimum wyniesie 25.
Oto wynik końcowy: tabela, w której masz liczbę pikseli każdej kategorii w każdym wielokącie. Jeśli chcesz przekonwertować go na obszar, pomnóż liczbę komórek przez rozdzielczość rastra.
źródło
Narzędzie Obszary tabelaryczne znajduje się w przyborniku Spatial Analyst w sekcji Narzędzia strefowe dla ArcGIS 10. Nie jesteś pewien, jakiej wersji używasz. Pamiętaj, że musisz mieć osobną licencję na rozszerzenie Spatial Analyst, aby korzystać z tych narzędzi. Wygląda na to, że to narzędzie zrobi to, czego szukasz.
Mam nadzieję że to pomoże
źródło