<?php
# $Id: class_gml2.php 3099 2008-10-02 15:29:23Z nimix $
# http://www.mapbender.org/index.php/class_gml2.php
# Copyright (C) 2002 CCGIS 
#
# This program is free software; you can redistribute it and/or modify
# it under the terms of the GNU General Public License as published by
# the Free Software Foundation; either version 2, or (at your option)
# any later version.
#
# This program is distributed in the hope that it will be useful,
# but WITHOUT ANY WARRANTY; without even the implied warranty of
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
# GNU General Public License for more details.
#
# You should have received a copy of the GNU General Public License
# along with this program; if not, write to the Free Software
# Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA.

require_once(dirname(__FILE__)."/../../core/globalSettings.php");
require_once(dirname(__FILE__)."/../classes/class_connector.php");
require_once(dirname(__FILE__)."/../classes/class_json.php");

class gml2 {
	var $geomtype_point = 'Point';					
	var $geomtype_polygon = 'Polygon';
	var $geomtype_line = 'LineString';
	var $geomtype_multipolygon = 'MultiPolygon'; 
	var $geomtype_multiline = 'MultiLine';
	var $geometries = array();	
	var $member = -1;
	var $geomtype = array();	
	var $keys = array();
	var $values = array();
	var $geometry = array();
	var $bbox = array();
	var $doc;
	var $geomFeaturetypeElement = null;
	
	
	function gml2(){
		$this->geometries = array($this->geomtype_point, $this->geomtype_polygon, $this->geomtype_line, $this->geomtype_multipolygon, $this->geomtype_multiline);
	}
	function getGml($req){
		$x = new connector($req);
		return $x->file;
	}

	function parseFile($req){
		#$data = implode("",file($req));
		$x = new connector($req);
		$data = $x->file;
		$data = $this->removeWhiteSpace($data);
		$this->parseXML($data);		
		#$e = new mb_exception("data = ".$data); 		
	}
	
	function parseXML($data) {
		if (func_num_args() == 2) {
			$this->geomFeaturetypeElement = func_get_arg(1);
		}
		$this->doc = $this->removeWhiteSpace($data);
		return $this->toGeoJSON();
	}

	function removeWhiteSpace ($string) {
		return preg_replace("/\>(\s)+\</", "><", trim($string));
	}
	
	function sepNameSpace($s){
		$c = mb_strpos($s,":"); 
		if($c>0){
			return mb_substr($s,$c+1);
		}
		else{
			return $s;
		}		
	}
	
	function toGeoJSON () {
		$gmlDoc = new SimpleXMLElement($this->doc);
		
		$gmlDoc->registerXPathNamespace('xls', 'http://www.opengis.net/xls');
		$gmlDoc->registerXPathNamespace('wfs', 'http://www.opengis.net/wfs');
		$gmlDoc->registerXPathNamespace('gml', 'http://www.opengis.net/gml');
		
		// build feature collection
		$featureCollection = new FeatureCollection();
		
		// segments of the featzreCollection
		$gmlFeatureMembers = $gmlDoc->xpath("//gml:featureMember");
		
		if(count($gmlFeatureMembers)>0){
			$cnt=0;
			foreach ($gmlFeatureMembers as $gmlFeatureMember) {
				$featureMember_dom = dom_import_simplexml($gmlFeatureMember);
				
				$feature = new Feature();
				if ($this->geomFeaturetypeElement != null) {
					$feature->parse($featureMember_dom, $this->geomFeaturetypeElement);
				}
				else {
					$feature->parse($featureMember_dom);
				}
				if (isset($feature->geometry)) {
					$featureCollection->addFeature($feature);
				}
				$cnt++;
			}
			
			return $featureCollection->toGeoJSON();
		}
		else{
			return "{}";
		}
	}

	
	/**
	 * Exports the file to SHP.
	 * 
	 * @param string $filenamePrefix the filename without an ending like .shp
	 */
	function toShape ($filenamePrefix) {
		$unique = TMPDIR . "/" . $filenamePrefix;
		$fShape = $unique.".shp";
		$fGml = $unique.".gml";
//		$pathOgr = '/appserver/postgresql/templates/postgresql824/bin/ogr2ogr';
		$pathOgr = '/usr/bin/ogr2ogr';
		$w = $this->toFile($fGml);

		/*
		 * @security_patch exec done
		 */
 		$exec = $pathOgr.' -f "ESRI Shapefile" "'.$fShape.'" '.$fGml;
 		exec(escapeshellcmd($exec));
 		
 		$exec = 'zip -j '.$unique.' '.$unique.'.shp '.$unique.'.dbf '.$unique.'.shx '.$unique.'.gfs '.$unique.'.gml ';
 		exec(escapeshellcmd($exec));

		$exec = 'rm -f '.$unique.' '.$unique.'.shp '.$unique.'.dbf '.$unique.'.shx '.$unique.'.gfs '.$unique.'.gml';
		exec(escapeshellcmd($exec));
		
		$exec = 'chmod 777 '.$unique.'.*';
		exec(escapeshellcmd($exec));
		//echo "<a href='../tmp/".$unique.".zip'>Download ".$prefix."<a>";
	}
	
