clear
clc
t=(1/16:1/8:1)'*2*pi;%八边形顶点
x1=sin(t);
x1=[x1;x1(1)];
y1=cos(t);
y1=[y1;y1(1)];
plot(x1,y1)
%漳州台八边形经纬度
zz_octogan=[ x1+117.63 y1+24.43];
==============gmt的psxy -Sg可以直接画图======
set ps=zzoctagon.ps
set J=M6i
set R=115/121.5/22/29
gmt gmtset FONT_TITLE 40p,36,black
gmt gmtset FORMAT_GEO_MAP ddd:mm
gmt gmtset MAP_FRAME_TYPE fancy
gmt gmtset MAP_FRAME_WIDTH 0.1c
gmt gmtset GRID_CROSS_SIZE_PRIMARY = 0c
gmt gmtset TICK_LENGTH -0.0c
gmt gmtset ANNOT_FONT_PRIMARY = Helvetica
rem 写入文件头
gmt psxy -J%J% -R%R% -Ba2f2WSne -T -K > %ps%
rem gmt grdimage china.grd -Itemp.grd -J -R -Ba2f2WSne -Cwiki-scotland.cpt -K -O>>%ps%
rem gmt grdcut topo30.grd -Gtemp.grd -R%R%
rem gmt grdgradient temp.grd -A50 -Gtemp.grd.int -Ne0.6
rem 画地形
rem gmt grdimage temp.grd -J%J% -R%R% -Ba2f2WSne -Cetopo1.cpt -Itemp.grd.int -V -P -K -O >> %ps%
gmt pscoast -J -R -B -Gwhite -Slightblue -K -O >>%ps%
gmt psxy shengjie.xy -R -J -Ba2f2WSne -K -O -W0.8/255/0/0 >> %ps%
gmt psxy fujian_xi.xy -R -J -K -O -W0.7/255/0/0 >>%ps%
gmt psxy deng_fault.xy -R -J -K -O -W0.8/255/0/0 >>%ps%
rem gmt psxy ..Data2zl.txt -R -J -K -O -Sg -Gblue >> %ps%
rem echo 118.5 26|psxy -R -JX -Sg0.4 -W0.1p/0 -G0/0/255 -K -O>> %ps%
rem gmt psxy zl.txt -R -J -Sc -W0.5,black, -Gyellow -K -O >> %ps%
gawk "{print $1,$2}" zl.txt |psxy -R -J -Sc0.4c -Gred -K -O >> %ps%
rem gmt psxy octagon.txt -R -J -Sc0.4c -Gred -K -O >> %ps%
gmt psxy octagon.txt -R -J -K -O -W0.8/255/0/0 >>%ps%
REM gawk "{print $5,$4,0.12*$3 }" hn6.txt |psxy -R -J -Sc -W0.5,black, -Gred -K -O >> %ps%
rem gmt psxy zl.txt -J%J% -R%R% -B%B% -Sc0.2c -Gred -P -K -O >> %ps%
REM gmt gmtset FONT_ANNOT_PRIMARY 10p,STSong-Light--UniGB-UTF8-H,black
REM gmt pslegend legend_hn.txt -J -R -Dx0.1i/-0.65i/4i/0.4i/BL -F+gwhite -O -K -L1.5 -V >> %ps%
gmt psxy -J%J% -R%R% -T -O >> %ps%
gmt psconvert %ps% -Tj -E400 -A -P -C-sFONTPATH=c:windowsFonts
del gmt.conf gmt.history