对绘制气象色斑图createShadedLayer()方法的简单理解（meteoinfo包）

已于 2024-09-02 11:40:06 修改
阅读量1.4k 收藏 32
点赞数 46
文章标签： java
于 2024-09-02 11:29:27 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/q1686821393/article/details/141806748
版权

目录

一、方法概述

  • 包名：org.meteoinfo.geo.meteodata
  • 作用：创建一个带有阴影效果的矢量图层。
  • 返回值：VectorLayer 矢量图层
  • 参数：
参数类型参数名描述
GridDatagridData气象数据
LegendSchemeaLS图例方案
StringlName创建的矢量图层名称
StringfieldName指定了 gridData 中的哪个字段数据将被用于图层渲染
booleanisSmooth是否应用平滑算法

二、源码理解

1. 取出入参的数据

//转换为多边形 POLYGON 图例方案
LegendScheme ls = aLS.convertTo(ShapeTypes.POLYGON);
//获取图例色阶数值范围
Object[] ccs = LegendManage.getContoursAndColors(ls);
double[] cValues = (double[])ccs[0];
//取出气象数据里的最大值、最小值
gridData.getMaxMinValue(maxmin);
double maxData = maxmin[0];
double minData = maxmin[1];
//取出气象数据里的经纬度数组
double[] xArray = gridData.getXArray();
double[] yArray = gridData.getYArray();

2. 生成等值线数组

int[][] S1 = new int[gridData.getYNum()][gridData.getXNum()];
/**
 *	gridData.getData(): 气象格点数据
 *	cValues：图例色阶数值范围
 *	xArray：经度数组
 *	yArray：纬度数组
 *	gridData.getDoubleMissingValue()：缺失值
 *	S1：用于处理中间数据的空数组？
 */
Object[] cbs = ContourDraw.tracingContourLines(gridData.getData(), cValues, xArray, yArray, gridData.getDoubleMissingValue(), S1);
//等值线列表
List<PolyLine> ContourLines = (List)cbs[0];
//等值线图层边界
List<Border> borders = (List)cbs[1];
if (isSmooth) {
  //平滑
  ContourLines = Contour.smoothLines(ContourLines);
}

这里的方法tracingContourLines()，调用了wcontour包里的类方法，来生成等值线结果。（以后分析《关于等值线的绘制方法》）

public static Object[] tracingContourLines(double[][] gridData, double[] cValues, double[] X, double[] Y, double noData, int[][] S1) {
    int nc = cValues.length;
    List<Border> borders = Contour.tracingBorders(gridData, X, Y, S1, noData);
    List<PolyLine> contourLines = Contour.tracingContourLines(gridData, X, Y, nc, cValues, noData, borders, S1);
    return new Object[]{contourLines, borders};
}

3. 创建等值线 POLYGON 列表，新建矢量图层aPolygonShape

//生成等值线多边形 POLYGON 列表
List<Polygon> ContourPolygons = ContourDraw.tracingPolygons(gridData.getData(), ContourLines, borders, cValues);
//创建新图层，图元类型为 POLYGON
VectorLayer aLayer= new VectorLayer(ShapeTypes.POLYGON);
//创建新字段，用于存储每个多边形范围内的最低值
Field aDC = new Field(fieldName + "_Low", DataType.DOUBLE);
aLayer.editAddField(aDC);
//创建新字段，用于存储每个多边形范围内的最高值
aDC = new Field(fieldName + "_High", DataType.DOUBLE);
aLayer.editAddField(aDC);

4. 遍历等值线 POLYGON 列表，将其添加到矢量图层

4.1 开始遍历 POLYGON 列表

Iterator var25 = ContourPolygons.iterator();//迭代器
while(var25.hasNext()) {
	......
}

4.2 获取 POLYGON 顶点列表 pList

//获取每一个Polygon对象
Polygon aPolygon = (Polygon)var25.next();
//获取aPolygon内部的最小值
double aValue = aPolygon.LowValue;
//用于存储多边形的顶点
List<PointD> pList = new ArrayList();
//遍历 aPolygon 的外轮廓点列表
Iterator var29 = aPolygon.OutLine.PointList.iterator();

//将每个点坐标复制到 PointD 对象，并添加到 pList 列表中。
PointD aPoint;
while(var29.hasNext()) {
    wcontour.global.PointD pointList = (wcontour.global.PointD)var29.next();
    aPoint = new PointD();
    aPoint.X = pointList.X;
    aPoint.Y = pointList.Y;
    pList.add(aPoint);
}

//检查点列表是否为顺时针顺序
if (!GeoComputation.isClockwise(pList)) {
    Collections.reverse(pList);
}

4.3 将 POLYGON 顶点列表 pList，添加到多边形对象 aPolygonShape 中

