A novel ultrasound software system for lumbar level identification in obstetric patients.

PURPOSE: Using machine learning, we developed a proprietary ultrasound software called the Spine Level Identification (SLIDE) system, which automatically identifies lumbar landmarks in real time as the operator slides the transducer over the lumber spine. Here, we assessed the agreement between SLIDE and manual palpation and traditional lumbar ultrasound (LUS) for determining the primary target L3-4 interspace. METHODS: Upon institutional ethics approval and informed consent, 76 healthy term parturients scheduled for elective Caesarean delivery were recruited. The L3-4 interspace was identified by manual palpation and then by the SLIDE method. The reference standard was located using traditional LUS by an experienced operator. The primary outcome was the L3-4 interspace identification agreement of manual palpation and SLIDE with the reference standard, as percentage agreement and Gwet's agreement coefficient (AC1). RESULTS: The raw agreement was 70% with Gwet's agreement coefficient (AC1) = 0.59 (95% confidence interval [CI], 0.41 to 0.77) for manual palpation and 84% with Gwet's AC1 = 0.82 (95% CI, 0.70 to 0.93) for SLIDE. When the levels differ from the reference, the manual palpation method identified L2-3 more often than L4-5 while the SLIDE method identified equally above or below L3-4. The SLIDE system had greater agreement than palpation in locating L3-4 and all other lumber interspaces after controlling for body mass index (adjusted odds ratio, 2.99; 95% CI, 1.21 to 8.7; P = 0.02). CONCLUSION: The SLIDE system had higher agreement with traditional ultrasound than manual palpation did in identifying L3-4 and all other lumber interspaces after adjusting for BMI in healthy term obstetric patients. Future studies should examine factors that affect agreement and ways to improve SLIDE for clinical integration. STUDY REGISTRATION: www. CLINICALTRIALS: gov (NCT02982317); registered 5 December 2016.

RéSUMé: OBJECTIF: À l'aide de l'apprentissage automatique, nous avons développé un logiciel d'échographie propriétaire appelé SLIDE (pour Spine Level Identification, c.-à-d. système d'identification du niveau vertébral), qui identifie automatiquement les points de repère lombaires en temps réel lorsque l'opérateur fait passer le transducteur sur la colonne lombaire. Ici, nous avons évalué l'agrément entre le SLIDE et la palpation manuelle et l'échographie lombaire traditionnelle pour déterminer l'espace intervertébral cible principal L3­L4. MéTHODE: Après avoir obtenu l'approbation du comité d'éthique de l'établissement et le consentement éclairé, 76 parturientes en bonne santé et à terme devant bénéficier d'un accouchement par césarienne programmée ont été recrutées. L'espace intervertébral L3­L4 a été identifié par palpation manuelle puis avec le logiciel SLIDE. L'étalon de référence a été localisé à l'aide d'une échographie lombaire traditionnelle par un opérateur expérimenté. Le critère d'évaluation principal était l'agrément entre l'identification de l'espace intervertébral L3­L4 par palpation manuelle et par logiciel SLIDE avec l'étalon de référence, en pourcentage d'agrément et coefficient d'agrément de Gwet (CA1). RéSULTATS: L'agrément brut était de 70 % avec le coefficient d'agrément de Gwet (CA1) = 0,59 (intervalle de confiance [IC] à 95 %, 0,41 à 0,77) pour la palpation manuelle et de 84 % avec le CA1 de Gwet = 0,82 (IC 95 %, 0,70 à 0,93) pour le logiciel SLIDE. Lorsque les niveaux lombaires différaient de la référence, la méthode de palpation manuelle a identifié L2­L3 plus souvent que L4­L5, tandis que la méthode SLIDE a identifié les vertèbres supérieures ou inférieures à L3­L4 de manière égale. Le système SLIDE a affiché un agrément plus important que la palpation pour localiser L3­L4 et tous les autres espaces intervertébraux lombaires après ajustement pour tenir compte de l'indice de masse corporelle (rapport de cotes ajusté, 2,99; IC 95 %, 1,21 à 8,7; P = 0,02). CONCLUSION: Le système SLIDE avait affiché un agrément plus élevé avec l'échographie traditionnelle que la palpation manuelle pour identifier le niveau L3­L4 et tous les autres espaces intervertébraux lombaires après ajustement pour tenir compte de l'IMC chez les patientes obstétricales à terme en bonne santé. Une étude future devrait examiner les facteurs qui affectent l'agrément et les moyens d'améliorer le logiciel SLIDE pour une intégration clinique. ENREGISTREMENT DE L'éTUDE: www.clinicaltrials.gov (NCT02982317); enregistrée le 5 décembre 2016.

