jwass · December 17, 2015 22:33
diff --git a/converted.geojson b/converted.geojson
diff --git a/jqtrde.py b/jqtrde.py
 import pandas as pd
 import geopandas as gpd

 import fiona.crs
 from shapely.geometry import Point

 # Open the file as a table. Name the two columns as x and y
 df = pd.read_csv('raw.txt', header=None, names=['x', 'y'], sep=' ')
 # Convert to shapely Points and to a GeoDataFrame with epsg 26919
 geoms = df.apply(lambda row: Point(row['x'], row['y']), axis=1)
 gdf = gpd.GeoDataFrame({'geometry': geoms}, crs=fiona.crs.from_epsg(26919))
 # Convert to lon/lat
 gdf.to_crs(epsg=4326, inplace=True)
 # Store GeoJSON
 with open('converted.geojson', 'w') as f:
    f.write(gdf.to_json())
diff --git a/raw.txt b/raw.txt
 613105 4975798
 612517 4975332
 618438 4976478
 618116 4976792
 612186 4974863
 607978 4982238
 617011 4980265
 617490 4980715
 573541 4965659
 570290 4977544
	import pandas as pd
	import geopandas as gpd

	import fiona.crs
	from shapely.geometry import Point

	# Open the file as a table. Name the two columns as x and y
	df = pd.read_csv('raw.txt', header=None, names=['x', 'y'], sep=' ')
	# Convert to shapely Points and to a GeoDataFrame with epsg 26919
	geoms = df.apply(lambda row: Point(row['x'], row['y']), axis=1)
	gdf = gpd.GeoDataFrame({'geometry': geoms}, crs=fiona.crs.from_epsg(26919))
	# Convert to lon/lat
	gdf.to_crs(epsg=4326, inplace=True)
	# Store GeoJSON
	with open('converted.geojson', 'w') as f:
	f.write(gdf.to_json())