# Converting Shapefile points to Lat/Long?

Quick Overview of the Project: I have 9000 records of individuals, each with a coordinate pair in a csv. I also have a shapefile/folder(dbf, prj, sbn, sbx, shp, shx) of municipalities within the range of coordinates. I'd like to know which municipality(shape) in the shapefile corresponds to each record in the csv

Problem: When I look at the shapefile in python, the grid isn't lat/long. I subscribed to a trial version of arcGIS to be able to load up the coordinates and look at them, and they map correctly, so I must just be missing something. Attached is my code and the map printed by Python.

How can I transpose my csv coordinates to the map or vice versa?

``````import shapefile as shp
import matplotlib.pyplot as plt
from shapely.geometry import Point

import csv

n = 0
m = 1

coorDict = {}
coor = ''
fileName = 'TestLatLong'

sfRec = sf.records()
plt.figure()
for shape in sf.shapeRecords():
x = [i[0] for i in shape.shape.points[:]]
y = [i[1] for i in shape.shape.points[:]]
plt.plot(x,y)
plt.show()
plt.figure()

n = 0
#Then create a dictionary of lat/longs from 'TestLatLong.csv'
``````

The output of all shape records and data in the first shape record above For comparison, all the data should be roughly in a rectangle from (40.881321, -73.926644) to (41.442483, -73.542874), and as you can see this map isn't.