transformações epsg:9603

Na tabela seguinte sintetizam-se os códigos de operações de transformação com 3 parâmetros registadas na base de dados EPSG para Portugal Continental. 

Para cada transformação indica-se também o ano de publicação dos parâmetros de transformação pelo IGP ou, nos casos aplicáveis, se se trata uma transformação "derivada", i.e., se se considera aplicáveis os valores de parâmetros estimados para um datum de destino diferente.

Todos os valores encontrados se baseiam em ajustes de transformações de Molodensky publicados pelo IGP (embora alguns já não estejam actualmente disponíveis no respectivo site). 
Como descrito anteriormente, as operações de transformação estão associadas ao método Geocentric Translations (geog2D domain) com o código epsg:9603.

  

	PARA ►	ETRS89	WGS 84
DE ▼	epsg	4258	4326
Datum 73	4274	epsg:5036 [ano 2009]	epsg:1983 [ano 2009]
Lisbon	4207	epsg:5038 [ano 2009]	epsg:1984 [ano 1998]
ED50	4230	epsg:5040 [copiada de espg:1985]	epsg:1985 [ano 1998]
Lisbon 1890	4666	epsg:5039 [copiada de espg:1986]	epsg:1986 [ano 1998]

Para referência futura, a tabela seguinte lista os valores dos parâmetros para cada transformação:

Parâmetro ▼			Código da operação de transformação ▼
Código	Nome	Unidade	5036	1983	5038	1984	5040	1985	5039	1986
8605	X-axis translation	m	-223.150	-223.237	-303.861	-304.046	-87.987	-87.987	508.088	508.088
8606	Y-axis translation	m	110.132	110.193	-60.693	-60.576	-108.639	-108.639	-191.042	-191.042
8607	Z-axis translation	m	36.711	36.649	103.607	103.640	-121.593	-121.593	565.223	565.223

Não se considera necessária qualquer actualização ou correcção à base de dados EPSG, dado que:

1. as transformações com 3 parâmetros são menos exactas que as transformações com 7 parâmetros; 
2. existem transformações com 7 parâmetros disponíveis para todas as combinações de datum de origem vs. datum de destino; 
3. a generalidade das aplicações SIG permite a utilização de transformações com 7 parâmetros. 

inconsistências nas transformações EPSG para PT

Conclusões da análise das (instâncias de) transformações registadas na base de dados do EPSG para Portugal Continental:

1. Estão registadas transformações cujos parâmetros foram determinados ou publicados pelo IGP (directamente ou através do Eurogeographics) nos anos de 1998, 2000, 2001 e 2009.
2. Nem todas as transformações cujos parâmetros foram determinados ou publicados pelo IGP (directamente ou através do Eurogeographics) nos anos de 1998, 2000, 2001 e 2009 se encontram registadas.
   Por exemplo, não estão registados os actuais parâmetros da transformação de Bursa-Wolf de Datum Lisboa (Hayford) para ETRS89.
   
   
3. As transformações de 7 parâmetros são registadas por vezes segundo a convenção Position Vector (epsg:9606), por vezes segundo a convenção Coordinate Frame Rotation (epsg:9607), e noutros casos encontram-se replicadas nas duas notações.
4. Num caso estranho, a transformação epsg:1792 (segundo o método epsg:9606) é assinalada como deprecated mas a transformação epsg:1992 (segundo o método epsg:9607) ainda é válida. Contudo, os valores dos parâmetros são os basicamente os mesmos e e são válidos (trata-se da transformação de Datum 73 para ETRS89 publicada pelo IGP em 2001). Admito que a epsg:1792 tenha sido retirada porque os valores não estão registados com o mesmo número de casas decimais.
   
   
5. Por vezes (mas nem sempre), uma transformação originalmente ajustada para WGS84 é "duplicada" noutra transformação para ETRS89. A inversa também ocorre. 
6. Não estão registadas transformações entre data locais.

Para tornar a base de dados consistente:

