c++ 手工做一个STL文件

【实例截图】

【核心代码】

// AuzDlg.cpp : implementation file
//

#include "stdafx.h"
#include "Auz.h"
#include "AuzDlg.h"
#include "math.h"

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

//#define pi 3.1415926535898

struct FXYZ{float x,y,z;};
struct FACET{FXYZ n;FXYZ dot[3];BYTE b[2];};
struct IXYZ{WORD d[3];};


FXYZ normal(FXYZ v0,FXYZ v1);
/////////////////////////////////////////////////////////////////////////////
// CAuzDlg dialog

CAuzDlg::CAuzDlg(CWnd* pParent /*=NULL*/)
	: CDialog(CAuzDlg::IDD, pParent)
{
	//{{AFX_DATA_INIT(CAuzDlg)
		// NOTE: the ClassWizard will add member initialization here
	//}}AFX_DATA_INIT
	m_hIcon = AfxGetApp()->LoadIcon(IDR_MAINFRAME);
}

void CAuzDlg::DoDataExchange(CDataExchange* pDX)
{
	CDialog::DoDataExchange(pDX);
	//{{AFX_DATA_MAP(CAuzDlg)
	DDX_Control(pDX, IDC_LIST1, m_lst);
	//}}AFX_DATA_MAP
}

BEGIN_MESSAGE_MAP(CAuzDlg, CDialog)
	//{{AFX_MSG_MAP(CAuzDlg)
	ON_WM_PAINT()
	ON_WM_QUERYDRAGICON()
	ON_BN_CLICKED(IDC_BUTTON1, OnButton1)
	ON_BN_CLICKED(IDC_BUTTON2, OnButton2)
	//}}AFX_MSG_MAP
END_MESSAGE_MAP()

/////////////////////////////////////////////////////////////////////////////
// CAuzDlg message handlers

BOOL CAuzDlg::OnInitDialog()
{
	CDialog::OnInitDialog();

	SetIcon(m_hIcon, TRUE);			// Set big icon
	SetIcon(m_hIcon, FALSE);		// Set small icon
	
	// TODO: Add extra initialization here
	
	return TRUE;  // return TRUE  unless you set the focus to a control
}

// If you add a minimize button to your dialog, you will need the code below
//  to draw the icon.  For MFC applications using the document/view model,
//  this is automatically done for you by the framework.

void CAuzDlg::OnPaint() 
{
	if (IsIconic())
	{
		CPaintDC dc(this); // device context for painting

		SendMessage(WM_ICONERASEBKGND, (WPARAM) dc.GetSafeHdc(), 0);

		// Center icon in client rectangle
		int cxIcon = GetSystemMetrics(SM_CXICON);
		int cyIcon = GetSystemMetrics(SM_CYICON);
		CRect rect;
		GetClientRect(&rect);
		int x = (rect.Width() - cxIcon   1) / 2;
		int y = (rect.Height() - cyIcon   1) / 2;

		// Draw the icon
		dc.DrawIcon(x, y, m_hIcon);
	}
	else
	{
		CDialog::OnPaint();
	}
}

HCURSOR CAuzDlg::OnQueryDragIcon()
{
	return (HCURSOR) m_hIcon;
}

