three.js/examples/js/VTKLoader.js

/**
 * @author Sebastien Valette [email protected]
 */

THREE.VTKLoader = function ( showStatus ) {

	THREE.Loader.call( this, showStatus );

};

THREE.VTKLoader.prototype = new THREE.Loader();
THREE.VTKLoader.prototype.constructor = THREE.VTKLoader;
THREE.VTKLoader.prototype.supr = THREE.Loader.prototype;

THREE.VTKLoader.prototype.load = function( url, callback) {

	this.onLoadStart();

	function  parseGeometry(text){
		var geometry=new THREE.Geometry();
		var lines=text.split("\n");

		function v( x, y, z ) {

			geometry.vertices.push( new THREE.Vector3( x, y, z ) );

		}

		function f3( a, b, c ) {

			geometry.faces.push( new THREE.Face3( a, b, c ) );

		}

		var lineIndex=0;

		var line=lines[0].split(" ");
		var lineLength=line.length;
		var columnIndex=-1;


		function readNextString ()
		{
			while (1)
			{
				var nextWord=line[columnIndex];
				columnIndex++;
				if (columnIndex==lineLength)
				{
					lineIndex++;
					columnIndex=0;
					if (lineIndex>lines.length)
						return ("");
					line=lines[lineIndex].split(" ");
					lineLength=line.length;
				}
				if (nextWord!=null)
					if (nextWord.length>0)
						return (nextWord);
			}
		}

		// read point data
		var found=false;
		while (!found)
		{
			var readString=readNextString();
			switch (readString.toUpperCase())
			{
				case "POINTS":
					found=true;
					break;
				case "":
					alert ("error while reading "+url+" : \n");
					return;
				default:
			}

		}

		var newIndex;
		var new2old;

		var numberOfPoints=parseInt(readNextString());

		if (numberOfPoints>5000000)
		{
			alert ("mesh is too big : "+numberOfPoints+" vertices");
			return;
		}

		var coord=[0,0,0];
		var index2=0;

		var number;
		var coordIndex;
		for (var j=0;j!=numberOfPoints;j++)
		{
			for (coordIndex=0;coordIndex<3;coordIndex++)
			{
				do
				{
					number=parseFloat(line[columnIndex]);
					columnIndex++;
					if (columnIndex==lineLength)
					{
						lineIndex++;
						columnIndex=0;
						if (lineIndex>lines.length)
						{
							alert ("error while reading "+url+" : \n");
							return;
						}
						line=lines[lineIndex].split(" ");
						lineLength=line.length;
					}
				} while (isNaN(number))
				coord[coordIndex]=number;
			}
			v(coord[0],coord[1],coord[2]);
		}

		found=false;
		while (!found)
		{
			var readString=readNextString();
			switch (readString)
			{
				case "POLYGONS":
					found=true;
					break;
				case "":
					alert ("error while reading "+url+" : \n");
					return;
				default:
			}
		}

		var numberOfPolygons=parseInt(readNextString());
		var numberOfpolygonElements=parseInt(readNextString());

		index2=0;
		var connectivity=[];
		for (var p=0;p!=numberOfpolygonElements;p++)
		{
			do
			{
				number=parseInt(line[columnIndex]);
				columnIndex++;
				if (columnIndex==lineLength)
				{
					lineIndex++;
					columnIndex=0;
					if (lineIndex>lines.length)
					{
						alert ("error while reading "+url+" : \n");
						return;
					}
					line=lines[lineIndex].split(" ");
					lineLength=line.length;
				}
			} while (isNaN(number))

			connectivity[index2]=number;
			index2++;

			if (index2==connectivity[0]+1)
			{
				var triangle=[0,0,0];
				index2=0;
				var numberOfTrianglesInCell=connectivity[0]-2;
				var vertex1=triangle[0]=connectivity[1];

				for (var i=0;i<numberOfTrianglesInCell;i++)
				{
					var vertex2=connectivity[i+2];
					var vertex3=triangle[2]=connectivity[i+3];
					f3(vertex1,vertex2,vertex3);
				}
			}
		}

		geometry.computeCentroids();
		geometry.computeFaceNormals();
		geometry.computeVertexNormals();
		geometry.computeBoundingSphere();
		return (geometry);
	}

	var xmlhttp=new XMLHttpRequest();
	xmlhttp.open("GET",url,true);

	function handler()
	{
		if(xmlhttp.readyState == 4)
		{
			if (xmlhttp.status!=200)
			{
				alert("Could not read file "+url+": error"+xmlhttp.status);
				return;
			}
			var myGeometry=parseGeometry(xmlhttp.responseText);
			callback(myGeometry);
		}
	}
	xmlhttp.onreadystatechange=handler;
	xmlhttp.send();

};