terrapop_models 0.1.0

data/db/migrate/20150729160936_fixingterrapop_reclassify_categorical_raster_to_binary_summariz_v3.rb

@@ -0,0 +1,68 @@
+class FixingterrapopReclassifyCategoricalRasterToBinarySummarizV3 < ActiveRecord::Migration
+  def change
+    sql =<<-SQL
+    CREATE OR REPLACE FUNCTION terrapop_reclassify_categorical_raster_to_binary_summariz_v3(IN sample_geog_lvl_id bigint, IN rasters_id bigint, IN bnd_num integer)
+        RETURNS TABLE(geog_instance bigint, geog_instance_label character varying, pixel_count bigint, binary_area double precision, percent_area double precision, total_area double precision) AS
+        $BODY$
+
+    DECLARE
+        area_raster_id integer;
+        categorical_raster_id integer;
+
+    BEGIN
+        SELECT rv.area_reference_id as area_id into area_raster_id FROM raster_variables rv WHERE rv.id =  rasters_id ;
+        SELECT rv.second_area_reference_id into categorical_raster_id FROM raster_variables rv WHERE rv.id = rasters_id ;
+
+        RETURN QUERY
+
+        WITH lookup AS
+            (
+            SELECT replace(replace(array_agg(classification::text || ':1')::text, '{', ''), '}', '') as exp
+            FROM raster_variables WHERE id IN (
+                    select raster_variable_classifications.mosaic_raster_variable_id 
+                    from raster_variable_classifications
+                    where raster_variable_classifications.raster_variable_id = rasters_id)
+            ), transformation as
+            (
+            SELECT ST_SRID(r.rast) as prj_value
+            FROM rasters r
+            WHERE r.raster_variable_id = categorical_raster_id
+            LIMIT 1
+            ), polygon as
+            (
+            SELECT sgl.id as sample_geog_level_id, gi.id as geog_instance_id, gi.label as geog_instance_label,
+            gi.code as geog_instance_code, ST_Transform(bound.geog::geometry, t.prj_value) as geom
+            FROM transformation t, sample_geog_levels sgl
+            inner join geog_instances gi on sgl.id = gi.sample_geog_level_id
+            inner join boundaries bound on bound.geog_instance_id = gi.id
+            WHERE sgl.id = sample_geog_lvl_id
+            ), r_table as
+            (
+            SELECT p.geog_instance_id, p.geog_instance_label, ST_union(ST_Reclass(ST_Clip(r.rast, bnd_num ,p.geom, TRUE),1,l.exp, '8BUI',0)) as rast
+            FROM lookup l, polygon p inner join rasters r on ST_Intersects(r.rast, p.geom)
+            WHERE r.raster_variable_id = categorical_raster_id
+            GROUP by p.sample_geog_level_id, p.geog_instance_id, p.geog_instance_label, p.geog_instance_code
+            ), a_table as
+            (
+            SELECT p.geog_instance_id, p.geog_instance_label, ST_union(ST_Clip(r.rast, p.geom)) as rast
+            FROM lookup l, polygon p inner join rasters r on ST_Intersects(r.rast, p.geom)
+            WHERE r.raster_variable_id = area_raster_id
+            GROUP by p.sample_geog_level_id, p.geog_instance_id, p.geog_instance_label, p.geog_instance_code
+            ), calc as
+            (
+            SELECT r.geog_instance_id, r.geog_instance_label, (ST_SummaryStats(r.rast)).*, (ST_SummaryStats(a.rast)).sum as total_area,
+            (ST_SummaryStats(ST_MapAlgebra(r.rast, 1, a.rast, 1, '[rast1]*[rast2]', '32BUI'))).sum as binary_area
+            FROM r_table r inner join a_table a on (r.geog_instance_id = a.geog_instance_id)
+            )
+            SELECT c.geog_instance_id, c.geog_instance_label, c.count as pixel_count, c.binary_area,
+            c.binary_area / c.total_area as percent_area, c.total_area as total_area
+            FROM calc c;
+
+    END;
+    $BODY$
+    LANGUAGE plpgsql VOLATILE;
+    SQL
+    
+    execute sql
+  end
+end