A comparative study on position and paramedian neuraxial access on healthy volunteers using three-dimensional models registered to lumbar spine ultrasound.

PURPOSE: Optimizing patient position and needle puncture site are important factors for successful neuraxial anesthesia. Two paramedian approaches are commonly utilized and we sought to determine whether variations of the seated position would increase the chance of puncture success. METHODS: We simulated paramedian needle passes on three-dimensional lumbar spine models registered to volumetric ultrasound data acquired from ten healthy volunteers in three different positions: 1) prone; 2) seated with thoracic and lumbar flexion; and 3) seated as in position 2, but with a 10° dorsal tilt. Simulated paramedian needle passes from the right side performed on validated models were used to determine L2-3 and L3-4 neuraxial target size and success. We selected two paramedian puncture sites according to standard anesthesia textbook descriptions: 10 mm lateral and 10 mm caudal from inferior edge of the superior spinous process as described by Miller, and 10 mm lateral from the superior edge of the inferior spinous process as described by Barash. RESULTS: A significant increase in the area available for dural puncture was found in the L2-3 (61-62 mm2) and L3-4 (76-79 mm2) vertebral levels for all seated positions relative to the prone position (P < 0.001). Similarly, a significant increase in the total number of successful punctures was found in the L2-3 (77-79) and L3-4 (119-120) vertebral levels for all seated positions relative to the prone position (P < 0.001). No differences were found between seated positions. The Barash puncture site achieved a higher number of successful punctures than the Miller puncture site in both the L2-3 (19) and L3-4 (84) vertebral levels (P < 0.001). CONCLUSION: An added dorsal table tilt did not increase puncture success in the seated position. The landmarks for puncture site described by Barash resulted in significantly more successful punctures compared with those described by Miller in all positions.

RéSUMé: OBJECTIF: L'optimisation de la position du patient et celle du site de ponction de l'aiguille sont des facteurs importants pour la réussite d'une anesthésie neuraxiale. Deux approches paramédianes sont fréquemment utilisées et nous avons tenté de déterminer si des variations de la position assise augmenteraient la probabilité de réussite de la ponction. MéTHODE: Nous avons simulé les passages paramédians de l'aiguille sur des modèles tridimensionnels de la colonne lombaire adaptés à partir de données d'échographie volumétriques acquises auprès de dix volontaires sains placés dans trois positions différentes : 1) couché sur le ventre; 2) assis en flexion thoraco-lombaire; et 3) assis comme en position 2, mais avec une inclinaison dorsale de 10°. Les passages paramédians simulés de l'aiguille du côté droit réalisés sur des modèles validés ont été utilisés pour déterminer la taille des cibles neuraxiales L2­3 et L3­4 ainsi que la réussite de la ponction. Nous avons sélectionné deux sites de ponction paramédians selon les descriptions de deux manuels d'anesthésie standard, soit 10 mm en latéral et 10 mm en caudal depuis le bord inférieur de l'apophyse épineuse supérieure tel que décrit par celui de Miller, et 10 mm en latéral depuis le bord supérieur de l'apophyse épineuse inférieure, tel que décrit par celui de Barash. RéSULTATS: Une augmentation significative de la surface disponible pour la ponction durale a été observée aux niveaux vertébraux L2­3 (61­62 mm2) et L3­4 (76­79 mm2) dans les deux positions assises par rapport à la position ventrale (P < 0,001). De la même manière, nous avons observé une augmentation significative du nombre total de ponctions durales réussies aux niveaux vertébraux L2­3 (77­79) et L3­4 (119­120) dans les deux positions assises par rapport à la position ventrale (P < 0,001). Aucune différence n'a été observée entre les deux positions assises. Le site de ponction selon le manuel de Barash a permis la réalisation d'un nombre plus élevé de ponctions réussies que le site de ponction selon celui de Miller, tant au niveau vertébral L2­3 (19) qu'au niveau L3­4 (84) (P < 0,001). CONCLUSION: L'ajout d'une inclinaison du plan dorsal n'a pas augmenté le taux de réussite de la ponction en position assise. Les repères utilisés pour le site de ponction décrits par le manuel de Barash ont entraîné un nombre significativement plus élevé de ponctions réussies que ceux décrits par celui de Miller, toutes positions confondues.

