Resample using rasterio

resample_rasterio.py

import rasterio
from rasterio import Affine
from rasterio.warp import reproject, Resampling


with rasterio.open('path/to/band') as src:
    scale = 2 # scaling amount of data so a scale of 2 would double the pixel size
    aff = src.transform
    
    # adjust the new affine transform to new scale
    newaff = Affine(aff.a / scale, aff.b, aff.c,
                aff.d, aff.e / scale, aff.f)

    meta = src.profile
    meta.update(affine = newaff)
    
    with rasterio.open('path/to/output/destination', 'w', **meta) as dst: 
         reproject(source=rasterio.band(src, 1), 
             destination=rasterio.band(dst, 1),
                  resampling = Resampling.bilinear)