OpenStreetMap道路绿色分数计算异常，求排查原因

核心问题：错误的聚合逻辑导致重复计算 #

你当前的SQL会对每个匹配的绿色多边形，单独计算它与道路缓冲区的交集占缓冲区面积的比例，再把所有比例累加。如果一条道路的缓冲区和多个绿色多边形重叠——不管这些多边形是否有区域重叠——每个多边形的占比都会被相加，最终总和必然超过1，甚至出现几十、几百的极端值：

• 若多个多边形无重叠，比如道路同时挨着森林和水域，各占缓冲区的0.4和0.5，总和会是0.9（这看起来合理，但本质逻辑错误）；
• 若多个多边形有重叠区域，比如道路被多个重叠的森林斑块覆盖，每个斑块的交集占比都会被重复计算，总和会急剧膨胀。

你的需求是道路周边30米范围内绿色区域的实际占比，应该计算所有绿色区域在缓冲区中的总面积，再除以缓冲区面积，而不是累加单个多边形的占比。

次要问题：连接条件不符合需求 #

当前JOIN条件ST_Intersects(l.way, ST_Buffer(p.way,30))判断的是「道路线本身是否在多边形的30米缓冲区内」，但你要的是「道路周边30米范围内的绿色区域」，应该改为判断道路的30米缓冲区是否与多边形相交，即ST_Intersects(ST_Buffer(l.way,30), p.way)，这样才能精准匹配符合范围的多边形。

修正后的SQL #

SELECT                                                                 
    l.osm_id,
    -- 合并所有交集区域计算总面积，再除以缓冲区面积；无匹配时返回0
    COALESCE(
        ST_Area(ST_Union(ST_Intersection(ST_Buffer(l.way, 30), p.way)))
        /
        ST_Area(ST_Buffer(l.way, 30)),
        0
    ) AS green_fraction
FROM planet_osm_line AS l
LEFT JOIN planet_osm_polygon AS p 
    ON ST_Intersects(ST_Buffer(l.way, 30), p.way)
    AND (p.natural = 'water' OR p.landuse = 'forest')
GROUP BY l.osm_id;

关键修正点： #

• 使用ST_Union合并所有多边形与道路缓冲区的交集，彻底避免重叠区域的重复计算，得到绿色区域的实际总面积；
• 改用LEFT JOIN确保没有匹配绿色多边形的道路也能被返回，此时分数为0；
• 将筛选条件（natural/landuse）移至JOIN子句，避免过滤掉无匹配的道路；
• 用COALESCE处理无匹配时的NULL值，确保分数始终在0-1范围内。

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