AstroNuclPhysics ® Nuclear Physics-Astrophysics-Cosmology OSTNUCLINE
Nuclear physics and physics of ionizing radiation Radionuclide scintigraphy

MATHEMATICAL ANALYSIS AND COMPREHENSIVE COMPUTER EVALUATION OF SCINTIGRAPHIC EXAMINATIONS


3.1. Radionuclide ventriculography

Radionuclide ventriculography is used for a detailed examination of the dynamics of the cardiac cycle, both by visual assessment of the motility of the heart walls , as well as by quantitative analysis and in the calculation of a number of parameters of central hemodynamics. Of the quantitative parameters, it is mainly the ejection fraction of the left (or right) ventricle, end-diastolic, pulse and residual volume , cardiac output, ejection and filling rates , time intervals of significant phases of the heart revolution. To assess the motility of the heart walls, contours of the heart walls in end-diastole and end-systole, parametric images of heart rate volume and paradoxical movement, Fourier images of phase and amplitude. Regional phase and amplitude analysis can be performed . If we perform the study in the LAO projection in connection with the study in the LL projection, we will get an overview of the motility of all segments of the left ventricular wall with this combination.

The program uses a completely automatic (mathematical) delimitation of the areas of interest of the heart chamber and background, which improves the accuracy and reproducibility of the quantification of hemodynamic parameters.

Data storage
After premedication with inactive Sn-pyrophosphate (to label erythrocytes
99mTc in vivo), the patient is given approximately 400-800 MBq of sodium pertechnetate 99m Tc. Connect the electrodes of the ECG-cardiomonitor connected to the acquisition computer and set the scintillation camera detector above the precordium in the LAO projection so that the left ventricle lies in the right middle part of the field of view. Select the zoom scale so that the left ventricle fills as much of the field of view as possible. If we want to quantify the right ventricle, it is necessary that this is displayed with a reserve in the left part of the field of view.

Recommended storage mode:

64 x 64 matrix , 16 bit, Study type: ECG -gated

Zoom approx. 1.5-2 ´
(otherwise individually according to the size of the heart and the field of view of the camera)

LAO projection , or and projection LL

Preset time 10 min. Pulse preselection 5,000,000

Preset number of cycles 700, Cycle elimination interval 10%

Before starting the study, care should be taken to stabilize the heart rhythm so that the average and current (last) length of the cycle are close; otherwise, cycles would be incorrectly excluded. Therefore, when preparing the study, we recommend switching to the monitoring storage mode only after the heart rhythm has stabilized.
Then we start our own storage of the phase dynamic study of the cardiac cycle synchronized by the ECG signal. If there is a significant change in heart rate during study storage, the study should be restarted and restarted with a new stabilized
rhythm.
The methodology of phase dynamic scintigraphy of the cardiac cycle gated by the ECG signal is described in detail in §4.4 "Gated
dynamic scintigraphy " of Chapter 4 "Scintigraphy" of the book "Nuclear physics and ionizing radiation physics ".

Evaluation of the study
After
invoking the required ventriculographic scintigraphic study using the basic OSTNUCLINE system, we start a comprehensive program VENTR - radionuclide ventriculography (its first part).

Images of end-diastole and end-systole

Images of the maximally filled ventricle (in end-diastole) and maximally emptied ventricle (in end-systole) are important both for visual assessment of chamber shape and size, for calculating ventricular volume and mainly for defining areas of interest of the ventricle and tissue background. With proper R-wave ECG synchronization, the initial frames up to about 1/4 cycle (and often several end frames) correspond to end-diastole, the end-systole is about 1/3 to 1/2 of the cycle length. However, this temporal location of ED and ES is highly variable individually, so that a flat-rate determination of the end-diastole and end-systolic phase "hard" would not always be accurate (for other systems, however, this is simply solved). In our program VENTR is the definition phase of the end-diastole and end-systole precisely adjusts the en of the actual situation at the particular heart.

The program first summarizes preliminary images of end-diastole and end-systole, calculates a preliminary heart rate image from them and constructs a preliminary volume curve from its most pulsating parts - this curve is always well "modulated", even in severe pathologies of ventricular motility. The upper and lower levels are determined automatically (with the possibility of manual shift) on this curve for the determination of end-diastole and end-systole points. The images corresponding to these curve points are then summed, normalized to each other to form end-diastole (ED) and end-systole images.(ES). This method allows completely accurate creation of these images even in the case of very poor functioning of the heart chamber, as the program finds relatively well pulsating sites on the chamber, from which the curve is sufficiently modulated to reliably determine end-diastole and end-systole.

Defining the regions of interest of the heart chamber and background

On images in end-diastole and end-systole, another important action can be taken for the complex processing of radionuclide ventriculography - delimitation of the areas of interest of the ventricles and background . This delimitation can be performed either manually or automatically in the VENTR program. For better reproducibility of the results, we recommend using automatic delineation , which based on geometric analysis of the radioactivity distribution image in most cases objectively delimits the left ventricular ROI (however, the automatically delimited ROI can be arbitrarily corrected as needed).

