include tqdm and reset projection

measprocess/geospatial/process.py

@@ -41,7 +41,7 @@ def project_onto_streets(
     street_series = street_series.to_crs(epsg)
 
     projected = []
-    for point in point_series:
+    for point in tqdm(point_series):
         street_ind = street_series.distance(point).argmin()
         projected.append(
             street_series.loc[street_ind].interpolate(
@@ -76,9 +76,9 @@ def project_onto_streets(
         plt.show()
 
     # Reset Street and Point Series CRS
-    street_series.to_crs(street_original_crs)
-    point_series.to_crs(point_original_crs)
-    projected.to_crs(point_original_crs)
+    street_series = street_series.to_crs(street_original_crs)
+    point_series = point_series.to_crs(point_original_crs)
+    projected = projected.to_crs(point_original_crs)
 
     return projected, deviation_projection