compress.py

import numpy as np
import rasterio

original = rasterio.open("original.tif").read()

# cmd:
#  convert -compress LZW -type TrueColor original.tif lossless.tif
lossless = rasterio.open("lossless.tif").read()

# cmd:
# convert original.tif -type TrueColor lossy.jpg
lossy = rasterio.open("lossy.jpg").read()

print original.shape, lossless.shape, lossy.shape

#passes yay
np.testing.assert_almost_equal(original, lossless)

# fails
np.testing.assert_almost_equal(original, lossy)