20200718新增:制作元素含量不同深度对比图
1.Surfer下制作不同深度对比图
重点掌握制作不同深度对比图的方法
掌握快速统一等值线图样式的方法
掌握快速统一图件大小的方法
2.voxler制作不同深度对比图
重点掌握制作不同深度图件的方法
掌握添加图框、坐标轴的方法
掌握添加静态文本、图例的方法
掌握添加三维文本的方法。
一、【Surfer篇】数据二维可视化
1.数据处理-坐标转换
我们原始的采样点数据最常见的是地理坐标系-度分秒格式的数据,如下图所示。这种数据格式(文本型)是不会被Surfer(或其他GIS)软件所能识别的,所以在成图前,我们需要把这种格式的数据转换成Surfer(或其他GIS)软件所能识别的数据格式(数值型)。
转换前:
X |
Y |
Z |
取样深度 |
氨氮 |
120°37′8.35″ |
29°54′18.51″ |
27.93048 |
1 |
0.07 |
120°37′7.69″ |
29°54′19.48″ |
25.43048 |
3.5 |
0.2 |
120°37′7.83″ |
29°54′18.9″ |
7.330478 |
21.6 |
0.13 |
120°37′9.05″ |
29°54′18.8″ |
30.02867 |
1 |
0.27 |
120°37′8.39″ |
29°54′17.89″ |
26.02867 |
5 |
0.03 |
120°37′7.08″ |
29°54′16.65″ |
12.02867 |
19 |
0.13 |
我们可以利用课件里的常用小工具.xls这个EXCEL文件实现度分秒与度的坐标的互转。AB列输入十进制单位的坐标,CD列即可自动生成度分秒格式的坐标
FG列输入度分秒格式的坐标,HI列即可自动生成度分秒格式的坐标十进制单位的坐标
常用小工具.xls这个EXCEL文件同时提供中央经线换算功能。可以根据所在地的经度快速计算计算出其三度带(六度带)情况下的带号和中央子午线
利用常用小工具.xls把度分秒格式坐标转换成十进制坐标后,可以利用surfer坐标转换功能,把地理坐标转换成高斯投影坐标。
Surfer没有的坐标系可以通过导入方式导入,比如本次课程导入中央经线为120度的三度带坐标系。
转换后成果如下表所示
工程号 |
X |
Y |
氨氮 |
经度 |
纬度 |
经度 |
纬度 |
CH3 |
559781 |
3309759.08 |
0.07 |
120.6189852 |
29.9051406 |
120°37′8.35″ |
29°54′18.51″ |
CH4 |
559763 |
3309789.1 |
0.07 |
120.6188023 |
29.90541233 |
120°37′7.69″ |
29°54′19.48″ |
CH5 |
559767 |
3309771.12 |
0.1 |
120.6188426 |
29.9052499 |
120°37′7.83″ |
29°54′18.9″ |
CH6 |
559800 |
3309768.31 |
13.57 |
120.6191817 |
29.90522296 |
120°37′9.05″ |
29°54′18.8″ |
CH7 |
559782 |
3309740.24 |
8.29 |
120.618998 |
29.90497058 |
120°37′8.39″ |
29°54′17.89″ |
CS10 |
559747 |
3309701.74 |
3.14 |
120.6186331 |
29.90462503 |
120°37′7.08″ |
29°54′16.65″ |
CS11 |
559790 |
3309713.22 |
0.03 |
120.6190779 |
29.90472647 |
120°37′8.68″ |
29°54′17.02″ |
2.制作采样点位图并转KML
当我们把坐标格式转换为数值型后就可以利用Surfer处理数据了。必备的一张图就是我们的设计点位图(施工前)、采样点位图(施工后)。创建一个post map即可。其中X是经度,Y是纬度,标注内容为采样点号(设计点号)。底图可以叠加上场地范围、谷歌影像等基础数据(不同数据源图形套合方法)。
当我们需要去场地现场调查时,我们可以把点位图转换为KMZ/KML格式,经优化显示后可在谷歌地球或者奥维等软件上显示采样点如图所示。
优化前
优化后
3.超标点位图
根据检测值,依据行业规范,可以用不同大小和颜色的点来显示超标点位
4.污染异常等值线图
根据原始数据网格化插值后可以生成污染等值线图,
也可以根据污染值显示污染范围
5.地下水流向图
需要准备的数据工程号、X、Y、高程(孔口),静止水位(深度)、标高(潜水面高程,高程- 静止水位所得)。高程(孔口)除了实测外也可以根据场调地形网格文件获取(方法?)
工程号 |
X |
Y |
高程 |
静止水位 |
标高 |
CH3 |
559781 |
3309759.08 |
30.55965 |
2.7 |
27.8597 |
CH4 |
559763 |
3309789.1 |
31.93267 |
2.3 |
29.6327 |
CH5 |
559767 |
3309771.12 |
31.65389 |
1.4 |
30.2539 |
CH6 |
559800 |
3309768.31 |
29.04307 |
1 |
28.0431 |
CH7 |
559782 |
3309740.24 |
29.83162 |
1.7 |
28.1316 |
CS10 |
559747 |
3309701.74 |
29.63637 |
2.2 |
27.4364 |
CS11 |
559790 |
3309713.22 |
28.15636 |
1.4 |
26.7564 |
CS12 |
559833 |
3309703.54 |
23.89906 |
2.4 |
21.4991 |
CS13 |
559872 |
3309702.27 |
21.6722 |
1.2 |
20.4722 |
CS14 |
559893 |
3309709.2 |
20.90439 |
2.2 |
18.7044 |
潜水面及流向(瞎编数据,所以看起来比较混乱)
6.三维样品分布图(Surfer新功能)
需要准备的数据。工程号、X、Y、高程(孔口),取样深度、标高(高程- 取样深度所得)
X |
Y |
高程(孔口) |
取样深度 |
Z |
559802.560 1、电脑:访问腾讯课堂官网 ke.qq.com 查看我的课表或下载win/mac客户端听课 2、手机/平板:下载腾讯课堂APP, 进入学习页面听课 |