Automatic delimitation of the area of interest of the ventricle
The center of gravity of the ventricle is automatically marked (with the possibility of manual correction) on the image of the end-diastole, which is taken as the origin of the polar coordinates in further calculations. In these polar coordinates, an image of the
local concavity of the radioactivity distribution image is then constructed (the calculation is performed simultaneously on the ED image and on the pulse-paradox image), on which the chamber boundaries as places with increased concavity should be sharply expressed. The program then searches this image and constructs the ROI as a junction of places with maximum concavity .

Method of computer delimitation of the area of ??interest of the heart chamber.
a) In the end-diastolic image, the center of gravity T is automatically defined as the origin of the polar coordinates. The center of gravity is to cut the sections (profiles) gradually at different angles.
b), c) Two typical shapes of these profiles. From the profile curves f (x) the function of local concavity K (x) is determined; the maxima of the local concavity lie at the boundary of the chamber.
d) Quantification of the local concavity of the profile using the size of the area between the function f (x) and its secant p (x).
e) The resulting image of local concavity has significant maxima at the border of the ventricle. The heart contour (ROI) can then be defined as the junction of points with maximum local concavity.


The course of automatic delineation of the ROI of the heart chamber
a) The center of gravity of the chamber is made incisions gradually at different angles, on which the local concavity is calculated and stored in the parametric image.
b) In the finished parametric image of local concavity, places of significant maxima are searched for and a line in ROI bits is guided through them - the ROI of the heart chamber (c) is created, which is then smoothed radially and is offered to possibly. manual adjustment.

This automatically defined ROI is then displayed on the background of the end-diastole image and at the same time on the background of the concavity image. This ROI can be modified manually as needed. ROI is further confronted with a combined paradoxical-pulse image (described below), where there is a very important gap separating the ventricle from the vestibule from the vestibule with the opposite dynamics. If it happens that the originally marked ROI partially encroaches into the atrium, it is necessary to delete the overlying part and lead the upper limit of the ROI through this gap between the ventricle and atrium (however, in the VENTR program this happens very rarely, because the automatic delineation procedure image into account). At this stage, we will also know whether the entire possible dyskinetic region is included and we can add it as needed.


The automatically generated ROI of the heart chamber can be definitively corrected not only according to the
parametric image of the local concavity, but also very precisely according to the combined
pulse-paradox image.

Once the ROI of the ventricle in the end-diastole is definitively marked, the ROI of the tissue background for that vent is automatically determined. The calculation is performed on the image of the end-systole and the result is a "crescent-shaped" area, the outer part of which is bounded by the end-diastolic ROI of the ventricle and the inner part by the image of the ventricle in the end-systole. The background ROI is avoided or structures with increased activity. If there is not enough space inside the end-diastolic ROI for a sufficient number of background points (this happens with very low chamber contractility), the program will gain the appropriate number of points from a close neighborhood outside the ROI. Automatically defined ROI background can also be manually modified, but the sentences Shiit is not necessary. The resulting ROI is representative of the average tissue background in the area of the ventricle in the ES-phase, as it closely abuts the ventricle in the ES. Background ROI is important for correction against tissue background of images, but mainly of curves for calculation of EF and other parameters of central hemodynamics (see below).
If we also want to quantify the activity of the
right ventricle , we mark (manually) the ROI of the right ventricle on the appropriate instruction. Again, we mark simultaneously the ED image and the combined heart rate-paradox image using well-visible dark transitions between the ventricle and the atrium. We try to avoid art. pulmonalis. The background ROI for the right ventricle is then defined automatically according to the same criteria as for the left ventricle.

Heart stroke and paradoxical image, heart contours

The images of the heart chamber in end-diastole and end-systole give us some idea of its contractility, but for a better assessment of the shape and mobility of the heart walls, it is appropriate to emphasize these walls in some way ("visibility") - to create contours of the heart chamber . The contours of the heart wall are constructed using mathematically generated isochromatic curves ("contour lines") on the image in end-diastole and end-systole and are displayed together in one image with significant color difference - the contour in end-diastole is blue, the contour in end-systole is red (Fig.3.1.1).
                      

Fig.3.1.1. Distinctive images of the heart revolution constructed by a program for assessing the motility of the heart walls.
Left: Images of end-diastole, end-systole, paradoxical and stroke image.
Right: Contours of the heart wall in end-diastole (blue) and end-systole (red).
This is a markedly pathological case with apical dyskinesis.