The Neyman Pearson detection of microsaccades with maximum likelihood estimation of parameters.

Despite the fact that the velocity threshold method is widely applied, the detection of microsaccades continues to be a challenging problem, due to gaze-tracking inaccuracy and the transient nature of microsaccades. Important parameters associated with a saccadic event, e.g., saccade duration, amplitude, and maximum velocity, are sometimes imprecisely estimated, which may lead to biases in inferring the roles of microsaccades in perception and cognition. To overcome the biases and have a better detection algorithm for microsaccades, we propose a novel statistical model for the tracked gaze positions during eye fixations. In this model, we incorporate a parametrization that has been previously applied to model saccades, which allows us to veridically capture the velocity profile of saccadic eye movements. Based on our model, we derive the Neyman Pearson Detector (NPD) for saccadic events. Implemented in conjunction with the maximum likelihood estimation method, our NPD can detect a saccadic event and estimate all parameters simultaneously. Because of its adaptive nature and its statistical optimality, our NPD method was able to better detect microsaccades in some datasets when compared with a recently proposed state-of-the-art method based on convolutional neural networks. NPD also yielded comparable performance with a recently developed Bayesian algorithm, with the added benefit of modeling a more biologically veridical velocity profile of the saccade. As opposed to these algorithms, NPD can lend itself better to online saccade detection, and thus has potential for human-computer interaction applications. Our algorithm is publicly available at https://github.com/hz-zhu/NPD-micro-saccade-detection.

Design of high energy laser pulse delivery in a multimode fiber for photoacoustic tomography.

In photoacoustic tomography (PAT), delivering high energy pulses through optical fiber is critical for achieving high quality imaging. A fiber coupling scheme with a beam homogenizer is demonstrated for coupling high energy pulses in a single multimode fiber. This scheme can benefit PAT applications that require miniaturized illumination or internal illumination with a small fiber. The beam homogenizer is achieved by using a cross cylindrical lens array, which provides a periodic spatial modulation on the phase of the input light. Thus the lens array acts as a phase grating which diffracts the beam into a 2D diffraction pattern. Both theoretical analysis and experiments demonstrate that the focused beam can be split into a 2D spot array that can reduce the peak power on the fiber tip surface and thus enhance the coupling performance. The theoretical analysis of the intensity distribution of the focused beam is carried out by Fourier optics. In experiments, coupled energy at 48 mJ/pulse and 60 mJ/pulse have been achieved and the corresponding coupling efficiency is 70% and 90% in a 1000-µm and a 1500-µm-core-diameter fiber, respectively. The high energy pulses delivered by the multimode fiber are further tested for PAT imaging in phantoms. PAT imaging of a printed dot array shows a large illumination area of 7 cm2 under 5 mm thick chicken breast tissue. In vivo imaging is also demonstrated on the human forearm. The large improvement in coupling energy can potentially benefit PAT with single fiber delivery to achieve large area imaging and deep penetration detection.

Fabrication of Circuits on Flexible Substrates Using Conductive SU-8 for Sensing Applications.

This article describes a new low-cost rapid microfabrication technology for high-density interconnects and passive devices on flexible substrates for sensing applications. Silver nanoparticles with an average size of 80 nm were used to create a conductive SU-8 mixture with a concentration of wt 25%. The patterned structures after hard baking have a sheet resistance of 11.17 Ω /□. This conductive SU-8 was used to pattern planar inductors, capacitors and interconnection lines on flexible Kapton film. The conductive SU-8 structures were used as a seed layer for a subsequent electroplating process to increase the conductivity of the devices. Examples of inductors, resistor-capacitor (RC) and inductor-capacitor (LC) circuits, interconnection lines and a near-field communication (NFC) antenna are presented as a demonstration. As an example of high-resolution miniaturization, we fabricated microinductors having line widths of 5 µ m. Mechanical bending tests were successful down to a 5 mm radius. To the best of the authors' knowledge, this is the first report of conductive SU-8 used to fabricate such planar devices and the first on flexible substrates. This is a proof of concept that this fabrication approach can be used as an alternative for microfabrication of planar passive devices on flexible substrates.

Broadening the detection view of 2D photoacoustic tomography using two linear array transducers.

