土方量计算，或者叫填挖方计算，体积计算，Skyline在很早的版本中就提供了这个的功能。
目前的软件版本，不仅仅可以对地形修改对象进行土方量计算，还可以在FLY工程中导入DEM数字高程模型数据，计算不同DEM数据之间的土方量变化情况。
这里摘取了针对地形修改对象的土方量计算代码，请大家鉴赏。
<!DOCTYPE html>
<html xmlns="http://www.w3.org/1999/xhtml">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8" />
<title>剖切指定位置的三维模型</title>
<script type="text/javascript">
// ------------------------------------------------------
// 若加载其他工程文件，可以修改下面一行中fly文件地址
var flyurl = "E:\\2018\\DEMO\\Default0726.fly";
// ------------------------------------------------------
function LoadFly() {
var sg = CreateSGObj();
sg.Open(flyurl);
}
window.setTimeout("LoadFly()", 1000); function StartCal() {
var SGWorld = CreateSGObj();
var oid = SGWorld.ProjectTree.FindItem("New Modify Terrain ##5175551");
CalculateList(oid, null);
} function CalculateList(value, geometry) {
var SGWorld = CreateSGObj();
var volumeGeometry;
var ModifyObj = SGWorld.ProjectTree.GetObject(value);
if (geometry == null)
volumeGeometry = ModifyObj.Geometry;
else
volumeGeometry = geometry; var origVisability = SGWorld.ProjectTree.GetVisibility(value); SGWorld.ProjectTree.SetVisibility(value, false);
calculateVolume(volumeGeometry, 2, 3, 0, ModifyObj.TreeItem.Name); SGWorld.ProjectTree.SetVisibility(value, true);
calculateVolume(volumeGeometry, 2, 3, 1, ModifyObj.TreeItem.Name); SGWorld.ProjectTree.SetVisibility(value, origVisability);
} var gAreaUnitFactor = 1;
var gVolumeUnitFactor = 1;
var gAreaUnit = "";
var gVolumeUnit = "";
var gQuaryResolution = 1024;
var gVolumeResultHTML = "";
//-------------
// calculateVolume
function calculateVolume(geometry, type, altitudeType, phase, ObjName) {
var SGWorld = CreateSGObj();
//debugger
//var htmlStr;
//htmlStr = SGLang.i18n("PrepareClaculation");
//SGWorld.Window.ShowMessageBarText(htmlStr, 1, -1); //var qualityStep = ResolutionsArray[$("#DensityID").val()];
var qualityStep = 10; var polygonGeometry = geometry.Clone(); // to enable edit geometry (ref) when geometry type is RelativeToPivot. Navon.
var envelope = polygonGeometry.Envelope; // Multi polygon evnelope (min/max)
var MinX = Math.min(envelope.Rings(0).Points(0).x, envelope.Rings(0).Points(2).x);
var MaxX = Math.max(envelope.Rings(0).Points(0).x, envelope.Rings(0).Points(2).x);
var MinY = Math.min(envelope.Rings(0).Points(0).y, envelope.Rings(0).Points(2).y);
var MaxY = Math.max(envelope.Rings(0).Points(0).y, envelope.Rings(0).Points(2).y); var stepX = Math.abs(MaxX - MinX) / gQuaryResolution;
var stepY = Math.abs(MaxY - MinY) / gQuaryResolution; // calculate cell 2D size
var cellPos0 = SGWorld.Creator.CreatePosition(MinX, MaxY, 0, 3, 0, 0, 0, 1);
var cellPos1 = SGWorld.Creator.CreatePosition(MinX, MaxY - stepY, 0, 3, 0, 0, 0, 1);
var cellPos2 = SGWorld.Creator.CreatePosition(MinX + stepX, MaxY, 0, 3, 0, 0, 0, 1);
var cellWidth = cellPos0.DistanceTo(cellPos1) * qualityStep;
var cellArea = cellPos0.DistanceTo(cellPos1) * cellPos0.DistanceTo(cellPos2) * qualityStep * qualityStep * gAreaUnitFactor; if (phase == 0) {
// For now we call this query several times because of an engine bug!!!
gElevationBufferBefore = SGWorld.Analysis.QueryElevationBuffer(MinX, MaxY, stepX, stepY, gQuaryResolution, gQuaryResolution).toArray(); // MaxX - MinX, MaxY - MinY);
gElevationBufferBefore = SGWorld.Analysis.QueryElevationBuffer(MinX, MaxY, stepX, stepY, gQuaryResolution, gQuaryResolution).toArray(); // MaxX - MinX, MaxY - MinY);
gElevationBufferBefore = SGWorld.Analysis.QueryElevationBuffer(MinX, MaxY, stepX, stepY, gQuaryResolution, gQuaryResolution).toArray(); // MaxX - MinX, MaxY - MinY);
}
else {
// For now we call this query several times because of an engine bug!!!
gElevationBufferAfter = SGWorld.Analysis.QueryElevationBuffer(MinX, MaxY, stepX, stepY, gQuaryResolution, gQuaryResolution).toArray(); // MaxX - MinX, MaxY - MinY);
gElevationBufferAfter = SGWorld.Analysis.QueryElevationBuffer(MinX, MaxY, stepX, stepY, gQuaryResolution, gQuaryResolution).toArray(); // MaxX - MinX, MaxY - MinY);
gElevationBufferAfter = SGWorld.Analysis.QueryElevationBuffer(MinX, MaxY, stepX, stepY, gQuaryResolution, gQuaryResolution).toArray(); // MaxX - MinX, MaxY - MinY); //// var group;
//var drawObjects = $("#showObjectsID").is(':checked') ? true : false;
//if (drawObjects) {
// CreateLayer(cellWidth * 100, ObjName);
// var result = featureLayerStyles["ImageLabel"](cellWidth / 40, abspath() + "/img/point.png"); // ,"[Color]");
//}
//group = SGWorld.ProjectTree.CreateGroup( SGLang.i18n("ToolName"), ""); var VolumeAdded = 0;
var VolumeRemoved = 0;
for (i = qualityStep / 2; i < gQuaryResolution; i = i + qualityStep) {
for (j = qualityStep / 2; j < gQuaryResolution; j = j + qualityStep) {
var x = MinX + i * stepX;
var y = MaxY - j * stepY;
var objColor = "#000000";
var volumeType = 0;
var pointGeometry = SGWorld.Creator.GeometryCreator.CreatePointGeometry([x, y, 0]);
if (polygonGeometry.SpatialRelation.Intersects(pointGeometry)) {
var elevationDiff = gElevationBufferAfter[j * gQuaryResolution + i] - gElevationBufferBefore[j * gQuaryResolution + i];
var volumeDiff = cellArea * elevationDiff * gVolumeUnitFactor;
if (elevationDiff > 0) { // Added
volumeType = 1;
VolumeAdded += Math.abs(volumeDiff);
objColor = 65280; // "#00ff00"; // green
}
else if (elevationDiff < 0) { // Removed
volumeType = 2;
VolumeRemoved += Math.abs(volumeDiff);
objColor = 255; // "#ff0000"; // red
}
else {
volumeType = 0;
objColor = 0; // "#000000"; // black
}
//if (drawObjects) { // && i % ObjInterval == 0 && j % ObjInterval == 0) {
// var pos = SGWorld.Creator.CreatePosition(x, y, gElevationBufferAfter[j * gQuaryResolution + i], 3, 0, 0, 0, 1);
// gLayer.FeatureGroups.Point.CreateFeature([pos.X, pos.Y, pos.Altitude], volumeType + ";" + cellArea.toFixed(3) + ";" + volumeDiff.toFixed(3) + ";" + objColor);
// //SGWorld.Creator.CreateSphere(pos, 2, 0, objColor, objColor, 2, group, "sphere");
//}
} // if intersect
} // for j
if (i % 10 == 0) {
var progress = (i / gQuaryResolution) * 100;
htmlStr = SGLang.i18n("calculatingVolume") + ObjName + ": " + progress.toFixed(0) + "%";
SGWorld.Window.ShowMessageBarText(htmlStr, 1, -1);
}
} // for i
SGWorld.Window.HideMessageBarText();
gVolumeResultHTML += "<br/><u>" + ObjName + ":</u><br/>" + "objectVolumeAdded" + "<b>" + VolumeAdded.toFixed(3) + "</b> " + gVolumeUnit + " <br/> " + "objectVolumeRemoved" + "<b>" + VolumeRemoved.toFixed(3) + "</b> " + gVolumeUnit + "<br/>";
//gTotalVolumeAdded += VolumeAdded;
//gTotalVolumeRemoved += VolumeRemoved;
alert(gVolumeResultHTML);
//if (drawObjects)
// gLayer.Save();
} // phase return true;
} // SGWorld70
function CreateSGObj() {
try {
var obj = document.getElementById("SGWorld");
if (obj == null) {
obj = document.createElement('object');
obj.setAttribute("name", "SGWorld");
obj.setAttribute("id", "SGWorld");
obj.style.height = "1px";
obj.style.width = "1px";
obj.setAttribute("classid", "CLSID:3A4F919C-65A8-11D5-85C1-0001023952C1");
document.body.appendChild(obj);
}
return obj;
}
catch (e) {
alert(e);
}
} </script>
</head>
<body style="margin: 0px; border: 0px;">
<div id="top">
<div>
<input id="Button1" type="button" value="分析" style="position:absolute;top:6px;right:100px;" onclick="StartCal()" /> </div>
</div>
<div style="position: absolute; top: 50px; bottom: 0px; width: 100%; background-color: black; margin: 0px; border: 0px; margin-top: 3px;">
<object id="TerraExplorerInformationWindow" classid="CLSID:3a4f9193-65a8-11d5-85c1-0001023952c1" style="margin-left: 3px; width: 19%; height: 99%;"></object>
<object id="TerraExplorer3DWindow" classid="CLSID:3a4f9192-65a8-11d5-85c1-0001023952c1" style="width: 80%; height: 99%;"></object>
</div>
</body> </html>