Le portail francophone de la géomatique

import pandas as pd
import geopandas as gpd
from shapely.geometry import Point


df = pd.read_excel('output/ZMT_1.xlsx')


def convert_to_decimal(coord):
    if isinstance(coord, str):
        
        parts = coord.split(' ')
        if len(parts) == 2:
            direction = parts[0]  # N, S, E, W
            degrees_minutes = parts[1].split('°')
            degrees = float(degrees_minutes[0])
            minutes = float(degrees_minutes[1][:-1])
            
            
            decimal = degrees + (minutes/60)
            
            if direction in ['S', 'W']:
                decimal = -decimal
            
            return decimal
        else:
            return None
    else:
        return None

df['Latitude_decimal'] = df['Latitude'].apply(convert_to_decimal)
df['Longitude_decimal'] = df['Longitude'].apply(convert_to_decimal)


df['geometry'] = df.apply(lambda row: Point(row['Longitude_decimal'], row['Latitude_decimal']), axis=1)


gdf = gpd.GeoDataFrame(df, geometry='geometry', crs='EPSG:4326')


output_path = 'output/zmt_kd.shp'
gdf.to_file(output_path)

print(f'Fichier shapefile enregistré avec succès à {output_path}')