1. Registar TODAS as transformações para WGS84 ou para ETRS89 cujos parâmetros tenham sido determinados ou publicados pelo IGP (directamente ou através do Eurogeographics) nos anos de 1998, 2000, 2001 e 2009.
   Mesmo que um determinado ajuste não seja o melhor ou mais recente, convém que permaneça um registo formal da sua existência.
   O conceito de deprecated no EPSG só se aplica a registos com erros formais ou materiais. Nas observações são assinalados textualmente os casos em que um registo substitui (supersedes) ou é substuído por outro.
   
2. No caso de transformações de 7 parâmetros, garantir que são registadas segundo a convenção Position Vector (epsg:9606).
   O que não invalida que se mantenham os registos que já existem e que seguem a convenção epsg:9607.
   
3. No caso de transformações para ETRF89, garantir que são "duplicadas" para WGS84.
   O objectivo é simplificar a publicação/partilhar de informação geográfica em aplicações de visualização que não requeiram grande exactidão (que também não seria garantida por este tipo de transformações...)
4. Se só estiver registada/publicada uma transformação para WGS84, garantir que é duplicada para ETRF89.
   Neste caso, o objectivo é garantir que é possível documentar a transformação utilizada para transformar dados em sistemas obsoletos para ETRF89.

Posto isto, o passo seguinte é listar todos as transformações que já estão no EPSG e as que faltam.

bursa-wolf

No labirinto das designações alternativas para métodos de transformação...

A designação "transformação de Bursa-Wolf" é utilizada para transformações de Helmert com 7 parâmetros: 3 parâmetros de translação geocêntrica, 3 parâmetros de rotação e um factor de escala.
Existem duas notações diferentes para os parâmetros de rotação, em função do que se convencione ser o sentido positivo de rotação:

• Dextrógiro - usado no método designado por Position Vector (comum na Europa, aconselhado pela ISO 19111, utilizado pelo CRS-EU, utilizado pelo IGP) e correspondente aos métodos epsg:9606 e epsg:1033 
• Levógiro - usado no método designado por Coordinate Frame ou Coordinate Frame Rotation (comum nos EUA e Austrália) e correspondente aos métodos epsg:9607 e epsg:1032.
  [esta convenção também é usada no método de Molodensky-Badekas].

Ver, por exemplo:

• Guidance Note Number 7, part 2, v.42, pp. 121 a 124 (OGP)
• Equation-based methods (ESRI)

molodensky

No labirinto das designações alternativas para métodos de transformação...

Sob a designação genérica de "transformação de Molodensky" encontram-se 3 métodos distintos (cf. Deakin (2004:2)):

It is now common in the literature to describe three Molodensky transformations:
(i) The Molodensky-Badekas transformation: a seven-parameter conformal transformation (or similarity transformation) [...]
(ii) The standard Molodensky transformation: a five-parameter transformation [...]
(iii) The abridged Molodensky transformation: a modified version of the standard Molodensky transformation [...]. The abridged Molodensky transformation equations do not contain the ellipsoidal heights of points to be transformed.

As transformações que interessam mais (neste caso particular)  são as referidas nas alíneas (ii) e (iii) e que serão discutidas mais abaixo.

Transformação de Molodensky-Badekas

Relativamente à transformação geocêntrica de Molodensky-Badekas (epsg:1034) vide Guidance Note Number 7, part 2, v.42, pp. 124:

To eliminate high correlation between the translations and rotations in the derivation of parameter values for the Helmert transformation methods, instead of the rotations being derived about the geocentric coordinate reference system origin they may be derived at a location within the points used in the determination. Three additional parameters, the coordinates of the rotation point, are then required, making 10 parameters in total.

Ou seja, para além dos 7-parâmetros de uma transformação conforme (a.k.a. Helmert ou Bursa-Wolf), é necessário dispôr das 3 coordenadas geocêntricas do ponto de avaliação da transformação (Xp, Yp, Zp), fornecidas no CRS de origem.
Visto de outra forma, a transformação epsg:1032 (Coordinate Frame Rotation) é um caso particular da transformação de Molodensky-Badekas em que o ponto de avaliação é a origem do CRS geocêntrico.

