se/SE/se-lib/Api/WfsQgisServer.php

<?php

Lib::loadClass('Api_WfsServerBase');
Lib::loadClass('Api_WfsException');
Lib::loadClass('Api_WfsGeomTypeConverter');

class Api_WfsQgisServer extends Api_WfsServerBase {

	public function parseXMLRequest() {
		$data = array();
		$reqContent = Request::getRequestBody();
		if (empty($reqContent)) {
			throw new Exception("Empty request");
		}
		$parserXml = xml_parser_create();
		xml_parser_set_option($parserXml, XML_OPTION_CASE_FOLDING, 0);
		xml_parser_set_option($parserXml, XML_OPTION_SKIP_WHITE, 1);
		if (0 == xml_parse_into_struct($parserXml, $reqContent, $tags)) {
			throw new Exception("Error parsing xml");
		}
		xml_parser_free($parserXml);
		if (empty($tags)) {
			throw new Exception("Empty structure from request");
		}

		$rootTagName = V::get('tag', '', $tags[0]);
		if ('Transaction' == $rootTagName) {
			return $this->_parseTransactionXmlStruct($reqContent, $tags);
		}

		throw new Api_WfsException("TODO ... L." . __LINE__, 501);

		$xml = new SimpleXMLElement($reqContent);
		$namespaces = $xml->getNameSpaces(true);

		if ('Transaction' == $xml->getName()) {
			$this->_parseTransactionXml($xml);
		}
		else {
			throw new Api_WfsException("Not Implemented " . htmlspecialchars($xml->getName()), 501);
		}
	}

	public function getFeatureAction() {
		$type = V::get('TYPENAME', '', $_REQUEST);
		$typeEx = explode(':', $type);
		$maxFeatures = V::get('MAXFEATURES', '10000', $_REQUEST, 'int');// TODO: Set Deafult Limit
		$srsname = V::get('SRSNAME', '', $_REQUEST);// eg. EPSG:4326
		if (count($typeEx) != 2) throw new HttpException("Wrong param TYPENAME", 400);
		return $this->getFeatures($typeEx[0], $typeEx[1], $maxFeatures, $srsname);
	}

