matlab生成Z为常数的.grd格式数据

当数据的Z值为常数时,surfer不能进行网格化而直接得到全为均值的.grd数据(即得到如下数据)。
Z值为常数
使用以下函数可生成上述形式的数据。
1、读入.grd数据

function [data,line,list,X_MIN,X_MAX,Y_MIN,Y_MAX,Z_MIN,Z_MAX] = ReadGrd(filename)
%ReadGrd
%读grd文件
%data 是输出的数据(矩阵)
%line,list 是行数和列数
%X_MIN,X_MAX 为x方向的最小值和最大值
%Y_MIN,Y_MAX 为y方向的最小值和最大值
%Z_MIN,Z_MAX 为z方向的最小值和最大值
%filename是文件名
%第一个数据为西南角的数据

fid = fopen(filename,'r');
ID=fscanf(fid,'%s',1);
list=fscanf(fid,'%d',1);%列数
line=fscanf(fid,'%d',1);%行数
Y_MIN=fscanf(fid,'%f',1);%列
Y_MAX=fscanf(fid,'%f',1);
X_MIN=fscanf(fid,'%f',1);%行
X_MAX=fscanf(fid,'%f',1);
Z_MIN=fscanf(fid,'%f',1);
Z_MAX=fscanf(fid,'%f',1);
data(line,list) = 0.0;
for k=1:1:line
    for t=1:list
        data(k,t) = fscanf(fid,'%f',1);
    end
end
fclose(fid);
end

2、输出.grd数据

function [] = WriteGrd(filename,data,list,line,X_MIN,X_MAX,Y_MIN,Y_MAX)
%writedGrd
%写grd文件
%data 是输入的数据(矩阵)
%line,list 是行数和列数
%X_MIN,X_MAX 为x方向的最小值和最大值
%Y_MIN,Y_MAX 为y方向的最小值和最大值
%Z_MIN,Z_MAX 为z方向的最小值和最大值
%filename是文件名

disp('数据写入文件-----------');

Z_MAX = max(max(data));
Z_MIN = min(min(data));
fid = fopen(filename,'wt');
fprintf(fid,'DSAA\n');
fprintf(fid,'%d ',list);%列数
fprintf(fid,'%d\n',line);%行数
fprintf(fid,'%f ',X_MIN);%x,y,z的变换范围
fprintf(fid,'%f\n',X_MAX);
fprintf(fid,'%f ',Y_MIN);
fprintf(fid,'%f\n',Y_MAX);
fprintf(fid,'%f ',Z_MIN);
fprintf(fid,'%f\n',Z_MAX);
j=1;
for k=1:1:line
    for t=1:list
        fprintf(fid,'%f ',data(j,1));
        j=j+1;
    end
    fprintf(fid,'\n');
end
fclose(fid);
end

3、设置Z为均值

clc;
clear;
filename_in = '*******************';
%输入一个与目标数据.grd格式参数相同的.grd数据
filename_out1 = '************************';
%目标.grd格式数据的保存路径
[JMDZ,Nx,Ny,X_MIN,X_MAX,Y_MIN,Y_MAX,DZ_MIN,DZ_MAX] = ReadGrd(filename_in);
data(Nx*Ny,1) = 0.0;
for i = 1:Nx
    for j = 1:Ny  
        data((i-1)*Ny + j ,1) = -2000;
        %均值为-2000
    end
end
WriteGrd_change(filename_out1,data,Ny,Nx,Y_MIN,Y_MAX,X_MIN,X_MAX);

使用时将三个.m文件置于同一路径内即可

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值