นี่คือคำสั่ง img2grdgmt ที่สามารถเรียกใช้ในผู้ให้บริการโฮสต์ฟรีของ OnWorks โดยใช้หนึ่งในเวิร์กสเตชันออนไลน์ฟรีของเรา เช่น Ubuntu Online, Fedora Online, โปรแกรมจำลองออนไลน์ของ Windows หรือโปรแกรมจำลองออนไลน์ของ MAC OS
โครงการ:
ชื่อ
img2grd - แยกส่วนย่อยของไฟล์ img ในรูปแบบ Mercator หรือ Geographic
เรื่องย่อ
img2grd img ไฟล์ ไฟล์กริด ภูมิภาค ชนิด -มินแลต/แม็กซ์แลต- นาที - นำทาง
-ขนาด-ระดับ- แม็กลอน ] [ -n]
หมายเหตุ ไม่อนุญาตให้มีช่องว่างระหว่างแฟล็กตัวเลือกและอาร์กิวเมนต์ที่เกี่ยวข้อง
DESCRIPTION
img2grd อ่านไฟล์รูปแบบ img แยกชุดย่อย และเขียนลงในไฟล์กริด ดิ -M
ตัวเลือกกำหนดว่าการฉายภาพ Spherical Mercator ของไฟล์ img เป็น .หรือไม่
รักษาไว้หรือถ้าควรเขียนตารางทางภูมิศาสตร์โดยยกเลิกการฉายภาพ Mercator ถ้า
เลือกตารางทางภูมิศาสตร์แล้วคุณยังสามารถขอการสุ่มตัวอย่างซ้ำได้อย่างแม่นยำ -R ให้
ที่จำเป็น อาร์กิวเมนต์
img ไฟล์
ไฟล์รูปแบบ Mercator img เช่น แรงโน้มถ่วงของทะเลหรือเขตข้อมูลภูมิประเทศของพื้นทะเล
ประมาณจากข้อมูลเครื่องวัดระยะสูงจากดาวเทียมโดย Sandwell and Smith หากผู้ใช้ตั้งค่า
ตัวแปรสภาพแวดล้อม $GMT_DATADIRแล้ว img2grd จะพยายามหา img ไฟล์ in
$GMT_DATADIR; ไม่งั้นจะพยายามเปิด img ไฟล์ โดยตรง.
-Gไฟล์กริด
ไฟล์กริด เป็นชื่อของไฟล์กริดเอาต์พุต
-NS[หน่วย]ทิศตะวันตก/ทางทิศตะวันออก/ใต้/ทางทิศเหนือ[/ซมิน/ซีแม็กซ์][NS]
ทิศตะวันตก, ทางทิศตะวันออก, ใต้และ ทางทิศเหนือ ระบุภูมิภาคที่สนใจและคุณอาจระบุ
เป็นองศาทศนิยมหรือในรูปแบบ [+-]dd:mm[:ss.xxx][W|E|S|N] ผนวก r ถ้าต่ำกว่า
พิกัดแผนที่ซ้ายและขวาบนจะได้รับแทน w/e/s/n ทั้งสอง
ชวเลข -ร และ -ถ ย่อมาจาก Global Domain (0/360 และ -180/+180 ในลองจิจูด)
ตามลำดับ โดยมีค่าละติจูด -90/+90) อีกทางหนึ่งสำหรับการสร้างกริด ให้
Rโค้ดลอน/ลาดพร้าว/nx/nyที่นี่มี รหัส เป็นการรวมกันของ 2 อักขระของ L, C, R (สำหรับซ้าย,
ตรงกลางหรือขวา) และ T, M, B สำหรับด้านบน ตรงกลาง หรือด้านล่าง เช่น BL สำหรับล่างซ้าย
นี่แสดงว่าจุดใดบนพื้นที่สี่เหลี่ยมที่ lon/ลาดพร้าว พิกัดอ้างอิง
ถึงและขนาดกริด nx และ ny ด้วยระยะห่างกริดผ่าน -I ใช้ในการสร้าง
ภูมิภาคที่เกี่ยวข้อง หรือระบุชื่อไฟล์กริดที่มีอยู่
และ -R การตั้งค่า (และระยะห่างกริด ถ้ามี) จะถูกคัดลอกจากกริด
การใช้ -Rหน่วย คาดว่าพิกัดที่คาดการณ์ไว้ (คาร์ทีเซียน) เข้ากันได้กับที่เลือก -J
และเราผกผันเพื่อกำหนดพื้นที่ทางภูมิศาสตร์รูปสี่เหลี่ยมผืนผ้าที่แท้จริง สำหรับ
มุมมองมุมมอง (-p) เลือกต่อท้าย /ซมิน/ซีแม็กซ์. ในกรณีของมุมมองมุมมอง
(-p) ช่วง z (ซมิน, ซีแม็กซ์) สามารถต่อท้ายเพื่อระบุมิติที่ XNUMX ได้ นี้
ต้องทำเมื่อใช้ตัว .เท่านั้น -เจซ ตัวเลือกไม่ใช่เมื่อใช้เฉพาะตัว -p ตัวเลือก
ในกรณีหลังนี้ จะมีการพล็อตมุมมองเปอร์สเปคทีฟของเครื่องบินโดยไม่มีส่วนที่สาม
มิติ.
ตัวเลือก อาร์กิวเมนต์
-C ตั้งค่าพิกัด x และ y Mercator ที่สัมพันธ์กับศูนย์ฉายภาพ [ค่าเริ่มต้นคือ
เทียบกับมุมล่างซ้ายของตาราง] ต้องใช้ -M.
-NS[มินแลต/แม็กซ์แลต]
ใช้ช่วงละติจูดที่ขยาย -80.738/+80.738 หรือ ผนวก
มินแลต/แม็กซ์แลต เป็นขอบเขตละติจูดของไฟล์ img อินพุต [ค่าเริ่มต้นคือ
-72.006/72.006]. ปกติไม่จำเป็นเพราะเราสามารถกำหนดขอบเขตจาก
การตรวจสอบขนาดไฟล์
-E สามารถใช้ได้เมื่อ -M ไม่ได้กำหนดให้บังคับให้กริดสุดท้ายมีเหมือนกันทุกประการ
ภูมิภาคตามที่ร้องขอด้วย -R. โดยค่าเริ่มต้น ขอบเขตสุดท้ายคือการประมาณการโดยตรงของ
ภูมิภาค Mercator ดั้งเดิมและโดยทั่วไปจะขยายเกิน .เล็กน้อย
ร้องขอช่วงละติจูดและนอกจากนี้ การเพิ่มกริดในละติจูดไม่ได้
ตรงกับการเพิ่มลองจิจูด อย่างไรก็ตาม การสุ่มตัวอย่างเพิ่มเติมแนะนำขนาดเล็ก
ข้อผิดพลาดในการแก้ไขและควรใช้เฉพาะเมื่อกริดเอาต์พุตต้องตรงกับ
ร้องขอภูมิภาคและมี x_inc = y_inc ในกรณีนี้ภูมิภาคที่กำหนดโดย -R จะต้องเป็น
กำหนดเป็นทวีคูณของการเพิ่มขึ้น (.เช่น -R0/45/45/72)
-I ระบุ นาที เป็นความกว้างของพิกเซล img อินพุตเป็นนาทีของลองจิจูด
[ค่าเริ่มต้นคือ 2.0] ปกติไม่จำเป็นเพราะเราสามารถกำหนดขนาดพิกเซลได้จาก
การตรวจสอบขนาด
-M ส่งออกตาราง Spherical Mercator [ค่าเริ่มต้นคือกริด lon/lat ทางภูมิศาสตร์] ดิ
การฉายภาพ Mercator แบบ Spherical ของไฟล์ img จะถูกรักษาไว้ เพื่อให้ภูมิภาค -R
กำหนดโดยผู้ใช้มีการปรับเปลี่ยนเล็กน้อย ขอบเขตที่แก้ไขสอดคล้องกับขอบ
ของพิกเซล [หรือกลุ่มของ นำทาง พิกเซล]. ส่วนหัวของไฟล์กริดถูกตั้งค่าเพื่อให้ x และ
ความยาวแกน y แทนระยะห่างจากขอบด้านตะวันตกและด้านใต้ของภาพ
วัดในหน่วยเริ่มต้นของผู้ใช้ด้วย -เจม1 และที่ปรับแล้ว -R. โดยการตั้งค่า
ผิดนัด โปรเจ_ วงรี = Sphere ผู้ใช้สามารถสร้างภาพซ้อนทับได้ด้วยการปรับ -R
เพื่อให้ตรงกัน ดู ตัวอย่าง ด้านล่าง. ที่ปรับแล้ว -R ถูกเขียนในตารางด้วย
ข้อสังเกตของส่วนหัวจึงสามารถพบได้ในภายหลัง ดู -C เพื่อกำหนดพิกัดสัมพันธ์กับ
ศูนย์ฉายภาพ
-Nนำทาง เฉลี่ยค่าในพิกเซล img อินพุตเป็น นำทาง by นำทาง สี่เหลี่ยมและสร้าง
หนึ่งพิกเซลเอาต์พุตสำหรับแต่ละตารางดังกล่าว ถ้าใช้กับ -T3 มันจะรายงานค่าเฉลี่ย
ข้อจำกัดระหว่าง 0 ถึง 1 หากใช้กับ -T2 ผลลัพธ์จะเป็นค่าข้อมูลเฉลี่ย
หรือ NaN โดยขึ้นอยู่กับว่าข้อจำกัดเฉลี่ยคือ > 0.5 นำทาง ต้องแบ่งเท่าๆกัน
ลงในขนาดของ imgfile ในหน่วยพิกเซล [ค่าเริ่มต้น 1 ไม่มีค่าเฉลี่ย].
-NS[ขนาด]
คูณค่าไฟล์ img ด้วย ขนาด ก่อนจัดเก็บไว้ในไฟล์กริด [ค่าเริ่มต้นคือ
1.0]. สำหรับไฟล์ img ล่าสุด: ไฟล์ img topo จะถูกเก็บไว้ใน (แก้ไข) เมตร [-S1];
ไฟล์แรงโน้มถ่วงของอากาศอิสระใน mGal*10 [-S0.1 เพื่อรับ mGal]; ไฟล์โก่งแนวตั้งใน
ไมโครเรเดียน*10 [-S0.1 เพื่อรับไมโครเรเดียน] ไฟล์การไล่ระดับแรงโน้มถ่วงแนวตั้งใน
ออทวอส*50 [-S0.02 เพื่อรับ Eotvos หรือ -S0.002 เพื่อรับ mGal/km]). ถ้าไม่ ขนาด ได้รับ
เราพยายามกำหนดมาตราส่วนโดยตรวจสอบชื่อไฟล์เพื่อหาเบาะแส
-Tชนิด ชนิด จัดการการเข้ารหัสข้อมูลข้อจำกัด ชนิด = 0 แสดงว่าไม่มี
ข้อมูลดังกล่าวถูกเข้ารหัสในไฟล์ img (ใช้สำหรับเวอร์ชันก่อนปี 1995 ของ
ข้อมูลแรงโน้มถ่วง) และรับข้อมูลทั้งหมด ชนิด > 0 แสดงว่าข้อมูลข้อจำกัดคือ
เข้ารหัส (เวอร์ชันปัจจุบัน) (1995 และใหม่กว่า (ปัจจุบัน) ของไฟล์ img) เพื่อให้สามารถ
สร้างไฟล์กริดดังนี้: -T1 รับค่าข้อมูลทุกจุด -T2 รับข้อมูล
ค่าที่จุดจำกัดและ NaN ที่จุดที่มีการสอดแทรก -T3 รับ 1 ที่
จุดจำกัดและ 0 ที่จุดที่มีการสอดแทรก [ค่าเริ่มต้นคือ 1]
-วี[ระดับ] (มากกว่า ... )
เลือกระดับการใช้คำฟุ่มเฟือย [c] แนะนำเป็นพิเศษที่นี่เพราะจะเป็นประโยชน์ในการดู
วิธีการปรับพิกัด
-Wแม็กลอน
ระบุ แม็กลอน เป็นขอบเขตลองจิจูดสูงสุดของไฟล์ img อินพุต รุ่น
ตั้งแต่ปี 1995 มี แม็กลอน = 360.0 ในขณะที่บางไฟล์ก่อนหน้านี้มี แม็กลอน = 390.0
[ค่าเริ่มต้นคือ 360.0]
-^ or เพียงแค่ -
พิมพ์ข้อความสั้น ๆ เกี่ยวกับไวยากรณ์ของคำสั่ง จากนั้นออก (หมายเหตุ: บน Windows
ใช้เพียง -).
-+ or เพียงแค่ +
พิมพ์ข้อความการใช้งาน (ช่วยเหลือ) อย่างกว้างขวาง รวมถึงคำอธิบายของทุก ๆ
ตัวเลือกเฉพาะโมดูล (แต่ไม่ใช่ตัวเลือกทั่วไปของ GMT) จากนั้นออก
-? or ไม่ ข้อโต้แย้ง
พิมพ์ข้อความการใช้งาน (วิธีใช้) ที่สมบูรณ์ รวมทั้งคำอธิบายตัวเลือก จากนั้น
ทางออก
--รุ่น
พิมพ์เวอร์ชัน GMT และออก
--show-datadir
พิมพ์เส้นทางแบบเต็มไปยังไดเรกทอรีแชร์ GMT และออก
ภูมิศาสตร์ ตัวอย่าง
พื้นที่ -M ควรยกเว้นตัวเลือกหากคุณต้องการให้ตารางผลลัพธ์อยู่ในพื้นที่ทางภูมิศาสตร์
พิกัด. เพื่อดึงข้อมูลในภูมิภาค -R-40/40/-70/-30 จาก world_grav.img.7.2 และ
ปฏิเสธการให้พิกัดทางภูมิศาสตร์คุณสามารถลอง
img2grd world_grav.img.16.1 -Gmerc_grav.nc -R-40/40/-70/-30 -V
เนื่องจากระยะห่างละติจูดในไฟล์ img นั้นเท่ากันในหน่วย Mercator ดังนั้น
ตารางผลลัพธ์จะไม่ตรงกับที่ระบุ -R อย่างแน่นอนและระยะห่างละติจูดจะไม่
เท่ากับระยะห่างลองจิจูด หากคุณต้องการการจับคู่แบบตรงทั้งหมดกับของคุณ -R และระยะห่างเท่ากัน
ในลองจิจูดและละติจูด ให้ใช้เครื่องหมาย -E ตัวเลือก:
img2grd world_grav.img.16.1 -Gmerc_grav.nc -R-40/40/-70/-30 -E -V
เมอร์คาเตอร์ ตัวอย่าง
เนื่องจากไฟล์ img อยู่ในภาพฉายของ Mercator คุณจึงไม่ควรแยกตารางทางภูมิศาสตร์
ถ้าแผนของคุณคือสร้างแผนที่ Mercator หากคุณทำอย่างนั้นแสดงว่าคุณสิ้นสุดการฉายภาพและ
การฉายภาพซ้ำของตารางทำให้สูญเสียรายละเอียดความยาวคลื่นสั้น ใช้ดีกว่า -M และพล็อตตาราง
ใช้การฉายภาพเชิงเส้นที่มีมาตราส่วนเดียวกันกับการฉายภาพ Mercator ที่ต้องการ (ดู GMT
ตัวอย่างที่ 29) เพื่อดึงข้อมูลในภูมิภาค -R-40/40/-70/-30 จาก world_grav.img.7.2, วิ่ง
gmt img2grd -M world_grav.img.7.2 -Gmerc_grav.nc -R-40/40/-70/-30 -V
โปรดทราบว่า -V ตัวเลือกบอกเราว่าช่วงถูกปรับเป็น
-R-40/40/-70.0004681551/-29.9945810754. เรายังใช้ ข้อมูลกริด จะพบว่าไฟล์กริด
ส่วนหัวแสดงภูมิภาคที่จะเป็น -R0/80/0/67.9666667 นี่คือเรนจ์ของ x,y เราจะได้
จากการฉายภาพ Spherical Mercator โดยใช้ -R-40/40/-70.0004681551/-29.9945810754 และ
-เจม1. ดังนั้น ในการนำ ship.lonlatgrav และใช้เพื่อสุ่มตัวอย่าง merc_grav.nc เราสามารถทำได้:
gmt ชุด PROJ_ELLIPSOID Sphere
gmt mapproject -R-40/40/-70.0004681551/-29.9945810754 -Jm1i ship.lonlatgrav | แผนที่ \
gmt grdtrack -Gmerc_grav.nc | โครงการแผนที่ gmt \
-R-40/40/-70.0004681551/-29.9945810754 -Jm1i -I > ship.lonlatgravsat
ขอแนะนำให้ใช้วิธีการข้างต้นในการฉายภาพและยกเลิกการฉายข้อมูลในลักษณะดังกล่าว
แอปพลิเคชันเพราะจากนั้นมีขั้นตอนการแก้ไขเพียงขั้นตอนเดียว (in ติดตาม). ถ้าอย่างใดอย่างหนึ่ง
ขั้นแรกให้พยายามแปลงไฟล์กริดเป็น lon,lat แล้วสุ่มตัวอย่าง มีสอง
ขั้นตอนการแก้ไข (ในการแปลงและการสุ่มตัวอย่าง)
ในการสร้าง lon,lat grid จากกริดด้านบน เราสามารถใช้
gmt grdproject merc_grav.nc -R-40/40/-70.0004681551/-29.9945810754 -Jm1i -I -D2m -Ggrav.nc
ในบางกรณีสิ่งนี้จะไม่ง่ายเหมือน -R ในระบบพิกัดทั้งสองอาจไม่สอดคล้องกัน
ดี. เมื่อสิ่งนี้เกิดขึ้น เรายังใช้ได้ (จริงๆ แล้วใช้จะดีกว่าเสมอ)
gmt grd2xyz merc_grav.nc | โครงการแผนที่ gmt \
-R-40/40/-70.0004681551/-29.994581075 -Jm1i -I | \
พื้นผิว gmt -R-40/40/-70/70 -I2m -Ggrav.nc
ในการสร้างแผนที่ Mercator ของภูมิภาคข้างต้น สมมติว่าค่า gmt.conf ของเราสำหรับ
PROJ_LENGTH_UNIT คือนิ้ว จากนั้นเนื่องจากไฟล์ merc_grav.nc ด้านบนถูกฉายด้วย -เจม1i
มีความกว้าง 80 นิ้ว เราสามารถสร้างแผนที่กว้าง 8 นิ้วโดยใช้ -เจเอ็กซ์0.1i บนโปรแกรมแผนที่ใด ๆ
นำไปใช้กับกริดนี้ (เช่น grdcontour, ภาพกริด, กริดวิว) แล้วสำหรับโอเวอร์เลย์ซึ่ง
ทำงานใน lon,lat (เช่น ใจสั่น, ปสก) เราสามารถใช้การปรับปรุงข้างต้นได้ -R และ -เจม0.1 ที่จะได้รับ
ทั้งสองระบบให้ตรงกัน
อย่างไรก็ตาม เราสามารถฉลาดกว่านี้ได้ โดยตระหนักว่าไฟล์ img อินพุตมีพิกเซล 2.0
นาทีกว้าง (หรือตรวจสอบ nx และ ny ด้วย grdinfo merc_grav.nc) เราตระหนักดีว่า
merc_grav.nc ใช้ความละเอียดสูงสุดของไฟล์ img และมี 2400 x 2039 พิกเซลและ
ที่ความกว้าง 8 นิ้ว นี่คือ 300 พิกเซลต่อนิ้ว เราตัดสินใจว่าเราไม่ต้องการอะไรมากมายและเรา
จะพอใจกับ 100 พิกเซลต่อนิ้ว ดังนั้นเราจึงต้องการเฉลี่ยข้อมูลเป็น 3 คูณ 3
สี่เหลี่ยม (ถ้าเราต้องการพล็อตเส้นชั้นความสูง เราก็อาจจะเลือกค่าเฉลี่ยของข้อมูลมากกว่านี้
(เช่น 6 คูณ 6) เพื่อให้ได้รูปทรงที่เรียบ) เนื่องจากปี 2039 หารด้วย 3 ไม่ลงตัว เราจึงจะได้
ปรับที่แตกต่างกัน -R เวลานี้:
gmt img2grd -M world_grav.img.7.2 -Gmerc_grav_2.nc -R-40/40/-70/-30 -N3 -V
คราวนี้เราพบว่าพื้นที่ที่ปรับแล้วคือ -R-40/40/-70.023256525/-29.9368261101 และ
เอาต์พุตคือ 800 x 601 พิกเซล ซึ่งเป็นขนาดที่ดีกว่าสำหรับเรา ตอนนี้เราสามารถประดิษฐ์ของเทียมได้
ไฟล์การส่องสว่างสำหรับสิ่งนี้โดยใช้ ไล่ระดับสี:
gmt grdgradient merc_grav_2.nc -Gillum.nc -A0/270 -Ne0.6
และถ้าเรามีไฟล์ CPT ชื่อ "grav.cpt" ด้วย เราก็สามารถสร้างแผนที่บรรเทาด้วยเฉดสีได้
อย่างนี้:
gmt grdimage merc_grav_2.nc -Iillum.nc -Cgrav.cpt -Jx0.1i -K > map.ps
gmt psbasemap -R-40/40/-70.023256525/-29.9368261101 -Jm0.1i -Ba10 -O >> map.ps
สมมติว่าคุณต้องการรับเฉพาะค่าข้อมูลที่จำกัดจากไฟล์ img ใน lat/lon
พิกัด. แล้ววิ่ง img2grd กับ -T2 ตัวเลือก ใช้ grd2xyz เพื่อดัมพ์ค่า pipe
ผ่าน grep -v NaN เพื่อกำจัด NaN และไปป์ผ่าน โครงการแผนที่ กับผกผัน
การฉายภาพดังกล่าวข้างต้น
ใช้ img2grdgmt ออนไลน์โดยใช้บริการ onworks.net