Screen land for projects near a grid connection

-- Land-use blocks >= 20 ha, within 2 km of a substation, outside protected areas, inside Kakheti
WITH kakheti AS (
  SELECT geom FROM osm.place_polygon WHERE osm_id = -2000091
),
substations AS (
  SELECT geom FROM osm.infrastructure_point
  WHERE  osm_type = 'power' AND osm_subtype = 'substation' AND geom && ST_MakeEnvelope(44.9, 41.0, 46.8, 42.6, 4326)
  UNION ALL
  SELECT ST_Centroid(geom) FROM osm.infrastructure_polygon
  WHERE  osm_type = 'power' AND osm_subtype = 'substation' AND geom && ST_MakeEnvelope(44.9, 41.0, 46.8, 42.6, 4326)
),
protected AS (
  SELECT geom FROM osm.place_polygon
  WHERE  boundary IN ('protected_area', 'national_park') AND geom && ST_MakeEnvelope(44.9, 41.0, 46.8, 42.6, 4326)
)
SELECT l.osm_id, l.osm_type AS land_use,
       ROUND((ST_Area(l.geom::geography) / 10000)::numeric, 1) AS hectares, l.geom
FROM   osm.landuse_polygon AS l, kakheti AS k
WHERE  l.osm_type IN ('farmland', 'meadow', 'vineyard', 'orchard')
  AND  l.geom && k.geom AND ST_Intersects(l.geom, k.geom)
  AND  ST_Area(l.geom::geography) >= 200000
  AND  EXISTS (
         SELECT 1 FROM substations AS s
         WHERE  s.geom && ST_Expand(l.geom, 0.03)
           AND  ST_DWithin(l.geom::geography, s.geom::geography, 2000)
       )
  AND  NOT EXISTS (SELECT 1 FROM protected AS p WHERE ST_Intersects(l.geom, p.geom));