Although commercial linear array transducers are widely used in clinical ultrasound, their application in photoacoustic tomography (PAT) is still limited due to the limited-view problem that restricts the image quality. In this paper, we propose a simple approach to address the limited-view problem in 2D by using two linear array transducers to receive PAT signal from different orientations. The positions of the two transducers can be adjusted to fit the specific geometry of an imaging site. This approach is made possible by using a new calibration method, where the relative position between the two transducers can be calibrated using ultrasound by transmitting ultrasound wave with one transducer while receiving with the other. The calibration results are then applied in the subsequent PAT imaging to incorporate the detected acoustic signals from both transducers and thereby increase the detection view. In this calibration method, no calibration phantom is required which largely simplifies and shortens the process. The efficacy of the calibration and improvement on the PAT image quality are demonstrated through phantom studies and in vivo imaging.

MR elastography and diffusion-weighted imaging of ex vivo prostate cancer: quantitative comparison to histopathology.

The purpose of this work was (1) to develop a magnetic resonance elastography (MRE) system for imaging of the ex vivo human prostate and (2) to assess the diagnostic power of mono-frequency and multi-frequency MRE and diffusion weighted imaging (DWI) alone and combined as correlated with histopathology in a patient study. An electromagnetic driver was designed specifically for MRE studies in small-bore MR scanners. Ex vivo prostate specimens (post-fixation) of 14 patients who underwent radical prostatectomy were imaged with MRE at 7 T (nine cases had DWI). In six patients, the MRE examination was performed at three frequencies (600, 800, 1000 Hz) to extract the power-law exponent Gamma. The images were registered to wholemount pathology slides marked with the Gleason score. The areas under the receiver-operator-characteristic curves (AUC) were calculated. The methods were validated in a phantom study and it was demonstrated that (i) the driver does not interfere with the acquisition process and (ii) the driver can generate amplitudes greater than 100 µm for frequencies less than 1 kHz. In the quantitative study, cancerous tissue with Gleason score at least 3 + 3 was distinguished from normal tissue in the peripheral zone (PZ) with an average AUC of 0.75 (Gd ), 0.75 (Gl ), 0.70 (Gamma-Gd ), 0.68 (apparent diffusion coefficient, ADC), and 0.82 (Gd + Gl + ADC). The differentiation between PZ and central gland was modest for Gd (p < 0.07), Gl (p < 0.06) but not significant for Gamma (p < 0.2). A correlation of 0.4 kPa/h was found between the fixation time of the prostate specimen and the stiffness of the tissue, which could affect the diagnostic power results. DWI and MRE may provide complementary information; in fact MRE performed better than ADC in distinguishing normal from cancerous tissue in some cases. Multi-frequency (Gamma) analysis did not appear to improve the results. However, in light of the effect of tissue fixation, the clinical implication of our results may be inconclusive and more experiments are needed.

MR elastography of prostate cancer: quantitative comparison with histopathology and repeatability of methods.

The purpose of this work was to assess trans-perineal prostate magnetic resonance elastography (MRE) for (1) repeatability in phantoms/volunteers and (2) diagnostic power as correlated with histopathology in prostate cancer patients. The three-dimensional (3D) displacement field was obtained using a fractionally encoded gradient echo sequence using a custom-made transducer. The repeatability of the method was assessed based on three repeat studies and by changing the driving frequency by 3% in studies on a phantom and six healthy volunteers. Subsequently, 11 patients were examined with MRE prior to radical prostatectomy. The areas under the receiver operating characteristic curves were calculated using a windowed voxel-to-voxel approach by comparing the 2D registered slides, masked with the Gleason score. For the repeatability study, the average intraclass correlation coefficient for elasticity images was 99% for repeat phantom studies, 98% for ±6 Hz phantom studies, 95% for volunteer repeat studies with 2 min acquisition time, 82% for ±2 Hz volunteer studies with 2 min acquisition time and 73% for repeat volunteer studies with 8 min acquisition time. For the patient study, the average elasticity was 8.2 ± 1.7 kPa in the prostate capsule, 7.5 ± 1.9 kPa in the peripheral zone (PZ), 9.7 ± 3.0 kPa in the central gland (CG) and 9.0 ± 3.4 kPa in the transition zone. In the patient study, cancerous tissue with Gleason score at least 3 + 3 was significantly (p < 0.05) different from normal tissue in 10 out of 11 cases with tumors in the PZ, and 6 out of 9 cases with tumors in the CG. However, the overall case-averaged area under the curve was 0.72 in the PZ and 0.67 in the CG. Cancerous tissue was not always stiffer than normal tissue. The inversion algorithm was sensitive to (i) vibration amplitude and displacement nodes and (ii) misalignment of the 3D wave field due to subject movement.

