TIN.sml

# TIN.SML
# Sample script for tutorial "Writing Scripts with SML".
clear();
PopupMessage("Select /LITEDATA/SURFMODL/SURFACE.RVC/ELEV_PTS");
# Get input vector with elevation points and output TIN object.
vector VecPts;
tin TINout;
GetInputVector(VecPts);
GetOutputTIN(TINout,"TINToolkit");	# initialize TIN for editing.
# Read number of vector points and get
# vector georeference information.
numeric numpts;
numpts = NumVectorPoints(VecPts);
class Georef georefV;
georefV = GetLastUsedGeorefObject(VecPts);
# Create implied georeference for TIN
class Georef georefT;
CreateImpliedGeoref(TINout,georefV.Projection);
georefT = GetLastUsedGeorefObject(TINout);
# Create arrays to hold coordinates of vector points.
array numeric xpoints [numpts];
array numeric ypoints [numpts];
array numeric zpoints [numpts];
# Loop to read point coordinates and store in arrays.
numeric i, x, y, z;
numeric mapX, mapY;
for i = 1 to numpts {
	x = VecPts.point[i].Internal.x;
	y = VecPts.point[i].Internal.y;
	z = VecPts.point[i].Internal.z;
	# Convert from object to map coordinates;
	ObjectToMap(VecPts,x,y,georefV,mapX,mapY);
	xpoints[i] = mapX;
	ypoints[i] = mapY;
	zpoints[i] = z;
}
# Create TIN object from coordinate arrays.
numeric zscale, zoffset, xytolerance;
zscale = 1;
zoffset = 0;
xytolerance = 1; # minimum allowed distance between nodes.
TINCreateFromNodes(TINout,numpts,xpoints,ypoints,zpoints,
			zscale,zoffset,xytolerance);
# Get number of TIN hulls, edges, and triangles and print to console.
numeric numhulls, numedges, numtri;
numhulls = TINNumberHulls(TINout);
numedges = TINNumberEdges(TINout);
numtri = TINNumberTriangles(TINout);
print("Number of TIN hulls = ", numhulls);
print("Number of TIN edges = ", numedges);
print("Number of TIN triangles = ", numtri);
CloseTIN(TINout);
# End of script.