pstext - 360文档中心

GMT使用手册

1
GMT 是用具有高可移植度的 ANSI C 语言编写，与 POSIX 标准相适应。它经过少许改动就可以在所有运行 UNIX 软件的硬件上运行。在 GMT 的编写中，我们遵循了标准的 UNIX 的编写原则：把原始数据--分析处理--图像显示的流程分解为一系列基本的步骤，每一个步骤都可以用一个 UNIX 或 GMT 工具来完成。这个方法带来了如下的好处：（1）只需要少量的程序，（2）每一个程序简洁，并且易于更新，（3）每一个步骤都是独立的，因此可以有多种应用（每一步骤都可以任意组合），（4）程序可以用 shell 或其他的管道（pipes）连接到一起，这样就可以根据需要建立自己想要的数据处理流程。因为每一个科研机构的数据库结构不同，这样的话，数据获取和后续的处理，作图之间的解耦就很重要。如果我们想用自定义（custom）的数据库，我们就可以自己编一个程序把这种数据库的格式转化为 GMT 可读的格式。接下来就可以用标准的 GMT 函数了。
GMT 中文攻略 β 版
北京大学地球与空间科学学院地球物理系佟啸鹏
目录
前言 ..................................................................................................................................................1 1 简介 ..............................................................................................................................................1 2 GMT 概览和快速索引.................................................................................................................3 3 GMT 的基本特征..........................................................................................................................5

pst码编码规则 -回复

pst码编码规则-回复pst码是由世界邮政联盟（UPU）制定的国际邮政编码系统。

它的全称是Postal Service Tracking (PST) Code，俗称“邮政编码”。

pst码的编码规则是根据地理区域划分，并且在每个地理区域内按照一定的编号规则进行编码。

下面将一步一步回答关于pst码编码规则的问题。

首先，pst码编码规则中的第一步是确定地理区域。

通常情况下，地理区域是按照国家或地区进行划分的。

不同的国家或地区都会有自己的pst码编码规则。

第二步是确定地理区域内的具体行政区划。

在一个国家或地区中，通常有不同的行政区划，例如省、州、县等。

每个行政区划都会被分配一个特定的pst码。

第三步是确定pst码的具体编码规则。

在每个行政区划内，pst码的编码规则通常是按照数字顺序进行分配的。

具体的编码规则可能随各个国家或地区的需要而有所不同。

以中国为例，pst码编码规则如下：1. 国家码：中国的国家码是86，代表中国这个国家。

2. 省级行政区划码：中国的省级行政区划码是根据汉字的拼音首字母进行编码的。

例如，广东省的行政区划码是GD，山东省的行政区划码是SD 等。

3. 地级行政单位码：在每个省级行政区划内，还会有各个地级行政单位，例如市、自治州、地区等。

地级行政单位码的编码规则是根据字母的顺序进行分配的。

例如，广州市的行政单位码是A，深圳市的行政单位码是B 等。

4. 县级行政单位码：在每个地级行政单位内，还会有各个县级行政单位，例如县、区、镇等。

县级行政单位码的编码规则通常是按照数字顺序进行分配的。

例如，广州市内的海珠区的编码是01，番禺区的编码是02等。

通过以上的编码规则，可以得到一个完整的pst码。

例如，广州市海珠区的pst码就是8601A01，其中86代表中国，01代表广东省，A代表广州市，01代表海珠区。

pst码的作用是方便邮政服务的快捷和准确。

通过pst码，邮件和包裹可以更快地送达目的地，并且可以准确地确定其所在地区。

win32 api-第四章

caseWM_PAINt (hwnd, &ps) ;
使用GDI函数
EndPaint (hwnd, &ps) ;
return 0 ;
在处理WM_PAINT消息时，必须成对地呼叫BeginPaint和EndPaint。如果窗口消息处理
这个区域称为「无效区域」或「更新区域」。正是显示区域内无效区域的存在，才会让Windows
将一个WM_PAINT消息放在应用程序的消息队列中。只有在显示区域的某一部分失效时，
窗口才会接受WM_PAINT消息。
Windows内部为每个窗口保存一个「绘图信息结构」，这个结构包含了包围无效区域的最小
可以保留绘制显示区域需要的所有信息，并且仅当「响应要求」－即Windows给窗口消息
处理程序发送WM_PAINT消息时才进行绘制。如果程序在其它时间需要更新其显示区域，
它可以强制Windows产生一个WM_PAINT消息。这看来似乎是在屏幕上显示内容的一种舍
近求远的方法。但您的程序结构可以从中受益。
的显示区域的方式，以及显示文字时Windows使用的字体。
当程序需要绘图时，它必须先取得设备内容句柄。在取得了该句柄后，Windows用内定的
属性值填入内部设备内容结构。在后面的章节中您会看到，可以通过呼叫不同的GDI函数改
变这些默认值。利用其它的GDI函数可以取得这些属性的目前值。当然，还有其它的GDI
程序不处理WM_PAINT消息，则它必须将WM_PAINT消息传递给Windows中
DefWindowProc（内定窗口消息处理程序）。DefWindowProc以下列代码处理WM_PAINT
消息：
case WM_PAINT:

windows编程输出文字

输出文字壹佰软件开发小组整理编输出文字译在前一章，您看到了一个简单的Windows 98程序，它在窗口中央，或者更准确地说，在显示区域中央显示一行文字。

正如我们学到的，显示区域是整个应用程序窗口中未被标题列、窗口边框，以及可选的菜单列、工具列、状态列和滚动条占据的部分。

简而言之，显示区域是窗口中可以由程序任意书写和传递视觉信息的部分。

对于程序的显示区域，您几乎可以为所欲为，只不过您不能假定窗口大小是某一特定尺寸，或者在程序执行时其大小会保持不变。

如果您不熟悉图形窗口环境的程序设计，这些限制可能会使您感到惊讶：不能再假设屏幕上的一行文字一定有80个字符了。

您的程序必须与其它Windows程序共享视讯显示器。

Windows使用者控制程序窗口在屏幕上显示的方式。

尽管可以建立固定大小的窗口（这对于计算器之类的应用是合理的），但在大多数情况下，使用者应该能够改变应用程序窗口的大小。

您的程序必须能够接受指定给它的大小，并且合理地利用这一空间。

这有两种可能的情况。

一种可能是，程序只有仅能显示「hello」的显示区域；还有另一种可能，即程序在一个大屏幕、高分辨率的系统上执行，其显示区域大得足以显示两整页文字。

灵活地处理这两种极端是Windows程序设计的要点之一。

这一章，我们将讲述程序在显示区域显示信息的方式，但比上一章说明的显示方式更加复杂。

当程序在显示区域显示文字或图形时，它经常要「绘制」它的显示区域。

本章着重讲述绘制的方法。

尽管Windows为显示图形提供了强大的图形设备接口（GDI）函数，但在这一章中，我只介绍简单文字行的显示。

我也将忽略Windows能够使用的不同字体外形及字体大小，仅使用Windows的内定系统字体。

这看起来似乎是一种限制，其实不然，本章涉及和解决的问题适用于所有Windows程序设计。

在混合显示文字和图形时，Windows内定字体的字符大小通常决定了图形的尺寸。

本章表面上是讨论绘图的方法，实际上是讨论与设备无关的程序设计基础。

pstack用法

pstack用法
pstack是一个用于查看正在运行的进程的函数堆栈的工具。

它通常在调试或分析程序崩溃或陷入无限循环的情况下使用。

pstack可以为给定进程的每个线程生成函数调用堆栈的列表，并显示函数名、地址和其他相关信息。

使用pstack的基本语法如下：
其中，pid是要查看堆栈的进程ID。

下面介绍一些常用的pstack选项：
1. -p: 输出显示pid和线程ID（tid）
2.-l:定义一个长列表，显示堆栈中的每个函数调用和其对应的源文件和行号
3. -F: 强制pstack执行，即使进程采取了coredumps默认的呈现机制
4.-T:包括线程ID的堆栈跟踪，可以查看每个线程的函数调用堆栈
5. -E: 反汇编版本的文件（只有在安装了addr2line和objdump工具时才可用）
6. -a: 可以检查指定进程的core文件，并显示其函数调用堆栈（需要core文件）
7. -c: 指定core文件的名称，并显示与core文件关联的函数调用堆栈
8. -e: 指定一个附加到程序的动态链接库，并显示堆栈信息（需要指定-so选项）
下面是一些pstack使用示例：
1. 使用pid查看进程的函数调用堆栈：
2.使用附加的共享库查看进程的函数调用堆栈：
3. 查看core文件中的函数调用堆栈：
4.使用长列表显示函数调用堆栈：
通过分析pstack的输出，可以获得正在运行的进程中函数调用的详细信息。

