นี่คือคำสั่ง t.rast.aggregate.dsgrass ที่สามารถเรียกใช้ในผู้ให้บริการโฮสติ้งฟรีของ OnWorks โดยใช้หนึ่งในเวิร์กสเตชันออนไลน์ฟรีของเรา เช่น Ubuntu Online, Fedora Online, โปรแกรมจำลองออนไลน์ของ Windows หรือโปรแกรมจำลองออนไลน์ของ MAC OS
โครงการ:
ชื่อ
t.rast.aggregate.ds - รวบรวมข้อมูลของชุดข้อมูลแรสเตอร์เวลาพื้นที่ที่มีอยู่โดยใช้
ช่วงเวลาของชุดข้อมูลเวลาช่องว่างที่สอง
KEYWORDS
ชั่วขณะ, การรวมตัว, แรสเตอร์, เวลา
เรื่องย่อ
t.rast.aggregate.ds
t.rast.aggregate.ds --ช่วยด้วย
t.rast.aggregate.ds [-ns] อินพุต=พร้อมชื่อ ตัวอย่าง=พร้อมชื่อ [ชนิด=พร้อมชื่อ] เอาท์พุต=พร้อมชื่อ ชื่อฐาน=เชือก
วิธี=เชือก [ชดเชย=จำนวนเต็ม] [Nprocs=จำนวนเต็ม] [การสุ่มตัวอย่าง=พร้อมชื่อ[,พร้อมชื่อ-
[ที่ไหน=sql_query-เขียนทับ-ช่วย-ละเอียด-เงียบสงบ-ui]
ธง:
-n
ลงทะเบียน Null maps
-s
ใช้เวลาเริ่มต้น - ตัดทอนตามความละเอียด - เป็นคำต่อท้าย (แทนที่ offset
ตัวเลือก)
--เขียนทับ
อนุญาตให้ไฟล์เอาต์พุตเขียนทับไฟล์ที่มีอยู่
--ช่วยด้วย
พิมพ์สรุปการใช้งาน
--รายละเอียด
เอาต์พุตโมดูล verbose
--เงียบ
เอาต์พุตโมดูลเงียบ
--UI
บังคับให้เปิดใช้กล่องโต้ตอบ GUI
พารามิเตอร์:
อินพุต=พร้อมชื่อ [ที่จำเป็น]
ชื่อของชุดข้อมูลแรสเตอร์เวลาพื้นที่อินพุต
ตัวอย่าง=พร้อมชื่อ [ที่จำเป็น]
ช่วงเวลาจากชุดข้อมูลเวลาของพื้นที่นี้ (raster, vector หรือ raster3d) ใช้สำหรับ
การคำนวณรวม
ชนิด=พร้อมชื่อ
ประเภทของชุดข้อมูลเวลาพื้นที่รวม
ตัวเลือก: สตริด, เอสทีวีดี, str3ds
ค่าเริ่มต้น: เส้น
เอาท์พุต=พร้อมชื่อ [ที่จำเป็น]
ชื่อของชุดข้อมูลแรสเตอร์เวลาพื้นที่เอาท์พุต
ชื่อฐาน=เชือก [ที่จำเป็น]
ชื่อพื้นฐานของแมปเอาท์พุตที่สร้างขึ้นใหม่
คำต่อท้ายตัวเลขที่คั่นด้วยขีดล่างจะถูกแนบเพื่อสร้างเอกลักษณ์
ระบุ
วิธี=เชือก [ที่จำเป็น]
การดำเนินการโดยรวมที่จะดำเนินการบนแผนที่แรสเตอร์
ตัวเลือก: เฉลี่ย, นับ, ค่ามัธยฐาน โหมด, ขั้นต่ำ min_raster, ขีดสุด, max_raster,
มาตรฐาน พิสัย, ผลรวม ความแปรปรวน ความหลากหลาย, ความลาดชัน ชดเชย เด็คโคออฟ, ควอร์ต1, ควอร์ต3,
เปอร์เซ็น90, ปริมาณ, ความเบ้ โด่ง
ค่าเริ่มต้น: เฉลี่ย
ชดเชย=จำนวนเต็ม
ออฟเซ็ตที่ใช้ในการสร้างรหัสแมปเอาท์พุต รหัสแมปเอาท์พุตจะถูกสร้างขึ้นเป็น:
basename_ (นับ + ออฟเซ็ต)
ค่าเริ่มต้น: 0
Nprocs=จำนวนเต็ม
จำนวนกระบวนการ r.mapcalc ที่จะรันแบบขนาน
ค่าเริ่มต้น: 1
การสุ่มตัวอย่าง=ชื่อ[,ชื่อ-
วิธีที่จะใช้สุ่มตัวอย่างชุดข้อมูลอินพุต
ตัวเลือก: เท่ากับ, ทับซ้อนกัน ทับซ้อนกัน เริ่ม, เริ่ม, เสร็จสิ้น เสร็จแล้ว ในระหว่าง,
มี
ค่าเริ่มต้น: มี
ที่ไหน=sql_query
WHERE เงื่อนไขของคำสั่ง SQL ที่ไม่มีคีย์เวิร์ด 'where' ที่ใช้ใน GIS ชั่วคราว
กรอบ
ตัวอย่าง: start_time > '2001-01-01 12:30:00'
DESCRIPTION
t.rast.aggregate.ds ทำงานเหมือน t.rast.aggregate แต่แทนที่จะกำหนดคงที่
ความละเอียดสำหรับการรวมชั่วขณะ ช่วงเวลาของแผนที่ทั้งหมดที่ลงทะเบียนเป็นวินาที
ชุดข้อมูลเวลาอวกาศ (สามารถเป็น STRDS, STR3DS หรือ STVDS) ใช้เพื่อรวมแผนที่ของ
ชุดข้อมูลแรสเตอร์เวลาพื้นที่อินพุต
หมายเหตุ
ต้องระบุวิธีการสุ่มตัวอย่างจากมุมมองของชุดข้อมูลตัวอย่าง มันกำหนด
สะโพกความสัมพันธ์ชั่วคราวระหว่างช่วงเวลาของชุดข้อมูลสุ่มตัวอย่างและช่องว่างอินพุต
ชุดข้อมูลแรสเตอร์เวลา
ตัวอย่าง
การเร่งรัด การรวมตัว
ในตัวอย่างนี้ เราสร้าง 7 แผนที่แรสเตอร์ที่จะลงทะเบียนในช่องว่างเดียว
ชุดข้อมูลแรสเตอร์ชื่อ ฝน_รายวัน โดยใช้ความละเอียดชั่วขณะรายวัน ชื่อของ
แผนที่แรสเตอร์จะถูกเก็บไว้ในไฟล์ข้อความที่ใช้สำหรับการลงทะเบียนแผนที่แรสเตอร์
ชุดข้อมูลเวกเตอร์เวลาอวกาศถูกสร้างขึ้นจากแผนที่เวกเตอร์สองแผนที่ที่มีเวลาต่างกัน
ปณิธาน. แผนที่ถูกสร้างขึ้นโดยใช้ v.random แผนที่แรกมีเม็ด 3 วัน
เม็ดที่สองของ 4 วัน
ชุดข้อมูลแรสเตอร์เวลาอวกาศ ฝน_รายวัน ด้วยความละเอียดชั่วขณะรายวันจะเป็น
รวมโดยใช้ชุดข้อมูลเวคเตอร์เวลาพื้นที่ส่งผลให้แรสเตอร์เวลาพื้นที่เอาต์พุต
ชุด ฝน_agg. วิธีการรวมถูกตั้งค่าเป็น รวม เพื่อสะสม
ค่าหยาดน้ำฟ้าของช่วงเวลาทั้งหมดในชุดข้อมูลเวกเตอร์เวลาอวกาศ การสุ่มตัวอย่าง
ตัวเลือกรับรองว่าเฉพาะแผนที่แรสเตอร์ที่เป็นชั่วคราวในช่วงเวลาของ
ชุดข้อมูลเวกเตอร์กาลอวกาศได้รับการพิจารณาสำหรับการคำนวณ ดังนั้นจึงตั้งค่าตัวเลือกเป็น
ประกอบด้วย (เลเยอร์แผนที่เวกเตอร์ประทับเวลาชั่วคราว บรรจุ เลเยอร์แผนที่แรสเตอร์):
MAPS="map_1 map_2 map_3 map_4 map_5 map_6 map_7"
สำหรับแผนที่ใน ${MAPS} ; ทำ
r.mapcalc expression="${map} = 1"
echo ${map} >> map_list.txt
ทำ
t.create type=strds temporaltype=absolute \
เอาท์พุท=ฝนตก_daily \
title="ปริมาณน้ำฝนรายวัน" \
Description="ทดสอบชุดข้อมูลพร้อมปริมาณน้ำฝนรายวัน"
t.register -i type=raster input=precipitation_daily \
file=map_list.txt start="2012-08-20" increment="1 days"
t.info type=strds อินพุต=precipitation_daily
+-------------------------- ชุดข้อมูล Space Time Raster ------------------------- ----+
| |
+-------------------------- ข้อมูลพื้นฐาน ---------------------------- -----------+
| รหัส: ........................ precipitation_daily@PERMANENT
| ชื่อ: ...................... precipitation_daily
| Mapset: ....................ถาวร
| ผู้สร้าง: .................... soeren
| ประเภทชั่วขณะ: ............. สัมบูรณ์
| เวลาสร้าง: ............ 2014-11-23 16:48:17.686979
| แก้ไขเวลา:.......... 2014-11-23 16:48:18.302978
| ประเภทความหมาย:.............. หมายถึง
+-------------------------- เวลาที่แน่นอน ---------------------------- --------------+
| เริ่มเวลา:.................. 2012-09-10 00:00:00
| เวลาสิ้นสุด:.................... 2012-09-17 00:00:00
| ความละเอียด:................ 1 วัน
| ประเภทชั่วคราวของแผนที่:...... ช่วง
+-------------------------- ขอบเขตเชิงพื้นที่ ---------------------------- --------------+
| เหนือ:................................ 80.0
| ใต้:...................... 0.0
| ตะวันออก:.. ................. 120.0
| ทิศตะวันตก:.......................... 0.0
| ด้านบน:........................ 0.0
| ล่าง:............................ 0.0
+-------------------------- ข้อมูลเมตาดาต้า ---------------------------- --------+
| ตารางการลงทะเบียนแรสเตอร์:...... raster_map_register_3225725979b14b5db343a00835b882c7
| ความละเอียดเหนือ-ใต้ นาที:. 10.0
| ความละเอียดเหนือ-ใต้ สูงสุด:. 10.0
| ความละเอียดตะวันออก-ตะวันตก ขั้นต่ำ:... 10.0
| ความละเอียดสูงสุด East-west:... 10.0
| ค่าต่ำสุดขั้นต่ำ:.......... 1.0
| ค่าต่ำสุดสูงสุด:.......... 1.0
| ค่าสูงสุดต่ำสุด:.......... 1.0
| ค่าสูงสุดสูงสุด:.......... 1.0
| ประเภทการรวม:........... ไม่มี
| จำนวนแผนที่ที่ลงทะเบียน:.. 7
|
| ชื่อ:
| ปริมาณน้ำฝนรายวัน
| คำอธิบาย:
| ทดสอบชุดข้อมูลพร้อมปริมาณน้ำฝนรายวัน
| ประวัติคำสั่ง:
| #2014-11-23 16:48:17น
| t.create type="strds" temporaltype="absolute"
| output="precipitation_daily" title="ปริมาณน้ำฝนรายวัน"
| Description="ทดสอบชุดข้อมูลพร้อมปริมาณน้ำฝนรายวัน"
| #2014-11-23 16:48:18น
| t.register -i type="rast" input="precipitation_daily"
| file="map_list.txt" start="2012-08-20" increment="1 days"
|
+------------------------------------------------- ----------------------------+
v.สุ่มเอาท์พุต=points_1 n=20
v.สุ่มเอาท์พุต=points_2 n=20
t.create type=stvds temporaltype=absolute \
เอาท์พุท=คะแนน \
title="คะแนน" \
Description="คะแนนรวม"
t.register -i type=อินพุตเวกเตอร์=จุด \
map=points_1 start="2012-08-20" increment="3 วัน"
t.register -i type=อินพุตเวกเตอร์=จุด \
map=points_2 start="2012-08-23" increment="4 วัน"
t.info type=stvds อินพุต=จุด
+-------------------------- ชุดข้อมูลเวกเตอร์ Space Time ------------------------ ----+
| |
+-------------------------- ข้อมูลพื้นฐาน ---------------------------- -----------+
| รหัส: ........................ points@PERMANENT
| ชื่อ: ...................... คะแนน
| Mapset: ....................ถาวร
| ผู้สร้าง: .................... soeren
| ประเภทชั่วขณะ: ............. สัมบูรณ์
| เวลาสร้าง: ............ 2014-11-23 16:48:49.193903
| แก้ไขเวลา:.......... 2014-11-23 16:48:50.185671
| ประเภทความหมาย:.............. หมายถึง
+-------------------------- เวลาที่แน่นอน ---------------------------- --------------+
| เริ่มเวลา:.................. 2012-08-20 00:00:00
| เวลาสิ้นสุด:.................... 2012-08-27 00:00:00
| ความละเอียด:................ 1 วัน
| ประเภทชั่วคราวของแผนที่:...... ช่วง
+-------------------------- ขอบเขตเชิงพื้นที่ ---------------------------- --------------+
| เหนือ:................................ 79.283411
| ใต้:...................... 5.724954
| ตะวันออก:.. ................. 118.881168
| ทิศตะวันตก:.......................... 0.016755
| ด้านบน:........................ 0.0
| ล่าง:............................ 0.0
+-------------------------- ข้อมูลเมตาดาต้า ---------------------------- --------+
| Vector register table:...... vector_map_register_6f02d33e0ee243d1a521aaaca39ecb31
| จำนวนคะแนน ........... 40
| จำนวนบรรทัด ............ 0
| จำนวนขอบเขต ....... 0
| จำนวนเซนทรอยด์ ........ 0
| จำนวนใบหน้า ............ 0
| จำนวนเมล็ด .......... 0
| จำนวนพื้นฐาน ....... 40
| จำนวนโหนด ............ 0
| จำนวนพื้นที่ ............ 0
| จำนวนเกาะ .......... 0
| จำนวนหลุม ............ 0
| จำนวนเล่ม .......... 0
| จำนวนแผนที่ที่ลงทะเบียน:.. 2
|
| ชื่อ:
| คะแนน
| คำอธิบาย:
| คะแนนรวม
| ประวัติคำสั่ง:
| #2014-11-23 16:48:49น
| t.create type="stvds" temporaltype="absolute"
| output="points" title="Points" description="คะแนนสำหรับการรวม"
| #2014-11-23 16:48:49น
| t.register -i type="vect" input="points"
| map="points_1" start="2012-08-20" increment="3 วัน"
| #2014-11-23 16:48:50น
| t.register -i type="vect" input="points"
| map="points_2" start="2012-08-23" increment="4 วัน"
|
+------------------------------------------------- ----------------------------+
t.rast.aggregate.ds อินพุต=precipitation_daily \
เอาต์พุต=การตกตะกอน_agg \
ตัวอย่าง=ประเภทคะแนน=stvds \
ชื่อฐาน=prec_agg \
วิธี=สุ่มตัวอย่าง=ประกอบด้วย
t.สนับสนุนอินพุต=precipitation_agg \
title="ปริมาณน้ำฝนรวม" \
Description="ชุดข้อมูลปริมาณน้ำฝนรวม"
t.info type=strds อินพุต=precipitation_agg
+-------------------------- ชุดข้อมูล Space Time Raster ------------------------- ----+
| |
+-------------------------- ข้อมูลพื้นฐาน ---------------------------- -----------+
| รหัส: ........................ precipitation_agg@PERMANENT
| ชื่อ: ...................... ฝน_agg
| Mapset: ....................ถาวร
| ผู้สร้าง: .................... soeren
| ประเภทชั่วขณะ: ............. สัมบูรณ์
| เวลาสร้าง: ............ 2014-11-23 16:53:23.488799
| แก้ไขเวลา:.......... 2014-11-23 16:53:28.714886
| ประเภทความหมาย:.............. หมายถึง
+-------------------------- เวลาที่แน่นอน ---------------------------- --------------+
| เริ่มเวลา:.................. 2012-08-20 00:00:00
| เวลาสิ้นสุด:.................... 2012-08-27 00:00:00
| ความละเอียด:................ 1 วัน
| ประเภทชั่วคราวของแผนที่:...... ช่วง
+-------------------------- ขอบเขตเชิงพื้นที่ ---------------------------- --------------+
| เหนือ:................................ 80.0
| ใต้:...................... 0.0
| ตะวันออก:.. ................. 120.0
| ทิศตะวันตก:.......................... 0.0
| ด้านบน:........................ 0.0
| ล่าง:............................ 0.0
+-------------------------- ข้อมูลเมตาดาต้า ---------------------------- --------+
| Raster register table:...... raster_map_register_7b025eb7431747c98c5c1ad971e8c282
| ความละเอียดเหนือ-ใต้ นาที:. 10.0
| ความละเอียดเหนือ-ใต้ สูงสุด:. 10.0
| ความละเอียดตะวันออก-ตะวันตก ขั้นต่ำ:... 10.0
| ความละเอียดสูงสุด East-west:... 10.0
| ค่าต่ำสุดขั้นต่ำ:.......... 3.0
| ค่าต่ำสุดสูงสุด:.......... 4.0
| ค่าสูงสุดต่ำสุด:.......... 3.0
| ค่าสูงสุดสูงสุด:.......... 4.0
| ประเภทของการรวม:........... sum
| จำนวนแผนที่ที่ลงทะเบียน:.. 2
|
| ชื่อ:
| ปริมาณน้ำฝนรวม
| คำอธิบาย:
| ชุดข้อมูลปริมาณน้ำฝนรวม
| ประวัติคำสั่ง:
| #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" การสุ่มตัวอย่าง = "มี"
| #2014-11-23 16:53:28น
| t.support input="precipitation_agg"
| title="ปริมาณน้ำฝนรวม"
| Description="ชุดข้อมูลปริมาณน้ำฝนรวม"
|
+------------------------------------------------- ----------------------------+
โมดิส ดาวเทียม เซ็นเซอร์ ประจำวัน ข้อมูล การรวมตัว ไปยัง 8 วัน
ในตัวอย่างนี้ แสดงการรวมจากข้อมูลรายวันถึงแปดวัน นี้ "แปดวัน
สัปดาห์" ใช้ในผลิตภัณฑ์เซ็นเซอร์ดาวเทียม MODIS บางตัว
# หมายเหตุ: ตัวอย่างเขียนด้วยภาษาเชลล์
# สร้างแผนที่ทุก ๆ 8 วันเป็นแผนที่เมล็ดพันธุ์
สำหรับปีใน `seq 2000 2001` ; ทำ
สำหรับ doy ใน `seq -w 1 8 365` ; ทำ
r.mapcalc -s expression="8day_${year}_${doy} = rand(0.0,40.0)"
ทำ
ทำ
# จากชื่อแต่ละแผนที่เราเอาปีปีและแปลงแล้ว
# ถึงวันที่เริ่มต้นและสิ้นสุด YYYY-MM-DD และสร้างไฟล์ด้วย
# ชื่อแผนที่ วันที่เริ่มต้นและวันที่สิ้นสุด
g.list type=raster pattern=8day_20??_* > names_list
สำหรับ NAME ใน `cat names_list` ; ทำ
#พาร์เซ
YEAR=`echo $NAME | ตัด -d'_' -f2`
DOY=`สะท้อน $NAME | ตัด -d'_' -f3`
# แปลง YYYY_DOY เป็น YYYY-MM-DD
DOY=`สะท้อน "$DOY" | sed 's/^0*//'`
doy_end=0
ถ้า [ $DOY -le "353" ] ; แล้ว
doy_end=$(( $DOY + 8 ))
elif [ $DOY -eq "361" ] ; แล้ว
ถ้า [ $[$YEAR % 4] -eq 0] && [ $[$YEAR % 100] -ne 0] || [ $[$YEAR % 400] -eq 0 ] ; แล้ว
doy_end=$(( $DOY + 6 ))
อื่น
doy_end=$(( $DOY + 5 ))
fi
fi
DATE_START=`date -d "${YEAR}-01-01 +$(( ${DOY} - 1 ))วัน" +%Y-%m-%d`
DATE_END=`date -d "${YEAR}-01-01 +$(( ${doy_end} -1 ))วัน" +%Y-%m-%d`
# ไฟล์ข้อความพร้อมชื่อแผนที่ วันที่เริ่มต้น และวันที่สิ้นสุด
echo "$NAME|$DATE_START|$DATE_END" >> list_map_start_end_time.txt
ทำ
#ตรวจสอบรายชื่อที่สร้าง
รายการแมว_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
#แผนที่ทั้งหมดยกเว้นแผนที่สุดท้ายในแต่ละปีคิดเป็น 8 วัน
#ช่วง. แต่การรวมตัวเริ่มต้นใหม่ทุกๆครั้ง
#วันที่ 1 มกราคม
# สร้าง strds เหมือน MODIS 8 วัน
t.create type=strds temporaltype=absolute \
output=8day_ts title="ชุดเวลา 8 วัน" \
description="STRDS พร้อม MODIS เช่นการรวม 8 วัน"
# ลงทะเบียนแผนที่
t.register type=อินพุตแรสเตอร์=8day_ts \
ไฟล์=list_map_start_end_time.txt
# ตรวจสอบ
t.info อินพุต=8day_ts
t.rast.list อินพุต = 8day_ts
# สุดท้าย คัดลอกการรวมเป็นอนุกรมเวลารายวัน
t.rast.aggregate.ds -s อินพุต=daily_ts ตัวอย่าง=8day_ts \
เอาต์พุต = 8day_agg ชื่อฐาน = 8day_agg \
method=ค่าเฉลี่ยการสุ่มตัวอย่าง=มี
# เพิ่มข้อมูลเมตา
t.สนับสนุนอินพุต=8day_agg \
title="8 วันรวม ts" \
description="ชุดข้อมูลรวมที่เหมือน MODIS 8 วัน"
# ตรวจสอบรายการแผนที่ใน strds รวมที่สร้างขึ้นใหม่
t.rast.list อินพุต=8day_agg
ชื่อ|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
ใช้ t.rast.aggregate.dsgrass ออนไลน์โดยใช้บริการ onworks.net