quarta-feira, 6 de abril de 2011

Correspondência entre zona UTM e SRID

Estou com um DVD do IBGE onde os shapes estão na projeção UTM / SAD69. Preciso indicar o SRID ao rodar o programa shp2pgsql pra conseguir transformar em coordenadas posteriormente.
A tabela abaixo tem o SRID (da EPSG) e a zona correspondente. Listei apenas as do Brasil.

29118 - zone 18N
29119 - zone 19N
29120 - zone 20N
29121 - zone 21N
29122 - zone 22N
29177 - zone 17S
29178 - zone 18S
29179 - zone 19S
29180 - zone 20S
29181 - zone 21S
29182 - zone 22S
29183 - zone 23S
29184 - zone 24S
29185 - zone 25S

A seguinte SQL faz a conversão com base na fórmula que coloquei no post Em qual zona UTM está um determinado ponto?:
select ibge, case when latitudese < 0 then 29160+trunc((longitudes + 186) / 6) else 29100+trunc((longitudes + 186) / 6) end srid from mapas.shp_munic

Outras correspondências eu busco aqui:
http://svn.osgeo.org/metacrs/csmap/trunk/CsMapDev/Dictionaries/coordsys.asc