这对于定位程序崩溃或陷入无限循环的问题非常有用。

04 指令解说

第四章、GMT的指令解說GMT一共有61個指令來控制圖形的繪製；但是在實際測試後，發現有些指令功能重疊，有些指令則並不常用。

故整理之後，只有9個；並將10個常用指令與其環境參數，整理分述如下。

在上一章講解過的共同環境參數，在此就不再贅言；上一章未講解到的環境參數部份，則只介紹最常用的環境參數。

實際上的應用則在下一章會提到。

指令解說：1. grdcontour–將網格檔(*.grd)繪出封閉曲線。

grdcontour grdfile - Ccont_int - Jparameters [ -Aanot_int[ffont_size][aangle][/r/g/b][o]] ] [ -Btickinfo ] [-Ddumpfile ] [ -Eazimuth/elevation ] [ -Fred/green/blue ] [- Ggap/width ] [ -K ] [ -Llow/high ] [ -M[flag] ] [ -N[unit]] [ -O ] [ -P ] [ -Qcut ] [ -Rwest/east/south/north[r] ] [ -Ssmoothfactor ] [ - T[+|-][gap/length][:LH] ] [ -U[/dx/dy/][label] ] [ -V ] [ -Wtype[pen] ] [ -Xx-shift ] [-Yy-shift ] [ -Zfactor[/shift] ] [ -ccopies ]‧grdfile：輸入的網格檔(*.grd)。

‧- Ccont_int：將網格檔套入等高線資訊；在此，我們使用色彩調色表(color palette table，以後簡稱*.cpt)來控制等高線的繪製。

cont_int為我們代入的*.cpt檔。

ex：grdcontour e117-125n18-27.grd -Cfig02.cpt -JM16 -Ba1f0.5/a1/f0.5 -R118/121/22/25 -X3 -Y3 –P2.grdimage–由網格檔來繪製我們所需要的影像。

GMT R包用户手册说明书

