1. 1. GMT 5.4 script

1. 1. ASTER grid 2019, near dyatlov pass, from earth explorer,
   2. translated to .nc file
2. gdal_translate -of netcdf -ot float32 ASTGTM2_N61E059_dem.tif aster.nc

grid=aster.nc grid2=viher5c.nc outi=out.ps proj=-JQ18 rect=-R59.41/59.47/61.73/61.79

1. rect=-R59.425/59.46/61.752/61.773

gmt makecpt -Ccopper -T000/1100/10 -Z > palette.cpt

echo "0 120/160/180 200 220/255/180" > palette2.cpt echo "200 220/255/180 1000 255/210/150" >> palette2.cpt echo "B black" >> palette2.cpt echo "F white" >> palette2.cpt echo "N 127.5" >> palette2.cpt

1. 1. create hillshade

1. gmt grdgradient aster.nc -Ghillshade-grad.nc -A270 -Ne0.6 -V
2. gmt grdhisteq hillshade-grad.nc -Ghillshade-hist.nc -N -V
3. gmt grdmath hillshade-hist.nc 5 DIV = hillshade-int.nc

1. 1. resample landsat satellite image

1. 1. cdo remapbil,hillshade-int.nc viher3.nc viher3c.nc

1. 1. ASTER image

1. gmt grdimage $proj $rect $grid -Cpalette.cpt -Ihillshade-int.nc -K > $outi
   1. landsat image

1. gmt grdimage $proj $rect band1.nc band2.nc band3.nc -K > $outi

1. 1. LANDSAT image

gmt grdimage $proj $rect -Cpalette2 -Ihillshade-int.nc viher5c.nc -K > $outi

gmt grdcontour $proj $rect -C20 -Wa0.8,127/0/0 -Wc0.3,150/0/0 -A100+f16p -B1m1m $grid -O -K >> $outi

1. 1. legend

gmt pscoast $proj $rect -W1 -N1 -S0 -Lf59.435/61.733/61.733/2 -O -K >> $outi

gmt psxy $proj $rect -W2.0p,0/0/32 -O -K<< EOF >>$outi 59.441314557592591,61.759064921390738 59.441982532345719,61.759290837580707 59.442422103032989,61.759667163104425 59.443146691171336,61.75991308684398 59.443658407546351,61.760124073132118 59.443995403063568,61.760497088033603 59.444296982921308,61.760558136817913 59.444428697732107,61.760694171952522 59.444667449146017,61.7607524006932 59.444938752938853,61.760831763632318 59.445146458343864,61.761023116660844 59.445419551132879,61.761311686096548 59.445939751565362,61.76153319994026 59.446367490676181,61.761763244955851 59.446790443785453,61.76203931549864 59.446864917478422,61.762291219425045 59.447281453581873,61.762615554390429 59.447625999543867,61.76271605886884 59.448055718164476,61.762729131747598 59.448078937819268,61.762814392520333 59.448344288358363,61.763077774503195 59.448700220471643,61.763353474798038 59.449407086089693,61.76362214137059 59.450152761536614,61.763837876724601 59.450500171736181,61.764111788680943 59.450995064797468,61.764342223351775 59.45164301419193,61.764471932961143 59.452064083115026,61.764478860006975 59.452138160505356,61.764574053189904 59.4526521901617,61.764891504405071 59.452822552036331,61.76509908635748 59.452821873638271,61.765171126657258 59.45291846123699,61.765271099091912 59.45308787813476,61.765320628243366 59.453272619735323,61.765294730356082 59.453427776796467,61.765297771711523 59.453567820249468,61.765394886824531 59.453676508269496,61.765446090648943 59.453797204676718,61.765447222359555 59.453822103848232,61.765548737836959 59.453790995432577,61.765615381407571 59.453825548356647,61.765676120897808 59.453869605600204,61.76577425271725 59.453739407455174,61.765868632413948 59.453710550918309,61.765920691218675 59.453702409566901,61.766005105044663 59.453638371268752,61.766083839547271 59.453486338947748,61.766149184540708 59.453599418724359,61.766343280391069 59.453954524431218,61.766399251446025 59.454184704861831,61.76644457976964 59.454260512248155,61.766527661514061 59.454152274608987,61.766642818170297 59.454028080024777,61.766745949142695 59.453918029275627,61.766837903885687 59.454003965944317,61.766893825119936 59.454607269774925,61.766800144682321 59.45451715455544,61.767130475980615 59.4545956591841,61.767215389611586 59.454748614307036,61.767282668763457 59.454773419365118,61.767393308255976 59.454919344833847,61.7674977121086 59.455195423020072,61.767472199880864 59.455394677751073,61.767459039733765 59.455643795650658,61.767474631943735 59.455731447237177,61.767585359165516 59.455726375353898,61.767651071373358 59.455846596297327,61.767746591933033 59.456033043706917,61.767698334151447 59.456250366145703,61.767730104176948 59.456164182356005,61.7679312118076 59.456230705720493,61.768024029409524 59.456212344837049,61.76814403545152 59.456003055211148,61.768314614517827 59.4557078277779,61.768472958382787 59.455331844879403,61.768617825070315 59.45523223679259,61.768843872040897 59.45534455450256,61.769013869328518 59.455270580940649,61.769268500888543 59.455641791499154,61.769554729890558 59.455535843082714,61.769713035198521 59.455099057205722,61.770030093133208 > 59.453535741666641,61.765312349409591 59.453604703093895,61.765190748522365 59.453786510493011,61.765083976615358 59.453936971788821,61.765039488211407 59.45414385607058,61.764837806639996 59.454162663732561,61.764742897207753 59.454231625159807,61.764659851214397 59.454413432558937,61.764579770937054 59.454425971000234,61.764484860708841 > 59.453911894906192,61.765457676668788 59.453943241009483,61.765380564331025 59.45392443334751,61.765294554187818 59.454043548540035,61.765244134336911 59.454250432821794,61.765226339075681 59.454375817234975,61.765140328501332 59.454444778662221,61.765009829239695 59.454538816972104,61.764950511210522 59.454601509178694,61.764849670298439 59.454726893591889,61.764748829055861 59.454896162549687,61.764731033508141 > 59.453899356464881,61.765775021562888 59.454150125291235,61.765724602499354 59.454281778925079,61.765712739178291 59.454469855544851,61.765605969082799 59.454526278530793,61.765490301061199 59.454645393723304,61.765282690700943 59.454814662681109,61.76520557792459 59.455071700728134,61.765155157927865 59.455197085141322,61.76509584017888 59.455234700465276,61.765045420002323 > EOF

