mapprojectgmt - ออนไลน์ในคลาวด์

โครงการ:

ชื่อ

mapproject - ทำการแปลงแผนที่ไปข้างหน้าและผกผัน การแปลง Datum และ geodesy

เรื่องย่อ

โครงการแผนที่ [ ตาราง ] พารามิเตอร์ ภูมิภาค [ b|B|f|F|o|O[Lon0/ละติจูด0-dx/dy- c|i|p ] [
[datums-หน่วย-x0/y0-+|-]หน่วย][+|-- ไลน์.xy-+|-]หน่วย][+-
[a|c|g|m-d|e -h]ราคาเริ่มต้นที่[/ไปยัง-ระดับ-w|h- -b] [
-d] [ -f] [ -g] [ -h] [ -i] [ -o] [
-p] [ -s] [ -:[i|o-

หมายเหตุ ไม่อนุญาตให้มีช่องว่างระหว่างแฟล็กตัวเลือกและอาร์กิวเมนต์ที่เกี่ยวข้อง

DESCRIPTION

โครงการแผนที่ อ่าน (ลองจิจูด, ละติจูด) ตำแหน่งจาก ตาราง [หรืออินพุตมาตรฐาน] และ
คำนวณ (x,y) พิกัดโดยใช้การฉายแผนที่และมาตราส่วนที่ระบุ เป็นทางเลือก มัน
สามารถอ่านตำแหน่ง (x,y) และคำนวณ (ลองจิจูด, ละติจูด) ค่าผกผัน
การเปลี่ยนแปลง สามารถใช้ในการแปลงจุดเชิงเส้น (x,y) ที่ได้จากการแปลงเป็นดิจิทัล
แผนที่การฉายภาพที่ทราบพิกัดทางภูมิศาสตร์ อาจคำนวณระยะทางด้วย
ติดตามไปยังจุดคงที่หรือเข้าใกล้เส้นมากที่สุด สุดท้ายก็ใช้ทำ .ได้
การแปลง Datum ต่างๆ อนุญาตให้ช่องข้อมูลเพิ่มเติมหลังจาก 2 คอลัมน์แรก
ซึ่งต้องมี (ลองจิจูด ละติจูด) หรือ (x,y) ดูตัวเลือก -: เกี่ยวกับวิธีการอ่าน
(ละติจูด,ลองจิจูด) ไฟล์

ที่จำเป็น อาร์กิวเมนต์

-Jพารามิเตอร์ (มากกว่า ... )
เลือกการฉายแผนที่

-NS[หน่วย]เอ็กซ์มิน/xmax/ยมิน/วายแม็กซ์[NS] (มากกว่า ... )
ระบุภูมิภาคที่สนใจ กรณีพิเศษสำหรับการฉาย UTM: If -C ถูกนำมาใช้
และ -R ไม่ได้รับ จากนั้นภูมิภาคจะถูกตั้งค่าให้ตรงกับโซน UTM ที่กำหนด ดังนั้น
เพื่อรักษาสารละลายทรงรีให้สมบูรณ์ (ดู ข้อจำกัด สำหรับ more
ข้อมูล).

ตัวเลือก อาร์กิวเมนต์

ตาราง ASCII หนึ่งรายการขึ้นไป (หรือเลขฐานสอง ดู -ไบ[ncols][ชนิด]) ไฟล์ตารางข้อมูลที่มี a
จำนวนคอลัมน์ข้อมูล หากไม่มีตารางให้ เราก็อ่านจากอินพุตมาตรฐาน

-Ab|B|f|F|o|O[Lon0/ละติจูด0]
-อัฟ คำนวณ (ไปข้างหน้า) ราบจากจุดคงที่ ลอน/ลาดพร้าว ไปยังจุดข้อมูลแต่ละจุด
ใช้ - อ เพื่อรับ back-azimuth จากจุดข้อมูลไปยังจุดคงที่ ใช้ -อ่าว จะได้รับการ
การวางแนว (-90/90) แทนที่จะเป็นแนวราบ (0/360) ตัวพิมพ์ใหญ่ F, B or O จะ
แปลงจากละติจูด geodetic เป็น geocentric และประมาณการ azimuth ของ geodesics
(สมมติว่าทรงรีปัจจุบันไม่ใช่ทรงกลม) หากไม่มีการกำหนดจุดคงที่
เราคำนวณ azimuth (หรือ back-azimuth) จากจุดก่อนหน้า

-ค[dx/dy]
ตั้งค่าศูนย์กลางของพิกัดที่ฉายให้อยู่ที่ศูนย์กลางการฉายแผนที่ [ค่าเริ่มต้นคือ
มุมล่างซ้าย]. ทางเลือก เพิ่มออฟเซ็ตในยูนิตที่คาดการณ์ไว้ที่จะเพิ่ม (หรือ
ลบเมื่อ -I ถูกกำหนด) เป็น (จาก) พิกัดที่คาดการณ์ไว้ เช่น false
ทิศตะวันออกและทิศเหนือสำหรับโซนฉายภาพเฉพาะ [0/0] หน่วยที่ใช้สำหรับ
ออฟเซ็ตคือหน่วยระยะการพล็อตที่มีผลใช้งาน (ดู PROJ_LENGTH_UNIT) เว้นแต่ -F is
ใช้ ซึ่งในกรณีนี้ ออฟเซ็ตเป็นเมตร

-Dc|i|หน้า
แทนที่ PROJ_LENGTH_UNIT ชั่วคราวและใช้ c (ซม.) i (นิ้ว) หรือ p (คะแนน)
แทนที่. ใช้กับ .ไม่ได้ -F.

-อี[datums]
แปลงจาก geodetic (lon, lat, height) เป็น Earth Centered Earth Fixed (ECEF)
(x,y,z) พิกัด (add -I สำหรับการแปลงผกผัน) ต่อท้าย Datum ID (ดู -Qd)
หรือให้ ทรงรี:dx,dy,dz ที่ไหน ทรงรี อาจเป็น ID ทรงรี (ดู -คิว) หรือ
ให้เป็น a[,*inv_f*] โดยที่ a เป็นกึ่งแกนเอกและ inv_f เป็นตัวผกผัน
แบน (0 ถ้าละเว้น) ถ้า datums คือ - หรือไม่ให้ถือว่า WGS-84

-NS[หน่วย]
บังคับ 1:1 มาตราส่วน เช่น เอาต์พุต (หรืออินพุต ดู -I) ข้อมูลอยู่ในการคาดการณ์จริง
เมตร หากต้องการระบุหน่วยอื่นๆ ให้ผนวกหน่วยที่ต้องการ (ดู UNITS) ปราศจาก -F,
เอาต์พุต (หรืออินพุต ดู -I) อยู่ในหน่วยที่กำหนดโดย PROJ_LENGTH_UNIT (แต่
เห็น -D).

-NS[x0/y0-หน่วย-
คำนวณระยะทางตามเส้นทาง or ไปยังจุดเสริมที่ตั้งไว้ด้วย -Gx0/y0. ผนวก
หน่วยระยะทาง (ดูหน่วย) รวมทั้ง c (ระยะทางคาร์ทีเซียนโดยใช้อินพุต
พิกัด) หรือ C (ระยะทางคาร์ทีเซียนโดยใช้พิกัดที่คาดการณ์ไว้) NS C หน่วย
ต้อง -R และ -J ที่จะตั้งค่า โดยไม่มีจุดคงที่เราคำนวณสะสม
ระยะทางตามเส้นทาง ผนวก - เพื่อให้ได้ระยะห่างที่เพิ่มขึ้นระหว่างต่อเนื่องกัน
คะแนน ผนวก + เพื่อระบุจุดที่ 2 ผ่านคอลัมน์พิเศษสองคอลัมน์ในไฟล์อินพุต

-I ทำการแปลงผกผัน เช่น รับ (ลองจิจูด ละติจูด) จากข้อมูล (x,y)

-Lไลน์.xy-หน่วย-
กำหนดระยะทางที่สั้นที่สุดจากจุดข้อมูลอินพุตไปยังเส้นที่ระบุใน
ไฟล์หลายส่วน ASCII ไลน์.xy. ระยะทางและพิกัดของ
จุดที่ใกล้ที่สุดจะถูกผนวกเข้ากับผลลัพธ์เป็นสามคอลัมน์ใหม่ ต่อท้าย
หน่วยระยะทาง (ดูหน่วย) รวมทั้ง c (ระยะทางคาร์ทีเซียนโดยใช้พิกัดอินพุต)
or C (ระยะทางคาร์ทีเซียนโดยใช้พิกัดที่คาดการณ์ไว้) NS C หน่วยที่ต้องการ -R และ
-J ที่จะตั้งค่า สุดท้ายผนวก + เพื่อรายงาน id ส่วนบรรทัดและเศษส่วน
เลขจุดแทน lon/lat ของจุดที่ใกล้ที่สุด

-N[a|c|g|m]
แปลงจากละติจูด geodetic (โดยใช้ทรงรีปัจจุบัน ดู PROJ_ELLIPSOID)
หนึ่งในสี่ละติจูดเสริมที่แตกต่างกัน (ลองจิจูดไม่ได้รับผลกระทบ) เลือก
ราคาเริ่มต้นที่ aยูทาลิก, cไม่เป็นทางการ, gศูนย์กลางและ mละติจูด eridional [geocentric] ใช้ -I
เพื่อแปลงจากละติจูดเสริมเป็นละติจูด geodetic

-Q[d|e แสดงรายการพารามิเตอร์การฉายภาพทั้งหมด หากต้องการแสดงเฉพาะ Datum ให้ใช้ -Qd. เฉพาะรายการ
ทรงรี ใช้ -คิว.

-S ระงับคะแนนที่อยู่นอกเขต

-ไทย]ราคาเริ่มต้นที่[/ไปยัง]
การแปลงพิกัดระหว่าง Datums ราคาเริ่มต้นที่ และ ไปยัง โดยใช้มาตรฐาน Molodensky
การเปลี่ยนแปลง ใช้ - ธ ถ้าคอลัมน์อินพุตที่ 3 มีความสูงเหนือวงรี [Default
ถือว่าความสูง = 0, คือบนทรงรี] ระบุ Datum โดยใช้ Datum ID (ดู
-Qd) หรือให้ ทรงรี:dx,dy,dz ที่ไหน ทรงรี อาจเป็น ID ทรงรี (ดู -คิว) หรือ
ให้เป็น a[,*inv_f*] โดยที่ a เป็นกึ่งแกนเอกและ inv_f เป็นตัวผกผัน
แบน (0 ถ้าละเว้น) ถ้า datums คือ - หรือไม่ให้ถือว่า WGS-84 -T อาจจะ
ใช้ร่วมกับ -R -J เพื่อเปลี่ยนจุดอ้างก่อนการฉายพิกัด
(เพิ่ม -I เพื่อใช้การแปลง Datum หลังจากการฉายภาพผกผัน) ทำให้เเน่นอน
การตั้งค่า PROJ_ELLIPSOID นั้นถูกต้องสำหรับกรณีของคุณ

-วี[ระดับ] (มากกว่า ... )
เลือกระดับการใช้คำฟุ่มเฟือย [c]

-ว[w|เอช]
พิมพ์ความกว้างและความสูงของแผนที่บนเอาต์พุตมาตรฐาน ไม่มีการอ่านไฟล์อินพุต เท่านั้น
เอาท์พุทความกว้างหรือความสูงผนวก w or hตามลำดับ หน่วยของ
ขนาดอาจเปลี่ยนแปลงได้ทาง -D.

-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 คือคอลัมน์แรก)

-oปลอกคอ- (มากกว่า ... )
เลือกคอลัมน์ผลลัพธ์ (0 คือคอลัมน์แรก)

-p[x|y|z]อาซิม/ยกระดับ[/ระดับ][+วLon0/ละติจูด0[/z0]][+วีx0/y0] (มากกว่า ... )
เลือกมุมมองเปอร์สเปคทีฟ

-NS[ปลอกคอ][ก|ร] (มากกว่า ... )
ตั้งค่าการจัดการระเบียน NaN

-:[ผม|o] (มากกว่า ... )
สลับคอลัมน์ที่ 1 และ 2 บนอินพุตและ/หรือเอาต์พุต

-^ 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

ตัวอย่าง

ในการแปลงพิกัด UTM เป็นเมตรเป็นตำแหน่งทางภูมิศาสตร์ ให้ไฟล์ utm.txt และ
รู้จักโซน UTM (และโซนหรือซีกโลก) ลอง

โครงการแผนที่ gmt utm.txt -Ju+11/1:1 -C -I -F

ในการแปลงไฟล์ที่มี (ลองจิจูด,ละติจูด) เป็น (x,y) ตำแหน่งในหน่วยซมบน Mercator
ตารางสำหรับมาตราส่วนที่กำหนด 0.5 ซม. ต่อองศา วิ่ง

gmt mapproject lonlatfile -R20/50/12/25 -Jm0.5c > xyfile

ในการแปลงไฟล์ 2 คอลัมน์, ไบนารี, ความแม่นยำสองเท่าหลายไฟล์ด้วย (ละติจูด, ลองจิจูด)
เป็น (x,y) ตำแหน่งในหน่วยนิ้วบนตาราง Transverse Mercator (ลองจิจูดกลาง 75W) สำหรับ
มาตราส่วน = 1:500000 และระงับจุดที่จะตกนอกพื้นที่แผนที่ ให้วิ่ง

แทร็ก gmt mapproject* -R-80/-70/20/40 -Jt-75/1:500000 -: -S -Di -bo -bi2 > tmfile.b

ในการแปลงพิกัดทางภูมิศาสตร์ (lon, lat, height) ในไฟล์ old.dat จาก NAD27
CONUS datum (Datum ID 131 ซึ่งใช้ทรงรี Clarke-1866) ถึง WGS 84 รัน

gmt mapproject old.dat -Th131 > new.dat

ในการคำนวณระยะทางที่ใกล้ที่สุด (เป็นกม.) ระหว่างแต่ละจุดในไฟล์อินพุต quakes.dat
และส่วนของเส้นตรงที่ระบุในไฟล์ ASCII แบบหลายส่วน coastline.xy, run

gmt mapproject quakes.dat -Lcoastline.xy/k > quake_dist.dat

ข้อ จำกัด

พื้นที่อินพุตสี่เหลี่ยมตั้งค่าด้วย -R โดยทั่วๆ ไป จะถูกแมปเป็นรูปไม่เหลี่ยม
ตาราง เว้นแต่ -C ถูกตั้งค่า จุดซ้ายสุดบนตารางนี้มี xvalue = 0.0 และ
จุดต่ำสุดจะมีค่า yvalue = 0.0 ดังนั้น ก่อนที่คุณจะแปลงแผนที่เป็นดิจิทัล ให้วิ่งสุดขั้ว
พิกัดแผนที่ผ่าน โครงการแผนที่ โดยใช้มาตราส่วนที่เหมาะสมและดูว่าค่าใด (x,y)
พวกมันถูกแมปไว้ ใช้ค่าเหล่านี้เมื่อตั้งค่าสำหรับการแปลงเป็นดิจิทัลเพื่อให้มี
การแปลงผกผันทำงานอย่างถูกต้องหรืออีกวิธีหนึ่งคือใช้ awk เพื่อปรับขนาดและเปลี่ยน
(x,y) ค่าก่อนการแปลง

สำหรับการฉายภาพบางส่วน อาจใช้สารละลายทรงกลมแม้ว่าผู้ใช้จะเลือก an
ทรงรี สิ่งนี้เกิดขึ้นเมื่อผู้ใช้ -R การตั้งค่าหมายถึงภูมิภาคที่เกินโดเมน
ซึ่งการขยายอนุกรมวงรีนั้นใช้ได้ โดยมีเงื่อนไขดังต่อไปนี้ (1)
แลมเบิร์ต Conformal Conic (-เจแอล) และ Albers Equal-Area (-เจบี) จะใช้สารละลายทรงกลม
เมื่อมาตราส่วนแผนที่เกิน 1.0E7 (2) Mercator ขวาง (-เจที) และ UTM (-จู) จะ
ใช้สารละลายทรงกลมเมื่อกำหนดเขตแดนตะวันตกหรือตะวันออกใน -R เป็นมากกว่า
10 องศาจากเส้นเมอริเดียนกลาง และ (3) เหมือนกันสำหรับ Cassini (-เจ.ซี) แต่มีขีดจำกัดของ
เพียง 4 องศา

วงรี AND สเฟียรอยด์

GMT จะใช้สูตรวงรีหากมีการใช้งานและผู้ใช้เลือก an
ทรงรีเป็นรูปร่างอ้างอิง (ดู PROJ_ELLIPSOID) ผู้ใช้จำเป็นต้องรู้เท่าทัน
ข้อผิดพลาดที่อาจเกิดขึ้น: (1) สำหรับการคาดการณ์บางอย่าง เช่น Transverse Mercator, Albers และ
รูปกรวยของแลมเบิร์ต เราใช้นิพจน์วงรีเมื่อพื้นที่ที่แมปคือ
เล็กและเปลี่ยนเป็นนิพจน์ทรงกลม (และแทนที่ตัวช่วยที่เหมาะสม
ละติจูด) สำหรับแผนที่ขนาดใหญ่ขึ้น สูตรทรงรีใช้ดังนี้ (ก) ขวาง
Mercator: เมื่อจุดทั้งหมดอยู่ภายใน 10 องศาของเส้นเมอริเดียนกลาง (b) การคาดคะเนรูปกรวย
เมื่อพิสัยตามยาวน้อยกว่า 90 องศา (ค) การฉายแคสสินีเมื่อจุดทั้งหมด
อยู่ภายใน 4 องศาของเส้นเมริเดียนกลาง (2) เมื่อคุณพยายามที่จะจับคู่ประวัติศาสตร์บางอย่าง
ข้อมูล (เช่น พิกัดที่ได้รับจากการฉายภาพและการอ้างอิงบางอย่าง
ทรงรี) คุณอาจพบว่า GMT ให้ผลลัพธ์ที่แตกต่างกันเล็กน้อย เป็นไปได้
แหล่งที่มาของความไม่ตรงกันนี้คือการคำนวณแบบเก่ามักใช้ตัวเลขที่มีนัยสำคัญน้อยกว่า สำหรับ
ตัวอย่างของ Snyder มักใช้ทรงรีคลาร์ก 1866 (กำหนดโดยเขาว่าเป็นมี
แฟบ f = 1/294.98) จาก f เราได้ความเยื้องศูนย์กำลังสองเป็น 0.00676862818
(นี่คือสิ่งที่ GMT ใช้) ในขณะที่ไนเดอร์ปัดเศษและใช้ 0.00676866 ความแตกต่างนี้สามารถ
ให้ความคลาดเคลื่อนหลายสิบซม. หากคุณต้องการสร้างพิกัดที่คาดการณ์ไว้
ด้วยความเยื้องศูนย์กลางที่แตกต่างกันเล็กน้อยนี้ คุณควรระบุทรงรีของคุณเองด้วย
พารามิเตอร์เดียวกับคลาร์ก 1866 แต่ด้วย f = 1/294.97861076 นอกจากนี้โปรดทราบว่าผู้สูงอายุ
ข้อมูลอาจถูกอ้างอิงไปยัง Datum ที่แตกต่างกัน และเว้นแต่คุณจะรู้ว่ามีการใช้ Datum ใดและ
แปลงข้อมูลทั้งหมดเป็นข้อมูลทั่วไป คุณอาจพบว่าไม่ตรงกันตั้งแต่หลักสิบถึงหลายร้อย
เมตร (3) สุดท้าย โปรดทราบว่า PROJ_SCALE_FACTOR มีค่าเริ่มต้นบางอย่างสำหรับบางค่า
ประมาณการดังนั้นคุณอาจต้องแทนที่การตั้งค่าเพื่อให้ตรงกับผลลัพธ์ที่ผลิต
กับการตั้งค่าอื่นๆ

ใช้ mapprojectgmt ออนไลน์โดยใช้บริการ onworks.net