	/**
	 * Writes a file containing the GML.
	 * 
	 * @param string $path the path to the file.
	 * @param string $path the path to the file.
	 * @return bool true if the file could be saved; else false.
	 */
	function toFile ($file) {
		$handle = fopen($file, "w");
		if (!$handle) {
			$e = new mb_exception("class_gml2.php: Filehandler error (" . $file . ").");
			return false;
		}
		if (!fwrite($handle,$this->__toString())) {
			$e = new mb_exception("class_gml2.php: Could not write file (" . $file . ").");
			fclose($handle);
			return false;
		}
		fclose($handle);
		return true;
	}
	
	function __toString () {
		return $this->doc;
	}


	/**
	 * @deprecated
	 */
	function parsegml($req){
		#$data = implode("",file($req));
		$x = new connector($req);
		$this->parse_xml($x->file);		 		
	}

	/**
	 * @deprecated
	 */
	function parse_xml($data){


		$this->doc = $data;

		$section = false;
		$geom = false;
		$boundedBy = false;
		$coordinates = false;
		$el = -1;
		$parser = xml_parser_create(CHARSET);
		xml_parser_set_option($parser,XML_OPTION_CASE_FOLDING,0);
		xml_parser_set_option($parser,XML_OPTION_SKIP_WHITE,1);
		xml_parser_set_option($parser,XML_OPTION_TARGET_ENCODING,CHARSET);
		xml_parse_into_struct($parser,$data,$values,$tags);
		$code = xml_get_error_code ($parser);
		if ($code) {
			$line = xml_get_current_line_number($parser); 
			$mb_exception = new mb_exception(xml_error_string($code) .  " in line " . $line);
		}
		xml_parser_free($parser);
		
		foreach ($values as $element) {
			#$e = new mb_exception($element[tag]);
			if(strtoupper($this->sepNameSpace($element[tag])) == strtoupper("boundedBy") && $element[type] == "open"){
				$boundedBy = true;
			}
			if ($boundedBy) {
				if(strtoupper($this->sepNameSpace($element[tag])) == strtoupper("coordinates")){
					$this->bbox = explode(",", str_replace(",,","",str_replace(" ",",",trim($element[value]))));
					$boundedBy=false;
				}
			}
			if(strtoupper($this->sepNameSpace($element[tag])) == strtoupper("featureMember") && $element[type] == "open"){
				$this->member++;
				$this->keys[$this->member] = array();
				$this->value[$this->member] = array();
				$this->geometry[$this->member] = array();
				$section = true;
				$cnt_geom = 0;
			}
			if($section == true){
				if( in_array($this->sepNameSpace($element[tag]),$this->geometries) && $element[type] == "open"){
					$geom = true;
					$this->geomtype[$this->member] = $this->sepNameSpace($element[tag]);
				}
				else if(in_array($this->sepNameSpace($element[tag]),$this->geometries) && $element[type] == "close"){
					$cnt_geom++;
					$geom = false;
				}
				if($geom == true){
					if (strtoupper($this->sepNameSpace($element[tag])) == strtoupper("coordinates")) {
						$this->geometry[$this->member][$cnt_geom] =  str_replace(",,","",str_replace(" ",",",trim($element[value])));
						$coordinates = true;
						// XXX: Increment counter to get all geometries of a feature member, 
						// comment it out to only show first geometry of featuremember
						$cnt_geom++;
					}
					else if (!$coordinates && trim($element[value])) {
						$coords = str_replace(",,","",str_replace(" ",",",trim($element[value])));
						$tmp_array = explode(",", $coords);
						if (count($tmp_array > 1)) {
							$this->geometry[$this->member][$cnt_geom] =  $coords;
							$coordinates = true;
						   // XXX: Increment counter to get all geometries of a feature member, 
						   // comment it out to only show first geometry of featuremember
							$cnt_geom++;
						}
					}
				}
				else if(strtoupper($this->sepNameSpace($element[tag])) == strtoupper("featureMember") && $element[type] == "close"){
					$section = false;	
					$el = -1;
				}
				else{
					$el++;
					$this->values[$this->member][$el] = $element[value];
					$this->keys[$this->member][$el] = $element[tag];	
				}
			}
		}
	}	

