用matlab怎么画视电阻率拟断面图,在MATLAB平台上实现可控源音频大地电磁反演数据三维可视化显示...

第29卷　增刊

物探化探计算技术

2007年10月

收稿日期6文章编号:1001—1749(2007)增刊(1)—0068—04

在MAT LAB 平台上实现可控源音频大地

电磁反演数据三维可视化显示

李晓昌

(中国地质科学院地球物理地球化学勘查研究所,河北廊坊　065000)

摘　要:这里运用MATLAB 的G U I 图形工具,实现了可控源音频大地电磁法(简称CS A MT)二维反演数据的三维可视化展示,从而提高了物探工作者成果解释的工作效率。关键词:MAT LAB;三维可视化;可控源音频大地电磁测深;反演数据;三维立体图中图分类号:P 63113+

25 文献标识码:A

0　前言

MAT LAB 原先是作为Matrix 实验室使用L I N 2PACK 和E I SP ACK 矩阵软件工具包的接口,后来逐渐发展成为集通用科学计算、图形交互、系统控制和程序语言设计为一体的商业软件[1～2]

。由交互式的图形(Figure )窗口工具,生成图形用户界面(G U I ),便于用户进行图形定制。国内物探工作者在绘制电磁测深反演数据的图形中,通常使用Sur fer 软件。但Surfer 软件是提供绘制二维图形的工具,为了制作三维可视化图形,需在Sur fer 界面上对每条测深数据的二维图形经过多次的旋转及变换后,再将这些多个图形拼凑成为三维可视化图形。其方法复杂,不便于资料解释。作者在本文通过使用MATLAB ,修改软件提供的色标,轻松地实现了对多条电磁测深反演数据的三维可视化图形展示,并可将图形旋转,便于从不同角度分析反演结果,显示不同深度(高程)切面的电阻率平面图。

1　三维可视化绘制原理

如图1所示,在三维坐标系统中,定义y 轴为测线,x 轴为测点,z 轴为深度或高程,这三个方向以米为单位。MATLAB 提供了多个函数来表示三

维数据。有的函数用三维来显示曲线,而有的函数

则负责绘制表面和构建框架[3～5]

。也可以用颜色

来表示第四维数据-对应某一个坐标(yxz )点的电阻率值。

三维图形与z =0平面相关的这个方向角度被

称为仰角,而与x =0这个平面有关的角度被称为方位角。默认的三维视角就是仰角为30°而方位角为37.5°。在三维图中可以通过改变视角,从不同角度观看三维图形。

2　三维可视化图形绘制

2.1　显示方法

电磁测深反演数据分布在三维空间的每个测

图1　三维坐标系统示意图

Fig .1　The sche m e of 3D coordinates

:2007-0-20