把surfer的.grd文件改写为.txt

What is this script for?
My tutor gave me a .grd file and ask me to deal with that with MATLAB, so I have to rewrite it into another form that I could make use of that afterward.

So here is the code.

filename: readte.m

clear,clc;
%This script is to read grd files and convert it into (lon,lat,value) form
%in .txt file, filename is output.txt.
%.grd files should be from surfer.
%Keith Jin, ZJU, 21/08/02, ver. 1.0
[matrix xmin xmax ymin ymax]=grd_read_v2('out.grd');
file = fopen('output.txt','w');
countx = 0;
fprintf(file,"lon lat te\n");
for x = xmin:0.10034188034188:xmaxcountx = countx+1;county = 0;for y = ymin:0.1:ymaxcounty = county+1;fprintf(file,"%.4f %.4f %f\n",x,y,matrix(county,countx));end %end y
end %end x
fclose('all');

and here’s another one grd read function, from

https://www.mathworks.com/matlabcentral/fileexchange/20880-surfer-grid-import-export?s_tid=FX_rc1_behav

function [matrix xmin xmax ymin ymax]=grd_read_v2(namefile)
% Function to read a GRD file
%                  (from Golden Software Surfer, ASCII format)
%
% [matrix xmin xmax ymin ymax]=grd_read_v2(name of file)
%
% Input:
%      nomarch = name of the file to be read, including ".grd" extension
% Output:
%      matrix =  matrix of the read data
%      xmin xmax ymin ymax = grid limits
%
% Coded by Alberto Avila Armella.
%          UPDATED & IMPROVED BY Jose Maria Garcia-Valdecasas
grdfile=fopen(namefile,'r');    % Open file
code=fgetl(grdfile);            % Reads surfer code 'DSAA'
% Grid dimensions (number of nodes)
aux=str2num(fgetl(grdfile)); nx=aux(1); ny=aux(2);
% X limits
aux=str2num(fgetl(grdfile)); xmin=aux(1); xmax=aux(2);
% Y limits
aux=str2num(fgetl(grdfile)); ymin=aux(1); ymax=aux(2);
% Z limits
aux=str2num(fgetl(grdfile)); zmin=aux(1); aux(2);
% Read matrix
[matrix,count] = fscanf(grdfile, '%f', [nx,ny]);
matrix=matrix';   % Trasposes matrix
fclose(grdfile);

把surfer的.grd文件改写为.txt相关推荐

VS配合VSS时，编译报错：未能向文件“....csproj.FileListAbsolute.txt”写入命令行,对路径的访问被拒绝。...
最近帮别人调试的时候,出现这个问题.环境:VS2008+VSS2005,多人配合完成一个项目. 在编译的时候报错:未能向文件"....csproj.FileListAbsolute.txt& ...
从指定文件夹里COPY指定的一批文件列表（TXT文件）
现实总是层出不穷的,像分形几何一样,一颗树永远在不断的相似变化 .... 现在的问题不是COPY,而是测试将一个包含近25万个文件的目录里COPY出上千个文件,放在另一个目录下作测试. 看看是不是因为 ...
python关闭读写的所有的文件-Python读写txt文本文件的操作方法全解析
一.文件的打开和创建 >>> f = open('/tmp/test.txt') >>> f.read() 'hello python! hello world! ...
python读取整个txt文件-Python读写txt文本文件
一.文件的打开和创建 1 2 3 4 5 >>> f= open('/tmp/test.txt') >>> f.read() 'hello python! hell ...
Android 读取assets文件下的txt文件
ndroid 读取assets文件下的txt文件,解决了读取txt文件的乱码问题: ? 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 ...
[20161107]关于数据文件位图区.txt
[20161107]关于数据文件位图区.txt --上个星期写一篇 [20161104]File Space Bitmap Block损坏能修复吗?.txt 链接:http://blog.itpub. ...
JAVA编程TXT文件_java读写txt文件的方法
java读写txt文件的方法发布时间:2020-06-26 15:54:02 来源:亿速云阅读:111 作者:Leah 本篇文章为大家展示了java读写txt文件的方法,代码简明扼要并且容易理解, ...
ML之FE：基于FE特征工程对RentListingInquries数据集进行预处理并导出为三种格式文件(csv格式/txt格式/libsvm稀疏txt格式)
ML之FE:基于FE特征工程对RentListingInquries数据集进行预处理并导出为三种格式文件(csv格式/txt格式/libsvm稀疏txt格式) 目录输出结果设计思路核心代码输出 ...
用python修改文件内容修改txt内容的3种方法
用python修改文件内容修改txt内容的3种方法方法一.修改原文件方式 def updateFile(file,old_str,new_str):"""替换文件中的字 ...
[20161031]rman备份与数据文件OS块.txt
[20161031]rman备份与数据文件OS块.txt --每个数据文件都有一个OS块,位于数据文件的第1块(也是0块).通过bbed无法访问: BBED> set dba 7,0 BBED- ...

把surfer的.grd文件改写为.txt

把surfer的.grd文件改写为.txt

把surfer的.grd文件改写为.txt相关推荐

最新文章

热门文章