请教：如何根据实测的xyz坐标求解土方体积？

hanqinhu

已经有数千个数据，但是与matlab常用的矩阵格式不同，这些数据表示的是土坡表面的三维坐标，即其x，y，z 的坐标，例如下面这些数据，三列分别表示不同点的xyz的坐标，这些点构成了一个土坡的表面。

4.793429e+005,3.942442e+006,2.835362e+001
4.793529e+005,3.942442e+006,2.819601e+001
4.793629e+005,3.942442e+006,2.811963e+001
4.793729e+005,3.942442e+006,2.799952e+001
4.793829e+005,3.942442e+006,2.767889e+001
4.793929e+005,3.942442e+006,2.739609e+001

问题是：如何根据这些数据求解在某个水平面以上土体的方量呢？


其实这些数据本可以导入surfer计算，但是自己想用matlab编一个程序，但是试了好多次却没能成功，上来请教一下诸位，还请不吝赐教！谢谢了！

lengyunfeng

你是用在建筑整平上还是用在其它上面？如果建筑或者线型构筑物的整平不仅仅涉及某个水平面，还涉及放坡面，用matlab来编的话比较麻烦。即便是你只要某个水平面以上的土方面积，这里面还会涉及某些地面高程点位于水平面以下，土方零线并不是连通，进而计算方法也会有很多差异。建议你还是用专业一点的软件来计算，用matlab编的没办法验证其准确性，即便编出来了也不敢放心的用的。有一个软件叫htcad你不妨试一下。

hanqinhu

呵呵，谢谢lengyunfeng了，是用于整平的，水平面，其实htcad和cass都有，其实我用surfer计算得到的结果和cass误差也极小，单纯从工作角度来说交差足够了，但是就想自己试一下用matlab行不行，纯粹个人爱好，只是水平有限，心有余而力不足，见笑了。

lengyunfeng

那你不妨把代码帖出来，大家研究一下说不定还真能帮到你。
3# hanqinhu