//新建多边形对象
PolygonShape aPolygonShape = new PolygonShape();
aPolygonShape.setPoints(pList); //设置顶点
aPolygonShape.setExtent(GeometryUtil.getPointsExtent(pList)); //设置范围为点列表的边界范围
aPolygonShape.lowValue = aValue;//设置最低值

//aPolygon 有洞（内部空区域）
if (aPolygon.HasHoles()) {
    Iterator var36 = aPolygon.HoleLines.iterator();
	//遍历洞线条列表 HoleLines，并将每个洞的点列表添加到 aPolygonShape 中
    label87:
    while(true) {
        PolyLine holeLine;
        do {
            if (!var36.hasNext()) {
                break label87;
            }

            holeLine = (PolyLine)var36.next();
        } while(holeLine.PointList.size() < 3);

        pList = new ArrayList();
        Iterator var32 = holeLine.PointList.iterator();

        while(var32.hasNext()) {
            wcontour.global.PointD pointList = (wcontour.global.PointD)var32.next();
            aPoint = new PointD();
            aPoint.X = pointList.X;
            aPoint.Y = pointList.Y;
            pList.add(aPoint);
        }

        aPolygonShape.addHole(pList, 0);
    }
}

4.4 处理一些特殊情况

如果 aPolygon 不是高值中心，且高值和低值相同，则将 aPolygonShape 的高值设置为 aValue，并根据 cValues 数组调整低值。

int valueIdx = Arrays.binarySearch(cValues, aValue);
if (valueIdx < 0) {
    valueIdx = -valueIdx - 1;
}

if (valueIdx >= cValues.length - 1) {
    aPolygonShape.highValue = maxData;
} else {
    aPolygonShape.highValue = cValues[valueIdx + 1];
}

if (!aPolygon.IsHighCenter && aPolygon.HighValue == aPolygon.LowValue) {
    aPolygonShape.highValue = aValue;
    if (valueIdx == 0) {
        aPolygonShape.lowValue = minData;
    } else {
        aPolygonShape.lowValue = cValues[valueIdx - 1];
    }
}

4.5 将多边形对象 aPolygonShape 插入到 aLayer 矢量图层中

//获取图层中的图元数量 shapeNum
int shapeNum = aLayer.getShapeNum();
try {
//将 aPolygonShape 插入到 aLayer 矢量图层中的指定位置 shapeNum
    if (aLayer.editInsertShape(aPolygonShape, shapeNum)) {
    //如果插入成功，更新图层中与 aPolygonShape 相关的字段值。
        aLayer.editCellValue(fieldName + "_Low", shapeNum, aPolygonShape.lowValue);
        aLayer.editCellValue(fieldName + "_High", shapeNum, aPolygonShape.highValue);
    }
} catch (Exception var34) {
    Logger.getLogger(DrawMeteoData.class.getName()).log(Level.SEVERE, (String)null, var34);
}

5. 组装矢量图层的数据

//图层名
aLayer.setLayerName(lName);
//图例字段名
ls.setFieldName(fieldName + "_Low");
aLayer.setLegendScheme(ls);
//图层绘制样式：SHADED 阴影填充
aLayer.setLayerDrawType(LayerDrawType.SHADED);

return aLayer;

小白理解，如有错误，请评论指出，欢迎批评。

三、上篇文章遗留问题

问题：在方法 createShadedLayer() 绘制阴影图层中，参数lName、fieldName 的运用逻辑？在参考案例(链接)中直接传递了空字符串。

解答：
lName 用于给生成的阴影图层命名，用来标识图层。在实际应用中，这个参数可以直接传递空字符串，因为它主要用于标识图层，有时候在代码中不需要特别指定
fieldName 图例方案里的字段名，在气象数据中，可能有多个变量（如温度、湿度等），这个参数就是用来选择其中一个变量来显示。如果仅有一个变量，不需要进行显示区分（图例区分），那么也可以直接传递空字符串。

代码：

VectorLayer aLayer = new VectorLayer(ShapeTypes.POLYGON);
aLayer.setLayerName(lName);//图层名
ls.setFieldName(fieldName + "_Low");
aLayer.setLegendScheme(ls);//图例字段名

四、源代码