data/db/migrate/20150730141036_adjust_jpeg_clip_for_categorical.rb

@@ -0,0 +1,67 @@
+class AdjustJpegClipForCategorical < ActiveRecord::Migration
+  def change
+    
+    sql = <<-SQL
+      CREATE OR REPLACE Function terrapop_jpeg_raster_clip_colormap_v2(sample_geog_lvl_id bigint, rasters_id bigint, raster_bnd integer, colormap text)
+      RETURNS TABLE (tiff bytea) AS
+      $BODY$
+
+        BEGIN
+        RETURN QUERY
+
+        WITH poly_table AS 
+        (
+        SELECT ST_Transform(bound.geog::geometry, (SELECT ST_SRID(r.rast) FROM rasters r WHERE r.raster_variable_id = rasters_id LIMIT 1)) as geom
+          FROM sample_geog_levels sgl
+          inner join geog_instances gi on sgl.id = gi.sample_geog_level_id
+          inner join boundaries bound on bound.geog_instance_id = gi.id
+          where gi.sample_geog_level_id = sample_geog_lvl_id
+        ),
+        new_rast AS
+        (
+          SELECT ST_Union(ST_Clip(r.rast, raster_bnd, p.geom, TRUE)) as rast
+          FROM poly_table p inner join rasters r on ST_Intersects(r.rast, p.geom)
+          where r.raster_variable_id = rasters_id
+        )
+        select ST_AsJPEG(ST_ColorMap(ST_Union(r.rast), 1, colormap), 1) as tiff from new_rast r;
+
+        END;
+      $BODY$
+      LANGUAGE 'plpgsql';
+    SQL
+    
+    sql0 = <<-SQL
+      CREATE OR REPLACE Function terrapop_jpeg_raster_clip_colormap_with_buffer_v1(sample_geog_lvl_id bigint, rasters_id bigint, raster_bnd integer, colormap text)
+      RETURNS TABLE (tiff bytea) AS
+      $BODY$
+
+        BEGIN
+        RETURN QUERY
+
+        WITH poly_table AS 
+        (
+        SELECT ST_Buffer(ST_Transform(bound.geog::geometry, (SELECT ST_SRID(r.rast) FROM rasters r WHERE r.raster_variable_id = rasters_id LIMIT 1)), 0.0000001) as geom
+          FROM sample_geog_levels sgl
+          inner join geog_instances gi on sgl.id = gi.sample_geog_level_id
+          inner join boundaries bound on bound.geog_instance_id = gi.id
+          where gi.sample_geog_level_id = sample_geog_lvl_id
+        ),
+        new_rast AS
+        (
+          SELECT ST_Union(ST_Clip(r.rast, raster_bnd, p.geom, TRUE)) as rast
+          FROM poly_table p inner join rasters r on ST_Intersects(r.rast, p.geom)
+          where r.raster_variable_id = rasters_id
+        )
+        select ST_AsJPEG(ST_ColorMap(ST_Union(r.rast), 1, colormap), 1) as tiff from new_rast r;
+
+        END;
+      $BODY$
+      LANGUAGE 'plpgsql';
+    SQL
+    
+    execute sql
+    
+    execute sql0
+    
+  end
+end

data/db/migrate/20150730153216_ststandardize_terrapop_reclassify_categorical_raster_to_binary_summariz_v3.rb

