Merge pull request #110 from OpenDRR/test_hexbin_unclipped

Test hexbin unclipped
OpenDRR · Mar 31, 2022 · 3423cf4 · 3423cf4
2 parents b5811ed + f279ddd
commit 3423cf4
Show file tree

Hide file tree

Showing 10 changed files with 1,602 additions and 428 deletions.
diff --git a/scripts/Create_hexbin_MH_risk_sauid_prioritization_aggregation_area.sql b/scripts/Create_hexbin_MH_risk_sauid_prioritization_aggregation_area.sql
@@ -1,7 +1,33 @@
 
 -- test aggregation to hexbin grids area proxy
 
---5km
+-- clipped
+-- 1km
+DROP VIEW IF EXISTS results_nhsl_hazard_threat.nhsl_hazard_threat_prioritization_hexgrid_1km;
+CREATE VIEW results_nhsl_hazard_threat.nhsl_hazard_threat_prioritization_hexgrid_1km AS 
+SELECT
+b.gridid_1,
+SUM("Et_BldgNum" * b.area_ratio) AS "Et_BldgNum",
+SUM("Et_AssetValue" * b.area_ratio) AS "Et_AssetValue",
+SUM("Et_PopNight" * b.area_ratio) AS "Et_PopNight",
+SUM(eq_shaking_score_abs * b.area_ratio) AS "eq_shaking_score_abs",
+SUM(eq_shaking_score_rel * b.area_ratio) AS "eq_shaking_score_rel",
+SUM(fld_priority_score_abs * b.area_ratio) AS "fld_priority_score_abs",
+SUM(fld_priority_score_rel * b.area_ratio) AS "fld_priority_score_rel",
+SUM(wildfire_priority_score_abs * b.area_ratio) AS "wildfire_priority_score_abs",
+SUM(wildfire_priority_score_rel * b.area_ratio) AS "wildfire_priority_score_rel",
+SUM(cy_priority_score_abs * b.area_ratio) AS "cy_priority_score_abs",
+SUM(cy_priority_score_rel * b.area_ratio) AS "cy_priority_score_rel",
+SUM(total_abs_score * b.area_ratio) AS "total_abs_score",
+SUM(total_rel_score * b.area_ratio) AS "total_rel_score",
+c.geom
+
+FROM results_nhsl_hazard_threat.nhsl_hazard_threat_indicators_s_tbl a
+LEFT JOIN boundaries."SAUID_HexGrid_1km_intersect" b ON a."Sauid" = b.sauid
+LEFT JOIN boundaries."HexGrid_1km" c ON b.gridid_1 = c.gridid_1
+GROUP BY b.gridid_1,c.geom;
+
+-- 5km
 DROP VIEW IF EXISTS results_nhsl_hazard_threat.nhsl_hazard_threat_prioritization_hexgrid_5km;
 CREATE VIEW results_nhsl_hazard_threat.nhsl_hazard_threat_prioritization_hexgrid_5km AS 
 SELECT
@@ -78,11 +104,11 @@ LEFT JOIN boundaries."HexGrid_25km" c ON b.gridid_25 = c.gridid_25
 GROUP BY b.gridid_25,c.geom;
 
 
--- 50km
-DROP VIEW IF EXISTS results_nhsl_hazard_threat.nhsl_hazard_threat_prioritization_hexgrid_50km;
-CREATE VIEW results_nhsl_hazard_threat.nhsl_hazard_threat_prioritization_hexgrid_50km AS 
+-- global fabric
+DROP VIEW IF EXISTS results_nhsl_hazard_threat.nhsl_hazard_threat_prioritization_hexgrid_global_fabric;
+CREATE VIEW results_nhsl_hazard_threat.nhsl_hazard_threat_prioritization_hexgrid_global_fabric AS 
 SELECT
-b.gridid_50,
+b.gridid,
 SUM("Et_BldgNum" * b.area_ratio) AS "Et_BldgNum",
 SUM("Et_AssetValue" * b.area_ratio) AS "Et_AssetValue",
 SUM("Et_PopNight" * b.area_ratio) AS "Et_PopNight",
