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

╰+攻爆jí腚メ 2022-06-09 14:53 664阅读 0赞

转自：[http://blog.csdn.net/kupepoem/article/details/20144697][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][http_blog.csdn.net_e5max_article_details_8854682]

[利用GDAL和QWT绘制图像直方图1--Gdal统计直方图][GDAL_QWT_1--Gdal]篇我们已经得到了图像的直方图，下面我们利用QWT工具来进行直方图的绘制（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][http_blog.csdn.net_e5max_article_details_8854682]。

二、直方图绘制类实现

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

**\[cpp\]** [view plain][] [copy][view plain]

1.  //Hisgram为直方图的数据类
2.  **class** Histogram: **public** QwtPlotHistogram 
3.  \{ 
4.  **public**: 
5.   Histogram(**const** QString &, **const** QColor &); 
6.  
7.  //设置直方图绘制的颜色和符号
8.  **void** setColor(**const** QColor &); 
9.  
10. //设置直方图的数据
11. //numValues为统计的间隔数，这在篇1中有介绍
12. //dfmin为统计图像波段的最小值
13. //dfmax为统计图像波段的最大值
14. //pArray为直方图数组，大小为numValues
15. **void** setValues(**double** dfmin,**double** dfmax,uint numValues, **const****double** \*pArray); 
16. \}; 
17. //实际界面绘制类
18. **class** CHisPlot: **public** QwtPlot 
19. \{ 
20.  Q\_OBJECT 
21. 
22. **public**: 
23.  CHisPlot(QWidget \* = NULL); 
24. 
25. //设置直方图数据
26. **void** setHis(**double** dfmin,**double** dfmax,**const** QVector<**double**> & vecHis ,**const** QString &strHisName); 
27. 
28. 
29. **private**: 
30. 
31.  Histogram \*m\_hisGram; 
32. \};

好了我们看实现吧

**\[cpp\]** [view plain][] [copy][view plain]

1.  //直方图数据类实现
2.  Histogram::Histogram(**const** QString &title, **const** QColor &symbolColor): 
3.   QwtPlotHistogram(title) 
4.  \{ 
5.  //设置柱状直方图
6.   setStyle(QwtPlotHistogram::Columns); 
7.  //设置颜色
8.   setColor(symbolColor); 
9.  \} 
10. 
11. **void** Histogram::setColor(**const** QColor &symbolColor) 
12. \{ 
13.  QColor color = symbolColor; 
14.  color.setAlpha(180); 
15. 
16.  setPen(QPen(Qt::black)); 
17.  setBrush(QBrush(color)); 
18. 
19.  QwtColumnSymbol \*symbol = **new** QwtColumnSymbol(QwtColumnSymbol::Box); 
20.  symbol->setFrameStyle(QwtColumnSymbol::Raised); 
21.  symbol->setLineWidth(2); 
22.  symbol->setPalette(QPalette(color)); 
23.  setSymbol(symbol); 
24. \} 
25. 
26. //设置数据
27. **void** Histogram::setValues(**double** dfmin,**double** dfmax,uint numValues, **const****double** \*values) 
28. \{ 
29.  QVector<QwtIntervalSample> samples(numValues); 
30. 
31. **double** dscale=(dfmax-dfmin)/numValues; 
32. 
33. **for** ( uint i = 0; i < numValues; i++ ) 
34.  \{ 
35.  QwtInterval interval(dfmin+i\*dscale, (i+1)\*dscale+dfmin);//设置统计间隔，这个感觉是QWT的设计的过人之处
36. 
37. 
38.  interval.setBorderFlags(QwtInterval::ExcludeMaximum); 
39. 
40.  samples\[i\] = QwtIntervalSample(values\[i\], interval); 
41.  \} 
42. 
43.  setData(**new** QwtIntervalSeriesData(samples)); 
44. \} 
45. 
46. CHisPlot::CHisPlot(QWidget \*parent): 
47.  QwtPlot(parent) 
48. \{ 
49. 
50. //设置背景
51.  setCanvasBackground(QColor(Qt::gray)); 
52.  plotLayout()->setAlignCanvasToScales(**true**); 
53. //设置坐标标题
54.  setAxisTitle(QwtPlot::yLeft, tr("Frequency")); 
55.  setAxisTitle(QwtPlot::xBottom, tr("Band")); 
56. 
57. 
58. \} 
59. 
60. //设置直方图数据
61. **void** CHisPlot::setHis(**double** dfmin,**double** dfmax,**const** QVector<**double**> & vecHis ,**const** QString & strHisName) 
62. \{ 
63. 
64. 
65.  m\_hisGram = **new** Histogram(strHisName, Qt::red); 
66. 
67. **int** ncount=vecHis.count(); 
68.  m\_hisGram->setValues( dfmin, dfmax,ncount, &vecHis\[0\]); 
69.  m\_hisGram->attach(**this**); 
70.  replot(); 
71. 
72.  etAutoReplot(**true**); 
73. 
74. 
75. \}

三、直方图的类使用

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

**\[cpp\]** [view plain][] [copy][view plain]

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

四、实现界面效果

![SouthEast][]

[http_blog.csdn.net_kupepoem_article_details_20144697]: http://blog.csdn.net/kupepoem/article/details/20144697
[http_blog.csdn.net_e5max_article_details_8854682]: http://blog.csdn.net/e5max/article/details/8854682
[GDAL_QWT_1--Gdal]: http://blog.csdn.net/kupepoem/article/details/20054773
[qwt]: http://blog.csdn.net/kupepoem/article/details/8235204
[view plain]: http://blog.csdn.net/kupepoem/article/details/20144697#
[SouthEast]: /images/20220609/f46aaf92234a4da1b37825e48c74782e.png