等值线算法（转载）

等值线算法

对于二维的规则网格，本文以104*104格点大小数据做为说明，数据下载：http://yunpan.cn/ccB7Y36eLXk35  访问密码 00c3 ，Qcloud.txt为原始数据，re01.txt为01代码

一.算法如下：

1.在该层的数据网格中求出所有四个相邻的数据点构成的正方形；

2.判断四个数据点数据与阈值之间的关系，生成01代码；

3.由上步生成的代码按照下图的关系求出等值线与四个数据点之间的拓扑关系；

          

4.由拓扑关系用线性插值方法求出等值线与正方形的交点；

5.顺序连接等值线段，即得到等值线。

二.源代码如下：

头文件atest.h：

#ifndef ATEST_H
#define ATEST_H
 
#include <QtGui/QtGui>
#include "ui_atest.h"
#include <iostream>
#include <stdlib.h>
#include <fstream>
#include <string>
#include <sstream>
#include <vector>
 
class ATest : public QMainWindow,public Ui_ATestClass
{
    Q_OBJECT
 
public:
    ATest(QWidget *parent = 0, Qt::WFlags flags = 0);
    ~ATest();
 
private:
    QString filename;
    std::vector<float> result;//存放所有矩形插值得到的位置数据
    bool drawtf;
    void hookupsignals();
    void paintEvent(QPaintEvent *);
private slots:
    void filesearch();
    void getParam();
};
 
#endif // ATEST_H

实现文件atest.cpp：

#include "atest.h"
using namespace std;
 
ATest::ATest(QWidget *parent, Qt::WFlags flags)
    : QMainWindow(parent, flags)
{
    setupUi(this);
    hookupsignals();
    drawtf=false;
}
 
ATest::~ATest()
{
 
}
 
//*************************************
// 时    间:  2015/7/17 10:44
// 权    限:  public
// 返    回:  std::vector<std::string>
// 方法说明:  将string转化成字符串数组
//*************************************
std::vector<std::string> split(std::string str,std::string pattern)
{ 
    std::string::size_type pos;  
    std::vector<std::string> result; 
    str+=pattern;//扩展字符串以方便操作  
    int size=str.size();    
    for(int i=0; i<size; i++)  
    {    
        pos=str.find(pattern,i);    
        if(pos<size)    
        {     
            std::string s=str.substr(i,pos-i);      
            result.push_back(s);     
            i=pos+pattern.size()-1;    
        }  
    }  
    return result;
}
 
//*************************************
// 时    间:  2015/7/18 21:22 
// 权    限:  public
// 返    回:  float
// 方法说明:  返回最终的插值结果
//*************************************
float execute(float yuzhi,float vmin,float vmax,float locmin,float locmax)
{
    float result=(locmax-locmin)*(yuzhi-vmin)/(vmax-vmin);
    return result;
}
 
//*************************************
// 时    间:  2015/7/17 15:00
// 权    限:  public
// 返    回:  vector<float>
// 方法说明:  插值函数，返回插值得到的点的坐标
//*************************************
vector<float> interpolate(vector<int> flag,vector<float> data,vector<float> locate)
{
    vector<float> result;
    float t;
    if(flag[0]==0)
    {
        if(flag[1]==0)
        {
            if(flag[2]==0)
            {
                if(flag[3]==0)//0000
                {              
                }  
                else//0001
                {
                    t=execute(data[4],data[2],data[3],locate[4],locate[6]);
                    result.push_back(locate[5]);
                    result.push_back(locate[4]+t);
                    t=execute(data[4],data[3],data[0],locate[7],locate[1]);
                    result.push_back(locate[7]+t);
                    result.push_back(locate[6]);   
                }
            }