@@ -0,0 +1,69 @@
+class StstandardizeTerrapopReclassifyCategoricalRasterToBinarySummarizV3 < ActiveRecord::Migration
+  def change
+    
+  sql =<<-SQL
+    CREATE OR REPLACE FUNCTION terrapop_reclassify_categorical_raster_to_binary_summariz_v4(IN sample_geog_lvl_id bigint, IN rasters_id bigint, IN bnd_num integer)
+        RETURNS TABLE(geog_instance_id bigint, geog_instance_label character varying, pixel_count bigint, binary_area double precision, percent_area double precision, total_area double precision) AS
+        $BODY$
+
+    DECLARE
+        area_raster_id integer;
+        categorical_raster_id integer;
+
+    BEGIN
+        SELECT rv.area_reference_id as area_id into area_raster_id FROM raster_variables rv WHERE rv.id =  rasters_id ;
+        SELECT rv.second_area_reference_id into categorical_raster_id FROM raster_variables rv WHERE rv.id = rasters_id ;
+
+        RETURN QUERY
+
+        WITH lookup AS
+            (
+            SELECT replace(replace(array_agg(classification::text || ':1')::text, '{', ''), '}', '') as exp
+            FROM raster_variables WHERE id IN (
+                    select raster_variable_classifications.mosaic_raster_variable_id 
+                    from raster_variable_classifications
+                    where raster_variable_classifications.raster_variable_id = rasters_id)
+            ), transformation as
+            (
+            SELECT ST_SRID(r.rast) as prj_value
+            FROM rasters r
+            WHERE r.raster_variable_id = categorical_raster_id
+            LIMIT 1
+            ), polygon as
+            (
+            SELECT sgl.id as sample_geog_level_id, gi.id as geog_instance_id, gi.label as geog_instance_label,
+            gi.code as geog_instance_code, ST_Transform(bound.geog::geometry, t.prj_value) as geom
+            FROM transformation t, sample_geog_levels sgl
+            inner join geog_instances gi on sgl.id = gi.sample_geog_level_id
+            inner join boundaries bound on bound.geog_instance_id = gi.id
+            WHERE sgl.id = sample_geog_lvl_id
+            ), r_table as
+            (
+            SELECT p.geog_instance_id, p.geog_instance_label, ST_union(ST_Reclass(ST_Clip(r.rast, bnd_num ,p.geom, TRUE),1,l.exp, '8BUI',0)) as rast
+            FROM lookup l, polygon p inner join rasters r on ST_Intersects(r.rast, p.geom)
+            WHERE r.raster_variable_id = categorical_raster_id
+            GROUP by p.sample_geog_level_id, p.geog_instance_id, p.geog_instance_label, p.geog_instance_code
+            ), a_table as
+            (
+            SELECT p.geog_instance_id, p.geog_instance_label, ST_union(ST_Clip(r.rast, p.geom)) as rast
+            FROM lookup l, polygon p inner join rasters r on ST_Intersects(r.rast, p.geom)
+            WHERE r.raster_variable_id = area_raster_id
+            GROUP by p.sample_geog_level_id, p.geog_instance_id, p.geog_instance_label, p.geog_instance_code
+            ), calc as
+            (
+            SELECT r.geog_instance_id, r.geog_instance_label, (ST_SummaryStats(r.rast)).*, (ST_SummaryStats(a.rast)).sum as total_area,
+            (ST_SummaryStats(ST_MapAlgebra(r.rast, 1, a.rast, 1, '[rast1]*[rast2]', '32BUI'))).sum as binary_area
+            FROM r_table r inner join a_table a on (r.geog_instance_id = a.geog_instance_id)
+            )
+            SELECT c.geog_instance_id, c.geog_instance_label, c.count as pixel_count, c.binary_area,
+            c.binary_area / c.total_area as percent_area, c.total_area as total_area
+            FROM calc c;
+
+    END;
+    $BODY$
+    LANGUAGE plpgsql VOLATILE;
+    SQL
+    
+    execute sql
+  end
+end

data/db/migrate/20150730163446_add_origin_to_extract_requests.rb

