%%?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