Day: October 27, 2025

What private property owner has the most acerage by county? I've been playing with these queries in R, trying to pull this data to make maps and analysis by using overly complex postgresql queries. At some point I'll upload the R code up into Github, but for now, here is the PostgresSQL query that actually ran relatively quickly -- 25 seconds or so on my old laptop.

      SELECT
        primary_owner,
        county_name,
        avg(ST_X(geom)) as x,
        avg(ST_Y(geom)) as y,
        sum(acres) as av
      INTO TEMPORARY propTot
      FROM
        nytax_ctr_point
      WHERE
        "roll_section"='1' AND
        "owner_type"='8' AND
        "acres" < 1000
      GROUP BY  
        primary_owner, county_name;

  SELECT 
    t1.*
  FROM
    propTot t1
  RIGHT JOIN (
    SELECT county_name, max(av) as av
    FROM propTot
    GROUP BY county_name
  ) t2
  ON t1.county_name = t2.county_name AND
  t1.av = t2.av;