地壳应变率计算的传统方式是将GPS站点组成三角形进行计算,因Matlab与GMT在构造Delaunay三角网时采用不同的算法,对于同一套数据生成结果有明显的差异。通过研究Matlab与GMT构建三角网结果输出的数据格式,利用Perl语言编写程序将两套结果进行转换,实现Matlab与GMT软件在绘制Delaunay三角网时的相互转化并构建同样的三角网结果,绘图时再结合GMT的地形底图,取得了较好的表达效果。
The traditional method for crustal strain calculation is constructing triangulation with all the GPS sites first,and then calculating the strain of every triangulation.If the different algorithms:Matlab and GMT are adopted,then for the same data set the different trianggulation outputs will be obtained.By the study of the data formats output from the constructed triangulation results with Matlab and GMT we have programmed the transform between both the results in the Peal language thus realized the transform to each other between the software Matlab and GMT in mapping Delaunay triangulation and constructed the same triangulation results.When in mapping we combined the geographic manuscriptmap,thus we have achieved a good expressing results.