Solucionado (ver solução)
Solucionado
(ver solução)
2
respostas

Conversão não funciona, muda o sistema de coordenadas mas a projeção continua a mesma...

De fato ele muda o sistema de coordenadas, mas continua com a mesma projeção ortográfica. Exemplo: Tenho um arquivo em WGS84 (Grau, Minuto e segundo), e quero converter em SAD1969 zona 22S, o shape vai para SAD1969 zona 22S, porém a projeção fica em Grau, Minuto e Segundo misturado com metros. Alguém já viu?

2 respostas
solução!

Oii Istefferson, como você está? Espero que esteja bem.

Sinto muito pela demora em obter um retorno.

Fazendo a simulação onde temos um CRS inicial do EPSG:4326/WGS 84 que por padrão possui a projeção lat/long e convertendo para o SAD 69 / UTM Zone 22 S, a projeção retornada é em UTM com a unidade padrão dessa projeção, que é em metros, como mostro no código abaixo:

Conversão inicial para o CRS EPSG:4326/WGS 84:

import geopandas as gpd
from shapely.geometry import Point

d = {'col1': ['name1', 'name2'], 'geometry': [Point(1, 2), Point(2, 1)]}
gdf = gpd.GeoDataFrame(d, crs=4326)
gdf

gdf.crs:

<Geographic 2D CRS: EPSG:4326>
Name: WGS 84
Axis Info [ellipsoidal]:
- Lat[north]: Geodetic latitude (degree)
- Lon[east]: Geodetic longitude (degree)
Area of Use:
- name: World.
- bounds: (-180.0, -90.0, 180.0, 90.0)
Datum: World Geodetic System 1984 ensemble
- Ellipsoid: WGS 84
- Prime Meridian: Greenwich

Resultado:

col1geometry
0name1POINT (1 2)
1name2POINT (2 1)

Conversão para o SAD 69 / UTM Zone 22 S:

gdf = gdf.to_crs(29192)
gdf

gdf.crs:

<Projected CRS: EPSG:29192>
Name: SAD69 / UTM zone 22S
Axis Info [cartesian]:
- E[east]: Easting (metre)
- N[north]: Northing (metre)
Area of Use:
- name: Brazil - onshore and offshore northern and southern hemispheres between 54°W and 48°W. In remainder of South America - between 54°W and 48°W, southern hemisphere, onshore.
- bounds: (-54.0, -35.71, -47.99, 7.04)
Coordinate Operation:
- name: UTM zone 22S
- method: Transverse Mercator
Datum: South American Datum 1969
- Ellipsoid: GRS 1967 Modified
- Prime Meridian: Greenwich

Resultado:

col1geometry
0name1POINT (7300226.213447465 10360879.91808068)
1name2POINT (7490067.727935915 10184793.57499825)

Caso possa disponibilizar seu código, base de dados e um print de como está a projeção, poderei fazer uma análise mais assertiva sobre o seu resultado.

Fico no aguardo e a disposição para quaisquer dúvidas.

Grande abraço!

Muito Obrigado!!