Package‘gmt’October13,2022Version2.0.3Date2022-06-06Title Interface Between GMT Map-Making Software and RImports utilsSystemRequirements gmtLazyData yesDescription Interface between the GMT map-making software and R,enabling the user to manipulate geographic data within R and call GMT commands to draw andannotate maps in postscript format.The gmt package is about interactive dataanalysis,rapidly visualizing subsets and summaries of geographic data,whileperforming statistical analysis in the R console.License GPL-3URL https://NeedsCompilation noAuthor Arni Magnusson[aut,cre]Maintainer Arni Magnusson<********************>Repository CRANDate/Publication2022-06-0620:10:08UTCR topics documented:gmt-package (2)deg2num (3)geodist (4)gmt (5)gmt.demo (6)gmt.system (7)num2deg (8)psbar (9)psclose (11)pscoast (12)12gmt-package pstext (13)psxy (14)r2gmt (15)Index17 gmt-package Interface between the GMT Map-Making Software and RDescriptionInterface between the GMT map-making software and R,enabling the user to manipulate geo-graphic data within R and call GMT commands to draw and annotate maps in postscript format. DetailsInitialize GMT session:gmt set graphical parameters and postscriptﬁleCreate a map:pscoast drawpsxy annotatepstext annotatepsbar annotatepscloseﬁnalizeConvert and calculate:deg2num degrees to numericgeodist distance between coordinatesnum2deg numeric to degreesExamples:demo.par,demo.coast,demo.xy,demo.text,demo.bar strings and data frames to create example map Internal:gmt.system,r2gmt helping functionsThis package provides functions to draw basic maps with GMT,along with helping functions that can be used to add more advanced features to a map.GMT users typically write shell scripts to draw maps.The gmt package is about interactive data analysis,rapidly visualizing subsets and summaries of geographic data,while performing statisticaldeg2num3 analysis in the R console.Author(s)Arni Magnusson.ReferencesWessel,P.,Smith,W.H.F.,Scharroo,R.,Luis,J.and Wobbe,F.The Generic Mapping Tools:GMT documentation.Available at https://.deg2num Convert Degrees to NumericDescriptionConvert deg:min:sec string to decimal number.Usagedeg2num(x)Argumentsx string or vector of strings.DetailsDegrees,minutes and seconds are separated by colons,and each can have a decimal point as well.First character must be a minus sign or number,and last character must be W,E,S,N,or number. ValueNumeric representation of the degree string(s).NoteThe string format is adopted from Appendix B.1.1in the GMT manual.See Alsoas.numeric converts strings to numbers when things are straightforward.deg2num is the opposite of num2deg.gmt-package gives an overview of the package.Examples#The opposite of num2deg()exampledeg2num(c("12:30:44.5W","17.5S","1:00:05","200:45E"))4geodistgeodist Distance Between Geographic CoordinatesDescriptionCalculate surface distance between geographic agegeodist(Nfrom,Efrom,Nto,Eto,units="km")ArgumentsNfrom latitude of origin.Efrom longitude of origin.Nto latitude of destination.Eto longitude of destination.units how distance is measured:"km"for kilometres,"nm"for nautical miles.DetailsLatitude and longitude are passed as decimal numbers,e.g.66.5for 66◦30 N.Vectors of coordinates are supported.ValueVector of distances.NoteThe surface distance between geographic coordinates is:D =cos −1 sin θ1sin θ2+cos θ1cos θ2cos(φ1−φ2)where distance and coordinates are expressed in radians.θ1and θ2is the latitude of origin and destination,and φ1and φ2is longitude.The calculations assume a perfect sphere and elevation differences are ignored.The SI deﬁnition of a nautical mile is exactly 1.852km.See Alsodiff ,Trig .gmt-package gives an overview of the package.Examplesgeodist(55.75,37.63,39.9,116.4)#Moscow -Beijing geodist(90,0,-90,0,"nm")#from pole to polegmt5 gmt Initialize GMT SessionDescriptionInitialize a GMT session by setting graphical parameters and current postscriptﬁle.Usagegmt(par=NULL,file="map.ps",style="s",quiet=TRUE)Argumentspar optional string of GMT parameters,"PARAMETER1=value1PARAMETER2=value2 ...",passed to GMT shell command gmtset.file postscriptﬁlename used in subsequent calls to GMT plot functions,passed to options.style default units and values:"s"for SI(international)or"u"for United States, passed to GMT shell command gmtset.quiet whether current settings should be displayed.DetailsThe file argument can be supplied with(recommended)or without a full directory path.Withouta path,the R working directory is used(see getwd and setwd).See the GMT documentation for details on graphical parameters,gmtdefaults,gmtset and other GMT commands.ValueList containing the current options("gmt.file").If par is NULL,no GMT parameters are changed,but the current parameter values and postscript ﬁlename can be reviewed.If par is a string(empty""will do),a GMT conﬁgurationﬁle is created in the current map directory, dirname(file).See Alsooptions could be used to set gmt.file directly.gmt,pscoast,psxy,pstext,psbar,and psclose work together to draw maps.gmt-package gives an overview of the package.6gmt.demoExamples##Not run:#Draw map and save as"map.ps"in R working directorygmt(demo.par)pscoast(demo.coast)psxy(demo.xy)pstext(demo.text,"-J-R-F+f+a+j-O-K")psbar(demo.bar,ref=66)psclose()#See directory gmt/example for details##End(Not run)gmt.demo GMT Example DataDescriptionTheseﬁve objects are provided to demonstrate the functionality of the gmt package.Usagedemo.pardemo.coastdemo.xydemo.textdemo.barFormatdemo.par and demo.coast are simple strings.demo.xy is a data frame containing:Lon LongitudeLat LatitudeSize Size of plot symboldemo.text is a data frame containing:Lon LongitudeLat LatitudeSize Size of plot symbolAngle Angle in degrees counter-clockwise from horizontalFont Font numberJustify Justiﬁcation codeText Text labelgmt.system7demo.bar is a data frame containing:Lon LongitudeLat LongitudeWidth Bar width in degreesHeight Bar height in degreesDetailsSee the GMT documentation for details on psxy,pstext and other GMT commands.See Alsogmt-package gives an overview of the package.Examples##Not run:#Draw map and save as"map.ps"in R working directorygmt(demo.par)pscoast(demo.coast)psxy(demo.xy)pstext(demo.text,"-J-R-F+f+a+j-O-K")psbar(demo.bar,ref=66)psclose()#See directory gmt/example for details##End(Not run)gmt.system Invoke shell commandDescriptionThis internal function invokes a shell command,possibly directing the output to aﬁle.Usagegmt.system(cmd,file=NULL,append=FALSE)Argumentscmd system command to be invoked,as a string.fileﬁlename to which output should be directed.append whether output should be appended to existingﬁle.8num2degValueCommand output as a vector of strings.Notegmt.system is a fast platform-independent wrapper for system,supporting redirection toﬁle.It is mainly called by other functions,but users mayﬁnd it useful for running various GMT com-mands.See Alsosystem,writeLines.gmt-package gives an overview of the package.Examples##Not run:#Assuming bermuda.nc is in R working directorygmt.system("gmt grdcontour bermuda.nc-JM7i-C250-A1000-B2",file="b.ps")##End(Not run)num2deg Convert Numeric to DegreesDescriptionConvert decimal number to deg:min:sec string.Usagenum2deg(x,lat=NA,dec=FALSE,digits=0,zero=FALSE)Argumentsx number or vector of numbers.lat whether x is latitude.dec whether to return decimal degrees instead of deg:min:sec.digits precision used when rounding decimal degrees or seconds.zero whether trailing:00zeros should be retained.DetailsElement-speciﬁc format is supported,using vectors for lat,digits,and precision.The resulting string ends with N or S when lat is TRUE,E or W when lat is FALSE,or a number when lat is NA.Valuedeg:min:sec string representation of the number(s).NoteThe string format is adopted from Appendix B.1.1in the GMT manual.See Alsoas.character converts plain numbers to strings.num2deg is the opposite of deg2num.gmt-package gives an overview of the package.Examples#The opposite of deg2num()examplenum2deg(c(-12.51236,-17.5,1.00139,200.75),lat=c(FALSE,TRUE,NA,FALSE),dec=c(FALSE,TRUE,FALSE,FALSE),digits=c(1,1,0,0))psbar Add Bars to GMT Mercator MapDescriptionCall GMT to add bars to a map and save in postscript format.Usagepsbar(x,cmd="-J-R-W1p-G180-O-K",file=getOption("gmt.file"), ref=0,digits=getOption("digits"))Argumentsx data frame,matrix,orﬁlename containing the data to be plotted.cmd string of arguments passed to GMT shell command psxy.fileﬁlename where the map is saved.ref reference latitude where height1renders a bar1degree high.digits precision used when rounding the geographic coordinates.DetailsThe data are arranged in four columns:Lon,Lat,Width,and Height,in that order.If x is aﬁlename,the data should be tabular with or without a header,separated by commas or whitespace.Theﬁrst line is interpreted as header if theﬁrst non-whitespace character is not minus, point,or number.This function provides an alternative to psxy-Sb and psxy-Sr for drawing bars on a Mercator map.See the GMT documentation for details on psxy and other GMT commands.ValueNull,but the map is annotated and saved in postscript format.The temporary GMT inputﬁle‘bar.gmt’is saved in directory dirname(tempdir()),for the user to view or edit.It is later removed by psclose().NoteThis function does the necessary calculations to draw bars in standard height given a Mercator-projected map.It is not intended for other projections.The derivative of the Mercator projection is used to standardize the bar height:f (θ)=12tanπ4+θ2cos2π4+θ2whereθis latitude in radians and f(θ)is the y-axis coordinate.The bar height at a given latitude is h×f (θref)/f (θ),where h is the height passed by the user andθref is a reference latitude where h=1renders a bar1degree high.See AlsoSimilar to barplot and postscript in native R graphics.gmt,pscoast,psxy,pstext,psbar,and psclose work together to draw maps.gmt-package gives an overview of the package.Examples##Not run:#Draw map and save as"map.ps"in R working directorygmt(demo.par)pscoast(demo.coast)psxy(demo.xy)pstext(demo.text,"-J-R-F+f+a+j-O-K")psbar(demo.bar,ref=66)psclose()#See directory gmt/example for details##End(Not run)psclose11 psclose Finalize GMT MapDescriptionCall GMT toﬁnalize a map and save in postscript format.Usagepsclose(file=getOption("gmt.file"),trailer=TRUE)Argumentsfileﬁlename where the map is saved.trailer whether a closing trailer should be appended to the postscriptﬁle.DetailsA closing trailer is required if the last plotting command included-K(default behaviour).ValueNULL,but the map isﬁnalized and saved in postscript format.NoteThis function performs two tasks:1.Appends a closing trailer to the postscriptﬁle(optional).2.Removes GMTﬁles in temporary directory.See AlsoAnalogous to dev.off and postscript in native R graphics.gmt,pscoast,psxy,pstext,psbar,and psclose work together to draw maps.gmt-package gives an overview of the package.Examples##Not run:#Draw map and save as"map.ps"in R working directorygmt(demo.par)pscoast(demo.coast)psxy(demo.xy)pstext(demo.text,"-J-R-F+f+a+j-O-K")psbar(demo.bar,ref=66)psclose()#See directory gmt/example for details12pscoast ##End(Not run)pscoast Draw GMT MapDescriptionCall GMT to draw a map(coastlines,borders,rivers)and save in postscript format.Usagepscoast(cmd,file=getOption("gmt.file"))Argumentscmd string of arguments passed to GMT shell command pscoast.fileﬁlename where the map will be saved.DetailsThe file argument can be supplied with(recommended)or without a full directory path.Withouta path,the current working directory is used(see getwd and setwd).See the GMT documentation for details on pscoast and other GMT commands.ValueNULL,but a map is drawn and saved in postscript format.See AlsoSimilar to plot and postscript in native R graphics.gmt,pscoast,psxy,pstext,psbar,and psclose work together to draw maps.gmt-package gives an overview of the package.Examples##Not run:#Draw map and save as"map.ps"in current working directorygmt(demo.par)pscoast(demo.coast)psxy(demo.xy)pstext(demo.text,"-J-R-F+f+a+j-O-K")psbar(demo.bar,ref=66)psclose()#See directory gmt/example for details#Map in one callpstext13 pscoast("-JM12c-R7E/38E/29N/48N-G100-B5","x.ps")##End(Not run)pstext Add Text/Symbols to GMT MapDescriptionCall GMT to add text/symbols to a map and save in postscript format.Usagepstext(x,cmd="-J-R-O-K",file=getOption("gmt.file"))Argumentsx data frame,matrix,orﬁlename containing the data to be plotted.cmd string of arguments passed to GMT shell command pstext.fileﬁlename where the map is saved.DetailsThe cmd argument can be used to specify the data columns:Format(cmd)Data columns Notes-J-R-O-K Lon Lat Text(default)-J-R-F+f+a+j-O-K Lon Lat Font Angle Justify Text(see example) If x is aﬁlename,the data should be tabular with or without a header,separated by commas or whitespace.Theﬁrst line is interpreted as header if theﬁrst non-whitespace character is not minus, point,or number.See the GMT documentation for details on pstext and other GMT commands.ValueNULL,but the map is annotated and saved in postscript format.The temporary GMT inputﬁle‘text.gmt’is saved in directory dirname(tempdir()),for the user to view or edit.It is later removed by psclose().See AlsoSimilar to text and postscript in native R graphics.gmt,pscoast,psxy,pstext,psbar,and psclose work together to draw maps.gmt-package gives an overview of the package.14psxy Examples##Not run:#Draw map and save as"map.ps"in current working directorygmt(demo.par)pscoast(demo.coast)psxy(demo.xy)pstext(demo.text,"-J-R-F+f+a+j-O-K")psbar(demo.bar,ref=66)psclose()#See directory gmt/example for details##End(Not run)psxy Add Lines/Symbols to GMT MapDescriptionCall GMT to add lines/symbols to a map and save in postscript format.Usagepsxy(x,cmd="-J-R-Scp-W2p-O-K",file=getOption("gmt.file"))Argumentsx data frame,matrix,orﬁlename containing the data to be plotted.cmd string of arguments passed to GMT shell command psxy.fileﬁlename where the map is saved.DetailsThe data are arranged in two(Lon,Lat)or more columns,depending on the-S argument.If x is aﬁlename,the data should be tabular with or without a header,separated by commas or whitespace.Theﬁrst line is interpreted as header if theﬁrst non-whitespace character is not minus, point,or number.See the GMT documentation for details on psxy and other GMT commands.ValueNULL,but the map is annotated and saved in postscript format.The temporary GMT inputﬁle‘xy.gmt’is saved in directory dirname(tempdir()),for the user to view or edit.It is later removed by psclose().See AlsoSimilar to points,lines,and postscript in native R graphics.gmt,pscoast,psxy,pstext,psbar,and psclose work together to draw maps.gmt-package gives an overview of the package.Examples##Not run:#Draw map and save as"map.ps"in current working directorygmt(demo.par)pscoast(demo.coast)psxy(demo.xy)pstext(demo.text,"-J-R-F+f+a+j-O-K")psbar(demo.bar,ref=66)psclose()#See directory gmt/example for details#Simple map and circlespscoast("-JM12c-R158/192/-42/-8-Di-Gdarkgreen-B10f5-A2000-K","quake.ps")psxy(data.frame(lon=quakes$lon,lat=quakes$lat,mag=10^quakes$mag/1e6),"-J-R-Scp-W0.5p,red-O","quake.ps")##End(Not run)r2gmt Prepare Data for GMTDescriptionThis internal function reads data,from aﬁlename or R object,and writes them to a GMT inputﬁle. Usager2gmt(x,datafile,append=FALSE)Argumentsx data frame,matrix,orﬁlename containing the data to be written to a temporary ﬁle.datafileﬁlename where the data will be written in GMT format,tab-separated without header.append whether data should be appended to an existingﬁle,separating segments with ">"lines.DetailsIf x is aﬁlename,the data should be tabular with or without a header,separated by commas or whitespace.Theﬁrst line is interpreted as header if theﬁrst non-whitespace character is not minus, point,or number.ValueThe data frame that was written to datafile.Noter2gmt is like write.table,except it allows x to be aﬁlename,and appends tables with the GMT> separator.It is mainly called by other functions,but users mayﬁnd it useful for writing input data for GMT commands.See Alsoscan,read.table,write,write.table.gmt-package gives an overview of the package.ExamplesLonLat1<-data.frame(Lon=1:3,Lat=4:6)LonLat2<-data.frame(Lon=7:8,Lat=9:10)##Not run:r2gmt(LonLat1,"temp.gmt")r2gmt(LonLat2,"temp.gmt",append=TRUE)##End(Not run)Indexas.character,9as.numeric,3barplot,10deg2num,2,3,9demo.bar,2demo.bar(gmt.demo),6 demo.coast,2demo.coast(gmt.demo),6 demo.par,2demo.par(gmt.demo),6 demo.text,2demo.text(gmt.demo),6 demo.xy,2demo.xy(gmt.demo),6dev.off,11diff,4geodist,2,4gmt,2,5,10–13,15gmt-package,2gmt.demo,6gmt.system,2,7lines,15num2deg,2,3,8options,5plot,12points,15postscript,10–13,15 psbar,2,5,9,11–13,15 psclose,2,5,10,11,12,13,15 pscoast,2,5,10,11,12,13,15 pstext,2,5,10–12,13,15 psxy,2,5,10–13,14r2gmt,2,15read.table,16scan,16system,8text,13Trig,4write,16write.table,16writeLines,8 17。

