img2grdgmt - Online sa Cloud

PROGRAMA:

NAME

img2grd - I-extract ang subset ng img file sa Mercator o Geographic na format

SINOPSIS

img2grd imgfile grdfile rehiyon uri [ ] [ [minlat/maxlat] ] [ ] [ minuto ] [] [ navg
] [ [sukatan] ] [ [antas] ] [ maxlon ] [ -n]

tandaan: Walang puwang ang pinapayagan sa pagitan ng flag ng opsyon at ng mga nauugnay na argumento.

DESCRIPTION

img2grd nagbabasa ng img format file, nag-extract ng subset, at nagsusulat nito sa grid file. Ang -M
ang opsyon ay nagdidikta kung ang Spherical Mercator projection ng img file ay
napanatili o kung ang isang Geographic na grid ay dapat isulat sa pamamagitan ng pag-undo sa Mercator projection. Kung
pinili ang geographic na grid maaari ka ring humiling ng resampling sa eksaktong -R binigay

KAILANGAN MGA PANGANGATWIRANG

imgfile
Isang Mercator img format na file gaya ng marine gravity o seafloor topography field
tinatantya mula sa satellite altimeter data nina Sandwell at Smith. Kung naitakda ng user
isang variable ng kapaligiran $GMT_DATADIR, Pagkatapos img2grd susubukang hanapin imgfile in
$GMT_DATADIR; kung hindi ay susubukan nitong buksan imgfile direkta.

-Ggrdfile
grdfile ay ang pangalan ng output grid file.

-R[yunit]kanluran/silangan/timog/hilaga[/zmin/zmax][r]
kanluran, silangan, timog, at hilaga tukuyin ang rehiyon ng interes, at maaari mong tukuyin
ang mga ito sa decimal degrees o sa [+-]dd:mm[:ss.xxx][W|E|S|N] na format. Idugtong r kung mas mababa
kaliwa at kanang itaas na mga coordinate ng mapa ang ibinibigay sa halip na w/e/s/n. Ang dalawa
mga shorthand -Rg at -Rd nakatayo para sa pandaigdigang domain (0/360 at -180/+180 sa longitude
ayon sa pagkakabanggit, na may -90/+90 sa latitude). Bilang kahalili para sa paglikha ng grid, bigyan
Rcodelon/lat/nx/ny, Kung saan code ay isang 2-character na kumbinasyon ng L, C, R (para sa kaliwa,
gitna, o kanan) at T, M, B para sa itaas, gitna, o ibaba. hal, BL para sa kaliwang ibaba.
Ito ay nagpapahiwatig kung aling punto sa isang hugis-parihaba na rehiyon ang lon/lat tumutukoy ang coordinate
sa, at ang mga sukat ng grid nx at ny na may mga grid spacing sa pamamagitan ng -I ay ginagamit upang lumikha
ang kaukulang rehiyon. Bilang kahalili, tukuyin ang pangalan ng isang umiiral na grid file
at ang -R Ang mga setting (at grid spacing, kung naaangkop) ay kinopya mula sa grid.
paggamit -Ryunit inaasahan ang inaasahang (Cartesian) na mga coordinate na katugma sa napili -J
at inversely project namin upang matukoy ang aktwal na rectangular geographic na rehiyon. Para sa
pagtingin sa ibang anggulo (-p), opsyonal na dugtungan /zmin/zmax. Sa kaso ng pananaw na pananaw
(-p), isang z-range (zmin, zmax) ay maaaring idagdag upang ipahiwatig ang ikatlong dimensyon. Ito
kailangang gawin lamang kapag ginagamit ang -Jz opsyon, hindi kapag ginagamit lamang ang -p pagpipilian.
Sa huling kaso, ang isang perspective view ng eroplano ay naka-plot, na walang pangatlo
sukat

OPSYONAL MGA PANGANGATWIRANG

-C Itakda ang x at y Mercator coordinate na may kaugnayan sa projection center [Default is
kaugnay sa ibabang kaliwang sulok ng grid]. Nangangailangan -M.

-D[minlat/maxlat]
Gamitin ang pinalawak na hanay ng latitude -80.738/+80.738. Bilang kahalili, idagdag
minlat/maxlat bilang lawak ng latitude ng input img file. [Ang default ay
-72.006/72.006]. Karaniwang hindi kinakailangan dahil matutukoy natin ang lawak mula sa
inspeksyon ng laki ng file.

