t.rast.aggregate.dsgrass - Online w chmurze

PROGRAM:

IMIĘ

t.rast.aggregate.ds - Agreguje dane z istniejącego rastrowego zbioru danych przestrzenno-czasowych za pomocą
przedziały czasowe drugiego zbioru danych przestrzenno-czasowych.

SŁOWA KLUCZOWE

czasowy, agregacja, raster, czas

STRESZCZENIE

t.rast.aggregate.ds
t.rast.aggregate.ds --help
t.rast.aggregate.ds [-ns] wkład=Nazwa próba=Nazwa [rodzaj=Nazwa] wydajność=Nazwa nazwa podstawowa=ciąg
metoda=ciąg [offset=liczba całkowita] [Nprocs=liczba całkowita] [próbkowania=Nazwa[,Nazwa,...]]
[gdzie=zapytanie_sql] [--przepisać] [--pomoc] [--gadatliwy] [--cichy] [--ui]

Flagi:
-n
Zarejestruj mapy Null

-s
Użyj czasu rozpoczęcia — obciętego zgodnie z szczegółowością — jako sufiksu (zastępuje przesunięcie
opcjonalnie)

--przepisać
Zezwalaj plikom wyjściowym na zastępowanie istniejących plików

--help
Wydrukuj podsumowanie wykorzystania

--gadatliwy
Pełne wyjście modułu

--cichy
Cichy moduł wyjściowy

--UI
Wymuś uruchomienie okna GUI

Parametry:
wkład=Nazwa [wymagany]
Nazwa wejściowego zbioru danych rastrowych czasoprzestrzennych

próba=Nazwa [wymagany]
Wykorzystywane są przedziały czasowe z tego zbioru danych czasoprzestrzennych (raster, wektor lub raster3d).
obliczenia agregacji

rodzaj=Nazwa
Typ zbioru danych przestrzenno-czasowych agregacji
Opcje: kroki, płyty, str3ds
Zaniedbanie: strd

wydajność=Nazwa [wymagany]
Nazwa wyjściowego zestawu danych rastrowych przestrzeni i czasu

nazwa podstawowa=ciąg [wymagany]
Nazwa podstawowa nowo wygenerowanych map wyjściowych
Sufiks liczbowy oddzielony znakiem podkreślenia zostanie dołączony, aby utworzyć unikat
identyfikator

metoda=ciąg [wymagany]
Operacja agregująca do wykonania na mapach rastrowych
Opcje: średni, liczyć, mediana, tryb, minimum, min_raster, maksymalny, max_raster,
standardowe odchylenie, Zakres, suma, zmienność, różnorodność, nachylenie, offsetowy, Detcoeff, kwarta1, kwarta3,
perc90, kwantyl, skośność, kurtoza
Zaniedbanie: średni

offset=liczba całkowita
Przesunięcie, które jest używane do tworzenia identyfikatorów map wyjściowych, identyfikator mapy wyjściowej jest generowany jako:
nazwa_bazowa_ (liczba + przesunięcie)
Zaniedbanie: 0

Nprocs=liczba całkowita
Liczba procesów r.mapcalc, które mają być uruchamiane równolegle
Zaniedbanie: 1

próbkowania=Nazwa nazwa,...]
Metoda, która ma być użyta do próbkowania wejściowego zbioru danych
Opcje: równy, nakładanie się, zachodzące na siebie, zaczyna się, zaczęła się, wykończenia, skończone, w trakcie,
zawiera
Zaniedbanie: zawiera

gdzie=zapytanie_sql
WHERE warunki instrukcji SQL bez słowa kluczowego „where” użyte w czasowym GIS

Przykład: czas_początku > '2001-01-01 12:30:00'

OPIS

t.rast.aggregate.ds działa jak t.rast.aggregate, ale zamiast definiować stałą
szczegółowość agregacji czasowej przedziałów czasowych wszystkich map zarejestrowanych w ciągu sekundy
zbiory danych czasoprzestrzennych (mogą to być STRDS, STR3DS lub STVDS) są wykorzystywane do agregowania map
wejściowy rastrowy zbiór danych przestrzenno-czasowych.