postinst 语法

postinst 语法postinst语法是Debian和Ubuntu等Linux操作系统中用于软件包安装脚本的一种格式。

在本文中，我们将详细介绍postinst语法，并提供一些示例来帮助读者更好地理解和使用该语法。

一、什么是postinst语法？postinst是一个由shell脚本编写的可执行文件，用于在软件包安装期间执行特定的操作。

它通常包含在Debian软件包的DEBIAN 目录中，并在软件包安装过程中自动执行。

postinst脚本可以用于完成各种任务，例如创建用户、设置环境变量、配置文件等。

二、postinst语法的基本结构postinst脚本的基本结构如下所示：#!/bin/shset -e# 配置变量和函数# 主要安装逻辑exit 0在上面的示例中，首先指定了脚本的解释器为/bin/sh，并通过set -e命令启用了错误检查。

这意味着如果脚本中的任何命令失败，脚本将立即退出。

接下来，可以定义一些变量和函数来帮助执行安装过程中的特定任务。

最后，在主要安装逻辑部分编写实际的安装代码。

最后，通过exit 0命令退出脚本。

三、postinst语法的常用命令和选项postinst脚本中可以使用许多命令和选项来执行各种任务。

下面列出了一些常用的命令和选项：1. dpkg-query命令：可以使用dpkg-query命令来查询已安装的软件包信息。

例如，dpkg-query -S <文件路径>将返回包含指定文件的软件包名称。

2. dpkg-divert命令：可以使用dpkg-divert命令来重定向文件。

例如，dpkg-divert --add --rename --divert <目标文件>.divert <目标文件>将重命名目标文件并创建一个新的文件，以便安装过程不会覆盖原始文件。

3. update-alternatives命令：可以使用update-alternatives命令来管理系统中的可替代项。

GMT4安装手册(windows)

GMT4.5.9绘图工具安装本软件运行需要安装4.5.9版本GMT（4.5.11版本、5.0版本及5以上版本可能会出现不兼容情况）；为了使用海岸线数据需要安装GSHHG；生成ps文件需要ghotstscript安装，查看ps文件需要安装gsView并设置中文显示。

1.1.1.GMT4.5.9安装（1）在安装文件“GMT”文件夹内，双击“gmt-4.5.9_install32.exe”。

在打开的GMT4安装向导中，选择“Next”。

（2）选择GMT安装的路径，推荐采用默认路径。

设置好后，点击“Next”。

（3）选择“Install”。

1.1.2.GSHHG安装（1）在安装文件GMT文件夹内，双击“gshhg-2.2.2_install.exe”。

（2）在打开的安装向导中，点击“Next”。

（3）将GSHHG安装在gmt同一路径下，推荐采用默认路径。

设置好后，点击“Next”。

（4）点击“Install”，开始安装。

1.1.3.Ghostscript安装查看ps文件，还需要安装支持各种文档编码格式的Ghostscipt软件。

若未能安装Ghostscipt，在打开GSView软件时会收到下载提醒。

本“GMT”文件夹内包含了gs9.07版本安装文件gs907w32.exe.点击“Next”，本软件的安装目录也不做要求。

1.1.4.GSView安装（1）在“GMT”文件夹内，双击“gsv49w32.exe”进行安装。

（2）在弹出的安装向导对话框中，点击“Setup”。

（3）选择安装语言，点击“English”。

（4）在弹出的安装提示中，选择“Next”。

（5）在弹出的安装提示中，选择“Next”。

（6）在弹出的安装提示中，选择“Next”。

（7）之后，选择软件安装路径，这里不作要求。

点击“Next”。

（8）确认路径，点击“Next”。

（9）完成安装，点击“Finish”。

1.1.5.GMT中文支持设置（1）添加中文字体到GMT字体列表在GMT安装目录下（默认为C:\programs\GMT4），\share\pslib下PS_font_info.d文件中按一下格式添加中文字体：结果如下图：（2）查看中文字体的编号在cmd命令窗口运行pstext -L，查看中文字体的编号。

GMT主要命令解释

The –B option–Bxinfo[/yinfo][:."title string":][W|w][E|e][S|s][N|n]该选项用选定的标志间隔指定了待绘制地图的边界。

xinfo 和 yinfo 的格式为[atick[m|c]][ftick[m|c]][gtick[m|c]][l|p][:"axis label":][:,"unit label":]这里a, f, 和g 分别为刻度标记、刻度和网格的间隔。

m|c 选择分（m）或秒（c）。

缺省设置是 4 个边界全部绘制（表示为W, E, S, N）。

通过重写代码来改变缺省设置（如, WSn）。

大写字母（如, W）说明同时画轴和轴的刻度标记。

轴标记的格式受你的.gmtdefaults 文件中的DEGREE_FORMAT 控制。

线性笛卡尔投影轴。

长刻度线为标记刻度，短刻度线为刻度间隔，轴名用―:‖括起，刻度单位―%‖用―: ,% :‖表示。

我们使用–R0/12/0/1 –JX3/0.4 –Ba4 f2 g1:Frequency::,%:.The –c option–c 选项指定绘图拷贝数。

[缺省为1]The –H option–H 选项让GMT 知道输入数据文件具有一个（缺省）或多个头纪录。

如果多于1 个头纪录，一定要在–H选项后指定数字。

如–H4。

见图1.9。

pscoast -R119/122/30.5/33 -Jm4c -B5The –J? options选择地图投影。

–J 后面的代码决定投影类型。

以用户单位指定地图宽度（或轴长）。

GMT 所允许的投影见图1.10墨卡托投影（Mercator: ）–JM width.正形投影（Orthographic: ）–JG lon0/lat0/width. lon0/lat0 指定投影中心Albers 圆锥投影（Albers conic: ）–JB lon0/lat0/lat1/lat2/width. 给出投影中心和2 条标准平行线Eckert 投影（Eckert IV and VI: ）–JK[f|s] lon0/width. 给出中央子午线线性投影（Linear: ）–JX width/height. 给出绘图宽度[和高度]。

GMT笔记-GMT常用画图等命令.

