[更新情報]

▽Go to footer

GMT で報告書用の説明図を作成する事例です. これまで説明図関係は guplot で作成していましたが,GMTを用いて説明図を作成する事例を作成してみました. パターンでの塗りつぶしと矢印の使い方がポイントです.

実行用バッチファイルと解説

rem *******************************************
rem * drawing example
rem *******************************************
del .*
rem ===========================================
rem * データ作成
rem ===========================================
rem 塗りつぶしを行う矩形領域を指定(_temph0.txt)(閉じてなくてよい)
rem ***Sand filling shape
echo -4.75 0.0 > _temph0.txt
echo  4.75 0.0 >> _temph0.txt
echo  4.75 1.7 >> _temph0.txt
echo -4.75 1.7 >> _temph0.txt
rem 塗りつぶしを行う矩形領域を指定(_temph1.txt)(閉じてなくてよい)
rem ***CSEM shape
echo -4.75 1.7 > _temph1.txt
echo  4.75 1.7 >> _temph1.txt
echo  4.75 5.2 >> _temph1.txt
echo -4.75 5.2 >> _temph1.txt
rem 塗りつぶしを行う円形領域を指定(_temph2.txt)(閉じてなくてよい)
rem ***Pipe A area
gawk "BEGIN{pi=3.141592654;for(i=0;i<=365;i++){print 1.5*cos(i/180*pi)-2.5,1.5*sin(i/180*pi)+1.7}}" > _temph2.txt
rem 塗りつぶしを行う円形領域を指定(_temph3.txt)(閉じてなくてよい)
rem ***Pipe B area
gawk "BEGIN{pi=3.141592654;for(i=0;i<=365;i++){print 1.5*cos(i/180*pi)+2.5,1.5*sin(i/180*pi)+1.7}}" > _temph3.txt
rem
gmtset FRAME_PEN 0
gmtset VECTOR_SHAPE 1
set fig_out=fig_gmt_pipe.eps
rem **************************************************************************
rem gmtset FRAME_PEN 0           psbasemapのframeのpen幅を0に指定
rem gmtset VECTOR_SHAPE 1        Svで描画するベクトルの矢印の形を1に指定
rem set fig_out=fig_gmt_pipe.eps 出力画像ファイル名を指定
rem **************************************************************************
rem ===========================================
rem * プロット実行
rem ===========================================
psbasemap -R-5/10/-2/6 -JX15/8 -B -P -K > %fig_out%
psxy _temph0.txt -R -J -Gp150/31 -W3 -O -K >> %fig_out%
psxy _temph1.txt -R -J -Gp150/8 -W3 -O -K >> %fig_out%
psxy _temph2.txt -R -J -G255 -W5 -O -K >> %fig_out%
psxy _temph3.txt -R -J -G255 -W5 -O -K >> %fig_out%
rem **************************************************************************
rem 2つの矩形領域の塗りつぶしと黒枠線を描画後、
rem 円形領域の塗りつぶしと黒枠線を描画
rem 上書きで塗りつぶされるので,塗りつぶす順番に注意する
rem -Gp150/31 150dpi  ハッチングパターン31で塗りつぶし
rem -W3               黒枠線を太さ3で描画
rem **************************************************************************
rem
echo -2.50  1.7 -2.50 -0.6 | psxy -R -J -SvS0.005/0.0/0.0 -G0 -O -K >> %fig_out%
echo  2.50  1.7  2.50 -0.6 | psxy -R -J -SvS0.005/0.0/0.0 -G0 -O -K >> %fig_out%
echo -2.50 -0.5  2.50 -0.5 | psxy -R -J -SvS0.005/0.2/0.1 -G0 -O -K >> %fig_out%
echo -4.75  0.0 -4.75 -1.1 | psxy -R -J -SvS0.005/0.0/0.0 -G0 -O -K >> %fig_out%
echo  4.75  0.0  4.75 -1.1 | psxy -R -J -SvS0.005/0.0/0.0 -G0 -O -K >> %fig_out%
echo -4.75 -1.0  4.75 -1.0 | psxy -R -J -SvS0.005/0.2/0.1 -G0 -O -K >> %fig_out%
echo  0.00 -0.7 10 0 0 CM 5000| pstext -R -J -N -O -K >> %fig_out%
echo  0.00 -1.2 10 0 0 CM LIMITS OF TRENCH EXCAVATION | pstext -R -J -N -O -K>> %fig_out%
rem **************************************************************************
rem -SvS ベクトルは始点と終点で指定,両矢印とする
rem 矢印形状は「arrowwidth/headlength/headwidth」で指定
rem 0.005/0.0/0.0 線の太さ0.005cmの線(矢印無し)
rem 0.005/0.2/0.1 線の太さ0.005cm,矢印頭の長さ0.2cm,矢印頭の幅0.1cm
rem pstextのデータは「x y size angle fontno justify text」で指定
rem **************************************************************************
rem
echo -3.56066 0.63934 -1.43934 2.76066 | psxy -R -J -SvS0.005/0.2/0.1 -G0 -O -K >> %fig_out%
echo  -2.64 1.84 8 45 0 CM 3000 FID MS PIPE| pstext -R -J -N -O -K >> %fig_out%
echo  1.43934 0.63934  3.56066 2.76066 | psxy -R -J -SvS0.005/0.2/0.1 -G0 -O -K >> %fig_out%
echo   2.36 1.84 8 45 0 CM 3000 FID MS PIPE| pstext -R -J -N -O -K >> %fig_out%
rem
echo  4.75  0.0  5.00 0.0 | psxy -R -J -SvS0.005/0.0/0.0 -G0 -O -K >> %fig_out%
echo  2.50  0.2  5.00 0.2 | psxy -R -J -SvS0.005/0.0/0.0 -G0 -O -K >> %fig_out%
echo  5.10  0.1 10 0 0 LM 200 MIN. BEDDING| pstext -R -J -N -O -K>> %fig_out%
rem
echo  4.75  5.2  5.30 5.2 | psxy -R -J -SvS0.005/0.0/0.0 -G0 -O -K >> %fig_out%
echo  2.50  3.2  5.30 3.2 | psxy -R -J -SvS0.005/0.0/0.0 -G0 -O -K >> %fig_out%
echo  5.20  3.2  5.20 5.2 | psxy -R -J -SvS0.005/0.2/0.1 -G0 -O -K >> %fig_out%
echo  5.30  4.4 10 0 0 LM Min. COVER| pstext -R -J -N -O -K>> %fig_out%
echo  5.30  4.0 10 0 0 LM 1000 @%%12%%\176@%%%% 2000| pstext -R -J -N -O -K>> %fig_out%
rem
echo   0.0 0.8 10 0 0 CM COMPACTED| pstext -R -J -N -O -K >> %fig_out%
echo   0.0 0.4 10 0 0 CM SAND FILL| pstext -R -J -N -O -K >> %fig_out%
echo   0.0 4.5 10 0 0 CM COMPACTED SUITABLE| pstext -R -J -N -O -K >> %fig_out%
echo   0.0 4.1 10 0 0 CM EXCAVATED MATERIAL| pstext -R -J -N -O -K >> %fig_out%
echo   0.0 5.3 10 0 0 CB FGL| pstext -R -J -N -O >> %fig_out%
rem
del _*
del .*

バッチファイルのダウンロードと画像確認

ファイル名概要
作図用バッチファイル
出力画像



inserted by FC2 system