netcdf3实现三四维数组数据中多个要素写入生成nc文件

最新推荐文章于 2024-08-08 22:16:05 发布
最新推荐文章于 2024-08-08 22:16:05 发布
阅读量929 收藏
点赞数 1
文章标签： java
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/m0_46682198/article/details/116786499
版权

需要一个jar包
E:\ServerFiles\Hebei\workspace\Demo1\WebContent\WEB-INF\lib\netcdfAll-4.6.10.jar


public static void createNCfile(String[] str1,String filepath1,int width,int height,float lonWest,float latSouth,float gridInterval,List<float[][]> Dewpoint_temperature ,List<float[][]> Pressure ,List<float[][]> Relative ,List<float[][]> Temperature, List<float[][]> Total_precipitation, List<float[][]> Wind ,List<float[][]> Wind_speed, List<float[][]> u_component_of_wind_height, List<float[][]> u_component_of_wind_gust ,List<float[][]> v_component_of_wind ,List<float[][]> v_component_of_wind_gust){
        //创建netcdf3写入文件对象
        NetcdfFileWriter fileWriter;
        try {
            fileWriter = NetcdfFileWriter.createNew(Version.netcdf3, filepath1);
            //分别创建经度和纬度方向的维度，lon的宽度为width,lat高度为height,并命名为lon和lat
            Dimension xdim=fileWriter.addDimension(null, "lon", width);
            Dimension ydim=fileWriter.addDimension(null, "lat", height);
            Dimension time=fileWriter.addDimension(null, "time", 1);
            Dimension grounp=fileWriter.addDimension(null, "height_above_ground", 1);
            Dimension grounp1=fileWriter.addDimension(null, "height_above_ground1", 1);
            //通过一维组合成二维array list
            ArrayList<Dimension> dims_Dewpoint_temperature=new ArrayList<Dimension>();
            ArrayList<Dimension> dims_Pressure=new ArrayList<Dimension>();
            ArrayList<Dimension> dims_Relative=new ArrayList<Dimension>();
            ArrayList<Dimension> dims_Temperature=new ArrayList<Dimension>();
            ArrayList<Dimension> dims_Total_precipitation=new ArrayList<Dimension>();
            ArrayList<Dimension> dims_Wind=new ArrayList<Dimension>();
            ArrayList<Dimension> dims_Wind_speed=new ArrayList<Dimension>();
            ArrayList<Dimension> dims_u_component_of_wind_height=new ArrayList<Dimension>();
            ArrayList<Dimension> dims_u_component_of_wind_gust=new ArrayList<Dimension>();
            ArrayList<Dimension> dims_v_component_of_wind=new ArrayList<Dimension>();
            ArrayList<Dimension> dims_v_component_of_wind_gust=new ArrayList<Dimension>();
            //注意这里add的顺序很重要,一定要和数据对应，如果把list顺序调换，那么二维的维度也会不同，相当于数学中的坐标系，第一个是x坐标，第二个是y坐标，x宽度为width,y高度为height
            dims_Dewpoint_temperature.add(time);
            dims_Dewpoint_temperature.add(grounp1);
            dims_Dewpoint_temperature.add(xdim);
            dims_Dewpoint_temperature.add(ydim);      
            dims_Pressure.add(time);
            dims_Pressure.add(xdim);
            dims_Pressure.add(ydim);
            dims_Relative.add(time);
            dims_Relative.add(grounp);
            dims_Relative.add(xdim);
            dims_Relative.add(ydim);
            dims_Temperature.add(time);
            dims_Temperature.add(grounp);
            dims_Temperature.add(xdim);
            dims_Temperature.add(ydim);
            dims_Total_precipitation.add(time);
            dims_Total_precipitation.add(xdim);
            dims_Total_precipitation.add(ydim);
            dims_Wind.add(time);
            dims_Wind.add(grounp);
            dims_Wind.add(xdim);
            dims_Wind.add(ydim);
            dims_Wind_speed.add(time);
            dims_Wind_speed.add(grounp);
            dims_Wind_speed.add(xdim);
            dims_Wind_speed.add(ydim);
            dims_u_component_of_wind_height.add(time);
            dims_u_component_of_wind_height.add(grounp);
            dims_u_component_of_wind_height.add(xdim);
            dims_u_component_of_wind_height.add(ydim);
            dims_u_component_of_wind_gust.add(time);
            dims_u_component_of_wind_gust.add(grounp);
            dims_u_component_of_wind_gust.add(xdim);
            dims_u_component_of_wind_gust.add(ydim);
            dims_v_component_of_wind.add(time);
            dims_v_component_of_wind.add(grounp);
            dims_v_component_of_wind.add(xdim);
            dims_v_component_of_wind.add(ydim);
            dims_v_component_of_wind_gust.add(time);
            dims_v_component_of_wind_gust.add(grounp);
            dims_v_component_of_wind_gust.add(xdim);
            dims_v_component_of_wind_gust.add(ydim);
            //创建名称为lon的变量,类型为folat，对应的维度为lon,对应Dimension里面定义的名称为"lon"的那个
            Variable vx=fileWriter.addVariable(null, "lon", DataType.FLOAT, "lon");
            //给vx变量加上数据说明，全程long_name
            fileWriter.addVariableAttribute(vx, new Attribute("long_name","longitude" ));
            //给vx变量加上units，对应的Degree_east，如果不加units为Degree_east(东经)，则不是地理坐标系，这个很重要
            fileWriter.addVariableAttribute(vx, new Attribute("units","Degrees_east" ));

            Variable vy=fileWriter.addVariable(null, "lat", DataType.FLOAT, "lat"); 
            fileWriter.addVariableAttribute(vy, new Attribute("long_name","latitude" ));
            //给vy变量加上units，对应的Degree_north，如果不加units为Degree_north(北纬)，则不是地理坐标系，这个很重要
            fileWriter.addVariableAttribute(vy, new Attribute("units","Degrees_north" ));
            
            Variable vtime=fileWriter.addVariable(null, "time", DataType.DOUBLE, "time");
            fileWriter.addVariableAttribute(vtime, new Attribute("long_name","time"));
            fileWriter.addVariableAttribute(vtime, new Attribute("units","Hour since 2021-04-26T04:00:00Z"));
            
            Variable vgrounp=fileWriter.addVariable(null, "height_above_ground", DataType.FLOAT, "height_above_ground");
            fileWriter.addVariableAttribute(vgrounp, new Attribute("long_name","height_above_ground"));
            fileWriter.addVariableAttribute(vgrounp, new Attribute("units","m"));
            
            Variable vgrounp1=fileWriter.addVariable(null, "height_above_ground1", DataType.FLOAT, "height_above_ground1");
            fileWriter.addVariableAttribute(vgrounp1, new Attribute("long_name","height_above_ground1"));
            fileWriter.addVariableAttribute(vgrounp1, new Attribute("units","m"));
            
            //创建变量名称为var的变量，对应的维度为dims,该dims为上面定义的一个ArrayList<Diminsion>，该list包含2个维度，经度和维度
            Variable v_dims_Dewpoint_temperature=fileWriter.addVariable(null, str1[0], DataType.FLOAT, dims_Dewpoint_temperature);
            fileWriter.addVariableAttribute(v_dims_Dewpoint_temperature, new Attribute("long_name","露点温度" ));
            fileWriter.addVariableAttribute(v_dims_Dewpoint_temperature, new Attribute("units","K" ));
            
            Variable v_dims_Pressure=fileWriter.addVariable(null, str1[1], DataType.FLOAT, dims_Pressure);
            fileWriter.addVariableAttribute(v_dims_Pressure, new Attribute("long_name","压力" ));
            fileWriter.addVariableAttribute(v_dims_Pressure, new Attribute("units","Pa" ));
            
            Variable v_dims_Relative=fileWriter.addVariable(null, str1[2], DataType.FLOAT, dims_Relative);
            fileWriter.addVariableAttribute(v_dims_Relative, new Attribute("long_name","相对湿度" ));
            fileWriter.addVariableAttribute(v_dims_Relative, new Attribute("units","%" ));
            
            Variable v_dims_Temperature=fileWriter.addVariable(null, str1[3], DataType.FLOAT, dims_Temperature);
            fileWriter.addVariableAttribute(v_dims_Temperature, new Attribute("long_name","地上温度" ));
            fileWriter.addVariableAttribute(v_dims_Temperature, new Attribute("units","K" ));
            
            Variable v_dims_Total_precipitation=fileWriter.addVariable(null, str1[4], DataType.FLOAT, dims_Total_precipitation);
            fileWriter.addVariableAttribute(v_dims_Total_precipitation, new Attribute("long_name","总降水量" ));
            fileWriter.addVariableAttribute(v_dims_Total_precipitation, new Attribute("units","kg.m-2" ));
            
            Variable v_dims_Wind=fileWriter.addVariable(null, str1[5], DataType.FLOAT, dims_Wind);
            fileWriter.addVariableAttribute(v_dims_Wind, new Attribute("long_name","风向" ));
            fileWriter.addVariableAttribute(v_dims_Wind, new Attribute("units","degree_true" ));
            
            Variable v_dims_Wind_speed=fileWriter.addVariable(null, str1[6], DataType.FLOAT, dims_Wind_speed);
            fileWriter.addVariableAttribute(v_dims_Wind_speed, new Attribute("long_name","风速" ));
            fileWriter.addVariableAttribute(v_dims_Wind_speed, new Attribute("units","m/s" ));
            
            Variable v_dims_u_component_of_wind_height=fileWriter.addVariable(null, str1[7], DataType.FLOAT, dims_u_component_of_wind_height);
            fileWriter.addVariableAttribute(v_dims_u_component_of_wind_height, new Attribute("long_name","地上风高度的u分量" ));
            fileWriter.addVariableAttribute(v_dims_u_component_of_wind_height, new Attribute("units","m/s" ));
            
            Variable v_dims_u_component_of_wind_gust=fileWriter.addVariable(null, str1[8], DataType.FLOAT, dims_u_component_of_wind_gust);
            fileWriter.addVariableAttribute(v_dims_u_component_of_wind_gust, new Attribute("long_name","风、阵风、离地高度的u分量" ));
            fileWriter.addVariableAttribute(v_dims_u_component_of_wind_gust, new Attribute("units","m/s" ));
            
            Variable v_dims_v_component_of_wind=fileWriter.addVariable(null, str1[9], DataType.FLOAT, dims_v_component_of_wind);
            fileWriter.addVariableAttribute(v_dims_v_component_of_wind, new Attribute("long_name","地上风力高度的v分量" ));
            fileWriter.addVariableAttribute(v_dims_v_component_of_wind, new Attribute("units","m/s" ));
            
            Variable v_dims_v_component_of_wind_gust=fileWriter.addVariable(null, str1[10], DataType.FLOAT, dims_v_component_of_wind_gust);
            fileWriter.addVariableAttribute(v_dims_v_component_of_wind_gust, new Attribute("long_name","风、阵风、地上高度的v分量" ));
            fileWriter.addVariableAttribute(v_dims_v_component_of_wind_gust, new Attribute("units","m/s" ));
            
            //使用create()来创建该文件，只有create之后才能在文件夹中显示该文件，并且写入数据
            fileWriter.create();
            //创建x和y方向上的变量。D1    代表一维
            ArrayFloat xvalues=new ArrayFloat.D1(width);
            ArrayFloat yvalues=new ArrayFloat.D1(height);
            ArrayFloat timevalues=new ArrayFloat.D1(1);
            ArrayFloat grounpvalues = new ArrayFloat.D1(1);
            ArrayFloat grounp1values = new ArrayFloat.D1(1);
            for (int i = 0; i < width; i++) {
                xvalues.setFloat(i, (float) (lonWest+gridInterval*i));
            }
            for (int j = 0; j < height; j++) {
                yvalues.setFloat(j, (float) (latSouth+gridInterval*j));
            }
            timevalues.setInt(0, 0);
            grounpvalues.setInt(0, 0);
            grounp1values.setInt(0, 0);
            //创建网格上的变量，D3代表3维，用(1,width, height)来描述他的定义域，这个顺序必须和ArrayList<Dimension> 加载的顺序一一对应，否则会出错
            ArrayFloat values_Dewpoint_temperature=new ArrayFloat.D4(1,1,width, height);
            ArrayFloat values_Pressure=new ArrayFloat.D3(1,width, height);
            ArrayFloat values_Relative=new ArrayFloat.D4(1,1,width, height);
            ArrayFloat values_Temperature=new ArrayFloat.D4(1,1,width, height);
            ArrayFloat values_Total_precipitation=new ArrayFloat.D3(1,width, height);
            ArrayFloat values_Wind=new ArrayFloat.D4(1,1,width, height);
            ArrayFloat values_Wind_speed=new ArrayFloat.D4(1,1,width, height);
            ArrayFloat values_u_component_of_wind_height=new ArrayFloat.D4(1,1,width, height);
            ArrayFloat values_u_component_of_wind_gust=new ArrayFloat.D4(1,1,width, height);
            ArrayFloat values_v_component_of_wind=new ArrayFloat.D4(1,1,width, height);
            ArrayFloat values_v_component_of_wind_gust=new ArrayFloat.D4(1,1,width, height);
            Index index_Dewpoint_temperature=values_Dewpoint_temperature.getIndex();
            Index index_Pressure=values_Pressure.getIndex();
            Index index_Relative=values_Relative.getIndex();
            Index index_Temperature=values_Temperature.getIndex();
            Index index_Total_precipitation=values_Total_precipitation.getIndex();
            Index index_Wind=values_Wind.getIndex();
            Index index_Wind_speed=values_Wind_speed.getIndex();
            Index index_u_component_of_wind_height=values_u_component_of_wind_height.getIndex();
            Index index_u_component_of_wind_gust=values_u_component_of_wind_gust.getIndex();
            Index index_v_component_of_wind=values_v_component_of_wind.getIndex();
            Index index_v_component_of_wind_gust=values_v_component_of_wind_gust.getIndex();
            for (int n = 0; n < 1; n++) {
            	for (int k = 0; k < 1; k++) {
                	for (int i = 0; i < width; i++) {
                        for (int j = 0; j < height; j++) {
                            //将二维数组[i][j]处的数据写入网格的i,j处
                        	values_Dewpoint_temperature.setFloat(index_Dewpoint_temperature.set(n,k,i, j), Dewpoint_temperature.get(k)[i][j]); 
                        	values_Relative.setFloat(index_Relative.set(n,k,i,j), Pressure.get(k)[i][j]);
                        	values_Temperature.setFloat(index_Temperature.set(n,k,i,j), Temperature.get(k)[i][j]);
                        	values_Wind.setFloat(index_Wind.set(n,k,i,j), Wind.get(k)[i][j]);
                        	values_Wind_speed.setFloat(index_Wind_speed.set(n,k,i,j), Wind_speed.get(k)[i][j]);
                        	values_u_component_of_wind_height.setFloat(index_u_component_of_wind_height.set(n,k,i,j), u_component_of_wind_height.get(k)[i][j]);
                        	values_u_component_of_wind_gust.setFloat(index_u_component_of_wind_gust.set(n,k,i,j), u_component_of_wind_gust.get(k)[i][j]);
                        	values_v_component_of_wind.setFloat(index_v_component_of_wind.set(n,k,i,j), v_component_of_wind.get(k)[i][j]);
                        	values_v_component_of_wind_gust.setFloat(index_v_component_of_wind_gust.set(n,k,i,j), v_component_of_wind_gust.get(k)[i][j]);                       	
                        }
                    }
    			}
			}
            for (int k = 0; k < 1; k++) {
            	for (int i = 0; i < width; i++) {
                    for (int j = 0; j < height; j++) {
                        //将二维数组[i][j]处的数据写入网格的i,j处
                    	values_Pressure.setFloat(index_Pressure.set(k,i,j), Pressure.get(k)[i][j]);
                    	values_Total_precipitation.setFloat(index_Total_precipitation.set(k,i,j), Total_precipitation.get(k)[i][j]);                   }
                }

			}
            fileWriter.write(vx, xvalues);
            fileWriter.write(vy, yvalues);
            fileWriter.write(vtime, timevalues);
            fileWriter.write(vgrounp, grounpvalues);
            fileWriter.write(vgrounp1, grounp1values);
            fileWriter.write(v_dims_Dewpoint_temperature, values_Dewpoint_temperature);
            fileWriter.write(v_dims_Pressure, values_Pressure);
            fileWriter.write(v_dims_Relative, values_Relative);
            fileWriter.write(v_dims_Temperature, values_Temperature);
            fileWriter.write(v_dims_Total_precipitation, values_Total_precipitation);
            fileWriter.write(v_dims_Wind,values_Wind);
            fileWriter.write(v_dims_Wind_speed, values_Wind_speed);
            fileWriter.write(v_dims_u_component_of_wind_height, values_u_component_of_wind_height);
            fileWriter.write(v_dims_u_component_of_wind_gust, values_u_component_of_wind_gust);
            fileWriter.write(v_dims_v_component_of_wind, values_v_component_of_wind);
            fileWriter.write(v_dims_v_component_of_wind_gust, values_v_component_of_wind_gust);
            fileWriter.close();
            System.out.println("文件已存入"+filepath1);
        } catch (IOException e) {
            // TODO Auto-generated catch block
            e.printStackTrace();
            System.out.println("文件创建失败，请检查路径");
        } catch (InvalidRangeException e) {
            // TODO Auto-generated catch block
            e.printStackTrace();
            System.out.println("文件写入错误，超出范围");
        }

    }
