diff --git a/util/removelines.py b/util/removelines.py
index b1c57bf66..3d4b4ede4 100755
--- a/util/removelines.py
+++ b/util/removelines.py
@@ -51,7 +51,7 @@ def removelines(infile, outfile, xcol='X', ycol='Y', ext=1, cut=None, **kwargs):
     #header.add_history('This xylist was filtered by the "removelines.py" program')
     #header.add_history('to remove horizontal and vertical lines of sources')
     #header['REMLINEN'] = (len(x) - len(xc), 'Number of sources removed by "removelines.py"')
-    T.writeto(outfile)
+    T.writeto(outfile, header=T.get_header())
     return 0
 
 def main():
diff --git a/util/uniformize.py b/util/uniformize.py
index 6a06e46fd..7258f3f80 100755
--- a/util/uniformize.py
+++ b/util/uniformize.py
@@ -14,7 +14,7 @@ def uniformize(infile, outfile, n, xcol='X', ycol='Y', ext=1, **kwargs):
     T = fits_table(infile, lower=False)
     if len(T) == 0:
         print('No sources')
-        T.writeto(outfile)
+        T.writeto(outfile, header=T.get_header())
         return
     x = T.get(xcol)
     y = T.get(ycol)
@@ -63,7 +63,7 @@ def uniformize(infile, outfile, n, xcol='X', ycol='Y', ext=1, **kwargs):
     J = np.array(J)
     #header.add_history('This xylist was filtered by the "uniformize.py" program')
     T.cut(J)
-    T.writeto(outfile)
+    T.writeto(outfile, header=T.get_header())
     return 0
 
 def main():