Under normal circumstances (with good and isotropic ventricular contraction), the entire end-systolic contour must lie deep and concentrically within the contour in the end-diastole. With an overall reduced ejection fraction, the "depth of immersion" of the end-systolic contour within the end-diastolic contour decreases. When there is a disorder in the local mobility of the heart wall (anisotropic pulsation), the concentration is disturbed - both contours approach in the appropriate place ( hypokinesis ), touch ( akinesia ), or even the end-systolic contour intersects and emerges from the end-diastolic contour ( dyskinesia ) - Fig. 3.1.2.

Fig.3.1.2. Typical shapes of heart wall contours in end-diastole ( ______ ) and in end-systole (...........) at different pathologies of left ventricular contractility.

From the images of end-diastole and end-systole, important parametric images mapping the regional distribution of the dynamics of systolic-diastolic function of the heart are constructed (Fig. 3.1.1):

The stroke image SV (image of Stroke Volume) arises as the difference between the image of end-diastole and end-systole:

SV = ED - ES.

It quantifies us (in relative numbers - pulses), how much blood decreases during the contraction of the ventricles in individual places of the image, ie the local intensity of pulsation of the chamber (in this it is somewhat similar to the amplitude Fourier image). The atria are not visible in the heart rate image (there is filling - negative values), the heart chambers should normally be displayed clearly and relatively homogeneously (a small decrease in the intensity of the heart rate image inside the ventricle at the inlet site is physiological). A dark defect in the left ventricular heartbeat is a sign of reduced contractility - hypokinesia or akinesia . A brighter island separated from the rest of the ventricle by a dark transition on the pulse image is (along with the corresponding finding on the paradoxical image) a manifestation of dyskinesia.

The paradox picture arises as the difference between the pictures in end-systole and in end-diastole:

PAR = ES - ED.

In normal contractility, in the paradoxical image, the chambers must be dark (there is emptying - negative values), the atria is light. The presence of a light "island" in the dark field of the left ventricle is an unmistakable indicator of the paradoxical (opposite, opposite ) movement of the heart wall in this place - dyskinesia (Fig. 3.1.1).

The combined stroke -paradoxical image of SP is given by a suitable composition of stroke and paradoxical image, which can be achieved by using absolute values of differences in the contents of individual elements of ED and ES images, thus turning negative values into positive:

SP = ÷ ED - ES ú = ÷ ES - ED ÷ = ÷ SV - PAR ÷ .

Chambers and atria are clearly shown in this image, which is very similar to an amplitude Fourier image. Very important is the sharp dark transition between the ventricles and atria in this image (it is the transition between ventricular dynamics and opposite atrial dynamics), which allows accurate the regions of interest of chambers to be specified so that they do not interfere with the atria (as described above).

Images of end-diastole and end-systole, contours of the heart wall in ED and ES, stroke and paradox images, are cleary simultaneous displayed in the VENTR program - or we can print them graphically.

Fourier amplitude-phase analysis

Mathematical methods of Fourier analysis can be used for a completely detailed analysis of the dynamics of the heart revolution in individual parts of the heart . From a general point of view, Fourier analysis is based on the theorem, according to which each function f (x) can be expressed as a superposition of harmonic functions A.cos ( w .x) and B.sin ( w .x) with different amplitudes A, B and circular frequencies. w . Cardiac activity is periodic at time t, so we can naturally model it with the function A.cos( w .t + j ) with a fixed circular frequency w = 2 p / T determined by the period T of the cardiac cycle, where A is the pulsation amplitude and j expresses the time (phase) shift of the beginning of the decrease in radioactivity after the arrival of the R-wave, ie the phase difference between the electrical and mechanical systole. It is customary to express the phase shift j in angular units - degrees o , while the whole heart cycle (period T) is assigned 360 o . The relationship between the phase shift j in degrees and the time shift D t in seconds between the electrical and mechanical systole is

D t = T. j / 360 .

