Zastępowanie wartości danych NoData w rastrze z wykorzystaniem informacji statystycznych obliczonych na podstawie wartości pikseli otaczających


Tomasz Bartuś



Problem

Mamy kilka plików ASCII grid NMT (w formacie ARC/INFO ASCII GRID). Wszystkie z tym samym 32-bitowym typem piksela z wartością zmiennoprzecinkową i rozdzielczością 5 m. Chcemy połączyć te pliki w jeden plik o rozdzielczości 5 m. Tworzymy plik mozaiki. Niestety na krawędziach niektórych elementów mozaiki znajdują się luki z pikselami o wartościach NoData.

Raster NoData Gaps
Fig. 1. Połączone rastry z widocznymi lukami z pikselami o wartościach NoData

Rozwiązanie

Rozwiązanie polega na zamianie wartości pikseli NoData wartościami średniej arytmetycznej obliczonymi na podstawie najbliższego sąsiedztwa pikseli luk. W tym przypadku wartości pikseli zostaną obliczone na podstawie średniej wyciągnietej z matrycy 5 × 5 pikseli. W tym celu wybierzmy:

ArcToolbox > Spatial Analyst Tools > Map Algebra > Raster Calculator

W kalkulatorze rastrów wprowadź wyrażenie:

Con(IsNull("raster"), FocalStatistics("raster", NbrRectangle(5,5, "CELL"), "MEAN"), "raster")

Uwaga!: słowo raster zamień nazwą zbioru rastrowego z tabeli zawartości

Źródła

How To: Remove and replace no data values within a raster using statistical information from the surrounding data values.

 
 

Informacje wstępne

Oprogramowanie, licencja, instalacja...

Dane

ArcGIS Desktop: 23MB (73MB)
ArcGIS Pro: 94MB (1,07GB)
 
Szczegółowy program zajęć

1. Wstęp do GIS

Umiejętności
Wykład 1: Wstęp do GIS
 
Ćwiczenie 1: Wycieczka po San Diego. Wstęp do ArcGIS
 
Wykład 2: Aplikacja ArcGIS
Ćwiczenie 2: Szacowanie szkód wywołanych przez tornado. Od metainformacji do wyników analizy
 
Wykład 3: Rozwiązywanie zadań przy użyciu GIS
Ćwiczenie 3: Wybór lokalizacji dla centrum młodzieżowego. Zapytania atrybutowe i przestrzenne
 
Test

2. Symbolizacja map

Umiejętności
Wykład 4: Symbole i adnotacje
Ćwiczenie 4: Lokalizacja ośrodka dziennego dla seniorów. Symbolizacja i etykietowanie obiektów na mapach
 
Wykład 5: Symbole bazujące na atrybutach
Ćwiczenie 5: Lokalizacja ośrodka dziennego dla seniorów. Symbole oparte na atrybutach
 
Wykład 6: Metody klasyfikacji danych
Ćwiczenie 6: Lokalizacja ośrodka dziennego dla seniorów. Klasyfikacja danych
 
Wykład 7: Mapy gęstości i proporcji
Ćwiczenie 7: Lokalizacja ośrodka dziennego dla seniorów. Mapy gęstości i normalizacja danych
 
Test

3. Układy współrzędnych i odwzorowania kartograficzne

Umiejętności
Wykład 8: Układy współrzędnych geograficznych
Ćwiczenie 8: Położenie kabla transatlantyckiego. Modyfikacje układów współrzędnych
 
Wykład 9: Odwzorowania kartograficzne
Ćwiczenie 9: Położenie kabla transatlantyckiego. Właściwości odwzorowań kartograficznych
 
Test

4. Organizacja danych geograficznych

Umiejętności
Wykład 10: Modele danych geograficznych
Ćwiczenie 10: Ocena zagrożenia powodziowego. Analiza danych geograficznych
 
Wykład 11: Geobazy
Ćwiczenie 11: Ocena zagrożenia powodziowego. Tworzenie przestrzennej bazy danych
 
Ćwiczenie 12: Salzburg. Geobaza wspierająca branżę turystyczną
 
 
Test

5. Tworzenie i edycja danych

