#!/usr/bin/env python3
# -*- coding: utf-8 -*-
"""
Created on Tue Feb 23 16:47:12 2021

@author: brian
"""

from pyproj import Proj, transform
import pandas as pd

pi = Proj(init='epsg:4326')
po = Proj(init='epsg:32617')

d = pd.read_csv('hwf-map.tiff.points.latlon')
for i in d.index:
    x,y = transform(pi,po,d.mapX[i],d.mapY[i])
    d.mapX[i] = x
    d.mapY[i] = y
d.to_csv('hwf-map.tiff.points',index=False)
    #print(x,',',y)