Documentação

Usuário Discussão:Peter/Entrega2024-06grid: mudanças entre as edições

Página do usuário Discussão

Usuário Discussão:Peter/Entrega2024-06grid (ver código-fonte)

Edição das 23h11min de 30 de junho de 2024

3 057 bytes adicionados ,  30 junho
m
→‎tentativa 3
(10 revisões intermediárias pelo mesmo usuário não estão sendo mostradas)
Linha 177: Linha 177:
; -- NADA DE MENSAGEM!
; -- NADA DE MENSAGEM!
-- nao roda o loop .. devolver a função e gerar texto e rodar na mão.
-- nao roda o loop .. devolver a função e gerar texto e rodar na mão.
</syntaxhighlight>


---------------============================
=== tentativa 3 ===
Tentando simplesmente devolver a string. Parece que agora funciona:
<syntaxhighlight lang="sql" style="font-size: 80%;">
DROP FUNCTION if exists osmc.grid_generate_all_levels_logistic2
;
CREATE or replace FUNCTION osmc.grid_generate_all_levels_logistic2(
    p_lev_max numeric,
    p_contry text,
p_contry_base integer,
    p_id_bitlen integer,  -- ex. 37
    p_lev0 numeric default 2.5,
    p_blocksize int default 100
) RETURNS text LANGUAGE plpgsql  AS $f$
DECLARE
  tpl text;
  cover_frag text;
  oid_list text[];
  lev numeric;
  s  text  :='';
  ret text :='';
  gg  text  :='geom4326';
BEGIN
tpl := $$
  SELECT cbits as gid_vbit, geom
  FROM osmc.decode_scientific_absolute_geoms_lixo1(
    natcod.parents_to_children_baseh(%1$s::real, %3$L::text[], 16, true, true)::text,
    %2$L,
    %4$s::int
  ) t
$$;
SELECT array_agg(x::text) INTO oid_list
FROM (
  SELECT array_agg(cod_b16h)::text as x
      FROM (
        select cod_b16h, row_number() over() ord
        from grid_br2.poc_caruaru_vw01logistic
        where length(gid_vbit)=grid_generate_all_levels_logistic2.p_id_bitlen
      ) t1
      GROUP BY ord/grid_generate_all_levels_logistic2.p_blocksize
) t2;
FOREACH cover_frag IN ARRAY oid_list LOOP
    s:='';
    FOR lev IN (select x from generate_series(p_lev0, p_lev_max, 2.5) t(x)) LOOP
      IF lev>p_lev0 THEN s := s || E'\n UNION ALL \n'; END IF;
      s := s || format(tpl, lev::text, p_contry, cover_frag, p_contry_base::text);
    END LOOP; -- lev
    s := s|| E'\n  ORDER BY 1';
    ret := ret || E'\n\nINSERT INTO grid_br2.tmp_cells \n'|| s --|| E';\n COMMIT; \n\n';
  END LOOP; -- cover_frag


  RETURN ret;
END;
$f$;


