利用GDAL和QWT绘制图像直方图2--QWT绘制直方图

 转自：http://blog.csdn.net/kupepoem/article/details/20144697 

系统环境：win7

  开发环境：vs2010，qt4.7，qwt6.01，GDAL1.9

参考：..\qwt-6.0.1\examples\tvplot工程

参考：http://blog.csdn.net/e5max/article/details/8854682

        
利用GDAL和QWT绘制图像直方图1--Gdal统计直方图篇我们已经得到了图像的直方图，下面我们利用QWT工具来进行直方图的绘制（QWT的编译和使用请查看我的博客使用qwt作曲线图）。下面我们看看怎么利用QWT绘制GDAL统计出来的直方图。

一、QWT框架简介

        Qwt框架提供了几种标准的图元实现。如网格（QwtPlotGrid类）、标签（QwtPlotMarker类）、标尺（QwtPlotScaleItem类）、显示SVG格式图形文件的图元（QwtPlotSvgItem类）、频谱图（QwtPlotSpectrogram类）、曲线（QwtPlotCurve类）、柱状图（QwtPlotHistogram类）、间隔区域（QwtPlotIntervalCurve类）、三维散点图（QwtPlotSpectroCurve类）等。用户可以继承QwtPlotItem实现符合自己需要的图元。具体请参考http://blog.csdn.net/e5max/article/details/8854682。

二、直方图绘制类实现
QWT继承QT的特点，数据和界面进行分离，具体的我还没深入研究。在此次直方图绘制中Histogram类为数据类，CHisPlot为界面绘制类。我们也不说闲话了，先看直方图绘制类的头文件吧。

//Hisgram为直方图的数据类
class Histogram: public QwtPlotHistogram
{
public:
Histogram(const QString &, const QColor &);

//设置直方图绘制的颜色和符号
void setColor(const QColor &);

//设置直方图的数据
//numValues为统计的间隔数，这在篇1中有介绍
//dfmin为统计图像波段的最小值
//dfmax为统计图像波段的最大值
//pArray为直方图数组，大小为numValues
void setValues(double dfmin,double dfmax,uint numValues, const double *pArray);
};
//实际界面绘制类
class CHisPlot: public QwtPlot
{
Q_OBJECT

public:
CHisPlot(QWidget * = NULL);

//设置直方图数据
void setHis(double dfmin,double dfmax,const QVector<double> & vecHis ,const QString &strHisName);

private:

Histogram *m_hisGram;
};

好了我们看实现吧

//直方图数据类实现
Histogram::Histogram(const QString &title, const QColor &symbolColor):
QwtPlotHistogram(title)
{
//设置柱状直方图
setStyle(QwtPlotHistogram::Columns);
//设置颜色
setColor(symbolColor);
}

void Histogram::setColor(const QColor &symbolColor)
{
QColor color = symbolColor;
color.setAlpha(180);

setPen(QPen(Qt::black));
setBrush(QBrush(color));

QwtColumnSymbol *symbol = new QwtColumnSymbol(QwtColumnSymbol::Box);
symbol->setFrameStyle(QwtColumnSymbol::Raised);
symbol->setLineWidth(2);
symbol->setPalette(QPalette(color));
setSymbol(symbol);
}

//设置数据
void Histogram::setValues(double dfmin,double dfmax,uint numValues, const double *values)
{
QVector<QwtIntervalSample> samples(numValues);

double dscale=(dfmax-dfmin)/numValues;

for ( uint i = 0; i < numValues; i++ )
{
QwtInterval interval(dfmin+i*dscale, (i+1)*dscale+dfmin);//设置统计间隔，这个感觉是QWT的设计的过人之处

interval.setBorderFlags(QwtInterval::ExcludeMaximum);

samples[i] = QwtIntervalSample(values[i], interval);
}

setData(new QwtIntervalSeriesData(samples));
}

CHisPlot::CHisPlot(QWidget *parent):
QwtPlot(parent)
{

//设置背景
setCanvasBackground(QColor(Qt::gray));
plotLayout()->setAlignCanvasToScales(true);
//设置坐标标题
setAxisTitle(QwtPlot::yLeft, tr("Frequency"));
setAxisTitle(QwtPlot::xBottom, tr("Band"));

}

//设置直方图数据
void CHisPlot::setHis(double dfmin,double dfmax,const QVector<double> & vecHis ,const QString & strHisName)
{

m_hisGram = new Histogram(strHisName, Qt::red);

int ncount=vecHis.count();
m_hisGram->setValues( dfmin, dfmax,ncount, &vecHis[0]);
m_hisGram->attach(this);
replot();

etAutoReplot(true);

}

三、直方图的类使用

实现的类写完了，我们看看怎么用CHisPlot类吧。不多说了，先放代码

CHisPlot*pHisPlot=new CHisPlot();
pHisPlot->setHis(pMinMax[0],pMinMax[1],vecHis,"his1");
pHisPlot->show();

四、实现界面效果