นี่คือคำสั่ง grdgradientgmt ที่สามารถเรียกใช้ในผู้ให้บริการโฮสติ้งฟรีของ OnWorks โดยใช้หนึ่งในเวิร์กสเตชันออนไลน์ฟรีของเรา เช่น Ubuntu Online, Fedora Online, โปรแกรมจำลองออนไลน์ของ Windows หรือโปรแกรมจำลองออนไลน์ของ MAC OS
โครงการ:
ชื่อ
grdgradient - คำนวณอนุพันธ์ของทิศทางหรือการไล่ระดับสีจาก grid
เรื่องย่อ
ไล่ระดับสี in_grdfile out_grdfile [ อาซิม[/อาซิม2-a][c][o][n-
[ส|พี]อะซิม/อีเลฟ[/ล้อมรอบ/กระจาย/กระจกเงา/แบบมีสไตล์- ธง -e][t][amp-ซิก[/ชดเชย]]
] [ ภูมิภาค ] [ สโลปไฟล์ -ระดับ- -ฉ ] [ -n]
หมายเหตุ ไม่อนุญาตให้มีช่องว่างระหว่างแฟล็กตัวเลือกและอาร์กิวเมนต์ที่เกี่ยวข้อง
DESCRIPTION
ไล่ระดับสี อาจใช้ในการคำนวณอนุพันธ์ของทิศทางในทิศทางที่กำหนด (-A),
หรือทิศทาง (-S) [และขนาด (-D)] ของการไล่ระดับสีเวกเตอร์ของข้อมูล
ค่าโดยประมาณในแถวแรก/สุดท้าย/คอลัมน์ของเอาต์พุตขึ้นอยู่กับเงื่อนไขขอบเขต (ดู
-L).
ที่จำเป็น อาร์กิวเมนต์
in_grdfile
ไฟล์กริด 2 มิติที่ใช้คำนวณอนุพันธ์ของทิศทาง (ดูรูปแบบไฟล์กริด
ด้านล่าง)
-Gout_grdfile
ชื่อของไฟล์กริดเอาต์พุตสำหรับอนุพันธ์ของทิศทาง (ดูรูปแบบไฟล์กริด
ด้านล่าง)
ตัวเลือก อาร์กิวเมนต์
-Aอาซิม[/อาซิม2]
ทิศทาง Azimuthal สำหรับอนุพันธ์ของทิศทาง อาซิม คือมุมใน x,y
ระนาบวัดเป็นองศาบวกตามเข็มนาฬิกาจากทิศเหนือ (ทิศทาง +y) ไปทาง
ทิศตะวันออก (ทิศทาง +x) ค่าลบของอนุพันธ์ทิศทาง
-[dz/dx*บาป(อาซิม) + dz/dy*cos(อาซิม)] พบ; การปฏิเสธให้ค่าบวก
เมื่อความชันของ z(x,y) ตกต่ำใน อาซิม ทิศทางความรู้สึกที่ถูกต้องสำหรับ
การแรเงาแสงของภาพ (ดู ภาพกริด และ กริดวิว) โดยแหล่งกำเนิดแสง
เหนือระนาบ x,y ที่ส่องแสงจาก อาซิม ทิศทาง. ทางเลือก จัดหาสอง
ราบ -Aอาซิม/อาซิม2ซึ่งในกรณีนี้การไล่ระดับสีในแต่ละทิศทางคือ
คำนวณและเก็บขนาดที่ใหญ่กว่าไว้ สิ่งนี้มีประโยชน์สำหรับ
ข้อมูลการส่องสว่างที่มีโครงสร้างเป็นเส้นตรงสองทิศทาง เช่น -A0/270
ส่องสว่างจากทิศเหนือ (บน) และทิศตะวันตก (ซ้าย)
-D[a][c][o][n]
หาทิศทางของการไล่ระดับบวก (ขึ้น-ลาด) ของข้อมูล ไปแทน
หาด้าน (ทิศทางลาดลง) ใช้ -ดา. โดยค่าเริ่มต้น เส้นทางคือ
วัดตามเข็มนาฬิกาจากทิศเหนือ เช่น อาซิม in -A ข้างต้น. ผนวก c ใช้ธรรมดา
มุมคาร์ทีเซียนวัดทวนเข็มนาฬิกาจากทิศทางบวก x (ตะวันออก)
ผนวก o เพื่อรายงานทิศทาง (0-180) แทนทิศทาง (0-360) ผนวก n
เพื่อเพิ่มมุมทั้งหมด 90 องศา (เช่น เพื่อให้เกิดการกระแทกเฉพาะที่ของพื้นผิว )
-E[s|พี]อะซิม/อีเลฟ[/ล้อมรอบ/กระจาย/กระจกเงา/แบบมีสไตล์]
คำนวณความกระจ่างใสของ Lambertian ที่เหมาะสมกับการใช้งาน ภาพกริด และ กริดวิวส่วน
Lambertian Reflection ถือว่าพื้นผิวในอุดมคติสะท้อนแสงทั้งหมดนั้น
กระทบกับพื้นผิวและพื้นผิวจะสว่างเท่ากันจากทุกทิศทางการรับชม อาซิม
และ ยกระดับ คือ แอซิมัทและระดับความสูงของเวกเตอร์แสง อุปทาน ล้อมรอบ
กระจาย กระจกเงา แบบมีสไตล์ ซึ่งเป็นพารามิเตอร์ที่ควบคุมคุณสมบัติการสะท้อนแสง
ของพื้นผิว ค่าเริ่มต้นคือ: 0.55/0.6/0.4/10 เพื่อฝากค่าบางอย่างไว้
ไม่ถูกแตะต้องระบุ = เป็นค่าใหม่ ตัวอย่างเช่น -E60/30/=/0.5 ตั้งค่า อาซิม ยกระดับ
และ กระจาย ถึง 60, 30 และ 0.5 และปล่อยให้พารามิเตอร์การสะท้อนแสงอื่น ๆ
ไม่ถูกแตะต้อง ผนวก s เพื่อใช้อัลกอริธึม Lambertian ที่ง่ายกว่า โปรดทราบว่าด้วยแบบฟอร์มนี้
คุณต้องระบุพารามิเตอร์ราบและระดับความสูงเท่านั้น ผนวก p ที่จะใช้
Peucker การประมาณเชิงเส้นทีละชิ้น (อัลกอริธึมที่ง่ายกว่า แต่เร็วกว่า ในกรณีนี้
อาซิม และ ยกระดับ เดินสายไปที่ 315 และ 45 องศา นี่หมายความว่าแม้ว่าคุณจะ
ให้ค่าอื่นๆ ที่พวกเขาจะถูกละเว้น)
-Lธง เงื่อนไขขอบเขต ธง อาจจะ x or y or xy แสดงว่าข้อมูลเป็นระยะในช่วงของ
x หรือ y หรือทั้งสองอย่าง หรือ ธง อาจจะ g ระบุสภาพทางภูมิศาสตร์ (x และ y คือ
lon และ lat) [ค่าเริ่มต้นใช้เงื่อนไข "ธรรมชาติ" (อนุพันธ์ย่อยบางส่วนปกติ
ถึงขอบเป็นศูนย์)]
-สุทธิ][amp-ซิก[/ชดเชย]]
การทำให้เป็นมาตรฐาน [ค่าเริ่มต้น: ไม่มีการทำให้เป็นมาตรฐาน] การไล่ระดับสีที่แท้จริง g ถูกชดเชยและ
ปรับขนาดเพื่อสร้างการไล่ระดับสีปกติ gn ด้วยขนาดเอาต์พุตสูงสุดของ amp.
If amp ไม่ได้รับ, ค่าเริ่มต้น amp = 1. ถ้า ชดเชย ไม่ได้ถูกกำหนดไว้เป็น
ค่าเฉลี่ยของ g. -N อัตราผลตอบแทน gn = amp * (g - ชดเชย)/สูงสุด(หน้าท้อง(g - ชดเชย)) -เน่
ทำให้เป็นมาตรฐานโดยใช้การแจกแจง Laplace สะสม gn = amp * (1.0 -
ประสบการณ์(ตร.ว(2) * (g - ชดเชย)/ ซิก)) ที่ไหน ซิก ถูกประมาณโดยใช้บรรทัดฐาน L1 ของ
(g - ชดเชย) ถ้าไม่ได้รับ -นท ทำให้เป็นมาตรฐานโดยใช้ Cauchy . สะสม
การกระจายผลผลิต gn = (2 * amp / PI) * แทน( (g - ชดเชย)/ ซิก) ที่ไหน ซิก
ถูกประมาณโดยใช้บรรทัดฐาน L2 ของ (g - ชดเชย) ถ้าไม่ได้รับ
-NS[หน่วย]เอ็กซ์มิน/xmax/ยมิน/วายแม็กซ์[NS] (มากกว่า ... )
ระบุภูมิภาคที่สนใจ ใช้ -R ตัวเลือกจะเลือกส่วนย่อยของ
in_grdfile ตาราง หากส่วนย่อยนี้เกินขอบเขตของกริด เฉพาะ
ภูมิภาคทั่วไปจะถูกดึงออกมา
-Sสโลปไฟล์
ชื่อของไฟล์กริดเอาต์พุตที่มีขนาดสเกลาร์ของเวกเตอร์การไล่ระดับสี ต้องใช้ -D
แต่ทำให้ -G ไม่จำเป็น.
-วี[ระดับ] (มากกว่า ... )
เลือกระดับการใช้คำฟุ่มเฟือย [c]
-ฉ กริดทางภูมิศาสตร์ (ขนาดลองจิจูด ละติจูด) จะถูกแปลงเป็นเมตร
ผ่านการประมาณ "Flat Earth" โดยใช้พารามิเตอร์ทรงรีปัจจุบัน
-n[b|c|l|n][+a][+bBC][+c][+tธรณีประตู] (มากกว่า ... )
เลือกโหมดการแก้ไขสำหรับกริด
-^ or เพียงแค่ -
พิมพ์ข้อความสั้น ๆ เกี่ยวกับไวยากรณ์ของคำสั่ง จากนั้นออก (หมายเหตุ: บน Windows
ใช้เพียง -).
-+ or เพียงแค่ +
พิมพ์ข้อความการใช้งาน (ช่วยเหลือ) อย่างกว้างขวาง รวมถึงคำอธิบายของทุก ๆ
ตัวเลือกเฉพาะโมดูล (แต่ไม่ใช่ตัวเลือกทั่วไปของ GMT) จากนั้นออก
-? or ไม่ ข้อโต้แย้ง
พิมพ์ข้อความการใช้งาน (วิธีใช้) ที่สมบูรณ์ รวมทั้งคำอธิบายตัวเลือก จากนั้น
ทางออก
--รุ่น
พิมพ์เวอร์ชัน GMT และออก
--show-datadir
พิมพ์เส้นทางแบบเต็มไปยังไดเรกทอรีแชร์ GMT และออก
GRID ระยะทาง UNITS
หากกริดไม่มีมิเตอร์เป็นหน่วยแนวนอน ให้ผนวก +uหน่วย ไปยังไฟล์อินพุต
ชื่อที่จะแปลงจากหน่วยที่ระบุเป็นเมตร หากกริดของคุณเป็นพื้นที่ทางภูมิศาสตร์ ให้แปลง
ระยะทางเป็นเมตรโดยการจัดหา -ฉ แทน.
ข้อเสนอแนะ
ถ้าไม่รู้ว่าอะไร -N ตัวเลือกที่จะใช้เพื่อสร้างไฟล์ความเข้มสำหรับ ภาพกริด or
กริดวิว, ครั้งแรกที่ดีคือ -เน่0.6.
โดยปกติแล้ว 255 เฉดก็เพียงพอแล้วสำหรับการสร้างภาพข้อมูล คุณสามารถบันทึกดิสก์ได้ 75%
เว้นวรรคโดยผนวก =nb/a ต่อท้ายชื่อไฟล์เอาต์พุต out_grdfile.
หากคุณต้องการสร้างแผนที่ส่องสว่างหลายส่วนของภูมิภาคย่อยของชุดข้อมูลขนาดใหญ่ และคุณ
ต้องการเอฟเฟกต์การส่องสว่างให้สอดคล้องกันในทุกแผนที่ ให้ใช้ -N ตัวเลือกและ
ให้ค่าเท่ากันของ ซิก และ ชดเชย ไปยัง ไล่ระดับสี สำหรับแต่ละแผนที่ การเดาที่ดีคือ
ชดเชย = 0 และ ซิก พบโดย ข้อมูลกริด -L2 or -L1 นำไปใช้กับ grd การไล่ระดับสีที่ไม่ปกติ
หากคุณต้องการเพียงแค่ x- หรือ y-อนุพันธ์ของกริด ใช้ กริด.
GRID ไฟล์ รูปแบบ
โดยค่าเริ่มต้น GMT จะเขียนกริดออกเป็นความแม่นยำเดียวใน netCDF . ร้องเรียน COARDS
รูปแบบไฟล์. อย่างไรก็ตาม GMT สามารถสร้างไฟล์กริดในกริดอื่น ๆ ที่ใช้กันทั่วไปได้
รูปแบบไฟล์และยังอำนวยความสะดวกที่เรียกว่า "การบรรจุ" ของกริดการเขียนจุดลอยตัว
ข้อมูลเป็นจำนวนเต็ม 1 หรือ 2 ไบต์ ในการระบุความแม่นยำ มาตราส่วน และออฟเซ็ต ผู้ใช้ควร
เพิ่มคำต่อท้าย =id[/ขนาด/ชดเชย[/น่าน]], ที่ไหน id เป็นตัวระบุสองตัวอักษรของกริด
ชนิดและความแม่นยำ และ ขนาด และ ชดเชย เป็นสเกลแฟกเตอร์เสริมและออฟเซ็ตเป็น
นำไปใช้กับค่ากริดทั้งหมดและ น่าน คือค่าที่ใช้ระบุข้อมูลที่ขาดหายไป เผื่อ
ตัวละครทั้งสอง id ไม่ได้ระบุไว้เช่นใน =/ขนาด กว่า a id=nf จะถือว่า เมื่อไหร่
อ่านตาราง รูปแบบโดยทั่วไปจะรู้จักโดยอัตโนมัติ ถ้าไม่ใช่ คำต่อท้ายเดียวกัน
สามารถเพิ่มลงในชื่อไฟล์กริดอินพุต ดู แปลงกริด และรูปแบบไฟล์กริดมาตราของ
ข้อมูลอ้างอิงทางเทคนิค GMT และตำราอาหารสำหรับข้อมูลเพิ่มเติม
เมื่ออ่านไฟล์ netCDF ที่มีหลายกริด GMT จะอ่านตามค่าเริ่มต้น
ตาราง 2 มิติแรกที่หาได้ในไฟล์นั้น เพื่อเกลี้ยกล่อม GMT ให้อ่านอีก
ตัวแปรหลายมิติในไฟล์กริด ผนวก ?นามสกุล ไปที่ชื่อไฟล์โดยที่
นามสกุล เป็นชื่อของตัวแปร โปรดทราบว่าคุณอาจต้องหนีความหมายพิเศษ
of ? ในโปรแกรมเชลล์ของคุณโดยใส่แบ็กสแลชไว้ข้างหน้า หรือโดยการวาง
ชื่อไฟล์และคำต่อท้ายระหว่างเครื่องหมายคำพูดหรือเครื่องหมายคำพูดคู่ NS ?นามสกุล ต่อท้ายก็ใช้ได้นะ
สำหรับกริดเอาต์พุตเพื่อระบุชื่อตัวแปรที่แตกต่างจากค่าเริ่มต้น: "z" ดู
แปลงกริด และส่วนดัดแปลงสำหรับ CF และรูปแบบไฟล์กริดของ GMT Technical
อ้างอิงและตำราอาหารสำหรับข้อมูลเพิ่มเติม โดยเฉพาะอย่างยิ่งในการอ่าน splices ของ 3-,
กริด 4 หรือ 5 มิติ
ตัวอย่าง
เพื่อสร้างไฟล์สำหรับการส่องสว่างข้อมูลใน geoid.nc โดยใช้การไล่ระดับสีแบบปกติใน
ช่วง [-0.6,0.6] เลียนแบบแหล่งกำเนิดแสงในทิศทางเหนือและตะวันตก:
gmt grdgradient geoid.nc -A0/270 -Ggradients.nc=nb/a -Ne0.6 -V
วิธีค้นหาแนวราบของผ้าพื้นทะเลในไฟล์ topo.nc:
gmt grdgradient topo.nc -Dno -Gazimuths.nc -V
ข้อมูลอ้างอิง
Horn, BKP, Hill-Shading and the Reflectance Map, การดำเนินการของ IEEE, Vol. 69 หมายเลข
1 มกราคม 1981 น. 14-47. (http://people.csail.mit.edu/bkph/papers/Hill-Shading.pdf)
ใช้ grdgradientgmt ออนไลน์โดยใช้บริการ onworks.net