The Kidneys Are Not All Normal: Transplanted Kidneys and Their Speckle Distributions.

OBJECTIVE: Modelling ultrasound speckle to characterise tissue properties has generated considerable interest. As speckle is dependent on the underlying tissue architecture, modelling it may aid in tasks such as segmentation or disease detection. For the transplanted kidney, where ultrasound is used to investigate dysfunction, it is unknown which statistical distribution best characterises such speckle. This applies to the regions of the transplanted kidney: the cortex, the medulla and the central echogenic complex. Furthermore, it is unclear how these distributions vary by patient variables such as age, sex, body mass index, primary disease or donor type. These traits may influence speckle modelling given their influence on kidney anatomy. We investigate these two aims. METHODS: B-mode images from n = 821 kidney transplant recipients (one image per recipient) were automatically segmented into the cortex, medulla and central echogenic complex using a neural network. Seven distinct probability distributions were fitted to each region's histogram, and statistical analysis was performed. DISCUSSION: The Rayleigh and Nakagami distributions had model parameters that differed significantly between the three regions (p ≤ 0.05). Although both had excellent goodness of fit, the Nakagami had higher Kullbeck-Leibler divergence. Recipient age correlated weakly with scale in the cortex (Ω: ρ = 0.11, p = 0.004), while body mass index correlated weakly with shape in the medulla (m: ρ = 0.08, p = 0.04). Neither sex, primary disease nor donor type exhibited any correlation. CONCLUSION: We propose the Nakagami distribution be used to characterize transplanted kidneys regionally independent of disease etiology and most patient characteristics.

Automatic measurement of kidney dimensions in two-dimensional ultrasonography is comparable to expert sonographers.

Purpose: Length and width measurements of the kidneys aid in the detection and monitoring of structural abnormalities and organ disease. Manual measurement results in intra- and inter-rater variability, is complex and time-consuming, and is fraught with error. We propose an automated approach based on machine learning for quantifying kidney dimensions from two-dimensional (2D) ultrasound images in both native and transplanted kidneys. Approach: An nnU-net machine learning model was trained on 514 images to segment the kidney capsule in standard longitudinal and transverse views. Two expert sonographers and three medical students manually measured the maximal kidney length and width in 132 ultrasound cines. The segmentation algorithm was then applied to the same cines, region fitting was performed, and the maximum kidney length and width were measured. Additionally, single kidney volume for 16 patients was estimated using either manual or automatic measurements. Results: The experts resulted in length of 84.8±26.4 mm [95% CI: 80.0, 89.6] and a width of 51.8±10.5 mm [49.9, 53.7]. The algorithm resulted a length of 86.3±24.4 [81.5, 91.1] and a width of 47.1±12.8 [43.6, 50.6]. Experts, novices, and the algorithm did not statistically significant differ from one another (p>0.05). Bland-Altman analysis showed the algorithm produced a mean difference of 2.6 mm (SD = 1.2) from experts, compared to novices who had a mean difference of 3.7 mm (SD = 2.9 mm). For volumes, mean absolute difference was 47 mL (31%) consistent with ∼1 mm error in all three dimensions. Conclusions: This pilot study demonstrates the feasibility of an automatic tool to measure in vivo kidney biometrics of length, width, and volume from standard 2D ultrasound views with comparable accuracy and reproducibility to expert sonographers. Such a tool may enhance workplace efficiency, assist novices, and aid in tracking disease progression.

3-D Ultrafast Shear Wave Absolute Vibro-Elastography Using a Matrix Array Transducer.