@@ -0,0 +1,5 @@
+class AddOriginToExtractRequests < ActiveRecord::Migration
+  def change
+    add_column :extract_requests, :origin, :string
+  end
+end

data/db/migrate/20150731151304_consistent_return_fieldsterrapop_gli_yield_areal_summarization.rb

@@ -0,0 +1,99 @@
+class ConsistentReturnFieldsterrapopGliYieldArealSummarization < ActiveRecord::Migration
+  def change
+    
+    sql1 =<<-SQL
+      CREATE OR REPLACE Function terrapop_gli_yield_areal_summarization_v2(sample_geog_lvl_id bigint, rasters_id bigint)
+      RETURNS TABLE (geog_instance_id bigint, geog_instance_label character varying, min double precision, max double precision, mean double precision, count bigint)
+          AS 
+          $BODY$
+          BEGIN
+              RETURN QUERY
+                  WITH bin_rast as 
+                  (
+                  SELECT sgl.id as sample_geog_level_id, gi.id as geog_instance_id, gi.label as geog_instance_label,
+                  gi.code as geog_instance_code, ST_Union(ST_Clip(r.rast, bound.geog::geometry)) as rast
+                  FROM sample_geog_levels sgl
+                  inner join geog_instances gi on sgl.id = gi.sample_geog_level_id
+                  inner join boundaries bound on bound.geog_instance_id = gi.id
+                  inner join rasters r on ST_Intersects(r.rast,bound.geog::geometry)
+                  where sgl.id = sample_geog_lvl_id and r.raster_variable_id = rasters_id
+                  group by sgl.id, gi.id, gi.label, gi.code
+                  ), area_ref_rast as
+                  (
+                  SELECT sgl.id as sample_geog_level_id, gi.id as geog_instance_id, gi.label as geog_instance_label,
+                  gi.code as geog_instance_code,ST_Union(ST_Clip(r.rast, bound.geog::geometry)) as rast
+                  FROM sample_geog_levels sgl inner join geog_instances gi on sgl.id = gi.sample_geog_level_id
+                  inner join boundaries bound on bound.geog_instance_id = gi.id
+                  inner join rasters r on ST_Intersects(r.rast,bound.geog::geometry)
+                  where sgl.id = sample_geog_lvl_id and r.raster_variable_id = ( select rv.area_reference_id from raster_variables rv where rv.id = rasters_id limit 1 )
+                  group by sgl.id, gi.id, gi.label, gi.code
+                  )
+                  SELECT bin_rast.geog_instance_id, bin_rast.geog_instance_label,
+                  sum((ST_SummaryStats(bin_rast.rast,1)).min) as min,
+                  sum((ST_SummaryStats(bin_rast.rast,1)).max) as max,
+                  sum((ST_SummaryStats(bin_rast.rast,1)).mean) as mean,
+                  sum((ST_SummaryStats(bin_rast.rast,1)).count)::bigint as yield_cell_count
+                  FROM bin_rast inner join area_ref_rast on ST_intersects(bin_rast.rast, area_ref_rast.rast)
+                  GROUP BY bin_rast.geog_instance_id, bin_rast.geog_instance_label;
+          END;
+          $BODY$
+      LANGUAGE 'plpgsql';
+    SQL
+    
+    sql2 = <<-SQL
+      CREATE OR REPLACE Function terrapop_gli_harvest_areal_summarization_v2(sample_geog_lvl_id bigint, rasters_id bigint)
+      RETURNS TABLE (geog_instance_id bigint, geog_instance_label character varying, percent double precision, harvest_area double precision)
+          AS 
+          $BODY$
+          BEGIN
+              RETURN QUERY
+                  WITH bin_rast as 
+                      (
+                      SELECT sgl.id as sample_geog_level_id, gi.id as geog_instance_id, gi.label as geog_instance_label,
+                      gi.code as geog_instance_code, ST_Union(ST_Clip(r.rast, bound.geog::geometry)) as rast
+                      FROM sample_geog_levels sgl
+                      inner join geog_instances gi on sgl.id = gi.sample_geog_level_id
+                      inner join boundaries bound on bound.geog_instance_id = gi.id
+                      inner join rasters r on ST_Intersects(r.rast,bound.geog::geometry)
+                      where sgl.id = sample_geog_lvl_id and r.raster_variable_id = rasters_id
+                      group by sgl.id, gi.id, gi.label, gi.code
+                      ), zero_rast as
+                      (
+                      SELECT sgl.id as sample_geog_level_id, gi.id as geog_instance_id, gi.label as geog_instance_label,
+                      gi.code as geog_instance_code,
+                      ST_Union(ST_Clip(ST_SetBandNoDataValue(r.rast,0), bound.geog::geometry)) as rast
+                      FROM sample_geog_levels sgl
+                      inner join geog_instances gi on sgl.id = gi.sample_geog_level_id
+                      inner join boundaries bound on bound.geog_instance_id = gi.id
+                      inner join rasters r on ST_Intersects(r.rast,bound.geog::geometry)
+                      where sgl.id = sample_geog_lvl_id and r.raster_variable_id = rasters_id
+                      group by sgl.id, gi.id, gi.label, gi.code 
+                      ), area_ref_rast as
+                      (
+                      SELECT sgl.id as sample_geog_level_id, gi.id as geog_instance_id, gi.label as geog_instance_label,
+                      gi.code as geog_instance_code,ST_Union(ST_Clip(r.rast, bound.geog::geometry)) as rast , sum((ST_SummaryStats(r.rast)).sum) as total_area
+                      FROM sample_geog_levels sgl inner join geog_instances gi on sgl.id = gi.sample_geog_level_id
+                      inner join boundaries bound on bound.geog_instance_id = gi.id
+                      inner join rasters r on ST_Intersects(r.rast,bound.geog::geometry)
+                      where sgl.id = sample_geog_lvl_id and r.raster_variable_id = ( select rv.area_reference_id from raster_variables rv where rv.id = rasters_id limit 1 )
+                      group by sgl.id, gi.id, gi.label, gi.code
+                      ), final_rast as
+                      (
+                      SELECT bin_rast.geog_instance_id, bin_rast.geog_instance_label,
+                      sum((ST_SummaryStats(ST_Intersection(ST_Intersection(bin_rast.rast,zero_rast.rast),area_ref_rast.rast,'band2'),1)).sum) as harvest_area
+                      FROM zero_rast inner join bin_rast on bin_rast.geog_instance_id = zero_rast.geog_instance_id
+                      inner join area_ref_rast on bin_rast.geog_instance_id = area_ref_rast.geog_instance_id
+                      GROUP BY bin_rast.geog_instance_id, bin_rast.geog_instance_label
+                      )
+                      select f.geog_instance_id, f.geog_instance_label, (f.harvest_area/a.total_area) as percent_area, f.harvest_area
+                      from final_rast f inner join area_ref_rast a on f.geog_instance_id = a.geog_instance_id;
+          END;
+          $BODY$
+      LANGUAGE 'plpgsql';
+    SQL
+   
+    execute(sql1)
+    execute(sql2)
+    
+  end
+end

