It looks like a similar question has been asked before, and there is some information here:

https://dicomiseasy.blogspot.com/2013/06/getting-oriented-using-image-plane..html

I have a case where it would be nice to be able to dynamically extract:

1. Table Tilt of MR scanner.
2. Plane of instance / series relative to the body.

The sequences in question are actually 'resl' sequences from an Esoate MR scanner. Based on another reference I tried this python script, and at least for the instances that I have it seems correct, but I suspect there is more to it than that with the body position and table tilt ?
_______

#!/usr/bin/env python
import pydicom
import numpy as np # python -m pip install numpy
import os

"""
This function takes IOP of an image and returns its plane (Sagittal, Coronal, Axial)

"""

# Change the Working Directory to the Scripts Working Directory

cwd = os.getcwd()
# Print the current working directory
print("Current working directory: {0}".format(cwd))

abspath = os.path.abspath(__file__)
dname = os.path.dirname(abspath)
os.chdir(dname)
print("Current working directory: {0}".format( os.getcwd()))

def file_plane(IOP):
IOP_round = [round(x) for x in IOP]
plane = np.cross(IOP_round[0:3], IOP_round[3:6])
plane = [abs(x) for x in plane]
if plane[0] == 1:
return "Sagittal"
elif plane[1] == 1:
return "Coronal"
elif plane[2] == 1:
return "Axial"

directory = 'testfiles'
# iterate over files in
# that directory
for filename in os.listdir(directory):
f = os.path.join(directory, filename)
# checking if it is a file
if os.path.isfile(f) and f.endswith('.dcm'):
print('Processing file: '+f)

a=pydicom.read_file(f)
IOP = a.ImageOrientationPatient
plane = file_plane(IOP)
print(f + ': ' + plane +'\n')

This is probably a better solution ?

Determining MR Image Planes, from Kaggle

https://www.kaggle.com/code/davidbroberts/determining-mr-image-planes

Stephen D. Scotti

This little Python method also seems to work with the samples that I have.

I am curious to know if there is a Dicom Tag for the table tilt / elevation for MRI scanners that can be supine to upright ?

called with get_image_plane(ImageOrientationPatient)

def get_image_plane(loc):

row_x = loc[0]
row_y = loc[1]
row_z = loc[2]
col_x = loc[3]
col_y = loc[4]
col_z = loc[5]

# Calculate the normal vector to the image slice
normal_x = row_y * col_z - row_z * col_y
normal_y = row_z * col_x - row_x * col_z
normal_z = row_x * col_y - row_y * col_x

normal = [abs(normal_x), abs(normal_y), abs(normal_z)]

max_component = max(normal)
if max_component == normal[0]:
return "Sagittal "+str(max_component)
elif max_component == normal[1]:
return "Coronal "+str(max_component)
else:
return "Axial "+str(max_component)

Stephen D. Scotti