Real-time ultrasound imaging plays an important role in ultrasound-guided interventions. The 3-D imaging provides more spatial information compared to conventional 2-D frames by considering the volumes of data. One of the main bottlenecks of 3-D imaging is the long data acquisition time, which reduces practicality and can introduce artifacts from unwanted patient or sonographer motion. This article introduces the first shear wave absolute vibro-elastography (S-WAVE) method with real-time volumetric acquisition using a matrix array transducer. In S-WAVE, an external vibration source generates mechanical vibrations inside the tissue. The tissue motion is then estimated and used in solving a wave equation inverse problem to provide the tissue elasticity. A matrix array transducer is used with a Verasonics ultrasound machine and a frame rate of 2000 volumes/s to acquire 100 radio frequency (RF) volumes in 0.05 s. Using plane wave (PW) and compounded diverging wave (CDW) imaging methods, we estimate axial, lateral, and elevational displacements over 3-D volumes. The curl of the displacements is used with local frequency estimation to estimate elasticity in the acquired volumes. Ultrafast acquisition extends substantially the possible S-WAVE excitation frequency range, now up to 800 Hz, enabling new tissue modeling and characterization. The method was validated on three homogeneous liver fibrosis phantoms and on four different inclusions within a heterogeneous phantom. The homogeneous phantom results show less than 8% (PW) and 5% (CDW) difference between the manufacturer values and the corresponding estimated values over a frequency range of 80-800 Hz. The estimated elasticity values for the heterogeneous phantom at 400-Hz excitation frequency show the average errors of 9% (PW) and 6% (CDW) compared to the provided average values by magnetic resonance elastography (MRE). Furthermore, both imaging methods were able to detect the inclusions within the elasticity volumes. An ex vivo study on a bovine liver sample shows less than 11% (PW) and 9% (CDW) difference between the estimated elasticity ranges by the proposed method and the elasticity ranges provided by MRE and acoustic radiation force impulse (ARFI).

3-D Transducer Mounted Shear Wave Absolute Vibro-Elastography: Proof of Concept.

Quantitative tissue stiffness characterization using ultrasound (US) has been shown to improve prostate cancer (PCa) detection in multiple studies. Shear wave absolute vibro-elastography (SWAVE) allows quantitative and volumetric assessment of tissue stiffness using external multifrequency excitation. This article presents a proof of concept of a first-of-a-kind 3-D hand-operated endorectal SWAVE system designed to be used during systematic prostate biopsy. The system is developed with a clinical US machine, requiring only an external exciter that can be mounted directly to the transducer. Subsector acquisition of radio frequency (RF) data allows imaging of shear waves with a high effective frame rate (up to 250 Hz). The system was characterized using eight different quality assurance phantoms. Due to the invasive nature of prostate imaging, at this early stage of development, validation of in vivo human tissue was instead carried out by intercostally scanning the livers of n = 7 healthy volunteers. The results are compared with 3-D magnetic resonance elastography (MRE) and an existing 3-D SWAVE system with a matrix array transducer (M-SWAVE). High correlations were found with MRE (99% in phantoms, 94% in liver data) and with M-SWAVE (99% in phantoms, 98% in liver data).

Multi-Frequency 3D Shear Wave Absolute Vibro-Elastography (S-WAVE) System for the Prostate.

This article describes a novel system for quantitative and volumetric measurement of tissue elasticity in the prostate using simultaneous multi-frequency tissue excitation. Elasticity is computed by using a local frequency estimator to measure the three-dimensional local wavelengths of steady-state shear waves within the prostate gland. The shear wave is created using a mechanical voice coil shaker which transmits simultaneous multi-frequency vibrations transperineally. Radio frequency data is streamed directly from a BK Medical 8848 transrectal ultrasound transducer to an external computer where tissue displacement due to the excitation is measured using a speckle tracking algorithm. Bandpass sampling is used that eliminates the need for an ultra-fast frame rate to track the tissue motion and allows for accurate reconstruction at a sampling frequency that is below the Nyquist rate. A roll motor with computer control is used to rotate the transducer and obtain 3D data. Two commercially available phantoms were used to validate both the accuracy of the elasticity measurements as well as the functional feasibility of using the system for in vivo prostate imaging. The phantom measurements were compared with 3D Magnetic Resonance Elastography (MRE), where a high correlation of 96% was achieved. In addition, the system has been used in two separate clinical studies as a method for cancer identification. Qualitative and quantitative results of 11 patients from these clinical studies are presented here. Furthermore, an AUC of 0.87±0.12 was achieved for malignant vs. benign classification using a binary support vector machine classifier trained with data from the latest clinical study with leave one patient out cross-validation.

Kidney Ultrasound for Nephrologists: A Review.

