MATLAB代码记录

数据整理、格式转换、计算

生成反演的单元体网格空间

%% 生成模型空间
clc;
clear;
%模型边界信息
X_min = 89.5;
X_max = 140.5;
dx = 1;
Y_min = -10.5;
Y_max = 35.5;
dy = 1;
Z_min = -20;
Z_max = -800;
dz = 20;
%导入上顶面起伏数据【注意数据Z轴正方向朝向，此处以向下为正】
%三列 【经度 纬度 深度】向下为正
data_top_depth = load('XXXX.dat');
%起伏层数据底部平面深度【该模型仅在第一层数据存在起伏厚度，下部数据均为等厚层】
horizontal_layers_depth = -20;
%生成的模型框架名称
output = fopen('XXXXX.dat','w');
%生成模型
for i = 1:length(data_top_depth)
    %模型参数：网格点中心位置【经度 纬度 深度 厚度】
    fprintf(output,'%f %f %f %f \n',data_top_depth(i,1),data_top_depth(i,2),-(data_top_depth(i,3)-horizontal_layers_depth)/2,abs((-horizontal_layers_depth)-data_top_depth(i,3)));
end
for z = Z_min-dz:-dz:Z_max
    for y = Y_min:dy:Y_max
        for x = X_min:dx:X_max
            %模型参数：网格点中心位置【经度 纬度 深度 厚度】
            fprintf(output,'%f %f %f %f \n',x,y,z+(dz/2),dz);
        end
    end
end
fclose(output);