selectdlg.cpp

显示地震剖面对地震专业的学生来说显示地震剖面特别重要本程序实现了地震剖面的绘制

// SelectDLG.cpp : implementation file
//

#include "stdafx.h"
#include "demo.h"
#include "SelectDLG.h"
#include "stdio.h"
#include "math.h"
#include "triangle.h"

#ifdef _DEBUG
#define new DEBUG_NEW
#undef THIS_FILE
static char THIS_FILE[] = __FILE__;
#endif

/////////////////////////////////////////////////////////////////////////////
// SelectDLG dialog


SelectDLG::SelectDLG(CWnd* pParent /*=NULL*/)
	: CDialog(SelectDLG::IDD, pParent)
{
	//{{AFX_DATA_INIT(SelectDLG)
	m_2dpos = _T("");
	m_3dpos = _T("");
	m_texpos = _T("");
	m_tripos = _T("");
	//}}AFX_DATA_INIT
}


void SelectDLG::DoDataExchange(CDataExchange* pDX)
{
	CDialog::DoDataExchange(pDX);
	//{{AFX_DATA_MAP(SelectDLG)
	DDX_Text(pDX, IDC_2DPOS, m_2dpos);
	DDX_Text(pDX, IDC_3DPOS, m_3dpos);
	DDX_Text(pDX, IDC_TEXPOS, m_texpos);
	DDX_Text(pDX, IDC_TRIPOS, m_tripos);
	//}}AFX_DATA_MAP
}


BEGIN_MESSAGE_MAP(SelectDLG, CDialog)
	//{{AFX_MSG_MAP(SelectDLG)
	ON_BN_CLICKED(IDC_LOAD2D, OnLoad2d)
	ON_BN_CLICKED(IDC_LOAD3D, OnLoad3d)
	ON_BN_CLICKED(IDC_LOADTEX, OnLoadtex)
	ON_BN_CLICKED(IDC_LOADTRI, OnLoadtri)
	//}}AFX_MSG_MAP
END_MESSAGE_MAP()

/////////////////////////////////////////////////////////////////////////////
// SelectDLG message handlers

void SelectDLG::OnLoad2d() 
{
	// TODO: Add your control notification handler code here
	static char szFilter[]="TXT文件(*.txt)|*.txt||";
	CFileDialog m_ldFile(TRUE,"txt","2Ddatafile",OFN_HIDEREADONLY|OFN_OVERWRITEPROMPT,szFilter);
 
    if (m_ldFile.DoModal()==IDOK)
	{   
		m_2dpos=m_ldFile.GetPathName();
		UpdateData(FALSE);		
	//}

	//读入txt文件
	n2dpoint=0;
	for(int i=0;i<MAXPOINT;i++)
	{
		m_point2d[i][0]=0;
		m_point2d[i][1]=0;
		m_point2d_old[i][0]=0;
		m_point2d_old[i][1]=0;
	}
	FILE* p2dfile;
	if((p2dfile=fopen(m_2dpos,"r"))==NULL)
		MessageBox("2D数据文件打开错误!",NULL,MB_ICONERROR|MB_SYSTEMMODAL);
	else
	{
		for(i=0;i<MAXPOINT;i++)
		{
			fscanf(p2dfile,"%f",&m_point2d[i][0]);
			fscanf(p2dfile,"%f",&m_point2d[i][1]);
			m_point2d_old[i][0]=m_point2d[i][0];
			m_point2d_old[i][1]=m_point2d[i][1];
			if((m_point2d[i][0]==0)&&(m_point2d[i][1]==0))
				break;
			n2dpoint++;
		}
		fclose(p2dfile);
		if(n2dpoint==MAXPOINT)
		{
			MessageBox("数据点数量超过程序最大可显示数量，可能导致图像显示不完整!",NULL,
			MB_ICONWARNING|MB_SYSTEMMODAL);
			return;
		}

		//对读入的二维点进行三角化
		int kk=0;
		int &numberoftriangle=kk;		
		int* vertexindex=new int[MAXPOINT*10];

		triangulation(m_point2d,n2dpoint,vertexindex,numberoftriangle);
		ntri=numberoftriangle;

		triangle=new int*[ntri];
		for(i=0;i<ntri;i++)
		{
			triangle[i]=new int[3];
			for(int j=0;j<3;j++)
				triangle[i][j]=vertexindex[i*3+j];		
		}
		delete vertexindex;
	}//else
	}//if (m_ldFile.DoModal()==IDOK)
}

