您现在的位置是:课程

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

-->