นี่คือคำสั่ง r.seriesgrass ที่สามารถเรียกใช้ในผู้ให้บริการโฮสติ้งฟรีของ OnWorks โดยใช้หนึ่งในเวิร์กสเตชันออนไลน์ฟรีของเรา เช่น Ubuntu Online, Fedora Online, โปรแกรมจำลองออนไลน์ของ Windows หรือโปรแกรมจำลองออนไลน์ของ MAC OS
โครงการ:
ชื่อ
r.series - ทำให้ค่าเซลล์เอาต์พุตแต่ละค่าเป็นฟังก์ชันของค่าที่กำหนดให้กับ
เซลล์ที่เกี่ยวข้องในเลเยอร์แผนที่แรสเตอร์อินพุต
KEYWORDS
แรสเตอร์ การรวม ซีรีส์
เรื่องย่อ
r.series
r.series --ช่วยด้วย
r.series [-nz] [อินพุต=พร้อมชื่อ[,พร้อมชื่อ-ไฟล์=พร้อมชื่อ] เอาท์พุต=พร้อมชื่อ[,พร้อมชื่อ-
วิธี=เชือก[,เชือก-ปริมาณ=ลอย[,ลอย-น้ำหนัก=ลอย[,ลอย-
[พิสัย=เฮ้-เขียนทับ-ช่วย-ละเอียด-เงียบสงบ-ui]
ธง:
-n
เผยแพร่ NULLs
-z
อย่าเปิดไฟล์ไว้
--เขียนทับ
อนุญาตให้ไฟล์เอาต์พุตเขียนทับไฟล์ที่มีอยู่
--ช่วยด้วย
พิมพ์สรุปการใช้งาน
--รายละเอียด
เอาต์พุตโมดูล verbose
--เงียบ
เอาต์พุตโมดูลเงียบ
--UI
บังคับให้เปิดใช้กล่องโต้ตอบ GUI
พารามิเตอร์:
อินพุต=ชื่อ[,ชื่อ-
ชื่อของแมปแรสเตอร์อินพุต
ไฟล์=พร้อมชื่อ
ไฟล์อินพุตที่มีชื่อแผนที่แรสเตอร์หนึ่งชื่อและตัวเลือกหนึ่งน้ำหนักต่อบรรทัด ตัวคั่นฟิลด์
ระหว่างชื่อและน้ำหนักคือ |
เอาท์พุต=ชื่อ[,ชื่อ- [ที่จำเป็น]
ชื่อแผนที่แรสเตอร์เอาต์พุต
วิธี=สตริง[,สตริง- [ที่จำเป็น]
การดำเนินการรวม
ตัวเลือก: เฉลี่ย, นับ, ค่ามัธยฐาน โหมด, ขั้นต่ำ min_raster, ขีดสุด, max_raster,
มาตรฐาน พิสัย, ผลรวม ความแปรปรวน ความหลากหลาย, ความลาดชัน ชดเชย เด็คโคออฟ, มูลค่า, ควอร์ต1,
ควอร์ต3, เปอร์เซ็น90, ปริมาณ, ความเบ้ โด่ง
ปริมาณ=ลอยตัว[,ลอยตัว-
ปริมาณที่จะคำนวณสำหรับวิธี = ปริมาณ
ตัวเลือก: 0.0-1.0
น้ำหนัก=ลอยตัว[,ลอยตัว-
ตัวคูณน้ำหนักสำหรับแต่ละแมปอินพุต ค่าเริ่มต้นคือ 1.0 สำหรับแต่ละแมปอินพุต
พิสัย=เฮ้
ละเว้นค่าที่อยู่นอกช่วงนี้
DESCRIPTION
r.series ทำให้ค่าเซลล์เอาต์พุตแต่ละค่าเป็นฟังก์ชันของค่าที่กำหนดให้กับ
เซลล์ที่เกี่ยวข้องในเลเยอร์แผนที่แรสเตอร์อินพุต
มีวิธีการดังต่อไปนี้:
· เฉลี่ย : ค่าเฉลี่ย
· นับ: จำนวนเซลล์ที่ไม่ใช่ NULL
· ค่ามัธยฐาน: ค่ามัธยฐาน
· โหมด: ค่าที่เกิดขึ้นบ่อยที่สุด
· ขั้นต่ำ: ค่าต่ำสุด
· สูงสุด: ค่าสูงสุด
· ช่วง: ช่วงของค่า (สูงสุด - นาที)
· stddev: ส่วนเบี่ยงเบนมาตรฐาน
· ผลรวม: ผลรวมของค่า
· ความแปรปรวน: ความแปรปรวนทางสถิติ
· ความหลากหลาย: จำนวนค่านิยมที่แตกต่างกัน
· ความชัน: ความชันการถดถอยเชิงเส้น
· ออฟเซ็ต: ออฟเซ็ตการถดถอยเชิงเส้น
· detcoeff: สัมประสิทธิ์การถดถอยเชิงเส้นของการกำหนด
· tvalue: การถดถอยเชิงเส้น t-value
· min_raster: หมายเลขแผนที่แรสเตอร์ที่มีค่าอนุกรมเวลาต่ำสุด
· max_raster: หมายเลขแผนที่แรสเตอร์ที่มีค่าอนุกรมเวลาสูงสุด
โปรดทราบว่าพารามิเตอร์ส่วนใหญ่ยอมรับหลายคำตอบ อนุญาตให้รวมหลายค่าได้
คำนวณในการรันครั้งเดียว เช่น:
อินพุต r.series=map1,...,mapN \
เอาต์พุต=map.mean,map.stddev \
วิธี=ค่าเฉลี่ย,stddev
หรือ:
อินพุต r.series=map1,...,mapN \
เอาต์พุต=map.p10,map.p50,map.p90 \
วิธีการ = ควอไทล์, ควอไทล์, ควอไทล์ \
ควอนไทล์=0.1,0.5,0.9
ต้องระบุจำนวนค่าเท่ากันสำหรับตัวเลือกทั้งหมด
หมายเหตุ
ไม่มีข้อมูล (โมฆะ) การจัดการ
ด้วยระบบเส้นทาง -n แฟล็ก เซลล์ใดๆ ที่เซลล์อินพุตใดๆ ที่เกี่ยวข้องเป็น NULL is
ตั้งค่าเป็น NULL โดยอัตโนมัติ (การกระจาย NULL) ฟังก์ชันการรวมไม่ได้ถูกเรียก ดังนั้น
วิธีการทั้งหมดประพฤติในลักษณะนี้เกี่ยวกับ -n ธง.
ไม่มี -n แฟล็ก รายการอินพุตทั้งหมดสำหรับแต่ละเซลล์ (รวมถึง NULL) จะถูกส่งต่อไปยัง
ฟังก์ชันรวม การรวมแต่ละรายการสามารถจัดการข้อมูลได้ตามต้องการ ส่วนใหญ่พวกเขา
เพียงคำนวณผลรวมของค่าที่ไม่ใช่ค่า NULL ให้ผลลัพธ์เป็น NULL ก็ต่อเมื่อทั้งหมด
อินพุตเป็น NULL
วางหลักประกันขั้นต่ำ และ สูงสุด การวิเคราะห์
การขอ min_raster และ max_raster เมธอดสร้างแผนที่ด้วยจำนวนแผนที่แรสเตอร์
ที่เก็บค่าต่ำสุด/สูงสุดของอนุกรมเวลา การนับเริ่มต้นที่ 0 รูปภาพ n
สำหรับแรสเตอร์ตัวแรกและตัวสุดท้ายที่อยู่ใน อินพุต=ตามลำดับ
พิสัย การวิเคราะห์
ถ้า ช่วง= ให้ตัวเลือก ค่าใด ๆ ที่อยู่นอกช่วงนั้นจะถือว่าเป็น
ถ้าพวกเขาเป็นโมฆะ NS พิสัย สามารถตั้งค่าพารามิเตอร์เป็น สูงต่ำ เกณฑ์: ค่าภายนอก
ของช่วงนี้จะถือว่าเป็น NULL (กล่าวคือ จะถูกละเว้นโดยผลรวมส่วนใหญ่ หรือ will
ทำให้ผลลัพธ์เป็น NULL ถ้าให้ -n) NS สูงต่ำ เกณฑ์เป็นจุดลอยตัว
ใช้เลย -inf or INF สำหรับธรณีประตูเดียว (เช่น ช่วง=0,inf เพื่อละเว้นค่าลบหรือ
ช่วง=-inf,-200.4 เพื่อละเว้นค่าที่สูงกว่า -200.4)
วัดเชิงเส้น ถดถอย
การถดถอยเชิงเส้น (ความชัน ออฟเซ็ต สัมประสิทธิ์การกำหนด ค่า t) ถือว่าเท่ากับ
ช่วงเวลา หากข้อมูลมีช่วงเวลาไม่ปกติ NULL raster map สามารถเป็น
แทรกลงในอนุกรมเวลาเพื่อให้ช่วงเวลาเท่ากัน (ดูตัวอย่าง)
ปริมาณ
r.series สามารถคำนวณปริมาณตามอำเภอใจ
หน่วยความจำ การบริโภค
การใช้หน่วยความจำไม่ใช่ปัญหาเนื่องจาก r.series ต้องเก็บหนึ่งแถวจากแต่ละแผนที่ที่a
เวลา
การจัดการ of เปิด ไฟล์ ขีด จำกัด
จำนวนแผนที่แรสเตอร์ที่จะประมวลผลนั้นกำหนดโดยขีดจำกัดของระบบปฏิบัติการ สำหรับ
ตัวอย่างเช่น ขีดจำกัดทั้งแบบแข็งและแบบอ่อนโดยปกติคือ 1024 ขีดจำกัดแบบอ่อนสามารถเปลี่ยนแปลงได้
ด้วยเช่น ulimit -n 1500 (ระบบปฏิบัติการที่ใช้ UNIX) แต่ไม่สูงกว่าฮาร์ด
ขีด จำกัด หากต่ำเกินไปคุณสามารถเพิ่มรายการใน .ในฐานะ superuser
/etc/security/limits.conf
#
your_username ฮาร์ด nofile 1500
สิ่งนี้จะเพิ่มขีดจำกัดฮาร์ดเป็น 1500 ไฟล์ ถูกเตือนว่าต้องการเปิดไฟล์มากขึ้น
แกะ. ดูเพิ่มเติมที่หน้า Wiki คำแนะนำสำหรับการประมวลผลข้อมูลแรสเตอร์ขนาดใหญ่
สำหรับแต่ละแผนที่ สามารถระบุตัวประกอบการถ่วงน้ำหนักได้โดยใช้เครื่องหมาย น้ำหนัก ตัวเลือก. การใช้ตุ้มน้ำหนัก
อาจมีความหมายเมื่อคำนวณผลรวมหรือค่าเฉลี่ยของแผนที่ที่มีขอบเขตเวลาต่างกัน
น้ำหนักเริ่มต้นคือ 1.0 จำนวนน้ำหนักต้องเท่ากันกับจำนวน
ใส่แผนที่และต้องมีลำดับเดียวกัน สามารถระบุน้ำหนักในไฟล์อินพุตได้
ใช้ ไฟล์ ตัวเลือกในการวิเคราะห์แผนที่แรสเตอร์จำนวนมากโดยไม่ต้องกดเปิดไฟล์
จำกัด และขีด จำกัด ขนาดของอาร์กิวเมนต์บรรทัดคำสั่ง การคำนวณช้ากว่า
อินพุต วิธีทางเลือก สำหรับทุกแถวเดี่ยวในแผนที่เอาท์พุต แมปอินพุตทั้งหมดจะเปิดขึ้น
และปิด จำนวน RAM จะเพิ่มขึ้นเป็นเส้นตรงตามจำนวนแผนที่อินพุตที่ระบุ
ตัวเลือกอินพุตและไฟล์จะไม่เกิดร่วมกัน อินพุตเป็นไฟล์ข้อความที่มีบรรทัดใหม่
แยกรายชื่อแผนที่แรสเตอร์และน้ำหนักที่เลือกได้ เป็นตัวคั่นระหว่างชื่อแผนที่
และน้ำหนักตัวอักษร "|" จะต้องใช้
ตัวอย่าง
การใช้ r.series ด้วยสัญลักษณ์แทน:
r.series input="`g.list pattern='insitu_data.*' sep=,`" \
output=insitu_data.stddev วิธี=stddev
หมายเหตุ g.รายการ สคริปต์ยังสนับสนุนนิพจน์ทั่วไปสำหรับการเลือกชื่อแผนที่
การใช้ r.series ด้วยแผนที่แรสเตอร์ NULL (เพื่อพิจารณาอนุกรมเวลา "สมบูรณ์"):
r.mapcalc "จำลอง = null ()"
r.series in=map2001,map2002,dummy,dummy,map2005,map2006,dummy,map2008 \
out=res_slope,res_offset,res_coeff meth=ความชัน,ออฟเซ็ต,detcoeff
ตัวอย่างสำหรับการคำนวณผลรวมหลายรายการในการรันครั้งเดียว (3 ผลรวมผลลัพธ์จาก two
อินพุตแผนที่):
r.series in=one,two out=result_avg,res_slope,result_count meth=sum,ลาด,นับ
ตัวอย่างการใช้ตัวเลือกไฟล์ของ r.series:
cat > input.txt << EOF
map1
map2
map3
EOF
ไฟล์ r.series=input.txt out=result_sum meth=sum
ตัวอย่างการใช้ตัวเลือกไฟล์ของ r.series รวมถึงน้ำหนัก น้ำหนัก 0.75 ควรเป็น
กำหนดให้ map2 เนื่องจากแผนที่อื่นๆ ไม่มีน้ำหนัก เราสามารถเว้นไว้ได้:
cat > input.txt << EOF
map1
แผนที่2|0.75
map3
EOF
ไฟล์ r.series=input.txt out=result_sum meth=sum
ตัวอย่างการนับจำนวนวันที่สูงกว่าอุณหภูมิที่กำหนดโดยใช้ค่าเฉลี่ยรายวัน
แผนที่ ('???' เป็นสัญลักษณ์แทน DOY):
# แนวทางสำหรับระบบที่ใช้เชลล์
r.series input=`g.list rast pattern="temp_2003_???_avg" sep=,` \
เอาต์พุต=temp_2003_days_over_25deg range=25.0,100.0 method=count
# วิธีการในสองขั้นตอน (เช่น สำหรับระบบ Windows)
g.list rast pattern="temp_2003_???_avg" output=mapnames.txt
ไฟล์ r.series=mapnames.txt \
เอาต์พุต=temp_2003_days_over_25deg range=25.0,100.0 method=count
ใช้ r.seriesgrass ออนไลน์โดยใช้บริการ onworks.net
