Geotools中实现NC转等值面

概述：

前面的文章有实现IDW插值并生成等值面的，本文在前文基础上实现气象NC数据生成等值面。

效果：

Arcgis预览图

实现后

3、代码

package com.lzugis.netcdf;import com.lzugis.CommonMethod;
import com.lzugis.geotools.EquiSurface;
import ucar.nc2.NetcdfFile;
import ucar.nc2.Variable;
import ucar.nc2.dataset.NetcdfDataset;
import wContour.Contour;
import wContour.Global.Border;
import wContour.Global.PolyLine;
import wContour.Global.Polygon;import java.io.File;
import java.io.IOException;
import java.util.ArrayList;
import java.util.HashMap;
import java.util.List;
import java.util.Map;public class Nc2EquiSurface {private static String rootPath = System.getProperty("user.dir");private EquiSurface equ = new EquiSurface();public Map getNcData(String ncpath) {Map map = new HashMap();NetcdfFile ncfile = null;try {ncfile = NetcdfDataset.open(ncpath);Variable varLon = ncfile.findVariable("lon");Variable varLat = ncfile.findVariable("lat");Variable varPre = ncfile.findVariable("TEM");float[] lon = (float[]) varLon.read().copyToNDJavaArray();float[] lat = (float[]) varLat.read().copyToNDJavaArray();float[][] pre = (float[][]) varPre.read().copyToNDJavaArray();double[] dLon = new double[lon.length], dLat = new double[lon.length];double[][] dPre = new double[pre.length][pre[0].length];for (int i = 0, len = lon.length; i < len; i++) {dLon[i] = Double.parseDouble(String.valueOf(lon[i]));}for (int i = 0, len = lat.length; i < len; i++) {dLat[i] = Double.parseDouble(String.valueOf(lat[i]));}for (int i = 0, len = pre.length; i < len; i++) {float[] _pre = pre[i];for (int j = 0, jlen = _pre.length; j < jlen; j++) {dPre[i][j] = Double.parseDouble(String.valueOf(_pre[j]));}}map.put("lon", dLon);map.put("lat", dLat);map.put("tem", dPre);} catch (IOException e1) {e1.printStackTrace();}return map;}public String nc2EquiSurface(Map ncData, double[] dataInterval) {String geojsonpogylon = "";List<PolyLine> cPolylineList = new ArrayList<PolyLine>();List<Polygon> cPolygonList = new ArrayList<Polygon>();double[][] _gridData = (double[][]) ncData.get("tem");int[][] S1 = new int[_gridData.length][_gridData[0].length];double[] _X = (double[]) ncData.get("lon"), _Y = (double[]) ncData.get("lat");double _undefData = -9999.0;List<Border> _borders = Contour.tracingBorders(_gridData, _X, _Y,S1, _undefData);int nc = dataInterval.length;cPolylineList = Contour.tracingContourLines(_gridData, _X, _Y, nc,dataInterval, _undefData, _borders, S1);// 生成等值线cPolylineList = Contour.smoothLines(cPolylineList);// 平滑cPolygonList = Contour.tracingPolygons(_gridData, cPolylineList,_borders, dataInterval);geojsonpogylon = equ.getPolygonGeoJson(cPolygonList);return geojsonpogylon;}public static void main(String[] args) {Nc2EquiSurface nc2equ = new Nc2EquiSurface();CommonMethod cm = new CommonMethod();long start = System.currentTimeMillis();String ncpath = rootPath + "/data/nc/gdfs1.nc";//获取NC的数据Map map = nc2equ.getNcData(ncpath);//根据NC生成等值面double[] dataInterval = new double[]{0, 5, 10, 15, 20, 25, 30, 35, 40, 45};String strGeojson = nc2equ.nc2EquiSurface(map, dataInterval);String strFile = rootPath + "/out/china_pre.json";cm.append2File(strFile, strGeojson);System.out.println("Total cost:" + (System.currentTimeMillis() - start));}
}

----------------------------------------------------------------------------------------------

如果我的文章对您有帮助，谢谢支持！

技术博客

CSDN：http://blog.csdn.NET/gisshixisheng

博客园：http://www.cnblogs.com/lzugis/

在线教程

http://edu.csdn.Net/course/detail/799

Github

https://github.com/lzugis/

联系方式

q q:1004740957

e-mail:niujp08@qq.com

公众号:lzugis15

Q Q 群:452117357(webgis)

337469080(Android)

Geotools中实现NC转等值面相关推荐

geotools中等值面的生成与OL3中的展示
概述: 本文讲述如何在geotools中IDW插值生成等值面,并根据给定shp进行裁剪,并生成geojson数据,以及Openlayers3中展示. 效果: 插值数据裁剪结果裁剪区域数据实现代码 ...
解决GeoTools中CQL解析中文字段名的问题
GeoTools中CQL无法解析中文字段名的过滤条件,会报异常错误,经过一个下午的努力,终于通过简单有效的方式解决啦 String filterCondition = "temp='&quo ...
Linux中安装nc（netcat）常见问题
Linux中安装nc(netcat)常见问题 Spark Streaming编写wordCount程序时,在Linux集群中需要安装nc,来对程序中使用到的端口进行开放. 在安装nc的时候,常出现的问 ...
GeoJson的生成与解析，JSON解析，Java读写geojson，geotools读取shp文件，Geotools中Geometry对象与GeoJson的相互转换
GeoJson的生成与解析一.wkt格式的geometry转成json格式二.json格式转wkt格式三.json格式的数据进行解析四.Java读写geojson 五.geotools读取sh ...
电路原理图中的“NC“是什么意思？
电路原理图中的"NC"是什么意思? 1.在看电路原理图的时候,电路原理图上有"NC",我查了下是表示"此处不贴任何电子器件",我看了下实际的 ...
Geotools中蜂巢的实现
概述本文讲述如何在geotools中实现蜂巢效果. 效果实现 1.扩展类IntersectionBuilder import com.vividsolutions.jts.geom.Coordin ...
GeoTools中的空间关系（Geometry Relationships）和空间操作（Geometry Operations）
参考原文:GeoTools中的空间关系和空间操作一.GeoTools包含的空间关系和空间操作: 空间关系(Geometry Relationships): 常见的空间关系(Geometry Rela ...
java操作geotiff_关于java使用geotools中的GeoTiffReader读取tif图像时报错
关于java使用geotools中的GeoTiffReader读取tif图像时报错:找不到类 javax.media.jai.PlanarImage 的错误代码放上: import java.io. ...
Linux中的nc测试端口是否开放
nc测试端口是否开放在Linux中有一个级强大的网络工具netcat,在默认情况下面都是没有安装的,现在介绍一下安装过程其实安装很简单一.安装使用 1.只需输入命令yum安装: [root@S ...

Geotools中实现NC转等值面

Geotools中实现NC转等值面相关推荐

最新文章

热门文章