Na base de dados EPSG - tal como para as transformações de 7-parâmetros -  também a transformação de Molodensky-Badekas pode ser declarada sob a forma duma concatenação de operações:

	Molodensky-Badekas (geog2D domain)	9636
1	Geographic 2D to Geographic 3D	9659
2	Geographic 3D to Geocentric	9602
3	Molodensky-Badekas (geocentric domain)	1034
4	Geocentric to Geographic 3D	9602
5	Geographic 3D to Geographic 2D	9659
	Molodensky-Badekas (geog3D domain)	1039
1	Geographic 3D to Geocentric	9602
2	Molodensky-Badekas (geocentric domain)	1034
3	Geocentric to Geographic 3D	9602

Transformação de Molodensky standard e abridged

Passando agora às transformações não-geocêntricas (Guidance Note Number 7, part 2, v.42, :128):

As an alternative to the computation of the new latitude, longitude and ellipsoid height by concatenation of three operations (geographic 3D to geocentric + geocentric to geocentric + geocentric to geographic 3D), the changes in these coordinates may be derived directly as geographic  coordinate offsets through formulas derived by Molodensky [...]

Deakin (2004) apresenta a derivação das fórmulas para as duas variantes:

• epsg:9604 - Molodensky [ver também a página do Instituto Geográfico Português]
• epsg:9605 - Abridged Molodensky [ver também Guidance Note Number 7, part 2, v.42, :128]

Estas transformações requerem 5 parâmetros, 2 dos quais podem ser obtidos a partir das definições dos elipsóides de origem e de destino, designadamente a diferença dos semi-eixos maiores e a diferença de achatamento dos elipsóides. Na base de dados EPSG, que já contem as definições dos vários elipsóides, basta indicar os 3 parâmetros restantes (dϕ, dλ, dh).


Translação geocêntrica

As transformações de Molodensky (epsg:9604 e epsg:9605) assumem que os eixos cartesianos dos dois elipsóides são paralelos.
Formalmente, a transformação de Molodensky standard é equivalente à translação geocêntrica (epsg:1031). Ou seja, desde que um determinado software de transformação de coordenadas realize a concatenação de operações correspondente ao método epsg:9603, os parâmetros de transformação são os mesmos.

Em termos numéricos, existem discrepâncias (que podem ser devidas a arredondamentos nos passos de cálculo intermédios). Segue-se uma experiência simples.
Testou-se a transformação de epsg:27493 (Datum 73/Hayford-Gauss) para epsg:3763 (ETRS89/PT-TM06), usando os parâmetros de transformação do IGP sobre uma amostra de 781 vértices geodésicos com coordenadas conhecidas nos dois CRS. Utilizaram-se os dois métodos - transformação de Molodensky standard e translação geocêntrica - implementados no software ArcGIS 9.3.1. As coordenadas foram armazenadas com precisão de 3 casas decimais.
Os pontos transformados têm sistematicamente uma diferença de 5, 6 ou 7 mm no terreno. É obviamente irrelevante em termos posicionais, mas topologicamente os pontos são diferentes.
É expectável que o mesmo tipo de situação surja em quaisquer implementações diferentes, devido às questões de arredondamentos.

quantas transformações?

quantas transformações estão registadas na base de dados EPSG para os CRS usados em Portugal?

	PARA ►	PTRA08	ETRS89	WGS 84
DE ▼	epsg	5013	4258	4326
Azores Occidental 1939	4182	3		2	5
Azores Central 1948	4183			1	1
Azores Central 1995	4665	6		10	16
Azores Oriental 1940	4184			1	1
Azores Oriental 1995	4664	3		4	7
Porto Santo	4615			1	1
Porto Santo 1995	4663	3		2	5
Selvagem Grande	4616			2	2
PTRA08	5013			1	1
Lisbon 1890	4666		1	2	3
Lisbon 1890 (Lisbon)	4904			2	2
Lisbon	4207		5	4	9
Lisbon (Lisbon)	4803			1	1
Datum 73	4274		6	4	10
ED50	4230		1	4	5
ETRS89	4258			1	1
		15	13	42	70

