計算三維網格的表面積

我有一個由網格和三角形定義的三維網格。我也有網格的法線。我想計算網格的面積，假設它總是封閉的。我發現一個有趣的實現3D體積的計算中this question，而且我申請它在C代碼構建一個由R.此調用的函數是代碼：計算三維網格的表面積

double SignedVolumeOfTriangle(double p1X, double p1Y, double p1Z, 
     double p2X, double p2Y, double p2Z, double p3X, double p3Y, double p3Z) { 
    double v321 = p3X*p2Y*p1Z; 
    double v231 = p2X*p3Y*p1Z; 
    double v312 = p3X*p1Y*p2Z; 
    double v132 = p1X*p3Y*p2Z; 
    double v213 = p2X*p1Y*p3Z; 
    double v123 = p1X*p2Y*p3Z; 
    return (double)(1.0/6.0)*(-v321 + v231 + v312 - v132 - v213 + v123); 
} 
void MeshVolume(double *X, double *Y, double *Z, int *numT, int *V1, int *V2, int *V3, double *Volume) { 
    int n;   
    *Volume=0;  
    for (n=0; n<*numT; n++) { 
     *Volume = *Volume + SignedVolumeOfTriangle(X[V1[n]], Y[V1[n]], Z[V1[n]], X[V2[n]], Y[V2[n]], Z[V2[n]], X[V3[n]], Y[V3[n]], Z[V3[n]]);  
    } 
    *Volume = fabs(*Volume); 
}

既不是問題也不在文章中我找到了計算網格面積的算法。有沒有人可以幫助我？

來源

2014-10-11 Nicola Dinapoli

你是否諮詢了其他人的消息？比如任何一種「圖書館」...... – alk 2014-10-11 08:36:09

不，我沒有。但是我想以最有效的方式解決問題，並且可能不使用外部庫。 – 2014-10-11 08:43:51

;-)有了「*圖書館*」，我的意思是他們收集知識的地方之一，例如在過去堆積書籍...... – alk 2014-10-11 08:50:15

你有一個封閉的體積，其表面由三角形組成。所有的三角形都對外表面起作用。對？

A = 0.5 * |PQ × PR| 
    = 0.5 * |PQ| * |PR| * sin(Ɵ)

其中

PQ = Q - P 
PR = R - P

和×表示cross product和Ɵ是向量之間的角度：

點P，Q和R之間的三角形的表面可以通過以下方式獲得。（交叉積的向量的大小是兩個原始向量之間的平行四邊形的面積，其中一半是三角形的面積）。

總和所有三角形的aeras。沒有必要取絕對值，因爲該區域只能是零或正值。所以：

double AreaOfTriangle(double p1X, double p1Y, double p1Z, 
     double p2X, double p2Y, double p2Z, 
     double p3X, double p3Y, double p3Z) 
{ 
    double ax = p2x - p1x; 
    double ay = p2y - p1y; 
    double az = p2z - p1z; 
    double bx = p3x - p1x; 
    double by = p3y - p1y; 
    double bz = p3z - p1z; 
    double cx = ay*bz - az*by; 
    double cy = az*bx - ax*bz; 
    double cz = ax*by - ay*bx; 

    return 0.5 * sqrt(cx*cx + cy*cy + cz*cz); 
}  

void MeshSurface(double *X, double *Y, double *Z, 
    int *numT, int *V1, int *V2, int *V3, double *Area) 
{ 
    int n; 

    *Area = 0.0; 

    for (n=0; n<*numT; n++) { 
     *Area += AreaOfTriangle(X[V1[n]], Y[V1[n]], Z[V1[n]], 
      X[V2[n]], Y[V2[n]], Z[V2[n]], 
      X[V3[n]], Y[V3[n]], Z[V3[n]]);  
    } 
}

來源

2014-10-11 10:25:35

謝謝M厄姆，我還沒有投票，但這是我正在尋找！ – 2014-10-11 12:27:58

計算三維網格的表面積

回答

相關問題