Hoe perceelgrenzen op te vragen via lattitude en longitude

[JvG]

Wie kan mij vertellen hoe ik de kadastrale grenzen van een perceel via een url kan opvragen op basis van een gegeven lattitude en longitude?

Inmiddels opgelost:

1. lat/lng vertalen naar EPSG 28992:
   http://epsg.io/trans?data=lng,lat&s_srs=4326&&t_srs=28992

2. Gevonden x en y - coördinaat gebruiken in CQL Filter:
   https://geodata.nationaalgeoregister.nl/kadastralekaartv3/wfs?service=WFS&request=getfeature&typename=kadastralekaartv3:perceel&SRS=EPSG:4326&CQL_FILTER=INTERSECTS(kadastralekaartv3:begrenzingperceel,Point(x y))

[wouter.visscher]

@JvG Yup, that will work

Misschien wel goed om te weten is dat CQL_FILTER niet onderdeel is van de WMS spec, maar een ‘vendor’ property van Geoserver. Officieel is het de query parameter FILTER degene die daarvoor geschikt is. Maar die is iets lastiger in elkaar te schroeven.

Je zou dan het volgende request kunnen doen:

https://geodata.nationaalgeoregister.nl/kadastralekaartv3/wfs?request=getfeature&service=wfs&version=1.0.0&outputFormat=GML2&typename=kadastralekaartv3:perceel&FILTER=<Filter><Contains><PropertyName>begrenzingperceel<%2FPropertyName><Point%20srsName%3D"http%3A%2F%2Fwww.opengis.net%2Fgml%2Fsrs%2Fepsg.xml%234326"><coordinates%20decimal%3D"."%20cs%3D"%2C">5.510354%2C52.871504<%2Fcoordinates><%2FPoint><%2FContains><%2FFilter>

De waarde achter de query parameter FILTER is een decoded xml waard:

<Filter>
	<Contains>
		<PropertyName>begrenzingperceel</PropertyName>
		<Point srsName="http://www.opengis.net/gml/srs/epsg.xml#4326">
			<coordinates decimal="." cs=",">5.510354,52.871504</coordinates>
		</Point>
	</Contains>
</Filter>

Je kan dan direct met de lat/lon waarde communiceren met de wfs, je krijgt wel in je response de geometry in RD terug. Voordeel is dat de parameter FILTER tegen verschillende WFS implementaties Geoserver, Mapserver, ESRI werkt

Vaak worden dit soort requesten die jij wil maken via een POST request gedaan:

https://geodata.nationaalgeoregister.nl/kadastralekaartv3/wfs?

<wfs:GetFeature service="WFS" version="1.0.0"
  outputFormat="GML2"
  xmlns:wfs="http://www.opengis.net/wfs"
  xmlns:ogc="http://www.opengis.net/ogc"
  xmlns:gml="http://www.opengis.net/gml"
  xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance"
  xsi:schemaLocation="http://www.opengis.net/wfs http://schemas.opengis.net/wfs/1.0.0/WFS-basic.xsd">
  <wfs:Query typeName="kadastralekaartv3:perceel">
		<ogc:Filter>
			<ogc:Contains>
				<ogc:PropertyName>begrenzingperceel</ogc:PropertyName>
				<gml:Point srsName="http://www.opengis.net/gml/srs/epsg.xml#28992">
					<gml:coordinates decimal="." cs=",">163298.25,542710.92</gml:coordinates>
				</gml:Point>
			</ogc:Contains>
		</ogc:Filter>
    </wfs:Query>
</wfs:GetFeature>

[JvG]

Super! Da’s inderdaad een stap minder. Dank voor je reactie Wouter.