命令：backtracker Forward and Backward flowlines and hotspot tracksbinlegs维护GMT–SYSTEM 索引文件blockmean L2 （x,y,z）table data filter/decimatorblockmedian L1 （x,y,z）table data filter/decimatorblockmode Mode estimate （x,y,z）table data filter/decimatorcps把shell scripts 和数据文件插入到Complete PostScrip（t CPS）文件cpsdecode由Complete PostScript（CPS）文件提取shell scripts 和数据文件cpsencode把shell scripts 和数据文件插入到Complete PostScript（CPS）文件dat2gmt把一个ASCII 文件转换为二进位gmt 文件filter1d Filter 1 –D table data sets（time series）fitcircle得到一组数据点的最佳拟合大圆或小圆gmt GMT 数据处理和显示软件包gmt2bin由gmt 文件创建二进制索引文件gmt2dat二进制gmt文件转换为ASCII文件gmtconvert ASCII 和二进制1–D表文件之间转换gmtconvert数据表文件格式转换gmtdefaults当前缺省设置列表gmtinfo得到单独航次的信息gmtlegs查明一个给定地区的航次gmtlist由<legid>.gmt 文件提取数据gmtmath使用Polish Notation 语法对表格数据进行代数操作gmtpath获得gmt 文件的完全路径gmtselect选择基于多空间标准的数据子集gmtset改变当前的.gmtdefaults 文件所选择的参数gmttrack一个航迹绘图程序grd2cpt由一个网格文件创建色谱表grd2xyz 2 –D 网格文件转换为表格数据grdclip限定网格数据组的z值范围grdcontour绘制2–D 网格数据组等值线grdcut由一个网格文件抽取一个子区就是裁剪掉一部分只要感兴趣的部分grdcut ingrid-G outgrid -R ranggrdedit改变一个2–D 网格文件的标头信息grdfft Perform operations on gridded files in the frequency domaingrdfilter Filter 2 –D gridded data sets in the space domaingrdgradient由网格文件计算方向梯度,得到另一个网格文件grdhisteq网格文件柱状图均化（equalization）grdimage由一个2–D 网格数据组创建影像图grdinfo得到网格文件的信息grdlandmask Create masking gridded files from shoreline data basegrdmask Reset grid nodes in/outside a clip path to constantsgrdmath对网格文件进行代数操作grdpaste沿一条公共边界合并多个网格文件grdproject把网格数据组投影到一个新的坐标系grdraster Extract subregion from a binary raster and write a grd file grdreformat把网格文件转换为另一种格式grdsample对一个2-D 网格文件重新采样创建一个新的网格grdtrack沿1-D 轨迹对一个2-D 网格数据组重新采样grdtrend对网格文件进行多项式趋势拟合grdvector绘制2-D 网格矢量区grdview由一个2-D 网格数据组创建3–D透视影像图grdvolume计算给定等值线以下的体积hotspotter Create CVA image from seamount flowlinesimg2mercgrd Extract region of img, preserving Mercator, save as grd makecpt创建色谱表文件makepattern Make GMT color pattern from b/w pattern or icon mapproject Transformation of coordinate systems for table datamgd77togmt Convert an MGD –77 ascii file to a binary gmt file minmax报告数据表文件极大极小值nearneighbor Nearest –neighbor gridding schemeoriginator Associate seamounts with hotspot point sourcesproject Project table data onto lines or great circlespsbasemap创建底图psclip Use polygon files to define clipping pathspscoast在图上绘制（和充填）海岸线、国界线和河流pscontour Contour or image raw table data by triangulationpscoupe Plot cross –sections of focal mechanisms.pshistogram绘制柱状图psimage Plot Sun rasterfiles on a mappslib v3.2 A PostScript based plotting librarypsmask Create overlay to mask out regions on mapspsmeca Plot focal mechanisms on mapspsmegaplot To create poster –size PostScript plots from page–size plot pspolar Plot polarities on the inferior focal half –sphere on maps psrose绘制扇形或玫瑰花图psscale在图上绘制灰度图例或彩色图例pssegy Create imagemasked postscript from SEGY filepssegyz Create imagemasked postscript from SEGY filepstext在图上绘制字串psvelo Plot velocity vectors, crosses, and wedges on mapspswiggle沿地图上一条规机绘制根据时间序列数据绘制曲线psxy在图上绘制符号、多边形和线条psxyz在图上绘制3–D符号、多边形和线条sample1d表格数据组重新采样spectrum1d Compute various spectral estimates from time –series splitxyz把xyz 文件分为几段surface一种连续弯曲的网格算法trend1d Fits polynomial or Fourier trends to y = f （x）seriestrend2d Fits polynomial trends to z = f （x,y）seriestriangulate Perform optimal Delauney triangulation and gridding参数：-B设定图边界尺标的间隔-H选项让GMT 知道输入数据文件具有一个（缺省）或多个头纪录。

MultiByteToWideChar

