mapprojectgmt - Bulutta Çevrimiçi

Program:

ADI

mapproject - İleri ve ters harita dönüşümleri, veri dönüşümleri ve jeodezi yapın

SİNOPSİS

harita projesi [ tablolar ] parametreler bölge [ b|B|f|F|o|O[boy0/lat0] ] [dx/dy] ] [ c|i|p ] [
[datumlar] ] [birim] ] [x0/y0/][[+|-]birim][+|-] ] [ ] [ satır.xy[/[+|-]birim][+] ] [
[a|c|g|m] ] [d|e ] [ ] [ [h]itibaren[/için] ] [seviye] ] [w|h] ] [ -b] [
-d] [ -f] [ -g] [ -h] [ -i] [ -o] [
-p] [ -s] [ -:[i|o]]

Not: Seçenek bayrağı ve ilişkili bağımsız değişkenler arasında boşluk bırakılamaz.

TANIM

harita projesi (boylam, enlem) konumlarını okur tablolar [veya standart girdi] ve
(x,y) koordinatlarını belirtilen harita projeksiyonu ve ölçekleri kullanarak hesaplar. İsteğe bağlı olarak,
(x,y) konumlarını okuyabilir ve tersini yaparak (boylam, enlem) değerleri hesaplayabilir
dönüşüm. Bu, sayısallaştırılarak elde edilen doğrusal (x,y) noktaları dönüştürmek için kullanılabilir.
coğrafi koordinatlara bilinen projeksiyon haritası. boyunca mesafeleri de hesaplayabilir
iz, sabit bir noktaya veya bir hatta en yakın yaklaşım. Son olarak, gerçekleştirmek için kullanılabilir
çeşitli veri dönüşümleri. İlk 2 sütundan sonra ek veri alanlarına izin verilir
(boylam, enlem) veya (x,y) olmalıdır. Seçeneğe bakın -: nasıl okunacağına dair
(enlem, boylam) dosyaları.

GEREKLİ ARGÜMANLAR

-Jparametreler (Daha ...)
Harita projeksiyonunu seçin.

-R[birim]x dakika/xmaks/imin/ymaks[D] (Daha ...)
İlgi alanı belirtin. UTM projeksiyonu için özel durum: Eğer -C kullanılır
ve -R verilmezse bölge, verilen UTM bölgesi ile çakışacak şekilde ayarlanır, bu nedenle
tam elipsoidal çözümü korumak için (daha fazlası için KISITLAMALAR'a bakın)
bilgi).

İSTEĞE BAĞLI ARGÜMANLAR

