您现在的位置是:课程
Surfer及Voxler软件在污染场地调查中的应用
2023-12-05 20:00课程 人已围观
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