java入门的小白
关注
Fortran写nc文件 f90 netcdf
ProMath的专栏
06-03 5068
netCDF的教程例子，用Fortran95把数据写到文件中 编译：gfortran  ****.f90 -I/usr/local/include -L/usr/local/lib -lnetcdff -o test1 运行：./test1 源代码： program simple_xy_wr   use netcdf   implicit none   ! 定义数据文件的名称   characte
python 创建netcdf_python写入二维netCDF文件（txt转nc文件）
weixin_39994461的博客
12-03 1985
netCDF文件后缀名是.nc，是一种二进制的网格文件，无法直接打开，下面给出基于python的netCDF文件生成的脚本。输入数据说明：数据为31°N-40°N，130°E-142°E区域表层的地震波速，网格分布大小不均匀。数据存储在txt文件中，分为三列，无表头。部分数据记录见下图，第一列是经度，第二列是维度，第三列是波速扰动值%过程与代码为了将其生成均匀网格的二维nc文件，我们需要1.文件读...
python 读取某文件夹下所有nc数据并写入
go_with_the_wind的博客
04-06 2562
import numpy as np import netCDF4 as nc import os path = '/media/ciro/B4FE-5315/polarRS_Test2data/polarRS_Test2data/CS2/Greenland/' files = os.listdir(path) x_list = np.array([]); y_list = np.array([]); h_list = np.array([]); for file in files: f = nc.
Java解析气象热力图NC文件：完整实现与代码解析
阿水的博客
08-08 5397
在气象和环境科学领域，NC（NetCDF）文件是一种广泛使用的数据格式，通常用于存储多维科学数据。在这篇博客中，我们将介绍如何使用Java来解析NC文件，并将提取的数据转换为JSON格式。本文将逐步解析代码，帮助你理解每一个关键步骤。
探索NetCDF世界的窗口：netcdfjs —— 高效读取与探索NetCDF v3文件的利器
gitblog_00062的博客
06-15 975
???? 探索NetCDF世界的窗口：netcdfjs —— 高效读取与探索NetCDF v3文件的利器 netcdfjsRead and explore NetCDF files项目地址:https://gitcode.com/gh_mirrors/ne/netcdfjs 在数据科学和地理信息系统的领域中，NetCDF（Network Common Data Form）是一种用于存储多维科学数据的标...
Fortran处理nc文件基于netcdf库
ZJ315296441的博客
08-10 4635
Fortran处理nc文件基于netcdf库 fortran90处理nc文件的学习笔记 编译器为 cygwin下的 gfortran !!! 一定要注意使用的netcdf.a和netcdff.a的区别 !!! 检查你指定的lib文件夹中有没有相对应的.a库 以及 include文件夹中有没有相应的头文件 !运行方式: gfortran.exe ReadAndWriteNC.f90 -I /usr/include -L /usr/lib -lnetcdff program ReadAndWriteNC
python-nc文件的创建和读写
weixin_67321393的博客
09-17 1610
【代码】python-nc文件的创建和读写。
Python实现将数据写入netCDF4中的方法示例
01-20
本文实例讲述了Python实现将数据写入netCDF4中的方法。分享给大家供大家参考，具体如下： nc文件为处理气象数据文件。用户可以去https://www.lfd.uci.edu/~gohlke/pythonlibs/ 搜索netCDF4,下载相应平台的whl文件，...
NCdo.rar_NC_nc文件_写入nc文件_根据矩阵创建一个NC文件
09-24
"NCdo.rar"是一个压缩包，包含了一个名为"NCdo.py"的Python脚本，它实现了从"fx_gridpoint.txt"文本文件中读取矩阵数据，并将其写入一个NC（NetCDF）文件的功能。以下是关于这个过程的一些关键知识点： 1. **NetCDF...
PanoplyWin-4.8.7.zip_NetCDF解析_nc文件_nc文件 显示_nc解析_panoply查看数据
09-15
《深入理解NetCDF：借助PanoplyWin进行nc文件解析与数据显示》 NetCDF（Network Common Data Form）是一种用于存储和处理科学数据的开放标准，它支持多维数组和元数据，广泛应用于气候、环境科学等领域。nc文件是...
python 创建netcdf_在Python中创建NetCDF文件时如何指定维度值？
weixin_39941859的博客
12-03 735
我正在创建一个NetCDF4文件，该文件当前有四个变量：1)地表温度(3D阵列-时间、纬度、经度)2)经度(1D-每个像素中心的坐标)3)纬度(1D-每个像素中心的坐标)图像采集时间1900-01:00我当前正在使用以下代码来执行此操作：#==========================WRITE THE NETCDF FILE==========================#newfil...
fortran使用netcdf 读写NC文件
h4x0r_007的专栏
07-15 2万+
先安装VS2010 再安装intel visual fortran2011 1、下载已编译的netcdf 去下载别人已编译好的netcdf库，http://bbs.lasg.ac.cn/bbs/viewthread.php?tid=955 NetCDF3.6.1Beta1.Win32.for.CVF.IVF.rar 为了简单操作，我将NetCDF3.6.1
【Python学习】——技巧分享（持续更新中）
elite0的博客
08-11 500
接上文，只需要对调经纬度即可。
NetCDF(nc)读写与格式转换介绍
weixin_44785184的博客
01-19 6134
本文介绍了NetCDF文件格式，并详细讲解了如何使用Python对NetCDF文件进行读写操作，进而介绍了NetCDF文件的地理参考，最后以两个数据为例讲解了怎么将NetCDF格式的数据转GeoTIFF格式的数据（.nc文件转为.tif文件）。
NC文件读写方式
热门推荐
痛苦+反省总结=进步
03-17 3万+
概述NetCDF全称为network Common Data Format，中文译法为“网络通用数据格式”，对程序员来说，它和zip、jpeg、bmp文件格式类似，都是一种文件格式的标准。netcdf文件开始的目的是用于存储气象科学中的数据，现在已经成为许多数据采集软件的生成文件的格式。 从数学上来说，netcdf存储的数据就是一个多自变量的单值函数。用公式来说就是f(x,y,z,…)=value
NetCDF(.nc) 文件读取报错：NetCDF: Not a valid ID 的解决
weixin_41671327的博客
01-08 6957
众所周知，NetCDF是一种很容易操作的数据集文件，常应用于气象学、图形学等学科中作为数据的存储。 今日，笔者遇到了一个读取该文件时的轻微报错（Runtime Error）：NetCDF: Not a valid ID 在网上查找解决办法，StackOverflow上一哥们也碰到了类似的报错，最后解决办法是：“Nc 的 操作放在了错误的位置上”（没错就这么模糊） 我最终的解决办法是：同一个N...
JavaScript 自增和自减的区别
木木木森林
10-14 513
自增和自减 自增 - 通过自增可以使变量在自身的基础上增加1 - a++; 使a自增1 - 对一个变量自增以后，原变量的值会立即自增1 - 自增分为两种：后++(a++) 和前++(++a) 无论是a++ 还是++a，都会立即使原变量的值自增1 不同的是a++ 和 ++a的值不同 a++的值等于
输出nc数据_NetCDF（NC）数据的使用、转换和分析
weixin_29475313的博客
01-13 5452
最近项目中需要处理和分析NC数据，所以我查了一下，百度百科的解释是：NetCDF(network Common Data Form)网络通用数据格式是由美国大学大气研究协会（University Corporation for Atmospheric Research，UCAR)的Unidata项目科学家针对科学数据的特点开发的，是一种面向数组型并适于网络共享的数据的描述和编码标准。目...
使用IDL进行netcdf文件读取，出图保存，转geotiff格式
weixin_42183989的博客
11-06 1860
Author:HD 上海某高校遥感专业 工科男 分享一些地理数据处理和数据挖掘的示例 相信积微者速成，相信分享的力量 CSDN@ 这是一个老师布置的小任务。中间包括nc文件的读取，出图样式的调整，色带的自定义花了很长时间去学习，查资料和实验。最终，效果还算满意。希望大家能少走一些弯路，这只是个简单的小示例，更详细的函数参数以及细微调整需要大家自己去动手。 初次写技术博客没有什么经验，如有不足，还请大家多多包涵，欢迎大家私信交流 ^ _ ^ !!!一起进步~ 这里也感谢要感谢 @ 傻灰点这里 @环境遥.
matlab中的多个数组变量如何输出到一个nc文件中
最新发布
09-15
在MATLAB中，将多个数组变量保存到一个.NETCDF (Network Common Data Form) 文件，也称为nc文件，通常需要使用`netcdf`工具箱。首先确保已经安装了这个工具箱。以下是基本步骤： 1. 导入所需库： ```matlab if ~isToolboxInstalled('netcdf') error('NetCDF toolbox not installed. Please install it first.'); end ``` 2. 创建或打开nc文件： ```matlab filename = 'my_data.nc'; % 指定输出文件名 ncid = netcdf(filename, 'w'); % 'w'表示创建新文件或覆盖已有文件 ``` 3. 添加维度（如果有必要）： ```matlab % 假设你有长度分别为m、n、p的三个数组 [m, n] = size(array1); p = size(array2, 3); % 如果array2是三维的 % 定义维度 dim1 = netcdf.addDimension(ncid, 'time', m); dim2 = netcdf.addDimension(ncid, 'space', n); dim3 = netcdf.addDimension(ncid, 'channel', p); ``` 4. 创建变量并设置其属性（包括数据类型和维度）： ```matlab % 对于每个数组，添加变量 var1 = netcdf.addVariable(ncid, 'variable1', 'double', [dim1, dim2]); var2 = netcdf.addVariable(ncid, 'variable2', 'double', [dim1, dim2, dim3]); ``` 5. 将数据写入nc文件： ```matlab % 写入array1的数据 netcdf.putVariable(ncid, var1, array1); % 写入array2的数据，注意调整第三个维度与dim3对应 netcdf.putVariable(ncid, var2, permute(array2, [dim1, dim2, dim3])); % permute可能取决于array2的原始形状 ``` 6. 关闭nc文件： ```matlab netcdf.close(ncid); ```
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值