@@ -99,15 +125,17 @@ SUM(total_rel_score * b.area_ratio) AS "total_rel_score",
 c.geom
 
 FROM results_nhsl_hazard_threat.nhsl_hazard_threat_indicators_s_tbl a
-LEFT JOIN boundaries."SAUID_HexGrid_50km_intersect" b ON a."Sauid" = b.sauid
-LEFT JOIN boundaries."HexGrid_50km" c ON b.gridid_50 = c.gridid_50
-GROUP BY b.gridid_50,c.geom;
+LEFT JOIN boundaries."SAUID_HexGrid_GlobalFabric_intersect" b ON a."Sauid" = b.sauid
+LEFT JOIN boundaries."HexGrid_GlobalFabric" c ON b.gridid = c.gridid
+GROUP BY b.gridid,c.geom;
 
--- 100km
-DROP VIEW IF EXISTS results_nhsl_hazard_threat.nhsl_hazard_threat_prioritization_hexgrid_100km;
-CREATE VIEW results_nhsl_hazard_threat.nhsl_hazard_threat_prioritization_hexgrid_100km AS 
+
+-- unclipped
+-- 1km
+DROP VIEW IF EXISTS results_nhsl_hazard_threat.nhsl_hazard_threat_prioritization_hexgrid_1km_uc;
+CREATE VIEW results_nhsl_hazard_threat.nhsl_hazard_threat_prioritization_hexgrid_1km_uc AS 
 SELECT
-b.gridid_100,
+b.gridid_1,
 SUM("Et_BldgNum" * b.area_ratio) AS "Et_BldgNum",
 SUM("Et_AssetValue" * b.area_ratio) AS "Et_AssetValue",
 SUM("Et_PopNight" * b.area_ratio) AS "Et_PopNight",
@@ -124,16 +152,41 @@ SUM(total_rel_score * b.area_ratio) AS "total_rel_score",
 c.geom
 
 FROM results_nhsl_hazard_threat.nhsl_hazard_threat_indicators_s_tbl a
-LEFT JOIN boundaries."SAUID_HexGrid_100km_intersect" b ON a."Sauid" = b.sauid
-LEFT JOIN boundaries."HexGrid_100km" c ON b.gridid_100 = c.gridid_100
-GROUP BY b.gridid_100,c.geom;
+LEFT JOIN boundaries."SAUID_HexGrid_1km_intersect_unclipped" b ON a."Sauid" = b.sauid
+LEFT JOIN boundaries."HexGrid_1km_unclipped" c ON b.gridid_1 = c.gridid_1
+GROUP BY b.gridid_1,c.geom;
 
+-- 5km
+DROP VIEW IF EXISTS results_nhsl_hazard_threat.nhsl_hazard_threat_prioritization_hexgrid_5km_uc;
+CREATE VIEW results_nhsl_hazard_threat.nhsl_hazard_threat_prioritization_hexgrid_5km_uc AS 
+SELECT
+b.gridid_5,
+SUM("Et_BldgNum" * b.area_ratio) AS "Et_BldgNum",
+SUM("Et_AssetValue" * b.area_ratio) AS "Et_AssetValue",
+SUM("Et_PopNight" * b.area_ratio) AS "Et_PopNight",
+SUM(eq_shaking_score_abs * b.area_ratio) AS "eq_shaking_score_abs",
+SUM(eq_shaking_score_rel * b.area_ratio) AS "eq_shaking_score_rel",
+SUM(fld_priority_score_abs * b.area_ratio) AS "fld_priority_score_abs",
+SUM(fld_priority_score_rel * b.area_ratio) AS "fld_priority_score_rel",
+SUM(wildfire_priority_score_abs * b.area_ratio) AS "wildfire_priority_score_abs",
+SUM(wildfire_priority_score_rel * b.area_ratio) AS "wildfire_priority_score_rel",
+SUM(cy_priority_score_abs * b.area_ratio) AS "cy_priority_score_abs",
+SUM(cy_priority_score_rel * b.area_ratio) AS "cy_priority_score_rel",
+SUM(total_abs_score * b.area_ratio) AS "total_abs_score",
+SUM(total_rel_score * b.area_ratio) AS "total_rel_score",
+c.geom
 
