I am trying to make a map of Chicago in Stata. I want this map to include borders for neighborhoods as well as major roads. I used the shp2dta command to convert the shp files into dta files.

Here's where I got the SHP file for the roads --> https://catalog.data.gov/dataset/tig...ased-shapefile

Here's where I got the SHP file for the neighborhood borders --> https://hub-cookcountyil.opendata.ar...e14a2eb234af_2

When I ran the shp2dta command on the roads shp file, I got a dta coordinate file with _X values ranging from -88.3 to -87.52 and _Y values ranging from 41.4 to 42.2. These values seem to align pretty closely with the latitude and longitude of roads in the city of chicago (assuming negative values for the longtiude means "Western Hemisphere").

However, when I ran the shp2dta command for the neighbhorhood borders shp file, I got a dta coordinate file with _X values ranging from 1,003,172 to 1,205,632 and _Y values ranging from 1,749,790 to 1,998,855. I have no idea what these values correspond to.

I was able to use spmap to make a map of the boundaries and a map of the roads. However, when I try to overlay the two, I get a (very small) map in the top-right of the output screen. I cannot zoom in enough to see if this map is a smaller version of the map I am trying to make.

I suspect that the problem lies with the different coordinate systems for the two files, but I am not sure. Any advice would be greatly appreciated.