Add CoY fields to PLUTO

Author: alexrichey

products/pluto/pluto_build/sql/create.sql

@@ -101,5 +101,9 @@ CREATE TABLE pluto (
     basempdate text,
     masdate text,
     polidate text,
-    edesigdate text
+    edesigdate text,
+    mnffar text,
+    affresfar text,
+    trnstzone text,
+    miharea text
 );

products/pluto/pluto_build/sql/export.sql

@@ -91,7 +91,11 @@ SELECT
     dcpedited::varchar(3),
     latitude::numeric(19, 7),
     longitude::numeric(19, 7),
-    notes::varchar(20)
+    notes::varchar(20),
+    mnffar::numeric(19, 11),
+    affresfar::numeric(19, 11),
+    trnstzone::varchar(1),
+    miharea::varchar(50)
 INTO export_pluto
 FROM pluto;
 
@@ -197,7 +201,9 @@ SELECT
     a.dcpedited AS "DCPEdited",
     a.latitude AS "Latitude",
     a.longitude AS "Longitude",
-    a.notes AS "Notes"
+    a.notes AS "Notes",
+    a.trnstzone AS "TrnstZone",
+    a.miharea AS "MIHArea"
 INTO unmapped
 FROM export_pluto AS a
 WHERE bbl::bigint IN (

products/pluto/pluto_build/sql/export_mappluto_shp.sql

@@ -92,6 +92,10 @@ SELECT
     a.latitude AS "Latitude",
     a.longitude AS "Longitude",
     a.notes AS "Notes",
+    a.mnffar AS "ManuFAR",
+    a.affresfar AS "AffResFAR",
+    a.trnstzone AS "TrnstZone",
+    a.miharea AS "MIHArea",
     round(st_length(b.:GEOM)::numeric, 11)::numeric(19, 7) AS "Shape_Leng",
     round(st_area(b.:GEOM)::numeric, 11)::numeric(19, 7) AS "Shape_Area",
     st_makevalid(b.:GEOM) AS geom

products/pluto/pluto_build/sql/far.sql

@@ -10,15 +10,19 @@ UPDATE pluto a
 SET
     residfar = b.residfar,
     commfar = b.commfar,
-    facilfar = b.facilfar
+    facilfar = b.facilfar,
+    mnffar = b.mnffar,
+    affresfar = b.affresfar
 FROM dcp_zoning_maxfar AS b
 WHERE a.zonedist1 = b.zonedist;
 -- zoning district 1 with / first p art 
 UPDATE pluto a
 SET
     residfar = (coalesce(a.residfar, b.residfar)),
     commfar = (coalesce(a.commfar, b.commfar)),
-    facilfar = (coalesce(a.facilfar, b.facilfar))
+    facilfar = (coalesce(a.facilfar, b.facilfar)),
+    mnffar = (coalesce(a.mnffar, b.mnffar)),
+    affresfar = (coalesce(a.affresfar, b.affresfar))
 FROM dcp_zoning_maxfar AS b
 WHERE split_part(a.zonedist1, '/', 1) = b.zonedist;
 
@@ -27,7 +31,9 @@ UPDATE pluto a
 SET
     residfar = (coalesce(a.residfar, b.residfar)),
     commfar = (coalesce(a.commfar, b.commfar)),
-    facilfar = (coalesce(a.facilfar, b.facilfar))
+    facilfar = (coalesce(a.facilfar, b.facilfar)),
+    mnffar = (coalesce(a.mnffar, b.mnffar)),
+    affresfar = (coalesce(a.affresfar, b.affresfar))
 FROM dcp_zoning_maxfar AS b
 WHERE split_part(a.zonedist1, '/', 2) = b.zonedist;
 
@@ -36,7 +42,9 @@ UPDATE pluto a
 SET
     residfar = (coalesce(a.residfar, b.residfar)),
     commfar = (coalesce(a.commfar, b.commfar)),
-    facilfar = (coalesce(a.facilfar, b.facilfar))
+    facilfar = (coalesce(a.facilfar, b.facilfar)),
+    mnffar = (coalesce(a.mnffar, b.mnffar)),
+    affresfar = (coalesce(a.affresfar, b.affresfar))
 FROM dcp_zoning_maxfar AS b
 WHERE a.zonedist2 = b.zonedist;
 
@@ -45,7 +53,9 @@ UPDATE pluto a
 SET
     residfar = (coalesce(a.residfar, b.residfar)),
     commfar = (coalesce(a.commfar, b.commfar)),
-    facilfar = (coalesce(a.facilfar, b.facilfar))
+    facilfar = (coalesce(a.facilfar, b.facilfar)),
+    mnffar = (coalesce(a.mnffar, b.mnffar)),
+    affresfar = (coalesce(a.affresfar, b.affresfar))
 FROM dcp_zoning_maxfar AS b
 WHERE split_part(a.zonedist2, '/', 1) = b.zonedist;
 
@@ -54,7 +64,9 @@ UPDATE pluto a
 SET
     residfar = (coalesce(a.residfar, b.residfar)),
     commfar = (coalesce(a.commfar, b.commfar)),
-    facilfar = (coalesce(a.facilfar, b.facilfar))
+    facilfar = (coalesce(a.facilfar, b.facilfar)),
+    mnffar = (coalesce(a.mnffar, b.mnffar)),
+    affresfar = (coalesce(a.affresfar, b.affresfar))
 FROM dcp_zoning_maxfar AS b
 WHERE split_part(a.zonedist2, '/', 2) = b.zonedist;
 
@@ -63,7 +75,9 @@ UPDATE pluto a
 SET
     residfar = (coalesce(a.residfar, b.residfar)),
     commfar = (coalesce(a.commfar, b.commfar)),
-    facilfar = (coalesce(a.facilfar, b.facilfar))
+    facilfar = (coalesce(a.facilfar, b.facilfar)),
+    mnffar = (coalesce(a.mnffar, b.mnffar)),
+    affresfar = (coalesce(a.affresfar, b.affresfar))
 FROM dcp_zoning_maxfar AS b
 WHERE a.zonedist3 = b.zonedist;
 
@@ -72,7 +86,9 @@ UPDATE pluto a
 SET
     residfar = (coalesce(a.residfar, b.residfar)),
     commfar = (coalesce(a.commfar, b.commfar)),
-    facilfar = (coalesce(a.facilfar, b.facilfar))
+    facilfar = (coalesce(a.facilfar, b.facilfar)),
+    mnffar = (coalesce(a.mnffar, b.mnffar)),
+    affresfar = (coalesce(a.affresfar, b.affresfar))
 FROM dcp_zoning_maxfar AS b
 WHERE split_part(a.zonedist3, '/', 1) = b.zonedist;
 
@@ -81,7 +97,9 @@ UPDATE pluto a
 SET
     residfar = (coalesce(a.residfar, b.residfar)),
     commfar = (coalesce(a.commfar, b.commfar)),
-    facilfar = (coalesce(a.facilfar, b.facilfar))
+    facilfar = (coalesce(a.facilfar, b.facilfar)),
+    mnffar = (coalesce(a.mnffar, b.mnffar)),
+    affresfar = (coalesce(a.affresfar, b.affresfar))
 FROM dcp_zoning_maxfar AS b
 WHERE split_part(a.zonedist3, '/', 2) = b.zonedist;
 
@@ -92,4 +110,9 @@ SET
     = (CASE WHEN a.residfar IS NULL OR a.residfar = '-' THEN 0::double precision ELSE a.residfar::double precision END),
     commfar
     = (CASE WHEN a.commfar IS NULL OR a.commfar = '-' THEN 0::double precision ELSE a.commfar::double precision END),
-    facilfar = (CASE WHEN a.facilfar IS NULL OR a.facilfar = '-' THEN 0 ELSE a.facilfar::double precision END);
+    facilfar
+    = (CASE WHEN a.facilfar IS NULL OR a.facilfar = '-' THEN 0 ELSE a.facilfar::double precision END),
+    mnffar
+    = (case when a.mnffar is null or a.mnffar = '-' then 0 else a.mnffar::double precision end),
+    affresfar
+    = (case when a.affresfar is null or a.affresfar = '-' then 0 else a.affresfar::double precision end);

products/pluto/pluto_build/sql/miharea.sql

@@ -0,0 +1,65 @@
+-- calculate how much (total area and percentage) of each lot is covered by MIH areas
+-- assign the MIH project ID to each tax lot based on which MIH area covers the majority of the lot
+-- a MIH area is only assigned if more than 10% of the lot is covered by the MIH area
+-- OR more than 50% of the MIH area overlaps with the lot
+
+DROP TABLE IF EXISTS mihperorder;
+CREATE TABLE mihperorder AS
+WITH mihper AS (
+    SELECT
+        p.id,
+        p.bbl,
+        n.project_id,
+        ST_AREA(
+            CASE
+                WHEN ST_COVEREDBY(p.geom, n.wkb_geometry) THEN p.geom
+                ELSE ST_MULTI(ST_INTERSECTION(p.geom, n.wkb_geometry))
+            END
+        ) AS segbblgeom,
+        ST_AREA(p.geom) AS allbblgeom,
+        ST_AREA(
+            CASE
+                WHEN ST_COVEREDBY(n.wkb_geometry, p.geom) THEN n.wkb_geometry
+                ELSE ST_MULTI(ST_INTERSECTION(n.wkb_geometry, p.geom))
+            END
+        ) AS segmihgeom,
+        ST_AREA(n.wkb_geometry) AS allmihgeom
+    FROM pluto AS p
+    INNER JOIN dcp_mih AS n
+        ON ST_INTERSECTS(p.geom, n.wkb_geometry)
+),
+
+grouped AS (
+    SELECT
+        id,
+        bbl,
+        project_id,
+        SUM(segbblgeom) AS segbblgeom,
+        SUM(segmihgeom) AS segmihgeom,
+        SUM(segbblgeom / allbblgeom) * 100 AS perbblgeom,
+        MAX(segmihgeom / allmihgeom) * 100 AS maxpermihgeom
+    FROM mihper
+    GROUP BY id, bbl, project_id
+)
+
+SELECT
+    id,
+    bbl,
+    project_id,
+    segbblgeom,
+    perbblgeom,
+    maxpermihgeom,
+    ROW_NUMBER() OVER (
+        PARTITION BY id
+        ORDER BY segbblgeom DESC, segmihgeom DESC
+    ) AS row_number
+FROM grouped
+WHERE perbblgeom >= 10 OR maxpermihgeom >= 50;
+
+-- assign the MIH project ID with the highest overlap to each lot
+UPDATE pluto a
+SET miharea = project_id
+FROM mihperorder AS b
+WHERE
+    a.id = b.id
+    AND row_number = 1;

products/pluto/pluto_build/sql/transitzone.sql

@@ -0,0 +1,8 @@
+-- assign transit zone designation based on spatial intersection
+-- with the greater transit zone area
+-- output is 'Y' if lot intersects with transit zone, NULL otherwise
+
+UPDATE pluto a
+SET trnstzone = 'Y'
+FROM dcp_greater_transit_zone AS b
+WHERE ST_INTERSECTS(a.geom, b.wkb_geometry);