title: GeoPandas Set CRS on Points
body: <p>Given the following GeoDataFrame:</p>
  
  <pre><code>h=pd.DataFrame({'zip':[19152,19047],
                 'Lat':[40.058841,40.202162],
                 'Lon':[-75.042164,-74.924594]})
  crs='none'
  geometry = [Point(xy) for xy in zip(h.Lon, h.Lat)]
  hg = GeoDataFrame(h, crs=crs, geometry=geometry)
  hg
  
         Lat          Lon     zip     geometry
  0   40.058841   -75.042164  19152   POINT (-75.042164 40.058841)
  1   40.202162   -74.924594  19047   POINT (-74.924594 40.202162)
  </code></pre>
  
  <p>I need to set the CRS as I did with another GeoDataFrame (like this):</p>
  
  <pre><code>c=c.to_crs("+init=epsg:3857 +ellps=GRS80 +datum=GGRS87 +units=mi +no_defs")
  </code></pre>
  
  <p>I've tried this: </p>
  
  <pre><code>crs={'init': 'epsg:3857'}
  </code></pre>
  
  <p>and this:</p>
  
  <pre><code>hg=hg.to_crs("+init=epsg:3857 +ellps=GRS80 +datum=GGRS87 +units=mi +no_defs")
  </code></pre>
  
  <p>...but no luck.</p>
  
  <p><strong>Some important notes:</strong></p>
  
  <ol>
  <li><p>The other GeoDataFrame for which the above .to_crs method worked was from a shape file and the geometry column was for polygons, not points.
  Its 'geometry' values looked like this after the .to_crs method was applied:</p>
  
  <p>POLYGON ((-5973.005380655156 3399.646267693398...   and when I try the above with the hg GeoDataFrame, they still look like regular lat/long coordinates.</p></li>
  <li><p>If/when this works out, I'll then concatenate these points with the polygon GeoDataFrame in order to plot both (points on top of polygons). </p></li>
  <li><p>When I try concatenating the GeoDataFrames first before using the .to_crs method, and then I use the method on both the point and polygon rows at once, I get the following error:</p>
  
  <p>ValueError: Cannot transform naive geometries.  Please set a crs on the object first.</p></li>
  </ol>
  
  <p>Thanks in advance!</p>