Umiejętności
Wykład 12: Tworzenie i edycja danych GIS
Ćwiczenie 13: Centrum handlowe Galeria. Modyfikacja śladów budynków
 
Wykład 13: Topologia
Ćwiczenie 14: Centrum handlowe Galeria. Modyfikacja obiektów z zachowaniem topologii
 
Wykład 14: Edycja atrybutów obiektów
 
Ćwiczenie 15: Centrum handlowe Galeria. Modyfikacja atrybutów
 
Wykład 15: Tworzenie nowych obiektów
 
Ćwiczenie 16: Park miejski. Tworzenie nowych klas obiektów i ich atrybutów
 
Test

6. Analizy GIS

Umiejętności
Wykład 16: Procedura analiz GIS
 
Ćwiczenie 17: Poszukiwania nieruchomości spełniającej kryteria. Definiowanie problemu i wybór danych do analiz
 
Wykład 17: Zapytania atrybutowe i przestrzenne, łączenie tabel
 
Ćwiczenie 18: Poszukiwania nieruchomości spełniającej kryteria. Analiza GIS
 
Wykład 18: Prezentacja wyników analiz
 
Ćwiczenie 19: Poszukiwania nieruchomości spełniającej kryteria. Przygotowanie mapy prezentacyjnej
 
Ćwiczenie 20: Poszukiwania nieruchomości spełniającej kryteria. Tworzenie raportów
 
Test

7. Geoprzetwarzanie i modelowanie

Umiejętności
Wykład 19: Geoprzetwarzanie
 
Ćwiczenie 21: Ocena zniszczeń pożarowych. Geoprzetwarzanie
 
Wykład 20: Modele i modelowanie
 
Ćwiczenie 22: Ocena zniszczeń pożarowych. Tworzenie i obsługa prostych modeli
 
Ćwiczenie 23: Przetarg na zakup drewna. Praca z istniejącymi modelami
 
Test

8. Tworzenie profesjonalnych map

Umiejętności
Wykład 21: Elementy mapy, praca w widoku układu
 
Ćwiczenie 24: Przetarg na zakup drewna. Praca z układami
 
Wykład 22: Tworzenie profesjonalnych map
 
Ćwiczenie 25: Siedliska nietoperzy w południowo zachodniej części Ameryki Północnej. Tworzenie profesjonalnych map"
 
Test

Tips and Tricks

Geomedia Professional

  1. Konwersja obiektów z klasy złożonej (Compound) do klas o geometrii punktowej, liniowej i powierzchniowej
  2. Wyszukiwanie duplikatów obiektów
  3. Łączenie obiektów w obrębie klasy
  4. Wycinanie fragmentu klasy (mapy, sceny)
  5. Obliczanie parametrów geometrycznych (np. łącznej powierzchni poligonów)
  6. Masowa modyfikacja atrybutów
  7. Liczba obiektów punktowych wewnątrz oczek siatki
  8. Liczba platów poligonowych wewnątrz oczek siatki
  9. Centroidy obiektów poligonowych
  10. Liczba kategorii płatów i wsp. SHDI
  11. Łączna długość linii wewnątrz ob. poligonowych

Geomedia Grid

  1. Import danych .xyz do GG
  2. Rasteryzacja
  3. Esport GG do ArcGis
  4. średnia wysokość położenia poligonów

G.S. Surfer

  1. Konwersja współrzędnych strony na współrzędne mapy
  2. Surfer Variogram Tutoria
  3. Anizotropia

ArcGIS Desktop

  1. Autokorelacja przestrzenna
  2. MAUP
  3. Łączenie zbiorów ASCII GRID w jedną całość
  4. Problem pikseli NoData w rastrach
  5. Konwersja zbiorów XYZ do rastra
  6. Konwersja zbiorów rastrowych do XYZ
  7. Ekstrakcja max wartości z rastra w obrębie poligonów
  8. Współrzędne geograficzne obiektów punktowych
  9. Statystyki fokalne
  10. Poster: Topologia geobazy - reguły topologiczne

BDOT

  1. Klasyfikacja obiektów BDOT10k

Inne

  1. Kalkulator Zadorskiego