利用geoserver地图服务分析数据---占地分析

• 具体可从以下链接查看geoserver的一些地图服务操作
• Web Feature Service (WFS)
• 使用的地图是：MapBox

一、占地分析

思路：

1. 从wms图层中根据图层名称拿到全部的图层属性信息
2. 获取到的图层类型是多边型【MultiPolygon】利用Turf中flatten方法减少层级的嵌套
3. 将得到的GeoJSON对象转换为WGS84投影
4. 使用空间查询的方法进行数据查询，分析出两个图层相交的部分

export const SPACE_QUERY_FUN_L = async (ShpName1, ShpName2) => {
	
	const queryInfo1 = await axios.get(
		`/geoserver/wfs?service=wfs&version=1.1.0&request=GetFeature&outputFormat=application/json&typeName=${ShpName1}`
	)


	
	const e = turf.flatten(queryInfo1.data.features[0].geometry)

	
	const allValues = e.features
	const allFeatures = []
	for (let i = 0; i < allValues.length; i++) {
		const converted = turf.toWgs84(allValues[i])
		allFeatures.push(converted)
	}
	
	const queryInfo = await axios.get(
		`/geoserver/wfs?SERVICE=WFS&REQUEST=GetFeature&VERSION=2.0.0&TYPENAMES=${ShpName1}&OUTPUTFORMAT=application/json&COUNT=10000000000000000&SRSNAME=EPSG:4326&CQL_FILTER=${`INTERSECTS(the_geom, SRID=4326;POLYGON((${String(
			_.chunk(_.flattenDeep(allFeatures[0].geometry.coordinates), 2)
				.map(item => item.join(' '))
				.join(', ')
		)})))`}`
	)

	if (queryInfo.data.features.length) {
		return queryInfo.data
	}
}

参数说明

1. 此次采取的几何方式是Intersects，即采用面与点相交，得到这个面里面的点数据
2. 多边形经纬度为一个数组，第一个点经纬度要与最后一个点一致，保证多边形的闭合
3. 过滤条件里有一个PropertyName属性，为必须字段，根据图层的属性来查看，具体查看方式为点击图层，然后点及要查询的图层名称，要素类型的第一个属性，即为该字段的值