-E Maaaring gamitin kapag -M ay hindi nakatakda upang pilitin ang huling grid na magkaroon ng eksaktong pareho
rehiyon gaya ng hiniling sa -R. Bilang default, ang huling rehiyon ay isang direktang projection ng
ang orihinal na rehiyon ng Mercator at karaniwang lalampas nang bahagya sa
hiniling na hanay ng latitude, at higit pa rito ang pagtaas ng grid sa latitude ay hindi
tumugma sa pagtaas ng longitude. Gayunpaman, ang sobrang resampling ay nagpapakilala ng maliit
interpolation error at dapat lamang gamitin kung ang output grid ay dapat tumugma sa
hiniling na rehiyon at mayroong x_inc = y_inc. Sa kasong ito ang rehiyon na itinakda ng -R dapat
ibinibigay sa multiple ng increment (.hal, -R0/45/45/72).

-I Ipahiwatig minuto bilang ang lapad ng isang input img pixel sa mga minuto ng longitude.
[Default ay 2.0]. Karaniwang hindi kinakailangan dahil matutukoy natin ang laki ng pixel mula sa
inspeksyon ng laki.

-M Mag-output ng Spherical Mercator grid [Ang default ay isang geographic na lon/lat grid]. Ang
Ang Spherical Mercator projection ng img file ay napanatili, upang ang rehiyon -R
itinakda ng gumagamit ay bahagyang binago; ang binagong rehiyon ay tumutugma sa mga gilid
ng mga pixel [o mga pangkat ng navg mga pixel]. Ang header ng grid file ay nakatakda upang ang x at
Ang mga haba ng y axis ay kumakatawan sa distansya mula sa kanluran at timog na mga gilid ng imahe,
sinusukat sa mga default na unit ng user, na may -Jm1 at ang inayos -R. Sa pamamagitan ng pagtatakda ng
default PROJ_ ELLIPSOID = Sphere, ang user ay maaaring gumawa ng mga overlay gamit ang adjusted -R
para magkatugma sila. Tingnan mo HALIMBAWA sa ibaba. Ang inayos -R ay nakasulat din sa grid
header remark, para mahanap ito sa ibang pagkakataon. Tingnan mo -C upang itakda ang mga coordinate na may kaugnayan sa
sentro ng projection.

-Nnavg I-average ang mga value sa input img pixels sa navg by navg mga parisukat, at lumikha
isang output pixel para sa bawat parisukat. Kung gagamitin sa -T3 mag-uulat ito ng average
hadlang sa pagitan ng 0 at 1. Kung ginamit sa -T2 ang output ay magiging average na halaga ng data
o NaN ayon sa kung ang average na pagpilit ay > 0.5. navg dapat pantay na hatiin
sa mga sukat ng imgfile sa mga pixel. [Default 1 walang pag-average].

-S[sukatan]
I-multiply ang mga halaga ng img file sa pamamagitan ng sukatan bago mag-imbak sa grid file. [Ang default ay
1.0]. Para sa mga kamakailang img file: img topo file ay naka-imbak sa (naitama) metro [-S1];
free-air gravity file sa mGal*10 [-S0.1 upang makakuha ng mGal]; vertical deflection files in
micro-radians*10 [-S0.1 para makakuha ng mga micro-radian], papasok ang mga file ng vertical gravity gradient
Eotvos*50 [-S0.02 para makakuha ng Eotvos, o -S0.002 upang makakuha ng mGal/km]). Kung hindi sukatan ay ibinigay
sinusubukan naming tukuyin ang sukat sa pamamagitan ng pagsusuri sa pangalan ng file para sa mga pahiwatig.

-Turi uri pinangangasiwaan ang pag-encode ng impormasyon ng hadlang. uri = 0 ay nagpapahiwatig na hindi
ang naturang impormasyon ay naka-encode sa img file (ginamit para sa pre-1995 na bersyon ng
gravity data) at nakukuha ang lahat ng data. uri Ang > 0 ay nagpapahiwatig na ang impormasyon ng hadlang ay
naka-encode (1995 at mas bago (kasalukuyang) bersyon ng mga img file) upang ang isa ay maaaring
gumawa ng grid file tulad ng sumusunod: -T1 nakakakuha ng mga halaga ng data sa lahat ng punto, -T2 nakakakuha ng data
mga halaga sa mga napiling punto at NaN sa mga interpolated na punto; -T3 nakakakuha ng 1 sa
constrained points at 0 sa interpolated points [Default is 1].

