Q: There is a way to obtain the value of the different DQ flag contained only in the postage stamp?
>>> from astropy.io import fits
>>> from reftools.interpretdq import DQParser, ImageDQ
>>> dqdata = fits.getdata('jbt7a3k7q_flc.fits', 3) # DQ, 1
>>> dqparser = DQParser.from_instrument('ACS')
>>> acsdq = ImageDQ(dqdata[:100, :100], dqparser=dqparser)
>>> acsdq.interpret_all()
Parsing DQ flag(s)...
Done!
Run time: 0.006 s
N_FLAGGED: 366/10000 (3.660%)
FLAG=1 : 0 (0.000%)
FLAG=2 : 0 (0.000%)
FLAG=4 : 0 (0.000%)
FLAG=8 : 0 (0.000%)
FLAG=16 : 111 (1.110%)
FLAG=32 : 21 (0.210%)
FLAG=64 : 115 (1.150%)
FLAG=128 : 0 (0.000%)
FLAG=256 : 0 (0.000%)
FLAG=512 : 100 (1.000%)
FLAG=1024 : 0 (0.000%)
FLAG=2048 : 0 (0.000%)
FLAG=4096 : 49 (0.490%)
FLAG=8192 : 0 (0.000%)
FLAG=16384: 0 (0.000%)
FLAG=32768: 0 (0.000%)
>>> acsdq.data
array([[ 512, 512, 512, ..., 512, 512, 512],
[ 0, 0, 0, ..., 4144, 0, 0],
[ 0, 0, 0, ..., 0, 0, 0],
...,
[ 0, 0, 0, ..., 0, 0, 0],
[ 0, 0, 0, ..., 0, 0, 0],
[ 0, 0, 0, ..., 0, 0, 0]], dtype=int16)
>>> acsdq.data.shape
(100, 100)