O quadro anterior baseia-se na versão 7.6 da base de dados EPSG em formato MS Access, inclui transformações deprecated (por motivos que serão clarificados mais tarde) e inclui transformações do tipo:

• epsg:9603 - Geocentric translations (geog2D domain)
• epsg:9606 - Position Vector transformation (geog2D domain)
• epsg:9607 - Coordinate Frame Rotation (geog2D domain)
• operações concatenadas, que se assinalam a itálico  (e que se  reportam apenas às definições de data com referência ao meridiano de Lisboa).

Critérios de criação de transformações

Algumas das transformações foram criadas pelo próprio OGP de acordo com os critérios expressos na pág. 24 do Guidance Note Number 7, v. 2.1, que em seguida se transcrevem:

Many transformations that are available are between a local or national  system and a system which is defined relative to the global standard, the International Earth Rotation Service Terrestrial Reference System (ITRS). At a geodetic level of significance these realisations of the ITRS are discrete. But at a level of accuracy  quite satisfactory for most GIS and oil industry applications  –  nominally submetre accuracy  –  the differences between these ITRS realisations are not significant. The EPSG dataset facilitates the use of these transformations in a hub implementation by:
•  choosing the GPS satellite navigation system's CRS WGS 84, which in recent years has been maintained within 0.1m of the ITRS, as a standard for use as a hub CRS.
•  creating transformations from national or regional realisations of the ITRS to WGS 84 which result in no change in coordinate values, with a transformation accuracy of 1m. 
•  duplicating selected transformations to national or regional realisations of the ITRS as a transformation to WGS 84, with a transformation accuracy  of the lesser of 1 metre or that of the original transformation.
[...]It must be recognised that these transformations may not necessarily be approved by the originator of the original transformation. It should also be recognised that CRSs other than WGS 84 may be more appropriate as a hub. For example, in Europe when there is a need  to retain positions with respect to the EuroAsian tectonic plate, because with time positions on this plate change with respect to a global best fit system such as WGS 84, the adoption of ETRS89 as a hub would be preferable.

Ou seja, para a exactidão requerida pela generalidade das aplicações SIG (e para a exactidão permitida pela utilização de transformações ajustadas de forma a minimizar o erro quadrático médio na área geográfica de aplicação), o OGP considera equivalentes os CRS ETRF89 (epsg:4258), PTRA08 (epsg:5013) e WGS84 (epsg:4326) .

Transformações para WGS84 

Dado que o WGS84 (epsg:4326) é o CRS hub para as concatenações ímplicitas de transformações e/ou conversões, é importante verificar estas transformações.
Seguir-se-á a seguinte chave de decisão:

(if)

Para um dado método e área geográfica de aplicação, existe mais de uma transformação para WGS84?

(then)

Caso exista, seleccionar a transformação mais exacta (i.e. com menor EMQ ou outro indicador agregado de erro, se comparáveis entre registos diferentes), ou mais recente (o indicador de erro pode "piorar" se forem usados mais pontos no ajuste, o que não obsta a que as estimativas dos parâmetros sejam melhores).
* Se a fonte da transformação seleccionada não for uma NMA, manter não obstante a transformação mais exacta reportada por uma NMA.

(elseif)

Para um dado método e área geográfica de aplicação, existe apenas 1 transformação para WGS84?

(then)

Caso exista, verificar se os valores dos parâmetros correspondem aos da transformação mais exacta ou recente reportada por uma NMA (que poderá ser uma transformação para um ITRS).

(else)

Se não existe qualquer transformação, criar uma com base nos valores mais exactos/recentes reportados por uma NMA (que poderá ser uma transformação para um ITRS).

Transformações para ETRS89 e PTRA08

Como resulta da sequência anterior, para poder verificar as transformações para WGS84 é necessário verificar primeiro as transformações para outras realizações do ITRS.