void SelectDLG::OnLoad3d() 
{
	// TODO: Add your control notification handler code here
	static char szFilter[]="TXT文件(*.txt)|*.txt||";
	CFileDialog m_ldFile(TRUE,"txt","3Ddatafile",OFN_HIDEREADONLY|OFN_OVERWRITEPROMPT,szFilter);

	if (m_ldFile.DoModal()==IDOK)
	{
	 	m_3dpos=m_ldFile.GetPathName();
		UpdateData(FALSE);
	//}

	//读入txt文件
	n3dpoint=0;
	for(long i=0;i<MAXPOINT;i++)
	{
		m_point3d[i][0]=0;
		m_point3d[i][1]=0;
		m_point3d[i][2]=0;
		m_point3d_old[i][0]=0;
		m_point3d_old[i][1]=0;
		m_point3d_old[i][2]=0;
	}
	FILE* p3dfile;
	if((p3dfile=fopen(m_3dpos,"r"))==NULL)
		MessageBox("3D数据文件打开错误!",NULL,MB_ICONERROR|MB_SYSTEMMODAL);
	else
	{
		for(i=0;i<MAXPOINT;i++)
		{
			fscanf(p3dfile,"%f",&m_point3d[i][0]);
			fscanf(p3dfile,"%f",&m_point3d[i][1]);	
			fscanf(p3dfile,"%f",&m_point3d[i][2]);
			m_point3d_old[i][0]=m_point3d[i][0];
			m_point3d_old[i][1]=m_point3d[i][1];
			m_point3d_old[i][2]=m_point3d[i][2];
			if((m_point3d[i][0]==0)&&(m_point3d[i][1]==0)&&(m_point3d[i][2]==0))
				break;
			n3dpoint++;
		}
		fclose(p3dfile);
		if(n3dpoint==MAXPOINT)
		{
			MessageBox("数据点数量超过程序最大可显示数量，可能导致图像显示不完整!",NULL,
			MB_ICONWARNING|MB_SYSTEMMODAL);
			return;
		}

		//归一化3D点坐标，使其形心在原点且各顶点到形心的平均距离为1
		//计算xyz三坐标的均值
		float mean[3]={0,0,0};
		for(i=0;i<n3dpoint;i++)
		{
			mean[0]+=m_point3d[i][0];
			mean[1]+=m_point3d[i][1];
			mean[2]+=m_point3d[i][2];
		}
		mean[0]=mean[0]/n3dpoint;
		mean[1]=mean[1]/n3dpoint;
		mean[2]=mean[2]/n3dpoint;
		//使点的形心在原点
		for(i=0;i<n3dpoint;i++)
		{
			m_point3d[i][0]=m_point3d[i][0]-mean[0];
			m_point3d[i][1]=m_point3d[i][1]-mean[1];
			m_point3d[i][2]=m_point3d[i][2]-mean[2];
		}
		//使它们到原点的平均距离的平方是1
		float meandis=0;
		for(i=0;i<n3dpoint;i++)
		{
			meandis+=(float)pow(m_point3d[i][0]*m_point3d[i][0]+m_point3d[i][1]*m_point3d[i][1]+m_point3d[i][2]*m_point3d[i][2],0.5);
		}
		meandis=meandis/n3dpoint;
		for(i=0;i<n3dpoint;i++)
		{
			m_point3d[i][0]=m_point3d[i][0]/meandis;
			m_point3d[i][1]=m_point3d[i][1]/meandis;
			m_point3d[i][2]=m_point3d[i][2]/meandis;
		}
	}//else
	}//if (m_ldFile.DoModal()==IDOK)
}

void SelectDLG::OnLoadtex() 
{
	// TODO: Add your control notification handler code here
	CFileDialog m_ldFile(TRUE,"*","texturefile",NULL,"Image Files (*.jpg;*.bmp)|*.jpg;*.bmp||");

	if (m_ldFile.DoModal()==IDOK)
	{
		m_texpos=m_ldFile.GetPathName();
		UpdateData(FALSE);
		if(!mDib.Load((char*)(LPCTSTR)m_texpos))
			MessageBox("纹理图像打开错误!",NULL,MB_ICONERROR|MB_SYSTEMMODAL);
		else
		{
			//调整纹理图像的大小使之为2的幂(此处均设为512x512)
			mDib2=mDib;
			ImageSize.cx=mDib.CX();
			ImageSize.cy=mDib.CY();
			mDib2.Change(512,512,24,NULL);
			GLubyte *pTextureBits=(GLubyte *)mDib.Buffer();
			GLubyte *pTextureBits2=(GLubyte *)mDib2.Buffer();
			gluScaleImage(GL_BGR_EXT,ImageSize.cx,ImageSize.cy,GL_UNSIGNED_BYTE,pTextureBits,512,512,GL_UNSIGNED_BYTE,pTextureBits2);
		}
	}	
}

