title: Changing the Coordinate Reference System (CRS) on Census Tract Data
body: <p><strong>How can I change the  Coordinate Reference System (CRS) on Census Tract Data to match the Coordinate Reference System (CRS) on the UCLA Congressional District Shapefile?</strong></p>
  <p>I am working with the <code>sf</code> package and Census Tract data called <code>Census_70.1</code>. The data looks like the following:</p>
  <pre><code>Simple feature collection with 73669 features and 15 fields
  geometry type:  MULTIPOLYGON
  dimension:      XY
  bbox:           xmin: -7115208 ymin: -1685018 xmax: 3321632 ymax: 4591848
  epsg (SRID):    NA
  proj4string:    +proj=aea +lat_1=29.5 +lat_2=45.5 +lat_0=37.5 +lon_0=-96 +x_0=0 +y_0=0 +datum=NAD83 +units=m +no_defs
  
  geometry
  1  MULTIPOLYGON (((777980.6 -6...
  2  MULTIPOLYGON (((797299.6 -6...
  3  MULTIPOLYGON (((781590.6 -7...
  4  MULTIPOLYGON (((788040.7 -7...
  5  MULTIPOLYGON (((782468.1 -7...
  6  MULTIPOLYGON (((783828.9 -7...
  7  MULTIPOLYGON (((776624.7 -7...
  8  MULTIPOLYGON (((775274.6 -7...
  9  MULTIPOLYGON (((791123.2 -7...
  10 MULTIPOLYGON (((797096.8 -7...
  </code></pre>
  <p>I need to change the CRS to match the following congressional district shapefile:</p>
  <pre><code>Simple feature collection with 436 features and 19 fields (with 1 geometry empty)
  geometry type:  MULTIPOLYGON
  dimension:      XY
  bbox:           xmin: -179.1473 ymin: 18.91383 xmax: 179.7785 ymax: 71.35256
  epsg (SRID):    4269
  proj4string:    +proj=longlat +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +no_defs
  </code></pre>
  <p>I've attempted to use the following code to re-set the CRS using the following code:</p>
  <pre><code>Census_70.2 &lt;- st_transform(st_set_crs(Census_70.1, &quot;+proj=longlat +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +no_defs&quot;), crs = 4269, check = T) 
  </code></pre>
  <p>However, while the <code>proj4string</code> changes on <code>Census_70.2</code>, the geometry stays the same.</p>
  <pre><code>Simple feature collection with 34492 features and 9 fields
  geometry type:  MULTIPOLYGON
  dimension:      XY
  bbox:           xmin: -7115713 ymin: -1295867 xmax: 2146924 ymax: 3476641
  epsg (SRID):    NA
  proj4string:    +proj=longlat +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +no_defs
  
                           geometry
  1  MULTIPOLYGON (((777980.6 -6...
  2  MULTIPOLYGON (((797299.6 -6...
  3  MULTIPOLYGON (((781590.6 -7...
  4  MULTIPOLYGON (((788040.7 -7...
  5  MULTIPOLYGON (((782468.1 -7...
  6  MULTIPOLYGON (((783828.9 -7...
  7  MULTIPOLYGON (((776624.7 -7...
  8  MULTIPOLYGON (((775274.6 -7...
  9  MULTIPOLYGON (((791123.2 -7...
  10 MULTIPOLYGON (((797096.8 -7...
  
  </code></pre>
  <p>Furthermore, by running the transformation code, it distorts the map from its original format in <code>Census_70.1</code>
  <a href="https://i.stack.imgur.com/B4Rhh.png" rel="nofollow noreferrer"><img src="https://i.stack.imgur.com/B4Rhh.png" alt="enter image description here" /></a></p>
  <p>To this distortion in <code>Census_70.2</code>
  <a href="https://i.stack.imgur.com/CVmE0.png" rel="nofollow noreferrer"><img src="https://i.stack.imgur.com/CVmE0.png" alt="enter image description here" /></a></p>
  <p>What am I doing incorrectly?</p>