t.rast.aggregate.dsgrass - Dalam Talian di Awan

JADUAL:

NAMA

t.rast.aggregate.ds - Mengagregatkan data set data raster ruang masa sedia ada menggunakan
selang masa set data masa ruang kedua.

Kata-kata berkaitan

temporal, pengagregatan, raster, masa

SINOPSIS

t.rast.aggregate.ds
t.rast.aggregate.ds - membantu
t.rast.aggregate.ds [-ns] input=nama sampel=nama [jenis=nama] output=nama nama asas=rentetan
kaedah=rentetan [mengimbangi=integer] [nprocs=integer] [pensampelan=nama[,nama,...]]
[di mana=sql_query] [--menimpa] [--membantu] [--kata kerja] [--tenang] [--ui]

Bendera:
-n
Daftar peta Null

-s
Gunakan masa mula - dipotong mengikut butiran - sebagai akhiran (mengatasi offset
pilihan)

--tulis ganti
Benarkan fail output untuk menulis ganti fail sedia ada

- membantu
Cetak ringkasan penggunaan

--verbose
Keluaran modul verbose

--senyap
Keluaran modul senyap

--ui
Paksa melancarkan dialog GUI

Parameter:
input=nama [diperlukan]
Nama set data raster masa ruang input

sampel=nama [diperlukan]
Selang masa daripada set data masa ruang ini (raster, vektor atau raster3d) digunakan untuk
pengiraan agregasi

jenis=nama
Jenis set data masa ruang pengagregatan
Pilihan: strds, stvds, str3ds
default: strds

output=nama [diperlukan]
Nama set data raster masa ruang keluaran

nama asas=rentetan [diperlukan]
Nama asas peta keluaran baharu yang dijana
Akhiran berangka yang dipisahkan dengan garis bawah akan dilampirkan untuk mencipta yang unik
pengecam

kaedah=rentetan [diperlukan]
Operasi agregat untuk dilakukan pada peta raster
Pilihan: purata, kira, median, mod, minimum, min_raster, maksimum, max_raster,
stddev, pelbagai, jumlah, kelainan, kepelbagaian, cerun, mengimbangi, detcoeff, sukuan1, sukuan3,
perc90, kuantil, kecondongan, kurtosis
default: purata

mengimbangi=integer
Offset yang digunakan untuk mencipta id peta output, id peta output dijana sebagai:
basename_ (kiraan + offset)
default: 0

nprocs=integer
Bilangan proses r.mapcalc untuk dijalankan secara selari
default: 1

pensampelan=nama [, nama,...]
Kaedah yang akan digunakan untuk pensampelan set data input
Pilihan: sama, bertindih, bertindih, bermula, bermula, selesai, selesai, semasa,
mengandungi
default: mengandungi

di mana=sql_query
WHERE syarat pernyataan SQL tanpa kata kunci 'di mana' digunakan dalam GIS temporal
rangka kerja
Contoh: masa_mula > '2001-01-01 12:30:00'

DESCRIPTION

t.rast.aggregate.ds berfungsi seperti t.rast.aggregate tetapi bukannya menentukan tetap
kebutiran untuk pengagregatan temporal selang masa semua peta yang didaftarkan dalam satu saat
set data masa ruang (boleh STRDS, STR3DS atau STVDS) digunakan untuk mengagregat peta bagi
masukkan dataset raster masa ruang.

NOTA

Kaedah pensampelan mesti dinyatakan dari sudut set data pensampel. Ia mentakrifkan
pinggul hubungan temporal antara selang set data pensampelan dan ruang input
set data raster masa.

CONTOH

Pemendakan pengagregatan
Dalam contoh ini kami mencipta 7 peta raster yang akan didaftarkan dalam satu ruang masa
dataset raster dinamakan hujan_harian menggunakan butiran temporal harian. Nama-nama bagi
peta raster disimpan dalam fail teks yang digunakan untuk pendaftaran peta raster.

Set data vektor ruang masa dicipta daripada dua peta vektor dengan temporal yang berbeza
resolusi. Peta dicipta menggunakan v.random. Peta pertama mempunyai butiran 3 hari
kedua butir 4 hari.

