lat / lon para utm para lat / lon é extremamente defeituoso, como pod
Eu tentei o seguinte, insira: dados lat / lon, então eu vou calcular uma caixa ao redor, digamos 50 m, então +/- 50 m no valor leste / norte.
Agora reconverti-o para lat / lon e com um script:
http: //robotics.ai.uiuc.edu/~hyoon24/LatLongUTMconversion.p Eu obtenho um resultado que simplesmente não pode ser, já antes é de cerca de 7, depois em torno de 2.
zone, easting, northing = LLtoUTM(23, location.get_lat(), location.get_lon())
topUTM = northing + error
bottomUTM = northing - error
leftUTM = easting - error
rightUTM = easting + error
left, top = UTMtoLL(23, leftUTM, topUTM, zone)
O erro no meu código ou o script pode ter falhas?
Então eu tentei usar pyproj, apenas lat / lon para utm para lat / lon para ver o que acontece
>>> p = pyproj.Proj(proj='utm', zone=32, ellps='WGS84')
>>> p
<pyproj.Proj object at 0x7ff9b8487dd0>
>>> x,y = p(47.9941214, 7.8509671)
>>> print x,y
5159550.36822 1114087.43925
>>> print p(x,y,inverse=True)
(47.971558538495991, 7.8546573140162605)
E aqui não está tão longe quanto no script acima, mas ainda parece suficientemente forte para não poder usá-lo. Por quê? O que posso fazer para obter resultados mais exatos?
EDITAR
Executei test () e todos os testes foram aprovado
no arquivo epsg não existe. O mais próximo que encontrei foi este:
<32632> +proj=utm +zone=32 +ellps=WGS84 +datum=WGS84 +units=m +no_defs <>
no tmerc. Também o que eu precisaria passar o towgs84 como parâmetros? Os acima?