2 384
edições
osmc:Metodologia/Algoritmo SQL/Issues: mudanças entre as edições
osmc:Metodologia/Algoritmo SQL/Issues (ver código-fonte)
Edição das 23h44min de 6 de junho de 2024
, 6 junho→Issue 02 - Função estimadora de cobertura interior e quadrado envolvente
| Linha 126: | Linha 126: | ||
A interseção da geometria da célula com o polígono do país tem a sua área indicada pela coluna ''area_km2''. A coluna ''is_contained'' é um flag, verdadeiro quando a célula tem interseção com as bordas do país. Analisando as áreas do Brasil, o ''flag'' foi falso apenas para os <math>1099512 \ km^2</math> de área, cuja raíz quadada é o tamanho do lado do quadrado, e uma das [[osmc:Metodologia#Projeção_e_cobertura_nacionais|potências de 2 utilizadas]], <math>\sqrt{1099512\ km^2} = 1048576 \ m = 2^{20} \ m = 1048,58 \ km</math>. Analogamente em Camarões (CM), <math>\sqrt{68719\ km^2} = 262143.1 \ m \approx 262144 \ m = 2^{18} \ m = 262.144 \ km</math>. | A interseção da geometria da célula com o polígono do país tem a sua área indicada pela coluna ''area_km2''. A coluna ''is_contained'' é um flag, verdadeiro quando a célula tem interseção com as bordas do país. Analisando as áreas do Brasil, o ''flag'' foi falso apenas para os <math>1099512 \ km^2</math> de área, cuja raíz quadada é o tamanho do lado do quadrado, e uma das [[osmc:Metodologia#Projeção_e_cobertura_nacionais|potências de 2 utilizadas]], <math>\sqrt{1099512\ km^2} = 1048576 \ m = 2^{20} \ m = 1048,58 \ km</math>. Analogamente em Camarões (CM), <math>\sqrt{68719\ km^2} = 262143.1 \ m \approx 262144 \ m = 2^{18} \ m = 262.144 \ km</math>. | ||
== Issue 03 - Melhora em == | |||
Conforme [https://gis.stackexchange.com/q/444441/7505 "How to reduce the clothesline-effect in ST_Transform?"], pode-se minimizar o efeito-varal. A ideia aqui nesta issue é usasar simplificação de curva, para controle e eventual complementação do algoritmo. | |||
<syntaxhighlight lang="sql"> | |||
select hlevel, count(*) n_cells, round(sqrt(max(st_area(geom)))/1000.0) side_km, | |||
round(avg(ST_NPoints(geom))) npts_geom, round(avg(ST_NPoints(geom4326))) npts_geom, | |||
round(avg(ST_NPoints(st_simplify(geom4326,0.00000005)))) npts_simpl | |||
from grid_cm.all_levels | |||
group by 1 order by 1; | |||
</syntaxhighlight> | |||
<pre> | |||
hlevel | n_cells | side_km | npts_geom | npts_geom | npts_simpl | |||
--------+---------+---------+-----------+-----------+------------ | |||
0 | 14 | 262 | 5 | 749 | 293 | |||
0.5 | 28 | 185 | 5 | 769 | 280 | |||
1 | 56 | 131 | 5 | 749 | 152 | |||
1.5 | 112 | 93 | 5 | 769 | 141 | |||
2 | 224 | 66 | 5 | 749 | 77 | |||
2.5 | 448 | 46 | 5 | 769 | 71 | |||
3 | 896 | 33 | 5 | 749 | 39 | |||
3.5 | 1792 | 23 | 5 | 769 | 36 | |||
4 | 3584 | 16 | 5 | 749 | 20 | |||
4.5 | 7168 | 12 | 5 | 769 | 19 | |||
5 | 14336 | 8 | 5 | 749 | 11 | |||
5.5 | 28672 | 6 | 5 | 769 | 11 | |||
</pre> | |||