Procent wielokąta w jednym pliku kształtu w obrębie wielokąta innego

Jestem nowicjuszem, przepraszam, jeśli to oczywiste / został już zapytany i odpowiedział, ale nic nie mogłem znaleźć.

Mam dwa kształty: 1. administracyjna warstwa graniczna dla hrabstwa w Wielkiej Brytanii znana jako granica LSOA, która ma 500 małych stref 2. strefa powodziowa.

Idealnie chcę dowiedzieć się, które z małych stref LSOA znajdują się w ≥50% w strefie powodziowej i kończą na tak / nie lub 1/0 dla każdej z 500 stref LSOA.

Ale nie wiem jak to zrobić. Pomyślałem, że mogę dołączyć do dwóch plików kształtów, ale nie ma między nimi wspólnej cechy. Potem pomyślałem, że mogę użyć funkcji Dołącz atrybut według lokalizacji, która działała i pokazuje mi, które LSOA znajdują się w strefie powodzi, ale to prawie wszystkie z nich (patrz zdjęcie 2).

Myślę, że to jest problem SQL, ale nie wiem. Jestem nowy w QGIS i nigdy nie korzystałem z PostgreSQL.

Każda pomoc byłaby bardzo mile widziana. Mogę podać wszelkie informacje, których potrzebujesz, aby mi pomóc.

Jest to stosunkowo proste zadanie przy użyciu narzędzi geoprzetwarzania zawartych w QGIS.

1. Oblicz obszar swoich stref LSOA.

   • Otwórz tabelę atrybutów warstwy LSOA.
   • Włącz tryb edycji.
   • Otwórz kalkulator pola.
   • Utwórz nowe pole typu „Liczba dziesiętna (rzeczywista)” z wyrażeniem „$ area”.
   • Wyłącz tryb edycji (zapisywanie zmian).

2. Scal warstwę strefy zalewowej w jeden element wieloczęściowy.

   • Vector > Geometry Tools > Singleparts to Multipart.
   • Wybierz „--- Scal wszystko ---” w polu Unikatowy identyfikator.

3. Przecinaj warstwę strefy LSOA z wieloczęściową warstwą strefy zalania.

   • Vector > Geoprocessing Tools > Intersect.
   • Warstwa wejściowa to strefy LSOA, warstwa przecinająca to strefy powodziowe.

4. Powstała warstwa będzie częściami stref LSOA (z atrybutami z warstwy stref LSOA), które nakładają się na warstwę stref powodziowych. Aby obliczyć udział każdej strefy LSOA w strefie powodziowej:

   • Następnie oblicz obszar przeciętych elementów (jak w kroku 1)
   • Dodaj kolejne pole, dzieląc oryginalny (całkowity) obszar przez przecięty obszar. Wynikiem jest liczba dziesiętna od 0 do 1. Pomnóż przez 100, aby podać wartość procentową.

5. Połącz oryginalną warstwę LSOA z warstwą przecinającą się, używając unikalnego identyfikatora współdzielonego przez obie warstwy.

6. Wyeksportuj połączoną warstwę jako nowy plik kształtu.

7. Usuń zduplikowane atrybuty.

Zrobione!

Bez kroku 2 dla każdej funkcji LSOA utworzono by oddzielną funkcję. Prawdopodobnie nie jest to, czego chcesz, jeśli interesuje Cię tylko całkowity zasięg dla każdej strefy LSOA. Jeśli chcesz rozróżnić zalew fluwialny / pływowy / pluwialny (a dane strefy powodziowej to obsługują), możesz przekonwertować pojedyncze elementy na wieloczęściowe, określając pole „TYPE” jako pole Unique ID.

Możesz użyć funkcji spatialite i niektórych przestrzennych funkcji SQL.

Select t1.geometry, t1.ID, area(t1.geometry), area(t2.geometry) ...... (anything you need to have in the table results)

(area(intersection(t1.geometry,t2.geometry))) as "Commun_AREA"

, ("Commun_AREA"*100/(area(t1.geometry))) as "Percent_AREA"

From lsoa as t1, flood_zone as t2

Where Intersects( t1.geometry,t2.geometry ) = 1

Wydaje się, że można to zrobić o wiele łatwiej niż odpowiedzi. Osobiście użyłbym prostego skryptu python:

floodName = "the layer name here"
boundryName = "the layer name here"
fieldName = "the name of the field to contain the output 1/0"
minCoverage = 0.5 # the minimum amount of area covered to write 1
updateMap = [] # this will store values to be written    

# get layers
floodLayer = QgsMapLayerRegistry.instance().mapLayersByName(floodName)[0]
boundryLayer = QgsMapLayerRegistry.instance().mapLayersByName(boundryName)[0]
fieldIndex = boundryLayer.dataProvider().fieldNameIndex(fieldName)    

# iterate through boundries
for b in boundryLayer.getFeatures():
    # get only flood features that intersect with this feature's bounding box
    # this will make the script go way faster than it would otherwise
    request = QgsFeatureRequest().setFilterRect(b.geometry().boundingBox())
    floodGeom = geometry()
    floodFeat = QgsFeature()
    iter = floodLayer.getFeatures(request)
    iter.nextFeature(feat)
    while iter.nextFeature(feat):
        floodGeom = floodGeom.combine(feat.geometry())
    intersectGeom = b.geometry().intersection(feat.geometry())
    if intersectGeom.area() > minCoverage * b.geometry().area():
        updateMap[b.id()] = {fieldIndex : 1}
    else:
        updateMap[b.id()] = {fieldIndex : 0}

boundryLayer.dataProvider().changeAttributeValues(updateMap)

ocenia to tylko wielokąty zalewowe, które przecinają się z obwiednią każdej warstwy granicznej, więc uruchomienie powinno być dość szybkie, a następnie aktualizuje tylko jedno pole w istniejącej warstwie (zamiast złożonej operacji tworzenia całej nowej warstwy i kopiowania starych wartości następnie usuwanie)

Miałem ten sam problem, co KJ, postępując zgodnie z instrukcjami Snorfalorpagusa przy użyciu metody „Przecięcie” w kroku 3. Obliczenie zajęło trochę czasu, a to, co mi pozostało, było puste.

Próbowałem wykonać te same kroki, z wyjątkiem użycia metody „Klip” w QGIS zamiast przecinania - tak więc w twoim przykładzie pozostałyby części obszarów, które NIE są objęte strefą powodzi. Wydawało się, że to działa z jakiegoś powodu i mogłem skorzystać z obliczenia pola „Obszar” z poprzedniego kroku, a także nowego obliczenia „Obszar” dla pozostałych części każdego wielokąta, aby obliczyć% każdego obszaru, który NIE był pokryty drugą warstwą wielokąta.

Technicznie jest to odwrotność tego, o co prosiłeś. Ale stamtąd to tylko kwestia odjęcie każdą wartość od 1 do dostać to, co jest objęte strefą przeciwpowodziowej.