	/**
	 * @deprecated
	 */
	function getMemberCount(){
		return ($this->member+1);	
	}
	
	/**
	 * @deprecated
	 */
	function getValueBySeparatedKey($memberCount,$keyName){
		for($i=0; $i<count($this->keys[$memberCount]); $i++){
			if($this->sepNameSpace($this->keys[$memberCount][$i]) == $keyName){
				return $this->values[$memberCount][$i];
			}	
		}
	}
	
	/**
	 * @deprecated
	 */
	function getValueByKey($memberCount,$keyName){
		for($i=0; $i<count($this->keys[$memberCount]); $i++){
			if($this->keys[$memberCount][$i] == $keyName){
				return $this->values[$memberCount][$i];
			}	
		}
	}
	
	/**
	 * @deprecated
	 */
	function getXfromMemberAsString($memberCount,$geomCount){
		$t = explode(",",$this->geometry[$memberCount][$geomCount]);
		$x = array();
		for($i=0; $i<(count($t)-1); $i=$i+2){
			array_push($x,$t[$i]);
		}
		return implode(",",$x);
	}
	
	/**
	 * @deprecated
	 */
	function getYfromMemberAsString($memberCount,$geomCount){
		$t = explode(",",$this->geometry[$memberCount][$geomCount]);
		$y = array();
		for($i=1; $i<=(count($t)-1); $i=$i+2){
			array_push($y,$t[$i]);
		}
		return implode(",",$y);
	}
	
	/**
	 * @deprecated
	 */
	function getGeometriesFromMember($memberCount){
		return $this->geometry[$memberCount];	
	}
	
	/**
	 * @deprecated
	 */
	function getGeometryTypeFromMember($memberCount){
		return 	$this->geomtype[$memberCount];
	}

	/**
	 * @deprecated
	 */
	function exportGeometriesToJS ($isInFrame) {
		$prefix = "";
		if ($isInFrame == true) {
			$prefix = "parent.";
		}
		$js = "";
		$js .= "var geom = new ".$prefix."GeometryArray();\n";
		for ($i=0; $i<count($this->geometry); $i++) {
			$js .= $this->exportMemberToJS($i, $isInFrame);
			$js .= "geom.addCopy(q);\n";
		}
		return $js;
	}
	