The Fourier analysis of radionuclide ventriculography can be imagined in the following steps :

  1. For each point of the image matrix we create a curve f (t) of the time course of the pulsation of the radio indicator - for the matrix 64x64 we get 4096 curves. The shape of each curve depends on the place in the heart from which it comes. Figure 3.1.1 shows four typical curves: a) A well-modulated curve from the normokinetic site of the left or right ventricle. b) The curve from the hypokinetic region is less modulated. c) The curve from the dyskinetic site has pathologically almost opposite dynamics than from the normokinetic chamber. d) The curve from the atrium has a physiologically opposite dynamics than the ventricles.
  2. Each of these curves is approximated by the corresponding cosine function f (t) » A.cos( w .t + j ) with a fixed frequency w = 2 p / T given a period T of the cardiac cycle, whose variable amplitude A expresses the pulsation amplitude at a given place and the phase variable j characterizes the time interval (delay) between the gating pulse (R-wave of the ECG) and the onset of ejection at a given site of the heart. In Fig.3.1.3 we see that curve a) from the normokinetic site of the left ventricle is approximated by a cosine A.cos( w .t + j ) with a relatively high amplitude A and a time (= phase) shift j less than 1/4 of the period T of the cardiac cycle. Curve b) from the hypokinetic region is expressed by a cosine of very small amplitude A. Curve c) from the dyskinetic site corresponds to a cosine function, the beginning of which (and thus phase j ) is delayed by about 1/2 of the period T behind the R-wave. An even larger absolute phase shift j is shown by the cosine function interspersed with a curve from the atrium - more than 1/2 of the period T, which is, however, physiological in the opposite pulsating atrium.
  3. We do this for each element i, j of the image matrix, so we have a total of i ´ j (= 4096 for the matrix 64 ´ 64) of the curves f ij and the cosine functions
    f
    cosij = A ij . cos ( w .t + j ij )
    with different amplitudes A
    ij and phases j ij .
  4. We will now create two new image matrices: in the first one we store the local values of the amplitudes A ij , in the second the local values of the phases j ij . This creates two important parametric images :
    An
    amplitude image , expressing how powerful the pulsation of the heart is at each point of the image. The usual brightness-color modulation scale is used for the amplitude image as for scintigraphic images - the higher the amplitude (better pulsation), the brighter and warmer the color of the given element of the image.
    Phase image
    , expressing the delay , with which the given place of the heart reacts to the synchronizing electric R-wave by mechanical contraction (systole). For the phase image, it is suitable to use a special color modulation scale , which significantly and unambiguously distinguishes normokinetic sites from sites with different degrees of pathological phase shift.
    In Fig.3.1.3 below we see a typical pathological case of Fourier amplitude and phase images, a similar case can be seen in color in Fig.3.1.4. Cosine functions from the normokinetic posterolateral region (curves of type a)) provide high values of amplitude (bright colors in the amplitude image) and low values of phase approx. (0.15-0.25) .T = 55
    ° - 90°, which is assigned a dark-light blue color in the phase image. Hypokinetic regions (type b curves) may also have normal phases (blue in the phase image), but low amplitudes expressed by dark colors in the amplitude image. Dyskinetic area (curves type c)), there apex, usually have somewhat n lower level of amplitude than normokinetic, but the absolute phase value is around about (0.5 to 0.6) .T = 140° - 220°, which corresponds to red to yellow to modulation phase image. Atrial images (type d curves) have a relatively high pulsation amplitude (bright colors in the amplitude image) and especially a very high phase j value (0.6 - 0.7) .T = 220° - 240° corresponding to the yellow to white color in the phase image (Fig. 3.1.4).
Fig.3.1.3. Principle of construction of amplitude and phase Fourier images illustrated on the example of pathological ventriculography.

We will now briefly describe how the Fourier phase analysis of cardiac cycle dynamics is implemented in the VENTR program. When asked whether this analysis should be performed, we will answer in the affirmative especially when, according to the image of heart rate volume or paradoxical image, we suspect a regional heart chamber motility disorder. Before the actual calculation, the time course of the radioactivity in the entire field of view of the camera is first displayed and the end point of the phase analysis is automatically defined on it (with the possibility of manual modification). The calculation then does not include images located near the end of the cycle behind this point, which may be affected by heart rate fluctuations. The result of Fourier analysis at each point of the image are parametric images of the phase and amplitude of cardiac pulsation - Fig.3.1.4.

Fig.3.1.4. Results of regional Fourier phase-amplitude analysis for a severely pathological case of extensive apical dyskinesia. (Physiological Fourier images are seen in the resulting protocol of a normal patient in Fig. 3.1.9.)

A phase distribution histogram is constructed from the phase image , where on the horizontal axis is the phase angle (0 to 360 o ) and on the vertical axis is the number of points (ie the size of the area) with this phase angle of pulsation. The color modulation of the phase image is chosen so that the blue color corresponds roughly to normokinesis (ie phase angles around 60 o ) and the red and yellow colors to asynchrony and dyskinesia with opposite dynamics (phase angle close to 360°) (phase angle approx. 120° and higher); the atriums are displayed in yellow to white. The corresponding verticals on the phase histogram are shown in the same colors - Fig.3.1.4. The amplitude image has the same color modulation as the scintigraphic images - well-pulsating areas with high amplitude are displayed in brighter, darker colors (blue - black). Dyskinetic (red - yellow - white), hypokinetic areas are the area in the amplitude image is expressed by an island (because the surroundings are necessarily akinetic) the brighter the color, the higher the amplitude of dyskinesia. Overall, it can be said that the amplitude image is very similar to the image of heart rate volume, respectively. combined heartbeat - paradoxical image.
If we want to perform
regional quantification phase and amplitude of the cardiac cycle, we will answer the relevant question in the affirmative. Then we gradually mark the respective structures as ROIs on the phase and amplitude images simultaneously (we end each marking with "Mark" and "Continue"). According to the instructions on the display, we first mark the normokinetic region (usually a posterolateral segment), which is then taken as a reference for determining the phase shift and relative amplitude of other regions. After marking each area of interest on the instruction of the program, we will insert its name . For each such marked region, a histogram of the phase distribution in this region is generated and displayed, and the size of the region relative to the whole chamber, the average absolute phase, the relative amplitude, and the phase shift relative to the reference normokinetic region *) are calculated.
*) Only non-zero ones are taken into account in the calculation cells defined by ROI, because zero values ??are created by truncation and have no real meaning. Therefore, when defining the ROI, we do not have to anxiously avoid or. zero notches or islands.
After completing the regional Fourier analysis (answer "no" to the question whether to continue with further quantification), a table of the results of the regional amplitude-phase analysis (Fig. 3.1.4) will be displayed in the lower right part of the display. areas can be printed for documentation (it is especially useful for pathological cases).

