Jestem względnie początkującym do korzystania z GIS i korzystam z QGIS 2.0.1 w Linuksie. Mam dwie historyczne mapy, które chcę przeanalizować, które pokazują wzorce użytkowania gruntów na tym samym obszarze w dwóch różnych momentach w czasie. Mam je zeskanowane i georeferencyjne jako warstwy w jednym pliku. Obok siebie wyglądają tak:
Najważniejsze, co mnie interesuje, to porównanie zasięgu jasnych i ciemnozielonych obszarów między dwiema mapami. Czy jest to możliwe, a jeśli tak, jakie jest najprostsze podejście? Czy można to zrobić na podstawie analizy rastrowej? A jeśli muszę zrobić plik kształtu, jaki jest najlepszy sposób?
Co już rozważałem:
Ręczne rysowanie plików kształtów jako wielokątów, jak opisano w tym samouczku . To byłoby BARDZO nudne.
Tworzenie uproszczonych obrazów rastrowych o wysokim kontraście za pomocą wyboru kolorów, filtrów itp. Metodą prób i błędów w Gimpie i przekształcanie go w plik kształtu. Wyniki były bardzo niechlujne.
Odpowiedzi:
Posteryzacja była świetnym początkiem: wyeliminowała większość artefaktów kompresji i uprościła kartografię na tyle, aby umożliwić dodatkowe czyszczenie.
Znaczna część czyszczenia kategorycznego rastra obejmuje tak zwane operacje „morfologiczne”. Obejmują one rozszerzenie jednej kategorii na sąsiadów, ponowne jej zmniejszenie i ponowne grupowanie przyległych komórek monokategoryjnych w ich własne kategorie.
Zwykle potrzebne są pewne eksperymenty, choćby dlatego, że artefakty, które mają zostać usunięte - litery, linie kreskowania itp. - będą różnić się rozmiarem pikseli między skanami. Aby rozpocząć, zilustruję, co te procedury mogą osiągnąć na przykładzie.
Oryginał po posteryzacji wygląda następująco. To siatka z trzema kategoriami pokazanymi w trzech kolorach. Naszym celem jest stworzenie siatki, w której ciemnozielone obszary są przekształcane w ciągłe fragmenty, bez nadmiaru, kropek lub nieistotnej linii, nadające się do późniejszej analizy za pomocą algebry rastrowej.
Po rozwinięciu ciemnozielonych obszarów o zaledwie jeden piksel do wszystkich otaczających obszarów ten obraz daje:
(Aby uzyskać bardziej precyzyjną kontrolę, możesz ograniczyć ekspansję tylko do czarnych obszarów, jeśli pozwala na to GIS).
Aby wyeliminować wiele cienkich linii zielonych artefaktów i małych wysepek, zmniejszmy zieleń z powrotem o dwa piksele
a następnie, aby zrównoważyć wszystkie rozwijające się i kurczące (aby zmniejszyć obciążenie), rozszerzymy go o jeszcze jeden piksel:
Grupowanie regionów identyfikuje następujące ciągłe obszary zieleni:
Każda łatka jest pokazana w innym kolorze.
Użyj operacji warunkowej lub operacji SetNull, aby wyeliminować małe poprawki. Jak malutki Przejrzałem tabelę atrybutów i odkryłem, że wiele łat zajmowało od 6 do 47 komórek; potem nastąpił skok do 422 komórek. Wybrałem próg w ramach tego skoku (100) i usunąłem wszystkie komórki z liczbą (nie wartościami!) Mniejszą niż ten próg. Oto, co pozostało, nałożone na oryginał do porównania:
Osiągnęliśmy dość precyzyjną reprezentację obszarów zainteresowania, odpowiednich do wykrywania i kwantyfikacji zmian w stosunku do podobnie przetwarzanych obrazów. Wziąłem trochę pracy, ale jest to o wiele mniej pracy niż ręczne digitalizowanie oryginalnego skanu i - pod warunkiem, że skany są wykonywane w stałej rozdzielczości - mogą być półautomatyczne. (Ponieważ oryginalne mapy używają różnych kolorów, na początku musi nastąpić inteligentna interwencja, aby wybrać odpowiednie kolory do powiększania i zmniejszania.) Każdy z kroków jest również dość szybkim obliczeniem, więc prawdopodobnie możesz sobie pozwolić na zeskanowanie oryginału mapy w bardzo wysokich rozdzielczościach dla największej precyzji.
źródło
Uzyskanie przybliżonej warstwy rastrowej w Gimpie i przekonwertowanie jej na wektor w QGIS prawdopodobnie zaoszczędziło mi trochę czasu, ale wygląda na to, że nie ma sposobu, aby uniknąć godzin czyszczenia wynikowych plików kształtów, wierzchołek po wierzchołku.
źródło
Ok, może to zadziała, może nie. w zależności od jakości skanu. możesz ustawić przezroczystość określonego koloru na wartość procentową lub użyć narzędzia wartości, aby wyizolować żądany kolor.
Nie zamierzam tego przypisywać, ponieważ zadałem pytanie, zanim ... mój rzeczywiście chciał wybrać domy z otwartym mapowaniem widoku ulicy. Daj mi więc znać, jeśli to pomoże.
Zidentyfikuj wielokąty na obrazie rastrowym
Daj mi znać, jeśli to pomoże .... mogę usunąć, jeśli nie jest to w porządku.
źródło
W Gimp masz kilka narzędzi wyboru, które mogą ułatwić ci pracę.
Myślę o niewyraźnym narzędziu do zaznaczania (wybieranie według obszaru koloru), które można ustawić tak, aby było mniej lub bardziej wrażliwe na zmiany kolorów (przy użyciu pierwszego przesłanego obrazu uzyskałem dobre wyniki z wartością progową 13,0). W ten sposób uzyskuję całkiem czysty wynik bez dużej ilości pasożytów, które następnie można łatwo i szybko usunąć wybierając narzędzie do wycierania
Po wybraniu możesz zmienić kolor tych obszarów za pomocą kontrastowych kolorów, ponownie zaimportować do GIS, a następnie wektoryzować?
Ale to nie rozwiąże problemu wskazanego przez gawędziarza o wykluciu terenu, ale skoro nie jest ich wiele, być może możesz je wektoryzować od zera bez poświęcania na to dużo czasu?
źródło