UWAGI

Metodę próbkowania należy określić z punktu widzenia zbioru danych próbnika. To definiuje
relacje czasowe są przesunięte między przedziałami zbioru danych próbkowania a przestrzenią wejściową
zbiór danych rastrowych czasu.

PRZYKŁADY

Osad zbiór
W tym przykładzie tworzymy 7 map rastrowych, które zostaną zarejestrowane w jednej czasoprzestrzeni
rastrowy zbiór danych o nazwie opady_dzienne przy użyciu codziennej granulacji czasowej. Nazwy
mapy rastrowe są przechowywane w pliku tekstowym, który służy do rejestracji mapy rastrowej.

Zbiór danych wektorowych czasoprzestrzennych jest tworzony z dwóch map wektorowych o różnym czasie
rezolucja. Mapy są tworzone przy użyciu v.random. Pierwsza mapa ma granulat 3-dniowy
drugi granulat 4 dni.

Rastrowy zbiór danych czasoprzestrzennych opady_dzienne z codzienną szczegółowością czasową
zagregowane przy użyciu zbioru danych wektorów przestrzenno-czasowych, w wyniku czego otrzymano wyjściowy raster przestrzenno-czasowy
zestaw danych opady_agg. Metoda agregacji jest ustawiona na suma gromadzić
wartości opadów we wszystkich przedziałach zbioru danych wektorów przestrzenno-czasowych. Próbkowanie
Opcja zapewnia, że ​​tylko mapy rastrowe, które są tymczasowe w przedziałach czasowych
Do obliczeń uwzględnia się zbiory danych wektorów czasoprzestrzennych. Dlatego opcja jest ustawiona na
zawiera (tymczasowo warstwy mapy wektorowej ze znacznikiem czasu zawierać warstwy mapy rastrowej):
MAPS="mapa_1 mapa_2 mapa_3 mapa_4 mapa_5 mapa_6 mapa_7"
dla mapy w ${MAPS} ; Do
r.mapcalc wyrażenie="${mapa} = 1"
echo ${map} >> map_list.txt
zrobić
t.create type=strds temporaltype=absolutne \
wyjście = opad_dzienny \
title="Dzienne opady" \
description="Testowy zestaw danych z dziennymi opadami"
t.register -i type=raster input=opad_dzienny \
file=map_list.txt start="2012-08-20" przyrost="1 dzień"
t.info type=strds input=opady_dzienne
+-------------------- Zbiór danych rastrowych czasoprzestrzennych ------------------------- ----+
| |
+-------------------- Podstawowe informacje ------------------------------------------ ----------+
| Identyfikator: ............................... opad_dzienny@PERMANENT
| Nazwa: ...........................opad_dzienny
| Zestaw map: ................... STAŁY
| Twórca: ................... soeren
| Typ czasowy: ............... absolutny
| Czas utworzenia: ............... 2014-11-23 16:48:17.686979
| Czas modyfikacji:.......... 2014-11-23 16:48:18.302978
| Typ semantyczny:............ średni
+-------------------- Czas bezwzględny ------------------------------------------ --------------+
| Godzina rozpoczęcia: .............. 2012-09-10 00:00:00
| Godzina zakończenia: .............. 2012-09-17 00:00:00
| Ziarnistość:............ 1 dzień
| Temporalny typ map:......interwałowy
+-------------------- Zasięg przestrzenny --------------------------- -------------+
| Północ:........................... 80.0
| Południe:........................... 0.0
| Wschód:............................. 120.0
| Zachód:........................... 0.0
| Góra:........................... 0.0
| Dół:........................... 0.0
+-------------------- Informacje o metadanych --------------------------- -------+
| Tabela rejestrów rastrowych:...... raster_map_register_3225725979b14b5db343a00835b882c7
| Rozdzielczość północ-południe min:. 10.0
| Maks. rozdzielczość północ-południe:. 10.0
| Rozdzielczość wschód-zachód min:... 10.0
| Rozdzielczość wschód-zachód max:... 10.0
| Minimalna wartość min:......... 1.0
| Minimalna wartość maks.:......... 1.0
| Maksymalna wartość min.:......... 1.0
| Maksymalna wartość maks.:......... 1.0
| Typ agregacji:........... Brak
| Liczba zarejestrowanych map:.. 7
|
| Tytuł:
| Opady dzienne
| Opis:
| Zestaw danych testowych z dziennymi opadami
| Historia poleceń:
| # 2014-11-23 16:48:17
| t.create type="strds" temporaltype="absolute"
| output="precipitation_daily" title="Dzienne opady"
| description="Testowy zestaw danych z dziennymi opadami"
| # 2014-11-23 16:48:18
| t.register -i type="rast" input="opady_dzienne"
| file="map_list.txt" start="2012-08-20" przyrost="1 dni"
|
+------------------------------------------------- ---------------------------+
v.losowe wyjście=punkty_1 n=20
v.losowe wyjście=punkty_2 n=20
t.create type=stvds temporaltype=absolutny \
wynik = punkty \
tytuł="Punkty" \
opis="Punkty za agregację"
t.register -i typ=wejście wektorowe=punkty \
map=points_1 start="2012-08-20" przyrost="3 dni"
t.register -i typ=wejście wektorowe=punkty \
map=points_2 start="2012-08-23" przyrost="4 dni"
t.info typ = wejście stvds = punkty
+-------------------- Zbiór danych wektorowych czasoprzestrzeni ------------------------- ----+
| |
+-------------------- Podstawowe informacje ------------------------------------------ ----------+
| Identyfikator: .......................punkty@PERMANENT
| Imię i nazwisko: ............... punktów
| Zestaw map: ................... STAŁY
| Twórca: ................... soeren
| Typ czasowy: ............... absolutny
| Czas utworzenia: ............... 2014-11-23 16:48:49.193903
| Czas modyfikacji:.......... 2014-11-23 16:48:50.185671
| Typ semantyczny:............ średni
+-------------------- Czas bezwzględny ------------------------------------------ --------------+
| Godzina rozpoczęcia: .............. 2012-08-20 00:00:00
| Godzina zakończenia: .............. 2012-08-27 00:00:00
| Ziarnistość:............ 1 dzień
| Temporalny typ map:......interwałowy
+-------------------- Zasięg przestrzenny --------------------------- -------------+
| Północ:........................... 79.283411
| Południe:........................... 5.724954
| Wschód:............................. 118.881168
| Zachód:........................... 0.016755
| Góra:........................... 0.0
| Dół:........................... 0.0
+-------------------- Informacje o metadanych --------------------------- -------+
| Vector register table:...... vector_map_register_6f02d33e0ee243d1a521aaaca39ecb31
| Liczba punktów .............. 40
| Liczba linii ............ 0
| Liczba granic ........ 0
| Liczba centroidów ........ 0
| Liczba twarzy .............. 0
| Liczba jąder ........ 0
| Liczba prymitywów ........ 40
| Liczba węzłów ............ 0
| Liczba obszarów .............. 0
| Liczba wysp ............... 0
| Liczba otworów .............. 0
| Liczba tomów ............... 0
| Liczba zarejestrowanych map:.. 2
|
| Tytuł:
| Zwrotnica
| Opis:
| Punkty za agregację
| Historia poleceń:
| # 2014-11-23 16:48:49
| t.create type="stvds" temporaltype="absolute"
| wyjście="punkty" title="Punkty" opis="Punkty za agregację"
| # 2014-11-23 16:48:49
| t.register -i type="vect" input="points"
| map="points_1" start="2012-08-20" przyrost="3 dni"
| # 2014-11-23 16:48:50
| t.register -i type="vect" input="points"
| map="points_2" start="2012-08-23" przyrost="4 dni"
|
+------------------------------------------------- ---------------------------+
t.rast.aggregate.ds input=opady_dzienne \
wyjście=opady_agg \
próbka=typ punktów=stvds \
nazwa_bazowa=prec_agg \
metoda=sumowanie próbkowania=zawiera
t.support wejście=opady_agg \
title="Zagregowane opady" \
opis="Zagregowany zbiór danych o opadach"
t.info typ=strds wejście=opady_agg
+-------------------- Zbiór danych rastrowych czasoprzestrzennych ------------------------- ----+
| |
+-------------------- Podstawowe informacje ------------------------------------------ ----------+
| Identyfikator:........................... opady_agg@PERMANENT
| Nazwa: ...................... opady_agg
| Zestaw map: ................... STAŁY
| Twórca: ................... soeren
| Typ czasowy: ............... absolutny
| Czas utworzenia: ............... 2014-11-23 16:53:23.488799
| Czas modyfikacji:.......... 2014-11-23 16:53:28.714886
| Typ semantyczny:............ średni
+-------------------- Czas bezwzględny ------------------------------------------ --------------+
| Godzina rozpoczęcia: .............. 2012-08-20 00:00:00
| Godzina zakończenia: .............. 2012-08-27 00:00:00
| Ziarnistość:............ 1 dzień
| Temporalny typ map:......interwałowy
+-------------------- Zasięg przestrzenny --------------------------- -------------+
| Północ:........................... 80.0
| Południe:........................... 0.0
| Wschód:............................. 120.0
| Zachód:........................... 0.0
| Góra:........................... 0.0
| Dół:........................... 0.0
+-------------------- Informacje o metadanych --------------------------- -------+
| Raster register table:...... raster_map_register_7b025eb7431747c98c5c1ad971e8c282
| Rozdzielczość północ-południe min:. 10.0
| Maks. rozdzielczość północ-południe:. 10.0
| Rozdzielczość wschód-zachód min:... 10.0
| Rozdzielczość wschód-zachód max:... 10.0
| Minimalna wartość min:......... 3.0
| Minimalna wartość maks.:......... 4.0
| Maksymalna wartość min.:......... 3.0
| Maksymalna wartość maks.:......... 4.0
| Typ agregacji:........... suma
| Liczba zarejestrowanych map:.. 2
|
| Tytuł:
| Zagregowane opady
| Opis:
| Zagregowany zbiór danych o opadach
| Historia poleceń:
| # 2014-11-23 16:53:23
| t.rast.aggregate.ds input="opady_dzienne"
| wyjście="precipitation_agg" sample="points" type="stvds" basename="prec_agg"
| metoda="suma" próbkowanie="zawiera"
| # 2014-11-23 16:53:28
| t.support input="precipitation_agg"
| title="Zagregowane opady"
| opis="Zagregowany zbiór danych o opadach"
|
+------------------------------------------------- ---------------------------+

