A novel MRI-based three-dimensional model of stomach volume, surface area and geometry in response to gastric filling and emptying ================================================================================================================================== * D Bertoli * EB Mark * D Liao * C Brock * JB Frøkjær * AM Drewes ## ABSTRACT **Background** Gastric motility and accommodation have a critical role in maintaining normal gastrointestinal homeostasis. Different modalities can be adopted to quantify those processes, i.e., scintigraphy to measure emptying time and intragastric Barostat for accommodation assessment. However, magnetic resonance imaging (MRI) can assess the same parameters non-invasively without ionizing radiation. Our study aimed to develop a detailed three-dimensional (3D) MRI model of the stomach to describe gastric volumes, surface areas, wall tension distribution, and inter-observer agreement. **Methods** Twelve healthy volunteers underwent an MRI protocol of six axial T2-weighted acquisitions. Each dataset was used to construct a 3D model of the stomach: Firstly, the volumes of the whole stomach, gastric liquid, and air were segmented. After landmark placing, a raw 3D model was generated from segmentation data. Subsequently, irregularities were removed, and the model was divided into compartments. Finally, surface area and 3D geometry parameters (inverse curvatures) were extracted. The inverse curvatures were used to as a proxy for wall tension distribution without measuring the intragastric pressure. **Key Results** The model was able to describe changes in volume and surface geometry for each compartment with a distinct pattern in response to filling and emptying. The surface tension was distributed non-homogeneously between compartments and showed dynamical changes at various time points. **Conclusion & Inferences** The presented model offers a detailed tool for evaluating gastric volumes, surface geometry, and wall tension in response to filling and emptying and will provide insights into gastric emptying and accommodation in diseases such as diabetic gastroparesis. **KEY POINTS** * - While MRI and ultrasound are getting progressively accepted as methods for evaluating gastric emptying and accommodation, they still do not provide essential insights into those processes. * - This study presents a three-dimensional stomach model able to report volume and surface data, and to describe the distribution of the gastric wall tension of gastric compartments by applying the Young-Laplace law. * - We observed that volumes and surface geometry showed distinct emptying patterns in each compartment and that the wall geometry distributed non-homogeneously in the stomach, showing different dynamical changes during the emptying phase. Our observations indicate the fundus as essential in the first phase of digestion, confirming its role as reservoir. * - The non-invasive model has the potential to give detailed information about the gastric volumes and surface geometry in response to filling and emptying, with the potential to understand the pathophysiology and improve treatment in patients with gastroparesis. ![Figure1](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/06/23/2022.06.21.22276694/F1.medium.gif) [Figure1](http://medrxiv.org/content/early/2022/06/23/2022.06.21.22276694/F1) Keywords * Abdomen * Algorithms * Gastric emptying * Magnetic resonance imaging * Reproducibility ## INTRODUCTION In modern practice, the interest in the function of the stomach has led to an increasing focus on the assessment of gastric emptying, motility, and accommodation, as these processes are closely associated with each other. Various examinations are commonly accepted such as: gamma scintigraphy,1,2 single-photon emission computed tomography,3,4 breath testing,5,6 intragastric Barostat,7 ultrasound,8,9 smartpill,10 and magnetic resonance imaging (MRI).11,12 Some of these methods are limited by invasiveness, the use of radioactive isotopes or by the inability to evaluate gastric emptying or anatomy. In contrast, others are time-consuming and offer only indirect measurements. Gamma scintigraphy is presently accepted as the gold standard in assessment of gastric emptying as it provides a reliable measure that does not suffer from confounding gastric secretions.13 Unfortunately, it does not offer a precise evaluation of gastric volumes, the contents’ intragastric distribution and wall geometry, resulting in a significant lack of insight into the mechanisms and pathophysiology behind diseases such as gastroparesis. Furthermore, its use is limited due to radiation exposure. The current standard for assessment of accommodation is the intragastric Barostat, an invasive technique not devoid of disadvantages and technical limitations.7,14 Because of these limitations, ultrasound and MRI are becoming progressively more attractive due to their ability to morphologically evaluate gastric emptying and motility, including intra-gastric volume distribution and geometry,15–17 without exposing the patient to unnecessary ionizing radiations or invasive procedures.18 However, ultrasound-based evaluation methods are prone to be observer-dependent and unable to visualize gastric air.19 MRIs, on the other hand, can visualize the entire stomach in one scan session providing three-dimensional (3D) morphological and functional information through repeated scans. While MRI examinations of the stomach and of the gastric function have been of interest in recent literature, the existing methods still have methodological limitations. For example, the intragastric pressure cannot be directly quantified and, therefore, the absolute wall tension value cannot be measured. Tension data are of utmost importance in the clinical application of these methods, as the gastric wall tension determines the perception of gastric distention.20–22 We hypothesized that an MRI 3D model would be able to provide essential insights into gastric geometry, motility, and accommodation. Furthermore, we hypothesized that through the application of the Young-Laplace law on the collected geometric data, the model would be able to describe the distribution of the wall tension across different compartments. Therefore, the primary aim of this study was to present a framework for a detailed analysis of gastric volumes, surface areas, and wall tension distribution in response to gastric filling and emptying using T2-weighted MR images before and after a liquid meal. The secondary aim was to apply the method to 12 healthy subjects, exploring the details of gastric filling and emptying processes. Finally, we aimed to measure the inter-observer agreement between two raters to further validate the observations on which our model was built. ## MATERIALS AND METHODS ### Study subject selection Data were obtained from 12 healthy subjects, normal-weighted, without prior history of gastric disorders or other diseases affecting the gastrointestinal (GI) function. None of the subjects was treated with any medication that could affect the GI system. Informed consent was obtained per national and local institutions’ ethical standards. The study was approved by the North Denmark Region Committee on Health Research Ethics (N-20090008). ### Study design Following a minimum of 6 hours of fasting (solids and liquids), a baseline scan (at t:-30 min) was obtained in a transversal plan with the subjects in supine position. Subsequently, subjects were asked to step outside the MRI scanner for 10 min and ingest a liquid meal of lightly heated 250 kcal tomato soup (500 ml “Karolines Køkken” (Arla Foods, Central Denmark Region, Denmark). Nutritional information per 100 g: fat 3g, carbohydrates 4.1g, proteins 1.1g. After ingesting the meal, the subjects were scanned at five additional time points (0, 15, 30, 45, and 60 minutes) with the same sequence as the baseline scan. ### MRI protocol MRI scans were performed using a 1.5T General Electric model Discovery MR450 (General Electric Medical Systems, Milwaukee, WI, USA). Axial steady-state gradient echo (FIESTA) T2-weighted image series covering the entire stomach was obtained with an echo time (TE) = 1.5 ms and a repetition time (TR) = 3.5 ms, in-plane resolution = 0.7422 x 0.7422 mm, flip angle = 45°, slice thickness = 5 mm, 30-35 slices, and no image gap or overlap. Fat signal was not suppressed. Images were obtained in approximately 20 seconds under a single breath-hold to minimize respiratory artifacts. As shown in previous similar studies,23,24 this sequence yielded a high signal of water and high contrast between gastric fluid content, comprehending meal and secretions, and gastric air as well as between the total gastric content and surrounding body tissues. No contrast-enhancing agents were deemed necessary to increase the accuracy of the measurement. An example of raw MRI images is shown in figure 1(A). ![Figure 1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/06/23/2022.06.21.22276694/F2.medium.gif) [Figure 1:](http://medrxiv.org/content/early/2022/06/23/2022.06.21.22276694/F2) Figure 1: *Three-dimensional model development*. Firstly (A) axial T2-weighted images of the stomach were obtained. Then (B), the volumes of liquid gastric content (blue), air content (white), and total gastric content (red) were segmented. Landmarks at the central part of the cardia, angular incisure, and the topmost part of the stomach were placed in this phase (C). Secondly, a raw 3-dimensional model was generated from the segmented total gastric volume images (D). Irregularities were then removed, and the model was divided into compartments. Finally, inverse curvature data could be analyzed on the surface (E). ### Model preparation #### Manual segmentation of gastric fluid content, gastric air, and total gastric volume The segmentation of high-signal gastric fluid content, low-signal gastric air, and total gastric volume was manually performed using a segmentation platform implemented in MATLAB (v.R2018b, MathWorks, Natick, MA, USA). An example of the segmentation steps is shown in figure 1(B). The total gastric volume was defined by the outer border of the stomach, including the stomach wall. The gastric fluid content volume was defined as the high-signal volume within the total gastric volume (including areas with low signal caused by residues in the stomach or signal artifacts), while the gastric air was defined as the low-signal volume within the total gastric volume. Segmentations were performed by two observers (EBM and CS). Consensus was reached before further analysis. Firstly, segmentation of gastric fluid content and gastric air was performed. Secondly, an approximation of the total gastric volume boundary was initially drawn automatically using the previously segmented volumes to avoid redundant work. Then the boundaries were manually edited to include the gastric wall and gastric areas not included in the other volumes. Segmentations of all three compartments took approximately 2 hours per subject’s dataset, comprehending all 6 scans. #### Compartments and geometry of the stomach The stomach volume and surface models were divided into antrum, corpus, and fundus. The antrum was obtained by dividing the stomach with a plane (Antrum-Corpus plane, AC-plane) passing the angular incisure and perpendicular to the less curvature in the coronal plane, see figure 1(C). Fundus and corpus were divided by a plane (Corpus-Fundus plane, CF-plane) passing the central point of the cardia and perpendicular to the line between the topmost point of the stomach and the middle point of the AC-plane line in the coronal plane, see figure 1(C).25,26 The angular incisures used to divide the stomach into antrum, corpus and fundus were defined and independently confirmed by two observers (LD and EBM), and if there was any disagreement, consensus was reached before further analysis. #### Three-dimensional volume and surface modelling of the stomach 3D raw stomach volume models were generated from the segmented total gastric volume data, see figure 1(D). The volume of the total stomach model was calculated in mL by multiplying the number of image pixels, pixel size, and slice thickness. Based on the volume of the 3D stomach model, a first raw surface was computed using the isosurface function in MATLAB with an isovalue of 0. The point cloud obtained from the isosurface was down-sampled using a box grid filter in a box size of 5 mm and then presented as the surface of the stomach, represented with triangular facets, see figure 1(D). In order to remove the high curvature variations while avoiding shrinkage of the original surface, a modified non-shrinking Gaussian smoothing method was used to remove the irregularities of the reconstructed stomach surfaces due to the discretization of the images.27 The surface area of the entire stomach was then calculated in cm2 by adding up the areas of all triangular facets.28 The same method was previously applied, e.g., in a work of Liao et al.29 At this step of our pipeline, geometric data were extracted from the model, and inverse curvature data was calculated. For details of surface computations, see Appendix 1: *surface smoothing*. The processing of geometric data will be discussed in the following sections. ### Data processing #### Gastric emptying half-times The change of gastric liquid, air, and total gastric volume during the emptying phase was curve-fitted to a linear-exponential model (LinExp),30 allowing to quantify the gastric emptying half-times as: ![Formula][1] This model can handle an initial volume increase (due to gastric secretions) using the coefficient kappa (κ) and the subsequent volume decrease described by the *tempt* coefficient.24 The *V* coefficient describes the start volume (t:0 min). This calculation was performed with the online tool’ apps.menne-biomed.de/gastempt/’. This tool has been previously utilized in other studies.31 In this work, the fitting method *nlme population fit* was utilized. #### Wall tension distribution In this work, the gastric wall tension was not directly measured or estimated but through the following principles and equations it was possible to describe how the tension was distributed in different compartments (called *Wall tension distribution* in this work). In medicine, the Young-Laplace equation (2) has been used to describe multiple phenomena like the formation of diverticula or the expansion of abdominal aorta aneurysms32: ![Formula][2] where ΔP is defined as Δ*P* = *Pi* − *P*o, where *P*o is the intra-abdominal pressure and *Pi* is the intragastric pressure; assuming isotropic tensions, *γ* is the tension on the gastric wall, and r1 and r2 are the principal radii of curvature. In differential geometry, the principal curvatures at a given point on a surface are the maximum and minimum values of all the possible curvatures obtained through the section of the surface by planes tangent to the surface containing the normal vector ![Graphic][3]. These two curvatures are denoted k1 and k2 and can be expressed as: ![Formula][4] Integrating (3) in (2): ![Formula][5] Resolving (4) for the ratio *γ* / Δ*P*, (5) is obtained: ![Formula][6] Introducing in (5) the mean curvature, an extrinsic measure of curvature defined in fluid mechanics as: ![Formula][7] We obtain that: ![Formula][8] Furthermore, assuming Po and Pi (and therefore Δ*P*) to be constant (*c*) within each timepoint, the tension would only be proportional to the inverse mean curvature, allowing us to describe the wall tension distribution in different gastric compartments within the same timepoint, given only geometric data: ![Formula][9] and therefore: ![Formula][10] To report this distribution, inverse curvature data were normalized to antral values at the same timepoint, as we expected the volumes to be accommodated mainly in the proximal stomach.26 As the constant *c* is of unknown value, this method does not allow us to estimate absolute wall tension values but allows us to evaluate its distribution on the gastric surface. For a visual representation of the forces involved, see figure 2. *k*1 and *k*2 (expressed as mm−1) were calculated from the 3D surface model of the stomach using the algorithm presented by Hartmann et al.28 See Appendix 1: *principal curvatures computation*. These analyses were at last applied to the previously generated 3D model. See figure 3 and figure 1(E). ![Figure 2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/06/23/2022.06.21.22276694/F3.medium.gif) [Figure 2:](http://medrxiv.org/content/early/2022/06/23/2022.06.21.22276694/F3) Figure 2: *Schematic representation of the physical forces involved in the accommodation process*. Po: abdominal pressure, P*i*: intragastric pressure, k1 and k2 are the principal curvatures with r1 and r2 their relative radii, N is the normal vector perpendicular to the surface, *γ* is the tension on the gastric wall. ![Figure 3:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/06/23/2022.06.21.22276694/F4.medium.gif) [Figure 3:](http://medrxiv.org/content/early/2022/06/23/2022.06.21.22276694/F4) Figure 3: *Tension distribution*. The image shows an example of an empty stomach at baseline, characterized by lower wall tensions (blue color due to low inverse mean curvature values), compared to an image of a full stomach at t:0 min, characterized by higher tensions, particularly in the fundus (yellow/red colors due to high inverse mean curvature values). Another important assumption in our model is that the stomach was represented as a thin-walled membrane. This assumption is based on the premise that to properly apply the Laplace’s law on a non-infinitesimal-walled structure, the following requirement should be met: ![Formula][11] where *t* is the wall thickness, *ri* the internal radius, and *P*o the intra-abdominal pressure. Studies reported mean intra-abdominal pressures *P*o of 1.7 (1.2) cmH20 and mean intra-gastric pressures *Pi* of 2.9 (1.7) cmH2O, resulting in ![Graphic][12] in the pre-swallow resting state.33,34 Resulting in: ![Formula][13] This requirement was met in our dataset. Wall thickness data were not reported as outside the scope of this work. Applying these equations to the geometric data collected allow to describe the distribution of the wall tension in different compartments without measuring or estimating the intragastric pressure. ### Data normalization Volume and surface area data were normalized to their total values to better convey the compartment dynamics during the emptying phase. The inverse curvature data were normalized to antral values at the same timepoint (see *Wall tension distribution* above). ### Inter-observer agreement As the initial part of this framework is heavily observer-dependent, inter-observer differences in the total gastric volume were explored through Dice similarity coefficients. For validation, intraclass correlation coefficient (ICC) assessments between the two raters (EBM and CS) were performed on the extracted geometric data (inverse curvatures) on which our model was built. Data were obtained by two different raters (EBM and CS). No limit of agreement was established a priori. ### Statistical analysis Data were analyzed in R (R Core Team, 2021, R Foundation for Statistical Computing, Vienna, Austria) with *rstatix* package (R package version 0.7.0), and figures were produced using the package ggplot2.35 After identifying and removing extreme outliers (defined as values above the third quartile + 3 interquartile ranges or below the first quartile – 3 interquartile ranges), volume, surface area, and inverse curvature data were inspected for normality through Q-Q plots and Shapiro-Wilk tests. Correlation tests were then performed with two-way repeated-measure analysis of variances (ANOVAs, *time* and *compartment* were used as independent variables, dependent variable: volume, surface areas, or inverse curvature, both non-normalized and normalized), tested for sphericity with Mauchly’s test, and corrected with Greenhouse-Geisser and Huynh-Feldt corrections. In the case of positive correlation, posthoc analyses were performed with pairwise t-tests adjusted with sequential Bonferroni tests (Bonferroni-Holm). A *p*-value < 0.05 was considered statistically significant. Data were reported as mean (SD). ICC estimates and their 95% confidence intervals were calculated using R packages *blandr* and *psych*.36,37 The model ICC1 (single-measure, absolute-agreement, two-way mixed effect) was utilized. Data was not transformed for analysis and were visually reported through a Bland-Altman plot. ## RESULTS All the recruited subjects fully complied with the study protocol, and a complete dataset was obtained from all twelve subjects (7 males), mean age 27.8 years (range 23-34 years). All the subjects well tolerated the meal. ### Gastric volumes Volumes of gastric fluid content, gastric air, and the total gastric volume can be seen in figure 4(A). ![Figure 4:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/06/23/2022.06.21.22276694/F5.medium.gif) [Figure 4:](http://medrxiv.org/content/early/2022/06/23/2022.06.21.22276694/F5) Figure 4: *Gastric volume data*. A: total gastric volume (red), total liquid volume (blue), and total gas volume (green) are shown. B: Normalized fundus volume (red), normalized corpus volume (green), and normalized antrum volume (blue) are shown. All data are reported as mean ± 1 standard error (error bars). The total gastric volume increased from baseline values 140.3±31.9ml to 669.1±41.4ml at t:0 min (*p* < .001) and decreased to 428.7±96.8ml at t:60 min (*p* < .001), remaining however over baseline values (*p* < .001). The initial increase in volume was mostly due to the ingested meal: the liquid volume increased from 38.7±23.2ml at baseline to 515.8±29.9ml at t:0 min (*p* < .001). The liquid volume showed a similar decreasing pattern with a minimal volume of 293.6±73.9ml at t:60 min, never reaching baseline values (*p* < .001). The gastric gas volume, after a limited initial increase from baseline values of 26.8±13.7ml to 108.7±54.9ml at t:0 min (*p* < .001), remained constant and over baseline values throughout the examination (*p* < .001), with a final volume of 101.7±62.9ml at t:60 min. All compartments’ total volumes (comprehending gastric wall, liquid, and air content) showed a similar dynamic with an initial increase (all *p* <.001) from baseline values to t:0 min and a successive decrease to t:60 min (all *p* <.05). No compartment volume returned to baseline values at t:60 min (all *p* <.05). Data showed a difference at baseline between the volumes of each compartment with a fundus volume of 13.4±6.0ml, corpus volume of 80.4±26.7ml, and antrum volume of 46.5±18.2ml (all *p* <.001). This observation was confirmed by a two-way repeated ANOVA that showed an influence of both time and compartment localization on the compartment volumes (all *p* < .001). Normalized data showed that fundus volumes increased from 9.6±3.6% at baseline to 19.3±7.5% at t:0 min and never returned to baseline values throughout the examination, with a volume of 14.6±5.3% at t:60 min. No change from baseline values of corpus (56.4±10.7%) or antrum (64.0±12.0%) was observed (all *p* > .05). The different dynamics can be observed in figure 4(B). The complete dataset can be found in Appendix 2, Table A1. ### Assessment of emptying half-times The emptying half-times of the liquid volume showed a mean value of 69.3 (SD 15.2) minutes. The emptying half-times of the total gastric volume showed a mean value of 86.0 (SD 18.0) minutes. ### Gastric surfaces The total gastric surface area increased from 220.4±27.2cm2 at baseline to 536.4±25.0cm2 at t:0 min (*p* < .001) and decreased to 413.4±54.2cm2 at t:60 min. See figure 5(A). ![Figure 5:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/06/23/2022.06.21.22276694/F6.medium.gif) [Figure 5:](http://medrxiv.org/content/early/2022/06/23/2022.06.21.22276694/F6) Figure 5: *Gastric surface area data*. A: Total gastric surface area is shown. B: Normalized fundus surface area (red), normalized corpus surface area (green), and normalized antrum surface area (blue) are shown. All data are reported as mean ± 1 standard error (error bars). All compartmental surface areas showed a similar dynamic with a starting increase from baseline to t:0 min (all *p* < .001) and a successive decrease toward t:60 min (all *p* < .05). Fundus surface area was smaller than the other compartmental areas (all *p* < .001) with a baseline value of 24.5±5.7cm2 for fundus, 111.8±22.6cm2 for corpus, and 84.2±30.8cm2 for antrum. This observation was confirmed by a two-way repeated ANOVA that showed an influence of both time and compartment localization on the compartment surface (all *p* < .001). Normalized data showed that the fundic surface area increased from 11.2±2.8% at baseline to 20.7±4.1% at t:15 min (*p* = .001), never returning to baseline values throughout the examination. This increment was mirrored by an initial decrease of corpus surface area from 50.9±9.6% at baseline to 41.8±8.9% at t:15 min (*p* < .001). Corpus surface area never returned to baseline throughout the study. No change was observed in antrum surface area, which remained around baseline (37.9±10.7%) during the examination. The different surface area dynamics can be seen in figure 5(B). The complete dataset can be found in Appendix 2, Table A2. ### Gastric wall tension distribution The wall tension distribution was different between compartments (*p* < .001), with mean inverse curvature values in the fundus being higher than those in the corpus (*p* < .001) and corpus values being higher than those in the antrum (*p* < .001) at every timepoint. The dynamical changes in the tension distribution can be seen in figure 6(A) for non-normalized data and figure 6(B) for normalized data. Data showed that the corpus tension was 26±26% higher than the antral tension (*p* = .017). At t:0 min, the tension in the fundus increased to values 34±34% higher than those in the antrum (*p* = .018), while corpus tension values still were 21±24% higher than antrum ones (*p* = .026). Fundus tension continued to increase to t:15 min to 51±32% of the antral tension (*p* < .001). From t:15 min, the tension distributed toward the antrum with values in fundus and corpus steadily diminishing throughout the examination to t:60min. ![Figure 6:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/06/23/2022.06.21.22276694/F7.medium.gif) [Figure 6:](http://medrxiv.org/content/early/2022/06/23/2022.06.21.22276694/F7) Figure 6: *Inverse curvature data*. A: the non-normalized values of fundus (red), corpus (green), and antrum mean inverse curvature value (blue) are shown. B: fundus (red) and corpus (green) inverse curvatures are shown normalized to antral (blue) values at the same timepoint. The complete dataset can be found in Appendix 2, Table A3. ### Inter-observer agreement Data showed a median Dice similarity coefficient for total gastric volume of 0.957 (IQR = 0.021) between the two raters. Dice similarity coefficients for liquid and gas can be seen in figure 7(A). Median values ranged from baseline values: 0.846 (0.047) for total gastric volume, 0.851 (0.101) for liquid volume, and 0.828 (0.020) for gas volume, to 0.965 (0.009) for total gastric volume, 0.975 (0.013) for liquid volume, and 0.902 (0.050) for gas volume at t:0 min. ![Figure 7:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/06/23/2022.06.21.22276694/F8.medium.gif) [Figure 7:](http://medrxiv.org/content/early/2022/06/23/2022.06.21.22276694/F8) Figure 7: *Dice similarity coefficient and intraclass correlation coefficient data*. A: Dice similarity coefficients for the liquid (blue) and gas (green) volumes are shown. B: Bland-Altman plot for the correlation agreement between the two inverse curvature datasets across every timepoint. Blue dots: observations, blue line: bias, red dashed lines: 95% limits of agreement. ICC validation results can be seen in figure 7(B). The estimated bias was −0.07mm−1 (−0.28, 0.13), with lower limit of agreement: −1.78mm−1 (−2.12, −1.43) and upper limit:1.63mm−1 (1.28, 1.98). The correlation coefficient was 0.87 (lower bound 0.79, upper bound 0.91, *p* < .005). ## DISCUSSION The present study presents a 3D MRI-based model for a detailed geometric assessment of gastric compartment volumes, surface areas, and wall tension distribution using MRI. The most important findings study were: 1) The model was able to describe changes in volume, surface areas, and wall tension distribution in fundus, corpus, and antrum; 2) The changes in the observed gastric variables in response to gastric filling and emptying showed a distinct pattern for each compartment; 3) The wall tension distributed non-homogeneously between different compartments and showed different dynamical changes during the emptying phase. Furthermore, during gastric filling, we observed that the most significant changes in volume and surface area occurred in the fundus, underlining its role as a reservoir of undigested food. The geometric data on which our model was built also gave us important insights into the dynamics of the emptying process, as, during the first 15 minutes after meal intake, the normalized fundus tension increased, representing its predominant role in the first phase of the emptying process. ### Gastric volumes and emptying As expected, the volume of the ingested liquids played a major role in the observed total gastric volume increase, while the gas volume was only accessory and generally constant throughout the examination. Furthermore, we observed that the fundus played a major role in accommodating newly ingested meals, as it was the only compartment to increase its normalized volume and surface area values during the first 15 min. The observed fasting volumes were in line with the literature,38,39 as were the dynamics of the emptying process.40 The volume and dynamics of total gas content were also in line with findings observed in the literature.26 With a careful approximation (due to differences in the compartment subdivision of the stomach), it appears that the antral volumes at the start of the filling phase were also in line with the literature.26 MRI gastric emptying studies have reported emptying half-times between 7 and 330 minutes in healthy subjects for liquid and mixed meals.41 Our half-time measurements of 69.3 min for liquid and 86.0 min for the total volume comply with these observations. Direct comparison with previous data is difficult due to differences between study protocols, such as different amounts and kinds of meals,23 MRI sequences,25 or different sample sizes or demographics.30 Moreover, some studies positioned the patients in a right lateral decubitus, which permits an easier outflow through the antrum,42 and it is not always stated which positions were used between each scan.43 However, the emptying phase dynamic of the fundus, corpus, and antrum in the present study showed a good agreement with a recent study done by Banerjee et al.26 ### Gastric surface areas analysis While no standard reference values of surface areas can be found in the literature, the dynamic of our observations is in line with other recent studies.44 In our work, the fundus appears to sustain the most relevant area increase after meal intake, as it was the only compartment to show an increase in its normalized surface area from baseline to t:0 min. Previously, we have done geometric surface analyses on 3D imaging data from the stomach,45 gallbladder,46 and rectum in humans to investigate the geometry and mechanical change of the models during distension.47 In our previous surface modelling analysis, we re-sliced the 3D volume model along the central line of the model to generate the surface, which was very time-consuming.47 The method presented in this study has a superior time efficiency compared to the previous modelling analysis as it implemented a surface smoothing which could be performed straightforwardly from the volume model reconstructed from the MR images. ### Gastric tension distribution In our observations, the gastric tension was distributed unevenly across different compartments, with the highest inverse curvature values in the fundus and the lowest in the antrum. Furthermore, we observed that during the earliest phase of digestion (15 min), the fundus reached tension values up to 51% higher than those observed in the antrum. These observations concords with the reported changes in volume and surface areas elicited in the fundus by the meal, indicating the fundus as the most influential compartment in the first phase of digestion, confirming its role as reservoir. Our data showed wide variations, as observable by relatively high standard deviations. As data were normally distributed and the inter-observer variability was excellent (see *Inter-observer agreement* below), these differences can be imputed to significant inter-subject variability, which is also often observed in other tracts of the GI system.10,48,49 ### Inter-observer agreement Our analysis showed an excellent agreement in total gastric volumes between the two observers, defined by an overall Dice similarity coefficient of 0.938. Our observed coefficient is significantly higher than 0.7, defined in the literature as a good overlap.50 ICC analysis also showed good reliability, defined as an ICC value between 0.75 and 0.9.51 ### Potential clinical implications Implementation of analytical tools like the currently described image acquisition and analysis framework in the clinical practice will require optimization of the analysis time, including optimized segmentation and landmark annotation. A more automatic and standardized method will rapidly assess gastric function, emptying, and other motility indexes to evaluate and diagnose motility disturbances. These disturbances range from delayed gastric emptying and gastroparesis to abnormally rapid transit, usually referred to as “dumping syndrome”.52,53 Furthermore, it will give access to the assessment of the gastric accommodation process, which plays an essential role in functional dyspepsia symptoms or in gastroesophageal reflux.15 These assessments are vital in clinical practice and pharmacological studies.54 Also, more automatic and standardized methods could be helpful in future research studies, giving them access to fast and straightforward 3D gastric models that could lead to a more robust understanding of the motility and the relative interferences of diseases or potential treatment approaches (drugs, neuromodulation, gastric pacing, and others). ### Study limitations In the present study, the subjects were placed in supine position after meal intake in the following image acquisitions. This position displaces the gastric content away from the antrum, potentially slowing gastric emptying.42 A possible improvement could be obtained by oblique positioning in Fowler’s position (semi- or low-) as that would offer a more realistic setting for the gastric environment.55 However, in the majority of similar studies, the study subjects were also placed in supine position as this is easier with standard MRI equipment.41 An unavoidable issue in the segmentation of gastric volumes is the lack of a precise anatomical boundary between the duodenum and pylorus, making it challenging to outline identical segmented compartments of the stomach during all six scans. Custom boundaries were adopted as previously described, but their overall impact on measurements is still unclear.56 Due to the lack of a direct measurement of the intragastric pressure, estimating exact absolute wall tension values was not possible in this study. This issue could be overcome by measuring the intragastric pressure before and after food intake.57 Furthermore, more insights into the stomach’s motility could be assessed through cine-imaging or through thinner slices, which could potentially be implemented in our framework. ### Conclusions As per our aims, the analysis framework presented in this study offered a tool for evaluating gastric volumes, surface areas, and a dynamic estimation of wall tension based on MRI, feasible with only minor preparation. The method was able to provide insights into both gastric accommodation and emptying. Compared to existing gastric emptying studies, our model is improved by estimating the gastric wall tension distribution changes in different stomach compartments and holds promise for future clinical implementations and research. ## DISCLOSURE No competing interest declared. ## AUTHOR CONTRIBUTIONS Davide Bertoli, Esben Bolvig Mark, Jens Brøndum Frøkjær, Donghua Liao, and Asbjørn Mohr Drewes designed the research; Davide Bertoli, Esben Bolvig Mark, and Donghua Liao were involved in data collection; Asbjørn Mohr Drewes, Davide Bertoli, Christina Brock, and Jens Brøndum Frøkjær directed the data interpretation; Davide Bertoli drafted and wrote the manuscript; All the authors revised the final manuscript. ## Supporting information Appendix 1 [[supplements/276694_file02.docx]](pending:yes) Appendix 2 [[supplements/276694_file03.docx]](pending:yes) ## Data Availability All data produced in the present study are available upon reasonable request to the authors ## **ACKNOWLEDGMENTS** We thank Kenneth Krogh Jensen (Aalborg University Hospital) for the assistance in the data collection and the Radiology Research Unit for providing MRI scanner facilities. We thank Christoffer Svinth for the assistance in the data collection. ## Footnotes * **Grant Support**: This study is supported by a Novo Nordisk Foundation grant. Grant No. #0052045. * Received June 21, 2022. * Revision received June 21, 2022. * Accepted June 23, 2022. * © 2022, Posted by Cold Spring Harbor Laboratory The copyright holder for this pre-print is the author. All rights reserved. The material may not be redistributed, re-used or adapted without the author's permission. ## REFERENCES 1. 1.Solnes LB, Sheikhbahaei S, Ziessman HA. Nuclear scintigraphy in practice: Gastrointestinal motility. Am J Roentgenol. 2018;211(2):260–266. 2. 2.Abell TL, Camilleri M, Donohoe K, et al. Consensus recommendations for gastric emptying scintigraphy: A joint report of the American neurogastroenterology and motility society and the society of nuclear medicine. Am J Gastroenterol. 2008;103(3):753–763. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1572-0241.2007.01636.x&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18028513&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F23%2F2022.06.21.22276694.atom) 3. 3.Henrichon S, Seltzer M, Siegel A. Use of SPECT/CT to confirm inconclusive gastric emptying scintigraphy results. Clin Nucl Med. 2015;40(6):e341–e342. 4. 4.Simonian HP, Kantor S, Linda C. K, et al. Simultaneous assessment of gastric accommodation and emptying of solid and liquid meals. Gastroenterology. 2006;124(4):A53. 5. 5.Bluemel S, Menne D, Fried M, Schwizer W, Steingoetter A. On the validity of the 13C-acetate breath test for comparing gastric emptying of different liquid test meals: A validation study using magnetic resonance imaging. Neurogastroenterol Motil. 2015;27(10):1487–1494. 6. 6.Goetze O, Wieczorek J, Mueller T, Przuntek H, Schmidt WE, Woitalla D. Impaired gastric emptying of a solid test meal in patients with Parkinson’s disease using 13 C-sodium octanoate breath test. Neurosci Lett. 2005;375(3):170–173. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neulet.2004.11.007&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15694254&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F23%2F2022.06.21.22276694.atom) 7. 7.Mundt MW, Samsom M. Fundal dysaccommodation in functional dyspepsia: head-to-head comparison between the barostat and three-dimensional ultrasonographic technique. Gut. 2006;55(12):1725–1730. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NjoiZ3V0am5sIjtzOjU6InJlc2lkIjtzOjEwOiI1NS8xMi8xNzI1IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjIvMDYvMjMvMjAyMi4wNi4yMS4yMjI3NjY5NC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 8. 8.Okabe T, Terashima H, Sakamoto A. Determinants of liquid gastric emptying: comparisons between milk and isocalorically adjusted clear fluids. Br J Anaesth. 2015;114(1):77–82. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bja/aeu338&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25260696&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F23%2F2022.06.21.22276694.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000347412000012&link_type=ISI) 9. 9.Perlas A, Mitsakakis N, Liu L, et al. Validation of a mathematical model for ultrasound assessment of gastric volume by gastroscopic examination. Anesth Analg. 2013;116(2):357–363. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1213/ANE.0b013e318274fc19&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23302981&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F23%2F2022.06.21.22276694.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000314078300016&link_type=ISI) 10. 10.Farmer AD, Pedersen AG, Brock B, et al. Type 1 diabetic patients with peripheral neuropathy have pan-enteric prolongation of gastrointestinal transit times and an altered caecal pH profile. Diabetologia. 2017;60(4):709–718. 11. 11.Bilecen D, Scheffler K, Seifritz E, Bongartz G, Steinbrich W. Hydro-MRI for the visualization of gastric wall motility using RARE magnetic resonance imaging sequences. Abdom Imaging. 2000;25(1):30–34. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s002619910005&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=10652917&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F23%2F2022.06.21.22276694.atom) 12. 12.Faas H, Steingoetter A, Feinle C, et al. Effects of meal consistency and ingested fluid volume on the intragastric distribution of a drug model in humans - A magnetic resonance imaging study. Aliment Pharmacol Ther. 2002;16(2):217–224. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11860404&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F23%2F2022.06.21.22276694.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000173864600006&link_type=ISI) 13. 13.Tomita T, Okugawa T, Yamasaki T, et al. Use of scintigraphy to evaluate gastric accommodation and emptying: Comparison with barostat. J Gastroenterol Hepatol. 2013;28(1):106–111. 14. 14.Drewes AM, Gregersen H, Arendt-Nielsen L. Experimental pain in gastroenterology: a reappraisal of human studies. Scand J Gastroenterol. 2003;38(11):1115–1130. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1080/00365520310004399&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=14686714&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F23%2F2022.06.21.22276694.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000186927400003&link_type=ISI) 15. 15.Buisman WJ, van Herwaarden-Lindeboom MYA, Mauritz FA, et al. Validation of a Novel 3-Dimensional Sonographic Method for Assessing Gastric Accommodation in Healthy Adults. J Ultrasound Med. 2016;35(7):1411–1418. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NToidWx0cmEiO3M6NToicmVzaWQiO3M6OToiMzUvNy8xNDExIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjIvMDYvMjMvMjAyMi4wNi4yMS4yMjI3NjY5NC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 16. 16.De Zwart IM, De Roos A. MRI for the evaluation of gastric physiology. Eur Radiol. 2010;20(11):2609–2616. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s00330-010-1850-3&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20585784&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F23%2F2022.06.21.22276694.atom) 17. 17.Buisman WJ, Mauritz FA, Westerhuis WE, Gilja OH, van der Zee DC, van Herwaarden-Lindeboom MYA. Evaluation of gastric volumes: Comparison of 3-D ultrasound and magnetic resonance imaging. Ultrasound Med Biol. 2016;42(7):1423–1430. 18. 18.Menys A, Keszthelyi D, Fitzke H, et al. A magnetic resonance imaging study of gastric motor function in patients with dyspepsia associated with Ehlers-Danlos Syndrome-Hypermobility Type: A feasibility study. Neurogastroenterol Motil. 2017;29(9):1–9. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/nmo.13089&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28418103&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F23%2F2022.06.21.22276694.atom) 19. 19.Holt S, Cervantes J, Wilkinson AA, Kirk Wallace JH. Measurement of gastric emptying rate in humans by real-time ultrasound. Gastroenterology. 1986;90(4):918–923. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=3512357&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F23%2F2022.06.21.22276694.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1986A488600014&link_type=ISI) 20. 20.Distrutti E, Azpiroz F, Soldevilla A, Malagelada J. Gastric wall tension determines perception of gastric distention. Gastroenterology. 1999;116(5):1035–1042. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0016-5085(99)70006-5&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=10220495&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F23%2F2022.06.21.22276694.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000079947900009&link_type=ISI) 21. 21.Piessevaux H, Tack J, Wilmer A, Coulie B, Geubel A, Janssens J. Perception of changes in wall tension of the proximal stomach in humans. Gut. 2001;49(2):203–208. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NjoiZ3V0am5sIjtzOjU6InJlc2lkIjtzOjg6IjQ5LzIvMjAzIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjIvMDYvMjMvMjAyMi4wNi4yMS4yMjI3NjY5NC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 22. 22.Tack J, Caenepeel P, Corsetti M, Janssens J. Role of tension receptors in dyspeptic patients with hypersensitivity to gastric distention. Gastroenterology. 2004;127(4):1058–1066. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1053/j.gastro.2004.07.005&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15480984&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F23%2F2022.06.21.22276694.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000224604400011&link_type=ISI) 23. 23.Hoad CL, Parker H, Hudders N, et al. Measurement of gastric meal and secretion volumes using magnetic resonance imaging. Phys Med Biol. 2015;60(3):1367–1383. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1088/0031-9155/60/3/1367&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25592405&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F23%2F2022.06.21.22276694.atom) 24. 24.Fruehauf H, Goetze O, Steingoetter A, et al. Intersubject and intrasubject variability of gastric volumes in response to isocaloric liquid meals in functional dyspepsia and health. Neurogastroenterol Motil. 2007;19(7):553–561. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1365-2982.2007.00904.x&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17593136&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F23%2F2022.06.21.22276694.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000247477500003&link_type=ISI) 25. 25.Min YW, Min BH, Kim S, Choi D, Rhee PL. Effect of DA-9701 on gastric motor function assessed by magnetic resonance imaging in healthy volunteers: A randomized, double-blind, placebo-controlled trial. PLoS One. 2015;10(9):1–13. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0131080&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25919488&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F23%2F2022.06.21.22276694.atom) 26. 26.Banerjee S, Pal A, Fox M. Volume and position change of the stomach during gastric accommodation and emptying: A detailed three-dimensional morphological analysis based on MRI. Neurogastroenterol Motil. 2020;32(8). 27. 27.Taubin G. Curve and surface smoothing without shrinkage. In: Proceedings of IEEE International Conference on Computer Vision. IEEE Comput. Soc. Press; :852–857. doi:10.1109/ICCV.1995.466848 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1109/ICCV.1995.466848&link_type=DOI) 28. 28.E. H. Blending and Implicit with a Parametric Surface. In: Computer Aided Geometric Design. ; 1995:825–835. 29. 29.Liao D, Zhao J, Gregersen H. Three-dimensional geometry analysis of the stomach in type II diabetic GK rats. Diabetes Res Clin Pract. 2006;71(1):1–13. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.diabres.2005.05.016&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16054265&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F23%2F2022.06.21.22276694.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000234822100001&link_type=ISI) 30. 30.Schmitz A, Kellenberger CJ, Liamlahi R, Studhalter M, Weiss M. Gastric emptying after overnight fasting and clear fluid intake: A prospective investigation using serial magnetic resonance imaging in healthy children. Br J Anaesth. 2011;107(3):425–429. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bja/aer167&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21676893&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F23%2F2022.06.21.22276694.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000293910400018&link_type=ISI) 31. 31.Fruehauf H, Menne D, Kwiatek MA, et al. Inter-observer reproducibility and analysis of gastric volume measurements and gastric emptying assessed with magnetic resonance imaging. Neurogastroenterol Motil. 2011;23(9):854–861. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1365-2982.2011.01743.x&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21740482&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F23%2F2022.06.21.22276694.atom) 32. 32.Costanzo F, Brasseur JG. The invalidity of the Laplace law for biological vessels and of estimating elastic modulus from total stress vs. strain: a new practical method. Math Med Biol. 2015;32(1):1–37. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/imammb/dqt020&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24071531&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F23%2F2022.06.21.22276694.atom) 33. 33.Shaw JM, Hamad NM, Coleman TJ, et al. Intra-abdominal pressures during activity in women using an intra-vaginal pressure transducer. J Sports Sci. 2014;32(12):1176–1185. 34. 34.El-Serag HB, Tran T, Richardson P, Ergun G. Anthropometric correlates of intragastric pressure. Scand J Gastroenterol. 2006;41(8):887–891. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1080/00365520500535402&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16803686&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F23%2F2022.06.21.22276694.atom) 35. 35.H W. Ggplot2: Elegant Graphics for Data Analysis. Springer-Verlag New York; 2016. [https://ggplot2.tidyverse.org](https://ggplot2.tidyverse.org). 36. 36.Deepankar Datta. blandr: a Bland-Altman Method Comparison package for R. 2017. 37. 37.Revelle W. psych: Procedures for Psychological, Psychometric, and Personality Research. 2022. [https://cran.r-project.org/package=psych](https://cran.r-project.org/package=psych). 38. 38.Chitkara DK, Camilleri M, Zinsmeister AR, et al. Gastric sensory and motor dysfunction in adolescents with functional dyspepsia. J Pediatr. 2005;146(4):500–505. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jpeds.2004.11.031&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15812453&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F23%2F2022.06.21.22276694.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000228414200012&link_type=ISI) 39. 39.Grimm M, Koziolek M, Kühn J-P, Weitschies W. Interindividual and intraindividual variability of fasted state gastric fluid volume and gastric emptying of water. Eur J Pharm Biopharm. 2018;127:309–317. 40. 40.Khalaf A, Hoad CL, Menys A, et al. MRI assessment of the postprandial gastrointestinal motility and peptide response in healthy humans. Neurogastroenterol Motil. 2018;30(1):1–9. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/nmo.13351&link_type=DOI) 41. 41.Bertoli D, Steinkohl E, Mark EB, Brock C, Drewes AM, Frøkjær JB. Quantification of gastric emptying with magnetic resonance imaging in healthy volunteers: A systematic review. Neurogastroenterol Motil. March 2022. 42. 42.Boulby P, Gowland P, Adams V, Spiller RC. Use of echo planar imaging to demonstrate the effect of posture on the intragastric distribution and emptying of an oil/water meal. Neurogastroenterol Motil. 1997;9(1):41–47. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1046/j.1365-2982.1997.d01-6.x&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=9058391&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F23%2F2022.06.21.22276694.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1997WJ88100008&link_type=ISI) 43. 43.Carbonea SF, Tanganelli I, Capodiventoc S, Ricci V, Volterrani L. Magnetic resonance imaging in the evaluation of the gastric emptying and antral motion: Feasibility and reproducibility of a fast not invasive technique. Eur J Radiol. 2010;75(2):212–214. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ejrad.2009.04.071&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19481400&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F23%2F2022.06.21.22276694.atom) 44. 44.Storlid EL, Hausken T, Lied GA, Gilja OH, Hatlebakk JG. Gastric accommodation in healthy subjects studied by ultrasound, manometry, and impedancemetry. Neurogastroenterol Motil. 2018;30(4):e13249. 45. 45.Liao D, Gregersen H, Hausken T, Gilja OH, Mundt M, Kassab G. Analysis of surface geometry of the human stomach using real-time 3-D ultrasonography in vivo. Neurogastroenterol Motil. 2004;16(3):315–324. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1365-2982.2004.00522.x&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15198654&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F23%2F2022.06.21.22276694.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000222072100008&link_type=ISI) 46. 46.Liao D, Duch BU, Stødkilde-Jørgensen H, Zeng YJ, Gregersen H, Kassab GS. Tension and Stress Calculations in a 3-D Fourier Model of Gall Bladder Geometry Obtained from MR Images. Ann Biomed Eng. 2004;32(5):744–755. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15171628&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F23%2F2022.06.21.22276694.atom) 47. 47.Frokjaer j. b., Liao D, Bergmann A, et al. Three-dimensional biomechanical properties of the human rectum evaluated with magnetic resonance imaging. Neurogastroenterol Motil. 2005;17(4):531–540. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1365-2982.2005.00647.x&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16078942&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F23%2F2022.06.21.22276694.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000230929900008&link_type=ISI) 48. 48. Velde S Vande, Notebaert A, Meersschaut V, Herregods N, Van Winckel M, Van Biervliet S. Colon transit time in healthy children and adolescents. Int J Colorectal Dis. 2013;28(12):1721–1724. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s00384-013-1750-5&link_type=DOI) 49. 49.Sugita M, Matsumoto M, Tsukano Y, Fukunaga C, Yamamoto T. Gastric emptying time after breakfast in healthy adult volunteers using ultrasonography. J Anesth. 2019;33(6):697–700. 50. 50.Zou KH, Warfield SK, Bharatha A, et al. Statistical validation of image segmentation quality based on a spatial overlap index1. Acad Radiol. 2004;11(2):178–189. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S1076-6332(03)00671-8&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=14974593&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F23%2F2022.06.21.22276694.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000188549900009&link_type=ISI) 51. 51.Koo TK, Li MY. A Guideline of Selecting and Reporting Intraclass Correlation Coefficients for Reliability Research. J Chiropr Med. 2016;15(2):155–163. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jcm.2016.02.012&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27330520&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F23%2F2022.06.21.22276694.atom) 52. 52.Rostas JW, Mai TT, Richards WO. Gastric Motility Physiology and Surgical Intervention. Surg Clin North Am. 2011;91(5):983–999. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21889025&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F23%2F2022.06.21.22276694.atom) 53. 53.Freeman R. Diabetic autonomic neuropathy. In: ; 2014:63–79. doi:10.1016/B978-0-444-53480-4.00006-0 54. 54.Van Den Abeele J, Brouwers J, Tack J, Augustijns P. Exploring the link between gastric motility and intragastric drug distribution in man. Eur J Pharm Biopharm. 2017;112:75–84. 55. 55.Treier R, Steingoetter A, Weishaupt D, et al. Gastric motor function and emptying in the right decubitus and seated body position as assessed by magnetic resonance imaging. J Magn Reson Imaging. 2006;23(3):331–338. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/jmri.20507&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16463302&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F06%2F23%2F2022.06.21.22276694.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000236058400010&link_type=ISI) 56. 56.Lu KH, Cao J, Oleson ST, Powley TL, Liu Z. Contrast-Enhanced Magnetic Resonance Imaging of Gastric Emptying and Motility in Rats. IEEE Trans Biomed Eng. 2017;64(11):2546–2554. 57. 57.de Zwart IM, Haans JJL, Verbeek P, Eilers PHC, de Roos A, Masclee AAM. Gastric accommodation and motility are influenced by the barostat device: assessment with magnetic resonance imaging. Am J Physiol Liver Physiol. 2007;292(1):G208–G214. [1]: /embed/graphic-3.gif [2]: /embed/graphic-4.gif [3]: /embed/inline-graphic-1.gif [4]: /embed/graphic-5.gif [5]: /embed/graphic-6.gif [6]: /embed/graphic-7.gif [7]: /embed/graphic-8.gif [8]: /embed/graphic-9.gif [9]: /embed/graphic-10.gif [10]: /embed/graphic-11.gif [11]: /embed/graphic-14.gif [12]: /embed/inline-graphic-2.gif [13]: /embed/graphic-15.gif