นี่คือคำสั่ง xyz2grdgmt ที่สามารถเรียกใช้ในผู้ให้บริการโฮสติ้งฟรีของ OnWorks โดยใช้หนึ่งในเวิร์กสเตชันออนไลน์ฟรีของเรา เช่น Ubuntu Online, Fedora Online, โปรแกรมจำลองออนไลน์ของ Windows หรือโปรแกรมจำลองออนไลน์ของ MAC OS
โครงการ:
ชื่อ
xyz2grd - แปลงตารางข้อมูลเป็นไฟล์กริด
เรื่องย่อ
xyz2grd [ ตาราง ] ไฟล์กริด เพิ่มขึ้น ภูมิภาค [[f|l|m|n|r|s|u|z-
xname/ชื่อ/ชื่อ/ขนาด/ชดเชย/โมฆะ/ชื่อเรื่อง/ข้อสังเกต -zfile.zfile-ระดับ-ธง-
[ -ไบ] [ -du] [ -f] [ -h] [ -i] [ -r ] [
-:[i|o-
หมายเหตุ ไม่อนุญาตให้มีช่องว่างระหว่างแฟล็กตัวเลือกและอาร์กิวเมนต์ที่เกี่ยวข้อง
DESCRIPTION
xyz2grd อ่านตาราง z หรือ xyz หนึ่งตารางขึ้นไป และสร้างไฟล์กริดไบนารี xyz2grd จะ
รายงานว่าบางโหนดไม่ได้กรอกข้อมูล โหนดที่ไม่มีข้อ จำกัด ดังกล่าวถูกตั้งค่า
เป็นค่าที่ระบุโดยผู้ใช้ [ค่าเริ่มต้นคือ NaN] โหนดที่มีค่ามากกว่าหนึ่งค่าจะเป็น
กำหนดให้เป็นค่ากลาง เป็นตัวเลือก (โดยใช้ -Z) ตาราง z 1 คอลัมน์อาจถูกอ่านโดยสมมติ
มีโหนดทั้งหมด (z-tables สามารถจัดระเบียบได้หลายรูปแบบ ดู -Z ด้านล่าง)
ที่จำเป็น อาร์กิวเมนต์
-Gไฟล์กริด
ไฟล์กริด เป็นชื่อของไฟล์กริดเอาต์พุตไบนารี (ดูรูปแบบไฟล์ 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] (มากกว่า ... )
ระบุภูมิภาคที่สนใจ
ตัวเลือก อาร์กิวเมนต์
ตาราง ASCII หนึ่งตัวขึ้นไป [หรือเลขฐานสอง ดู -ไบ] ไฟล์ที่มีค่า z หรือ (x,y,z) เดอะ xyz
แฝดสามไม่จำเป็นต้องเรียงลำดับ ต้องจัดเรียงตาราง z หนึ่งคอลัมน์และ -Z
ต้องตั้งค่า
-A[f|l|m|n|r|s|u|z]
โดยค่าเริ่มต้น เราจะคำนวณค่าเฉลี่ยหากรายการหลายรายการอยู่ในโหนดเดียวกัน
ใช้ -A เพื่อเปลี่ยนพฤติกรรมนี้ เว้นแต่จะถูกละเว้น if -Z จะได้รับ ผนวก f or s
เพื่อเก็บจุดข้อมูลแรกหรือจุดสุดท้ายที่กำหนดให้กับแต่ละโหนด ผนวก
l or u เพื่อหาค่าต่ำสุด (ต่ำสุด) หรือสูงสุด (สูงสุด) ในแต่ละโหนด
ตามลำดับ ผนวก m or r เพื่อคำนวณค่าเฉลี่ยหรือค่า RMS ที่แต่ละโหนด
ตามลำดับ ผนวก n เพียงนับจำนวนจุดข้อมูลที่ได้รับมอบหมาย
ไปยังแต่ละโหนด (ต้องการเพียงสองคอลัมน์อินพุต x และ y as z ไม่ปรึกษา)
ผนวก z เพื่อรวมค่าหลายค่าที่เป็นของโหนดเดียวกัน
-Dxname/ชื่อ/ชื่อ/ขนาด/ชดเชย/โมฆะ/ชื่อเรื่อง/ข้อสังเกต
ให้ค่าสำหรับ xname, ชื่อ, ชื่อ (ให้ชื่อของตัวแปรเหล่านั้นและใน
วงเล็บเหลี่ยม หน่วยของพวกมัน เช่น "ระยะทาง [km]") ขนาด (เพื่อคูณค่ากริด
หลังจากอ่าน [ปกติ 1]) ชดเชย (เพื่อเพิ่มลงในกริดหลังจากปรับขนาด [ปกติ 0])
โมฆะ (ค่าที่ใช้แทนข้อมูลที่ขาดหายไป [NaN]), ชื่อเรื่อง (อะไรก็ได้ที่คุณชอบ) และ
ข้อสังเกต (อะไรก็ได้ที่คุณชอบ). หากต้องการปล่อยให้ค่าเหล่านี้ไม่ถูกแตะต้อง ให้ปล่อยฟิลด์
ว่างเปล่า. ฟิลด์ที่ว่างเปล่าในตอนท้ายอาจถูกข้ามไป อีกทางหนึ่ง เพื่อให้ "/" เป็น
ส่วนหนึ่งของค่าใดค่าหนึ่ง ให้ใช้อักขระที่ไม่ใช่ตัวเลขและตัวอักษร (และไม่เท่ากับ
sign) เป็นตัวคั่นทั้งที่ขึ้นต้นและลงท้ายด้วย ตัวอย่างเช่น:
-D:xname:ชื่อ:ชื่อ:ขนาด:ชดเชย:โมฆะ:ชื่อเรื่อง:ข้อสังเกต: ใช้คำพูดเพื่อจัดกลุ่มข้อความ
ที่มีมากกว่าหนึ่งคำ โปรดทราบว่าสำหรับกริดทางภูมิศาสตร์ (-ฉ) xname และ ชื่อ เป็น
ตั้งค่าโดยอัตโนมัติ
-NS[zfile.zfile]
สลับลำดับไบต์ของอินพุตเท่านั้น ไม่มีการสร้างไฟล์กริด ยังต้อง
จัดหา -Z ตัวเลือก. ผลลัพธ์ถูกเขียนถึง zfile.zfile (หรือ stdout ถ้าไม่ได้ให้มา)
-วี[ระดับ] (มากกว่า ... )
เลือกระดับการใช้คำฟุ่มเฟือย [c]
-Z[ธง]
อ่านตาราง ASCII [หรือไบนารี] 1 คอลัมน์ นี่ถือว่าโหนดทั้งหมดเป็น
นำเสนอและจัดเรียงตามแบบแผนการจัดลำดับที่ระบุไว้ใน ธง.
หากข้อมูลที่เข้ามาแทนแถว ให้สร้าง ธง เริ่มกับ T(op) ถ้าแถวแรกคือ y =
ymax หรือ B(ottom) ถ้าแถวแรกคือ y = ymin จากนั้นผนวก L or R เพื่อระบุว่า
องค์ประกอบแรกอยู่ที่ปลายแถวด้านซ้ายหรือด้านขวา ในทำนองเดียวกันสำหรับรูปแบบคอลัมน์: start
กับ L or R เพื่อวางตำแหน่งคอลัมน์แรกแล้วผนวก T or B ขึ้นตำแหน่งก่อน
องค์ประกอบในแถว หมายเหตุ: ตัวบ่งชี้สองแถว/คอลัมน์นี้จำเป็นสำหรับ .เท่านั้น
กริด; สำหรับตารางอื่น ๆ ที่พวกเขาไม่ได้ใช้ สำหรับกริดไลน์ที่ลงทะเบียนแล้ว: ถ้า data
เป็นระยะใน x แต่ข้อมูลที่เข้ามาไม่มีคอลัมน์ (ซ้ำซ้อน) ที่ x
= xmax ต่อท้าย x. สำหรับข้อมูลเป็นระยะใน y โดยไม่มีแถวซ้ำซ้อนที่ y = ymax ให้ผนวก
y. ผนวก sn ข้ามไปก่อน n จำนวนไบต์ (อาจเป็นส่วนหัว) ถ้า
byte-order หรือคำที่ต้องการจะสลับ, ผนวก w. เลือกหนึ่งในหลายข้อมูล
ประเภท (ไบนารีทั้งหมดยกเว้น a):
A การแสดง ASCII ของค่าจุดลอยตัวตั้งแต่หนึ่งค่าขึ้นไปต่อบันทึก
a การแสดง ASCII ของรายการเดียวต่อบันทึก
c int8_t เซ็นชื่ออักขระ 1 ไบต์
u uint8_t อักขระ 1 ไบต์ที่ไม่ได้ลงชื่อ
h int16_t, ลงนามจำนวนเต็ม 2 ไบต์
H uint16_t, จำนวนเต็ม 2 ไบต์ที่ไม่ได้ลงชื่อ
i int32_t, ลงนามจำนวนเต็ม 4 ไบต์
I uint32_t, จำนวนเต็ม 4 ไบต์ที่ไม่ได้ลงชื่อ
l int64_t, จำนวนเต็มยาว (8 ไบต์)
L uint64_t, จำนวนเต็มยาว (8 ไบต์) ที่ไม่ได้ลงนาม
f ความแม่นยำจุดทศนิยม 4 ไบต์เดียว
d ความแม่นยำสองเท่าของทศนิยม 8 ไบต์
รูปแบบเริ่มต้นคือการวางแนวเส้นสแกนของตัวเลข ASCII: -ZTLa. สังเกตได้ว่า -Z เพียง
ใช้กับอินพุต 1 คอลัมน์ ความแตกต่างระหว่าง A และ a คือหลังสามารถ
ถอดรหัสทั้งคู่ ข้อมูลTนาฬิกา และ วว:mm:ss[.xx] รูปแบบในขณะที่อดีตเป็นอย่างเคร่งครัดสำหรับ
ค่าจุดลอยตัวปกติ
-bi[ncols][NS] (มากกว่า ... )
เลือกอินพุตไบนารีดั้งเดิม [ค่าเริ่มต้นคือ 3 คอลัมน์อินพุต] ตัวเลือกนี้ใช้ได้เท่านั้น
ถึงไฟล์อินพุต xyz; ดู -Z สำหรับตาราง z
-duไม่มีข้อมูล (มากกว่า ... )
แทนที่คอลัมน์อินพุตที่เท่ากับ ไม่มีข้อมูล กับแนน ยังตั้งค่าโหนดโดยไม่มีอินพุต xyz
แฝดสามเป็นค่านี้ [ค่าเริ่มต้นคือ NaN]
-ฉ[ฉัน|o]คอลอินโฟ (มากกว่า ... )
ระบุชนิดข้อมูลของคอลัมน์อินพุตและ/หรือเอาต์พุต
-h[i|o][n][+c][+d][+rข้อสังเกต][+รชื่อเรื่อง] (มากกว่า ... )
ข้ามหรือสร้างเรกคอร์ดส่วนหัว ไม่ใช้กับข้อมูลไบนารี
-iปลอกคอ[ล][สขนาด][โอชดเชย-...] (มากกว่า ... )
เลือกคอลัมน์อินพุต (0 คือคอลัมน์แรก)
-r (มากกว่า ... )
ตั้งค่าการลงทะเบียนโหนดพิกเซล [เส้นตาราง]
-:[ผม|o] (มากกว่า ... )
สลับคอลัมน์ที่ 1 และ 2 บนอินพุตและ/หรือเอาต์พุต
-^ 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 เป็นตัวระบุสองตัวอักษรของกริด
ชนิดและความแม่นยำ และ ขนาด และ ชดเชย เป็นสเกลแฟกเตอร์เสริมและออฟเซ็ตเป็น
นำไปใช้กับค่ากริดทั้งหมดและ น่าน คือค่าที่ใช้ระบุข้อมูลที่ขาดหายไป ดู
แปลงกริด และรูปแบบไฟล์กริดมาตราของการอ้างอิงทางเทคนิค GMT และตำราสำหรับ
ข้อมูลมากกว่านี้.
เมื่อเขียนไฟล์ netCDF ตารางจะถูกจัดเก็บตามค่าเริ่มต้นโดยใช้ชื่อตัวแปร "z" ถึง
ระบุชื่อตัวแปรอื่น นามสกุล, ผนวก ?นามสกุล ไปที่ชื่อไฟล์ โปรดทราบว่าคุณอาจ
ต้องหนีความหมายพิเศษของ ? ในโปรแกรมเชลล์ของคุณโดยใส่แบ็กสแลชใน
ข้างหน้า หรือโดยการวางชื่อไฟล์และคำต่อท้ายระหว่างเครื่องหมายคำพูดหรือเครื่องหมายคำพูดคู่
ภูมิศาสตร์ AND เวลา พิกัด
เมื่อประเภทกริดเอาท์พุตเป็น netCDF พิกัดจะมีป้ายกำกับว่า "ลองจิจูด"
"ละติจูด" หรือ "เวลา" ตามคุณสมบัติของข้อมูลที่ป้อนหรือกริด (ถ้ามี) หรือตาม
-f or -R ตัวเลือก. ตัวอย่างเช่น ทั้งสอง -f0x -f1t และ -R90w/90e/0t/3t จะส่งผลให้ a
เส้นลองจิจูด/เวลา เมื่อพิกัด x, y หรือ z เป็นเวลา มันจะถูกเก็บไว้ในกริด
เป็นเวลาสัมพัทธ์ตั้งแต่ยุคตามที่ระบุโดย TIME_UNIT และ TIME_EPOCH ใน gmt.conf ไฟล์
หรือบนบรรทัดคำสั่ง นอกจากนี้ หน่วย คุณลักษณะของตัวแปรเวลาจะระบุ
ทั้งหน่วยนี้และยุค
การแลกเปลี่ยน ข้อ จำกัด
สามารถอ่านข้อมูลทุกประเภท แม้กระทั่งจำนวนเต็ม 64 บิต แต่กริดภายในจะถูกจัดเก็บโดยใช้
ลอย ดังนั้นค่าจำนวนเต็มที่เกิน mantissa 23 บิตของประเภท float อาจไม่
เป็นตัวแทนอย่างแน่นอน เมื่อไหร่ -S ถูกใช้ไม่มีกริดโดยนัยและเราอ่านข้อมูลเป็น an
ภาชนะคู่กลาง ซึ่งหมายความว่าสามารถแสดงจำนวนเต็ม 64 บิตได้ทั้งหมดโดยใช้
mantissa 53 บิตของประเภทสองเท่า
ตัวอย่าง
ในการสร้างไฟล์กริดจากข้อมูล ASCII ใน hawaii_grv.xyz ให้ใช้
GMT xyz2grd hawaii_grv.xyz -Ddegree/degree/mGal/1/0//"ฮาวาย แรงโน้มถ่วง"/"GRS-80
ทรงรี ใช้แล้ว"
-Ghawaii_grv_new.nc -R198/208/18/25 -I5m -V
เพื่อสร้างไฟล์กริดจากไบนารีดิบ (3 คอลัมน์, สแกนไลน์แบบแม่นยำเดียว)
data raw.b ใช้
gmt xyz2grd raw.b -Dm/m/m/1/0 -Graw.nc -R0/100/0/100 -I1 -V -Z -bi3f
เพื่อสร้างไฟล์กริดจากไบนารี่ดิบ USGS DEM (ข้อมูลเชิง scanline เชิงเลขจำนวนเต็มสั้น
topo30.b บน NGDC global relief Data CD-ROM โดยมีค่า -9999 บ่งชี้ว่าหายไป
data หนึ่งต้องในเครื่องบางเครื่องกลับลำดับไบต์ บนเครื่องดังกล่าว (เช่น Sun) ให้ใช้
gmt xyz2grd topo30.b -Dm/m/m/1/0 -Gustopo.nc -R234/294/24/50 -I30s -di-9999 -ZTLhw
สมมติว่าคุณได้รับไฟล์ไบนารีที่มีทศนิยม 4 ไบต์ที่เขียนบน a
เครื่องที่มีลำดับไบต์ที่แตกต่างจากของคุณ คุณสามารถสลับลำดับไบต์ด้วย
gmt xyz2grd float.bin -Snew_floats.bin -V -Zf
ใช้ xyz2grdgmt ออนไลน์โดยใช้บริการ onworks.net