Visual evaluation of heart wall motility

The following is a simultaneous display of all significant images (ED, ES, SV, paradox, amplitude and phase images), according to which we can visually assess the regional motility of the heart walls. At the request of the program, we can choose three ways to insert the text of the visual evaluation (Fig. 3.1.5) :

1. Normal finding (standard text)
This method automatically generates a standard formulation of normal evaluation, eg:
"We do not observe a regional ventricular wall motility disorder
on phase scintigraphic images of the cardiac cycle. "
In addition, a standard formulation of normal conclusion is automatically generated:
" Conclusion : Visual evaluation of images of individual phases of the cardiac cycle as well as quantitative analysis of the volume curve indicate good global and local contractility of the walls of the left ventricle.
".

2. Pathological finding (semi-standard text)
In this method, the following text is automatically preselected:
" Phase scintigraphic images show:
segment hypokinesia:
segment akinesia:
segment dyskinesia:
segment asynchrony:
"

In the lines of individual types of disorders, we then enter the names of the segments of the heart wall in which we observe the respective disorder and, if necessary other data (eg about the extent of the fault).

3. Non-standard free text
If none of the previous methods suits the given case, we choose the third method, in which we insert any free text into the free frame using the keyboard.

No normal conclusion is generated in pathological evaluation or non-standard text. All preset texts can be easily edited, supplemented and modified as needed using the mouse and keyboard.

Calculation of absolute ventricular volume

An important parameter in nuclear cardiology is the absolute volume of the heart chamber in milliliters. Absolute volume is important both in itself (it may indicate, for example, dilatation of the heart chamber), and serves as a "conversion factor" between relative and absolute quantities. Knowledge of end-distolic volume allows to create a volume curve of the ventricle from the time course of radioactivity in the chamber, to determine the minute volume of heart and heart rate and residual volume from the value of ejection fraction and heart rate, to determine immediate and maximum ejection and filling rates in milliliters per second.

To calculate the absolute end-diastolic volume of the heart chamber , in principle 4 options can be selected at the request of the program (however, we recommend method 1) :

1. Geometric ratio method
This method calculates the end-diastolic volume of the heart chamber on the basis of geometric analysis of the distribution of radioactivity in the heart chamber image, specifically on the basis of the ratio between the maximum in the ventricular image and the area integral of radioactivity in the whole ROI - Fig.3.1.6. The only external parameter that this method needs - the camera display scale [mm / cell] - is stored in the CAMERA.PAR file. This method was developed at our workplace within the research task of 1985-1990 and tested on a number of phantom measurements and clinical examinations. Based on this experience, we can recommend it to all users as the simplest, most durable and most reproducible. The result of the calculation here practically does not depend on the depth of the ventricular placement and is only slightly dependent on the exact marking of the ROI of the heart chamber.

Fig.3.1.6. Relative geometric-analytical method for determining the absolute volume of the heart chamber.
(was developed in Department of Nuclear Medicine in Ostrava)

To put this method into practice, we do not need to perform any series of phantom calibration measurements, just a simple method (two point sources at a distance of 10 cm) to measure the scale of the camera display and save it in the CAMERA.PAR file.
Note: With good zoom coding (1, 1.2, 2 or 50%, 75%, 100%), it is sufficient to set the display scale in the full-size field of view (zoom 1 or 100%) and the actual display scale for other zoom values. the program calculates itself.

2. Geometric ellipsoidal method
This method calculates the end-diastolic volume of the heart chamber based on a geometric analysis of the size of the ROI of the heart chamber. The shape of the chamber is approximated by a rotating ellipsoid, the major axis of which is given by the two farthest points in the ROI and the minor axis being counted as the axis of an ellipse which has a major axis determined as before and whose area is equal to the ROI of the chamber. To put this method into practice, it is first necessary to perform a series of calibration measurements using a phantom (eg a rubber balloon) for different volumes in the range of about 60 - 300 ml. and the coefficients regres
e vl thus obtained
to add to the CAMERA.PAR file. This method was developed at our workplace at the same time as the first version of the VENTR program on the CLINCOM device (around 1977) and is still used occasionally. When carefully calibrated and performed, it provides good results, but is very sensitive to the accuracy of the chamber ROI marking. We therefore recommend using a rather proportional geometric method.

