นี่คือคำสั่ง grdtrackgmt ที่สามารถเรียกใช้ในผู้ให้บริการโฮสต์ฟรีของ OnWorks โดยใช้หนึ่งในเวิร์กสเตชันออนไลน์ฟรีของเรา เช่น Ubuntu Online, Fedora Online, โปรแกรมจำลองออนไลน์ของ Windows หรือโปรแกรมจำลองออนไลน์ของ MAC OS
โครงการ:
ชื่อ
grdtrack - ตัวอย่างกริดที่ตำแหน่ง (x,y) ที่ระบุ
เรื่องย่อ
ติดตาม [ ไซไฟล์ ] กริด1 กริด2 - f|p|m|r|R[+l- ความยาว[u]/ds[การเว้นวรรค][+a-ดีไฟล์
] [ เส้น[,เส้น-+aaz][+d][+iInc[u-+lความยาว[u-+nnp][+oaz][+rรัศมี[u-
ภูมิภาค ] [ วิธี/การปรับเปลี่ยน -รัศมี[u-+e|p-ระดับ- -b] [
-d] [ -f] [ -g] [ -h] [ -i] [ -n] [
-o] [ -s] [ -:[i|o-
หมายเหตุ ไม่อนุญาตให้มีช่องว่างระหว่างแฟล็กตัวเลือกและอาร์กิวเมนต์ที่เกี่ยวข้อง
DESCRIPTION
ติดตาม อ่านไฟล์กริดอย่างน้อยหนึ่งไฟล์ (หรือไฟล์ Sandwell/Smith IMG) และตาราง (จาก
ไฟล์หรืออินพุตมาตรฐาน แต่เห็น -E สำหรับข้อยกเว้น) ด้วยตำแหน่ง (x,y) [หรือ (lon,lat)] ใน
สองคอลัมน์แรก (อาจมีคอลัมน์เพิ่มเติม) มันสอดแทรกกริดที่
ตำแหน่งในตารางและเขียนตารางโดยเพิ่มค่า interpolated เป็น (one
หรือมากกว่า) คอลัมน์ใหม่ อีกทางหนึ่ง (-C) ข้อมูลที่ป้อนเข้าถือเป็นส่วนบรรทัดและ
เราสร้าง cross-profiles มุมฉากที่จุดข้อมูลแต่ละจุดหรือด้วยการแยกที่เท่ากัน
และสุ่มตัวอย่างกริดตามโปรไฟล์เหล่านี้ bicubic [ค่าเริ่มต้น], bilinear, B-spline หรือ
เพื่อนบ้านที่ใกล้ที่สุด (ดู -n) มีการใช้การประมาณการโดยต้องมีเงื่อนไขขอบเขตที่
ขอบเขตของภูมิภาค (ดู -n; ค่าเริ่มต้นใช้เงื่อนไข "ธรรมชาติ" (อนุพันธ์ย่อยบางส่วนที่สอง
ปกติถึงขอบเป็นศูนย์) เว้นแต่ระบบจะรับรู้กริดเป็นระยะโดยอัตโนมัติ)
ที่จำเป็น อาร์กิวเมนต์
-Gไฟล์กริด
ไฟล์กริด เป็นไฟล์กริดไบนารี 2 มิติที่มีฟังก์ชัน f(x,y) หากตารางที่กำหนด
อยู่ในรูปแบบ Sandwell/Smith Mercator คุณต้องผนวกรายการ . ที่คั่นด้วยเครื่องหมายจุลภาค
อาร์กิวเมนต์ที่มีมาตราส่วนเพื่อคูณข้อมูล (โดยปกติคือ 1 หรือ 0.1) โหมด
ซึ่งย่อมาจากสิ่งต่อไปนี้: (0) ไฟล์ Img ที่ไม่มีรหัสข้อ จำกัด ส่งคืน data
ทุกจุด (1) ไฟล์ Img ที่มีรหัสข้อจำกัด ส่งคืนข้อมูลทุกจุด (2)
ไฟล์ Img ที่มีการเข้ารหัสข้อจำกัด ส่งคืนข้อมูลเฉพาะที่จุดจำกัดและ NaN
ที่อื่นและ (3) ไฟล์ Img ที่มีการเข้ารหัสข้อ จำกัด ส่งคืน 1 ที่ข้อ จำกัด และ0
ที่อื่น และอาจเลือกละติจูดสูงสุดในไฟล์ IMG [80.738] คุณอาจทำซ้ำ
-G มากที่สุดเท่าที่คุณมีกริดที่คุณต้องการสุ่มตัวอย่าง หรือใช้ -G+ลรายการ
เพื่อส่งรายชื่อไฟล์ กริดจะถูกสุ่มตัวอย่างและผลลัพธ์จะถูกส่งออกใน
คำสั่งที่ได้รับ (ดูรูปแบบไฟล์ GRID ด้านล่าง)
ตัวเลือก อาร์กิวเมนต์
ไซไฟล์ นี่คือ ASCII (หรือเลขฐานสอง ดู -ไบ) ไฟล์ที่ 2 คอลัมน์แรกถือ (x,y)
ตำแหน่งที่ผู้ใช้ต้องการสุ่มตัวอย่างชุดข้อมูล 2 มิติ
-Af|pm|r|ร
สำหรับการสุ่มตัวอย่างแทร็ก (if -C เป็นการตั้งค่า) เราสามารถเลือกวิธีการดำเนินการนี้ได้
ผนวก f เพื่อรักษาจุดเดิม แต่เพิ่มจุดกลางหากจำเป็น [ค่าเริ่มต้น]
m as fแต่ก่อนอื่นตามเส้นเมอริเดียน (ตาม y) แล้วขนานกัน (ตาม x) p as fแต่
อันดับแรก ตามขนาน (ตาม y) ตามด้วยเส้นเมอริเดียน (ตาม x) r เพื่อสุ่มตัวอย่างที่
ตำแหน่งที่เท่ากัน จุดอินพุตไม่จำเป็นต้องรวมอยู่ในเอาต์พุตและ
R as rแต่ปรับระยะห่างที่กำหนดเพื่อให้พอดีกับความยาวของแทร็ก สุดท้ายผนวก
+l ถ้าควรวัดระยะทางตามเส้นรูปสี่เหลี่ยมขนมเปียกปูน (loxodromes) ละเว้นเว้นแต่
-C ถูกนำมาใช้.
-Cความยาว[ยู]/ds[การเว้นวรรค][+ก]
ใช้ส่วนของบรรทัดอินพุตเพื่อสร้างระยะเท่ากันและ (เป็นทางเลือก) ที่มีระยะห่างเท่ากัน
ชุดของโปรไฟล์ข้ามตามที่เราสุ่มตัวอย่างกริด [ค่าเริ่มต้นเพียงตัวอย่าง
ตารางที่ตำแหน่งอินพุต] ระบุสเกลความยาวสองสเกลที่ควบคุมวิธีการ
การสุ่มตัวอย่างเสร็จสิ้น: ความยาว กำหนดความยาวเต็มของแต่ละ cross-profile ในขณะที่ ds
คือระยะการสุ่มตัวอย่างตามแต่ละ cross-profile หรือไม่ก็ได้ ต่อท้าย /การเว้นวรรค เป็นเวลา
ระยะห่างที่เท่ากันระหว่าง cross-profiles [ค่าเริ่มต้นสร้าง cross-profiles ที่
พิกัดอินพุต]. โดยค่าเริ่มต้น ข้ามโปรไฟล์ทั้งหมดมีทิศทางเดียวกัน ผนวก
+a เพื่อสลับทิศทางของ cross-profiles ผนวกหน่วยที่เหมาะสมกับ ความยาว;
มันกำหนดหน่วยที่ใช้สำหรับ ds [และ การเว้นวรรค] (ดูหน่วยด้านล่าง) หน่วยเริ่มต้นสำหรับ
กริดทางภูมิศาสตร์คือเมตรในขณะที่กริดคาร์ทีเซียนหมายถึงหน่วยผู้ใช้ ผลลัพธ์
คอลัมน์จะเป็น lon, ลาดพร้าว, อ., ทิศทางของดาววัดบนพื้นโลก, z1, z2, ... , zn ( zi เป็นตัวอย่าง
ค่าสำหรับแต่ละ n กริด)
-Dดีไฟล์
ในคอนเสิร์ตกับ -C เราสามารถบันทึกบรรทัดต้นฉบับ (อาจสุ่มตัวอย่างใหม่) ลงในไฟล์ได้
ดีไฟล์ [ค่าเริ่มต้นจะบันทึกเฉพาะ cross-profiles] คอลัมน์จะเป็น lon, ลาดพร้าว, อ.,
ทิศทางของดาววัดบนพื้นโลก, z1, z2, ... (ค่าสุ่มตัวอย่างสำหรับแต่ละกริด)
-Eเส้น[,เส้น,...][+กaz][+d][+iInc[u]][+ลความยาว[u]][+นnp][+โอaz][+รรัศมี[หรือ]
แทนที่จะอ่านพิกัดของแทร็กอินพุต ให้ระบุโปรไฟล์ผ่านพิกัดและ
ตัวดัดแปลง รูปแบบของแต่ละคน เส้น is เริ่มต้น/หยุดที่นี่มี เริ่มต้น or หยุด เป็นอย่างใดอย่างหนึ่ง
lon/ลาดพร้าว (x/y สำหรับข้อมูลคาร์ทีเซียน) หรือคีย์ XY 2 อักขระที่ใช้
รูปแบบรูปแบบการให้เหตุผลแบบ "pstext" เพื่อระบุจุดบนแผนที่เป็น
[LCR][BMT]. นอกจากนี้ คุณสามารถใช้ Z-, Z+ เพื่อหมายถึงค่าต่ำสุดและสูงสุดทั่วโลกได้
ตำแหน่งในกริด (ใช้ได้เฉพาะเมื่อให้กริดเดียวเท่านั้น) แทนที่จะเป็นสอง
พิกัดที่คุณสามารถระบุที่มาและหนึ่งใน +a, +o,หรือ +r. คุณสามารถผนวก
+iInc[u] เพื่อกำหนดช่วงเวลาสุ่มตัวอย่าง; ถ้าไม่ได้ให้ เราก็ผิดนัดครึ่งหนึ่ง
ช่วงกริดขั้นต่ำ NS +a ตั้งค่าราบของโปรไฟล์ของความยาวที่กำหนด
เริ่มต้นที่แหล่งกำเนิดที่กำหนดในขณะที่ +o กำหนดโปรไฟล์ไว้ที่จุดเริ่มต้น ทั้งสอง
ต้องการ +l. สำหรับการสุ่มตัวอย่างแบบวงกลม ให้ระบุ +r เพื่อกำหนดวงกลมรัศมีที่กำหนด
เน้นที่แหล่งกำเนิด ตัวเลือกนี้ต้องใช้อย่างใดอย่างหนึ่ง +n or +i. +nnp ตั้งค่า
จำนวนคะแนนที่ต้องการในขณะที่ +lความยาว ให้ความยาวทั้งหมดของโปรไฟล์
ผนวก +d เพื่อแสดงผลระยะทางตามเส้นทางหลังพิกัด หมายเหตุ: ไม่
ไฟล์แทร็กจะถูกอ่าน โปรดทราบว่าคุณสามารถเลือกหน่วยระยะทางได้เพียงหน่วยเดียวเท่านั้น
การให้หน่วยต่าง ๆ จะส่งผลให้เกิดข้อผิดพลาด หากไม่มีการระบุหน่วยเรา
ค่าเริ่มต้นคือระยะทางวงกลมใหญ่ในหน่วยกม. (หากเป็นภูมิศาสตร์) ถ้าทำงานด้วย
ข้อมูลทางภูมิศาสตร์ที่คุณสามารถเพิ่ม - (Flat Earth) หรือ + (Geodesic) ถึง Inc, ความยาว,หรือ
รัศมี เพื่อเปลี่ยนโหมดการคำนวณระยะทาง [Great Circle]
-N Do ไม่ ข้ามจุดที่อยู่นอกโดเมนของกริด [Default only output
จุดภายในโดเมนกริด]
-NS[หน่วย]เอ็กซ์มิน/xmax/ยมิน/วายแม็กซ์[NS] (มากกว่า ... )
ระบุภูมิภาคที่สนใจ
-Sวิธี/การปรับเปลี่ยน
ร่วมกับ -Cคำนวณโปรไฟล์แบบเรียงซ้อนเดียวจากทุกโปรไฟล์ทั่ว
แต่ละส่วน ผนวกวิธีการคำนวณการซ้อน: a = ค่าเฉลี่ย (เฉลี่ย), m =
ค่ามัธยฐาน p = โหมด (โอกาสสูงสุด), l = ต่ำกว่า L = ต่ำกว่าแต่พิจารณาเท่านั้น
ค่าบวก u = บน, U = บน แต่พิจารณาเฉพาะค่าลบ [a]
การปรับเปลี่ยน ควบคุมการส่งออก; เลือกหนึ่งตัวเลือกหรือมากกว่าจากตัวเลือกเหล่านี้: +a : ผนวก
ค่าที่ซ้อนกันให้กับ cross-profiles ทั้งหมด +d : ผนวกส่วนเบี่ยงเบนสแต็กต่อทั้งหมด
ข้ามโปรไฟล์ +d : ผนวกข้อมูลที่เหลือ (data - stack) เข้ากับ cross-profiles ทั้งหมด
+s[ไฟล์] : บันทึกโปรไฟล์ที่ซ้อนกันไปที่ ไฟล์ [grdtrack_stacked_profile.txt] +cความจริง :
คำนวณซองจดหมายบนโปรไฟล์ที่ซ้อนกันเป็น +/- ความจริง *การเบี่ยงเบน [2]. หมายเหตุ: (1)
การเบี่ยงเบนขึ้นอยู่กับ วิธี และเป็น st.dev (a), สเกล L1 (m และ p) หรือครึ่งช่วง
(บน-ล่าง)/2. (2) ไฟล์โปรไฟล์ที่ซ้อนกันประกอบด้วย 1 กลุ่มบวก 4-6
คอลัมน์ หนึ่งกลุ่มสำหรับแต่ละตารางตัวอย่าง คอลัมน์แรกถือระยะข้าม
ในขณะที่ 4 ตัวแรกในกลุ่มเก็บค่าซ้อน ส่วนเบี่ยงเบน ค่าต่ำสุด และ max
มูลค่า. ถ้า วิธี เป็นหนึ่งใน a|m|p จากนั้นเราก็เขียนความมั่นใจล่างและบน
ขอบเขต (ดู +c). เมื่อหนึ่งหรือมากกว่าของ +a, +dและ +r ใช้แล้วเราต่อท้าย
ผลลัพธ์ไปยังจุดสิ้นสุดของแต่ละแถวสำหรับ cross-profiles ทั้งหมด ออเดอร์ซ้อนเสมอ
ค่า (+a) ตามด้วยส่วนเบี่ยงเบน (+d) และเศษเหลือ (+r). เมื่อมีมากกว่าหนึ่ง
ตารางถูกสุ่มตัวอย่าง ลำดับของคอลัมน์ 1-3 นี้จะถูกทำซ้ำสำหรับแต่ละตาราง
-NS[รัศมี[u]][+e|p]
เพื่อใช้กับการสุ่มตัวอย่างกริดแบบปกติ และจำกัดไว้เพียงกริดเดียวที่ไม่ใช่ IMG ถ้า
โหนดที่ใกล้ที่สุดไปยังจุดอินพุตคือ NaN ค้นหาภายนอกจนกว่าเราจะพบ
โหนดที่ไม่ใช่ NaN ที่ใกล้ที่สุดและรายงานค่านั้นแทน เลือกระบุการค้นหา
รัศมีซึ่งจำกัดการพิจารณาให้อยู่ในระยะนี้จากอินพุต
จุด. เพื่อรายงานตำแหน่งของโหนดที่ใกล้ที่สุดและระยะห่างจากอินพุต
จุด ต่อท้าย +e. เพื่อแทนที่จุดป้อนข้อมูลด้วยพิกัดของ
โหนดที่ใกล้ที่สุด ผนวก +p.
-วี[ระดับ] (มากกว่า ... )
เลือกระดับการใช้คำฟุ่มเฟือย [c]
-Z เขียนเฉพาะค่า z ตัวอย่างเท่านั้น [ค่าเริ่มต้นจะเขียนคอลัมน์ทั้งหมด]
-: สลับระหว่างอินพุต/เอาต์พุต (ลองจิจูด ละติจูด) และ (ละติจูด ลองจิจูด)
[ค่าเริ่มต้นคือ (ลองจิจูด ละติจูด)]
-bi[ncols][NS] (มากกว่า ... )
เลือกอินพุตไบนารีดั้งเดิม [ค่าเริ่มต้นคือ 2 คอลัมน์อินพุต]
-โบ[ncols][ชนิด] (มากกว่า ... )
เลือกเอาต์พุตไบนารีดั้งเดิม [ค่าเริ่มต้นคือมากกว่าหนึ่งอินพุต]
-d[ฉัน|o]ไม่มีข้อมูล (มากกว่า ... )
แทนที่คอลัมน์อินพุตที่เท่ากับ ไม่มีข้อมูล ด้วย NaN และทำย้อนกลับในเอาต์พุต
-ฉ[ฉัน|o]คอลอินโฟ (มากกว่า ... )
ระบุชนิดข้อมูลของคอลัมน์อินพุตและ/หรือเอาต์พุต
-g[a]x|y|d|X|Y|D|[Col]z[+|-]ช่องว่าง[หรือ] (มากกว่า ... )
กำหนดช่องว่างข้อมูลและตัวแบ่งบรรทัด
-h[i|o][n][+c][+d][+rข้อสังเกต][+รชื่อเรื่อง] (มากกว่า ... )
ข้ามหรือสร้างเรกคอร์ดส่วนหัว
-iปลอกคอ[ล][สขนาด][โอชดเชย-...] (มากกว่า ... )
เลือกคอลัมน์อินพุต (0 คือคอลัมน์แรก)
-n[b|c|l|n][+a][+bBC][+c][+tธรณีประตู] (มากกว่า ... )
เลือกโหมดการแก้ไขสำหรับกริด
-oปลอกคอ- (มากกว่า ... )
เลือกคอลัมน์ผลลัพธ์ (0 คือคอลัมน์แรก)
-NS[ปลอกคอ][ก|ร] (มากกว่า ... )
ตั้งค่าการจัดการระเบียน NaN
-^ or เพียงแค่ -
พิมพ์ข้อความสั้น ๆ เกี่ยวกับไวยากรณ์ของคำสั่ง จากนั้นออก (หมายเหตุ: บน Windows
ใช้เพียง -).
-+ or เพียงแค่ +
พิมพ์ข้อความการใช้งาน (ช่วยเหลือ) อย่างกว้างขวาง รวมถึงคำอธิบายของทุก ๆ
ตัวเลือกเฉพาะโมดูล (แต่ไม่ใช่ตัวเลือกทั่วไปของ GMT) จากนั้นออก
-? or ไม่ ข้อโต้แย้ง
พิมพ์ข้อความการใช้งาน (วิธีใช้) ที่สมบูรณ์ รวมทั้งคำอธิบายตัวเลือก จากนั้น
ทางออก
--รุ่น
พิมพ์เวอร์ชัน GMT และออก
--show-datadir
พิมพ์เส้นทางแบบเต็มไปยังไดเรกทอรีแชร์ GMT และออก
UNITS
สำหรับหน่วยระยะทางของแผนที่ ให้ต่อท้าย หน่วย d สำหรับองศาอาร์ค m สำหรับอาร์คนาทีและ s สำหรับอาร์ค
ที่สองหรือ e สำหรับมิเตอร์ [ค่าเริ่มต้น] f สำหรับเท้า k สำหรับกม. M สำหรับไมล์กฎเกณฑ์ n สำหรับการเดินเรือ
ไมล์และ u สำหรับเท้าสำรวจของสหรัฐฯ โดยค่าเริ่มต้น เราจะคำนวณระยะทางดังกล่าวโดยใช้ทรงกลม
การประมาณด้วยวงกลมขนาดใหญ่ นำหน้า - เป็นระยะทาง (หรือหน่วยไม่มีระยะทางคือ
กำหนด) เพื่อทำการคำนวณ "Flat Earth" (เร็วกว่า แต่แม่นยำน้อยกว่า) หรือ prepend + ไปยัง
ทำการคำนวณ geodesic ที่แน่นอน (ช้ากว่า แต่แม่นยำกว่า)
ASCII FORMAT PRECISION
รูปแบบเอาต์พุต ASCII ของข้อมูลตัวเลขถูกควบคุมโดยพารามิเตอร์ใน your gmt.conf
ไฟล์. ลองจิจูดและละติจูดถูกจัดรูปแบบตาม FORMAT_GEO_OUT ในขณะที่ other
ค่าต่างๆ จะถูกจัดรูปแบบตาม FORMAT_FLOAT_OUT โปรดทราบว่ารูปแบบที่มีผลสามารถ
ทำให้สูญเสียความแม่นยำในการส่งออก ซึ่งอาจนำไปสู่ปัญหาต่าง ๆ ที่ปลายน้ำ ถ้า
คุณพบว่าผลลัพธ์ไม่ได้เขียนด้วยความแม่นยำเพียงพอ ให้ลองเปลี่ยนเป็นไบนารี
เอาท์พุท (-โบ หากมี) หรือระบุทศนิยมเพิ่มเติมโดยใช้การตั้งค่า FORMAT_FLOAT_OUT
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 หรือต่ำกว่า ธรณีประตู.
ตัวอย่าง
เพื่อสุ่มตัวอย่างไฟล์ hawaii_topo.nc ตามแทร็ก SEASAT track_4.xyg (ตาราง ASCII
ประกอบด้วยเส้นแวง ละติจูด และแรงโน้มถ่วงจาก SEASAT นำหน้าด้วยส่วนหัว
บันทึก):
grdtrack track_4.xyg -Ghawaii_topo.nc -h > track_4.xygt
เพื่อสุ่มตัวอย่างไฟล์รูปแบบ Sandwell/Smith IMG topo.8.2.img (การตรวจวัดความลึก 2 นาที)
บนตาราง Mercator) และตาราง Muller et al age.3.2.nc ตามพิกัด lon,lat
ให้ไว้ในไฟล์ cruise_track.xy ลอง
grdtrack cruise_track.xy -Gtopo.8.2.img,1,1 -Gage.3.2.nc > ความลึก-age.d
เพื่อสุ่มตัวอย่างไฟล์รูปแบบ Sandwell/Smith IMG Grav.18.1.img (เปิดความผิดปกติของอากาศอิสระ 1 นาที
ตาราง Mercator) ตามแนวขวางยาว 100 กม. ซึ่งตั้งฉากกับส่วนของเส้นตรง
ให้ไว้ในไฟล์ track.xy สร้าง cross-profiles ทุก ๆ 25 กม. และสุ่มตัวอย่าง grid
ทุก 3 กม. ลอง
grdtrack track.xy -Ggrav.18.1.img,0.1,1 -C100k/3/25 -Ar > xprofiles.txt
ในการสุ่มตัวอย่าง data.nc กริดตามบรรทัดจากมุมล่างซ้ายไปมุมขวาบน
โดยใช้ระยะห่างกริด 1 กม. และระยะทางเอาต์พุตด้วย ลอง
grdtrack -ELB/RT+i1k+d -Gdata.nc > profiles.txt
ใช้ grdtrackgmt ออนไลน์โดยใช้บริการ onworks.net