DNGS/Proj: mudanças entre as edições

← Edição anterior

DNGS/Proj (ver código-fonte)

Edição das 11h28min de 20 de novembro de 2023

4 405 bytes adicionados , 20 de novembro de 2023

→‎Proposta ajuste fino

Peter

2 402

edições

@@ Linha 1: / Linha 1: @@
 A seguir as projeções e o seletor de jurisdição, para entregar a projeção correta de cada jurisdição.
+Apesar do requisto de igual-área, o país tem a liberdade de optar por projeções "quasi igual-área", como foi o caso da Colômbia (por estar próxima do Equador a projeção UTM tem distorção de área pequena e satisfatória para o censo).
 == Igual-área ==
@@ Linha 21: / Linha 23: @@
 [[Arquivo:GLANCE-ContinentsGrid1.png|thumb|420px|Grades continentais do GLANCE.]]
-Podemos nos inspirar no [[DGGS/Proj/GLANCE]] da NASA, e sua estratégia de "conectar levemente peças de um mosaico irregular". Usa [[wikipedia:Lambert_Azimuthal_Equal_Area]], como grande parte das projeções DGGS. Para definir as projeções da sua cobertura L0 usa um simples ''template'', de modo que podemos fazer o mesmo por país, mudando apenas os centros (''longitude_of_center'' e ''latitude_of_center'') e calibrando conforme altitude mediana.
+Podemos nos inspirar no [[DGGS/Proj/GLANCE]] da NASA, e sua estratégia de "conectar levemente peças de um mosaico irregular". Usa [[wikipedia:Lambert_Azimuthal_Equal_Area]], como grande parte das projeções DGGS. Para definir as projeções da sua cobertura L0 usa um simples ''template'', de modo que podemos fazer o mesmo por país, mudando apenas os centros (''longitude_of_center'' e ''latitude_of_center'') e [https://gis.stackexchange.com/q/469283/7505 calibrando conforme altitude mediana].
 <pre>
@@ Linha 39: / Linha 41: @@
 Para avaliar melhor os recursos de conexão entre as peças do mosaico, talvez interessante conferir algoritmos mais atualizados do [https://github.com/TUW-GEO/Equi7Grid Equi7Grid] (que não é DGGS  como o GLANCE mas usa mesma estratégia de mosaico L0).
+===Proposta ajuste fino ===
+O ajuste da projeção ao relevo local (altitude mediana dos locais mais povoados e sem picos extremos). Infelizmente o ajuste local requer sistema diferente do WGS84, o que pode gerar "efeitos colaterais", e portanto requer um profundo estudo antes de ser adotado.
+Exemplo do Brasil,  projeção Albers oficial do IBGE:
+* PROJ-string <code>+proj=aea +lat_0=-12 +lon_0=-54 +lat_1=-2 +lat_2=-22 +x_0=5000000 +y_0=10000000 +ellps=WGS84 +units=m</code>.<br/>Nela está implicito o Datum "D_SIRGAS2000", similar ao WGS84 com <code>[6378137,298.2572221009113]</code> (podemos aproximar a segunda para 298.257222101) e <code>UNIT["Degree",0.017453292519943295]</code>.
+* <code>lat_0</code>=latitude_of_origin; <code>lon_0</code>=central_meridian; <code>lat_1</code>=standard_parallel_1; <code>lat_2</code>=standard_parallel_2; <code>x_0</code>=false_easting; <code>y_0</code>=false_northing.
+Projeção com origem [https://gis.stackexchange.com/a/469286/7505 deslocada para 600 metro acima do nível do mar]:
+* PROJ-string <code>+proj=aea +lat_0=-12 +lon_0=-54 +lat_1=-2 +lat_2=-22 +x_0=5000000 +y_0=10000000 +a=6378737 +rf=298.257223563 +towgs84=0,0,0,0,0,0,-100 +units=m +no_defs +type=crs</code>.<br/>Nela está implicito o WGS84 pelo <code>towgs84</code>, e o raio <code>+a=6378737</code> é o mesmo, mas o segundo valor, talvez seja melhor adotar do SIRGAS (298.257222101) ao invés do <code>rf=298.257223563</code> (-298.257222101=-0.0000015).
+* A chave para a elevação está em <code>towgs84</code>, "Assuming a sphere with a radius of 6 million meters, 100 parts per million is 600 meters".
+:<small>PS: as projeções  Albers dos [https://epsg.io/5070 EUA] e Brasil foram definidas com relação ao GRS80 ("is an ellipsoid/spheroid, not a true geodetic datum / geographic coordinate system").  Na prática faz uso do DATUM NAD83, que por sua vez difere apenas uma fração de metros (0.0000015) do WGS84, sendo equivalentes para praticamente 100% das aplicações. Ver [[wikipedia:Geodetic_datum#World_Geodetic_System_1984_(WGS_84)|Geodetic datum na Wikipedia]].</small>
+<small>NOTA sobre ajuste sem impacto: uma alternativa para manter o WGS84 é tratar apenas as aplicações sensíveis. Por exemplo na medida de área de lote, apĺicar  correção do valor de área, com base na mediana das altitudes e ignorando zonas de altitudes extremas e densidade populacional zero (picos e fossas não-habitadas).</small>
+=== Métrica das limitações de uso ===
+Projeções oficiais (ex. Albers do IBGE) já estabeleceram um "uso oficial", independente do erro (distorção de área) da projeção na sua métrica de área de terrenos locais.
+Ainda assim, é importante saber quais outras aplicações podem aderir à mesma projeção: é necessário expressar o erro de medida implícito na projeção, e seu perfil de distorção ao longo do território nacional.
+Projeções não-oficiais precisam ser avaliadas. A metodologia AFAcodes sugere que as alternativas sejam  comparadas entre si, para se decidir qual usar: comparar através da métrica de distorção. <br/>A cada país pode ser necessária também uma prévia seleção das aplicações e suas exigências, como requisito mínimo.
+Dentro da metodologia AFAcodes estabelecemos os seguintes critérios e métricas:
+# Projeções locais de referência: tipicamente projeções municipais (ex. UTM oficial do Plano Diretor) utilizadas para agrimensura ou aplicações imobiliárias, e dados oficiais com respectivos polígonos de lote e medida de área oficial.  Sugere-se avaliar 1 a 3 municípios mais populosos do país, supondo que representem amostra de mais de 50% da população.
+# Comparar área dos lotes calculada pela projeção AFAcodes.
+# Supor que "erro de área válido" é o implícito pela [[Subpavimentação|cobertura interior do lote]], que será a fonte de medida com erro sistemático.
+# ''Fit'' de altitude pela ''altitude [[wikipedia:median|mediana]] das regiões povoadas'' (ignorar picos e vales desabitados), por serem os locais onde haverá maior demanda por múltiplas aplicações.
 == Quasi-igual-área ==