Ekstrakcja rowów za pomocą LIDAR DTM (rozdzielczość 1m)

9

Próbuję wyodrębnić rowy za pomocą LIDAR DTM za pomocą oprogramowania eCognition, ale tak naprawdę mam trudność od czasu, gdy oparłem swój przepływ pracy zatytułowany „Ekstrakcja sieci rowów i charakterystyka hydrogeomorfologiczna przy użyciu DTM pochodzącego z LiDAR na terenach podmokłych” t mają szczegółowe zestawy reguł wprowadzone w eCognition.

Chcę tylko zapytać, kto z was tutaj ma pomysł, jak wydobywać rowy w krajobrazie rolniczym, a może w otoczeniu miejskim, na przykład przy drogach?

Lub jeśli kiedykolwiek nie wykonałeś takiej pracy lub projektu, kto inny może mieć świetny pomysł, który pozwoli mi wydobyć rowy.

Dróżnik
źródło

Odpowiedzi:

9

Wcześniej musiałem mapować rowy z 1 m DEM pochodzących z krajobrazu rolniczego DEM. Opracowanie odpowiedniego przepływu pracy z pewnością jest trudnym zadaniem. Twoja zdolność do skutecznego wyodrębnienia sieci rowów będzie zależeć od wielu czynników. Na przykład, czy interesują Cię tylko rowy przydrożne? Jeśli tak, to czy drogi znajdują się na nasypach (jak to zwykle bywa w warunkach rolniczych) i czy masz dokładny plik wektorowy drogi (może to mieć kluczowe znaczenie dla tego zadania)? A może interesują Cię również rowy melioracyjne? Czy masz surową chmurę punktów LiDAR, czy tylko interpolowaną DEM? Jeśli masz później, bardzo ważne jest, aby nie stosować żadnego filtru wygładzającego, co niestety jest zwykle wykonywane z LiDAR DEM z powodu ich nadmiernej chropowatości powierzchni przed analizą. Niestety, nisko położony przydrożny rów można łatwo usunąć za pomocą filtra średniego 3 x 3 w 1 m DEM. Jeśli masz chmurę punktów LiDAR, zaleciłbym interpolację siatki DEM za pomocą schematu interpolacji Nearest Neighbor (przy założeniu wysokiej gęstości punktów), ponieważ chociaż zwiększy to chropowatość powierzchni w porównaniu z czymś takim jak IDW, lepiej zachowa rowy .

Teraz, zakładając, że masz wektor drogi i że jesteś zainteresowany tylko rowami przydrożnymi, oto przepływ pracy, którego możesz użyć:

  1. Wykonaj filtr różnicy wysokości od średniej. W otwartych źródłach GIS Whitebox Geospatial Analysis Tools , których użyłem do tego zadania, istnieje dosłownie narzędzie o nazwie „Różnica od średniej wysokości” (DFME), które jest idealne do tego przepływu pracy. Jeśli jednak z jakiegoś powodu nie możesz korzystać z Whitebox, wykonaj tradycyjny filtr średni (filtr skrzynkowy) i odejmij wynik od oryginalnego DEM (można również zastosować filtr „górnoprzepustowy”). Być może będziesz musiał poeksperymentować z wielkością filtra, która będzie zależeć od szerokości twoich cech rowu, ale powinna być nieco szersza niż rowy. Dla moich danych ustawiam parametr „Search Neighborhood Size” narzędzia DFME na 5 komórek, co tworzy filtr 11 x 11.

  2. Musisz wyzerować raster DFME, aby wyodrębnić wszystkie komórki siatki o „niskiej” wartości DFME. Znowu będzie to zależeć od twoich danych, a zwłaszcza od głębokości rowów w twoim DEM. Po prostu użyłem do tego kalkulatora rastrowego Whitebox, z wyrażeniem [niskie komórki] = [DFME] <(- 0,15). Jednostki tego parametru „0,15” w progu są takie same jak jednostki Z twojego DEM. To skutecznie mówi, daj mi wszystkie komórki siatki, które leżą co najmniej 15 cm (mój DEM jest w metrach) poniżej ich otoczenia.

  3. Buforuj wektor drogi, tak aby tworzył maskę drogi, która jest wystarczająco szeroka, aby obejmowała drogę i rowy przydrożne. Będzie to zależeć od szerokości twojej drogi. Jeśli masz rozległy teren, może istnieć wiele szerokości dróg, np. Główne drogi są zwykle szersze niż drogi drugorzędne. W moim przypadku bufor drogowy o długości 10 m działał dobrze.

  4. Konwertuj ten wielokąt bufora drogowego na raster, używając DFME lub DEM jako obrazu podstawowego, z którego raster wyjściowy wyprowadzi swoją rozdzielczość i zasięg. Możesz wcześniej przyciąć bufor drogowy do śladu warstwy siatki rastrowej, aby przyspieszyć ten proces, jeśli masz bardziej rozległą wektorową sieć dróg niż strona DEM, co miało miejsce w moim przypadku. W zależności od tego, jak działa konwersja wektora do rastra, być może zmieniono wartość w buforze drogowym raster na 1 dla dróg i 0 dla wszystkiego innego. Ponownie kalkulator rastrowy może być do tego przydatny.

  5. Pomnóż swój końcowy raster buforu drogowego przez swój progowy obraz DFME.

  6. Jeśli jesteś naprawdę elegancki, możesz zastosować algorytm przerzedzania linii do rastra wynikający z kroku 5, aby stworzyć ładną sieć cienkich linii rowów przydrożnych.