Ultrasound imaging is a key investigatory step in the evaluation of chronic kidney disease and kidney transplantation. It uses nonionizing radiation, is noninvasive, and generates real-time images, making it the ideal initial radiographic test for patients with abnormal kidney function. Ultrasound enables the assessment of both structural (form and size) and functional (perfusion and patency) aspects of kidneys, both of which are especially important as the disease progresses. Ultrasound and its derivatives have been studied for their diagnostic and prognostic significance in chronic kidney disease and kidney transplantation. Ultrasound is rapidly growing more widely accessible and is now available even in handheld formats that allow for bedside ultrasound examinations. Given the trend toward ubiquity, the current use of kidney ultrasound demands a full understanding of its breadth as it and its variants become available. We described the current applications and future directions of ultrasound imaging and its variants in the context of chronic kidney disease and transplantation in this review.

A Unified Representation of Control Logic in Human-Ultrasound Machine Interaction.

Advances in human-computer interaction (HCI) technologies have granted sonographers and radiologists a much improved user experience when operating different ultrasound (US) machines. Continued HCI improvements in US would benefit from a systematic study of the HCI control logic used in this domain. Such a study has not been presented previously and is the subject of this paper. We surveyed sonographers to determine the most frequently used controls in US machines. We standardized the representation of the US machine HCI control logic by using the unified modelling language (UML). We used UML diagrams to analyze the HCI control logic of 10 different cart-based US machines from several major manufacturers, and we discovered that the control logic for the most frequently used functions are identical. While this control logic does not follow an established standard, it has been commonly adopted. Using the UML for the visualization and formulation of control logic, we can target logically optimal interactions (whose operation steps cannot be further reduced), e.g., adjustment of B-mode gain, frequency and depth, and can derive methods to simplify logically sub-optimal interactions, e.g., the pointing and selecting operation, as well as image measurements.Our study provides insights into existing HCI approaches used in US machines and establishes a rigorous UML-based framework for future US machine design to improve interoperability, efficiency and ease-of-use.

Model-Based Quantitative Elasticity Reconstruction Using ADMM.

We introduce two model-based iterative methods to obtain shear modulus images of tissue using magnetic resonance elastography. The first method jointly finds the displacement field that best fits tissue displacement data and the corresponding shear modulus. The displacement satisfies a viscoelastic wave equation constraint, discretized using the finite element method. Sparsifying regularization terms in both shear modulus and displacement are used in the cost function minimized for the best fit. The second method extends the first method for multifrequency tissue displacement data. The formulated problems are bi-convex. Their solution can be obtained iteratively by using the alternating direction method of multipliers. Sparsifying regularizations and the wave equation constraint filter out sensor noise and compressional waves. Our methods do not require bandpass filtering as a preprocessing step and converge fast irrespective of the initialization. We evaluate our new methods in multiple in silico and phantom experiments, with comparisons with existing methods, and we show improvements in contrast to noise and signal-to-noise ratios. Results from an in vivo liver imaging study show elastograms with mean elasticity comparable to other values reported in the literature.

Breast Cancer Detection Using Multimodal Time Series Features From Ultrasound Shear Wave Absolute Vibro-Elastography.

In shear wave absolute vibro-elastography (S-WAVE), a steady-state multi-frequency external mechanical excitation is applied to tissue, while a time-series of ultrasound radio-frequency (RF) data are acquired. Our objective is to determine the potential of S-WAVE to classify breast tissue lesions as malignant or benign. We present a new processing pipeline for feature-based classification of breast cancer using S-WAVE data, and we evaluate it on a new data set collected from 40 patients. Novel bi-spectral and Wigner spectrum features are computed directly from the RF time series and are combined with textural and spectral features from B-mode and elasticity images. The Random Forest permutation importance ranking and the Quadratic Mutual Information methods are used to reduce the number of features from 377 to 20. Support Vector Machines and Random Forest classifiers are used with leave-one-patient-out and Monte Carlo cross-validations. Classification results obtained for different feature sets are presented. Our best results (95% confidence interval, Area Under Curve = 95%±1.45%, sensitivity = 95%, and specificity = 93%) outperform the state-of-the-art reported S-WAVE breast cancer classification performance. The effect of feature selection and the sensitivity of the above classification results to changes in breast lesion contours is also studied. We demonstrate that time-series analysis of externally vibrated tissue as an elastography technique, even if the elasticity is not explicitly computed, has promise and should be pursued with larger patient datasets. Our study proposes novel directions in the field of elasticity imaging for tissue classification.

Project SWAVE 2.0: An overview of the study design for multimodal placental image acquisition and alignment.

