R语言二进制grd文件读取合并导出

news/2024/12/22 1:11:30/

数据格式:1* 31* 10,共17个文件,存在D:\cyh\zpmi\路径下

setwd("D:/cyh/zpmi/")
path <- "D:/cyh/zpmi"
fileNames <- dir(path)
filePath <- sapply(fileNames, function(x){paste(path,x,sep='/')})
wind<-array(0,dim = c(1,31,10,17))
for (i in 1:17) {filename<-file(filePath[i],"rb")data<-readBin(filename,numeric(),1*31*10,size = 4)wind[,,,i]<-array(data=data,dim = c(1,31,10))close(filename)
}
write.filename = file("D:/cyh/zpmi/cyh.grd", "wb")
# 写为二进制文件
writeBin( as.vector(wind),write.filename,size = 4)
# 关闭文件
close(write.filename)
#######读取合并后文件#######
filename<-file("D:\\cyh\\zpmi\\cyh.grd","rb")
data<-readBin(filename,numeric(),1*31*10*17,size = 4)
bestdata<-array(data=data,dim = c(1,31,10,17))
close(filename)

未导入路径图片
导入后路径图片


http://www.ppmy.cn/news/359361.html

相关文章

linux下gmt4绘制站点分布图,GMT4 脚本风格指南

这篇博文会介绍写 GMT 脚本时需要遵循的一些风格与习惯。本文的目的是&#xff0c;希望遵循本风格指南的 GMT 脚本能够更易读、易改、更健壮、可移植性更高。 使用脚本来执行 GMT 命令 GMT 遵循了 UNIX 的设计思想&#xff0c;将不同的功能分别放在不同的命令中&#xff0c;因而…

Matlab读取grd文件,转nc文件存储

1.grd文件读取&#xff08;多维信息提取&#xff09; 例子1&#xff1a;“ clc;clear all; fidfopen(‘G:\work\4\pre\rain6.grd’,‘rb’); [ssh,count]fread(fid,inf,‘float32’); %inf表示读取全部数据float32为浮点型 fclose(fid); %size(a) sshreshape(ssh,27,17,1,3,366…

c语言读grd文件,基于GDAL库,读取.grd文件(以海洋地形数据为例)C++版

技术背景 海洋地形数据主要是通过美国全球地形起伏数据(GMT)获得&#xff0c;数据格式为grd(GSBG)二进制数据&#xff0c;打开软件通过是Surfer软件&#xff0c;surfer软件可进行数据的编辑处理&#xff0c;以及进一步的可视化表达等功能操作&#xff1b;由于Surfer软件不支持二…

Grd文件格式说明

绘图软件Golden Surfer采用的网格文件(grd)主要有Surfer 6 Text、Surfer 6 Binary以及Surfer 7三种存储格式&#xff0c;目前常用前两种数据格式。Surfer 6 Text格式以Ascii码存储&#xff0c;可直接用文本编辑器查看&#xff0c;提取及修改数据均较为方便&#xff1b;另外两种…

基于GDAL库读取GRD数据文件的C++代码

项目需要解析etop数据&#xff0c;都是grd格式文件&#xff0c;可以用surfer软件打开&#xff0c;但是并不适合开发&#xff0c;以及对数据的解析 这里基于gdal库提供的接口对grd文件进行读取和解析&#xff0c;代码如下 #include <iostream> #include <gdal_priv.h…

GDAL 读取和保存 Grd 文件

文章目录 GDAL 读取和保存 Grd 文件Sufer 6 Binary 文件格式读取 Grd 文件Grd 文件和 GDAL tiff 文件的区别手动实现读取 Grd 文件使用GDAL直接读取Grd文件 保存为 grd 文件测试测试结果测试代码 GDAL 读取和保存 Grd 文件 绘图软件 Golden Surfer 的网格文件&#xff08;Grd&…

将Surfer grd格网文件转成ArcGIS asc格式

今天有同事问&#xff1a;将Surfer grd格网文件转成ArcGIS的栅格文件 博主记得Surfer就有这个功能&#xff0c;但是使用Surfer20的时候&#xff0c;怎么也找不到那个功能。 后来&#xff0c;找遍了grid相关的工具栏&#xff0c;终于找到了&#xff0c;记录一下&#xff0c;免…

python读取grd数据_一个关于grads显示grd格式数据的问题。

是合成差值场的t检验程序。 program main implicit none integer i,j,t,irec integer,parameter :: x=145 integer,parameter :: y=73 integer,parameter :: n=58 !总体样本数 real,parameter :: n1=10.0 !偏早年份 real,parameter :: n2=10.0 !偏晚年份 real :: dat(…