--- global fabric
-DROP VIEW IF EXISTS results_nhsl_hazard_threat.nhsl_hazard_threat_prioritization_hexgrid_global_fabric;
-CREATE VIEW results_nhsl_hazard_threat.nhsl_hazard_threat_prioritization_hexgrid_global_fabric AS 
+FROM results_nhsl_hazard_threat.nhsl_hazard_threat_indicators_s_tbl a
+LEFT JOIN boundaries."SAUID_HexGrid_5km_intersect_unclipped" b ON a."Sauid" = b.sauid
+LEFT JOIN boundaries."HexGrid_5km_unclipped" c ON b.gridid_5 = c.gridid_5
+GROUP BY b.gridid_5,c.geom;
+
+
+-- 10km
+DROP VIEW IF EXISTS results_nhsl_hazard_threat.nhsl_hazard_threat_prioritization_hexgrid_10km_uc;
+CREATE VIEW results_nhsl_hazard_threat.nhsl_hazard_threat_prioritization_hexgrid_10km_uc AS 
 SELECT
-b.gridid,
+b.gridid_10,
 SUM("Et_BldgNum" * b.area_ratio) AS "Et_BldgNum",
 SUM("Et_AssetValue" * b.area_ratio) AS "Et_AssetValue",
 SUM("Et_PopNight" * b.area_ratio) AS "Et_PopNight",
@@ -150,6 +203,82 @@ SUM(total_rel_score * b.area_ratio) AS "total_rel_score",
 c.geom
 
 FROM results_nhsl_hazard_threat.nhsl_hazard_threat_indicators_s_tbl a