1. 1. djatlov pass

gmt psxy $proj $rect -Sc2.0 -W3p,180/0/0,- -O -K<< EOF >>$outi 59.462778 61.754722 EOF

1. 1. tent

gmt psxy $proj $rect -St1.5 -G127/127/200 -W3p,0/0/180,- -O -K<< EOF >>$outi 59.43123 61.75877 EOF

1. 1. cedar

gmt psxy $proj $rect -St1.8 -G40/120/40 -W3p,0/42/0,-. -O -K<< EOF >>$outi 59.4549787 61.76477 EOF

gmt psxy $proj $rect -S+0.75 -W4.0 -G0/0/0 -O -K<< EOF >>$outi 59.44493692006586,61.762129336012116,k 59.44836167718233,61.76220595524,s 59.4498444493604,61.763239014956255,d 59.454716619139845,61.76479863131844,kr2 59.454716619139846,61.76479863131845,kr2b 59.45401968997703,61.76506217040845,zo4 59.45401968997704,61.76506217040846,zo4 59.45401968997705,61.76506217040847,zo4 59.45401968997706,61.76506217040848,zo4 EOF

echo "59.43123 61.7610 30 0 1 rm Tent" |gmt pstext $proj $rect -F+f30p,Sans-Serif,DarkBlue+jBL -K -O >> $outi echo "59.465 61.7627 30 0 1 rm Cedar" |gmt pstext $proj $rect -F+f30p,Sans-Serif,DarkBlue+jBL -K -O >> $outi echo "59.453 61.7655 20 0 1 rm Ravine" |gmt pstext $proj $rect -F+f30p,Sans-Serif,DarkBlue+jBL -K -O >> $outi

echo "59.462778 61.754722 30 1 4 rm Dyatlov Pass" |gmt pstext $proj $rect -F+f30p,Sans-Serif,0/32/0+jBL -K -O >> $outi echo "59.435 61.75 30 0 5 rm Holat-Sayakhyl" |gmt pstext $proj $rect -F+f30p,Italic,20/20/20+jBL -K -O >> $outi

ps2pdf out.ps out.pdf pdf2svg out.pdf out.svg

1. svgo out.svg