SELECT volat_file_write(  -- com 37 e 500 dá 5Mb
  '/tmp/inserts_logistic2.sql'
  , osmc.grid_generate_all_levels_logistic2( 2.5, 'BR', 18, 27, 2.5, 300 )
);
-- deu certo!
chmod 777
psql < sql
---
select count(*) from grid_br2.tmp_cells; -- 16352
select count(*)
from grid_br2.tmp_cells
where ST_Intersects(geom,(select geom from grid_br2.municipio where isolabel_ext='BR-PE-Caruaru'))
; --14489
select count(*) from grid_br2.tmp_cells t inner join grid_br2.poc_caruaru c on c.gid_vbit=t.gid_vbit; -- 14489
---- -------------------------------------------
SELECT volat_file_write(
  '/tmp/inserts_logistic2.sql'
  , osmc.grid_generate_all_levels_logistic2( 2.5, 'BR', 18, 37, 2.5, 400 )
);  -- 4828425 byts = 4.7M. Em 17 segundos faz 20*12800. Faz 14 mil celulas por segundo. celulas de 8 metros de lado são 64m2.
-- depois de rodar
delete from grid_br2.tmp_cells where not(ST_Intersects(geom,(select geom from grid_br2.municipio where isolabel_ext='BR-PE-Caruaru')) );
-- DELETE 70068
insert into  grid_br2.poc_caruaru SELECT * FROM grid_br2.tmp_cells ORDER BY gid_vbit;
-- INSERT 0 14420108
delete from grid_br2.tmp_cells;
</syntaxhighlight>
Como Caruaru tem da ordem de 920 km2, são 14406456 (14 milhões) de células, portanto 1 mil segundos ou 16 minutos.
=== tentativa 4 ===
Finalmente reescrevendo como procedure:
<syntaxhighlight lang="sql" style="font-size: 80%;">
DROP FUNCTION if exists osmc.grid_generate_all_levels_logistic2
DROP FUNCTION if exists osmc.grid_generate_all_levels_logistic2
;
;
CREATE or replace FUNCTION osmc.grid_generate_all_levels_logistic2(
DROP PROCEDURE if exists osmc.grid_generate_all_levels_logistic2
;
CREATE PROCEDURE osmc.grid_generate_all_levels_logistic2(
     p_lev_max numeric, p_contry text,
     p_lev_max numeric, p_contry text,
p_contry_base integer,
p_contry_base integer,
     p_id_bitlen integer,  -- ex. 37
     p_id_bitlen integer,  -- ex. 37
     p_lev0 numeric default 2.5
     p_lev0 numeric default 2.5
) RETURNS text LANGUAGE plpgsql  AS $f$
) LANGUAGE plpgsql  AS $f$
DECLARE
DECLARE
   tpl text;
   tpl text;
Linha 206: Linha 292:
  $$;
  $$;
  RAISE NOTICE 'entrei...';
  RAISE NOTICE 'entrei...';
select array_agg(x) FROM (
select array_agg(x::text) INTO oid_list
SELECT array_agg(cod_b16h)::text as x  
FROM (
    FROM (
  SELECT array_agg(cod_b16h)::text as x  
      select cod_b16h, row_number() over() ord
      FROM (
      from grid_br2.poc_caruaru_vw01logistic
        select cod_b16h, row_number() over() ord
      where length(gid_vbit)=p_id_bitlen  -- ! control here
        from grid_br2.poc_caruaru_vw01logistic
    ) t1
        where length(gid_vbit)=27
    GROUP BY ord/1000
      ) t1
  ) t2     INTO oid_list;
      GROUP BY ord/100
  ) t2;
  FOREACH cover_frag IN ARRAY oid_list LOOP
  FOREACH cover_frag IN ARRAY oid_list LOOP
     s:='';
     s:='';
Linha 223: Linha 310:
     END LOOP; -- lev
     END LOOP; -- lev
     s := s|| E'\n  ORDER BY 1';
     s := s|| E'\n  ORDER BY 1';
     s := E'INSERT INTO grid_b2.tmp_cells \n' || s || E';\n COMMIT;';
     s := E'INSERT INTO grid_br2.tmp_cells \n' || s || E';\n COMMIT;';
     RAISE NOTICE 'RODANDO %', s;
     RAISE NOTICE '.. RODANDO';
     s := s || E'\n\n';
    PERFORM s;
     s := '';
   END LOOP; -- cover_frag
   END LOOP; -- cover_frag
   RAISE NOTICE 'sai...';
   RAISE NOTICE 'sai...';
  RETURN s;
END;
END;
$f$;
$f$;


select  osmc.grid_generate_all_levels_logistic2( 2.5, 'BR', 37, 18 );
CALL osmc.grid_generate_all_levels_logistic2( 2.5, 'BR', 27, 18 );
 
</syntaxhighlight>
</syntaxhighlight>
Peter
2 402

edições

Disponível em “https://wiki.addressforall.org/doc/Especial:MobileDiff/4141...4179”