void SelectDLG::OnLoadtri() 
{
	// TODO: Add your control notification handler code here
	static char szFilter[]="TXT文件(*.txt)|*.txt||";
	CFileDialog m_ldFile(TRUE,"txt","tridatafile",OFN_HIDEREADONLY|OFN_OVERWRITEPROMPT,szFilter);
 
    if (m_ldFile.DoModal()==IDOK)
	{   
		m_tripos=m_ldFile.GetPathName();
		UpdateData(FALSE);
		
		//读入txt文件
		ntri=0;
		delete triangle;
		for(int i=0;i<MAXPOINT;i++)
		{
			tritemp[i][0]=0;
			tritemp[i][1]=0;
			tritemp[i][2]=0;
		}
		FILE* ptrifile;
		if((ptrifile=fopen(m_tripos,"r"))==NULL)
			MessageBox("三角化索引数据打开错误!",NULL,MB_ICONERROR|MB_SYSTEMMODAL);
		else
		{
			for(i=0;i<MAXPOINT;i++)
			{
				fscanf(ptrifile,"%d",&tritemp[i][0]);
				fscanf(ptrifile,"%d",&tritemp[i][1]);
				fscanf(ptrifile,"%d",&tritemp[i][2]);
				if((tritemp[i][0]==0)&&(tritemp[i][1]==0)&&(tritemp[i][2]==0))
					break;
				ntri++;
			}
			fclose(ptrifile);
			triangle=new int*[ntri];
			for(i=0;i<ntri;i++)
			{
				triangle[i]=new int[3];
				for(int j=0;j<3;j++)
					triangle[i][j]=tritemp[i][j]-1;		
			}
		}//else
	}//if (m_ldFile.DoModal()==IDOK)	
}

//二维点三角化函数/////////////////////////////////////////////////////////////////////
void SelectDLG::triangulation(float m_point2d[MAXPOINT][2], int numberofcorner, int vertexindex[], int & numberoftriangle)
{
  triangulateio_s in,mid,vorout;
  in.numberofpoints = numberofcorner;
  in.numberofpointattributes = 0;
  in.pointlist = (float  *) malloc(in.numberofpoints * 2 * sizeof(float ));
  int i;
  for(i=0;i<in.numberofpoints;i++)
  {
  	  in.pointlist[i*2]=m_point2d[i][0];
	  in.pointlist[i*2+1]=m_point2d[i][1];
  }

  in.numberofsegments = 0;
  in.numberofholes = 0;
  in.pointmarkerlist = (int *) NULL; 
  mid.pointlist = (float  *) NULL;            
  mid.pointattributelist = (float  *) NULL;
  mid.pointmarkerlist = (int *) NULL; 
  mid.trianglelist = (int *) NULL;        
  mid.triangleattributelist = (float  *) NULL;
  mid.neighborlist = (int *) NULL;         
  mid.segmentlist = (int *) NULL;
  mid.segmentmarkerlist = (int *) NULL;
  mid.edgelist = (int *) NULL;            
  mid.edgemarkerlist = (int *) NULL; 
  vorout.pointlist = (float  *) NULL;        
  vorout.pointattributelist = (float  *) NULL;
  vorout.edgelist = (int *) NULL;         
  vorout.normlist = (float  *) NULL;         
  triangulate("ezv", &in, &mid, &vorout);
  numberoftriangle=mid.numberoftriangles;

  for (i = 0; i < mid.numberoftriangles; i++) 
  {	 
	  for (int j = 0; j < mid.numberofcorners; j++)
	  {		  
		  vertexindex[i * mid.numberofcorners + j]=
			  mid.trianglelist[i * mid.numberofcorners + j];
	  }	 
  }
 // fclose(angle);
  free(in.pointlist);
  free(mid.pointlist);
  free(mid.pointattributelist);
  free(mid.pointmarkerlist);
  free(mid.trianglelist);
  free(mid.neighborlist);
  free(mid.segmentlist);
  free(mid.segmentmarkerlist);
  free(mid.edgelist);
  free(mid.edgemarkerlist);
  free(vorout.pointlist);
  free(vorout.pointattributelist);
  free(vorout.edgelist);
  free(vorout.normlist);
}