标准SEG-Y数据读取及绘图

【实例简介】

本实例基于Qt，设计了用于读取标准SEG-Y数据的segy类、地震道trace类、地震剖面绘制的seismicview类，实现了SEG-Y数据的读取以及地震剖面绘制。

【实例截图】

【核心代码】

/*----------------------SeismicView CLASS----------------------*/
#ifndef SEISMICVIEW_H
#define SEISMICVIEW_H

#include <QGraphicsView>
class SegyClass;
class MainWindow;
class SeismicView : public QGraphicsView
{
     Q_OBJECT

public:
    SeismicView(QWidget *parent=0);
    void SetMainW(MainWindow *pMainWindow);


    bool m_ShowMode;

    bool m_bLBtDown;

    QPoint m_MouseDownPoint;
    QPoint m_MouseCurrPoint;

    float m_startData;
    float m_endData;
    int m_nData;

    int m_nStart;
    int m_nEnd;
    int m_nTrace;
private:
    void destroy();
    void paintEvent(QPaintEvent *  event );
    void mouseMoveEvent(QMouseEvent *mouseEvent);
    void mousePressEvent(QMouseEvent *mouseEvent);
    void mouseReleaseEvent(QMouseEvent *mouseEvent);
private:
    QPixmap pix;
    SegyClass *m_pSegy;
    MainWindow *m_pMainWindow;

private slots:
    void OpenFile();
    void SaveAsBmp();
    void SetShowMode();
    void SaveSgyFile();


};

#endif // SEISMICVIEW_H




#define X_INTER 8
#define Y_INTER 2
#define SAMPLE_INTER 1
#define AMP 10
#define X_START 20
#define Y_START 20

SeismicView::SeismicView(QWidget *parent)
   : QGraphicsView(parent)
{
    m_pSegy=NULL;
    m_pMainWindow=NULL;
    m_ShowMode=true;
    m_bLBtDown=false;

    m_startData=0.0;
    m_endData=0.0;
    m_nData=0;
}
void SeismicView::destroy()
{
    //delete m_pSegy;
    //QGraphicsView::destroy();
}
void SeismicView::SetMainW(MainWindow *pMainWindow)
{
    m_pMainWindow=pMainWindow;


}
void SeismicView::SetShowMode()
{
    m_ShowMode=!m_ShowMode;
    this->viewport()->update();

}
void SeismicView::OpenFile()
{
    m_pSegy=new  SegyClass;
    //m_pMainWindow->ui->actionDoubleCache->setCheckable(1);
    //m_pMainWindow->ui->actionDoubleCache->setChecked(1);

    QPainter paint(viewport());
    QString fileName=QFileDialog::getOpenFileName(this,"Open File","E:\\","sgy(*.sgy *.segy)");

    if (!fileName.isEmpty())
    {
        m_pSegy->ReadSegyFile(fileName);
       // /*
        if(m_ShowMode)
        {
            QRect ur(0,0,m_pSegy->GetWidth(X_INTER),m_pSegy->GetHeight(Y_INTER));
            QPixmap pix1(m_pSegy->GetWidth(X_INTER),m_pSegy->GetHeight(Y_INTER));
            pix=pix1;

            pix.fill(this,ur.topLeft());
            QPainter p(&pix);

            //QPainterPath path;

            p.setPen(QPen(Qt::black,1));
            p.setBrush(QColor(0, 0, 0));

            m_pSegy->DrawMem(&p,X_INTER,20,20,Y_INTER,AMP,1);

            setMinimumSize(m_pSegy->GetWidth(X_INTER),m_pSegy->GetHeight(Y_INTER));
        }

       // */
        this->viewport()->update();

    }
}
void SeismicView::mousePressEvent(QMouseEvent *mouseEvent)
{
    if(mouseEvent->button()==Qt::LeftButton)
    {
        m_bLBtDown = true;
        m_MouseDownPoint=mouseEvent->pos();
        m_nStart=(mouseEvent->y()-Y_START)/Y_INTER;
    }
}
void SeismicView::mouseReleaseEvent(QMouseEvent *mouseEvent)
{


}
void SeismicView::mouseMoveEvent(QMouseEvent *mouseEvent)
{
    if(m_bLBtDown)
    {
       m_MouseCurrPoint=mouseEvent->pos();

    }

    if(m_pSegy!=NULL)
    {
        QString str=QString("Value=%1 TraceN=%2,SampleN=%3").arg(m_pSegy->GetDataValue(mouseEvent->x(),mouseEvent->y(),X_INTER,Y_INTER,20,20,AMP))
                    .arg(m_pSegy->GetTraceN(mouseEvent->x(),mouseEvent->y(),X_INTER,Y_INTER,20,20,AMP))
                    .arg(m_pSegy->GetSampleN(mouseEvent->x(),mouseEvent->y(),X_INTER,Y_INTER,20,20,AMP));
        ((MainWindow*)m_pMainWindow)->statusBar()->showMessage(str);
    }

}
void SeismicView::paintEvent(QPaintEvent *event)
{
    QPainter paint(viewport());



    if(m_pSegy)
    {
        if(m_ShowMode)
        {
            QRect ur(0,0,m_pSegy->GetWidth(X_INTER),m_pSegy->GetHeight(Y_INTER));
            paint.drawPixmap(ur,pix);
        }
        else
        {

            paint.setBrush(QColor(0, 0, 0));
            m_pSegy->DrawMem(&paint,X_INTER,20,20,Y_INTER,AMP,1);
        }
    }

    QGraphicsView::paintEvent(event);
}
void SeismicView::SaveAsBmp()
{

    QString fileName=QFileDialog::getSaveFileName(this,"Open File","E:\\","jpb(*.jpg);;bmp(*.bmp)");
     if(!fileName.isEmpty())
     {
         pix.save(fileName,fileName.right(3).toStdString().data(),1);
     }

}
void SeismicView::SaveSgyFile()
{

    QString fileName=QFileDialog::getSaveFileName(this,"Open File","E:\\","segy(*.sgy);;bmp(*.segy)");
     if(!fileName.isEmpty())
     {
         //pix.save(fileName,fileName.right(3).toStdString().data(),1);
         m_pSegy->SaveSegy(fileName.toStdString().c_str());
     }

}