	/**
	 * @deprecated
	 */
	function exportMemberToJS ($i, $isInFrame) {
		$prefix = "";
		if ($isInFrame == true) {
			$prefix = "parent.";
		}
		$js = "";
		if ($this->getGeometryTypeFromMember($i) == $this->geomtype_point) {
			$js .= "var current_geomtype = ".$prefix."geomType.point;\n";
		}
		elseif ($this->getGeometryTypeFromMember($i) == $this->geomtype_line || $this->getGeometryTypeFromMember($i) == $this->geomtype_multiline) {
			$js .= "var current_geomtype = ".$prefix."geomType.line;\n";
		}
		elseif ($this->getGeometryTypeFromMember($i) == $this->geomtype_polygon || $this->getGeometryTypeFromMember($i) == $this->geomtype_multipolygon) {
			$js .= "var current_geomtype = ".$prefix."geomType.polygon;\n";
		}
		else{
			$e = new mb_notice("unknown geometry type: '".$this->getGeometryTypeFromMember($i)."' or no geometry existing");
			return "";
		}	
		
		$js .= "var q = new ".$prefix."MultiGeometry(current_geomtype);\n";
		
		for ($j=0; $j<count($this->geometry[$i]); $j++) {
			$js .= "q.addGeometry(current_geomtype);\n";
			
			$x_array = explode(",", $this->getXfromMemberAsString($i, $j));
			$y_array = explode(",", $this->getYfromMemberAsString($i, $j));
			
			for ($k=0; $k<count($x_array); $k++) {
				$js .= "q.get(-1).addPointByCoordinates(parseFloat(".$x_array[$k]."), parseFloat(".$y_array[$k]."));\n";
				//$js .= "alert(parseFloat(".$x_array[$k]."), parseFloat(".$y_array[$k]."));";
			}
			$js .= "q.get(-1).close();\n";
		}
//		$js .= "alert(q);\n";
		return $js;
	}

}


class FeatureCollection {
	var $type = "FeatureCollection";
	var $featureArray = array();
	
	public function __construct() {
		
	}
	
	public function addFeature ($aFeature) {
		array_push($this->featureArray, $aFeature);
	}
	
	public function toGeoJSON () {
		$str = "";
		$len = count($this->featureArray); 
		if ($len > 0) {
			$str .= "{\"type\": \"FeatureCollection\", \"features\": [";
			for ($i=0; $i < $len; $i++) {
				if ($i > 0) {
					$str .= ",";
				}	
				$str .= $this->featureArray[$i]->toGeoJSON();
			}
			$str .= "]}";
		}
		return $str;
	}
}

class Feature {
	var $type = "Feature";
	var $fid;
	var $geometry = false;
	var $properties = array();
	var $geomFeaturetypeElement = null;
	
	public function __construct() {
	}
	
	function sepNameSpace($s){
		list($ns,$FeaturePropertyName) = split(":",$s);
		$nodeName = array('ns' => $ns, 'value' => $FeaturePropertyName);
		return $nodeName;
	}
	