Geometric method for determination of absolute ventricular volume by radionuclide ventriculography

3. Simple sample method
The sample method determines the end-diastolic volume of the heart chamber based on the ratio of the integral of the number of pulses in the ROI of the heart chamber and the number of pulses measured in the blood sample taken. We take a blood sample of a defined volume during a radioventriculographic study, then place it in the field of view of the camera and measure the number of pulses for a defined time (or take a static scintigraphic image and determine the integral number of accumulated pulses). When evaluating with the VENTR program, after selecting a simple sample method, we first enter the measured activity of 1 ml. blood sample in pulse count./sec. The following is a question of how to make a correction for the absorption of gamma radiation from the left ventricle. The absorption depth of the chamber can be specified explicitly, or the program will determine it approximately from the height and weight of the patient. The sampling method should theoretically provide more accurate values of ventricular volume, but is more difficult to practice (blood sampling and measurement) and uncertainty in the absorption of radiation from the ventricle may also introduce certain errors.

Sample method for determination of absolute ventricular volume by radionuclide ventriculography

4. Sample method combined with the integral of the first flow
This relatively complicated method is based on the relationship between the applied bolus activity and the integral of the first bolus flow through the heart chamber in bolus angiocardiography, which precedes just before equilibrium ventriculography. Before application, the bolus activity is measured in a suitable arrangement either under the camera (in a suitable absorbing case) or another detector or activity meter (we must have a conversion factor between this measurement and direct measurement on the camera). Then we apply the bolus and store the angiocardiographic study, followed by a radioventriculographic study. When evaluating a first-pass cardiographic study, we determine the integral of the first bolus flow through the left ventricle (using exponential extrapolation). When calculating the EDV method 4 in the VENTR program, we enter the value of the applied bolus activity in imp./sec., integral of the first flow in imp. and heart rate at first transit. This method has the advantage ofthat it does not need to determine the depth of chamber placement, but is generally methodologically very demanding and is applicable only in combined RKG + ventriculographic studies (details are given in the book Ullmann V., Kuba J., Kuchař O., Mrhač L., Dudzik J. : Computer processing of dynamic scintigraphic studies, final report of research task No. 30-02-03 of the Ministry of Health, 1985.

5. Do not calculate ED volume
The calculation of end-diastolic volume is
optional . If we select the “Do not count” option, the EDV does not count, which is expressed on the display and in the final protocol by the statement " End-diastolic volume not counted ". We choose this option when, for example, due to inappropriate projection, the heart chamber overlaps with other structures or part of the chamber is missing and the EDV value would not be correct.

For each method, the calculated end-diastole value appears on the display. volume with the question whether we agree with the calculation. If we answer in the negative, the table of methods for determining EDV will reappear and we can repeat the calculation with the same or another method (with possibly modified parameters).

Mathematical analysis of the volume curve; ejection fraction

To calculate the quantitative parameters of central hemodynamics, curves of the time course of radioactivity in the marked ROIs (ie the left ventricle and background, or even the right ventricle and its background) are created . The following is a mathematical processing of these curves, in which some corrections of the curve must first be made :

Event correction. distortion (artificial reduction) of the end section. This distortion is caused by fluctuations in the period of the cardiac cycle during the study, which leads to a slightly smaller number of cycles contributing to the last few images than to the other images. The program corrects this distortion by multiplying the respective endpoints of the curve by a coefficient created as the ratio of the number of cycles accumulated in the initial section of the curve and the number of cycles accumulated for the respective endpoint section.

The time course curve of radioactivity in the heart chamber shows statistical fluctuations, which can have an adverse effect at some stages of mathematical processing (especially in derivation). Therefore, it is useful to suppress these interfering statistical fluctuations by appropriate smoothing of the m curve. However, conventional smoothing methods (eg multiple 3-point smoothing), if they are to be sufficiently effective, also distort the dynamics of the curve - it can lead to a reduction in the pulsation amplitude and thus to a reduction in the value of the ejection fraction. In the VENTR program, a special algorithm is used to smooth the curve, which smooths the curve strong enough without causing its distortion - the values of dynamic parameters such as ejection fraction are not affected.

The time course curve of radioactivity in the ROI of the heart chamber captures the time dependence of the instantaneous radioactivity not only in the chamber, but is superimposed with the radioactivity of the tissues above and below the chamber during planar camera projection. In order to obtain a "pure" curve of the course of radioactivity in the heart chamber itself, it is necessary to correct it against the tissue background . The exact correction of this kind is fromproblematic because the tissue background is inhomogeneous and it is not known what part of the background is caused by Compton-scattered radiation from the ventricle, blood supply to the heart wall and blood supply to other tissues - the proportion of constant and variable background components is not known. In the VENTR program, a curve from the above-described "moon-shaped" ROI between the contours of the ventricle in the ED and ES is used for background correction, from which the value at the points corresponding to the end-systole is taken for background subtraction.

After correction for the tissue background, an undistorted curve of the time dependence of the instantaneous radioactivity (proportional to the instantaneous amount of blood) in the heart chamber during the cycle *) is formed, from which we can objectively calculate dynamic parameters. To calculate the quantitative parameters, the points of end-diastole ED, end-systole ES, beginning and end of ejection and filling are automatically defined on the curve (with the possibility of manual modification).
*) In general, the dynamics of the curve could be distorted by two other factors - dead time and self-absorption of gamma radiation in a pulsating chamber. Fortunately, none of these adverse factors apply to equilibrium ventriculography for the following reasons.
The loss of pulses through the dead time is determined by the total frequency of pulses in the entire field of view of the camera. In ventriculography, both the ventricles and the atria, which pulse in the opposite direction, are usually in the visual field - in systole, the decrease in pulse frequency from the ventricle is roughly compensated by the increase in atrial radioactivity, and in diastole the growth of ventricular activity is compensated by the pulse rate. Total radioactivity and therefore the pulse frequency in the whole visual field
j e are roughly constant, so is the constant or. percentage of loss of impulses due to dead time - therefore it will have no effect on the quantification of dynamics.
As for self-absorption, if the heart chamber pulsed under the camera in free space
without the surrounding environment, then when the chamber is filled, the self-absorption of gamma radiation in the blood volume of the chamber itself would be somewhat greater than when the chamber is emptied. In reality, however, the heart chamber pulsates in an environment with an absorption coefficient approximately the same as that of the blood in the chamber, which changes the situation. The relevant calculation and phantom measurements (both were performed as part of our research task ..........) show that if a chamber of approximately spherical or elliptical shape is isotropically pulsating in an environment with the same absorption coefficient, not only is there no increase in loss by absorbing gamma radiation in the filled chamber, but on the contrary, there is a slightly lower total absorption in the filled chamber than in the empty chamber. In practice, there is no appreciable effect on the pulsation dynamics by the absorption of gamma radiation, so that the corresponding curve V (t) can be considered as the actual volume curve of the ventricle.

