นี่คือคำสั่ง grdsamplegmt ที่สามารถเรียกใช้ในผู้ให้บริการโฮสติ้งฟรีของ OnWorks โดยใช้หนึ่งในเวิร์กสเตชันออนไลน์ฟรีของเรา เช่น Ubuntu Online, Fedora Online, โปรแกรมจำลองออนไลน์ของ Windows หรือโปรแกรมจำลองออนไลน์ของ MAC OS
โครงการ:
ชื่อ
grdsample - สุ่มตัวอย่างกริดไปยัง lattice ใหม่
เรื่องย่อ
เกรดตัวอย่าง in_grdfile out_grdfile [ เพิ่มขึ้น ] [ ภูมิภาค -ระดับ- -f] [
-n] [ -r ] [ -x-n-
หมายเหตุ ไม่อนุญาตให้มีช่องว่างระหว่างแฟล็กตัวเลือกและอาร์กิวเมนต์ที่เกี่ยวข้อง
DESCRIPTION
เกรดตัวอย่าง อ่านไฟล์กริดและสอดแทรกเพื่อสร้างไฟล์กริดใหม่ด้วย: a
การลงทะเบียนที่แตกต่างกัน (-r or -T); หรือระยะห่างกริดใหม่หรือจำนวนโหนด (-I), และ
อาจเป็นภูมิภาคย่อยใหม่ (-R). bicubic [ค่าเริ่มต้น], bilinear, B-spline หรือ
ใช้การแก้ไขเพื่อนบ้านที่ใกล้ที่สุด ดู -n สำหรับการตั้งค่า โปรดทราบว่าการใช้ -R เท่านั้นคือ
เทียบเท่ากับ ตัด or กริดดิต -S. เกรดตัวอย่าง ได้อย่างปลอดภัยสร้างตาข่ายที่ดีจากหยาบ
หนึ่ง; การสนทนาอาจใช้นามแฝงเว้นแต่ว่าข้อมูลจะถูกกรองโดยใช้ กรี้ดดด or
ตัวกรองกริด.
เมื่อ -R ถูกละเว้น ตารางเอาท์พุตจะครอบคลุมพื้นที่เดียวกันกับกริดอินพุท เมื่อไหร่ -I
ถูกละเว้น ระยะห่างกริดของกริดผลลัพธ์จะเหมือนกับกริดอินพุต ทั้ง
-r or -T สามารถใช้เพื่อเปลี่ยนการลงทะเบียนกริด เมื่อละเว้น ตารางเอาต์พุตจะ
มีการลงทะเบียนเดียวกันกับตารางอินพุต
ที่จำเป็น อาร์กิวเมนต์
in_grdfile
ชื่อของไฟล์ตารางไบนารี 2 มิติอินพุต (ดูรูปแบบไฟล์ GRID ด้านล่าง)
-Gout_grdfile
ชื่อของไฟล์กริดเอาต์พุต (ดูรูปแบบไฟล์ GRID ด้านล่าง)
ตัวเลือก อาร์กิวเมนต์
-Iซินซี[หน่วย-หยิน[หน่วย-
x_inc [และเป็นทางเลือก y_inc] คือระยะห่างกริด หรือ ต่อท้ายคำต่อท้าย
ตัวแก้ไข ตามภูมิศาสตร์ (องศา) พิกัด: ผนวก m เพื่อระบุอาร์คนาทีหรือ s
เพื่อระบุอาร์ควินาที ถ้าหน่วยใดหน่วยหนึ่ง e, f, k, M, n or u ถูกต่อท้าย
แทน การเพิ่มขึ้นจะถือว่าเป็นเมตร ฟุต กม. ไมล์ ทะเล
ไมล์หรือฟุตสำรวจของสหรัฐฯ ตามลำดับ และจะถูกแปลงให้เทียบเท่า
องศาลองจิจูดที่ละติจูดกลางของภูมิภาค (การแปลงขึ้นอยู่กับ
PROJ_ELLIPSOID) ถ้า /y_inc ให้ไว้แต่ตั้งค่าเป็น 0 จะถูกรีเซ็ตเท่ากับ x_inc;
มิฉะนั้นจะถูกแปลงเป็นองศาละติจูด ทั้งหมด พิกัด: ถ้า = is
ต่อท้าย max . ที่สอดคล้องกัน x (ทางทิศตะวันออก) หรือ y (ทางทิศเหนือ) อาจจะปรับเล็กน้อย
เพื่อให้พอดีกับการเพิ่มขึ้นที่กำหนด [โดยค่าเริ่มต้น การเพิ่มขึ้นอาจถูกปรับ
เล็กน้อยเพื่อให้พอดีกับโดเมนที่กำหนด] สุดท้าย แทนที่จะให้เพิ่มขึ้น คุณอาจ
ระบุไฟล์ จำนวน of โหนด ที่ต้องการโดยการต่อท้าย + เป็นจำนวนเต็มที่ให้มา
การโต้แย้ง; การเพิ่มขึ้นจะถูกคำนวณใหม่จากจำนวนโหนดและ
โดเมน. ค่าที่เพิ่มขึ้นที่ได้ขึ้นอยู่กับว่าคุณได้เลือก a . หรือไม่
กริดไลน์ที่ลงทะเบียนหรือกริดที่ลงทะเบียนพิกเซล ดูรูปแบบไฟล์แอปสำหรับรายละเอียด
หมายเหตุ: if -Rไฟล์กริด ใช้แล้วระยะห่างกริดได้รับการเริ่มต้นแล้ว ใช้
-I เพื่อแทนที่ค่า
-NS[หน่วย]เอ็กซ์มิน/xmax/ยมิน/วายแม็กซ์[NS] (มากกว่า ... )
ระบุภูมิภาคที่สนใจ
-T แปลระหว่างการลงทะเบียนกริดและพิกเซล หากอินพุตถูกลงทะเบียนกริด
เอาต์พุตจะถูกลงทะเบียนพิกเซลและในทางกลับกัน
-วี[ระดับ] (มากกว่า ... )
เลือกระดับการใช้คำฟุ่มเฟือย [c]
-ฉ[ฉัน|o]คอลอินโฟ (มากกว่า ... )
ระบุชนิดข้อมูลของคอลัมน์อินพุตและ/หรือเอาต์พุต
-n[b|c|l|n][+a][+bBC][+c][+tธรณีประตู] (มากกว่า ... )
เลือกโหมดการแก้ไขสำหรับกริด
-r (มากกว่า ... )
ตั้งค่าการลงทะเบียนโหนดพิกเซล [เส้นตาราง]
-NS[[-]n] (มากกว่า ... )
จำกัดจำนวนคอร์ที่ใช้ในอัลกอริธึมแบบมัลติเธรด (ต้องใช้ OpenMP)
-^ or เพียงแค่ -
พิมพ์ข้อความสั้น ๆ เกี่ยวกับไวยากรณ์ของคำสั่ง จากนั้นออก (หมายเหตุ: บน Windows
ใช้เพียง -).
-+ or เพียงแค่ +
พิมพ์ข้อความการใช้งาน (ช่วยเหลือ) อย่างกว้างขวาง รวมถึงคำอธิบายของทุก ๆ
ตัวเลือกเฉพาะโมดูล (แต่ไม่ใช่ตัวเลือกทั่วไปของ GMT) จากนั้นออก
-? or ไม่ ข้อโต้แย้ง
พิมพ์ข้อความการใช้งาน (วิธีใช้) ที่สมบูรณ์ รวมทั้งคำอธิบายตัวเลือก จากนั้น
ทางออก
--รุ่น
พิมพ์เวอร์ชัน GMT และออก
--show-datadir
พิมพ์เส้นทางแบบเต็มไปยังไดเรกทอรีแชร์ GMT และออก
GRID VALUES PRECISION
โดยไม่คำนึงถึงความแม่นยำของข้อมูลที่ป้อน โปรแกรม GMT ที่สร้างไฟล์กริดจะ
ภายในถือกริดในอาร์เรย์ทศนิยม 4 ไบต์ ทำเพื่อเก็บความทรงจำ
และยิ่งไปกว่านั้น ส่วนใหญ่ถ้าไม่ใช่ข้อมูลจริงทั้งหมดสามารถจัดเก็บได้โดยใช้ทศนิยม 4 ไบต์
ค่านิยม ข้อมูลที่มีความแม่นยำสูงกว่า (เช่น ค่าความแม่นยำสองเท่า) จะสูญเสียข้อมูลนั้นไป
แม่นยำเมื่อ GMT ทำงานบนกริดหรือเขียนกริดใหม่ เพื่อจำกัดการสูญเสีย
ความแม่นยำในการประมวลผลข้อมูล คุณควรพิจารณาปรับข้อมูลให้เป็นมาตรฐานก่อนเสมอ
การประมวลผล
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 มิติ
ผลกระทบที่เกิด OF GRID สุ่มตัวอย่างใหม่
การสุ่มตัวอย่างหรือการสุ่มตัวอย่างกริดจะใช้อัลกอริธึมต่างๆ (ดู -n) ที่อาจนำไปสู่
การบิดเบือนที่เป็นไปได้หรือผลลัพธ์ที่ไม่คาดคิดในค่าที่สุ่มตัวอย่าง หนึ่งผลที่คาดหวัง
ของการสุ่มตัวอย่างด้วย splines มีแนวโน้มว่าค่าที่สุ่มตัวอย่างใหม่จะเกินเล็กน้อย
ขีดจำกัดต่ำสุด/สูงสุดทั่วโลกของกริดดั้งเดิม หากสิ่งนี้ไม่เป็นที่ยอมรับ คุณสามารถกำหนด
การตัดค่าที่สุ่มตัวอย่างใหม่เพื่อไม่ให้เกินค่าต่ำสุด/สูงสุดที่ป้อนเข้าโดย
เพิ่ม +c เพื่อคุณ -n ตัวเลือก
ข้อเสนอแนะ
หากจุดแก้ไขไม่ได้อยู่บนโหนดของกริดอินพุต แสดงว่า NaN อยู่ที่โหนดใดๆ ใน
พื้นที่ใกล้เคียงโดยรอบจุดจะให้ค่า NaN ที่สอดแทรก บีคิวบิก
interpolation [default] ให้ผลอนุพันธ์อันดับ XNUMX อย่างต่อเนื่อง แต่ต้องการย่านใกล้เคียงของ
4 โหนดโดย 4 โหนด การแก้ไขแบบไบลิเนียร์ [-n] ใช้เพียงย่าน 2 คูณ 2 แต่
ให้ผลต่อเนื่องเป็นศูนย์เท่านั้น ใช้ bicubic เมื่อความนุ่มนวลเป็นสิ่งสำคัญ ใช้ bilinear
เพื่อลดการแพร่กระจายของ NaNs
ตัวอย่าง
ในการสุ่มตัวอย่างตาราง 5 x 5 นาทีใน hawaii_5by5_topo.nc เป็นตาราง 1 นาที:
gmt grdsample hawaii_5by5_topo.nc -I1m -Ghawaii_1by1_topo.nc
ในการแปลไฟล์ Surface.nc ที่ลงทะเบียน gridline เป็นการลงทะเบียนพิกเซลในขณะที่ยังคงรักษา
ภูมิภาคและช่วงกริดเดียวกัน:
gmt grdsample surface.nc -T -Gpixel.nc
ใช้ grdsamplegmt ออนไลน์โดยใช้บริการ onworks.net