基于聚类分析的ModEM三维反演复杂地形网格剖分
胡士晖, 闵刚, 孙浥钦, 陈春江, 李春婷, 张志豪

Gridding of complex terrains based on cluster analysis for ModEM 3D inversion
HU Shi-Hui, MIN Gang, SUN Yi-Qin, CHEN Chun-Jiang, LI Chun-Ting, ZHANG Zhi-Hao
numsMatrix2cov.m
function numsMatrix2cov(Data,Layer)
elev_layer = Data;
Layer = sort(Layer); %正向排序簇中心数组
max_depest_layer = length(Layer);
for i = 1:max_depest_layer %循环获得该层协方差矩阵并输出
temp = elev_layer';
test_value = max_depest_layer - i +1;
temp(temp < test_value) = 0;
temp(temp >= test_value) = 1;
temp_space = [i,i];
writematrix(temp_space,'cov部分.xls','WriteMode','append');
writematrix(temp,'cov部分.xls','WriteMode','append');
end
fclose('all');
end