public static VectorLayer createShadedLayer(GridData gridData, LegendScheme aLS, String lName, String fieldName, boolean isSmooth) {
        LegendScheme ls = aLS.convertTo(ShapeTypes.POLYGON);
        Object[] ccs = LegendManage.getContoursAndColors(ls);
        double[] cValues = (double[])ccs[0];
        double[] maxmin = new double[2];
        gridData.getMaxMinValue(maxmin);
        double maxData = maxmin[0];
        double minData = maxmin[1];
        int[][] S1 = new int[gridData.getYNum()][gridData.getXNum()];
        double[] xArray = gridData.getXArray();
        double[] yArray = gridData.getYArray();
        if (gridData.getXDelta() < 0.0D) {
            ArrayUtils.reverse(xArray);
            gridData.xReverse();
        }

        if (gridData.getYDelta() < 0.0D) {
            ArrayUtils.reverse(yArray);
            gridData.yReverse();
        }

        Object[] cbs = ContourDraw.tracingContourLines(gridData.getData(), cValues, xArray, yArray, gridData.getDoubleMissingValue(), S1);
        List<PolyLine> ContourLines = (List)cbs[0];
        List<Border> borders = (List)cbs[1];
        if (isSmooth) {
            ContourLines = Contour.smoothLines(ContourLines);
        }

        List<Polygon> ContourPolygons = ContourDraw.tracingPolygons(gridData.getData(), ContourLines, borders, cValues);
        VectorLayer aLayer = new VectorLayer(ShapeTypes.POLYGON);
        Field aDC = new Field(fieldName + "_Low", DataType.DOUBLE);
        aLayer.editAddField(aDC);
        aDC = new Field(fieldName + "_High", DataType.DOUBLE);
        aLayer.editAddField(aDC);
        Iterator var25 = ContourPolygons.iterator();

        while(var25.hasNext()) {
            Polygon aPolygon = (Polygon)var25.next();
            double aValue = aPolygon.LowValue;
            List<PointD> pList = new ArrayList();
            Iterator var29 = aPolygon.OutLine.PointList.iterator();

            PointD aPoint;
            while(var29.hasNext()) {
                wcontour.global.PointD pointList = (wcontour.global.PointD)var29.next();
                aPoint = new PointD();
                aPoint.X = pointList.X;
                aPoint.Y = pointList.Y;
                pList.add(aPoint);
            }

            if (!GeoComputation.isClockwise(pList)) {
                Collections.reverse(pList);
            }

            PolygonShape aPolygonShape = new PolygonShape();
            aPolygonShape.setPoints(pList);
            aPolygonShape.setExtent(GeometryUtil.getPointsExtent(pList));
            aPolygonShape.lowValue = aValue;
            if (aPolygon.HasHoles()) {
                Iterator var36 = aPolygon.HoleLines.iterator();

                label87:
                while(true) {
                    PolyLine holeLine;
                    do {
                        if (!var36.hasNext()) {
                            break label87;
                        }

                        holeLine = (PolyLine)var36.next();
                    } while(holeLine.PointList.size() < 3);

                    pList = new ArrayList();
                    Iterator var32 = holeLine.PointList.iterator();

                    while(var32.hasNext()) {
                        wcontour.global.PointD pointList = (wcontour.global.PointD)var32.next();
                        aPoint = new PointD();
                        aPoint.X = pointList.X;
                        aPoint.Y = pointList.Y;
                        pList.add(aPoint);
                    }

                    aPolygonShape.addHole(pList, 0);
                }
            }

            int valueIdx = Arrays.binarySearch(cValues, aValue);
            if (valueIdx < 0) {
                valueIdx = -valueIdx - 1;
            }

            if (valueIdx >= cValues.length - 1) {
                aPolygonShape.highValue = maxData;
            } else {
                aPolygonShape.highValue = cValues[valueIdx + 1];
            }

            if (!aPolygon.IsHighCenter && aPolygon.HighValue == aPolygon.LowValue) {
                aPolygonShape.highValue = aValue;
                if (valueIdx == 0) {
                    aPolygonShape.lowValue = minData;
                } else {
                    aPolygonShape.lowValue = cValues[valueIdx - 1];
                }
            }

            int shapeNum = aLayer.getShapeNum();

            try {
                if (aLayer.editInsertShape(aPolygonShape, shapeNum)) {
                    aLayer.editCellValue(fieldName + "_Low", shapeNum, aPolygonShape.lowValue);
                    aLayer.editCellValue(fieldName + "_High", shapeNum, aPolygonShape.highValue);
                }
            } catch (Exception var34) {
                Logger.getLogger(DrawMeteoData.class.getName()).log(Level.SEVERE, (String)null, var34);
            }
        }

        aLayer.setLayerName(lName);
        ls.setFieldName(fieldName + "_Low");
        aLayer.setLegendScheme(ls);
        aLayer.setLayerDrawType(LayerDrawType.SHADED);
        return aLayer;
    }
