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

当数据的Z值为常数时，surfer不能进行网格化而直接得到全为均值的.grd数据（即得到如下数据）。

使用以下函数可生成上述形式的数据。
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文件置于同一路径内即可