//Description:
// This function maps a character string to a wide-character (Unicode) string
//
//Parameters:
// lpcszStr: [in] Pointer to the character string to be converted
BOOL WCharToMByte(LPCWSTR lpcwszStr, LPSTR lpszStr, DWORD dwSize)
{
DWORD dwMinSize;
dwMinSize = WideCharToMultiByte(CP_OEMCP,NULL,lpcwszStr,-1,NULL,0,NULL,FALSE);
// lpszStr: [out] Pointer to a buffer that receives the translated string.
// dwSize: [in] Size of the buffer
//
//Return Values:
// TRUE: Succeed
// FALSE: Failed
BOOL MByteToWChar(LPCSTR lpcszStr, LPWSTR lpwszStr, DWORD dwSize)
{
// Get the required size of the buffer that receives the Unicode
// string.
wchar_t *pwText;
pwText = new wchar_t[dwNum];
if(!pwText)

设计师必备常用网站

1、设计在线（新闻绝对是业界最及时的）2、视觉中国（新新设计师的代表性网站，有些杂，速度也不快）3、淘图网(高质量素材网站，更新很快) 4、设计师家园（建筑与室内设计师非去不可）5、亚洲CI网（网站本身不怎么样，但是论坛非常不错，高手云集）6、V6DP（论坛，年轻设计师的天下）7、（云集非常多大师作品以及他们的网站导航）8、Arting365（有些不错的独门信息及作品）9、台湾创意设计中心（了解台湾设计的必去站点）2、（云集非常多大师作品以及他们的网站导航，2004年底成立）3、鲜创意（综合性网站，2005年初成立）4、数字艺术（搜狐数码频道的分栏目，2005年初成立）5、涂鸦王国（自由开方性质的网上涂鸦站点，2004年底成立）6、中国美术家（非常全面，但是速度非常慢，2005年初成立）7、中国CG联盟（数字内容创作、电影电视制作，视觉娱乐设计，2004年底成立）8、火星时代动画网（不说了，看题目就知道，2004年成立）9、视觉同盟网（综合性网站，2004年成立）10、中国广告协会（不说了，看题目就知道，2005年初成立）设计师应该看得40个英文网站/工业设计/dmi/html/index.htm设计管理http://www.ddc.dk/丹麦工业设计/index-en.htm设计资讯/最新产品设计介绍/最新卫浴产品设计介绍http://www.form.de著名的设计杂志/现代艺术中心/default.asp著名的设计杂志/艺术设计博物馆官方/维多利亚和阿尔波特博物馆/design/statement/现代设计史方面的东东/阿莱西官方网站/通过网络走红的俄罗斯设计团体卡尔刘，内有大量工业设计草图/directory.htm家具设计史/teach/gothic/17，18世纪设计史/design设计博物馆，设计史.au/design/http://www.danskdesign.nu/index.php/language/en丹麦设计产品在线商店/设计与情感，设计交流/content/design/设计资讯http://www.droogdesign.nl荷兰的设计组织http://www.forma.co.hu/shop/设计商店，图多/这个要是不知道就拖出去枪毙1个小时/index.htm美国商业设计组织/home.htm设计与技术/著名的工业设计师/links/丹麦家具设计/当红设计明星http://www.uusimuoto.fi/内容很丰富可惜不是英文的//汽车设计/瑞士设计/blog/en/高科技趋势http://www.jidpo.or.jp/en/日本工业设计组织/设计公司/斯堪的纳维亚设计/index.jsp建筑环艺设计/视觉传达/yanagi/indexe.shtml日本当代设计大师柳宗理的个人网站/ideo.图像谷：CG-SKY：火星时代：视觉中国：设计联盟：点燃灵感：遥远理想：5D互动多媒体：国外CG门户：插画中国：AI俱乐部：蓝色理想：形色主义：闪客帝国：闪吧论坛：/bbs闪盟在线:中国UI设计网：素材精品屋：爱图网: /中国设计在线：365设计在线：/bbs数码艺术设计：/bbs七色鸟设计空间：中国photoshop联盟：http://www.ph→素材/教程资源类←天极美工教室：五维教程网：天使梦工厂：滤镜论坛：/forum_show.cgi 在线论坛：/bbsIT解题小组：网页教学网：冰火谷论坛：/bbs中国WEB开发者网站：设计师Rokey个人站：理想帝国：/www专注于PS技术:/pstext.htm 照片处理网：我爱策划网：占线图库：百万图库：/photo图酷：酷字千寻：/index.asp有风的日子：壁纸酷：3Lian素材网：好素材：素材词典：网页制作大宝库：承志模板站：中国GIF论坛：临沂信息港网页素材：/sucaiku 七色风网页家园：模板天下：天下秀：温州女孩素材网：柠檬素材部屋：三维资源：/eden/bbs/715497博熙设计：非常图酷：构图网：矢量LOGO：/find1.htm美图：/s8搜图网：/cgi-bin/imagefolio.cgi 伊美姬网：桌面城市：中国图片库：夜叉图库：网页制作大宝库：流星雨矢量图库：中华图库：/photo热线素材精品：壁纸好莱坞：Wincustomize：axialis：Eboy Home：ICON List：/iconICON中国：图标集中营：馨雅壁纸：/cn/index.php海龙文化：/civilization/index.asp非鸟矢量图库：/ai闪吧矢量图库：/index_vector.asp流星雨矢量图库：梧桐树下矢量图库：/flash/2000slt广之源矢量素材：/Subject.asp素材网：猫猫图片酷：/ai国外一个站点（wmf格式）：http://www.clip-art.ru闪客兵团矢量图库：/index6.htm闪客浪心：/flash.htm萤窗小语---矢量图库：/siliang.htm桌面城市：/scatalog.asp?scatalog=54海阔天空：/html/material-2.htm闪光玻璃：/index/slt.htm娱乐媒质传播网：/vectorsflash爱好者：/ziyuan/shiliang/index.htm闪客堡部分矢量资源：/download.php?sort=1corel公司（cdr格式）：/freebies/clipart.htm搜狐网素材版：/webcourse/webmon...aiku/index.html闪客天地Flash矢量资源（fla格式）：/vector.asp【推荐】电脑美术设计中的色彩表现：/443/1928943.shtml蓝精灵flash互动网矢量资源：/flashfla/flashvl.htm闪盟矢量图库（swf格式）：/home/read.php?qid=vector21点工作室矢量资源：/resource/bitmap/bmp.html平面设计与桌面出版--矢量素材区：http://202.103.233.162/cd/downhtml/down-1.htm asp深泽直→国外Photoshp教程资源←Text ，一流的网页设计师,教程和设计作品比得不看,还有PSD下载啊!/index.htm 教程和资源都不用说了,VERY GOOD!!!/tutorialsphotoshop.php 教程和链接都专家级水平老牌的PHOTOSHOP站点,现在教程正在更新,又有惊喜了!!/index1.htm专家级的网页设计作品,教程不得不看,网页作品是精品中的精品,设计师必去的网站!http://misery.subnet.at/v10/words/news_2k01.html个性化的设计师,平面作品和教程都个性化,好不好,看看便会给它迷倒!想找PSD素材的必到地方,教程不错但更新慢,但PHOTOSHOP资源是一般资源网站不能比的!/testa/tutorials.html年轻的设计师和那与年龄不相符的设计作品和教程,简直专家级水平!http://www.planet-tutorials.de/tutorials_photoshop.php很多的教程....老牌的PHOTOSHOP站点,常有惊喜的作品!老牌的PHOTOSHOP站点,学习PHOTOSHOP的必去的网站!比较多的作品/files/freebies/index.html很多的教程和资源!/vbrush/main.html一个专门的PHOTOSHOP笔刷资源网,就它就够了!一个用色很讲究的设计师网站,看过便知什么是用色的协调,知道我们的水平和他们有多远,强烈推荐!!/down/soft/636.htm -----Adobe photoshop cs 8.0完美破解汉化版人与ideo 都应该知道吧。

GMT笔记-GMT常用画图等命令

命令：backtracker Forward and Backward flowlines and hotspot tracksbinlegs维护GMT–SYSTEM 索引文件blockmean L2 （x,y,z）table data filter/decimatorblockmedian L1 （x,y,z）table data filter/decimatorblockmode Mode estimate （x,y,z）table data filter/decimatorcps把shell scripts 和数据文件插入到Complete PostScrip（t CPS）文件cpsdecode由Complete PostScript（CPS）文件提取shell scripts 和数据文件cpsencode把shell scripts 和数据文件插入到Complete PostScript（CPS）文件dat2gmt把一个ASCII 文件转换为二进位gmt 文件filter1d Filter 1 –D table data sets（time series）fitcircle得到一组数据点的最佳拟合大圆或小圆gmt GMT 数据处理和显示软件包gmt2bin由gmt 文件创建二进制索引文件gmt2dat二进制gmt文件转换为ASCII文件gmtconvert ASCII 和二进制1–D表文件之间转换gmtconvert数据表文件格式转换gmtdefaults当前缺省设置列表gmtinfo得到单独航次的信息gmtlegs查明一个给定地区的航次gmtlist由<legid>.gmt 文件提取数据gmtmath使用Polish Notation 语法对表格数据进行代数操作gmtpath获得gmt 文件的完全路径gmtselect选择基于多空间标准的数据子集gmtset改变当前的.gmtdefaults 文件所选择的参数gmttrack一个航迹绘图程序grd2cpt由一个网格文件创建色谱表grd2xyz 2 –D 网格文件转换为表格数据grdclip限定网格数据组的z值范围grdcontour绘制2–D 网格数据组等值线grdcut由一个网格文件抽取一个子区就是裁剪掉一部分只要感兴趣的部分grdcut ingrid-G outgrid -R ranggrdedit改变一个2–D 网格文件的标头信息grdfft Perform operations on gridded files in the frequency domaingrdfilter Filter 2 –D gridded data sets in the space domaingrdgradient由网格文件计算方向梯度,得到另一个网格文件grdhisteq网格文件柱状图均化（equalization）grdimage由一个2–D 网格数据组创建影像图grdinfo得到网格文件的信息grdlandmask Create masking gridded files from shoreline data basegrdmask Reset grid nodes in/outside a clip path to constantsgrdmath对网格文件进行代数操作grdpaste沿一条公共边界合并多个网格文件grdproject把网格数据组投影到一个新的坐标系grdraster Extract subregion from a binary raster and write a grd file grdreformat把网格文件转换为另一种格式grdsample对一个2-D 网格文件重新采样创建一个新的网格grdtrack沿1-D 轨迹对一个2-D 网格数据组重新采样grdtrend对网格文件进行多项式趋势拟合grdvector绘制2-D 网格矢量区grdview由一个2-D 网格数据组创建3–D透视影像图grdvolume计算给定等值线以下的体积hotspotter Create CVA image from seamount flowlinesimg2mercgrd Extract region of img, preserving Mercator, save as grd makecpt创建色谱表文件makepattern Make GMT color pattern from b/w pattern or icon mapproject Transformation of coordinate systems for table datamgd77togmt Convert an MGD –77 ascii file to a binary gmt file minmax报告数据表文件极大极小值nearneighbor Nearest –neighbor gridding schemeoriginator Associate seamounts with hotspot point sourcesproject Project table data onto lines or great circlespsbasemap创建底图psclip Use polygon files to define clipping pathspscoast在图上绘制（和充填）海岸线、国界线和河流pscontour Contour or image raw table data by triangulationpscoupe Plot cross –sections of focal mechanisms.pshistogram绘制柱状图psimage Plot Sun rasterfiles on a mappslib v3.2 A PostScript based plotting librarypsmask Create overlay to mask out regions on mapspsmeca Plot focal mechanisms on mapspsmegaplot To create poster –size PostScript plots from page–size plot pspolar Plot polarities on the inferior focal half –sphere on maps psrose绘制扇形或玫瑰花图psscale在图上绘制灰度图例或彩色图例pssegy Create imagemasked postscript from SEGY filepssegyz Create imagemasked postscript from SEGY filepstext在图上绘制字串psvelo Plot velocity vectors, crosses, and wedges on mapspswiggle沿地图上一条规机绘制根据时间序列数据绘制曲线psxy在图上绘制符号、多边形和线条psxyz在图上绘制3–D符号、多边形和线条sample1d表格数据组重新采样spectrum1d Compute various spectral estimates from time –series splitxyz把xyz 文件分为几段surface一种连续弯曲的网格算法trend1d Fits polynomial or Fourier trends to y = f （x）seriestrend2d Fits polynomial trends to z = f （x,y）seriestriangulate Perform optimal Delauney triangulation and gridding参数：-B设定图边界尺标的间隔-H选项让GMT 知道输入数据文件具有一个（缺省）或多个头纪录。

GMT笔记-GMT常用画图等命令

命令：backtracker Forward and Backward flowlines and hotspot tracksbinlegs维护GMT–SYSTEM 索引文件blockmean L2 （x,y,z）table data filter/decimatorblockmedian L1 （x,y,z）table data filter/decimatorblockmode Mode estimate （x,y,z）table data filter/decimatorcps把shell scripts 和数据文件插入到Complete PostScrip（t CPS）文件cpsdecode由Complete PostScript（CPS）文件提取shell scripts 和数据文件cpsencode把shell scripts 和数据文件插入到Complete PostScript（CPS）文件dat2gmt把一个ASCII 文件转换为二进位gmt 文件filter1d Filter 1 –D table data sets（time series）fitcircle得到一组数据点的最佳拟合大圆或小圆gmt GMT 数据处理和显示软件包gmt2bin由gmt 文件创建二进制索引文件gmt2dat二进制gmt文件转换为ASCII文件gmtconvert ASCII 和二进制1–D表文件之间转换gmtconvert数据表文件格式转换gmtdefaults当前缺省设置列表gmtinfo得到单独航次的信息gmtlegs查明一个给定地区的航次gmtlist由<legid>.gmt 文件提取数据gmtmath使用Polish Notation 语法对表格数据进行代数操作gmtpath获得gmt 文件的完全路径gmtselect选择基于多空间标准的数据子集gmtset改变当前的 .gmtdefaults 文件所选择的参数gmttrack一个航迹绘图程序grd2cpt由一个网格文件创建色谱表grd2xyz 2 –D 网格文件转换为表格数据grdclip限定网格数据组的z值范围grdcontour绘制2–D 网格数据组等值线grdcut由一个网格文件抽取一个子区就是裁剪掉一部分只要感兴趣的部分grdcut ingrid-G outgrid -R ranggrdedit改变一个2–D 网格文件的标头信息grdfft Perform operations on gridded files in the frequency domain grdfilter Filter 2 –D gridded data sets in the space domaingrdgradient由网格文件计算方向梯度,得到另一个网格文件grdhisteq网格文件柱状图均化（equalization）grdimage由一个2–D 网格数据组创建影像图grdinfo得到网格文件的信息grdlandmask Create masking gridded files from shoreline data base grdmask Reset grid nodes in/outside a clip path to constants grdmath对网格文件进行代数操作grdpaste沿一条公共边界合并多个网格文件grdproject把网格数据组投影到一个新的坐标系grdraster Extract subregion from a binary raster and write a grd file grdreformat把网格文件转换为另一种格式grdsample对一个2-D 网格文件重新采样创建一个新的网格grdtrack沿1-D 轨迹对一个2-D 网格数据组重新采样grdtrend对网格文件进行多项式趋势拟合grdvector绘制2-D 网格矢量区grdview由一个2-D 网格数据组创建3–D透视影像图grdvolume计算给定等值线以下的体积hotspotter Create CVA image from seamount flowlinesimg2mercgrd Extract region of img, preserving Mercator, save as grd makecpt创建色谱表文件makepattern Make GMT color pattern from b/w pattern or icon mapproject Transformation of coordinate systems for table data mgd77togmt Convert an MGD –77 ascii file to a binary gmt file minmax报告数据表文件极大极小值nearneighbor Nearest –neighbor gridding schemeoriginator Associate seamounts with hotspot point sourcesproject Project table data onto lines or great circlespsbasemap创建底图psclip Use polygon files to define clipping pathspscoast在图上绘制（和充填）海岸线、国界线和河流pscontour Contour or image raw table data by triangulation pscoupe Plot cross –sections of focal mechanisms.pshistogram绘制柱状图psimage Plot Sun rasterfiles on a mappslib v3.2 A PostScript based plotting librarypsmask Create overlay to mask out regions on mapspsmeca Plot focal mechanisms on mapspsmegaplot To create poster –size PostScript plots from page–size plot pspolar Plot polarities on the inferior focal half –sphere on mapspsrose绘制扇形或玫瑰花图psscale在图上绘制灰度图例或彩色图例pssegy Create imagemasked postscript from SEGY filepssegyz Create imagemasked postscript from SEGY filepstext在图上绘制字串psvelo Plot velocity vectors, crosses, and wedges on mapspswiggle沿地图上一条规机绘制根据时间序列数据绘制曲线psxy在图上绘制符号、多边形和线条psxyz在图上绘制3–D符号、多边形和线条sample1d表格数据组重新采样spectrum1d Compute various spectral estimates from time –series splitxyz把xyz 文件分为几段surface一种连续弯曲的网格算法trend1d Fits polynomial or Fourier trends to y = f （x）seriestrend2d Fits polynomial trends to z = f （x,y）seriestriangulate Perform optimal Delauney triangulation and gridding参数：-B设定图边界尺标的间隔-H选项让GMT 知道输入数据文件具有一个（缺省）或多个头纪录。

GMT4安装手册(windows)

GMT4.5.9绘图工具安装本软件运行需要安装4.5.9版本GMT（4.5.11版本、5.0版本及5以上版本可能会出现不兼容情况）；为了使用海岸线数据需要安装GSHHG；生成ps文件需要ghotstscript安装，查看ps文件需要安装gsView并设置中文显示。

1.1.1.GMT4.5.9安装（1）在安装文件“GMT”文件夹内，双击“gmt-4.5.9_install32.exe”。

在打开的GMT4安装向导中，选择“Next”。

（2）选择GMT安装的路径，推荐采用默认路径。

设置好后，点击“Next”。

（3）选择“Install”。

1.1.2.GSHHG安装（1）在安装文件GMT文件夹内，双击“gshhg-2.2.2_install.exe”。

（2）在打开的安装向导中，点击“Next”。

（3）将GSHHG安装在gmt同一路径下，推荐采用默认路径。

设置好后，点击“Next”。

（4）点击“Install”，开始安装。

1.1.3.Ghostscript安装查看ps文件，还需要安装支持各种文档编码格式的Ghostscipt软件。

若未能安装Ghostscipt，在打开GSView软件时会收到下载提醒。

本“GMT”文件夹内包含了gs9.07版本安装文件gs907w32.exe.点击“Next”，本软件的安装目录也不做要求。

1.1.4.GSView安装（1）在“GMT”文件夹内，双击“gsv49w32.exe”进行安装。

（2）在弹出的安装向导对话框中，点击“Setup”。

（3）选择安装语言，点击“English”。

（4）在弹出的安装提示中，选择“Next”。

（5）在弹出的安装提示中，选择“Next”。

（6）在弹出的安装提示中，选择“Next”。

（7）之后，选择软件安装路径，这里不作要求。

点击“Next”。

（8）确认路径，点击“Next”。

（9）完成安装，点击“Finish”。

1.1.5.GMT中文支持设置（1）添加中文字体到GMT字体列表在GMT安装目录下（默认为C:\programs\GMT4），\share\pslib下PS_font_info.d文件中按一下格式添加中文字体：结果如下图：（2）查看中文字体的编号在cmd命令窗口运行pstext -L，查看中文字体的编号。

各种字符编码间的转换方法：MultiByteToWideChar和MultiByteToWideChar(DOC)

MultiByteToWideChar和MultiByteToWideChar对各种字符编码间的转换函数原型：int WideCharToMultiByte(UINT CodePage,DWORD dwFlags,LPCWSTR lpWideCharStr,int cchWideChar,LPSTR lpMultiByteStr,int cbMultiByte,LPCSTR lpDefaultChar,LPBOOL lpUsedDefaultChar);此函数把宽字符串转换成指定的新的字符串，如ANSI，UTF8等，新字符串不必是多字节字符集。

参数：CodePage：指定要转换成的字符集代码页，它可以是任何已经安装的或系统自带的字符集，你也可以使用如下所示代码页之一。

参数说明：1、CodePage——指定要转换成的字符集代码页，它可以是任何已经安装的或系统自带的字符集，可选择以下代码页：CP_ACP //当前系统ANSI代码页CP_MACCP //当前系统Macintosh代码页CP_OEMCP //当前系统OEM代码页，一种原始设备制造商硬件扫描码CP_SYMBOL //Symbol代码页，用于Windows 2000及以后版本CP_THREAD_ACP //当前线程ANSI代码页，用于Windows 2000及以后版本CP_UTF7 //UTF-7，设置此值时lpDefaultChar和lpUsedDefaultChar都必须为NULLCP_UTF8 //UTF-8，设置此值时lpDefaultChar和lpUsedDefaultChar都必须为NULL用GetLocaleInfo 函数获取当前系统的代码页，936: 简体中文, 950: 繁体中文，949：韩文2、dwFlags——一般用0 就可以了指定是否转换成预制字符或合成的宽字符，对控制字符是否使用像形文字，以及怎样处理无效字符：MB_PRECOMPOSED //总是使用预制字符，即有单个预制字符时，就不会使用分解的基字符和不占空间字符。

GMT使用手册解析

datafiles
➢ CONVERT OR EXTRACT SUBSETS OF DATA: grd2xyz Convert 2-D gridded data to ASCII
table grdcut Cut a sub-region from a grd file grdpaste Paste together grdfiles along common
1. 在数字后直接加上单位符号，如：－X4c 2. 通过设置MEASURE_UN IT参数来实现。
GMT有许多预设的设定，要看GMT预设参数，在Linux环境下打 gmtdefaults -D。若要看使用者現正使用的参数，打gmtdefaults -L。
預在job內改参数，gmtset + 要改的参数，例如： gmtset ANOT_FONT_SIZE 12 （字型大小改成12） gmtset PAPER_MEDIA A4 (纸张大小设成A4，若要画A0或A1 的大图，则改为A0或A1) gmtset BASEMAP_TYPE FANCY (画出的地形图边框为一黑一白间隔，此为GMT预设值)
主要包括下面的几个文件：
netcdf.tar.Z GMT_suppl.tar.gz
bzip2.zip GMT_pdf.tar.gz GMT_scripts.tar.gz GMT_man.tar.gz GMT_high.tar.gz GMT_progs.tar.gz
triangle.tar.gz install_gmt GMT_tut.tar.gz GMT_ps.tar.gz GMT_progs.tar.gz GMT_web.tar.gz GMT_full.tar.gz
pswiggle Draw anomalies along track psxy Plot symbols, polygons, and lines in 2-D （二维折线图） psxyz Plot symbols, polygons, and lines in 3-D （三维折线图）

GMT笔记-GMT常用画图等命令

命令：backtracker Forward and Backward flowlines and hotspot tracksbinlegs维护GMT–SYSTEM 索引文件blockmean L2 （x,y,z）table data filter/decimatorblockmedian L1 （x,y,z）table data filter/decimatorblockmode Mode estimate （x,y,z）table data filter/decimatorcps把shell scripts 和数据文件插入到Complete PostScrip（t CPS）文件cpsdecode由Complete PostScript（CPS）文件提取shell scripts 和数据文件cpsencode把shell scripts 和数据文件插入到Complete PostScript（CPS）文件dat2gmt把一个ASCII 文件转换为二进位gmt 文件filter1d Filter 1 –D table data sets（time series）fitcircle得到一组数据点的最佳拟合大圆或小圆gmt GMT 数据处理和显示软件包gmt2bin由gmt 文件创建二进制索引文件gmt2dat二进制gmt文件转换为ASCII文件gmtconvert ASCII 和二进制1–D表文件之间转换gmtconvert数据表文件格式转换gmtdefaults当前缺省设置列表gmtinfo得到单独航次的信息gmtlegs查明一个给定地区的航次gmtlist由<legid>.gmt 文件提取数据gmtmath使用Polish Notation 语法对表格数据进行代数操作gmtpath获得gmt 文件的完全路径gmtselect选择基于多空间标准的数据子集gmtset改变当前的.gmtdefaults 文件所选择的参数gmttrack一个航迹绘图程序grd2cpt由一个网格文件创建色谱表grd2xyz 2 –D 网格文件转换为表格数据grdclip限定网格数据组的z值范围grdcontour绘制2–D 网格数据组等值线grdcut由一个网格文件抽取一个子区就是裁剪掉一部分只要感兴趣的部分grdcut ingrid-G outgrid -R ranggrdedit改变一个2–D 网格文件的标头信息grdfft Perform operations on gridded files in the frequency domaingrdfilter Filter 2 –D gridded data sets in the space domaingrdgradient由网格文件计算方向梯度,得到另一个网格文件grdhisteq网格文件柱状图均化（equalization）grdimage由一个2–D 网格数据组创建影像图grdinfo得到网格文件的信息grdlandmask Create masking gridded files from shoreline data basegrdmask Reset grid nodes in/outside a clip path to constantsgrdmath对网格文件进行代数操作grdpaste沿一条公共边界合并多个网格文件grdproject把网格数据组投影到一个新的坐标系grdraster Extract subregion from a binary raster and write a grd file grdreformat把网格文件转换为另一种格式grdsample对一个2-D 网格文件重新采样创建一个新的网格grdtrack沿1-D 轨迹对一个2-D 网格数据组重新采样grdtrend对网格文件进行多项式趋势拟合grdvector绘制2-D 网格矢量区grdview由一个2-D 网格数据组创建3–D透视影像图grdvolume计算给定等值线以下的体积hotspotter Create CVA image from seamount flowlinesimg2mercgrd Extract region of img, preserving Mercator, save as grd makecpt创建色谱表文件makepattern Make GMT color pattern from b/w pattern or icon mapproject Transformation of coordinate systems for table datamgd77togmt Convert an MGD –77 ascii file to a binary gmt file minmax报告数据表文件极大极小值nearneighbor Nearest –neighbor gridding schemeoriginator Associate seamounts with hotspot point sourcesproject Project table data onto lines or great circlespsbasemap创建底图psclip Use polygon files to define clipping pathspscoast在图上绘制（和充填）海岸线、国界线和河流pscontour Contour or image raw table data by triangulationpscoupe Plot cross –sections of focal mechanisms.pshistogram绘制柱状图psimage Plot Sun rasterfiles on a mappslib v3.2 A PostScript based plotting librarypsmask Create overlay to mask out regions on mapspsmeca Plot focal mechanisms on mapspsmegaplot To create poster –size PostScript plots from page–size plot pspolar Plot polarities on the inferior focal half –sphere on maps psrose绘制扇形或玫瑰花图psscale在图上绘制灰度图例或彩色图例pssegy Create imagemasked postscript from SEGY filepssegyz Create imagemasked postscript from SEGY filepstext在图上绘制字串psvelo Plot velocity vectors, crosses, and wedges on mapspswiggle沿地图上一条规机绘制根据时间序列数据绘制曲线psxy在图上绘制符号、多边形和线条psxyz在图上绘制3–D符号、多边形和线条sample1d表格数据组重新采样spectrum1d Compute various spectral estimates from time –series splitxyz把xyz 文件分为几段surface一种连续弯曲的网格算法trend1d Fits polynomial or Fourier trends to y = f （x）seriestrend2d Fits polynomial trends to z = f （x,y）seriestriangulate Perform optimal Delauney triangulation and gridding参数：-B设定图边界尺标的间隔-H选项让GMT 知道输入数据文件具有一个（缺省）或多个头纪录。

在地图上绘制波形

在地图上绘制波形将 SAC 波形数据绘制在地图上也算是常见的需求之⼀。

最常见的情况是，在地图上有⼀些台站，需要将波形画在台站的附近。

此问题的难点在于，在画地图时，-R 指定的是区域范围，即横轴是经度，纵轴是纬度，⽽pssac 在绘制波形时，横轴是时间，纵轴是振幅。

因⽽想要将波形放在特定的位置，不是⼀件容易的事情。

⽣成⽰例⽤ SAC ⾃带的数据，⽣成⼏个⽰例数据，供接下来使⽤:1 2SAC> dg sub tele ntkl.z nykl.z onkl.z sdkl.z SAC> w ntkl.z nykl.z onkl.z sdkl.z⽤saclst查看⼀下这⼏个波形数据的相关信息:1 2 3 4 5$ saclst b e stlo stla f *.zntkl.z 199.965 1600.95 -114.592 62.4797 nykl.z 199.962 1600.97 -74.53 44.5483 onkl.z 199.618 1600.62 -93.7022 50.8589 sdkl.z 199.098 1600.11 -104.036 44.1204先绘制地图在考虑如何把波形画在地图上之前，先得有⼀个地图。

⽤如下代码绘制⼀个地图：1 2 3 4 5 6 7 8 9 10 11 12 13#!/bin/bashJ=M15cR=-120/-60/30/65PS=map.pspsxy -J$J -R$R -T -K > $PS# 绘制海岸线pscoast -J$J -R$R -B10/10 -Ggray -K -O -A1000 >> $PS# 绘制台站位置saclst stlo stla f *.z | awk '{print $2, $3}' | psxy -J$J -R$R -Sa0.5c -Gblack -K -O >> $PS# 绘制台站名saclst stlo stla f *.z | awk '{print $2, $3,"15 0 0 TR", $1}' | pstext -J$J -R$R -D-0.1c/-0.1c -K -O >> $PS psxy -J$J -R$R -T -O >> $PS代码中使⽤了saclst和awk提取了每个波形所对应的台站经纬度，并通过管道传递给psxy和pstext命令。