Na poniższym zdjęciu sieć rowów przydrożnych jest pokazana na czarno, nałożona na obraz DFME, a odcień DEM jest przezroczysty. Myślę, że w tym przypadku działało to całkiem dobrze, ale znowu wymaga nieco finezji i zabawy z różnymi parametrami.

wprowadź opis zdjęcia tutaj

Jeśli okaże się, że nie interesują Cię tylko rowy przydrożne, to w Whitebox znajduje się narzędzie o nazwie Map Gully Depth, które może być wykorzystane do tego zadania w zależności od twoich danych i krajobrazu. Jest trudny w użyciu, więc daj mi e-mail, jeśli zdecydujesz się wybrać tę trasę, a ja chętnie udzielę wskazówek. Głównym problemem jest to, że trudno jest zastosować wzorce przepływu powierzchniowego (np. Obrazy akumulacji przepływu) do mapowania rowów, ponieważ w przeciwieństwie do strumieni, rowy w warunkach rolniczych są wykorzystywane zarówno do tymczasowego magazynowania wody, jak i do transportu wody. Na ogół nie mają gradientów, które można znaleźć w naturalnych strumieniach. Niemniej jednak opracowałem algorytm naruszenia depresji w Whitebox, który może być użyty do poprawy drenażu przez rowy, co może być również przydatne do mapowania rowów w polu.

WhiteboxDev
źródło
Bardzo dziękuję za informacje, które udostępniłeś na temat wydobycia rowu @WhiteboxDev. Będę bardzo szczęśliwy, mogąc być z tobą połączony, ponieważ uważam, że masz wszystkie niezbędne pomysły, które są mi potrzebne do realizacji mojego celu: wydobycia rowów. Właściwie interesują mnie zarówno rowy przydrożne, jak i rowy widoczne w krajobrazach rolniczych, które tworzą sieci. Spróbuję wyodrębnić przydrożne rowy przy użyciu przepływu pracy, o którym wspomniałeś powyżej, ale nie mam jeszcze oprogramowania Whitebox GAT. Czy to za darmo? Gdzie mogę pobrać instalator? Tymczasem może ArcGIS zrobi dla
1
Witaj @DitchExtractor tak, Whitebox to oprogramowanie typu open source i bezpłatne, działające w systemach Windows, Mac i Linux. Możesz go uzyskać, klikając link w mojej odpowiedzi, który przeniesie Cię na stronę główną projektu. Następnie wybierz pobierz. Są tam również moje dane kontaktowe, jeśli chcesz się ze mną skontaktować bezpośrednio.
WhiteboxDev,