tablo Bir veya daha fazla ASCII (veya ikili, bkz. -bi[ncol'ler][tip]) tutan veri tablosu dosya(lar)ı
veri sütunlarının sayısı. Tablo verilmezse standart girdiden okuruz.

-Ab|B|f|F|o|O[boy0/lat0]
-Af (ileri) azimutu sabit noktadan hesaplar boylam/enlem her veri noktasına
kullanım -Uzak veri noktalarından sabit noktaya geri azimut almak için. Kullanmak -Ao almak
azimutlar (90/90) yerine yönler (-0/360). büyük harf F, B or O irade
jeodezikten jeosentrik enlemlere dönüştürün ve jeodeziklerin azimutunu tahmin edin
(mevcut elipsoidin bir küre olmadığı varsayılarak). Sabit bir nokta verilmezse, o zaman
azimutu (veya arka azimutu) önceki noktadan hesaplıyoruz.

-C[dx/dy]
Tahmini koordinatların merkezini harita projeksiyon merkezinde olacak şekilde ayarlayın [Varsayılan
sol alt köşe]. İsteğe bağlı olarak, eklenecek öngörülen birimlere ofsetler ekleyin (veya
ne zaman çıkarılır -I false gibi öngörülen koordinatlara (dan) ayarlanır
belirli izdüşüm bölgeleri [0/0] için doğu yönü ve kuzey yönü. için kullanılan birim
ofsetler, geçerli olan çizim mesafesi birimidir (bkz. PROJ_LENGTH_UNIT) -F is
kullanılır, bu durumda ofsetler metre cinsindendir.

-Dc|i|p
PROJ_LENGTH_UNIT'i geçici olarak geçersiz kıl ve kullan c (santimetre), i (inç) veya p (puan)
Bunun yerine. ile kullanılamaz -F.

-E[datumlar]
Jeodeziden (boylam, enlem, yükseklik) Dünya Merkezli Yer Sabiti'ne (ECEF) dönüştürme
(x,y,z) koordinatları (ekleyin -I ters dönüşüm için). Veri kimliğini ekleyin (bkz. -Qd)
veya Ver elipsoit:dx,dy,dz nerede elipsoit elipsoid kimliği olabilir (bkz. -Qe) Veya
olarak verildi a[,*inv_f*], nerede a yarı ana eksendir ve inv_f tersi
düzleştirme (atlanmışsa 0). Eğer datumlar - veya verilmemişse WGS-84'ü varsayıyoruz.

-F[birim]
1:1 ölçeklendirmeyi zorla, yani çıktı (veya giriş, bkz. -I) veriler gerçek tahmindedir
metre. Diğer birimleri belirtmek için istediğiniz birimi ekleyin (bkz. BİRİMLER). Olmadan -F,
çıktı (veya girdi, bkz. -I) PROJ_LENGTH_UNIT tarafından belirtilen birimlerdedir (ancak
görmek -D).

-G[x0/y0/][[+|-]birim][+|-]
Pist boyunca mesafeleri hesaplayın or ile ayarlanan isteğe bağlı noktaya -Gx0/y0. Ekle
mesafe birimi (bkz. BİRİMLER), dahil c (Girdi kullanılarak kartezyen uzaklık
koordinatlar) veya C (Öngörülen koordinatları kullanarak kartezyen mesafesi). NS C birim
gerektirir -R ve -J ayarlanacak. Sabit bir nokta verilmeden kümülatı hesaplıyoruz
yol boyunca mesafeler. Ekle - ardışık arasındaki artan mesafeyi elde etmek için
puan. Ekle + giriş dosyasındaki iki ekstra sütun aracılığıyla 2. noktayı belirtmek için.

-I Ters dönüşümü yapın, yani (x,y) verisinden (boylam,enlem) alın.

-Lsatır.xy[/[+|-]birim][+]
Giriş veri noktalarından verilen hat(lar)a olan en kısa mesafeyi belirleyin.
ASCII çok parçalı dosya satır.xy. uzaklığı ve koordinatları
en yakın nokta çıktıya üç yeni sütun olarak eklenecektir. Ekle
mesafe birimi (bkz. BİRİMLER), dahil c (Giriş koordinatlarını kullanarak kartezyen mesafesi)
or C (Öngörülen koordinatları kullanarak kartezyen mesafesi). NS C birim gerektirir -R ve
-J ayarlanacak. Son olarak, ekle + satır segmenti kimliğini ve kesirliyi bildirmek için
en yakın noktanın boylam/lat yerine nokta numarası.

-N[a|c|g|m]
Jeodezik enlemlerden dönüştürme (geçerli elipsoidi kullanarak; bkz. PROJ_ELLIPSOID)
dört farklı yardımcı enlemden birine (boylamlar etkilenmez). Seçmek
itibaren autalik, cresmi, geosentrik ve meridional enlemler [jeosentrik]. Kullanmak -I
yardımcı enlemlerden jeodezik enlemlere dönüştürmek için.

-Q[d|e Tüm projeksiyon parametrelerini listeleyin. Yalnızca verileri listelemek için şunu kullanın: -Qd. Sadece listelemek için
elipsoidler, kullanım -Qe.

-S Bölgenin dışına düşen noktaları bastırın.

-NS]itibaren[/için]
Veriler arasındaki koordinat dönüşümleri itibaren ve için standart Molodensky kullanarak
dönüşüm. Kullanmak -T 3. giriş sütununun yüksekliği elipsoidin üzerindeyse [Varsayılan
yükseklik = 0, yani elipsoid üzerinde varsayılır]. Veri kimliğini kullanarak verileri belirtin (bkz.
-Qd) yada ver elipsoit:dx,dy,dz nerede elipsoit elipsoid kimliği olabilir (bkz. -Qe) Veya
olarak verildi a[,*inv_f*], nerede a yarı ana eksendir ve inv_f tersi
düzleştirme (atlanmışsa 0). Eğer datumlar - veya verilmemişse WGS-84'ü varsayıyoruz. -T olabilir
ile birlikte kullanılır -R -J koordinat projeksiyonundan önce veriyi değiştirmek için
(Ekle -I ters projeksiyondan sonra veri dönüşümünü uygulamak için). Emin olun
PROJ_ELLIPSOID ayarı sizin durumunuz için doğrudur.

-V[seviye] (Daha ...)
Ayrıntı düzeyini seçin [c].

-W[w|h]
Standart çıktıda harita genişliğini ve yüksekliğini yazdırır. Hiçbir girdi dosyası okunmaz. sadece
genişliği veya yüksekliği çıktılayın, ekleyin w or h, sırasıyla. birimleri
boyutlar aracılığıyla değiştirilebilir -D.

-bi[ncol'ler][T] (Daha ...)
Yerel ikili girişi seçin. [Varsayılan 2 giriş sütunudur].

-Bö[ncol'ler][tip] (Daha ...)
Yerel ikili çıktıyı seçin. [Varsayılan girişle aynıdır].

-d[i|o]veri yok (Daha ...)
Eşit olan giriş sütunlarını değiştirin veri yok NaN ile ve çıkışta tersini yapın.

-f[i|o]kolinfo (Daha ...)
Giriş ve/veya çıkış sütunlarının veri türlerini belirtin.

-g[a]x|y|d|X|Y|D|[col]z[+|-]boşluk[A] (Daha ...)
Veri boşluklarını ve satır sonlarını belirleyin.

-h[i|o][n][+c][+d][+rdüşünce][+rbaşlık] (Daha ...)
Başlık kayıtlarını atlayın veya oluşturun.

-iyaka[l][lerölçek][Öofset][,...] (Daha ...)
Giriş sütunlarını seçin (0 ilk sütundur).

-oyaka[,...] (Daha ...)
Çıktı sütunlarını seçin (0 ilk sütundur).

-p[x|y|z]azizim/elev[/seviye][+wboy0/lat0[/z0]][+vx0/y0] (Daha ...)
Perspektif görünümünü seçin.

-S[yaka][a|r] (Daha ...)
NaN kayıtlarının işlenmesini ayarlayın.

-:[i|o] (Daha ...)
Giriş ve/veya çıkışta 1. ve 2. sütunu değiştirin.

-^ or sadece -
Komutun sözdizimi hakkında kısa bir mesaj yazdırın, ardından çıkar (NOT: Windows'ta
sadece kullan -).

-+ or sadece +
Açıklama da dahil olmak üzere kapsamlı bir kullanım (yardım) mesajı yazdırın.
modüle özgü seçenek (ancak GMT ortak seçenekleri değil), ardından çıkar.

-? or yok hayır argümanlar
Seçeneklerin açıklamasını da içeren eksiksiz bir kullanım (yardım) mesajı yazdırın, ardından
çıkışlar.

--versiyon
GMT sürümünü yazdırın ve çıkın.

--gösteri-veri dizini
GMT paylaşım dizininin tam yolunu yazdırın ve çıkın.

BİRİMLERİ

Harita mesafe birimi için, ekleyin birim d ark derecesi için, m ark dakikası için ve s ark için
ikinci veya e metre için [Varsayılan], f ayak için, k km için, M kanun mili için, n denizcilik için
mil ve u ABD anket ayağı için. Varsayılan olarak, bu tür mesafeleri bir küresel kullanarak hesaplarız.
büyük daireler ile yaklaşım. başa ekle - bir mesafeye (veya birim mesafe yok ise)
verilen) "Düz Dünya" hesaplamaları yapmak (daha hızlı ancak daha az doğru) veya başa eklemek için + için
kesin jeodezik hesaplamalar yapın (daha yavaş ama daha doğru).

ASCII FORMAT HASSAS

Sayısal verilerin ASCII çıktı biçimleri, bilgisayarınızdaki parametreler tarafından kontrol edilir. gmt.conf
dosya. Boylam ve enlem FORMAT_GEO_OUT'a göre biçimlendirilirken, diğer
değerler FORMAT_FLOAT_OUT'a göre biçimlendirilir. Geçerli olan formatın
çıktıda hassasiyet kaybına yol açar, bu da akış yönünde çeşitli sorunlara yol açabilir. Eğer
çıktının yeterli hassasiyetle yazılmadığını fark ederseniz, ikiliye geçmeyi düşünün
çıktı (-Bö varsa) veya FORMAT_FLOAT_OUT ayarını kullanarak daha fazla ondalık sayı belirtin.

ÖRNEKLER

UTM koordinatlarını metre cinsinden coğrafi konumlara dönüştürmek için utm.txt dosyası ve
UTM bölgesini (ve bölge veya yarım küreyi) bilerek, deneyin

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

(Boylam,enlem) olan bir dosyayı bir Mercator üzerinde cm cinsinden (x,y) konumlarına dönüştürmek için
derece başına 0.5 cm'lik belirli bir ölçek için ızgara, çalıştırın

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

Birkaç 2 sütunlu, ikili, çift duyarlıklı dosyayı (enlem, boylam) ile dönüştürmek için
için bir Enine Merkatör ızgarasında (merkezi boylam 75W) inç cinsinden (x,y) konumlarına
ölçek = 1:500000 ve harita alanının dışına düşecek noktaları bastırın, çalıştırın

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

NAD27'den old.dat dosyasındaki jeodezik koordinatları (boylam, enlem, yükseklik) dönüştürmek için
CONUS verisi (Clarke-131 elipsoidini kullanan Veri Kimliği 1866) WGS 84'e, çalıştırın

gmt harita projesi eski.dat -Th131 > yeni.dat

quakes.dat girdi dosyasındaki her nokta arasındaki en yakın mesafeyi (km olarak) hesaplamak için
ve multisegment ASCII dosyası coastline.xy'de verilen çizgi parçaları, çalıştırın

gmt harita projesi quakes.dat -Lcoastline.xy/k > quake_dist.dat

KISITLAMALAR

ile ayarlanan dikdörtgen giriş bölgesi -R genel olarak dikdörtgen olmayan bir şekilde eşlenecek
Kafes. Meğer ki -C ayarlandığında, bu ızgaradaki en soldaki nokta x değeri = 0.0'a sahiptir ve
en alt noktanın y değeri = 0.0 olacaktır. Bu nedenle, bir haritayı dijitalleştirmeden önce, aşırı
harita koordinatları harita projesi uygun ölçeği kullanarak ve (x,y) değerlerinin ne olduğunu görün
üzerine haritalanırlar. olması için sayısallaştırmayı ayarlarken bu değerleri kullanın.
ters dönüşüm doğru çalışıyor veya alternatif olarak awk ölçeklendirmek ve kaydırmak için
(x,y) dönüştürmeden önceki değerler.

Bazı izdüşümler için, kullanıcının bir seçim yapmasına rağmen küresel bir çözüm kullanılabilir.
elipsoid. Bu, kullanıcılar -R ayar, etki alanını aşan bir bölge anlamına gelir
burada elipsoidal seri açılımları geçerlidir. Bunlar koşullar: (1)
Lambert Konformal Konik (-JL) ve Albers Eşit-Alan (-JB) küresel çözümü kullanacak
harita ölçeği 1.0E7'yi aştığında. (2) Enine Merkatör (-JT) ve UTM (-JU) olacak
batı veya doğu sınırı verildiğinde küresel çözümü kullanın -R fazla
Merkezi meridyenden 10 derece ve (3) Cassini için aynı (-JC) ancak bir limit ile
sadece 4 derece.

ELİPSOİDLER VE KÜRESELLER

GMT, uygulandıkları ve kullanıcı bir
referans şekli olarak elipsoid (bkz. PROJ_ELLIPSOID). Kullanıcının birkaç şeyden haberdar olması gerekir.
olası tuzaklar: (1) Enine Mercator, Albers ve
Lambert'in uyumlu koniği, haritalanan alanlar olduğunda elipsoid ifadeleri kullanırız.
küçük ve küresel ifadelere geçin (ve uygun yardımcı
enlemler) daha büyük haritalar için. Elipsoid formüller aşağıdaki gibi kullanılır: (a) Enine
Merkatör: Tüm noktalar merkezi meridyenin 10 derece içinde olduğunda, (b) Konik çıkıntılar
uzunlamasına menzil 90 dereceden az olduğunda, (c) Tüm noktalar olduğunda Cassini projeksiyonu
4 derece merkezi meridyen içindedir. (2) Bazı tarihsel bilgileri eşleştirmeye çalıştığınızda
veriler (örneğin, belirli bir izdüşüm ve belirli bir referansla elde edilen koordinatlar)
elipsoid) GMT'nin biraz farklı sonuçlar verdiğini görebilirsiniz. bir ihtimal
Bu uyumsuzluğun kaynağı, eski hesaplamaların genellikle daha az anlamlı basamaklar kullanmasıdır. İçin
Örneğin, Snyder'ın örnekleri genellikle Clarke 1866 elipsoidini kullanır (kendisi tarafından
bir düzleştirme f = 1/294.98). f'den eksantrikliğin karesini 0.00676862818 olarak alıyoruz
(GMT bunu kullanır), Snyder ise yuvarlar ve 0.00676866 kullanır. Bu fark,
birkaç on cm'lik tutarsızlıklar verin. Yansıtılan koordinatları yeniden oluşturmanız gerekirse
bu biraz farklı eksantriklikle, kendi elipsoidinizi
Clarke 1866 ile aynı parametreler, ancak f = 1/294.97861076 ile. Ayrıca, yaşlı olduğunu unutmayın
veriler farklı verilere atıfta bulunulabilir ve hangi verilerin kullanıldığını bilmiyorsanız ve
tüm verileri ortak bir veriye dönüştürün, onlarca yüzlerce uyumsuzluk yaşayabilirsiniz
metre. (3) Son olarak, PROJ_SCALE_FACTOR'un bazıları için belirli varsayılan değerlere sahip olduğunu unutmayın.
projeksiyonlar, üretilen sonuçları eşleştirmek için ayarı geçersiz kılmanız gerekebilir.
diğer ayarlarla.

onworks.net hizmetlerini kullanarak mapprojectgmt'yi çevrimiçi kullanın