void CAuzDlg::OnButton1() 
{
	// TODO: Add your control notification handler code here
	char st[80],head[120];
	int i,j;
	double x,y,z,r;
	FXYZ v0,v1,v2;
	long TtlFct=29;

	FACET fc[29]={0};
//struct FACET{FXYZ n;FXYZ dot[3];BYTE b[2];};	//法向量，三节点，保留
	FXYZ jd[16]={-50,50,25,-25,50,25,-25,25,25,25,25,25,25,50,25,50,50,25,50,0 ,25,-50,0,25,
		-50,50,-25,-25,50,-25,-25,25,-25,25,25,-25,25,50,-25,50,50,-25,50,0 ,-25,-50,0,-25};

	IXYZ N3[29]={0,1,2,0,2,7,2,3,6,2,6,7,3,4,5,3,5,6,8,10,9,8,15,10,10,14,11,10,15,14,10,15,14,11,13,12,11,14,13,
		5,4,12,5,12,13,1,0,8,1,8,9,0,7,15,0,15,8,6,5,13,6,13,14,7,6,14,7,14,15,
		3,2,10,3,10,11,2,1,9,2,9,10,4,3,11,4,11,12};

//	for(i=0;i<16;i  ){sprintf(st,"%02d %9.3f %9.3f %9.3f",i,jd[i].x,jd[i].y,jd[i].z);m_lst.AddString(st);}m_lst.AddString("");
//	for(i=0;i<TtlFct;i  ){sprintf(st,"%02d %2d %2d %2d",i,N3[i].d[0],N3[i].d[1],N3[i].d[2]);m_lst.AddString(st);}m_lst.AddString("");
	FILE *fp;
	fp=fopen("凹.stl","wb");
	sprintf(head,"这是一个手工编制的STL格式文件俯视图是一个汉字凹字");
	fwrite(head,1,80,fp);
	fwrite(&TtlFct,4,1,fp);

	for(i=0;i<TtlFct;i  )
	{
//		m_lst.AddString("");for(j=0;j<3;j  ){sprintf(st,"%02d %9.3f %9.3f %9.3f",i,jd[N3[i].d[j]].x,jd[N3[i].d[j]].y,jd[N3[i].d[j]].z);m_lst.AddString(st);}
		x=jd[N3[i].d[1]].x-jd[N3[i].d[0]].x;
		y=jd[N3[i].d[1]].y-jd[N3[i].d[0]].y;
		z=jd[N3[i].d[1]].z-jd[N3[i].d[0]].z;
		r=sqrt(x*x y*y z*z);v0.x=x/r;v0.y=y/r;v0.z=z/r;
		x=jd[N3[i].d[2]].x-jd[N3[i].d[0]].x;
		y=jd[N3[i].d[2]].y-jd[N3[i].d[0]].y;
		z=jd[N3[i].d[2]].z-jd[N3[i].d[0]].z;
		r=sqrt(x*x y*y z*z);v1.x=x/r;v1.y=y/r;v1.z=z/r;
		fc[i].n=normal(v0,v1);
		fc[i].dot[0]=jd[N3[i].d[0]];
		fc[i].dot[1]=jd[N3[i].d[1]];
		fc[i].dot[2]=jd[N3[i].d[2]];

		m_lst.AddString("");sprintf(st,"n: i=%02d %9.6f %9.6f %9.6f",i,fc[i].n.x,fc[i].n.y,fc[i].n.z);m_lst.AddString(st);
		for(j=0;j<3;j  ){sprintf(st,"j=%2d %9.3f %9.3f %9.3f",j,fc[i].dot[j].x,fc[i].dot[j].y,fc[i].dot[j].z);m_lst.AddString(st);}
	}
//	m_lst.AddString("");
	fwrite(fc,sizeof(FACET),TtlFct,fp);fclose(fp);


}


void CAuzDlg::OnButton2() 
{
	// TODO: Add your control notification handler code here
	char st[120];
	long TtlFct,k;
	int i;

	FILE *fp;
	if(!(fp=fopen("凹.stl","rb"))){m_lst.AddString("找不到凹.stl");return;}
	fread(st,1,80,fp);m_lst.AddString(st);
	fread(&TtlFct,4,1,fp);sprintf(st,"面片数：%d",TtlFct);m_lst.AddString(st);
	FACET *fc=(FACET *)new FACET[TtlFct];
	fread(fc,sizeof(FACET),TtlFct,fp);fclose(fp);

	m_lst.ResetContent();
	for(k=0;k<TtlFct;k  )
	{
		m_lst.AddString("");
		sprintf(st,"n: k=%02d %9.6f %9.6f %9.6f",k,fc[k].n.x,fc[k].n.y,fc[k].n.z);m_lst.AddString(st);
		for(i=0;i<3;i  )
		{
			sprintf(st,"i=%2d %9.3f %9.3f %9.3f",i,fc[k].dot[i].x,fc[k].dot[i].y,fc[k].dot[i].z);m_lst.AddString(st);
		}


	}
	


	delete [] fc;

}

FXYZ normal(FXYZ v0,FXYZ v1){
	float wk;
	FXYZ nrl;
	nrl.x=v0.y*v1.z - v1.y*v0.z;
	nrl.y=v0.z*v1.x - v0.x*v1.z;
	nrl.z=v0.x*v1.y - v1.x*v0.y;
	wk=float(sqrt(nrl.x*nrl.x nrl.y*nrl.y nrl.z*nrl.z));
	nrl.x=nrl.x/wk;
	nrl.y=nrl.y/wk;
	nrl.z=nrl.z/wk;
	return nrl;
}