que transformações?

passagem à parte mais complicada das transformações

que tipo de transformações estão registadas na base de dados EPSG para os CRS usados em Portugal?

• epsg:9603 - Geocentric translations (geog2D domain)
• epsg:9606 - Position Vector transformation (geog2D domain)
• epsg:9607 - Coordinate Frame Rotation (geog2D domain)

O nome das transformações é estranho, porque são transformações geocêntricas, mas associadas a sistemas de coordenadas geodésicas 2D. Esta situação é clarificada no Guidance Note Number 7 v. 2: mesmo no caso de CRS geocêntricos, a transformação é sempre declarada em termos de CRS geographic2D.

[...] most  transformations  in the EPSG Dataset  that operate in the geographic or geocentric coordinate domain continue to be associated with the CRS subtype geographic 2D even when they operate between geocentric or geographic 3D CRSs. Applications using this data may need to make appropriate adjustments and use the implicit concatenated operation technique [...]

Ou seja, na prática as transformações 9603, 9606 e 9607 são concatenações de 4 conversões e 1 transformação central, geocêntrica, como se detalha no Guidance Note Number 7, part 2, v.42
As sequências de passos são as seguintes:

	Geocentric Translations (geog2D domain)	9603
1	Geographic 2D to Geographic 3D	9659
2	Geographic 3D to Geocentric	9602
3	Geocentric Translations (geocentric domain)	1031
4	Geocentric to Geographic 3D	9602
5	Geographic 3D to Geographic 2D	9659
	Position Vector transformation (geog2D domain)	9606
1	Geographic 2D to Geographic 3D	9659
2	Geographic 3D to Geocentric	9602
3	Position Vector transformation (geocentric domain)	1033
4	Geocentric to Geographic 3D	9602
5	Geographic 3D to Geographic 2D	9659
	Coordinate Frame Rotation (geog2D domain)	9607
1	Geographic 2D to Geographic 3D	9659
2	Geographic 3D to Geocentric	9602
3	Coordinate Frame Rotation (geocentric domain)	1032
4	Geocentric to Geographic 3D	9602
5	Geographic 3D to Geographic 2D	9659

  

transformações epsg:9601

Esta parte vai ser mais complicada:
quais as transformações registadas na base de dados EPSG para os CRS locais usados em Portugal?

Primeiro as triviais, i.e., as de tipo epsg:9601 (Longitude rotation) aplicadas a CRS geográficos 2D.

Com interesse para PT, estão registadas as:

• epsg:1991

De: epsg:4904 [geographic2D CRS] < epsg:6904 [datum] < epsg:7004 [ellipsoid]

Para: epsg:4666 [geographic2D CRS] < epsg:6666 [datum] < epsg:7004 [ellipsoid] 

Ou seja, o CRS é o "mesmo", mas na definição do datum epsg:6666  Lisbon 1890, associado ao elipsóide de Bessel, é usado o meridiano de Greenwich como referência (e não o meridiano de Lisboa).

Já não me lembro em que software era, mas as definições por referência ao meridiano de Lisboa, embora formalmente válidas como no datum epsg:6904  Lisbon 1890 (Lisbon), originavam erros nas transformações on-the-fly.

Por isso, mais vale prevenir e usar só as definições baseadas no datum epsg:6666.

• epsg:1759

De: epsg:4803 [geographic2D CRS] < epsg:6803 [datum] < epsg:7022 [ellipsoid]

Para: epsg:4207 [geographic2D CRS] < epsg:6207 [datum] < epsg:7022 [ellipsoid] 

Ou seja, o CRS é o "mesmo", mas na definição do datum epsg:6207  Lisbon 1937, associado ao elipsóide de Internacional 1924 (a.k.a. Hayford), é usado o meridiano de Greenwich como referência (e não o meridiano de Lisboa).

Como no caso anterior, mais vale simplificar e usar só as definições baseadas no datum epsg:6207.

Nota: existia uma transformação epsg:1259, que foi descontinuada devido a um erro no registo do valor do parâmetro associado.