-V[antas] (higit pa ...)
Piliin ang antas ng verbosity [c]. Partikular na inirerekomenda dito, dahil ito ay kapaki-pakinabang upang makita
kung paano inaayos ang mga coordinate.

-Wmaxlon
Ipahiwatig maxlon bilang ang maximum longitude na lawak ng input img file. Mga bersyon
mula noong 1995 ay nagkaroon maxlon = 360.0, habang mayroon ang ilang naunang mga file maxlon = 390.0.
[Ang default ay 360.0].

-^ or m -
Mag-print ng maikling mensahe tungkol sa syntax ng command, pagkatapos ay lumabas (TANDAAN: sa Windows
gamitin lang -).

-+ or m +
Mag-print ng malawak na paggamit (tulong) na mensahe, kasama ang paliwanag ng anuman
opsyong tukoy sa module (ngunit hindi ang mga karaniwang opsyon sa GMT), pagkatapos ay lalabas.

-? or hindi argumento
Mag-print ng kumpletong mensahe ng paggamit (tulong), kasama ang paliwanag ng mga opsyon, pagkatapos
labasan.

--bersyon
I-print ang bersyon ng GMT at lumabas.

--show-datadir
I-print ang buong path sa GMT share directory at lumabas.

HEOGRAPHIC HALIMBAWA

Ang -M dapat na hindi kasama ang opsyon kung kailangan mong nasa geographic ang output grid
mga coordinate. Upang kunin ang data sa rehiyon -R-40/40/-70/-30 mula sa world_grav.img.7.2 at
reproject upang magbunga ng mga geographic na coordinate, maaari mong subukan

img2grd world_grav.img.16.1 -Gmerc_grav.nc -R-40/40/-70/-30 -V

Dahil ang latitude spacing sa img file ay katumbas ng distansya sa mga unit ng Mercator, ang
ang resultang grid ay hindi tutugma sa tinukoy -R eksakto, at ang latitude spacing ay hindi
katumbas ng longitude spacing. Kung kailangan mo ng eksaktong tugma sa iyong -R at ang parehong espasyo
sa longitude at latitude, gamitin ang -E opsiyon:

img2grd world_grav.img.16.1 -Gmerc_grav.nc -R-40/40/-70/-30 -E -V

MERCATOR HALIMBAWA

Dahil ang mga img file ay nasa isang Mercator projection, HINDI mo dapat i-extract ang isang geographic na grid
kung ang plano mo ay gumawa ng mapa ng Mercator. Kung ginawa mo iyon tinapos mo ang pag-project at
reprojection ang grid, nawawala ang maikling-wavelength na detalye. Mas magandang gamitin -M at i-plot ang grid
gamit ang isang linear na projection na may parehong sukat ng nais na Mercator projection (tingnan ang GMT
Halimbawa 29). Upang kunin ang data sa rehiyon -R-40/40/-70/-30 mula sa world_grav.img.7.2, tumakbo

gmt img2grd -M world_grav.img.7.2 -Gmerc_grav.nc -R-40/40/-70/-30 -V

Tandaan na ang -V ang opsyon ay nagsasabi sa amin na ang hanay ay naayos sa
-R-40/40/-70.0004681551/-29.9945810754. Magagamit din natin grdinfo upang mahanap na ang grid file
ipinapakita ng header ang rehiyon nito -R0/80/0/67.9666667 Ito ang hanay ng x,y na makukuha natin
mula sa isang Spherical Mercator projection gamit ang -R-40/40/-70.0004681551/-29.9945810754 at
-Jm1. Kaya, para kunin ang ship.lonlatgrav at gamitin ito para sampolan ang merc_grav.nc, magagawa natin ito:

gmt set PROJ_ELLIPSOID Sphere

gmt mapproject -R-40/40/-70.0004681551/-29.9945810754 -Jm1i ship.lonlatgrav | \
gmt grdtrack -Gmerc_grav.nc | gmt mapproject \
-R-40/40/-70.0004681551/-29.9945810754 -Jm1i -I > ship.lonlatgravsat

