代码分享：面波数据快速成图

代码分享：面波数据快速成图

前言

目前，物探数据主要用surfer软件成图，surfer软件具有强大的插值和绘图功能，成图比较美观。但是，生产过程中大量的物探数据，依靠excel和surfer来成图耗费人力时间成本。本博文在MATLAB平台上开发了一套用于面波数据快速成图的小程序，仅供大家借鉴。

1、成图效果展示

1.1 原始图像

对面波数据采用geogiga软件反演，导出视横波数据，在matlab中编辑克里金插值算法的代码，输出图像。

1.2 高程转换

将地表GPS测量的高程，与面波探测的深度进行转换，得到真实的高程。

1.3 里程换算

将地表GPS测量的里程与高程，与面波探测的深度和水平距离进行换算，由于面波测点在地表不等距分布，因此里程也是不等间距分布，换算之后得到真实的高程与里程。

1.4 图像加工

为了得到比较美观的图像，在MATLAB中对图像进行加工。

2、数据读取与图像保存

2.1 读取面波视横波速度数据

选择数据文件夹。

% 读取面波数据
[FileName,PathName] = uigetfile('*.txt','请选择视横波速度文件',...
    'MultiSelect','on');
filename = strcat(PathName,FileName);
data = importdata(filename);
fprintf('\n读取视横波速度完成!\n请按任意键继续...\n');

提取数据，自编函数dealMBdata。

% 初始参数设置
% 最大深度
depth_max = 80;
% 插值点数
num_points = 40;

% 面波数据预处理
[points,vs_value,xlocation] = dealMBdata(data);

2.2 数据与图像保存

% 保存数据