Dataset raster masa ruang hujan_harian dengan butiran temporal harian akan
diagregatkan menggunakan dataset vektor masa ruang yang menghasilkan raster masa ruang keluaran
dataset hujan_agg. Kaedah pengagregatan ditetapkan kepada jumlah untuk mengumpul
nilai kerpasan semua selang dalam set data vektor masa ruang. persampelan
pilihan memastikan bahawa hanya peta raster yang bersifat sementara semasa selang masa
set data vektor masa ruang dipertimbangkan untuk pengiraan. Oleh itu pilihan ditetapkan kepada
mengandungi (lapisan peta vektor dicap masa secara sementara mengandungi lapisan peta raster):
MAPS="peta_1 peta_2 peta_3 peta_4 peta_5 peta_6 peta_7"
untuk peta dalam ${MAPS} ; buat
r.mapcalc expression="${map} = 1"
echo ${map} >> map_list.txt
dilakukan
t.create type=strds temporaltype=absolute \
output=kerpasan_harian \
title="Pemendakan harian" \
description="Uji set data dengan kerpasan harian"
t.register -i type=raster input=precipitation_daily \
file=map_list.txt start="2012-08-20" kenaikan="1 hari"
t.info type=strds input=precipitation_daily
+-------------------- Set Data Raster Masa Ruang ------------------------- ----+
| |
+-------------------- Maklumat asas ---------------------------------------- ----------+
| Id: ........................ kerpasan_harian@KEKAL
| Nama: ...................... curah hujan_harian
| Mapset: .................... KEKAL
| Pencipta: ................... soeren
| Jenis sementara: ............. mutlak
| Masa penciptaan: ............. 2014-11-23 16:48:17.686979
| Masa pengubahsuaian:.......... 2014-11-23 16:48:18.302978
| Jenis semantik:.............. min
+-------------------- Masa mutlak ---------------------------------------- --------------+
| Masa mula:................. 2012-09-10 00:00:00
| Masa tamat:................... 2012-09-17 00:00:00
| Kebutiran:................ 1 hari
| Jenis peta sementara:...... selang
+--------------------- Tahap ruang ---------------------------------------- -------------+
| Utara:...................... 80.0
| Selatan:...................... 0.0
| Timur:.. .................... 120.0
| Barat:....................... 0.0
| Atas:........................ 0.0
| Bawah:..................... 0.0
+-------------------- Maklumat metadata ---------------------------------------- -------+
| Jadual daftar raster:...... raster_map_register_3225725979b14b5db343a00835b882c7
| Resolusi Utara-Selatan min:. 10.0
| Resolusi Utara-Selatan maks:. 10.0
| Resolusi timur-barat min:... 10.0
| Maksimum resolusi timur-barat:... 10.0
| Nilai minimum min:.......... 1.0
| Nilai minimum maks:.......... 1.0
| Nilai maksimum min:.......... 1.0
| Nilai maksimum maks:.......... 1.0
| Jenis pengagregatan:........... Tiada
| Bilangan peta berdaftar:.. 7
|
| Tajuk:
| Hujan harian
| Penerangan:
| Uji set data dengan kerpasan harian
| Sejarah arahan:
| # 2014-11-23 16:48:17
| t.create type="strds" temporaltype="absolute"
| output="kerpasan_harian" title="kerpasan harian"
| description="Uji set data dengan kerpasan harian"
| # 2014-11-23 16:48:18
| t.register -i type="rast" input="precipitation_daily"
| file="map_list.txt" start="2012-08-20" kenaikan="1 hari"
|
+------------------------------------------------ ---------------------------+
v.keluaran rawak=mata_1 n=20
v.keluaran rawak=mata_2 n=20
t.create type=stvds temporaltype=absolute \
keluaran=mata \
title="Mata" \
description="Mata untuk pengagregatan"
t.register -i type=input vektor=titik \
map=points_1 start="2012-08-20" increment="3 hari"
t.register -i type=input vektor=titik \
map=points_2 start="2012-08-23" increment="4 hari"
t.info type=stvds input=points
+-------------------- Set Data Vektor Masa Angkasa ------------------------- ----+
| |
+-------------------- Maklumat asas ---------------------------------------- ----------+
| Id: ........................ mata@TETAP
| Nama: ...................... mata
| Mapset: .................... KEKAL
| Pencipta: ................... soeren
| Jenis sementara: ............. mutlak
| Masa penciptaan: ............. 2014-11-23 16:48:49.193903
| Masa pengubahsuaian:.......... 2014-11-23 16:48:50.185671
| Jenis semantik:.............. min
+-------------------- Masa mutlak ---------------------------------------- --------------+
| Masa mula:................. 2012-08-20 00:00:00
| Masa tamat:................... 2012-08-27 00:00:00
| Kebutiran:................ 1 hari
| Jenis peta sementara:...... selang
+--------------------- Tahap ruang ---------------------------------------- -------------+
| Utara:...................... 79.283411
| Selatan:...................... 5.724954
| Timur:.. .................... 118.881168
| Barat:....................... 0.016755
| Atas:........................ 0.0
| Bawah:..................... 0.0
+-------------------- Maklumat metadata ---------------------------------------- -------+
| Vector register table:...... vector_map_register_6f02d33e0ee243d1a521aaaca39ecb31
| Bilangan mata ........... 40
| Bilangan baris ............ 0
| Bilangan sempadan ....... 0
| Bilangan centroid ........ 0
| Bilangan muka ............ 0
| Bilangan biji .......... 0
| Bilangan primitif ....... 40
| Bilangan nod ............ 0
| Bilangan kawasan ............ 0
| Bilangan pulau .......... 0
| Bilangan lubang ............ 0
| Bilangan jilid .......... 0
| Bilangan peta berdaftar:.. 2
|
| Tajuk:
| mata
| Penerangan:
| Mata untuk pengagregatan
| Sejarah arahan:
| # 2014-11-23 16:48:49
| t.create type="stvds" temporaltype="absolute"
| output="points" title="Points" description="Mata untuk pengagregatan"
| # 2014-11-23 16:48:49
| t.register -i type="vect" input="points"
| map="points_1" start="2012-08-20" increment="3 hari"
| # 2014-11-23 16:48:50
| t.register -i type="vect" input="points"
| map="points_2" start="2012-08-23" increment="4 hari"
|
+------------------------------------------------ ---------------------------+
t.rast.aggregate.ds input=precipitation_daily \
keluaran=kerpasan_agg \
sampel=jenis mata=stvds \
basename=prec_agg \
kaedah=jumlah persampelan=mengandungi
input t.support=precipitation_agg \
title="Aggregated precipitation" \
description="Agregat precipitation dataset"
t.info type=strds input=precipitation_agg
+-------------------- Set Data Raster Masa Ruang ------------------------- ----+
| |
+-------------------- Maklumat asas ---------------------------------------- ----------+
| Id: ........................ kerpasan_agg@KEKAL
| Nama: ...................... hujan_agg
| Mapset: .................... KEKAL
| Pencipta: ................... soeren
| Jenis sementara: ............. mutlak
| Masa penciptaan: ............. 2014-11-23 16:53:23.488799
| Masa pengubahsuaian:.......... 2014-11-23 16:53:28.714886
| Jenis semantik:.............. min
+-------------------- Masa mutlak ---------------------------------------- --------------+
| Masa mula:................. 2012-08-20 00:00:00
| Masa tamat:................... 2012-08-27 00:00:00
| Kebutiran:................ 1 hari
| Jenis peta sementara:...... selang
+--------------------- Tahap ruang ---------------------------------------- -------------+
| Utara:...................... 80.0
| Selatan:...................... 0.0
| Timur:.. .................... 120.0
| Barat:....................... 0.0
| Atas:........................ 0.0
| Bawah:..................... 0.0
+-------------------- Maklumat metadata ---------------------------------------- -------+
| Raster register table:...... raster_map_register_7b025eb7431747c98c5c1ad971e8c282
| Resolusi Utara-Selatan min:. 10.0
| Resolusi Utara-Selatan maks:. 10.0
| Resolusi timur-barat min:... 10.0
| Maksimum resolusi timur-barat:... 10.0
| Nilai minimum min:.......... 3.0
| Nilai minimum maks:.......... 4.0
| Nilai maksimum min:.......... 3.0
| Nilai maksimum maks:.......... 4.0
| Jenis pengagregatan:........... jumlah
| Bilangan peta berdaftar:.. 2
|
| Tajuk:
| Kerpasan agregat
| Penerangan:
| Set data kerpasan agregat
| Sejarah arahan:
| # 2014-11-23 16:53:23
| t.rast.aggregate.ds input="precipitation_daily"
| output="precipitation_agg" sample="points" type="stvds" basename="prec_agg"
| method="sum" sampling="contains"
| # 2014-11-23 16:53:28
| t.support input="precipitation_agg"
| title="Agregat kerpasan"
| description="Agregat precipitation dataset"
|
+------------------------------------------------ ---------------------------+