Development of non-invasive and in utero placenta imaging techniques can potentially identify biomarkers of placental health. Correlative imaging using multiple multiscale modalities is particularly important to advance the understanding of placenta structure, function and their relationship. The objective of the project SWAVE 2.0 was to understand human placental structure and function and thereby identify quantifiable measures of placental health using a multimodal correlative approach. In this paper, we present a multimodal image acquisition protocol designed to acquire and align data from ex vivo placenta specimens derived from both healthy and complicated pregnancies. Qualitative and quantitative validation of the alignment method were performed. The qualitative analysis showed good correlation between findings in the MRI, ultrasound and histopathology images. The proposed protocol would enable future studies on comprehensive analysis of placental anatomy, function and their relationship. ● An overview of a novel multimodal placental image acquisition protocol is presented. ● A co-registration method using surface markers and external fiducials is described. ● A preliminary correlative imaging analysis for a placenta specimen is presented.

SWAVE 2.0 Imaging of Placental Elasticity and Viscosity: Potential Biomarkers for Placenta-Mediated Disease Detection.

Pregnancy complications such as pre-eclampsia (PE) and intrauterine growth restriction (IUGR) are associated with structural and functional changes in the placenta. Different elastography techniques with an ability to assess the mechanical properties of tissue can identify and monitor the pathological state of the placenta. Currently available elastography techniques have been used with promising results to detect placenta abnormalities; however, limitations include inadequate measurement depth and safety concerns from high negative pressure pulses. Previously, we described a shear wave absolute vibro-elastography (SWAVE) method by applying external low-frequency mechanical vibrations to generate shear waves and studied 61 post-delivery clinically normal placentas to explore the feasibility of SWAVE for placental assessment and establish a measurement baseline. This next phase of the study, namely, SWAVE 2.0, improves the previous system and elasticity reconstruction by incorporating a multi-frequency acquisition system and using a 3-D local frequency estimation (LFE) method. Compared with its 2-D counterpart, the proposed system using 3-D LFE was found to reduce the bias and variance in elasticity measurements in tissue-mimicking phantoms. In the aim of investigating the potential of improved SWAVE 2.0 measurements to identify placental abnormalities, we studied 46 post-delivery placentas, including 26 diseased (16 IUGR and 10 PE) and 20 normal control placentas. By use of a 3.33-MHz motorized curved-array transducer, multi-frequency (80,100 and 120 Hz) elasticity measures were obtained with 3-D LFE, and both IUGR (15.30 ± 2.96 kPa, p = 3.35e-5) and PE (12.33 ± 4.88 kPa, p = 0.017) placentas were found to be significantly stiffer compared with the control placentas (8.32 ± 3.67 kPa). A linear discriminant analysis (LDA) classifier was able to classify between healthy and diseased placentas with a sensitivity, specificity and accuracy of 87%, 78% and 83% and an area under the receiver operating curve of 0.90 (95% confidence interval: 0.8-0.99). Further, the pregnancy outcome in terms of neonatal intensive care unit admission was predicted with a sensitivity, specificity and accuracy of 70%, 71%, 71%, respectively, and area under the receiver operating curve of 0.78 (confidence interval: 0.62-0.93). A viscoelastic characterization of placentas using a fractional rheological model revealed that the viscosity measures in terms of viscosity parameter n were significantly higher in IUGR (2.3 ± 0.21) and PE (2.11 ± 0.52) placentas than in normal placentas (1.45 ± 0.65). This work illustrates the potential relevance of elasticity and viscosity imaging using SWAVE 2.0 as a non-invasive technology for detection of placental abnormalities and the prediction of pregnancy outcomes.

U-LanD: Uncertainty-Driven Video Landmark Detection.

This paper presents U-LanD, a framework for automatic detection of landmarks on key frames of the video by leveraging the uncertainty of landmark prediction. We tackle a specifically challenging problem, where training labels are noisy and highly sparse. U-LanD builds upon a pivotal observation: a deep Bayesian landmark detector solely trained on key video frames, has significantly lower predictive uncertainty on those frames vs. other frames in videos. We use this observation as an unsupervised signal to automatically recognize key frames on which we detect landmarks. As a test-bed for our framework, we use ultrasound imaging videos of the heart, where sparse and noisy clinical labels are only available for a single frame in each video. Using data from 4,493 patients, we demonstrate that U-LanD can exceedingly outperform the state-of-the-art non-Bayesian counterpart by a noticeable absolute margin of 42% in R2 score, with almost no overhead imposed on the model size.