专利名称：一种等值线插值方法
技术领域：
本发明涉及等值线实时生成过程中的插值方法。
背景技术：
等值线(Contour)来源于测绘学上的等高线，是用一组光滑曲线来表示连续面状分布的制图现象数量特征渐变的方法。在工程应用领域，等值线被广泛应用于用一组监测数据反应大自然实际变化状态的最直观的工具，如降雨范围、温度变化、土壤墒情变化、病虫害分布等。目前能够进行等值线分析及绘制的方法如Surfer、Matlab、ArcGis软件中涉及的多种方法，或如中国发明专利申请CN 200710016657. 7、中国发明专利ZL 200510066109. 6 中公开的方法。这些方法中，等值线计算一般采用插值法(内插法)，即将计算区域划分为网格，然后向各个网格投影已有数据的散点，每个网格上得到一组数据，将具有相同数据的网格点连接起来，再通过圆滑处理获得等值线。把等值线分割出的区域用不同颜色进行渲染而得到等值面。实际工作中，这些软件或工具生成的等值线一般都存在以下问题1.不能在可视化界面实时修正等值线。2.生成等值线不能进行非锯齿化无缝边界裁剪或者速度较慢。3.生成的等值线图输出非矢量化，或者矢量化但没有属性信息。
发明内容
(一)要解决的技术问题本发明的目的是要克服上述不足，提供一种在可视化地图中，在线改变离散点，即时插值，修改等值线的方法。( 二 )技术方案为了达到上述目的，本发明提供了一种等值线生成方法，所述方法包括给定平面中的η个散点，在指定区域内对该区域进行规则矩形网格剖分，利用带权二次多项式插值算法获得各网格节点的值，然后根据网格节点坐标及值，跟踪绘制出等值线，具体步骤如下(1)指定等值线区域，为一矩形范围[MinX，MinY]，[MaxX, MaxY]；(2)把该矩形区域定义为(m-1) X (n_l)个网格，各网格点的序号排列为横向j =1，2，3. ..η;纵向i = 1，2，3. ..m，每一个网格的横边长为Dx，纵边长为Dy;(3)使用带权的二次多项式算法，插值出各网格点的值；(4)找出各等值点并连接生成该矩形区域内的等值线；(5)给定实际工作中有效的区域多边形，对矩形区域内的等值线进行裁剪，最终得到有效区域内的等值线图；以上步骤生成的等值线，都具有矢量信息，每条线都是具有实际意义的坐标所组成，并且具有等值线数值，可将该矢量信息进行输出。
其中，所述带权的二次多项式插值算法步骤为获取所述等值线图片中平面点内 η个观测点的坐标，P0 (x0, y0), Pl(xi; yi)…Plri (Xlri，y^)，所述η个点(η彡7)的函数值为
权利要求
1 一种等值线生成方法，所述方法包括给定平面中的η个散点，在指定区域内对该区域进行规则矩形网格剖分，利用带权二次多项式插值算法获得各网格节点的值，然后根据网格节点坐标及值，跟踪绘制出等值线，具体步骤如下(1)指定等值线区域，为一矩形范围[MinX，MinY]，[MaxX,MaxY]；(2)把该矩形区域定义为(m-1)X (n-1)个网格，各网格点的序号排列为横向j = 1， 2，3. ..η;纵向i = 1，2，3. ..m，每一个网格的横边长为Dx，纵边长为Dy;(3)插值获得各网格点的值；(4)找出各等值点并连接生成该矩形区域内的等值线；(5)给定实际工作中有效的区域多边形，对矩形区域内的等值线进行裁剪，最终得到有效区域内的等值线图；其特征在于步骤(3)中，使用带权的二次多项式算法进行插值，步骤为获取所述等值线图片中平面点内η个观测点的坐标，η彡7，P0 (x0, y0)，Ρι(χ1 Yl)……Plri (Xlri，yn_i)，所述 η个点的函数值为Ztl，Z1, -Zn^1 ；对待插值点F(a，b)进行双曲面插值，双曲面插值多项式为
2.根据权利要求1所述的等值线生成方法，其特征在于通过在步骤3中改变η个点的值，进行重新插值。
3.根据权利要求2所述的等值线生成方法，其特征在于在地图中改变离散点包括新增离散点、删除离散点、修改离散点的值、移动离散点。
全文摘要
本发明涉及一种在等值线图片中改变等值线插值点点后修改等值线的方法，通过获取所述等值线图片中平面点内n个点的坐标，对待插值点p(x，y)进行双曲面插值。本发明的方法由于在多项式插值算法中引用了权，能有效地改善由于数据存在明显的趋势和数据空间分布不以致造成的插值误差。本发明还涉及建立了非退化直线方法，并在此基础上，开发了点与多边形关系(点在多边形内)新的算法，该算法用于等值线边界裁减中，降低运算时间。
文档编号G06F19/00GK102289605SQ20111026031
公开日2011年12月21日 申请日期2011年9月5日 优先权日2011年9月5日
发明者焦中志, 翟彩霞, 薛伟 申请人:北京清流技术发展有限公司