DNGS/Proj
A seguir as projeções e o seletor de jurisdição, para entregar a projeção correta de cada jurisdição.
Igual-área
Projeções nacionais:
- Albers Brasil: https://github.com/AddressForAll/pg_pubLib-v1/blob/main/src/pubLib05pgis-extraSRID.sql
- Albers EUA
- Albers Canadá
- ...
Projeções continentais:
- Qualquer DNGS é válida, mas aqui o foco é na baixa complexidade, reusando projeções tradicionais e prevendo um ajuste linear conforme altitude mediana de cada país.
- Albers Africa: https://epsg.io/102022
- Lambert cilindrica equatorial: a matematica mais simples https://wiki.addressforall.org/doc/DGGS/Proj/rHEALPix
- Europa: wikipedia:European grid adota wikipedia:Lambert azimuthal equal-area projection com centro em 52° N, 10° E.
- Glance... padrão continental da NASA para fotos de satélite.
PS: o ajuste nacional por fator de escala não altera o geocódigo, mas o ajuste por planos oblíquos altera a área e portanto invalidam geocódigos continentais.
Proposta Glance
Em 2018 o projeto measures-glance/glance-grids (inspirado no primeiro paper do Equi7 Equi7Grid) ... Usa wikipedia:Lambert_Azimuthal_Equal_Area
Ver template Python abaixo e resultado em https://github.com/measures-glance/glance-grids/blob/master/data/crs_wkt/GLANCE_GRIDS_V01_EU_CRS.wkt
PROJCS["BU MEaSUREs Lambert Azimuthal Equal Area - {continent} - {version}", GEOGCS["GCS_WGS_1984", DATUM["D_WGS_1984", SPHEROID["WGS_1984",6378137.0,298.257223563]], PRIMEM["Greenwich",0.0], UNIT["degree",0.0174532925199433]], PROJECTION["Lambert_Azimuthal_Equal_Area"], PARAMETER["false_easting",0.0], PARAMETER["false_northing",0.0], PARAMETER["longitude_of_center",{lon_of_center}], PARAMETER["latitude_of_center",{lat_of_center}], UNIT["meter",1.0]]
Quasi-igual-área
Nacional:
- UTM Colômbia: https://github.com/AddressForAll/pg_pubLib-v1/blob/main/src/pubLib05pgis-extraSRID.sql
Continental:
- Equi7 ver wikipedia:Azimuthal_Equidistant. Ver https://cartography.tuwien.ac.at/eurocarto/wp-content/uploads/2015/09/3_6_ppt.pdf e https://doi.org/10.1016/j.cageo.2014.07.005 . Ver https://github.com/TUW-GEO/Equi7Grid
Seletor de jurisdições
O padrão DNGS também prevê as convenções necessárias para se decidir qual país usar, ou seja, selecionar com rapidez e baixo custo de CPU a jurisdição presente em cada área do globo. Outra característica é que a projeção do país pode receber uma correção valor de área, com base na mediana das altitudes e ignorando zonas extremas com densidade populacional zero (picos e fossas não-habitadas).
A decisão sobre qual jurisdição usar depende do contexto (já dado a priori) ou de uma decisão algorítmica baseada nos BBOXes de cada país. Abaixo o exemplo dos vizinhos Brasil (BR), Colômbia (CO), Equador (EC) e Uruguai (UY). A interseção de fato é apenas um pequeno retângulo, onde destacamos em colorido a área de decisão baseada nas fronteiras. Nas demais áreas há um único país por BBOX, de modo que a decisão (baseada em BBOXes apenas) é simples e rápida.
Uma vez identificado o país, a sua projeção é adotada, e a grade pode ser evidenciada localmente.
A BBOX original de cada país é decomposta em BBOXes "puras" e "de fronteira". As puras permitem rápida decisão (alta performance), enquanto as de fronteira requerem avaliação de pertinência ao polígono de fronteira (menor performance). Abaixo o resultado da decomposição quando analisando apenas Brasil e Colômbia.