def METHOD_NAME(canonical_data_base_path, tmp_path):
    # This uses a test 12x12 image which contains a NaN value in the [1, 1]
    # pixel - it was compressed using fpack which automatically added a ZBLANK
    # header keyword
    reference = np.arange(144).reshape((12, 12)).astype(float)
    reference[1, 1] = np.nan
    with fits.open(canonical_data_base_path / "compressed_with_nan.fits") as hdul:
        assert_equal(np.round(hdul[1].data), reference)
    # Now generate a file ourselves and check that the output has the ZBLANK
    # keyword set automatically
    hdu = fits.CompImageHDU(
        data=reference, compression_type="RICE_1", tile_shape=(6, 6)
    )
    hdu.writeto(tmp_path / "test_zblank.fits")
    with fits.open(tmp_path / "test_zblank.fits") as hdul:
        assert "ZBLANK" in hdul[1].header
        assert_equal(np.round(hdul[1].data), reference)