The EDV value, determined above, makes it possible to create a volumetric curve V (t) from the normal curve A (t) of the time course of radioactivity in the heart chamber , each point of which indicates the instantaneous volume V of the heart chamber at the respective time t of the heart cycle. Important parameters of central hemodynamics are calculated from this curve (Fig. 3.1.7) :

Fig.3.1.7. Important parameters of central hemodynamics calculated by volume curve analysis.

Some time intervals can also be calculated from the volume scream, but a more accurate determination of time intervals is performed on a derived (velocity) curve - see below. The calculated volume parameters of the left ventricle are displayed together with the volume curve (time is marked on the horizontal axis, volume in [ ml. ] Is calibrated on the vertical axis ) and we can print them out for documentation - Fig. 3.1.7. If the ROI of the right ventricle has also been marked, this quantification will also be performed for the right ventricle.

Quantification of speed parameters

The time derivative of the volume curve V (t) gives rise to the velocity curve R (t) - Fig.3.1.8

R (t) = d V(t) / dt

which at each of its points indicates the instantaneous rate of ejection or ventricular filling during the cardiac cycle. The smoothing of the velocity curve so that it does not affect the dynamics is performed analogously to the volume curve. The points of the maximum ejection and filling speed, the beginning and end of the ejection, the beginning and the end of the fast filling phase are defined by the program on the speed curve. The relevant speed parameters are calculated :

Further, from the velocity curve, a number of time intervals characterizing the time schedule of significant phases of the heart revolution are determined :

Fig.3.1.8. Velocity and time parameters of ejection and filling calculated by velocity curve analysis

The displayed velocity curve (time is on the horizontal axis, ejection and filling speed in [ ml./sec. ] Is calibrated on the vertical axis ) together with a number of speed and time parameters can be printed out again for documentation (Fig.3.1.8). If the ROI of the right ventricle has also been marked, the described calculations will also be performed for the right ventricle.

The resulting protocol

The following is a summary of the most important evaluation results: images of end-diastole and end-systole, contours of the heart wall, paradoxical and pulse image, Fourier images of phase and amplitude including histogram of phase distribution, volume and velocity curve, below them a clear table of the most important calculated quantitative parameters left and right ventricle. At the bottom of the screen, a verbal assessment is offered for editing. These verbal evaluation texts (including the conclusion) can be freely modified or supplemented as required. Among other things, we can reconcile verbal evaluation with the results of quantitative parameters. After the final value, I will appreciate the "Signature:" section, where we will insert the name of the evaluating doctor.