	public function getFeatures($nsPrefix, $type, $maxFeatures, $srsname) {
		$DBG = (V::get('DBG_GEO', '', $_GET) > 0);// TODO: Profiler
		$typeName = "{$nsPrefix}:{$type}";
		$acl = $this->getAclFromTypeName($typeName);
		$fldList = $this->_getFieldListFromAcl($acl);

		$baseNsUri = $this->getBaseNamespaceUri();
		//$wfsNs = 'p5_default_db_' . $type;//$nsPrefix;
		$wfsNs = 'p5_default_db';//$nsPrefix;
		$wfsNsUri = "{$baseNsUri}/" . substr($nsPrefix, 3);
		$featureTypeUri = $this->getBaseUri() . "?SERVICE=WFS&VERSION=1.0.0&TYPENAME={$typeName}&REQUEST=DescribeFeatureType";
		// https://biuro.biall-net.pl/dev-pl/se-master/wfs-qgis.php/default_db/
		// https://biuro.biall-net.pl/dev-pl/se-master/wfs-data.php/default_db/TEST_PERMS/?SERVICE=WFS&VERSION=1.0.0&TYPENAME=p5_default_db:TEST_PERMS&REQUEST=DescribeFeatureType

		// get BBox from geom_field (only one geom fld is allowed)
		$geomFld = null;
		{
			foreach ($fldList as $idZasob => $fldName) {
				if ($acl->isGeomField($fldName)) {
					$geomFld = $fldName;
				}
			}
		}

		$dbGeomType = $acl->getGeomFieldType($geomFld);

		$searchParams = array();
		$searchParams['limit'] = $maxFeatures;
		$searchParams['order_by'] = $acl->getPrimaryKeyField();
		$searchParams['order_dir'] = 'DESC';
		{// BBOX
			// 54.26931096743426,18.48242909824306,54.26738118403914,18.478738378639246
			$bbox = V::get('BBOX', '', $_GET);
			if (!empty($bbox)) {
				if (preg_match("/^\d+(.\d+)?,\d+(.\d+)?,\d+(.\d+)?,\d+(.\d+)?$/", $bbox, $matches)) {
					$searchParams['f_the_geom'] = "BBOX:{$bbox}";
				} else {
					// throw new Exception("Error Processing Request", 1);// ?
					$this->DBG("Error BBOX wrong format", __LINE__, __FUNCTION__, __CLASS__);
				}
			}
		}
		//if ($geomFld) $searchParams["f_{$geomFld}"] = 'IS NOT NULL';
		//if ($geomFld) $searchParams["f_{$geomFld}"] = 'GeometryType=' . strtoupper($dbGeomType);
		$geomType = strtoupper($dbGeomType);
		if ($geomFld) $searchParams["ogc:Filter"] = <<<OGC_FILTER
<ogc:Filter>
	<ogc:And>
		<ogc:Not>
			<ogc:PropertyIsNull>
				<ogc:PropertyName>{$geomFld}</ogc:PropertyName>
			</ogc:PropertyIsNull>
		</ogc:Not>
		<ogc:PropertyIsEqualTo>
			<ogc:Function name="GeometryType">
				 <ogc:PropertyName>{$geomFld}</ogc:PropertyName>
			</ogc:Function>
			<ogc:Literal>{$geomType}</ogc:Literal>
		</ogc:PropertyIsEqualTo>
	</ogc:And>
</ogc:Filter>
OGC_FILTER;
if($DBG){echo 'getItems:';print_r($searchParams);echo "\n";}
		$this->DBG("getItems:" . json_encode($searchParams), __LINE__, __FUNCTION__, __CLASS__);
		$items = $acl->getItems($searchParams);
		$this->DBG("items(" . count($items) . ")", __LINE__, __FUNCTION__, __CLASS__);

		header('Content-type: application/xml; charset=utf-8');
		$xmlWriter = new XMLWriter();
		$xmlWriter->openUri('php://output');
		// $xmlWriter->openMemory();// DBG
		$xmlWriter->setIndent(true);
		if (!$xmlWriter) throw new HttpException("Error no XMLWriter", 404);
		$xmlWriter->startDocument('1.0','UTF-8');
		//$xmlWriter->startElementNS('wfs', 'FeatureCollection', 'http://www.opengis.net/wfs');
		$xmlWriter->startElement('wfs:FeatureCollection');
//		$xmlWriter->writeAttributeNS('xmlns', 'wfs', 'http://www.w3.org/2000/xmlns/', 'http://www.opengis.net/wfs');
		$xmlWriter->writeAttribute('xmlns:wfs', 'http://www.opengis.net/wfs');
		$xmlWriter->writeAttribute('xmlns', 'http://www.opengis.net/wfs');
//		$xmlWriter->writeAttributeNS('xmlns', 'gml', 'http://www.w3.org/2000/xmlns/', 'http://www.opengis.net/gml');
//		$xmlWriter->writeAttributeNS('xmlns', 'xsi', 'http://www.w3.org/2000/xmlns/', 'http://www.w3.org/2001/XMLSchema-instance');
//		$xmlWriter->writeAttributeNS('xmlns', $wfsNs, 'http://www.w3.org/2000/xmlns/', $wfsNsUri);
		$xmlWriter->writeAttribute('xmlns:gml', 'http://www.opengis.net/gml');
		$xmlWriter->writeAttribute('xmlns:xsi', 'http://www.w3.org/2001/XMLSchema-instance');
		$xmlWriter->writeAttribute("xmlns:{$wfsNs}", $wfsNsUri);
		$xmlWriter->writeAttribute('xsi:schemaLocation', "{$wfsNsUri} {$featureTypeUri}");

if($DBG){echo '(geomFld: '.$geomFld.'):';print_r($acl->getFieldType($geomFld));echo "\n";}
		if (empty($items)) {
			$pKeyField = $acl->getPrimaryKeyField();
			$fakeItem = new stdClass();
			$fakeItem->{$pKeyField} = 0;
			if ('polygon' == $dbGeomType) {
				$fakeItem->the_geom = "POLYGON(())";
			} else if ('linestring' == $dbGeomType) {
				$fakeItem->the_geom = "LINESTRING()";
			} else if ('point' == $dbGeomType) {
				$fakeItem->the_geom = "POINT(0,0)";
			}
			$items[0] = $fakeItem;
		}
		$dbgLoop = 0;
		$this->DBG("before loop...", __LINE__, __FUNCTION__, __CLASS__);
		foreach ($items as $itemKey => $item) {
			if (0 == (++$dbgLoop) % 500) $this->DBG("items loop:{$dbgLoop}", __LINE__, __FUNCTION__, __CLASS__);
			if($DBG){echo 'item['.$itemKey.'] ('.$geomFld.')isEmpty('.empty($item->{$geomFld}).'):';print_r($item->{$geomFld});echo "\n";}
			if ($geomFld) {
				if (empty($item->{$geomFld})) {
					continue;// QGIS crash when WFS contain features with empty geom field
				}
			}
			$xmlWriter->startElement('gml:featureMember');
				$xmlWriter->startElement("{$wfsNs}:{$type}");
					$xmlWriter->writeAttribute('fid', "{$type}.{$itemKey}");
					foreach ($fldList as $idZasob => $fldName) {
						// if (!$acl->isAllowed($idZasob, 'R', $item)) {echo '<!-- '."[{$idZasob}]({$fldName}) not allowed to read: ";print_r($item);echo' -->';}
						if (!$acl->isAllowed($idZasob, 'R', $item)) continue;
						// if ($acl->isGeomField($fldName)) // BUG: wolno
						if ($geomFld != null && $fldName == $geomFld) {
							$xmlWriter->startElement("{$wfsNs}:{$fldName}");
								$this->_typeConverter->createGmlFromWkt_xmlWriter($item->{$fldName}, $xmlWriter);
							$xmlWriter->endElement();// {$wfsNs}:{$fldName}
						} else {
							$value = str_replace('&', '&amp;', $item->{$fldName});
							if (empty($value) && '0' !== $value) {
								continue;
							} else {
								$xmlWriter->startElement("{$wfsNs}:{$fldName}");
								$xmlWriter->text($value);
								$xmlWriter->endElement();// {$wfsNs}:{$fldName}
							}
						}
					}
				$xmlWriter->endElement();// {$wfsNs}:{$type}
			$xmlWriter->endElement();// gml:featureMember
		}
		$xmlWriter->endElement();// wfs:FeatureCollection
		$xmlWriter->endDocument();
		$this->DBG("items loop END", __LINE__, __FUNCTION__, __CLASS__);
		exit;
	}

