Prepare new indicators

orbisgis · Feb 6, 2024 · 3728ad8 · 3728ad8
1 parent eb9de18
commit 3728ad8
Show file tree

Hide file tree

Showing 2 changed files with 65 additions and 24 deletions.
diff --git a/geoindicators/src/main/groovy/org/orbisgis/geoclimate/geoindicators/LCZIndicators.groovy b/geoindicators/src/main/groovy/org/orbisgis/geoclimate/geoindicators/LCZIndicators.groovy
@@ -212,29 +212,59 @@ String grid_distances(JdbcDataSource datasource, String sprawl, String grid_indi
 /**
  * @author Erwan Bocher (CNRS)
  */
-String scaling_grid(JdbcDataSource datasource, String rsu_lcz) {
+String scaling_grid(JdbcDataSource datasource, String grid_indicators) {
 
-   datasource.execute("""
-    DROP TABLE IF EXISTS GRID_100, GRID_300, GRID_900;
-    CREATE TABLE  GRID_100 AS SELECT * FROM ST_MakeGrid('$rsu_lcz', 100,100);
-    CREATE TABLE  GRID_300 AS SELECT * FROM ST_MakeGrid('$rsu_lcz', 300,300);
-    CREATE TABLE  GRID_900 AS SELECT * FROM ST_MakeGrid('$rsu_lcz', 900,900);
+    if (!grid_indicators) {
+        error("No grid_indicators table to aggregate the LCZ fraction")
+        return
+    }
+    def gridCols = datasource.getColumnNames(grid_indicators)
+
+    def lcz_columns_cool_areas = ["LCZ_PRIMARY_101", "LCZ_PRIMARY_102",
+                             "LCZ_PRIMARY_103", "LCZ_PRIMARY_104","LCZ_PRIMARY_107"]
+    def lcz_fractions = gridCols.intersect(lcz_columns_cool_areas)
+
+    def grid_scaling_indices = postfix("grid_scaling_indices")
+    def grid_fraction_lod_1 = postfix("grid_fraction_lod_1")
+    def grid_fraction_lod_2 = postfix("grid_fraction_lod_2")
+    datasource.execute("""DROP TABLE IF EXISTS $grid_scaling_indices,$grid_fraction_lod_1,$grid_fraction_lod_2;
+    CREATE TABLE $grid_scaling_indices as SELECT the_geom,
+    ID_GRID, ID_COL as ID_COL_LOD_0, ID_ROW as ID_ROW_LOD_0, (CAST (ID_COL/3 AS INT)+1) AS ID_COL_LOD_1,
+    (CAST (ID_ROW/3 AS INT)+1) AS ID_ROW_LOD_1,
+    (CAST (ID_COL/9 AS INT)+1) AS ID_COL_LOD_2,
+    (CAST (ID_ROW/9 AS INT)+1) AS ID_ROW_LOD_2,
+    ${lcz_fractions.join(",")}, ${lcz_fractions.join("+")} AS SUM_LCZ_COOL FROM $grid_indicators;    
+    CREATE TABLE $grid_fraction_lod_1 as SELECT ID_COL_LOD_1, ID_ROW_LOD_1, SUM(SUM_LCZ_COOL) AS LOD_1_SUM_LCZ_COOL, 
+    ST_UNION(ST_ACCUM(the_geom)) as the_geom, count(*) as count
+    FROM $grid_scaling_indices GROUP BY  ID_COL_LOD_1, ID_ROW_LOD_1;
+    CREATE TABLE $grid_fraction_lod_2 as SELECT ID_COL_LOD_2, ID_ROW_LOD_2, SUM(SUM_LCZ_COOL) AS LOD_2_SUM_LCZ_COOL,
+    ST_UNION(ST_ACCUM(the_geom)) as the_geom,count(*) as count
+    FROM $grid_scaling_indices GROUP BY  ID_COL_LOD_2, ID_ROW_LOD_2;
     """.toString())
+    datasource.createIndex(grid_scaling_indices, "ID_COL_LOD_1")
+    datasource.createIndex(grid_scaling_indices, "ID_ROW_LOD_1")
+
+    datasource.save(grid_fraction_lod_1, "/tmp/grid_lod1.geojson", true)
+    datasource.save(grid_fraction_lod_2, "/tmp/grid_lod2.geojson", true)
+
+    def grid_fractions = postfix("grid_fractions")
+
+    //"neighbors"
+
+    /*int level =1
+    //Collect neighbors
+
+     datasource.execute( """
+DROP TABLE IF EXISTS neighbors;
+create table neighbors as
+        SELECT the_geom, LOD_${level}_SUM_LCZ_COOL as  LOD_${level}_LCZ_COOL_2 from $grid_fraction_lod_1
+            where ID_COL_LOD_$level=ID_COL_LOD_$level AND ID_ROW_LOD_$level=ID_ROW_LOD_$level-1""".toString())
+
+
+    datasource.save("neighbors", "/tmp/neighbors.geojson", true)*/
+
 
-    datasource.save("GRID_100", "/tmp/grid_100.geojson", true)
-    datasource.save("GRID_300", "/tmp/grid_300.geojson", true)
-    datasource.save("GRID_900", "/tmp/grid_900.geojson", true)
 
-    def values =[:]
-    int row_index= 0
-    datasource.eachRow("select * from grid_100 order by id_row"){row ->
-        if(row_index<=3){
-            row_index++
-        }
-        else {
-            row_index=0
-        }
-    }
     return null
 
 }

diff --git a/...cators/src/test/groovy/org/orbisgis/geoclimate/geoindicators/SprawlIndicatorsTests.groovy b/...cators/src/test/groovy/org/orbisgis/geoclimate/geoindicators/SprawlIndicatorsTests.groovy
@@ -47,14 +47,25 @@ class SprawlIndicatorsTests {
         h2GIS.save(buffering_cool_areas, "/tmp/buffering_cool_areas.geojson", true)*/
 
 
+        String grid_indicators = h2GIS.load("/home/ebocher/Autres/data/geoclimate/uhi_lcz/Dijon/grid_indicators.geojson", true)
 
-        String grid_indicators = h2GIS.load("/home/ebocher/Autres/data/geoclimate/uhi_lcz/Lorient/grid_indicators.geojson", true)
-
-        String scale_grid = Geoindicators.SprawlIndicators.scaling_grid(h2GIS, grid_indicators)
+        //String scale_grid = Geoindicators.SprawlIndicators.scaling_grid(h2GIS, grid_indicators)
 
        // h2GIS.save(scale_grid, "/tmp/scale_grid.csv", true)
 
-        /*
+        //GRID SCALES
+
+        h2GIS.execute("""DROP TABLE IF EXISTS lod_1,lod_2, lod_3;
+        CREATE TABLE lod_1 as select * from ST_MakeGrid('$grid_indicators', 300,300);
+        CREATE TABLE lod_2 as select * from ST_MakeGrid('$grid_indicators', 900,900);
+        CREATE TABLE lod_3 as select * from ST_MakeGrid('$grid_indicators', 2700,2700);""".toString())
+
+        h2GIS.save("lod_1", "/tmp/lod1.geojson", true)
+        h2GIS.save("lod_2", "/tmp/lod2.geojson", true)
+        h2GIS.save("lod_3", "/tmp/lod3.geojson", true)
+
+
+
         def sprawlLayer = Geoindicators.SprawlIndicators.compute_sprawl_areas(h2GIS, grid_indicators)
 
         h2GIS.save(sprawlLayer, "/tmp/sprawl_areas.geojson", true)
@@ -84,7 +95,7 @@ class SprawlIndicatorsTests {
 
         String cool_areas_distances = Geoindicators.SprawlIndicators.grid_distances(h2GIS, cool_areas_inverse, grid_indicators)
 
-        h2GIS.save(cool_areas_distances, "/tmp/cool_areas_distances.geojson", true)*/
+        h2GIS.save(cool_areas_distances, "/tmp/cool_areas_distances.geojson", true)
 
     }