Recompute Feature Matrix according to two GRanges objects

[Gab-23]

Is there a way to recompute a Feature Matrix according to two Granges objects?
For example, I want to consider only those fragments overlapping with set of regions gr1 AND gr2,
and discard other fragments, not overlapping both regions

Thanks in advance!

[timoast]

Hi, I think this would be the same as intersecting the two genomic regions and quantifying that right?