Department of Nuclear Medicine,  University Hospital  Ostrava                   Date: ..............
Name of patient: ....................... .                   Birth certificate number: .........................
Mathematical analysis and complex evaluation of radionuclide ventriculography
Evaluation:
We do not observe regional disorders of heart wall motility on phase scintigraphic images of the cardiac cycle, nor on Fourier images of phase and amplitude.

Conclusion:
Visual evaluation of images of individual phases of the cardiac cycle and quantitative analysis of cardiac dynamics indicate good global and local contractility of the walls of the left ventricle.
                                                                                       
Signature: MUDr. Jozef Kubinyi

Fig.3.1.9. The resulting protocol of complex evaluation of radionuclide ventriculography - an example of a normal finding.

The comprehensive evaluation of ventriculography ends with the printing of a clear report containing administrative data (name of workplace, name and birth number of the patient, name of examination, date, etc.), prominent images, curves, calculated quantitative parameters, verbal assessment of heart wall contractility and final evaluation of the whole examination with by the doctor's signature - fig.3.1.9 (normal) and fig.3.1.10 (markedly pathological case ).

Department of Nuclear Medicine,  University Hospital  Ostrava                 Date: ..............
Name of patient: ....................... .                   Birth certificate number: .........................
Mathematical analysis and complex evaluation of radionuclide ventriculography
Evaluation:
The following regional disorders of left ventricular wall motility are observed in phase scintigraphic images of the cardiac cycle and in Fourier images of phase and amplitude:

   segment hypokinesia:  Virtually all except posterolat.
         
segment akinesia:
   segment asynchrony:
    segment dyskinesia:  A p i c a l - very extensive !       

Conclusion:
Visual evaluation of images of individual phases of the cardiac cycle and quantitative analysis of cardiac dynamics indicate a severe disorder of left ventricular wall contractility with extensive hemodynamically significant apical dyskinesia.
Extremely reduced ejection fraction dilated LV
                                                                                       
Signature: MUDr. Jozef Kubinyi

Fig.3.1.10. The resulting protocol for a comprehensive evaluation of radionuclide ventriculography
- an example of an extremely pathological finding.

Separate Fourier analysis
Fourier phase analysis can be performed not only as part of a comprehensive evaluation within the VENTR program, but also separately. This separate design is advantageous, for example, in ventriculography in LL projection, where the quantification of ejection fraction and other dynamic parameters would not be correct due to the mutual overlap of cardiac structures, but we are interested in the local contractility of the respective heart wall segments. The FOURI program, which start,is used for this purposewe take from the menu ClinProg. The only difference from phase analysis in the VENTR program is that the FOURI program first generates a time course curve of radioactivity in the entire field of view in order to determine the endpoint of the Fourier analysis. Images of phase and amplitude and histograms of phase distribution together with a table of quantitative data of regional amplitude-phase analysis of evaluated areas can be printed on a printer. Fourier phase analysis program Four projections at j different s than LAO recommended to perform before a comprehensive evaluation program LAO projection ventry to any knowledge of motility disturbances respective segments of the heart wall (invisible in LAO projection) may include in the final word evaluation program VENTR.

Program structure
The VENTR program consists of the following parts (capable of independent function):

VENTR1 - precise determination of end-diastole and end-systole phases, summation of their images, construction of paradoxical image, heart rate image and combined paradoxical-heart rate image
VENTR2 - automatic delimitation of heart chamber and background
ROI VENTR3 - photography of ED, ES, paradoxical and heart rate images, wall contours in ED and ES, Fourier analysis
VENTR4 - verbal evaluation of heart wall contractility, creation of curves
VENTR5 - mathematical processing of curves and calculation of quantitative parameters of central hemodynamics, editing and printing of the final protocol

At the same time, this structure shows how to proceed when the program is interrupted or restarted in order to repeat some part of the calculation. If we want to repeat the mathematical processing of curves or modify the wording of the word evaluation and print a new protocol, just run the last part of the program - VENTR5. If we want to perform Fourier phase analysis (whether to repeat it or we have not performed it before), we run VENTR3, or we will use a separate FOURI program (as mentioned above).

Cast of SAVE AREA after the end of the program:
SA 3 - heart rate image
SA 4 - paradoxical image
SA 5 - end-diastole image including ROI
SA 6 - end-systole image including ROI
SA 7 - curves from ventricular and background areas
SA 8 - combined paradoxical +
heart rate image SA 9 - heart wall contours in ED and ES
SA10 - Fourier image of amplitude
SA11 - Fourier image of phase

2. Incorporation of complex programs   3.2. Bolus radiocardiography

Nuclear physics and physics of ionizing radiation   |  OSTNUCLINE
Nuclear and radiation physics Radiation detection and spectrometry Radiation applications
With cintigraphy Computer evaluation of scintigraphy Radiation protection
Gravity, black holes and spacetime physics ;Anthropic principle or cosmic God
AstroNuclPhysics ® Nuclear Physics - Astrophysics - Cosmology - Philosophy

Vojtech Ullmann