区域台站数据场图生成——学习版

 0. 数据下载

DEM数据：地理空间数据云 (gscloud.cn)

全国基础地理数据：全国地理信息资源目录服务系统 (webmap.cn)

1. 地理底图的制作

1.1数据下载

在“地理空间数据云”下载研究区的DEM数据，网站：地理空间数据云 (gscloud.cn)。

选择“高级检索”，选择所需要的数据集以及研究区范围，本次以经纬度来限定。

下载数据解压并添加至Arcmap中，接下来进行数据镶嵌。

1.2数据镶嵌

数据镶嵌的目的是消除图幅数据之间色彩的差异。

打开ArcToolbox，选择“数据管理工具”-“栅格”-“栅格数据集”-“镶嵌”，输入栅格选择所有图幅，选取一个栅格图幅作为目标栅格。本次选取58-06图幅，如下：

点击确定，进行数据镶嵌。

先前完成后可以将其他图幅移除，以免数据过多。

1.3数据裁剪

方法一：采用四角顶点的点数据生成面要素进行裁剪

本次数据处理需要的是一个经纬度范围内的图形。

经纬度转平面坐标步骤参考：arcgis经纬度转平面坐标教程 - 知乎 (zhihu.com)

图像裁剪参考：ArcGIS | 根据角点经纬度裁剪shp文件_裁剪线.shp-CSDN博客

方法二：采用经纬度范围进行裁剪

ArcToolbox-数据管理工具-栅格-栅格处理-裁剪，输入数据及经纬度范围进行裁剪，本次下载的数据为STRM_90m分辨率的dem数据，可以直接按照经纬度范围进行裁剪。

1.4底图绘制

将裁剪后的dem数据生成山体阴影图，若图像噪点过多，修改Z值，默认为1，但本次修改Z值为0.000005后生成的山阴图较为合适。

将裁剪的dem原始数据，修改色带以及透明度，叠加到山体阴影图之上，就能够较为直观的显示出地形的变化。

最终图（部分）如下：

可进行高程配色的修改调整，使图片更加美观。

1.5掩膜工具提取矢量数据

ArcToolbox-制图工具-掩膜工具-交叉图层掩膜

其中，掩膜图层选择裁剪范围，被掩膜图层选择需要裁剪的矢量图层。

2.GPS场图绘制

2.1前置操作

筛选经纬度范围内的GPS台站数据、excel处理数据得到相对的E&N分量。

2.2 ArcGIS处理

导入arcgis，右键“显示XY数据”；

将生成的图层右键-数据-导出，生成.shp格式的文件；

2.2.1计算矢量大小

右键-打开属性表，计算矢量大小：新建字段-命名为Field-改为双精度并设置小数位数-右键字段-字段计算器。

2.2.2箭头设置

右键台站点图层-属性-符号系统数量-分级符号-字段值选择Field。

模板-编辑符号-类型-字符标记符号-字体选择ESRI IGL Font20-选择箭头。生成无方向的箭头符号。

2.2.3箭头方向设置

右键台站点图层-属性-符号系统数量-分级符号-高级-旋转。

按表达式并构建表达式：( 180 /3.14 ) * Atn2( [R_Ve_l], [R_Vn_l] )。

根据需要调整符号系统。

可根据实际绘制的图形要求，右键-属性-符号系统-模板，来修改符号的初始旋转角度，本次设置90°，得到最终的GPS_velocity图（未绘制置信椭圆）。

2.2.4符号系统不随地图缩放变化

右键“图层”-属性-参考比例：使用当前比例尺-确定。

2.3GPSprofile绘制

2.3.1Matlab计算profile范围四角经纬度坐标

运行corner_position.m文件，在文件中更改需要绘制field提取范围的中心点纬度、中心点经度、数据宽度、长度及断层走向。运行文件即得到四点经纬度坐标结果data。

2.3.2ArcMap由点生成面

将data的四个点坐标输入到.xls的Excel文件中；

表格导入ArcMap，显示XY数据，转换为.shp格式的文件；

方法一：ArcToolbox-数据管理工具-要素-点集转线，首先将点要素转换为线要素，选定闭合线；第二步，使用要素转面将闭合线要素转为面。

方法二：在ArcCatalog创建面.shp文件，点击开始编辑-创建要素-绘制封闭面图形。

2.3.3计算范围内平行及垂直于断层的速度投影

需要的数据为：GPS_data，每列分别为经度、纬度、E向速度分量、N向速度分量、E向速度误差、N向速度误差以及相关度。

该步的目的是计算出平行及垂直与断层走向的E、N速度投影，筛选出剖面范围内的GPS数据点，生成figure：表示GPS数据点位与矩形范围的位置关系散点图。

2.3.4平行、垂直于断层速度投影的图件绘制

LPF_gps_profile进行平行、垂直于断层速度误差图的绘制，最终得到result表格数据，每列分别代表：GPS点到断层的距离（km）、平行于断层的速度投影（mm/yr）、平行于断层的速度投影误差（mm/yr）、垂直于断层的速度投影（mm/yr）、垂直于断层的速度投影误差（mm/yr）。

2.3.5误差处理及计算_Origin——平均值参考线及置信区间

双Y轴图的绘制：参考：Origin做双y轴图 - 知乎 (zhihu.com)

调整绘图比例：双击表格空白处，弹出窗口处调节

插入参考线：参考：OriginPro：如何在origin中画一条Y=0.3的水平线？方法3：你想要的方法【数据绘图】 - 知乎 (zhihu.com)

添加误差条带：双击坐标轴-参照线-左侧选择“水平”-细节-追加，选择参照线生成方法及填充，本次选择的统计均值及标准差。
显示设置：上个步骤中设置参照线对话框-标签-勾选显示-自定义-文本，进行自定义文本设置，例如：均值+标准差、Mean+SD…。

在上部及右侧添加坐标轴：右键图件中间空白区域-新图层轴-打开对话框-自定义-勾选需要添加的坐标轴-关联坐标轴刻度-确定。

一些绘图参考（有些不相关）：

Origin：误差带图的绘制 - 知乎 (zhihu.com)

Origin绘图（4）——带误差棒散点图 - 知乎 (zhihu.com)

Origin绘图|带填充区误差线图_origin画误差带图-CSDN博客

【SPSS绘图】误差图 - 知乎 (zhihu.com)     散点误差条图的绘制

最终话：待学习

01 GTM绘制velo

velo — GMT 中文手册 (gmt-china.org)

02 置信椭圆的概念及绘制

概念：置信椭圆（误差椭圆）详解-CSDN博客

03 python关于矢量场图的绘制

python 矢量场的构造、计算；通过 matplotlib 作图、分析、筛选符合要求的像素点_计算矢量场的散度 python plot-CSDN博客

这次的记录就先到这，如果有新东西下篇再开