data/db/migrate/20150731192327_add_boundary_only_field_to_terrapop_sample.rb

@@ -0,0 +1,5 @@
+class AddBoundaryOnlyFieldToTerrapopSample < ActiveRecord::Migration
+  def change
+    add_column :terrapop_samples, :boundaries_only, :boolean, default: false
+  end
+end

data/db/migrate/20150803213048_add_function_to_detect_boundaries_that_wrap_globe.rb

@@ -0,0 +1,50 @@
+class AddFunctionToDetectBoundariesThatWrapGlobe < ActiveRecord::Migration
+  def change
+    
+    sql=<<-SQL
+    CREATE OR REPLACE FUNCTION terrapop_boundary_wrap(sample_geog_lvl_id bigint) RETURNS Boolean AS
+    $BODY$
+      DECLARE
+
+        far_west double precision;
+        far_east double precision;
+
+      BEGIN
+
+        with boundbox as
+        (
+        select ST_X((ST_Dumppoints(ST_convexhull(geom))).geom) as coordinates
+          from sample_geog_levels sgl inner join geog_instances gi on sgl.id = gi.sample_geog_level_id inner join boundaries bound on bound.geog_instance_id = gi.id WHERE sgl.id = sample_geog_lvl_id
+        )
+        select min(coordinates)
+        into far_west
+        from boundbox;
+    
+
+        with boundbox as
+        (
+        select ST_X((ST_Dumppoints(ST_convexhull(geom))).geom) as coordinates
+          from sample_geog_levels sgl inner join geog_instances gi on sgl.id = gi.sample_geog_level_id inner join boundaries bound on bound.geog_instance_id = gi.id WHERE sgl.id = sample_geog_lvl_id
+        )
+        select max(coordinates)
+        into far_east
+        from boundbox;
+
+        RAISE NOTICE  ' % ', far_west;
+        RAISE NOTICE  ' % ', far_east;
+
+        IF far_west < -178 AND far_east > 178 THEN
+          RETURN True;
+        ELSE
+          RETURN False;
+        END IF;
+
+      END;
+    $BODY$
+      LANGUAGE plpgsql VOLATILE
+      COST 100;
+    SQL
+    
+    execute sql
+  end
+end