-LEFT JOIN boundaries."SAUID_HexGrid_GlobalFabric_intersect" b ON a."Sauid" = b.sauid
-LEFT JOIN boundaries."HexGrid_GlobalFabric" c ON b.gridid = c.gridid
-GROUP BY b.gridid,c.geom;
+LEFT JOIN boundaries."SAUID_HexGrid_10km_intersect_unclipped" b ON a."Sauid" = b.sauid
+LEFT JOIN boundaries."HexGrid_10km_unclipped" c ON b.gridid_10 = c.gridid_10
+GROUP BY b.gridid_10,c.geom;
+
+-- 25km
+DROP VIEW IF EXISTS results_nhsl_hazard_threat.nhsl_hazard_threat_prioritization_hexgrid_25km_uc;
+CREATE VIEW results_nhsl_hazard_threat.nhsl_hazard_threat_prioritization_hexgrid_25km_uc AS 
+SELECT
+b.gridid_25,
+SUM("Et_BldgNum" * b.area_ratio) AS "Et_BldgNum",
+SUM("Et_AssetValue" * b.area_ratio) AS "Et_AssetValue",
+SUM("Et_PopNight" * b.area_ratio) AS "Et_PopNight",
+SUM(eq_shaking_score_abs * b.area_ratio) AS "eq_shaking_score_abs",
+SUM(eq_shaking_score_rel * b.area_ratio) AS "eq_shaking_score_rel",
+SUM(fld_priority_score_abs * b.area_ratio) AS "fld_priority_score_abs",
+SUM(fld_priority_score_rel * b.area_ratio) AS "fld_priority_score_rel",
+SUM(wildfire_priority_score_abs * b.area_ratio) AS "wildfire_priority_score_abs",
+SUM(wildfire_priority_score_rel * b.area_ratio) AS "wildfire_priority_score_rel",
+SUM(cy_priority_score_abs * b.area_ratio) AS "cy_priority_score_abs",
+SUM(cy_priority_score_rel * b.area_ratio) AS "cy_priority_score_rel",
+SUM(total_abs_score * b.area_ratio) AS "total_abs_score",
+SUM(total_rel_score * b.area_ratio) AS "total_rel_score",
+c.geom
+
+FROM results_nhsl_hazard_threat.nhsl_hazard_threat_indicators_s_tbl a
+LEFT JOIN boundaries."SAUID_HexGrid_25km_intersect_unclipped" b ON a."Sauid" = b.sauid
+LEFT JOIN boundaries."HexGrid_25km_unclipped" c ON b.gridid_25 = c.gridid_25
+GROUP BY b.gridid_25,c.geom;
+
+
+-- 50km
+DROP VIEW IF EXISTS results_nhsl_hazard_threat.nhsl_hazard_threat_prioritization_hexgrid_50km_uc;
+CREATE VIEW results_nhsl_hazard_threat.nhsl_hazard_threat_prioritization_hexgrid_50km_uc AS 
+SELECT
+b.gridid_50,
+SUM("Et_BldgNum" * b.area_ratio) AS "Et_BldgNum",
+SUM("Et_AssetValue" * b.area_ratio) AS "Et_AssetValue",
+SUM("Et_PopNight" * b.area_ratio) AS "Et_PopNight",
+SUM(eq_shaking_score_abs * b.area_ratio) AS "eq_shaking_score_abs",
+SUM(eq_shaking_score_rel * b.area_ratio) AS "eq_shaking_score_rel",
+SUM(fld_priority_score_abs * b.area_ratio) AS "fld_priority_score_abs",
+SUM(fld_priority_score_rel * b.area_ratio) AS "fld_priority_score_rel",
+SUM(wildfire_priority_score_abs * b.area_ratio) AS "wildfire_priority_score_abs",
+SUM(wildfire_priority_score_rel * b.area_ratio) AS "wildfire_priority_score_rel",
+SUM(cy_priority_score_abs * b.area_ratio) AS "cy_priority_score_abs",
+SUM(cy_priority_score_rel * b.area_ratio) AS "cy_priority_score_rel",
+SUM(total_abs_score * b.area_ratio) AS "total_abs_score",
+SUM(total_rel_score * b.area_ratio) AS "total_rel_score",
+c.geom
+
+FROM results_nhsl_hazard_threat.nhsl_hazard_threat_indicators_s_tbl a
+LEFT JOIN boundaries."SAUID_HexGrid_50km_intersect_unclipped" b ON a."Sauid" = b.sauid
+LEFT JOIN boundaries."HexGrid_50km_unclipped" c ON b.gridid_50 = c.gridid_50
+GROUP BY b.gridid_50,c.geom;
+
+-- 100km
+DROP VIEW IF EXISTS results_nhsl_hazard_threat.nhsl_hazard_threat_prioritization_hexgrid_100km_uc;
+CREATE VIEW results_nhsl_hazard_threat.nhsl_hazard_threat_prioritization_hexgrid_100km_uc AS 
+SELECT
+b.gridid_100,
+SUM("Et_BldgNum" * b.area_ratio) AS "Et_BldgNum",
+SUM("Et_AssetValue" * b.area_ratio) AS "Et_AssetValue",
+SUM("Et_PopNight" * b.area_ratio) AS "Et_PopNight",
+SUM(eq_shaking_score_abs * b.area_ratio) AS "eq_shaking_score_abs",
+SUM(eq_shaking_score_rel * b.area_ratio) AS "eq_shaking_score_rel",
+SUM(fld_priority_score_abs * b.area_ratio) AS "fld_priority_score_abs",
+SUM(fld_priority_score_rel * b.area_ratio) AS "fld_priority_score_rel",
+SUM(wildfire_priority_score_abs * b.area_ratio) AS "wildfire_priority_score_abs",
+SUM(wildfire_priority_score_rel * b.area_ratio) AS "wildfire_priority_score_rel",
+SUM(cy_priority_score_abs * b.area_ratio) AS "cy_priority_score_abs",
+SUM(cy_priority_score_rel * b.area_ratio) AS "cy_priority_score_rel",
+SUM(total_abs_score * b.area_ratio) AS "total_abs_score",
+SUM(total_rel_score * b.area_ratio) AS "total_rel_score",
+c.geom
+
+FROM results_nhsl_hazard_threat.nhsl_hazard_threat_indicators_s_tbl a
+LEFT JOIN boundaries."SAUID_HexGrid_100km_intersect_unclipped" b ON a."Sauid" = b.sauid
+LEFT JOIN boundaries."HexGrid_100km_unclipped" c ON b.gridid_100 = c.gridid_100
+GROUP BY b.gridid_100,c.geom;