Download 本文表示 - 寒地土木研究所
Transcript
技術ニュース GMT を用いたデータマッピング手法について 村上 泰啓* 㧝㧚ߪߓߦ ࡂ ࡢ ࠗ ᄢ ቇ ߩ WEB ࠨ ࠗ ࠻ ߢ ή ᢱ 㐿 ߐ ࠇ ߡ ࠆ GMT1)(Generic Mapping Tools)ࠍ↪ߔࠆߣޔ࿑ ߩ᭽ߥ 㜞ຠߥ࿑ࠍࠆߎߣ߇ߢ߈߹ߔޕGMT ߪ⇇ਛߢᄙߊ ߩ⎇ⓥ⠪߿ᐭ㑐ଥ⠪߇↪ߒߡࠆ࠷࡞ߢޔౝᄖߩ WEB ࠨࠗ࠻ߢ↪ᣇᴺ߇⸃⺑ߐࠇࠆߥߤޔ ᵴ⊒ߦ↪ߐ ࠇߡࠆࡈߩ࠷࡞ߩ৻ߟߢߔࠦޔࠄ߇ߥߒ߆ߒޕ ࡑࡦ࠼ࡊࡠࡦࡊ࠻ߢ㚟േߔࠆߚዋޘൎᚻ߇ᖡߊޔ ߒ߆߽ WINDOWS ߩ⸃⺑߇ᄙߊߥߚ╩ޔ⠪ߪ߃ ࠆࠃ߁ߦߥࠆ߹ߢᤨ㑆߇߆߆ࠅ߹ߒߚޕ GMT ࠍዉߔࠆ ࿑ )/6 ߢឬ↹ߒߚᣣᧄㄝߩᒻ㧔ᒻ࠺࠲ߪ ߹ߢߪޔᵩ᳓ᤨߩ㒠㔎ಽᏓ࿑߿⺞ᩏ⟎࿑ࠍ⧰ഭߒߡ )'$%1 ߩ ⑽⸃ᐲߩ߽ߩࠍ↪㧕 ࿑ߒߡ߹ߒߚ߇ޔGMT ߢߎ߁ߒߚ࿑߇◲නߦឬߌࠆߎ ߣࠍ⍮ࠅޔ ⺰ᢥ߿⊒⾗ᢱߦᄙ↪ߔࠆࠃ߁ߦߥࠅ߹ߒߚޕ Ꮢ⽼ߩ GIS ߦߪᢜ߹ߖࠎ߇ޔGMT ࠍ↪ߔࠇ߫࿑ߛߌߢߥߊ߫߃ޔ᳓᷷ಽᏓ߿ᔕജಽᏓߥߤᄌ㊂ߩ ੑᰴర⊛ߥಽᏓ߿ࡌࠢ࠻࡞࿑ߥߤ߽⍴ᤨ㑆ߦ࿑ߢ߈߹ߔߩߢޔ㜞ଔߥ࠰ࡈ࠻࠙ࠚࠕࠍዉߖߕߣ߽ᣢሽߩ PC ⅣႺߢ߽චಽߦ⋡ᮡࠍ㆐ᚑߢ߈ࠆࠊߌߢߔޕGMT ߩࠗࡦࠬ࠻࡞ᣇᴺ߿◲ᤃߥࠬࠢࡊ࠻ߪޔWEB ߢ߽ᣢߦ 㐿()ߐࠇߡ߹ߔ߇ޔWINDOWS ߩ gawk ߦࠃࠆࡈࠔࠗ࡞ᠲߩᣇᴺ߿ᣂߒᯏ⢻(ᣣઃゲޔಠ)ߦ㑐ߔ ࠆ⺑ߪ߶ߣࠎߤήߩ߇ታᖱߢߔޔߪߢߎߎޕೋᔃ⠪ะߌ߆ࠄᵩ᳓ᤨߩ╬㔎㊂✢࿑߿᳓ᷓᣇะߩ᳓⾰㗄⋡ߩᄌ ൻࠍᤨ㑆ゲߢߔࠆ႐วߥߤޔ႐ߩᛛⴚ⠪߇ᔕ↪ߒ߿ߔࠃ߁◲නߥ⺑ߣࠨࡦࡊ࡞ࠬࠢࡊ࠻ࠍ૬⸥ߒ ߹ߒߚߩߢޔᓮᵴ↪㗂ߌࠇ߫ᐘ↟ߢߔޕ 㧞㧚GMT ߩ⸳ቯ ࡂࡢࠗᄢቇߩ WEB ࠨࠗ࠻ http://gmt.soest.hawaii.edu/ߦߡࠍ࡞ࠗࠔࡈߩ ޔᚻߒ߹ߔ ߪߢߎߎޕWindows ߩ߽ߩࠍㆬᛯߒ߹ߔޕԘ㨪Ԛࠍ㗅ߦታⴕߒߡߊߣޔc:¥programs¥GMT ߩਅߦࡈࠜ࡞࠳߇ࠄࠇޔGMT ᧄ ߩ࠶࠻ࠕ࠶ࡊ߇ቢੌߒ߹ߔ ࡦ࡚ࠫࡃޕ3 એ೨ߩ GMT ߢߪࡄࠬߩ⸳ቯࠍᚻേߢⴕ߁ᔅⷐ߇ࠅ߹ߒߚ߇ࡃޔ ࡚ࠫࡦ㧠એ㒠ߩ GMT ߢߪࡄࠬ߇⥄േ⊛ߦ⸳ቯߐࠇࠆࠃ߁ߦߥࠅ߹ߒߚޕᰴߦԛߩ gawk316.zip ࠍ⸃ಓߒߡᓧࠄ ࠇߚ gawk.exe ࠍ c:¥gawk ࡈࠜ࡞࠳ߦᩰ⚊ߒޔ ࡑࠗࠦࡦࡇࡘ࠲ࠍฝࠢ࠶ࠢߒߡࠪࠬ࠹ࡓߩࡊࡠࡄ࠹ࠖࠍ㐿߈ޔ ⚦⸳ቯ࠲ࡉߩⅣႺᄌᢙ࠲ࡉࠍࠢ࠶ࠢߒࡓ࠹ࠬࠪޔⅣႺᄌᢙߩ Path ࠍㆬᛯߒߡ✬㓸ࡏ࠲ࡦࠍࠢ࠶ࠢߔࠆߣޔ ࠪࠬ࠹ࡓᄌᢙߩ✬㓸࠙ࠖࡦ࠼࠙߇㐿ߊߩߢޔᄌᢙ୯ߩᨒߩᦨᓟߩㇱಽߦ ;c:¥gawk ࠍㅊടߒޔOK ࠍࠢ࠶ࠢߒ߹ ߔ ޕnetcdf-3.6.1-win32.zip ࠍ⸃ಓᓟ⚿ޔᨐࠍ c:¥netcdf ࡈࠜ࡞࠳ߦᩰ⚊ߒߡ߅߈߹ߔ ߽ࠇߎޕgawk ߣห᭽ߦࡄࠬ 㧝 ࡂࡢࠗᄢቇ 9'$ ࠨࠗ࠻ߢᚻน⢻ߥ )/6 㑐ㅪߩࡄ࠶ࠤࠫ৻ⷩ Ԙ GMT_basic_install.exe GMT ߩࠦࡑࡦ࠼ࡄ࠶ࠤࠫ⟲ߩࠗࡦࠬ࠻ ԙ GSHHS_ highfull_install.exe ో⇇ߩᶏጯ✢ޔ᳓ၞߩࠗࡦ࠺࠲ߩࠗࡦࠬ࠻ Ԛ GMT_pdf_install.exe ขᛒ⺑ᦠ ԛ gawk316.zip awk ߩ windows ߩࠗࡦࠬ࠻ Ԝ netcdf-3.6.1-win32.zip Netcdf ᒻᑼߩ࠺࠲ࠍ⺒ߺㄟߺὑߩࠗࡉࠗࡦࠬ࠻ 寒地土木研究所月報 №685 2010年6月 45 set waku=139/146/41/46 # set outfile=hok_rd.ps # set bp=a2.0f1.0g0.0 # set rd=..¥topodata¥hok_rd.dat # #1 set r=2 ( / (a / ( ) f / ) ( ) g ( )) ¥topodata¥ %rd% psbasemap -Jm%r% -R%waku% -B%bp% -G200 -P -K > %outfile% # -G200 gawk "{print $1,$2}" %rd% | psxy -Jm%r% -R%waku% -B%bp% -W1/255/0/0 -M -P -V -O -K >> %outfile% gawk %rd% 1,2 | psxy pscoast -Jm%r% -R%waku% -W1 -Df -P -V -S255 -O >> %outfile% NETCDF Ghostscript GMT 255 ) -S255 c:¥netcdf¥bin gawk ;c:¥netcdf¥bin GSview ( GSview 2010 2 GSview GMT US # (0 ps eps) Windows Ghostscript GSview WEB gs871w32.exe gsv49w32.exe URL (WEBサイトのURL http://auemath.aichi-edu.ac.jp/~khotta/ghost/#disk ) C:¥programs¥GMT¥share¥conf¥gmt.conf 2 SI GMT GMT hok_coast.ps ( ) ࡊ ࡠࡦࡊ࠻ GMT hok_rd.bat GMT 1 set waku set 2 3 5 gawk 46 1 ( ) 2 ( pscoast 4 psbasemap ) outfile psxy 寒地土木研究所月報 №685 2010年6月 # ᷹ⷰᚲฬ ⚻ᐲ ✲ᐲ ✚㔎㊂(mm) 141.67 42.85 173 ੑୀ 143.88 42.97 154 ⩵㕙ᴛ 142.02 43.21 137 ⟤ᵹᷰ 141.91 43.16 124 ⧃࠳ࡓ 142.1 43.28 122 ৼᧂጊ 141.98 43.08 117 ᳓ญ 142.08 43.26 116 ᩵ᴛ࠳ࡓ 142 43.24 113 ↱ੳ 141.78 43.03 113 㔎ᾍ 141.82 43.08 112 ධ㐳ᴧ 141.72 42.95 106 # ⚻ᐲ 139.992 139.993 > 139.986 139.986 139.986 139.986 139.986 139.986 139.987 139.987 139.987 > ✲ᐲ 41.551 41.549 41.573 41.572 41.572 41.571 41.571 41.571 41.571 41.571 41.570 ࿑ ✚㔎㊂࠺࠲ Ꮐޔᴡ✂࠺࠲ ฝߩ৻ㇱ 㧪ߪࠦࡑࡦ࠼ࠪࠚ࡞ߩ࠳ࠗࠢ࠻Ṷ▚ሶߣ߫ࠇ⚿ޔ ᨐࠍᣂⷙߦࡈࠔࠗ࡞ߦᦠ߈ㄟߢߔޕpsbasemap ߢ ߞߚ࿑ߦ〝߿ᶏጯ✢ߩࠗࡦࠍ㊀ߨߡߊߩߢߔ߇ޔ ឬ↹ࠦࡑࡦ࠼߇ታⴕߐࠇࠆㇺᐲࠗࡔࠫࠍㅊ⸥ߒߡߊ ࿑ ᐕ ⽕㔎ᤨߩ✚㔎㊂ಽᏓ࿑ ᔅⷐ߇ࠆߩߢޔpsbasemap ࠦࡑࡦ࠼ߩᰴ߆ࠄߪ࠳ޔ ࠗࠢ࠻Ṷ▚ሶߩ>>ࠍߞߡ㊀ߨᦠ߈ߒߡ߹ߔޕ4 ⴕ ⋡ߩ set ᢥߢ〝✂ߩࠗࡦ࠺࠲ࠍ⸥㍳ߒߚ hok_rd.dat ߣ߁ࡈࠔࠗ࡞ฬߩ೨ߦ..¥topodata¥ߣ߁⸥ภ߇ߞߡ ߹ߔ࠻࠶࠼ޕ㧞ߟߩ⸥ภ..ߪࡃ࠶࠴ࡈࠔࠗ࡞ߩࠆࡈࠜ࡞࠳ߩࡈࠜ࡞࠳ࠍ␜ߔ⸥ภߢޔᰴߩ¥topodata¥ߪ ࡈࠜ࡞࠳ߩਅߦࠆࡈࠜ࡞࠳ࠍ␜ߒ߹ߔ࠼ࡦࡑࠦޕਛޔ%ߢ߹ࠇߚᄌᢙ߇ࠅ߹ߔ߇ ߪࠇߎޔset ᢥߢቯ⟵ߐ ࠇߚᄌᢙฬߢߔޕset ᢥߢࡈࠔࠗ࡞ฬ߿᧦ઙࠍ೨ߦ⸳ቯߒߡ߅ߊߣޔset ᢥߩౝኈࠍᄌᦝߔࠆߛߌߢޔએ㒠ߩࠬ ࠢࡊ࠻ࠍߜߜୃᱜߔࠆᔅⷐ߇ήߊߥࠆߣ߁ࡔ࠶࠻߇ࠅ߹ߔޕ-Jޔ-Rޔ-Bޔ-W ߥߤߩ⸥ภߪઁߩ GMT ࠦࡑࡦ࠼ߢ߽ࠊࠇࠆㅢࠦࡑࡦ࠼ࠗࡦࠝࡊ࡚ࠪࡦߣ߫ࠇޔ-J ߪᛩᓇᴺޔ-R ߪឬ↹ߔࠆゲߩ▸࿐ޔ-B ߪ࿑ࠍ ࿐ࡈࡓ߿⚻✲✢ߩ᭽ޔ-W ߪ✢ߩᄥߐ߿⦡ࠍ⸳ቯߔࠆߚߩ߽ߩߢߔ⺑ߥ⚦ߩࡦ࡚ࠪࡊࠝߪߢߎߎޕ ߪ⋭߈߹ߔޕ 㧔㧞㧕ᵩ᳓ᤨߩ✚㔎㊂ಽᏓ࿑ࠍឬߊ ࿑ ߦ␜ߔᒻᑼߢ㔎㊂࠺࠲㧔᷹ⷰᚲฬ✲ޔᐲ⚻ᐲᐳᮡ✚ޔ㒠㔎㊂㧕 ޔᴡ✂(✲ᐲ⚻ޔᐲ)ߥߤߩ࠺࠲ࠍḰ ߒ߹ߔޕ࿖ᢙ୯ᖱႎߦߪో࿖ߩ〝✂߿ᴡᎹ✂ߩ࠺࠲߇߹ࠇߡ߹ߔ߇✲ࠍࠄࠇߎޔᐲ⚻ᐲਗ߮ߦᄌ឵ߔ set hokriv=..¥topodata¥hokkaido.riv # ർᶏౝߩᴡᎹࠗࡦ࠺࠲㧔⚻ᐲ✲ᐲਗ߮㧕 set grddata=raingrid.grd # ✚㔎㊂ಽᏓߩࠣ࠶࠼࠺࠲ࡈࠔࠗ࡞ set rdat=train.txt # ᷹ⷰὐᲤߩ࠺࠲ࡈࠔࠗ࡞㧔⚻ᐲ✲ޔᐲ✚ޔ㔎㊂ਗ߮ߩ࠺࠲㧕ᜰቯ set of=hok_200308.ps # ജࡈࠔࠗ࡞ฬߩ⸳ቯ set cp=GMT_saru_c.cpt # 㔎㊂ಽᏓߩࠞࡄ࠶࠻ set waku=139/146/41/46 # ឬ↹㗔ၞߩ⸳ቯ set bp=a2f1g0 # ࿑ᨒߩᢙ୯㑆㓒߇ 2 ᐲࡓࡈޔ㑆㓒߇ 1 ᐲߪ✢⚻✲ޔήߒޕ set r=2 # ⚻ᐲ 1 ᐲߚࠅߩឬ↹ࠨࠗ࠭߇ 2cm psbasemap -Jm%r% -R%waku% -B%bp% -P -K > %of% # ၮ࿑࿑ makecpt -CGMT_seis.cpt -T0/450/10 -I > %cp% # ᣢሽࠞࡄ࠶࠻ࠍరߦ✚㔎㊂ᦨᄢ୯ߦวࠊߖߡᣂߒࡄ࠶࠻ࠍ⸳ቯޕ surface %rdat% -G%grddata% -I11.25c/7.5c -R%waku% -T0.25 -Ll0.0 # ✚㒠㔎㊂ߩὐᖱႎࠍ㕙(ࠣ࠶࠼)ߦᄌ឵ޕ grdimage %grddata% -Jm%r% -R%waku% -B%bp% -C%cp% -K -O >> %of% # ✚㔎㊂ಽᏓࠍᣂߒࠞࡄ࠶࠻ߢឬ↹ gawk "{print $1,$2}" %hokriv% | psxy -Jm%r% -W1/0/0/255 -M -R -P -V -K -O >> %of% # ർᶏౝߩᴡᎹ✂ឬ↹ gawk "{print $1,$2}" %rdat% | psxy -R%waku% -Jm%r% -Sc0.05 -W1/0/0/0 -G255/0/0 -O -K >> %of% # ർᶏౝߩ㔎㊂᷹ⷰὐឬ↹ #╬㔎㊂✢ߩឬ↹ grdcontour %grddata% -Jm%r% -R%waku% -Wa1/255/0/0 -Wc0/255/0/0 -C10 -A50f8t -L10/500 -G5/5 -K -O >> %of% pscoast -Jm%r% -R -B -W2 -Dh -P -V -S255 -K -O >> %of% # ᶏጯ✢ߩឬ↹ –Dh ߪᦨ㜞⸃ᐲޔ-S255 ߪ᳓ㇱߩႣࠅߟ߱ߒޕ psscale -D4/-1/10/0.3h -C%cp% -Ba50f50 -O >> %of% # 㒠㔎㊂ಽᏓߩ⦡ࠬࠤ࡞ឬ↹ ࿑ ࿑ ࠍឬ↹ߔࠆߚߩࠬࠢࡊ࠻㧔ࠣ✂ដߌㇱಽߪ⺑↪ߢࠅޔታ㓙ߦߪਇⷐ㧕 寒地土木研究所月報 №685 2010年6月 47 ࠆߎߣߢޔGMT ߢ↪ߢ߈ࠆࠃ߁ߦߥࠅ߹ߔޕ࿑ ฝ ߩᴡ✂࠺࠲ߩ࠺࠲ਗ߮ߪⴕޔᲤߦ⚻ᐲ✲ޔᐲߢᢙ ୯߇ਗ߮߇࠲࠺ߢࠢࡑ ޔಾࠄࠇߡ߹ߔޕᴡߢᢙ ߪⶄᢙߩࠗࡦߢ᭴ᚑߐࠇߡ߹ߔߩߢࠗࠆߥ⇣ߢ ޔ ࡦߩᖱႎߢࠆߣ߁ࠍ␜ߒߡ߹ߔޕ೨ㅀߩ〝 ࠗࡦ߽ห᭽ߥ࠺࠲᭴ㅧߦߥߞߡ߹ߔ ޔߡߐޕ࿑ ߩࡃ࠶࠴ࡈࠔࠗ࡞ࠍታⴕߔࠆߣޔ10~20 ⑽⒟ᐲߢ࿑ ߩ ࠃ߁ߥ╬㔎㊂✢࿑ࠍឬߊߎߣ߇ߢ߈߹ߔߩߎޕ࿑ߦ pstext ࠦࡑࡦ࠼ߢ⺑ࠍട߃ࠆߎߣ߽น⢻ߢߔ߇ޔ GSVIEW ߢ jpg ߥߤߩ↹ࡈࠔࠗ࡞ߦᄌ឵ߒޔᢥᦠࡈࠔ ࠗ࡞ߦ⾍ࠅઃߌߚᓟޔ ಠ߿ゲߩ⺑ޔ ࠺࠲ߩᚻవޔ 㓸⸘ᦼ㑆ߥߤߩㅊടᖱႎࠍ✬㓸ߒߚ߶߁߇ଢ߆߽ߒࠇ ߹ߖࠎޕ 㧔㧟㧕㔡ᄩಽᏓ࿑ߩឬ↹ ☨࿖⾰⺞ᩏᚲߩ WEB ࠨࠗ࠻ߢߪޔ1973 ᐕએ㒠ߦ⊒ ↢ߒߚࡑࠣ࠾࠴ࡘ࠼㧝એߩ㔡࠺࠲߇㐿ߐࠇߡ ߹ߔ(࠲࠺ߩߎޕ࿑)ߦߪ↢⊒ޔ⟎ࡘ࠴࠾ࠣࡑޔ ࠼ޔ㔡Ḯߩᷓߐߥߤ߇߹ࠇߡ߹ߔߩߢޔGMT ࠍ ߞߡ㔡ᄩಽᏓ࿑ࠍឬߡߺ߹ߔޔߪߢߎߎޕනߦ㔡Ḯ ⟎ࠍ␜ߔߛߌߢߥߊߩ࠼ࡘ࠴࠾ࠣࡑޔᄢ߈ߐߦࠃߞߡ 㔡Ḯ⟎ߩࡑࠢߩᄢ߈ߐ߿ޔ㔡Ḯᷓߐߦࠃߞߡ⦡ࠍᄌ ߃࠼ࡘ࠴࠾ࠣࡑޔᄢ߈ߐࠍ␜ߔಠ߽࿑ਛߦࠇߡߺ ߹ߔޕ࿁ߪ pslegend ߣ߁ࠦࡑࡦ࠼߇ߞߡ߹ߔޕ PDE PDE PDE PDE PDE PDE 1973 1973 1973 1973 1973 01 06 153931 01 07 121712.60 01 09 022114.80 01 09 170755.50 01 10 135240.60 38.00 46.43 5.27 36.85 37.81 141.69 51.40 -178.21 36.84 137.36 24 34 59 52 33 4.2 mbGS 4.9 mbGS 3.5 mbGS 5.1 mbGS 4.4 mbGS Preliminary Determination of Epicenters set waku=120/150/20/50 # set outfile=seis73-10.ps # set bp=a5.0f2.5g0 # set seis=..¥seisdata¥seis.txt # set r=0.5 # set cpt=depth.cpt # makecpt -T0/700/20 -I -Z > %cpt% # N1 S 0.5c c 0.10c 255/0/0 0.5p 0.9c mag2. S 0.5c c 0.15c 255/0/0 0.5p 0.9c mag3. S 0.5c c 0.20c 255/0/0 0.5p 0.9c mag4. S 0.5c c 0.25c 255/0/0 0.5p 0.9c mag5. S 0.5c c 0.30c 255/0/0 0.5p 0.9c mag6. S 0.5c c 0.35c 255/0/0 0.5p 0.9c mag7. S 0.5c c 0.40c 255/0/0 0.5p 0.9c mag8. S 0.5c c 0.45c 255/0/0 0.5p 0.9c mag9. 5F. ....... 3F. ....... 2F. ....... 2F. ....... 2F. ....... ( ) 1 0.5cm 0~700km 20km psbasemap -Jm%r% -R%waku% -Y5 -B%bp% -P -K > %outfile% # gawk "{print $7,$6,$8,$9/20 }" %seis% | psxy -J -R -B -C%cpt% -Sc -W1 -M -P -V -O -K >> %outfile% pscoast -J -R -W1 -Df -P -V -O -K >> %outfile% # psscale -D5/-1/8/0.2h -Ba100f50g50:Depth(km): -C%cpt% -O -K -V >> %outfile% # pslegend legend.txt -R -J -Dx0.5/18/3c/5c/TL -F -G255 -O -V -S > legend3.bat call legend3.bat >> %outfile% # 48 # ( -8) # 寒地土木研究所月報 №685 2010年6月 set of=japan_view30sec.ps # ജࡈࠔࠗ࡞ฬ set gg=grad30.grd # 㒶ᓇ࿑ࠣ࠶࠼ࡈࠔࠗ࡞ฬߩ⸳ቯ set cpo=GMT_globe.cpt set set set # ࠞࡄ࠶࠻ฬ waku=125/150/25/50/-20000/5000 # ឬ↹㗔ၞ # GEBCO ߩ 30 ⑽ⷺᒻ࠺࠲ࡈࠔࠗ࡞ߩෳᾖ grddata=GEBCO_08.nc bp=a5f2.5g0.0 grdgradient %grddata% -G%gg% -A45 -Ne0.4 # ゲߩ⸳ቯ # 㒶ᓇࠣ࠶࠼ᚑޕ-A45 ߪർ᧲ᣇะ߆ࠄߩశޔ-Ne ߪᓇߩᒝߐޕ grdview %grddata% -R%waku% -JZ4 -Y5 -JX15c -C%cpo% -I%gg% -E170/40 -B5:"Longitude":/5:"Latitude":/10000:"Elevation":SEZ -N-20000/lightgray -Qi300 -V -P -K > %of% psscale -D7/-1/10/0.3h -C%cpo% -Ba5000f2500g500 -P -O >> %of% # 㠽⍑࿑ߩᚑ # ࠬࠤ࡞ߩᚑ ࿑ ࿑ ࠍឬ↹ߔࠆߚߩࠬࠢࡊ࠻㧔ࠣ✂ដߌㇱಽߪ⺑↪ߢࠅޔታ㓙ߩࠬࠢࡊ࠻ߦߪਇⷐ㧕 pslegend ߪ GMT ߩᦨㄭߩࡃ࡚ࠫࡦߦട߃ࠄࠇߚᯏ⢻ߢޔಠࠍᝌߢ߈߹ߔޕpslegend ߢฯ߈ߔ࠺࠲ߪ GMT ߩࡃ࠶࠴ࡈࠔࠗ࡞ߢ࠻ࡊࠢࠬޔਛߢߘߩࡃ࠶࠴ࡈࠔࠗ࡞ࠍ call ᢥߢ⺒ߺㄟࠎߢታⴕߐߖߡಠࠍജߒߡ ߹ߔޕಠ↪ߩ࠺࠲ߩౝኈࠍ࿑ ߦ␜ߒߡ߅ࠅޔߢߎߎޔN 1 ߩ㧝ߪಠߩᢙࠍ␜ߒ߹ߔޕS ߪࠪࡦࡏ࡞ ࠍ⸥ㅀߔࠆߣ߁ት⸒ߢޔ0.5c ߪᏀߩಠᨒ߆ࠄ 0.5cm 㔌ࠇߚ⟎߇ࠪࡦࡏ࡞ߩਛᔃࠍ␜ߒ߹ߔޕc ߪࠍࡊ ࡠ࠶࠻ߔࠆߢߩ࡞ࡏࡦࠪޔᄢ߈ߐߪ 0.10cm⿒ޔ㧔255/0/0㧕ߢႣࠅߟ߱ߒߡ߹ߔޕಠᨒᏀ߆ࠄ 0.9cm ߩ ⟎ߦޔᢥሼ mag1 ࠍࡊࡠ࠶࠻ߒߡ߹ߔޕ࿑ Ꮐߩಠߢߪࡑࠣ࠾࠴ࡘ࠼ߦᔕߓߡߩࠨࠗ࠭߇ᄌൻߒߡ ߹ߔޕታ㓙ߩ㔡ᄩߩᄢ߈ߐ߽ಠߦḰߓߡ߹ߔޔߚ߹ޕ㔡ᄩߩ㔡ᐲ߇ᵻᣇ߆ࠄ⿒ޔ✛ޔ㕍ߦ⚡ޔᄌൻߒߡ ߅ࠅޔᶏḴઃㄭ߆ࠄᄢ㒽ߦะߌޔ㔡ᄩߩ⊒↢⟎߇ᷓߊ ߥࠆะ߇ߡขࠇ߹ߔޕ 㧔㧠㧕GEBCO ࠺࠲ࠍ↪ߚ㠽⍑࿑㧔ㅘⷞ࿑㧕 GEBCO(The General Bathymetric Chart of the Oceans)ߪᄢ ᵗ᳓ᷓ✚࿑ߣ߫ࠇޔޔWEB ߢ⸃ᐲ߇ 1 ಽⷺ(⚂ 1km)ߣ 30 ⑽ⷺ(⚂ 500㨙)ߩᒻ࠺࠲ࠍᚻߔࠆߎߣ߇ߢ ߈߹ߔޕ࿑ ߩࠬࠢࡊ࠻ࠍ↪ߒߡᣣᧄㄝߩ㒽ၞߣ ᳓ၞߩᒻࠍ㠽⍑࿑ߢ࿑␜ߒߚߩ߇࿑ ߢߔ࠻ࡊࠢࠬޕ ߢߪ㋦⋥ᣇะࠍ 4 ߦᒝ⺞ߒߡ߹ߔ߇ޔᣣᧄㄝߩᶏḴ ߿ἫጊፉߩಽᏓߥߤޔᒻߩ᭽ሶ߇⦟ߊಽ߆ࠅ߹ߔޕ 㧔㧡㧕᳓᷷ߩ᳓ᷓᣇะಽᏓߩ⚻ᐕᄌൻ࿑ ߎࠇ߹ߢߪޔ࿑ߩឬ↹߇ࡔࠗࡦߢߒߚ߇ᦨޔᓟߦ᳓ᷓ ᣇะߩ᳓᷹᷷ⷰ⚿ᨐߩᤨ㑆⊛ߥᄌേะࠍឬ↹ߔࠆᣇᴺߦ ߟߡ⺑ߒ߹ߔ⾂ࡓ࠳ޕ᳓ᳰ߿ᶏၞߦ⸳⟎ߒߚࡦࠨߢ ᳓ᷓᣇะߩ᳓⾰߫߃ޔ᳓᷷߿ṁሽ㉄⚛㊂ߥߤߩ᷹ⷰࠍⴕ ߁ߎߣ߇ࠅ߹ߔޕ᭽ⷐߥޘ࿃ߢ᳓ⅣႺ߇ᄌൻߔࠆ႐ว߇ ࿑ ࠳ࡓ⾂᳓ᳰߩ᳓᷷ಽᏓߩ⚻ᐕᄌൻ ࠅ߹ߔ߇ࠍ࠻ࡊࠢࠬࠆߔ⚫ߢߎߎޔ↪ߔࠆߣߎޔ # ᐕᣣ ᳓ 1980-05-09 180 1980-05-09 165.14 1980-05-30 172.47 1980-06-27 165.71 2001-09-04 162.17 2001-09-14 161.65 2001-09-18 161.35 2001-09-18 180 # ᐕᣣ ḓᐩ 1980-05-09 140 1980-05-09 147.54 1980-05-30 147.67 1980-06-27 148.31 2001-09-04 148.87 2001-09-14 152.55 2001-09-18 149.25 2001-09-18 140 # ᐕᣣ 1980-05-09 1980-05-09 1980-05-09 1980-05-09 1980-05-09 1980-05-09 1980-05-09 1980-05-09 1980-05-09 1980-05-09 ᤨೞ ᐕᣣ ᳓ ᳓ᷓ ᳓᷷ 00:00 29350 165.14 165.04 7.6 00:00 29350 165.14 163.14 7.2 00:00 29350 165.14 161.14 6.9 00:00 29350 165.14 159.14 6.2 00:00 29350 165.14 157.14 6 00:00 29350 165.14 155.14 5.8 00:00 29350 165.14 153.14 5.5 00:00 29350 165.14 151.14 5.4 00:00 29350 165.14 149.14 5.3 00:00 29350 165.14 147.14 5.1 ࿑ ࠺࠲ߩ ࿑ ߩࠬࠢࡊ࠻ਛߢෳᾖߔࠆࡈࠔࠗ࡞ЈᏀ߆ࠄ UWTARQNVRTPޔDVOARNQVRTPޔK\CTKYV6:6 寒地土木研究所月報 №685 2010年6月 49 gmtset INPUT_DATE_FORMAT yyyy-mm-dd # gmtset PLOT_DATE_FORMAT yy-o # 99-July gmtset ANNOT_FONT_SIZE_PRIMARY +10p # set twaku=1999-01-01T/2001-09-18T/145/165 # / / set size=15cT/10c # 15cm 10cm set tmpdat=..¥topodata¥Izari-wt.TXT # set tmpxyz=tmp.dat # set surpol=..¥topodata¥sur_polt.prn # set btmpol=..¥topodata¥btm_polt.prn # set outfile=iz_wt_pm.ps # set xaxis=pa6Of1o:"Date": #O set yaxis=a5f1:"Elevation(m)": #a5f1 set title="Water temperature profile":WSne # set bp=%xaxis%/%yaxis%:.%title% #xaxis,yaxis,title set cp=izwt.cpt # makecpt -CGMT_seis.cpt -T-2/25/1 -I > %cp% # gawk "{print $1 , $5 , $6 }" %tmpdat% > %tmpxyz% ex)1980-01-25 X yy -2 # 5m o 2 / T PlotDateFormat Y 2 99-July 6 1 bp 25 1 1 5 GMT_seis.cpt %tmpxyz% 6 # pscontour %tmpxyz% -R%twaku% -JX%size% -C%cp% -B%bp% -W0 -A2 -I -P -O -K >> %outfile% psbasemap -JX%size% -R%twaku% -B%bp% -Y5 -P -K > %outfile% # gawk "{print $1 , $2 }" %surpol% | psxy -JX -R -W1 -G190 -P -O -K >> %outfile% # gawk "{print $1 , $2 }" %btmpol% | psxy -JX -R -W1 -G190 -P -O -K >> %outfile% psscale -D7/-2/12/0.3h -Bsa5f5g5 -C%cp% -P -V -O >> %outfile% # # ߁ߒߚ᳓ⅣႺߩᤨ㑆⊛ߥᄌㆫ߽ኈᤃߦนⷞൻߢ߈߹ߔޕ࿑ ߪ⾂ࡓ࠳ࠆߣޔ᳓ᳰߩ᳓ᷓᣇะߩ᳓᷷ᄌൻࠍ⚻ ᐕ⊛ߦㅊߞߚ߽ߩߢߔߩߎޕ࿑ࠍᚑߔࠆߚߦߪᣣᤨޔ᳓ᷓޔ᳓᷷ਗ߮ߩ࠺࠲(࿑)ࠍḰߔࠆᔅⷐ߇ࠅޔ ࠬࠢࡊ࠻(࿑)ߢ߽ᤨ㑆ゲ࠺࠲ࠍ↪ࠆߚߩ․ᱶߥ߇ࠊࠇߡ߹ߔ⚦ߩ࠻ࡊࠢࠬߪߢߎߎޕ ߪ߃ߡ⺑ߒ߹ߖࠎ߇ߦ․ߪߢߎߎޔ࿑ ߩࠦࡦ࠲࿑ߩਅࠍࡑࠬࠢߔࠆ࠺࠲࠶࠻ࠍ⸳ቯߒߡ߹ߔޕ GMT ߩࡑ࠾ࡘࠕ࡞ߢߪ psmask ߣ߁ࠦࡑࡦ࠼߇ห᭽ߩᯏ⢻ࠍᨐߚߔߎߣߦߥߞߡࠆߩߢߔ߇ߊ߹߁ޔᯏ⢻ߒ ߥ߆ߞߚߚ✕ޔᕆㆱ㔍⊛ߦߎߩࠃ߁ߥࠬࠢࡊ࠻ߢኻಣߒߡ߹ߔ߁߽ߣࠆߔ߆ߒ߽ޕዋߒࠛࠟࡦ࠻ߥᣇᴺ ߇ࠆߩ߆߽ߒࠇ߹ߖࠎޕᓟߩᬌ⸛⺖㗴ߢߔߪߢߎߎޔ߅ߥޕ᳓᷷ߩᤨ㑆ಽᏓࠍᛒ߹ߒߚ߇ޔᮡ㜞୯ࠍ↪ ߒߚ╬㜞✢ಽᏓߥߤ߽ห᭽ߦឬߊߎߣ߇น⢻ߢߔޕ ߅ࠊࠅߦ ㄭᐕޔGoogle Map ߿ Google Earth ߥߤ WEB ߢ࿑ⴡᤊ↹ᖱႎ߇㐿ߐࠇޔኈᤃߦ࿑↹ᖱႎࠍ ᚻߢ߈ࠆଢߥᤨઍߦߥࠅ߹ߒߚ↹ߚߒ߁ߎޔࠄ߇ߥߒ߆ߒޕߦߪ↪㒢߇߆߆ߞߡࠆ႐ว߇ᄙߊ⺰ޔᢥ ߥߤ߳ߩᷝઃߦߪᵈᗧ߇ᔅ㗇ߢߔ ޕGMT ߪ࿑ߦ․ൻߒߚࡈߩ࠺࠲ࡑ࠶ࡇࡦࠣᚻᴺߢޔ⇇ߢ 15,000ੱએߩ↪⠪߇ࠆߎߣ߆ࠄ߽ޔ⸥ߩࠃ߁ߥ㒢ߪዋߥ߽ߩߣ߃߹ߔޔࠄ߇ߥߒ߆ߒޕGMT ߩ ↪ߦᒰߚߞߡߪߎࠇ߹ߢ UNIX ࡌࠬߩ⸃⺑߇ᄙߊޔWINDOWS ߩ߽ߩ߇ዋߥߩ߇⁁ߢߒߚޔߚ߹ޕGMT ߩࠦࡑࡦ࠼એᄖߩࠦࡑࡦ࠼ࠪࠚ࡞߿ gawk ߣߞߚࠦࡑࡦ࠼ߩ⍮⼂߽ᔅⷐߣߥࠆߚޔWINDOWS ࡙ࠩߦߣ ߞߡߪ৻ጀขߞઃ߈ߠࠄ߽ߩߦߥߞߡߚߣᕁ߹ߔ ߪߦ࠻ࡊࠢࠬߚߒ⚫ߢߎߎޕGMT ߩࠦࡑࡦ࠼એᄖ ߩ⺑߽⧯ᐓߢߔ߇ട߃߹ߒߚߩߢ߃ߐ࠲࠺ޔឥ߃߫ߔߋߦ↪ߢ߈ࠆࠃ߁ߦߒߡࠅ߹ߔ⥝ޕߩࠆᣇߪ ߎࠇࠍᯏߦ GMT ࠍߞߡߺߡߪ߆߇ߢߒࠂ߁ޕㄭᐕߩ᳇ᄌേߛߌߢߥߊޔⅣႺ߽ᄌൻߒޔౝኈ߇ ߹ߔ߹ߔᄙ᭽ൻߒߡ߹ߔޕᛛⴚ⠪⥄り߇࠺࠲ࠍนⷞൻߒߡޔ႐ߢ߇߈ߡࠆߩ߆ࠍℂ⸃ߒޔㄦㅦߦኻ ╷ࠍ⻠ߓࠆᔅⷐᕈ߽㜞߹ߞߡࠆߩߢߪߥߢߒࠂ߁߆ޕGMT ߪਁ⢻࠷࡞ߢߪࠅ߹ߖࠎ߇ޔ᭽ߥޘ႐㕙ߢᬺ ോߦᓎ┙ߟ࠷࡞ߢࠆߎߣߪ⏕߆ߢߔ╩ޔߛߚޕ⠪߽߹ߛߘࠇ߶ߤᷓߊ GMT ࠍℂ⸃ߒߡࠆ⸶ߢߪߥߩߢޔ ߎߎߢ⚫ߒߚࠬࠢࡊ࠻߇ᔅߕߒ߽ᱜߒߣߪ㒢ࠅ߹ߖࠎߣߞ߽ޕല₸⊛ߥࠬࠢࡊ࠻ߩࠨࡦࡊ࡞߇ࠅ߹ߒ ߚࠄߏᢎ␜㗂ߌࠇ߫ሜߒߢߔ⺑ߪߢ✬ᧄޔߚ߹ޕߒ߈ࠇߥ߆ߞߚ߽ᄙߩߢࠄߚߒ߹ࠅ╬⇼⾰ߒ߽ޔ㆙ ᘦߥߊࠦࡦ࠲ࠢ࠻㗂ߌࠇ߫ᐘ↟ߢߔޕ 50 寒地土木研究所月報 №685 2010年6月 ෳ⠨ᢥ₂ 㧝㧕9GUUGN2CPF9*(5OKVJ(TGGUQHVYCTGJGNRUOCRCPFFKURNC[FCVC'QU6TCPU#)7 㧞㧕)'$%1 6JG)GPGTCN$CVJ[OGVTKE%JCTVQHVJG1EGCPUJVVRYYYIGDEQPGV 村上 泰啓* Yasuhiro Murakami ᵏ㧖 寒地土木研究所 寒地水圏研究グループ ኙᧁ⎇ⓥᚲ 寒地河川チーム・景観ユニット (兼務) 主任研究員 ኙ᳓⎇ⓥࠣ࡞ࡊ 工学博士 Yasuhiro Murakami ኙᴡᎹ࠴ࡓ᥊࡙ⷰ࠾࠶࠻㧔ോ㧕 ਥછ⎇ⓥຬ Ꮏቇඳ჻ )/6 ߦ㑐ߔࠆ 9'$ ࠨࠗ࠻৻ⷩ 9'$ ࡍࠫ࠲ࠗ࠻࡞ ⠪ )/6 ߩᣇ ᚭᄢቇℂቇㇱ ╽ᭉ㤚ߐࠎ )/6 ࠃࠃࠟࠗ࠼㧔╙ 㧦)/6XGT ↪㧕 )/6 ߚ߹ߏࠢࡉ ⁛↥ᬺᛛⴚ✚ว⎇ⓥᚲ ᵴᢿጀ㔡 ⎇ⓥࡦ࠲ ၳᎹ᥍ᄩߐࠎ )/6 ߩㇱደ ᶏᵗ⎇ⓥ㐿⊒ᯏ᭴ ⅣႺࡈࡠࡦ࠹ ࠖࠕ⎇ⓥࡦ࠲ ᒻᐽߐࠎ )/6 ᔓ㍳ -CUJKOCߐࠎ *KIJ SWCNKV[ OCR ࠞࡈࠜ࡞࠾ࠕᄢ࠺ࡆࠬᩞფ⾗Ḯ RTQFWEVKQP ⎇ⓥᚲ 6JG)GPGTKE/CRRKPI6QQNU ࡂࡢࠗᄢቇ 51'56 ᩞ㧔)/6 ߩᧄኅ㧕 74. JVVRYYYUGKURNCPGVUEKMQDGWCELR`MCMGJ K)/6)/6*1961JVON JVVRUVCHHCKUVIQLRJJQTKMCYC)/6)/6I WKFG[QK[QKJVON JVVRYYYICKCJM[QVQWCELR`OMCVQ)/6)/6 VCOCIQJVON JVVRJ[FTQKKUWVQM[QCELR`CICVCCTEJKXG )/6 JVVRYYYRNCNCQTLRMCUJKOCIOVIOVJVON JVVRECUQKNTGUQWTEGNCYTWEFCXKUGFWFTWRCN DQQMGZRQTVJVON JVVRIOVUQGUVJCYCKKGFWIOVIOVAYKPFQYUJV ON )/6 8GTUKQP ࡂࡢࠗᄢቇ 2CWN9GUUGN9CNVGT*( JVVRYYYUQGUVJCYCKKGFW)/6IOVFQEIOVJV ON)/6A&QEU)/6A&QEUJVON 6GEJPKECN 4GHGTGPEG CPF 5OKVJ %QQMDQQM 寒地土木研究所月報 №685 2010年6月 51