data/db/migrate/20150814154127_add_processing_flag_to_extract_request.rb

@@ -0,0 +1,5 @@
+class AddProcessingFlagToExtractRequest < ActiveRecord::Migration
+  def change
+    add_column :extract_requests, :processing, :boolean, default: false
+  end
+end

data/db/migrate/20150814165605_add_colormapped_categorical_to_binary_jpeg_cookie_cutteringizing_function.rb

@@ -0,0 +1,49 @@
+class AddColormappedCategoricalToBinaryJpegCookieCutteringizingFunction < ActiveRecord::Migration
+  def change
+    sql1 =<<-SQL
+    CREATE OR REPLACE FUNCTION terrapop_categorical_to_binary_as_jpeg_colormap_v1(IN sample_geog_lvl_id bigint, IN rasters_id bigint, IN bnd_num integer DEFAULT 1, colormap TEXT DEFAULT 'greyscale')
+    RETURNS  bytea AS
+    $BODY$
+      DECLARE ret_tiff bytea;
+        BEGIN
+        
+        WITH lookup AS
+        (
+        SELECT replace(replace(array_agg(classification::text || ':1')::text, '{', ''), '}', '') as exp
+        FROM raster_variables WHERE id IN (
+                select raster_variable_classifications.mosaic_raster_variable_id 
+                from raster_variable_classifications
+                where raster_variable_classifications.raster_variable_id = rasters_id)
+        ), cat_rast as
+        (
+        SELECT rv.second_area_reference_id as cat_id
+        FROM raster_variables rv
+        WHERE rv.id = rasters_id
+        ),transformation as
+        (
+        SELECT ST_SRID(r.rast) as prj_value
+        FROM cat_rast c, rasters r
+        WHERE r.raster_variable_id = c.cat_id 
+        LIMIT 1
+        ), polygon as
+        (
+        SELECT sgl.id as sample_geog_level_id, gi.id as geog_instance_id, gi.label as geog_instance_label,
+        gi.code as geog_instance_code, ST_Transform(bound.geog::geometry, t.prj_value) as geom
+        FROM transformation t, sample_geog_levels sgl
+        inner join geog_instances gi on sgl.id = gi.sample_geog_level_id
+        inner join boundaries bound on bound.geog_instance_id = gi.id
+        WHERE sgl.id = sample_geog_lvl_id
+        )
+        SELECT ST_AsJPEG(ST_ColorMap(ST_union(ST_Reclass(ST_Clip(r.rast, bnd_num,p.geom, TRUE),1,l.exp, '8BUI',0))), 1) into ret_tiff
+        FROM lookup l, cat_rast c,  polygon p inner join rasters r on ST_Intersects(r.rast,p.geom)
+        WHERE r.raster_variable_id = (c.cat_id);
+        RETURN ret_tiff;
+        END;
+
+        $BODY$
+    LANGUAGE 'plpgsql';
+    SQL
+        
+    execute(sql1)
+  end
+end