	public function describeFeatureTypeAction() {
		$type = V::get('TYPENAME', '', $_REQUEST);
		if (empty($type)) {
			$reqContent = Request::getRequestBody();
			if (!empty($reqContent)) {
				return $this->_parseDescribeFeatureTypeRequest($reqContent);
			} else {
				return $this->_getDescribeFeatureAllTypes();
			}
			//throw new HttpException("Wrong param TYPENAME", 400);
		}
		$typeEx = explode(':', $type);
		if (count($typeEx) != 2) {
			throw new HttpException("Wrong param TYPENAME", 400);
		}
		return $this->_getDescribeFeatureType($typeEx[0], $typeEx[1]);
	}

	public function getCapabilitiesAction() {
		$wfsServerUrl = $this->getBaseUri();
		$serviceTitle = "Web Feature Service for QGIS";
		$serviceDescription = "This is the reference implementation of WFS 1.0.0 and WFS 1.1.0, supports all WFS operations including Transaction.";

		//header('Content-type: application/xml; charset="utf-8"');
		header('Content-type: application/xml');
		$this->_getCapabilities($wfsServerUrl, $serviceTitle, $serviceDescription);
		exit;
	}

	public function _getTableAclList() {// Use only Tables from default_db with the_geom field
		$tblAclList = array();
		$db = DB::getDB();
		$idDefaultDB = $db->_zasob_id;
		$fullTblAclList = $this->_usrAcl->getTablesAcl();
		foreach ($fullTblAclList as $tblAcl) {
			$dataSourceName = 'default_db';// TODO: getSourceName
			$tblName = $tblAcl->getName();
			if ($idDefaultDB != $tblAcl->getDB()) {// hide non default_db tables
				continue;
			}
			try {
				$acl = $this->getAclFromTypeName($typeName = "p5_{$dataSourceName}:{$tblName}");
			} catch (Exception $e) {
				// TODO: error log $e->getMessage();
			}
			if (!$acl) {
				// TODO: error log msg
				continue;
			}
			$fldList = $acl->getRealFieldList();
			if (!in_array('the_geom', $fldList)) {
				continue;
			}
			$tblAclList[] = $tblAcl;
		}
		return $tblAclList;
	}

	public function _getFieldListFromAcl($acl) {
		$fldList = $acl->getRealFieldListByIdZasob();
		// mv the_geom to the first place
		$orderedFldList = array();
		foreach ($fldList as $idZasob => $fldName) {
			if ('the_geom' == $fldName) $orderedFldList[$idZasob] = $fldName;
		}
		foreach ($fldList as $idZasob => $fldName) {
			if ('the_geom' == $fldName) continue;
			$orderedFldList[$idZasob] = $fldName;
		}
		return $orderedFldList;
	}

}