基于聚类分析的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

numsMatrix2ws.m function numsMatrix2ws(Data,Layer,res_object,res_air) elev_layer = Data; Layer = sort(Layer); max_depest_layer = length(Layer); % Z为间隔,用于输出z方向网格间距 Z = Layer(2:end) - Layer(1:end-1); Z = Z(end:-1:1); res_object = log(res_object); res_air = log(res_air); out_ws = fopen('ws部分.txt','w+'); fprintf(out_ws,"%.2f\t",Z); fprintf(out_ws,"\n\n"); %循环获得该层模型电阻率数组并输出 for i = 1:max_depest_layer temp = elev_layer; test_value = max_depest_layer - i +1; % 标识层位的变量 temp(temp >= test_value) = res_object +