	/**
	 * Parses the feature segment of a GML and stores the geometry in the
	 * $geometry variable of the class.
	 * 	
	 * Example of a feature segment of a GML. 
	 * 	<gml:featureMember>
	 * 		<ms:ROUTE fid="ROUTE.228168">
	 * 			<gml:boundedBy>
	 * 				<gml:Box srsName="EPSG:31466">
	 * 					<gml:coordinates>2557381.0,5562371.1 2557653.7,5562526.0</gml:coordinates>
	 * 				</gml:Box>
	 * 			</gml:boundedBy>
	 * 			<ms:geometry>
	 * 				<gml:LineString>
	 * 					<gml:coordinates>
	 * 						2557380.97,5562526 2557390.96,
	 * 						5562523.22 2557404.03,5562518.2 2557422.31,
	 * 						5562512 2557437.16,5562508.37 2557441.79,
	 * 						5562507.49 2557454.31,5562505.1 2557464.27,
	 * 						5562503.97 2557473.24,5562502.97 2557491.67,
	 * 						5562502.12 2557505.65,5562502.43 2557513.78,
	 * 						5562501.12 2557520.89,5562498.79 2557528.5,
	 * 						5562495.07 2557538.9,5562488.91 2557549.5,
	 * 						5562483.83 2557558.55,5562476.61 2557569.07,
	 * 						5562469.82 2557576.61,5562462.72 2557582.75,
	 * 						5562457.92 2557588.57,5562452.56 2557590.38,
	 * 						5562449.69 2557593.57,5562445.07 2557596.17,
	 * 						5562441.31 2557601.71,5562433.93 2557612.97,
	 * 						5562421.03 2557626,5562405.33 2557639.66,
	 * 						5562389.75 2557653.69,5562371.12 
	 * 					</gml:coordinates>
	 * 				</gml:LineString>
	 * 			</ms:geometry>
	 * 			<code>354</code>
	 * 			<Verkehr>0</Verkehr>
	 * 			<rlp>t</rlp>
	 * 		</ms:ROUTE>
	 * 	</gml:featureMember>
	 * 
	 * @return void
	 * @param $domNode DOMNodeObject the feature tag of the GML 
	 * 								(<gml:featureMember> in the above example)
	 */
	public function parse($domNode) {
		if (func_num_args() == 2) {
			$this->geomFeaturetypeElement = func_get_arg(1);
		}

		$currentSibling = $domNode->firstChild;
		
		$this->fid = $currentSibling->getAttribute("fid");
		
		$currentSibling = $currentSibling->firstChild;
		
		while ($currentSibling) {
		
			$name = $currentSibling->nodeName;
			$value = $currentSibling->nodeValue;
			
			$namespace = $this->sepNameSpace($name);
			$ns = $namespace['ns'];
			$columnName = $namespace['value'];
			$isGeomColumn = ($this->geomFeaturetypeElement == null || $columnName == $this->geomFeaturetypeElement);
			
			// check if this node is a geometry node.
			// however, even if it is a property node, 
			// it has a child node, the text node!
			// So we might need to do something more 
			// sophisticated here...
			if ($currentSibling->hasChildNodes() && $isGeomColumn){
				$geomNode = $currentSibling->firstChild; 
					$geomType = $geomNode->nodeName;
					switch ($geomType) {
						case "gml:Polygon" :
							$this->geometry = new GMLPolygon();
							$this->geometry->parsePolygon($geomNode);
							break;
						case "gml:LineString" :
							$this->geometry = new GMLLine();
							$this->geometry->parseLine($geomNode);
							break;
						case "gml:Point" :
							$this->geometry = new GMLPoint();
							$this->geometry->parsePoint($geomNode);
							break;
						case "gml:MultiLineString" :
							$this->geometry = new GMLMultiLine();
							$this->geometry->parseMultiLine($geomNode);
							break;
						case "gml:MultiPolygon" :
							$this->geometry = new GMLMultiPolygon();
							$this->geometry->parseMultiPolygon($geomNode);
							break;
						case "gml:Envelope" :
							$this->geometry = new GMLEnvelope();
							$this->geometry->parseEnvelope($geomNode);
							break;
						default:
							$this->properties[$columnName] = $value;
							break;
					}
			} 
			else {
					$this->properties[$columnName] = $value;
			}
			
			$currentSibling = $currentSibling->nextSibling;
		}
	}
	
	public function toGeoJSON () {
		$str = "";
		$str .= "{\"type\":\"Feature\", \"id\":\"".$this->fid."\", \"geometry\": ";
		if ($this->geometry) {
			$str .= $this->geometry->toGeoJSON();
		}
		else {
			$str .= "\"\"";
		}

		
		$prop = array();
		
		$str .= ", \"properties\": ";
		$cnt = 0;
		foreach ($this->properties as $key => $value) {
				$prop[$key] = $value;
				$cnt ++;
		}

		$json = new Mapbender_JSON();
		$str .= $json->encode($prop); 
		$str .= "}";
		
		return $str;
	}
}

class GMLLine {

	var $pointArray = array();

	public function __construct() {
		
	}