MODIS satelita czujnik codziennie dane zbiór do 8 dni
W tym przykładzie pokazano agregację danych dziennych do ośmiu dni. To „osiem dni
tydzień” jest stosowany w niektórych czujnikach satelitarnych MODIS.
# UWAGA: przykład jest napisany w języku powłoki
# twórz mapy co 8 dni jako mapy początkowe
dla roku w „nast. 2000 2001”; Do
dla doy w `seq -w 1 8 365` ; Do
r.mapcalc -s wyrażenie="8day_${rok}_${doy} = rand(0.0,40.0)"
zrobić
zrobić
# Z nazwy każdej mapy bierzemy rok i doba i konwertujemy je
# na datę RRRR-MM-DD początku i końca i utwórz plik z
# nazwy map, data początkowa i data końcowa
g.list type=wzorzec rastrowy=8day_20??_* > lista_nazw
dla IMIĘ w `liście_imion kotów`; Do
# Przeanalizuj
ROK=`echo $NAZWA | cięcie -d'_' -f2`
DOY=`echo $NAZWA | cięcie -d'_' -f3`
# przekonwertuj YYYY_DOY na YYYY-MM-DD
DOY=`echo "$DOY" | sed 's/^0*//'`
doy_end=0
if [$DOY -le "353"] ; Następnie
doy_end=$(( $DOY + 8 ))
elif [ $DOY -eq "361"] ; Następnie
if [ $[$ROK % 4] -eq 0 ] && [ $[$ROK % 100] -ne 0 ] || [ $[$ROK % 400] -eq 0 ] ; Następnie
doy_end=$(( $DOY + 6 ))
więcej
doy_end=$(( $DOY + 5 ))
fi
fi
DATE_START=`data -d "${YEAR}-01-01 +$(( ${DOY} - 1 ))dni" +%Y-%m-%d`
DATE_END=`data -d "${YEAR}-01-01 +$(( ${doy_end} -1 ))dni" +%Y-%m-%d`
# plik tekstowy z nazwami map, datą początkową i końcową
echo "$NAME|$DATE_START|$DATE_END" >> list_map_start_end_time.txt
zrobić
# sprawdź utworzoną listę.
lista kotów_map_start_end_time.txt
8day_2000_001|2000-01-01|2000-01-09
8day_2000_009|2000-01-09|2000-01-17
...
8day_2000_353|2000-12-18|2000-12-26
8day_2000_361|2000-12-26|2001-01-01
8day_2001_001|2001-01-01|2001-01-09
8day_2001_009|2001-01-09|2001-01-17
...
8day_2001_345|2001-12-11|2001-12-19
8day_2001_353|2001-12-19|2001-12-27
8day_2001_361|2001-12-27|2002-01-01
# wszystkie mapy z wyjątkiem ostatniej mapy w każdym roku reprezentują 8 dni
# interwały. Ale agregacja zaczyna się od nowa
# 1 stycznia.
# utwórz 8-dniowe strd-podobne do MODIS
t.create type=strds temporaltype=absolutne \
wyjście=8day_ts title="8-dniowy szereg czasowy" \
opis="STRDS z MODIS jak agregacja 8-dniowa"
# zarejestruj mapy
t.register type=rastrowe wejście=8day_ts \
plik=list_map_start_end_time.txt
# sprawdzać
t.info wejście = 8day_ts
Wejście t.rast.list=8day_ts
# na koniec skopiuj agregację do dziennego szeregu czasowego
t.rast.aggregate.ds -s wejście=dzienne_ts próbka=8dni_ts \
wyjście=8day_agg nazwabazowa=8day_agg \
metoda=średnia próbka=zawiera
# dodaj metadane
t. wsparcie input=8day_agg \
title="Zagregowany ts z 8 dni" \
opis="8-dniowy zagregowany zbiór danych podobny do MODIS"
# sprawdź listę map w nowo utworzonych zagregowanych plikach strd
t.rast.list wejście=8day_agg
nazwa|mapset|start_time|end_time
8day_agg_2000_01_01|modis|2000-01-01 00:00:00|2000-01-09 00:00:00
8day_agg_2000_01_09|modis|2000-01-09 00:00:00|2000-01-17 00:00:00
8day_agg_2000_01_17|modis|2000-01-17 00:00:00|2000-01-25 00:00:00
...
8day_agg_2000_12_18|modis|2000-12-18 00:00:00|2000-12-26 00:00:00
8day_agg_2000_12_26|modis|2000-12-26 00:00:00|2001-01-01 00:00:00
8day_agg_2001_01_01|modis|2001-01-01 00:00:00|2001-01-09 00:00:00
...
8day_agg_2001_12_11|modis|2001-12-11 00:00:00|2001-12-19 00:00:00
8day_agg_2001_12_19|modis|2001-12-19 00:00:00|2001-12-27 00:00:00
8day_agg_2001_12_27|modis|2001-12-27 00:00:00|2002-01-01 00:00:00

Użyj t.rast.aggregate.dsgrass online, korzystając z usług onworks.net