data/db/migrate/20150814183708_around_the_world_wrap_a_raster_function.rb

@@ -0,0 +1,157 @@
+class AroundTheWorldWrapARasterFunction < ActiveRecord::Migration
+  def change
+    sql =<<-SQL
+    CREATE OR REPLACE FUNCTION terrapop_wrap_global_raster_v1(IN sample_geog_lvl_id bigint, IN rast_id bigint, IN rast_band_num integer)
+      RETURNS SETOF raster AS
+    $BODY$
+    DECLARE
+
+        one_raster text := 'hello';
+        query text := '';
+
+        BEGIN
+
+        WITH lookup AS
+        (
+        SELECT replace(replace(array_agg(classification::text || ':1')::text, '{', ''), '}', '') as exp
+        FROM raster_variables WHERE id IN (
+                select raster_variable_classifications.mosaic_raster_variable_id
+                from raster_variable_classifications
+                where raster_variable_classifications.raster_variable_id =  rast_id )
+        ), 
+        cat_rast as
+        (
+        SELECT rv.second_area_reference_id as cat_id
+        FROM raster_variables rv
+        WHERE rv.id =  rast_id
+        )
+        SELECT r_table_schema || '.' || table_name as tablename
+        INTO one_raster
+        FROM cat_rast inner join new_rasters on cat_rast.cat_id = new_rasters.raster_variable_id ;
+
+        RAISE NOTICE  ' % ', one_raster;
+
+        query := $$ WITH lookup AS
+        (
+        SELECT replace(replace(array_agg(classification::text || ':1')::text, '{', ''), '}', '') as exp
+          FROM raster_variables WHERE id IN (
+                select raster_variable_classifications.mosaic_raster_variable_id
+                from raster_variable_classifications
+                where raster_variable_classifications.raster_variable_id = $$ || rast_id || $$)
+        )
+        , transformation as
+        (
+        SELECT ST_SRID(r.rast) as prj_value
+        FROM $$ || one_raster || $$ r
+        LIMIT 1
+        ), polygon as
+        (
+        SELECT sgl.id as sample_geog_level_id, gi.id as geog_instance_id, gi.label as geog_instance_label,
+        gi.code as geog_instance_code, ST_Buffer(ST_Transform(bound.geog::geometry, t.prj_value), 0.0001) as geom
+        FROM transformation t, sample_geog_levels sgl
+        inner join geog_instances gi on sgl.id = gi.sample_geog_level_id
+        inner join boundaries bound on bound.geog_instance_id = gi.id
+        WHERE sgl.id = $$ || sample_geog_lvl_id || $$
+        ), polydump as
+        (
+        SELECT (ST_dump(geom)).path[1] as id, (ST_dump(geom)).geom
+        FROM polygon
+        ), polygon_x as
+        (
+        Select p.id, ST_X(ST_Centroid(p.geom)) as x, p.geom
+        FROM polydump p 
+        ), polygon_west as
+        (
+        select id, x, geom
+        from polygon_x
+        where x > 0
+        ), 
+        polygon_east as
+        (
+        select id, x, geom
+        from polygon_x
+        where x < 0
+        ), raster_west as
+        (
+        SELECT ST_union(ST_Reclass(ST_Clip(r.rast, 1,p.geom, TRUE),1,l.exp, '8BUI',0)) as rast
+        FROM lookup l, polygon_west p inner join $$ || one_raster || $$ r on ST_Intersects(r.rast,p.geom)
+        ), raster_east as
+        (
+        SELECT ST_union(ST_Reclass(ST_Clip(r.rast, 1,p.geom, TRUE),1,l.exp, '8BUI',0)) as rast
+        FROM lookup l, polygon_east p inner join $$ || one_raster || $$ r on ST_Intersects(r.rast,p.geom) 
+        )
+        select ST_MapAlgebra(e.rast, 1, w.rast,1, '[rast1.val] + [rast2.val]' ,'32BUI', 'UNION') as rast
+        from raster_east e, raster_west w $$ ;
+
+        RAISE NOTICE  ' % ', query;
+
+        RETURN QUERY execute query;
+
+    END;
+    $BODY$
+      LANGUAGE plpgsql VOLATILE
+      COST 100
+      ROWS 1000;
+    SQL
+    
+    execute(sql)
+    
+    sql1 =<<-SQL
+    CREATE OR REPLACE FUNCTION terrapop_categorical_to_binary_as_jpeg_colormap_v2(IN sample_geog_lvl_id bigint, IN rasters_id bigint, IN bnd_num integer DEFAULT 1, colormap TEXT DEFAULT 'greyscale')
+    RETURNS  bytea AS
+    $BODY$
+      DECLARE ret_tiff bytea;
+        BEGIN
+        
+        WITH lookup AS
+        (
+        SELECT replace(replace(array_agg(classification::text || ':1')::text, '{', ''), '}', '') as exp
+        FROM raster_variables WHERE id IN (
+                select raster_variable_classifications.mosaic_raster_variable_id 
+                from raster_variable_classifications
+                where raster_variable_classifications.raster_variable_id = rasters_id)
+        ), cat_rast as
+        (
+        SELECT rv.second_area_reference_id as cat_id
+        FROM raster_variables rv
+        WHERE rv.id = rasters_id
+        ),transformation as
+        (
+        SELECT ST_SRID(r.rast) as prj_value
+        FROM cat_rast c, rasters r
+        WHERE r.raster_variable_id = c.cat_id 
+        LIMIT 1
+        ), polygon as
+        (
+        SELECT sgl.id as sample_geog_level_id, gi.id as geog_instance_id, gi.label as geog_instance_label,
+        gi.code as geog_instance_code, ST_Transform(bound.geog::geometry, t.prj_value) as geom
+        FROM transformation t, sample_geog_levels sgl
+        inner join geog_instances gi on sgl.id = gi.sample_geog_level_id
+        inner join boundaries bound on bound.geog_instance_id = gi.id
+        WHERE sgl.id = sample_geog_lvl_id
+        ),
+        rast_rast as 
+        (
+          SELECT ST_union(ST_Reclass(ST_Clip(r.rast, bnd_num,p.geom, TRUE),1,l.exp, '8BUI',0)) AS rast 
+          FROM lookup l, cat_rast c,  polygon p inner join rasters r on ST_Intersects(r.rast,p.geom)
+          WHERE r.raster_variable_id = (c.cat_id)
+        ),
+        rast_x as 
+        (
+          SELECT ST_Width(rast) AS x FROM rast_rast
+        ),
+        rast_y as 
+        (
+          SELECT ST_Height(rast) AS y FROM rast_rast
+        )
+        SELECT ST_AsJPEG(ST_ColorMap(ST_Resample(rr.rast, x.x, y.y), 1, 'fire'), 1, 100) into ret_tiff FROM rast_rast rr, rast_x x, rast_y y;
+        RETURN ret_tiff;
+        END;
+
+        $BODY$
+    LANGUAGE 'plpgsql';
+    SQL
+        
+    execute(sql1)
+  end
+end