	public function parseLine ($domNode) {
		$currentSibling = $domNode->firstChild;
		while ($currentSibling) {
			
			foreach(explode(' ',$currentSibling->nodeValue) as $cords){
				list($x,$y,$z) = explode(',',$cords);
				$this->addPoint($x, $y);
			}
			$currentSibling = $currentSibling->nextSibling;
		}
	}
	
	protected function addPoint ($x, $y) {
		array_push($this->pointArray, array("x" => $x, "y" => $y));
	}
	
	public function toGeoJSON () {
		$numberOfPoints = count($this->pointArray);
		$str = "";
		if ($numberOfPoints > 0) {
			$str .= "{\"type\": \"LineString\", \"coordinates\":[";
			for ($i=0; $i < $numberOfPoints; $i++) {
				if ($i > 0) {
					$str .= ",";
				}
				$str .= "[".$this->pointArray[$i]["x"].",".$this->pointArray[$i]["y"]."]";
			}
			$str .= "]}";
		}
		else {
			$e = new mb_exception("GMLLine: toGeoJSON: this point is null.");
		}
		return $str;
	}
}

class GMLPoint {

	var $point;

	public function __construct() {
		
	}

	public function parsePoint ($domNode) {
		$currentSibling = $domNode->firstChild;
		while ($currentSibling) {
			list($x, $y, $z) = explode(",", $currentSibling->nodeValue);
			$this->setPoint($x, $y);
			$currentSibling = $currentSibling->nextSibling;
		}
	}
	
	protected function setPoint ($x, $y) {
#		echo "x: " . $x . " y: " . $y . "\n";
		$this->point = array("x" => $x, "y" => $y);
	}
	
	public function toGeoJSON () {
		$str = "";
		if ($this->point) {
			$str .= "{\"type\": \"Point\", \"coordinates\":";
			$str .= "[".$this->point["x"].",".$this->point["y"]."]";
			$str .= "}";
		}
		else {
			$e = new mb_exception("GMLPoint: toGeoJSON: this point is null.");
		}
		return $str;
	}
}

class GMLPolygon {

	var $pointArray = array();
	var $innerRingArray = array();

	public function __construct() {
		
	}

	public function parsePolygon ($domNode) {
		$simpleXMLNode = simplexml_import_dom($domNode);

		$simpleXMLNode->registerXPathNamespace('gml', 'http://www.opengis.net/gml');
		
		$allCoords = $simpleXMLNode->xpath("gml:outerBoundaryIs/gml:LinearRing/gml:coordinates");
			
		$cnt=0;
		foreach ($allCoords as $Coords) {
			$coordsDom = dom_import_simplexml($Coords);
				
//			$name = $coordsDom->nodeName;
//			$value = $coordsDom->nodeValue;				
//			echo "===> name: ".$name. ", Value: ".$value."<br>";
			
			foreach(explode(' ',$coordsDom->nodeValue) as $pointCoords){

				list($x,$y,$z) = explode(',',$pointCoords);
				$this->addPoint($x, $y);
				}
			
			$cnt++;
		}
		
		$innerRingNodeArray = $simpleXMLNode->xpath("gml:innerBoundaryIs/gml:LinearRing");
		if ($innerRingNodeArray) {
			$ringCount = 0;
			foreach ($innerRingNodeArray as $ringNode) {
				$coordinates = $ringNode->xpath("gml:coordinates");
				foreach ($coordinates as $coordinate) {
					$coordsDom = dom_import_simplexml($coordinate);
						
					foreach(explode(' ',$coordsDom->nodeValue) as $pointCoords){
		
						list($x,$y,$z) = explode(',',$pointCoords);
						$this->addPointToRing($ringCount, $x, $y);
					}
				}
				$ringCount++;
			}
		}
	}
	
	protected function addPoint ($x, $y) {
		array_push($this->pointArray, array("x" => $x, "y" => $y));
	}
	
	protected function addPointToRing ($i, $x, $y) {
		if (count($this->innerRingArray) <= $i) {
			array_push($this->innerRingArray, array());
		}
		$index = count($this->innerRingArray);
		$currentIndex = ($i < $index ? $i : $index);
		array_push($this->innerRingArray[$currentIndex], array("x" => $x, "y" => $y));
	}
	