MODIS satelit sensor harian data pengagregatan kepada 8 hari
Dalam contoh ini pengagregatan daripada data harian kepada lapan hari ditunjukkan. Ini "lapan hari
minggu" digunakan dalam beberapa produk penderia satelit MODIS.
# NOTA: contoh ditulis dalam bahasa shell
# buat peta setiap 8 hari sebagai peta benih
untuk tahun dalam `seq 2000 2001`; buat
untuk doy dalam `seq -w 1 8 365` ; buat
r.mapcalc -s expression="8day_${year}_${doy} = rand(0.0,40.0)"
dilakukan
dilakukan
# Daripada nama setiap peta, kami mengambil tahun dan hari, dan menukarnya
# kepada tarikh YYYY-MM-DD untuk mula dan tamat, dan buat fail dengan
# nama peta, tarikh mula dan tarikh tamat
g.list type=raster pattern=8day_20??_* > names_list
untuk NAME dalam `cat names_list` ; buat
# Menghuraikan
YEAR=`echo $NAME | potong -d'_' -f2`
DOY=`echo $NAME | potong -d'_' -f3`
# tukar YYYY_DOY kepada YYYY-MM-DD
DOY=`echo "$DOY" | sed 's/^0*//'`
doy_end=0
jika [ $DOY -le "353" ] ; kemudian
doy_end=$(( $DOY + 8 ))
elif [ $DOY -eq "361" ] ; kemudian
jika [ $[$TAHUN % 4] -eq 0 ] && [ $[$TAHUN % 100] -ne 0 ] || [ $[$TAHUN % 400] -eq 0 ] ; kemudian
doy_end=$(( $DOY + 6 ))
lagi
doy_end=$(( $DOY + 5 ))
fi
fi
DATE_START=`tarikh -d "${YEAR}-01-01 +$(( ${DOY} - 1 ))hari" +%Y-%m-%d`
DATE_END=`tarikh -d "${YEAR}-01-01 +$(( ${doy_end} -1 ))hari" +%Y-%m-%d`
# fail teks dengan nama peta, tarikh mula dan tarikh tamat
echo "$NAME|$DATE_START|$DATE_END" >> list_map_start_end_time.txt
dilakukan
# semak senarai yang dibuat.
cat list_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
# semua peta kecuali peta terakhir dalam setiap tahun mewakili 8 hari
# selang waktu. Tetapi pengagregatan bermula sekali lagi setiap
# 1 Januari.
# buat strds seperti MODIS 8 hari
t.create type=strds temporaltype=absolute \
output=8day_ts title="siri masa 8 hari" \
description="STRDS dengan MODIS seperti pengagregatan 8 hari"
# daftar peta
t.register type=raster input=8day_ts \
file=list_map_start_end_time.txt
# semak
input t.info=8hari_ts
t.rast.list input=8day_ts
# akhirnya, salin pengagregatan ke siri masa harian
t.rast.aggregate.ds -s input=daily_ts sample=8day_ts \
output=8day_agg basename=8day_agg \
kaedah=persampelan purata=mengandungi
# tambah metadata
input t.support=8day_agg \
title="8 hari agregat ts" \
description="set data agregat seperti MODIS 8 hari"
# semak senarai peta dalam strds agregat yang baru dibuat
t.rast.list input=8day_agg
name|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

Gunakan t.rast.aggregate.dsgrass dalam talian menggunakan perkhidmatan onworks.net