I have two rasters, and one of them has one more column than the other. So I want to remove the last column in the larger one.
How can I change the raster width in place without having to do  a with open(f, 'r') and then another with open(f, 'w')?
I tried using r+, but when trying to do ds.width = 890, I got AttributeError: attribute 'width' of 'rasterio._base.DatasetBase' objects is not writable

with rio.open(f, 'r+') as ds:

        print(ds.profile) 

 

        ds_arr = ds.read(1) # Do I need to actually modify the np array here and then overwrite the file later?

        print(ds_arr.shape)

        

        profile = ds.profile

        

        width = profile['width']

        height = profile['height']

        

        if width == 891:

            ds.width = 890 # error

            
Profile: {'driver': 'GTiff', 'dtype': 'int16', 'nodata': -1.0, 'width': 891, 'height': 3081, 'count': 1, 'crs': CRS.from_epsg(32629), 'transform': Affine(15.0, 0.0, 232053.930561042,

       0.0, -15.0, 1438436.08069127), 'tiled': False, 'interleave': 'band'}