	public function toGeoJSON () {
		$numberOfPoints = count($this->pointArray);
		$str = "";
		if ($numberOfPoints > 0) {
			$str .= "{\"type\": \"Polygon\", \"coordinates\":[[";
			for ($i=0; $i < $numberOfPoints; $i++) {
				if ($i > 0) {
					$str .= ",";
				}
				$str .= "[".$this->pointArray[$i]["x"].",".$this->pointArray[$i]["y"]."]";
			}
			$str .= "]";
			
			for ($i=0; $i < count($this->innerRingArray); $i++) {
				$str .= ",[";
				for ($j=0; $j < count($this->innerRingArray[$i]); $j++) {
					if ($j > 0) {
						$str .= ",";
					}
					$str .= "[".$this->innerRingArray[$i][$j]["x"].",".$this->innerRingArray[$i][$j]["y"]."]";
				}
				$str .= "]";
			}
			$str .= "]}";
		}
		else {
			$e = new mb_exception("GMLPolygon: toGeoJSON: this point is null.");
		}
		return $str;
	}
}

class GMLEnvelope extends GMLPolygon{
/*      <gml:Envelope>
         <gml:lowerCorner>42.943 -71.032</gml:lowerCorner>
         <gml:upperCorner>43.039 -69.856</gml:upperCorner>
      </gml:Envelope>
*/
	public function parseEnvelope ($domNode) {
		$corner1 = $domNode->firstChild;
		$corner2 = $corner1->nextSibling;
		
		list($y1,$x1) = explode(' ',$corner1->nodeValue);
		list($y2,$x2) = explode(' ',$corner2->nodeValue);

		$this->addPoint($x1, $y1);
		$this->addPoint($x1, $y2);
		$this->addPoint($x2, $y2);
		$this->addPoint($x2, $y1);
		$this->addPoint($x1, $y1);
	}
}


class GMLMultiLine {

	var $lineArray = array();

	public function __construct() {
		
	}

	public function parseMultiLine ($domNode) {
		$simpleXMLNode = simplexml_import_dom($domNode);

		$simpleXMLNode->registerXPathNamespace('gml', 'http://www.opengis.net/gml');
		
		$allCoords = $simpleXMLNode->xpath("gml:lineStringMember/gml:LineString/gml:coordinates");
			
		$cnt=0;
		foreach ($allCoords as $Coords) {
			
			$this->lineArray[$cnt] = array();
			
			$coordsDom = dom_import_simplexml($Coords);
				
//			$name = $coordsDom->nodeName;
//			$value = $coordsDom->nodeValue;				
//			echo "===> name: ".$name. ", Value: ".$value."<br>";
			
			foreach(explode(' ',$coordsDom->nodeValue) as $pointCoords){
				list($x,$y,$z) = explode(',',$pointCoords);
				$this->addPoint($x, $y, $cnt);
				}
			
			$cnt++;
		}
		
	}
	
	protected function addPoint ($x, $y, $i) {
		array_push($this->lineArray[$i], array("x" => $x, "y" => $y));
	}
	
	public function toGeoJSON () {
		$numberlineArray = count($this->lineArray);
		$str = "";
		if ($numberlineArray > 0) {
			$str .= "{\"type\": \"MultiLineString\", \"coordinates\":[";
			
			for ($cnt =0; $cnt < $numberlineArray; $cnt++){
				if ($cnt > 0) {
						$str .= ",";
					}
					$str .="[";
			
				for ($i=0; $i < count($this->lineArray[$cnt]); $i++) {
					if ($i > 0) {
						$str .= ",";
					}
					$str .= "[".$this->lineArray[$cnt][$i]["x"].",".$this->lineArray[$cnt][$i]["y"]."]";
				}
				$str .="]";
			}
			$str .= "]}";
			
		}
		else {
			$e = new mb_exception("GMLMultiLine: toGeoJSON: this multiLine is null.");
		}
		return $str;
	}
}

