mapprojectgmt - Trực tuyến trên đám mây

CHƯƠNG TRÌNH:

TÊN

mapproject - Thực hiện các phép biến đổi bản đồ chuyển tiếp và đảo ngược, chuyển đổi dữ liệu và đo đạc

SYNOPSIS

dự án bản đồ [ bảng ] thông số khu [ b|B|f|F|o|O[lớn0/vĩ độ0]] [[dx/dy]] [ c|i|p ] [
[datums]] [[đơn vị]] [[x0/y0/] [[+|-]đơn vị][+|-]] [] [ dòng.xy[/ [+|-]đơn vị][+]] [
[a|c|g|m]] [[d|e ] [] [[h]từ[/đến]] [[cấp]] [[w|h]] [ -b] [
-d] [ -f] [ -g] [ -h] [ -i] [ -o] [
-p] [ -s] [ -:[i|o]]

Lưu ý: Không được phép có khoảng trống giữa cờ tùy chọn và các đối số liên quan.

MÔ TẢ

dự án bản đồ đọc vị trí (kinh độ, vĩ độ) từ bảng [hoặc đầu vào tiêu chuẩn] và
tính toán (x, y) tọa độ bằng cách sử dụng tỷ lệ và phép chiếu bản đồ đã chỉ định. Tùy ý, nó
có thể đọc các vị trí (x, y) và tính toán các giá trị (kinh độ, vĩ độ) làm ngược lại
sự biến đổi. Điều này có thể được sử dụng để biến đổi các điểm tuyến tính (x, y) thu được bằng cách số hóa
một bản đồ của phép chiếu đã biết đến tọa độ địa lý. Cũng có thể tính toán khoảng cách cùng
theo dõi, đến một điểm cố định hoặc cách tiếp cận gần nhất với một đường. Cuối cùng, có thể được sử dụng để thực hiện
các chuyển đổi dữ liệu khác nhau. Các trường dữ liệu bổ sung được phép sau 2 cột đầu tiên
mà phải có (kinh độ, vĩ độ) hoặc (x, y). Xem tùy chọn -: về cách đọc
(vĩ độ, kinh độ) tệp.

YÊU CẦU TRANH LUẬN

-Jthông số (hơn ...)
Chọn phép chiếu bản đồ.

-NS[đơn vị]xmin/xmax/ymin/ymax[NS] (hơn ...)
Chỉ định khu vực quan tâm. Trường hợp đặc biệt cho phép chiếu UTM: Nếu -C được sử dụng
và -R không được đưa ra thì vùng này được đặt trùng với vùng UTM đã cho để
để bảo toàn dung dịch ellipsoidal đầy đủ (Xem HẠN CHẾ để biết thêm
thông tin).

CHỌN TRANH LUẬN

bàn Một hoặc nhiều ASCII (hoặc nhị phân, xem -bi[ncols][kiểu]) (các) tệp bảng dữ liệu chứa một
số cột dữ liệu. Nếu không có bảng nào được đưa ra thì chúng tôi đọc từ đầu vào chuẩn.

-Ab | B | f | F | o | O [lớn0/vĩ độ0]
-Nếu tính toán góc phương vị (về phía trước) từ điểm cố định lon / lat đến từng điểm dữ liệu.
Sử dụng -ab để lấy lại góc phương vị từ điểm dữ liệu đến điểm cố định. Sử dụng -Không để có được
định hướng (-90/90) thay vì phương vị (0/360). Chữ hoa F, B or O sẽ
chuyển đổi từ vĩ độ trắc địa sang vĩ độ địa tâm và ước tính phương vị của trắc địa
(giả sử ellipsoid hiện tại không phải là hình cầu). Nếu không có điểm cố định nào được đưa ra thì
chúng tôi tính góc phương vị (hoặc phương vị quay lại) từ điểm trước đó.

-NS[dx/dy]
Đặt tâm của các tọa độ được chiếu thành tâm chiếu của bản đồ [Mặc định là
góc dưới bên trái]. Theo tùy chọn, thêm các hiệu số trong các đơn vị dự kiến ​​sẽ được thêm vào (hoặc
trừ khi -I được đặt) thành (từ) tọa độ dự kiến, chẳng hạn như false
các hướng mở rộng và hướng đông bắc cho các vùng chiếu cụ thể [0/0]. Đơn vị được sử dụng cho
hiệu số là đơn vị khoảng cách âm mưu có hiệu lực (xem PROJ_LENGTH_UNIT) trừ khi -F is
được sử dụng, trong trường hợp đó hiệu số tính bằng mét.

-Dc | i | p
Tạm thời ghi đè PROJ_LENGTH_UNIT và sử dụng c (cm), i (inch), hoặc p (điểm)
thay thế. Không thể được sử dụng với -F.

-E [datums]
Chuyển đổi từ trắc địa (vĩ độ, vĩ độ, độ cao) sang Trái đất có tâm điểm cố định (ECEF)
(x, y, z) tọa độ (thêm -I cho phép chuyển đổi nghịch đảo). Nối ID dữ liệu (xem -QĐ)
Hoặc cho dấu chấm lửng:dx,dy,dz Ở đâu dấu chấm lửng có thể là một ID ellipsoid (xem -Qe) hoặc
đưa ra như a[, * inv_f *], ở đâu a là trục bán chính và inv_f là nghịch đảo
làm phẳng (0 nếu bỏ qua). Nếu như datums là - hoặc không được cung cấp, chúng tôi giả định là WGS-84.

-NS[đơn vị]
Buộc chia tỷ lệ 1: 1, tức là, đầu ra (hoặc đầu vào, xem -I) dữ liệu trong thực tế dự kiến
mét. Để chỉ định các đơn vị khác, hãy nối đơn vị mong muốn (xem ĐƠN VỊ). Không có -F,
đầu ra (hoặc đầu vào, xem -I) nằm trong các đơn vị được chỉ định bởi PROJ_LENGTH_UNIT (nhưng
xem -D).

-NS[x0/y0/] [[+ | -]đơn vị] [+ | -]
Tính toán khoảng cách dọc theo đường or đến điểm tùy chọn được đặt với -Gx0 / y0. Nối
đơn vị khoảng cách (xem UNITS), bao gồm c (Khoảng cách Descartes sử dụng đầu vào
tọa độ) hoặc C (Khoảng cách Descartes sử dụng tọa độ dự phóng). Các C đơn vị
đòi hỏi -R và -J được thiết lập. Không có điểm cố định nào được đưa ra, chúng tôi tính tích lũy
khoảng cách dọc theo đường. Nối - để có được khoảng cách gia tăng giữa các
điểm. Nối + để chỉ định điểm thứ 2 thông qua hai cột phụ trong tệp đầu vào.

-I Thực hiện phép biến đổi Nghịch đảo, tức là lấy (kinh độ, vĩ độ) từ dữ liệu (x, y).

-Ldòng.xy[/ [+ | -]đơn vị] [+]
Xác định khoảng cách ngắn nhất từ ​​các điểm dữ liệu đầu vào đến (các) dòng được cho trong
tệp đa phân đoạn ASCII dòng.xy. Khoảng cách và tọa độ của
điểm gần nhất sẽ được thêm vào đầu ra dưới dạng ba cột mới. Nối
đơn vị khoảng cách (xem UNITS), bao gồm c (Khoảng cách Descartes sử dụng tọa độ đầu vào)
or C (Khoảng cách Descartes sử dụng tọa độ dự phóng). Các C đơn vị yêu cầu -R và
-J được thiết lập. Cuối cùng, nối thêm + để báo cáo id phân đoạn dòng và phân số
số điểm thay vì lon / vĩ độ của điểm gần nhất.

-N [a | c | g | m]
Chuyển đổi từ vĩ độ trắc địa (sử dụng ellipsoid hiện tại; xem PROJ_ELLIPSOID)
đến một trong bốn vĩ độ phụ khác nhau (các kinh độ không bị ảnh hưởng). Lựa chọn
từ akhông tưởng, cchính quy, geocentric, và mvĩ độ kinh tuyến [địa tâm]. Sử dụng -I
để chuyển đổi từ vĩ độ phụ sang vĩ độ trắc địa.

-Q [d | e Liệt kê tất cả các thông số hình chiếu. Để chỉ liệt kê các mức dữ liệu, hãy sử dụng -QĐ. Để chỉ danh sách
ellipsoids, sử dụng -Qe.

-S Áp chế các điểm nằm ngoài vùng.

-Thứ tự]từ[/đến]
Phối hợp chuyển đổi giữa các mức dữ liệu từ và đến sử dụng Molodensky tiêu chuẩn
sự biến đổi. Sử dụng -Thứ tự nếu cột đầu vào thứ 3 có chiều cao trên ellipsoid [Mặc định
giả sử chiều cao = 0, tức là, trên ellipsoid]. Chỉ định các mức dữ liệu bằng cách sử dụng ID dữ liệu (xem
-QĐ) Hoặc cho dấu chấm lửng:dx,dy,dz Ở đâu dấu chấm lửng có thể là một ID ellipsoid (xem -Qe) hoặc
đưa ra như a[, * inv_f *], ở đâu a là trục bán chính và inv_f là nghịch đảo
làm phẳng (0 nếu bỏ qua). Nếu như datums là - hoặc không được cung cấp, chúng tôi giả định là WGS-84. -T có lẽ
được sử dụng kết hợp với -R -J để thay đổi dữ liệu trước khi chiếu tọa độ
(thêm vào -I để áp dụng chuyển đổi dữ liệu sau phép chiếu nghịch đảo). Đảm bảo rằng
cài đặt PROJ_ELLIPSOID chính xác cho trường hợp của bạn.

-V [cấp] (hơn ...)
Chọn mức độ chi tiết [c].

-W [w | h]
In chiều rộng và chiều cao của bản đồ trên đầu ra tiêu chuẩn. Không có tệp đầu vào nào được đọc. Chỉ
xuất ra chiều rộng hoặc chiều cao, nối thêm w or h, tương ứng. Các đơn vị của
kích thước có thể được thay đổi thông qua -D.

-bi [ncols][NS] (hơn ...)
Chọn đầu vào nhị phân gốc. [Mặc định là 2 cột đầu vào].

-bo [ncols][kiểu] (hơn ...)
Chọn đầu ra nhị phân gốc. [Mặc định giống như đầu vào].

-d [i | o]không có dữ liệu (hơn ...)
Thay thế các cột đầu vào bằng không có dữ liệu với NaN và làm ngược lại trên đầu ra.

-f [i | o]thông tin về colinfo (hơn ...)
Chỉ định kiểu dữ liệu của cột đầu vào và / hoặc cột đầu ra.

-g [a] x | y | d | X | Y | D | [col] z [+ | -]khoảng cách[hoặc là] (hơn ...)
Xác định khoảng trống dữ liệu và ngắt dòng.

-h [i | o] [n] [+ c] [+ d] [+ rnhận xét] [+ rtiêu đề] (hơn ...)
Bỏ qua hoặc tạo (các) bản ghi tiêu đề.

-ivòng cổ[l] [stỉ lệ] [obù đắp] [,hữu ích. Cảm ơn !] (hơn ...)
Chọn các cột đầu vào (0 là cột đầu tiên).

-ovòng cổ[, ...] (hơn ...)
Chọn các cột đầu ra (0 là cột đầu tiên).

-p [x | y | z]phương vị/độ cao[/cấp độ z] [+ wlớn0/vĩ độ0[/z0]] [+ vx0/y0] (hơn ...)
Chọn hình chiếu phối cảnh.

-NS[vòng cổ] [a | r] (hơn ...)
Đặt xử lý các bản ghi NaN.

-: [i | o] (hơn ...)
Hoán đổi cột thứ nhất và thứ hai trên đầu vào và / hoặc đầu ra.

-^ or chỉ -
In một thông báo ngắn về cú pháp của lệnh, sau đó thoát ra (LƯU Ý: trên Windows
chỉ sử dụng -).

-+ or chỉ +
In thông báo sử dụng rộng rãi (trợ giúp), bao gồm giải thích về bất kỳ
tùy chọn dành riêng cho mô-đun (nhưng không phải tùy chọn chung theo GMT), sau đó thoát.

-? or Không đối số
In thông báo sử dụng hoàn chỉnh (trợ giúp), bao gồm giải thích về các tùy chọn, sau đó
lối thoát hiểm.

--phiên bản
In phiên bản GMT và thoát.

--show-datadir
In đường dẫn đầy đủ đến thư mục chia sẻ GMT và thoát.

CÁC ĐƠN VỊ

Đối với đơn vị khoảng cách bản đồ, hãy nối thêm đơn vị d đối với độ cung, m cho phút cung, và s cho vòng cung
thứ hai, hoặc e cho mét [Mặc định], f cho chân, k cho km, M cho số dặm quy chế, n cho hải lý
dặm, và u cho chân khảo sát Hoa Kỳ. Theo mặc định, chúng tôi tính toán các khoảng cách như vậy bằng cách sử dụng một hình cầu
xấp xỉ với các vòng tròn lớn. Thêm trước - đến một khoảng cách (hoặc đơn vị là không có khoảng cách là
cho trước) để thực hiện các phép tính "Trái đất phẳng" (nhanh hơn nhưng kém chính xác hơn) hoặc thêm vào trước + đến
thực hiện các phép tính trắc địa chính xác (chậm hơn nhưng chính xác hơn).

ASCII FORMAT ĐỘ CHÍNH XÁC

Định dạng đầu ra ASCII của dữ liệu số được kiểm soát bởi các tham số trong gmt.conf
tập tin. Kinh độ và vĩ độ được định dạng theo FORMAT_GEO_OUT, ngược lại
giá trị được định dạng theo FORMAT_FLOAT_OUT. Lưu ý rằng định dạng có hiệu lực có thể
dẫn đến mất độ chính xác trong đầu ra, có thể dẫn đến các vấn đề khác nhau ở hạ nguồn. Nếu như
bạn thấy đầu ra không được viết với đủ độ chính xác, hãy cân nhắc chuyển sang hệ nhị phân
đầu ra (-bo nếu có) hoặc chỉ định nhiều số thập phân hơn bằng cách sử dụng cài đặt FORMAT_FLOAT_OUT.

VÍ DỤ

Để chuyển đổi tọa độ UTM tính bằng mét thành vị trí địa lý, hãy cung cấp một tệp utm.txt và
biết vùng UTM (và vùng hoặc bán cầu), hãy thử

gmt mapproject utm.txt -Ju + 11/1: 1 -C -I -F

Để chuyển đổi tệp có (kinh độ, vĩ độ) thành (x, y) vị trí tính bằng cm trên Mercator
lưới cho một tỷ lệ nhất định là 0.5 cm mỗi độ, chạy

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

Để chuyển đổi một số tệp 2 cột, nhị phân, chính xác kép với (vĩ độ, kinh độ)
vào vị trí (x, y) tính bằng inch trên lưới Mercator ngang (kinh độ trung tâm 75W) cho
scale = 1: 500000 và chặn những điểm nằm ngoài khu vực bản đồ, chạy

gmt mapproject track. * -R-80 / -70 / 20/40 -Jt-75/1: 500000 -: -S -Di -bo -bi2> tmfile.b

Để chuyển đổi tọa độ trắc địa (vĩ độ, vĩ độ, độ cao) trong tệp old.dat từ NAD27
Dữ liệu CONUS (ID dữ liệu 131 sử dụng ellipsoid Clarke-1866) đến WGS 84, chạy

gmt mapproject old.dat -Th131> new.dat

Để tính khoảng cách gần nhất (tính bằng km) giữa mỗi điểm trong tệp đầu vào quakes.dat
và các phân đoạn dòng được đưa ra trong tệp Coastline.xy đa phân đoạn ASCII, chạy

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

NHỮNG HẠN CHẾ

Vùng đầu vào hình chữ nhật được đặt với -R nói chung sẽ được ánh xạ vào một không phải hình chữ nhật
lưới điện. Trừ phi -C được đặt, điểm ngoài cùng bên trái trên lưới này có xvalue = 0.0 và
điểm thấp nhất sẽ có yvalue = 0.0. Do đó, trước khi bạn số hóa bản đồ, hãy chạy cực
bản đồ tọa độ thông qua dự án bản đồ sử dụng thang đo thích hợp và xem giá trị (x, y) là gì
chúng được ánh xạ vào. Sử dụng các giá trị này khi thiết lập số hóa để có
biến đổi nghịch đảo hoạt động chính xác, hoặc cách khác, sử dụng ôi để mở rộng quy mô và thay đổi
(x, y) các giá trị trước khi biến đổi.

Đối với một số phép chiếu, giải pháp hình cầu có thể được sử dụng mặc dù người dùng đã chọn
ellipsoid. Điều này xảy ra khi người dùng -R cài đặt ngụ ý một vùng vượt quá tên miền
trong đó các phần mở rộng của chuỗi hình elip là hợp lệ. Đây là các điều kiện: (1)
Lambert Conic phù hợp (-JL) và Albers Equal-Area (-JB) sẽ sử dụng giải pháp hình cầu
khi tỷ lệ bản đồ vượt quá 1.0E7. (2) Mercator ngang (-JT) và UTM (-JU) sẽ sẽ
sử dụng giải pháp hình cầu khi ranh giới phía tây hoặc phía đông được đưa vào -R nhiều hơn
10 độ từ kinh tuyến trung tâm, và (3) như vậy đối với Cassini (-JC) nhưng với giới hạn là
chỉ 4 độ.

ELLIPSOID VÀ HÌNH CẦU

GMT sẽ sử dụng công thức hình elip nếu chúng được triển khai và người dùng đã chọn
ellipsoid làm hình dạng tham chiếu (xem PROJ_ELLIPSOID). Người dùng cần lưu ý một số
cạm bẫy tiềm ẩn: (1) Đối với một số phép chiếu, chẳng hạn như Transverse Mercator, Albers, và
Conic hình cầu của Lambert, chúng tôi sử dụng các biểu thức hình elip khi các khu vực được ánh xạ là
nhỏ và chuyển sang các biểu thức hình cầu (và thay thế phụ trợ thích hợp
vĩ độ) cho các bản đồ lớn hơn. Các công thức hình elip được sử dụng như sau: (a) Ngang
Mercator: Khi tất cả các điểm nằm trong phạm vi 10 độ của kinh tuyến trung tâm, (b) Phép chiếu hình nón
khi phạm vi dọc nhỏ hơn 90 độ, (c) Phép chiếu Cassini khi tất cả các điểm
nằm trong phạm vi 4 độ của kinh tuyến trung tâm. (2) Khi bạn đang cố gắng so khớp một số
dữ liệu (ví dụ: tọa độ thu được với một phép chiếu nhất định và một tham chiếu nhất định
ellipsoid) bạn có thể thấy rằng GMT cho kết quả hơi khác một chút. Một khả năng
nguồn gốc của sự không khớp này là do các phép tính cũ hơn thường sử dụng các chữ số có nghĩa ít hơn. Vì
ví dụ, các ví dụ của Snyder thường sử dụng ellipsoid Clarke 1866 (được ông định nghĩa là có
a phẳng f = 1 / 294.98). Từ f, chúng ta nhận được bình phương độ lệch tâm là 0.00676862818
(đây là những gì GMT sử dụng), trong khi Snyder làm tròn và sử dụng 0.00676866. Sự khác biệt này có thể
cho những sai lệch vài chục cm. Nếu bạn cần tạo lại các tọa độ được chiếu
với độ lệch tâm hơi khác này, bạn nên chỉ định ellipsoid của riêng mình với
thông số tương tự như Clarke 1866, nhưng với f = 1 / 294.97861076. Ngoài ra, hãy lưu ý rằng cũ hơn
dữ liệu có thể được tham chiếu đến các vùng dữ liệu khác nhau và trừ khi bạn biết mức dữ liệu nào đã được sử dụng và
chuyển đổi tất cả dữ liệu sang một mức dữ liệu chung mà bạn có thể gặp phải sự không khớp từ hàng chục đến hàng trăm
mét. (3) Cuối cùng, hãy lưu ý rằng PROJ_SCALE_FACTOR có các giá trị mặc định nhất định đối với một số
các phép chiếu, do đó bạn có thể phải ghi đè cài đặt để khớp với kết quả được tạo ra
với các cài đặt khác.

Sử dụng mapprojectgmt trực tuyến bằng các dịch vụ onworks.net