如何用SPARQL查询判断坐标点是否属于指定地理区域？

问题背景 #

我有带sf:Point类型的坐标数据（例如"POINT(51.95656412820512 5.880331774358974)"^^<http://www.opengis.net/ont/geosparql#wktLiteral>），希望通过查询Wikidata的wdt:P3896属性，确定该点所属的国家或行政区域。

现有尝试 #

我已能通过联邦SPARQL从Wikidata获取国家边界数据，查询语句如下：

PREFIX rdf: <http://www.w3.org/1999/02/22-rdf-syntax-ns#>
PREFIX wdt: <http://www.wikidata.org/prop/direct/>
PREFIX wd: <http://www.wikidata.org/entity/>
PREFIX rdfs: <http://www.w3.org/2000/01/rdf-schema#>
SELECT * WHERE {
    SERVICE <https://query.wikidata.org/sparql> {
        ?country wdt:P31 wd:Q6256.
        ?country wdt:P3896 ?boundary .
        FILTER(?country = wd:Q55) .
        ?country rdfs:label  ?label .
        FILTER(lang(?label)="en") .
    }
}

返回结果：

现在需要解决的核心问题：能否直接用SPARQL判断给定坐标点是否在荷兰境内？目前未找到能判断点是否在wdt:P3896对应的地理形状内的SPARQL函数，也无其他可行思路。

可行解决方案 #

Wikidata的SPARQL端点本身不支持直接解析wdt:P3896返回的外部地图文件并执行点-in-polygon判断，需通过以下两种方式实现需求：

1. 借助支持GeoSPARQL的三元组存储+预导入边界数据 #

• 操作步骤：
  • 下载wdt:P3896指向的荷兰边界数据，转换为GeoSPARQL兼容的sf:Polygon/sf:MultiPolygon类型WKT字面量
  • 将转换后的边界几何数据导入自己的三元组存储，关联到对应国家实体（比如新增属性wdt:P3896_geom存储几何数据）
  • 使用GeoSPARQL的sf:within函数编写查询判断点的归属
• 示例查询：

  PREFIX sf: <http://www.opengis.net/ont/sf#>
  PREFIX wdt: <http://www.wikidata.org/prop/direct/>
  PREFIX wd: <http://www.wikidata.org/entity/>
  PREFIX rdfs: <http://www.w3.org/2000/01/rdf-schema#>
  
  SELECT ?country_label WHERE {
      # 定义待判断的坐标点
      BIND("POINT(51.95656412820512 5.880331774358974)"^^<http://www.opengis.net/ont/geosparql#wktLiteral> AS ?point)
      # 关联本地存储的荷兰边界几何数据
      wd:Q55 wdt:P3896_geom ?country_polygon.
      wd:Q55 rdfs:label ?country_label.
      FILTER(lang(?country_label) = "en")
      # 执行点在多边形内的判断
      FILTER(sf:within(?point, ?country_polygon))
  }
  

2. 在SPARQL查询外处理几何判断 #

• 操作步骤：
  • 通过联邦SPARQL获取目标国家的wdt:P3896链接
  • 用脚本语言（Python/JavaScript等）下载并解析该链接对应的边界数据，转换为可处理的几何对象
  • 使用地理处理库执行点-in-polygon判断
• 示例Python代码片段：

  from shapely.wkt import loads
  import requests
  
  # 待判断的坐标点WKT
  point_wkt = "POINT(51.95656412820512 5.880331774358974)"
  point = loads(point_wkt)
  
  # 从Wikidata获取的荷兰边界数据链接
  boundary_url = "http://commons.wikimedia.org/data/main/Data:Nederland.map"
  # 下载并解析边界数据（需根据实际数据格式调整解析逻辑）
  response = requests.get(boundary_url)
  # 假设解析后得到多边形对象
  country_polygon = loads(response.text)
  
  # 判断点是否在境内
  print("该点位于荷兰境内" if point.within(country_polygon) else "该点不在荷兰境内")
  

关键说明 #

• Wikidata的SPARQL端点仅能返回外部数据链接，无法直接解析并处理链接指向的地图文件，因此纯SPARQL无法完成点-in-polygon判断
• GeoSPARQL的空间函数仅支持处理存储在三元组中的几何字面量，无法直接操作外部文件

内容的提问来源于stack exchange，提问作者Sibbe Bakker