Netcdf对数据进行裁剪

对三维数据进行裁剪

List<Range> tyxRanges = new ArrayList<>();
tyxRanges.add(new Range(null,0,12));
tyxRanges.add(new Range(null,0,12));
tyxRanges.add(new Range(null,0,12));
Array varData = varObject.read(tyxRanges);
ncWrite.write(varVar, varData);

只显示了左下角的区域。

增加对时间变量的裁剪

public String varExport(@ModelAttribute FileSaveVars params ) throws IOException, InvalidRangeException, ParseException {//      0° 设置时间String timeStart = "2018-01-01 00:30:00";String timeEnd = "2018-01-01 00:40:00";
//      1° 设置原始文件路径+转存文件路径String filePath = "/Users/caowei/workspace/jszone/大气系统/ncdemo/src/main/resources/static/MERRA2_400.tavg1_2d_flx_Nx.20180101.SUB.nc";String savePath = "/Users/caowei/workspace/test.nc";
//      2° 读取文件NetcdfFile ncFile = NetcdfFile.open(filePath, null);
//      3° 设置需要转存的变量String varName = "QLML";
//      4° 设置时间变量Variable timeObject = ncFile.findVariable("time");String timeAt = timeObject.getUnitsString();String regex = "\\d{4}-\\d{2}-\\d{2} \\d{2}:\\d{2}:\\d{2}";Pattern p = Pattern.compile(regex);Matcher matcher = p.matcher(timeAt);if (matcher.find()) {timeAt = matcher.group(0);}
//      4-1° 获取时间差int timeDiff1 = DateHelp.dateDifference(timeAt,timeStart);int timeDiff2 = DateHelp.dateDifference(timeAt,timeEnd);
//      4-2° 获取时间步长Array timeArray =  timeObject.read();int timeStep = (int)(timeArray.getFloat(1) - timeArray.getFloat(0)) ;int startIndex = timeDiff1/timeStep;int endIndex = timeDiff2/timeStep;
//      5° 获取需要转存的变量对象Variable varObject = ncFile.findVariable(varName);
//      6° 构造输出保存文件NetcdfFileWriter ncWrite = NetcdfFileWriter.createNew(NetcdfFileWriter.Version.netcdf4, savePath);
//      6° 设置文件的一些参数
//      7°-1纬度for(Dimension dimension : varObject.getDimensions()){System.out.println(dimension.getShortName());System.out.println(dimension.getLength());ncWrite.addDimension(null, dimension.getShortName(), dimension.getLength());}
//      7°-2变量Variable varVar = ncWrite.addVariable(null, varObject.getShortName(), varObject.getDataType(), varObject.getDimensionsString());System.out.println(varVar.getDimensions());
//      7°-3设置属性List<Attribute> AttrList = varObject.getAttributes();for (Attribute att : AttrList) {varVar.addAttribute(att);}
//      8°创建ncWrite.create();
//      9°将数据写入List<Range> tyxRanges = new ArrayList<>();tyxRanges.add(new Range(null,startIndex,endIndex));tyxRanges.add(new Range(null,0,120));tyxRanges.add(new Range(null,0,120));Array varData = varObject.read(tyxRanges);ncWrite.write(varVar, varData);
//      10°关闭保存ncWrite.close();return  "OK";
}

Netcdf对数据进行裁剪相关推荐

使用SHP数据批量裁剪TIF图像并转为JPG格式
问题描述已有数据:1张大范围的遥感影像和1张裁剪为同一大小的shp数据使用工具:ArcPy 实现目标: 1.shp数据按照属性分块输出 2.使用shp数据批量裁剪遥感影像,将影像裁剪为小块的tif ...
如何从省级地图数据中裁剪出州级地图数据
如何从省级地图数据中裁剪出州级地图数据开发工具与关键技术:SuperMap iDesktop 9D 作者:蒙伟志撰写时间: 先从省级地图裁剪出你所需要的地图,比如说我裁剪的是从化区,如图下粉红色区 ...
使用shp数据批量裁剪栅格数据并统计均值
批量实现利用shp文件,裁剪相应的多个栅格数据,并统计裁剪后栅格的均值至TXT文件,亲测可用. #使用shp数据批量裁剪栅格数据并统计均值 # -*- coding: cp936 -*- #使用shp ...
gis 数据框裁剪_BIM+GIS的八大挑战！大挑战，见未来
BIM与GIS的跨界融合,使微观领域的BIM信息与宏观领域的GIS信息实现交换和互操作,提升了BIM应用深度,将BIM的应用从单体延伸到建筑群甚至城市级,为GIS行业发展带来了新的契机,同时也带来了一 ...
gis 数据框裁剪_ArcGis入门学习教程带练习数据
百度网盘链接:https://pan.baidu.com/s/18iTEYMygozKZT2vuzLY8Dg 提取码:pwp8 目录实验一.使用 ARCMAP浏览地理数据 1 一.实验目的 ...
【ArcGIS遇上Python】窗体版Python批量处理地理数据--栅格裁剪
一.说明之前写过用Python代码块处理地理数据的例子,如裁剪.投影.格式转换等,曾想过用C#制作个窗体,嵌入Python代码,今天无意发现Python也能用tKinter构建简单的窗体,眼前一亮, ...
目标检测矩形框与polygon数据增加--裁剪，拓展，旋转
1.裁剪 import torch from torchvision import transforms import cv2 import numpy as np import types from ...
【ADNI】对 nii 3D 数据进行裁剪（Matlab）待完善...
crop_ADNI_GM.m function num_list = crop_ADNI_GM(filelistname, save_path, label) %%%% %% 脚本说明 % 1)对AD ...
gis 数据框裁剪_【GIS实操】如何用GIS提取“三调”数据？
做国土空间规划时,最重要的资料之一就是三调GIS数据.当资料到手,会发现一个问题:我只需要做一个镇/村的规划,但是整个市/镇数据在手,怎样才能将我需要的内容提取出来呢?这就需要用到GIS的剪裁工具!剪 ...

Netcdf对数据进行裁剪

Netcdf对数据进行裁剪相关推荐

最新文章

热门文章