这些年做项目的时候,碰到等值面基本都是arcgis server来支撑的,通过构建GP服务,一般的都能满足做等值线和等值面的需求。可是突然有一天,我发现如果没有arcgis server 的话,我既然没法生成等值面等值线了。况且,还有许多别的要求:

  1. 没有arcgis server支持,arcgis server毕竟是很大一笔开销,个人基本无法负担;
  2. 跨平台,有的服务器是linux,有的是windows,看来,只能是java的类库了;
  3. 免费,生成等值线和等值面的程序不能有费用;
  4. 输出多种格式,既能输出图片格式,又能输出矢量格式,尤其是 kml之类的矢量数据;

     在气象家园http://bbs.06climate.com/找到一个很不错的解决方案,MeteoInfo,作者是一个很厉害的人,具体的不说了,大家可以去看他的网站和博客,下面我说说如何使用他的类库来进行等值面的生成。如果大家有很好的类似的库,欢迎推荐给我,谢谢。

1:引用他的项目中的lib目录中的jar包,生成窗口对象:

  1 MapView mapView = new MapView();
  2 mapView.setBackground(new Color(255, 255, 255, 0));
  3 mapView.setBounds(0, 0, 200, 200);
  4 VectorLayer clipLayer = MapDataManage.readMapFile_ShapeFile(“d:/chengdu.shp”); //剪切图层,就是生成等值面的形状范围,
    //使用shp最方便了,这儿有一个坑,就是你是用的shp可以通过他自带的软件加载显示才行,否则程序会异常,至于为啥有的shp不能加载,我也没搞清楚
    //反正这个坑让我趟过去了。

2:添加站点信息,用于插值

  1 StationData stationData = new StationData();
  2 //
  3 for(int i=0;i<10;i++)
  4 {
  5 stationData.addData("st"+i, 114+0.1*i, 35+0.1*i, i*10); //站点名称,经度,维度,值
  6 }

3:设定插值参数

  1 GridDataSetting gridDataSetting = new GridDataSetting();
  2 gridDataSetting.dataExtent = clipLayer.getExtent();
  3 stationData.projInfo = clipLayer.getProjInfo();
  4 gridDataSetting.xNum = contourconfig.getGridx();// 格点点数
  5 gridDataSetting.yNum = contourconfig.getGridy();// 格点点数
  6 
  7 InterpolationSetting interSet = new InterpolationSetting();
  8 
  9 interSet.setGridDataSetting(gridDataSetting);
 10 interSet.setInterpolationMethod(InterpolationMethods.IDW_Radius);
 11 interSet.setRadius(5);
 12 interSet.setMinPointNum(1);
 13 GridData gridData = stationData.interpolateData(interSet);
 14 
 15 LegendScheme legendScheme =LegendManage.createLegendSchemeFromGridData(gridData, LegendType.UniqueValue,ShapeTypes.Polygon);
4:插值生成图层:
  1 VectorLayer contourLayer = DrawMeteoData.createShadedLayer(gridData, legendScheme, "ContourLayer", "Data",
  2 					true);
  3 VectorLayer lastLayer = contourLayer.clip(clipLayer);
  4 mapView.addLayer(lastLayer);

5:导出数据,可以是kml或者png之类的图片格式

  1 mapView.exportToPicture(“d:/ddd.png”); //地图导出为图片
  2 //地图导出为kml
  3 lastLayer.saveAsKMLFile(“d:/ddd.kml”);

这是使用这个做出的效果图之一:

2016-06-02 1658

浏览数据秀(dataxiu.com)网站,了解更多数据可视化方法技术。