Inirerekomenda na gamitin ang paraan sa itaas ng pag-project at pag-unproject ng data sa naturang
isang aplikasyon, dahil mayroon lamang isang interpolation na hakbang (sa grdtrack). Kung isa
unang sumusubok na i-convert ang grid file sa lon,lat at pagkatapos ay i-sample ito, mayroong dalawa
mga hakbang sa interpolation (sa conversion at sa sampling).

Upang makagawa ng lon,lat grid mula sa itaas na grid na magagamit namin

gmt grdproject merc_grav.nc -R-40/40/-70.0004681551/-29.9945810754 -Jm1i -I -D2m -Ggrav.nc

Sa ilang mga kaso hindi ito magiging madali gaya ng -R sa dalawang coordinate system ay maaaring hindi magkatugma
mabuti. Kapag nangyari ito, maaari rin naming gamitin (sa katunayan, maaaring palaging mas mahusay na gamitin)

gmt grd2xyz merc_grav.nc | gmt mapproject \
-R-40/40/-70.0004681551/-29.994581075 -Jm1i -Ako | \
gmt surface -R-40/40/-70/70 -I2m -Ggrav.nc

Upang makagawa ng isang mapa ng Mercator ng rehiyon sa itaas, ipagpalagay na ang halaga ng aming gmt.conf ay para sa
Ang PROJ_LENGTH_UNIT ay pulgada. At dahil ang nasa itaas na merc_grav.nc file ay inaasahang may -Jm1i
ito ay 80 pulgada ang lapad. Maaari tayong gumawa ng mapa na 8 pulgada ang lapad sa pamamagitan ng paggamit -Jx0.1i sa anumang mga programa sa mapa
inilapat sa grid na ito (hal, grdcontour, grdimage, grdview), at pagkatapos ay para sa mga overlay na
magtrabaho sa lon,lat (hal., psxy, pscoast) maaari nating gamitin ang nababagay sa itaas -R at -Jm0.1 upang makuha
magkatugma ang dalawang sistema.

Gayunpaman, maaari tayong maging mas matalino kaysa dito. Napagtatanto na ang input img file ay may mga pixel 2.0
minuto ang lapad (o sinusuri ang nx at ny gamit ang grdinfo merc_grav.nc) napagtanto namin iyon
Ginamit ng merc_grav.nc ang buong resolution ng img file at mayroon itong 2400 by 2039 pixels, at
sa 8 pulgada ang lapad ito ay 300 pixels bawat pulgada. Nagpasya kami na hindi namin kailangan ng ganoon karami at kami
ay masisiyahan sa 100 pixels per inch, kaya gusto naming i-average ang data sa 3 by 3
mga parisukat. (Kung gusto natin ng contour plot, malamang na pipiliin nating i-average ang data nang higit pa
(hal., 6 by 6) para makakuha ng makinis na contours.) Dahil ang 2039 ay hindi nahahati ng 3 makakakuha tayo ng
iba ang inayos -R sa pagkakataong ito:

gmt img2grd -M world_grav.img.7.2 -Gmerc_grav_2.nc -R-40/40/-70/-30 -N3 -V

Sa pagkakataong ito, nakita namin ang naayos na rehiyon ay -R-40/40/-70.023256525/-29.9368261101 at ang
ang output ay 800 by 601 pixels, mas magandang sukat para sa amin. Ngayon ay maaari tayong lumikha ng isang artipisyal
file ng pag-iilaw para sa paggamit nito grdgradient:

gmt grdgradient merc_grav_2.nc -Gillum.nc -A0/270 -Ne0.6

at kung mayroon din tayong CPT file na tinatawag na "grav.cpt" maaari tayong gumawa ng color shaded na relief map
ganito:

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

Ipagpalagay na gusto mong makuha lamang ang mga napiling halaga ng data mula sa isang img file, sa lat/lon
mga coordinate. Tapos tumakbo img2grd sa -T2 opsyon, gamitin grd2xyz upang dump ang mga halaga, pipe
sa pamamagitan ng grep -v NaN upang maalis ang mga NaN, at i-pipe through mapa proyekto kasama ang kabaligtaran
projection tulad ng nasa itaas.

Gumamit ng img2grdgmt online gamit ang mga serbisyo ng onworks.net