import nibabel as nb
import bids
layout = bids.BIDSLayout('/data/bids/ds003130-download/')
imgs = layout.get(extension='.nii.gz', datatype=['fmap', 'func'])
dirs = {"A": {"j": "PA", "j-": "AP"}} # Truncated for simplicity
for bfile in imgs:
img = bfile.get_image()
ornt = ''.join(nb.aff2axcodes(img.affine))
md = bfile.get_metadata()
peaxis = md['PhaseEncodingDirection']
pedir = dirs[ornt[1]][peaxis]
print(f"{bfile.filename:<46} {ornt} {peaxis:<3} {pedir}")