From raw images and footprint shapefiles to overlapping sections only in R

Question

I have two *.tif images, which have no projection but for which I know the footprint (available as shapefiles). I am looking for a way to extract the part of each images tha