class GMLMultiPolygon {

	var $polygonArray = array();
	var $innerRingArray = array();

	public function __construct() {
		
	}

	protected function addPointToRing ($i, $j, $x, $y) {
		if (count($this->innerRingArray[$i]) <= $j) {
			array_push($this->innerRingArray[$i], array());
		}
		array_push($this->innerRingArray[$i][$j], array("x" => $x, "y" => $y));
	}
	

	public function parseMultiPolygon ($domNode) {
//		echo $domNode->nodeName."<br>";
		$simpleXMLNode = simplexml_import_dom($domNode);

		$simpleXMLNode->registerXPathNamespace('gml', 'http://www.opengis.net/gml');

		$allPolygons = $simpleXMLNode->xpath("gml:polygonMember/gml:Polygon");
		
		$cnt=0;
		foreach ($allPolygons as $polygon) {
			$allCoords = $polygon->xpath("gml:outerBoundaryIs/gml:LinearRing/gml:coordinates");
				
			$this->polygonArray[$cnt] = array();
			foreach ($allCoords as $Coords) {
				
				$coordsDom = dom_import_simplexml($Coords);
					
				foreach (explode(' ',$coordsDom->nodeValue) as $pointCoords) {
					list($x,$y,$z) = explode(',',$pointCoords);
					$this->addPoint($x, $y, $cnt);
				}
			}
			
			$this->innerRingArray[$cnt] = array();
			$innerRingNodeArray = $polygon->xpath("gml:innerBoundaryIs");
			if ($innerRingNodeArray) {
				$ringCount = 0;
				foreach ($innerRingNodeArray as $ringNode) {
					$currentRingNode = $ringNode->xpath("gml:LinearRing");
					foreach ($currentRingNode as $node) {
						$coordinates = $node->xpath("gml:coordinates");
						foreach ($coordinates as $coordinate) {
							$coordsDom = dom_import_simplexml($coordinate);
								
							foreach(explode(' ',$coordsDom->nodeValue) as $pointCoords){
				
								list($x,$y,$z) = explode(',',$pointCoords);
								$this->addPointToRing($cnt, $ringCount, $x, $y);
							}
						}
						$ringCount++;
						
					}
				}
			}
			$cnt++;
			new mb_exception("create multipolygon " . serialize($this->innerRingArray));
		}		
	}
	
	protected function addPoint ($x, $y, $i) {

		array_push($this->polygonArray[$i], array("x" => $x, "y" => $y));
	}
	
	public function toGeoJSON () {
		$numberPolygonArray = count($this->polygonArray);
		$str = "";
		if ($numberPolygonArray > 0) {
			$str .= "{\"type\": \"MultiPolygon\", \"coordinates\":[";
			
			for ($cnt =0; $cnt < $numberPolygonArray; $cnt++){
				if ($cnt > 0) {
					$str .= ",";
				}
				$str .= "[";

				$str .= "[";
				for ($i=0; $i < count($this->polygonArray[$cnt]); $i++) {
					if ($i > 0) {
						$str .= ",";
					}
					$str .= "[".$this->polygonArray[$cnt][$i]["x"].",".$this->polygonArray[$cnt][$i]["y"]."]";
				}
				$str .= "]";
				
				for ($i=0; $i < count($this->innerRingArray[$cnt]); $i++) {
					$str .= ",[";
					for ($j=0; $j < count($this->innerRingArray[$cnt][$i]); $j++) {
						if ($j > 0) {
							$str .= ",";
						}
						$str .= "[".$this->innerRingArray[$cnt][$i][$j]["x"].",".$this->innerRingArray[$cnt][$i][$j]["y"]."]";
					}
					$str .= "]";
				}
				$str .= "]";
			}
			$str .= "]}";
			
		}
		else {
			$e = new mb_exception("GMLMultiPolygon: toGeoJSON: this multiLine is null.");
		}
		return $str;
	}
}
?>