资源简介
使用matlab利用水平导数法进行边界识别，识别地质体的边界，识别构造单元，并与其他方法进行比对，比较各自的效果。有一定的地球物理意义
代码片段和文件信息
%% AS
clear;clc;
%% inputing magnetic anomaly data
[RC xR yR DT] = ReadSurferGrd（‘D:\Program Files\edge_detect\ModelData\complex.grd‘）;
dx=（xR（2）-xR（1））/（RC（1）-1）;
dy=（yR（2）-yR（1））/（RC（2）-1）;
%% adding noise
% noise_level=0;
% noise=2*rand（RC（2）RC（1））-1;
% noise=noise.*abs（DT）*noise_level/100;
% DT=DT+noise;
%% the 3 directions derivative of the magnetic anomaly data
[DT_xDT_y]=gradient（DTdxdy）;
RC_max=max（[RC（1） RC（2）]）;
npts=2*2^nextpow2（RC_max）;
DT_z=vertical（DTnptsRC（1）RC（2）dx）;
%% AS
AS=sqrt（DT_x.^2+DT_y.^2+DT_z.^2）;
%% MATLAB pl
川公网安备 51152502000135号
评论
共有 条评论