- 1、本文档共15页,可阅读全部内容。
- 2、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
多面体面积与体积.
海伦公式:
double a,b,c,p;
/*double m,n,e,p,q,f;*/
ZCoord3D Min,Max,point;
VCoords3D::iterator iter = Coords.begin();
/*Min.x=(*iter).x;*/
Min.y=(*iter).y;
Max.x=(*iter).x;
Max.y=(*iter).y;
Min.z=(*iter).z;
Max.z=(*iter).z;
if((*iter).xMin.x (*iter).xMax.x || (*iter).yMin.y (*iter).yMax.y || (*iter).zMin.z (*iter).zMax.z)
point = *iter;
for(iter += 2; iter != Coords.end(); iter++)
{
/*m = (*(iter - 1)).x -point.x;
n = (*(iter - 1)).y - point.y;
e = (*(iter - 1)).z - point.z;
p = (*iter).x - point.x;
q = (*iter).y - point.y;
f = (*iter).z - point.z;
Area = 0.5 * sqrt((m * q - n * p) * (m * q - n * p)
+ (m * f - e * p) * (m * f - e * p)
+ (n * f - e * q) * (n * f-e * q));*/
a = sqrt(((iter - 1)-x - point.x) * ((iter - 1)-x - point.x)
+ ((iter - 1)-y - point.y) * ((iter - 1)-y - point.y)
+ ((iter - 1)-z - point.z) * ((iter - 1)-z - point.z));
c = sqrt((iter-x - (iter - 1)-x) * (iter-x - (iter - 1)-x)
+ (iter-y - (iter -1 )-y) * (iter-y - (iter - 1)-y)
+ (iter-z - (iter - 1)-z) * (iter-z - (iter - 1)-z));
b = sqrt((iter-x - point.x) * (iter-x - point.x)
+ (iter-y - point.y) * (iter-y - point.y)
+ (iter-z - point.z) * (iter-z - point.z));
p = 0.5 * (a + b + c);
Area = sqrt(p * (p - a) * (p - b) * (p - c));
s += Area;
只能算出凸多边形的面积,不能算出凹多边形面积。
网格法:
double Area = 0;
double Area1 = 0;
double Area2 = 0;
double f1,f2;
double SArea = 0;
ZCoord3D p,pFront,pRight,pDiagonal;
double A[2],B[2],C[2],D[2];
double AB,BC,CD,AC,BD;
ZCoord2D Min,Max;
//算点组成的面投影到平面上的面积
VCoords3D::iterator iter = Coords.begin();
Min.x = (*iter).x;
Min.y = (*iter).y;
Max.x = (*iter).x;
Max.y = (*iter).y;
for(; iter != Coords.end(); ++iter)
{
if((*iter).x Min.x)
Min.x = (*iter).x;
if((*iter).y Min.y)
Min.y = (*iter).y;
if((*iter).x Max.x)
Max.x = (*iter).x;
if((*iter).y Max.y)
Max.y = (*iter).y;
}
int Numx = Max.x - Min.x;
int Numy = Max.y - Min.y;
for (int i = 0;
文档评论(0)