宰宰陪我学Java
关注
Leaflet、kriging生成色斑图
04-27 1046
Leaflet、kriging生成色斑图渲染至canvas上，再将canvas对象转换成image图片加载到地图上，此方法放大时效果不是太好，范围连线、等值线不够平滑，有锯齿。
MeteoInfo绘图工具包，wContour绘图包，netcdfAll等绘图相关所有jar，MeteoInfoMap.exe等可执行文件，wContour中包
07-09
MeteoInfo绘图工具包，wContour绘图包，netcdfAll等绘图相关所有jar，MeteoInfoMap.exe等可执行文件，wContour中包含反距离加权算法（IDW），克里金插值算法(Kriging)
对于使用meteoinfo包绘制气象色斑图的理解
最新发布
无人区发疯的博客
08-30 1842
*** @Param grid 气象数据* @Param vectorLayer 矢量图层* @Param outPath 图片输出路径* @Param conf 图片配置信息*///读取色阶数值断点范围//1.用于创建一个渐变色图例方案//2.绘制图层//3.叠加图层，剪裁到合适地区//创建视图和布局//创建视图，添加底层地图shp文件和气象色斑图图层//创建布局，设置页面大小、图层排列等//返回视图的地理边界范围。
MeteoInfo-Java解析与绘图教程
张志强的博客
12-08 992
MeteoInfo-Java解析与绘图教程
色斑图程序
04-26
降雨量色斑图程序
meteoinfo必须要用java_发布MeteoInfo Java 1.2.2
weixin_33586063的博客
02-26 1832
主要更新了MeteoInfoLab至0.2版，已经有一些实用功能了，这里做些简单的介绍。MeteoInfoLab是以脚本编写和命令行交互为主的软件，用Jython语言对MeteoInfo库进行了封装，提供科学计算和绘图的功能，函数参照MatLab、Numpy、Matplotlib实现。目前只实现了很有限的功能。双击MeteoInfoLab.exe启动(Linux下用milab.sh启动)，界面中包...
制作气象格点数据NC文件 制作色斑图根据NC文件生成色班图 制作天气预报图
01-11
制作气象格点数据NC文件 根据NC文件生成色班图 制作天气预报图 netcdf openlayer java python GeoJSON basemap
MeteoInfo(气象制图软件)安装包
08-31
MeteoInfo是一个开源的 GIS 和科学计算环境软件（LGPL许可），其中包括两个主要应用：MeteoInfoMap 和 MeteoInfoLab 。 MeteoInfoMap 旨在交互地查看和分析气象和空间数据，它有两个版本：Java 和 C＃ 。 ...
MeteoInfo_Help_1.0.1_Chinese.rar_MeteoInfo_meteoinfo学习_气象_气象软件
07-15
meteoinfo的中文使用说明，很强大的一个气象软件，学气象的亲可以学习使用这个软件
25、气象填色图绘制
小兔子乖乖
07-20 1160
气象填色图（也称为气象色斑图）是气象学中常用的一种图表，用于展示温度、湿度、气压等气象要素的空间分布。这种图表通过不同的颜色或色阶来表示不同的气象数据值，使得观察者能够直观地理解数据的分布情况。
java做的天气预报图形界面程序
10-26
包含源代码与用exe4j制作的exe程序，自带jre，jre已经做过初步的删减。
气象云图数据解析代码之Java实现
04-14
Java语言实现解析气象云图数据，获得良好的效果，能够实现对AWX格式的气象云图数据精准解析，测试效果良好！
降雨量色斑图程序
01-05
非常专业的降雨量色斑图程序，值得拥有和收藏
热力图生成代码，适用于各种色斑图生成
10-21
热力图代码，适用各种色斑图生成，导入vs直接运行
MeteoInfo软件
07-07
很好用的一个气象绘图软件，能够很方便的绘制各种阴影图、站点图等。
Java图案绘制
qq_36159785的博客
11-20 1343
图案绘制主要是通过在JPanel类的子类中，重写paint()方法，并在该方法中使用Graphics2D类中的方法来实现的。 先上个效果图： —ps：其实很简单，真不知道该怎么描述了。。。在main方法中绘制一个显示图案的窗口：JFrame window = new DrawFrame(); window.setTitle("绘制图形"); window.setDefaultC
MeteoInfo 生成等值面
随风落叶
08-02 1783
以上是效果图，下面是测试代码： /** * 站点Bean * @author Administrator * */ public class StationBean { private String id ; //ID private double x; //X坐标 private double y; //y坐标 private double v; //值 publ...
MeteoInfo-Micaps绘制色斑图
m0_61759012的博客
12-09 1070
meteoInfo处理色斑图
MeteoInfo：免费气象数据分析软件
"MeteoInfo用户手册 - 免费气象数据显示和分析软件" MeteoInfo是一款专为气象数据处理设计的免费软件，它提供强大的数据图形显示和分析功能，适用于各种气象数据格式，并带有基本的GIS（地理信息系统）功能。软件的...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值