Preoperative prediction of lymphovascular and perineural invasion in gastric cancer using spectral computed tomography imaging and machine learning

BACKGROUND Lymphovascular invasion(LVI)and perineural invasion(PNI)are important prognostic factors for gastric cancer(GC)that indicate an increased risk of metastasis and poor outcomes.Accurate preoperative prediction of LVI/PNI status could help clinicians identify high-risk patients and guide treatment deci-sions.However,prior models using conventional computed tomography(CT)images to predict LVI or PNI separately have had limited accuracy.Spectral CT provides quantitative enhancement parameters that may better capture tumor invasion.We hypothesized that a predictive model combining clinical and spectral CT parameters would accurately preoperatively predict LVI/PNI status in GC patients.AIM To develop and test a machine learning model that fuses spectral CT parameters and clinical indicators to predict LVI/PNI status accurately.METHODS This study used a retrospective dataset involving 257 GC patients(training cohort,n=172;validation cohort,n=85).First,several clinical indicators,including serum tumor markers,CT-TN stages and CT-detected extramural vein invasion(CT-EMVI),were extracted,as were quantitative spectral CT parameters from the delineated tumor regions.Next,a two-step feature selection approach using correlation-based methods and information gain ranking inside a 10-fold cross-validation loop was utilized to select informative clinical and spectral CT parameters.A logistic regression(LR)-based nomogram model was subsequently constructed to predict LVI/PNI status,and its performance was evaluated using the area under the receiver operating characteristic curve(AUC).RESULTS In both the training and validation cohorts,CT T3-4 stage,CT-N positive status,and CT-EMVI positive status are more prevalent in the LVI/PNI-positive group and these differences are statistically significant(P<0.05).LR analysis of the training group showed preoperative CT-T stage,CT-EMVI,single-energy CT values of 70 keV of venous phase(VP-70 keV),and the ratio of standardized iodine concentration of equilibrium phase(EP-NIC)were independent influencing factors.The AUCs of VP-70 keV and EP-NIC were 0.888 and 0.824,respectively,which were slightly greater than those of CT-T and CT-EMVI(AUC=0.793,0.762).The nomogram combining CT-T stage,CT-EMVI,VP-70 keV and EP-NIC yielded AUCs of 0.918(0.866-0.954)and 0.874(0.784-0.936)in the training and validation cohorts,which are significantly higher than using each of single independent factors(P<0.05).CONCLUSION The study found that using portal venous and EP spectral CT parameters allows effective preoperative detection of LVI/PNI in GC,with accuracy boosted by integrating clinical markers.

Multispectral smart window:Dynamic light modulation and electromagnetic microwave shielding

A novel multispectral smart window has been proposed,which features dynamic modulation of light transmittance and effective shielding against electromagnetic microwave radiation.This design integrates liquid crystal dynamic scattering and dye doping techniques,enabling the dual regulation of transmittance and scattering within a singlelayer smart window.Additionally,the precise control of conductive film thickness ensures the attainment of robust microwave signal shielding.We present a theoretical model for ion movement in the presence of an alternating electric field,along with a novel approach to manipulate negative dielectric constant.The proposed model successfully enables a rapid transition between light transparent,absorbing and haze states,with an optimum drive frequency adjustable to approximately 300 Hz.Furthermore,the resistive design of the conductive layer effectively mitigates microwave radiation within the 2−18 GHz range.These findings offer an innovative perspective for future advancements in environmental construction.

Rare Bilateral Radius and Ulna Fractures: A Case Report

Introduction: Forearm fractures involving both the radius and ulna, present distinctive challenges in orthopedic trauma management. This case report explores the complexities of a rare case of bilateral fractures, emphasizing the importance of a comprehensive approach for optimal outcomes. Case Presentation: This report presents the case of a 40-year-old motorcyclist who was involved in a road accident and arrived at the emergency department shortly after the incident. Clinical examination revealed closed fractures in both forearms with mild swelling and severe tenderness. Fortunately, no neurovascular issues or compartment syndrome were detected. The patient was administered intravenous fluids and effective pain relief. Forearm splints were applied and the arms were elevated using pillows. Radiographs showed fractures in the right distal radius and ulna, as well as the left mid-shaft radius and ulna. Under general anesthesia, open reduction and internal fixation were per-formed using plates and screws for all four fractures. Postoperatively, the patient’s arms were immobilized and elevated, with gradual recovery during follow-up appointments over six months, eventually achieving full function without complaints. Conclusion: This article highlights the possibility, although infrequent, of bilateral shaft fractures of the ulna and radius. With thorough clinical assessment and radiological investigation, such fractures can be well-defined. The required definitive treatment, which usually involves open reduction and internal fixation, can be performed.

Enhancing silicon spectral emission in LIBS using Tesla coil discharge

Laser-induced breakdown spectroscopy(LIBS)is a powerful technique for elemental analysis,offering rapid analysis,minimal sample preparation,wide elemental coverage,and portability.To enhance the detection sensitivity of LIBS,increasing the spectral emission intensity is crucial.This paper explores the use of Tesla coil(TC)discharge as an alternative to spark discharge in silicon LIBS.The study examines the influence of TC discharge on both time-integrated and timeresolved spectra,with and without TC discharge;the corresponding electron temperature and density are obtained.The results show that TC discharge significantly amplifies the spectral intensity,improving signal sensitivity in LIBS analysis.Specifically,in the laser energy range from 7.4 to 24.0 mJ,TC discharge increased the average spectral line intensities of Si(II)385.60 nm and Si(I)390.55 nm by factors of 8.4 and 5.1,respectively.Additionally,the average electron temperature and density were enhanced by approximately 3.2%and 4.2%,respectively,under TC discharge.The advantages of TC discharge include higher energy deposition,extended discharge duration,reduced electrode erosion,and enhanced safety.This research contributes to advancing LIBS technology and expanding its applications in various fields.

Wind-sand flow simulation and radius optimization of highway embankment under different vertical curve radius

It is of great practical value to explore the correlation between the vertical curve radius of desert highway and the increase of sand accumulation in local lines,and to select the appropriate vertical curve radius for reducing the risk of sand accumulation.In this study,three-dimensional models of desert highway embankments with different vertical curve radii were constructed,and Fluent software was used to simulate the wind-sand flow field and sand accumulation distribution of vertical curve embankments.The results show that:(1)Along the direction of the road,the concave and the convex vertical curve embankments have the effect of collecting and diverging the wind-sand flow,respectively.When the radius of the concave vertical curve is 3000 m,5000 m,8000 m,10000 m and 20000 m,the wind velocity in the middle of the vertical curve is 31.76%,22.58%,10.78%,10.53%and 10.44%,higher than that at both ends.When the radius of the convex vertical curve is 6500 m,8000 m,10000 m,20000 m and 30000 m,the wind velocity at both ends of the vertical curve is 14.06%,9.99%,6.14%,3.22%and 2.41%,higher than that in the middle.The diversion effect also decreases with the increase of the radius.(2)The conductivity of the concave and convex vertical curve embankments with different radii is greater than 1,which is the sediment transport roadbed.The conductivity increases with the increase of radius and gradually tends to be stable.When the radius of the concave and convex vertical curves reaches 8000 m and 20000 m respectively,the phenomenon of sand accumulation is no longer serious.Under the same radius condition,the concave vertical curve embankment is more prone to sand accumulation than the convex one.(3)Considering the strength of the collection and diversion of the vertical curve embankment with different radii,and the sand accumulation of the vertical curve embankment in the desert section of Wuma Expressway,the radius of the concave vertical curve is not less than 8000 m,and the radius of the convex vertical curve is not less than 20000 m,which can effectively reduce the sand accumulation of the vertical curve embankment.In the desert highway area,the research results of this paper can provide reference for the design of vertical curve to ensure the safe operation of desert highway.

Experimental investigation of spectral evolution in flash radiation by hypervelocity impact on aluminum plates

In this study,a series of hypervelocity impact tests were carried out based on a two-stage light gas gun,and the sequence spectrum and radiation evolution data of the impact products under different impact conditions were obtained.The diameter of the projectile is 3-5 mm,the impact velocity is 3.13-6.58 km/s,and the chamber pressure is 0.56-990 Pa.The spectrum of ejected debris cloud in the 250-310 nm band were obtained using a transient spectral measurement system and a multi-channel radiometer measurement system.The test results reveal that the flash radiation intensity increases as a power function with the kinetic energy of the impact.Furthermore,the peak value of the line spectrum decreases as the chamber vacuum degree increases,while the radiation width gradually expands.The line spectrum in the spectral characterization curve corresponds to the ejected debris clouds splitting phase,which does not produce significant line spectrum during material fragmentation and is dominated by the continuum spectrum produced by blackbody radiation.There will appear one or three characteristic peaks in the flash radiation time curve,the first and second peaks correspond to the penetration phase and the third peak corresponds to the expansion phase of the ejected debris clouds on the time scale,the first and second peaks are more sensitive to the chamber vacuum degree,and when the pressure is higher than 99 Pa,the first and second characteristic peaks will disappear.The radiant heat attenuation of the flash under different impact conditions is significantly different,the attenuation exponent has a power function relationship with the impact velocity and the chamber vacuum degree,while the attenuation exponent has a linear relationship with the diameter of the projectile,the specific expression of the attenuation exponent is obtained by fitting.The findings from this research can serve as a valuable reference for remote diagnostic technologies based on flash radiation characteristics.

S^(2)ANet:Combining local spectral and spatial point grouping for point cloud processing

Background Despite the recent progress in 3D point cloud processing using deep convolutional neural networks,the inability to extract local features remains a challenging problem.In addition,existing methods consider only the spatial domain in the feature extraction process.Methods In this paper,we propose a spectral and spatial aggregation convolutional network(S^(2)ANet),which combines spectral and spatial features for point cloud processing.First,we calculate the local frequency of the point cloud in the spectral domain.Then,we use the local frequency to group points and provide a spectral aggregation convolution module to extract the features of the points grouped by the local frequency.We simultaneously extract the local features in the spatial domain to supplement the final features.Results S^(2)ANet was applied in several point cloud analysis tasks;it achieved stateof-the-art classification accuracies of 93.8%,88.0%,and 83.1%on the ModelNet40,ShapeNetCore,and ScanObjectNN datasets,respectively.For indoor scene segmentation,training and testing were performed on the S3DIS dataset,and the mean intersection over union was 62.4%.Conclusions The proposed S^(2)ANet can effectively capture the local geometric information of point clouds,thereby improving accuracy on various tasks.

Nuclear charge radius predictions by kernel ridge regression with odd-even effects

The extended kernel ridge regression(EKRR)method with odd-even effects was adopted to improve the description of the nuclear charge radius using five commonly used nuclear models.These are:(i)the isospin-dependent A^(1∕3) formula,(ii)relativistic continuum Hartree-Bogoliubov(RCHB)theory,(iii)Hartree-Fock-Bogoliubov(HFB)model HFB25,(iv)the Weizsacker-Skyrme(WS)model WS*,and(v)HFB25*model.In the last two models,the charge radii were calculated using a five-parameter formula with the nuclear shell corrections and deformations obtained from the WS and HFB25 models,respectively.For each model,the resultant root-mean-square deviation for the 1014 nuclei with proton number Z≥8 can be significantly reduced to 0.009-0.013 fm after considering the modification with the EKRR method.The best among them was the RCHB model,with a root-mean-square deviation of 0.0092 fm.The extrapolation abilities of the KRR and EKRR methods for the neutron-rich region were examined,and it was found that after considering the odd-even effects,the extrapolation power was improved compared with that of the original KRR method.The strong odd-even staggering of nuclear charge radii of Ca and Cu isotopes and the abrupt kinks across the neutron N=126 and 82 shell closures were also calculated and could be reproduced quite well by calculations using the EKRR method.

Revisiting the Electron Radius in Light of Length Contraction

In previous work, the electron radius was identified as the “actual electron radius.” However, this is more accurately described as the electron radius at rest. This study reexamines the electron with an emphasis on the electron radius under motion, incorporating the effects of length contraction. The findings suggest that the radius is subject to Lorentz contraction, which has interesting implications for relativistic effects at the subatomic level.

Decoupling of temporal/spatial broadening effects in Doppler wind LiDAR by 2D spectral analysis

Pulse echo accumulation is commonly employed in coherent Doppler wind LiDAR(light detection and ranging)under the assumption of steady wind.Here,the measured spectral data are analyzed in the time dimension and frequency dimension to cope with the temporal wind shear and achieve the optimal accumulation time.A hardware-efficient algorithm combining the interpolation and cross-correlation is used to enhance the wind retrieval accuracy by reducing the frequency sampling interval and then reduce the spectral width calculation error.Moreover,the temporal broadening effect and spatial broadening effect are decoupled according to the strategy we developed.

Estimating wheat spike-leaf composite indicator(SLI)dynamics by coupling spectral indices and machine learning

The contribution of spike photosynthesis to grain yield(GY)has been overlooked in the accurate spectral prediction of yield.Thus,it’s essential to construct and estimate a yield-related phenotypic trait considering spike photosynthesis.Based on field and spectral reflectance data from 19 wheat cultivars under two nitrogen fertilization conditions in two years,our objectives were to(i)construct a yield-related phenotypic trait(spike–leaf composite indicator,SLI)accounting for the contribution of the spike to photosynthesis,(ii)develop a novel spectral index(enhanced triangle vegetation index,ETVI3)sensitive to SLI,and(iii)establish and evaluate SLI estimation models by integrating spectral indices and machine learning algorithms.The results showed that SLI was sensitive to nitrogen fertilizer and wheat cultivar variation as well as a better predictor of yield than the leaf area index.ETVI3 maintained a strong correlation with SLI throughout the growth stage,whereas the correlations of other spectral indices with SLI were poor after spike emergence.Integrating spectral indices and machine learning algorithms improved the estimation accuracy of SLI,with the most accurate estimates of SLI showing coefficient of determination,root mean square error(RMSE),and relative RMSE values of 0.71,0.047,and 26.93%,respectively.These results provide new insights into the role of fruiting organs for the accurate spectral prediction of GY.This high-throughput SLI estimation approach can be applied for wheat yield prediction at whole growth stages and may be assisted with agronomical practices and variety selection.

Whole-volume histogram analysis of spectral-computed tomography iodine maps characterizes HER2 expression in gastric cancer

BACKGROUND Although surgery remains the primary treatment for gastric cancer(GC),the identification of effective alternative treatments for individuals for whom surgery is unsuitable holds significance.HER2 overexpression occurs in approximately 15%-20%of advanced GC cases,directly affecting treatment-related decisions.Spectral-computed tomography(sCT)enables the quantification of material compositions,and sCT iodine concentration parameters have been demonstrated to be useful for the diagnosis of GC and prediction of its invasion depth,angioge-nesis,and response to systemic chemotherapy.No existing report describes the prediction of GC HER2 status through histogram analysis based on sCT iodine maps(IMs).AIM To investigate whether whole-volume histogram analysis of sCT IMs enables the prediction of the GC HER2 status.METHODS This study was performed with data from 101 patients with pathologically confirmed GC who underwent preoperative sCT examinations.Nineteen parameters were extracted via sCT IM histogram analysis:The minimum,maximum,mean,standard deviation,variance,coefficient of variation,skewness,kurtosis,entropy,percentiles(1st,5th,10th,25th,50th,75th,90th,95th,and 99th),and lesion volume.Spearman correlations of the parameters with the HER2 status and clinicopathological parameters were assessed.Receiver operating characteristic curves were used to evaluate the parameters’diagnostic performance.RESULTS Values for the histogram parameters of the maximum,mean,standard deviation,variance,entropy,and percentiles were significantly lower in the HER2+group than in the HER2–group(all P<0.05).The GC differentiation and Lauren classification correlated significantly with the HER2 status of tumor tissue(P=0.001 and 0.023,respectively).The 99th percentile had the largest area under the curve for GC HER2 status identification(0.740),with 76.2%,sensitivity,65.0%specificity,and 67.3%accuracy.All sCT IM histogram parameters correlated positively with the GC HER2 status(r=0.237-0.337,P=0.001-0.017).CONCLUSION Whole-lesion histogram parameters derived from sCT IM analysis,and especially the 99th percentile,can serve as imaging biomarkers of HER2 overexpression in GC.

Anisotropy of Trabecular Bone from Ultra-Distal Radius Digital X-Ray Imaging: Effects on Bone Mineral Density and Age

Background: When applied to trabecular bone X-ray images, the anisotropic properties of trabeculae located at ultra-distal radius were investigated by using the trabecular bone scores (TBS) calculated along directions parallel and perpendicular to the forearm. Methodology: Data from more than two hundred subjects were studied retrospectively. A DXA (GE Lunar Prodigy) scan of the forearm was performed on each subject to measure the bone mineral density (BMD) value at the location of ultra-distal radius, and an X-ray digital image of the same forearm was taken on the same day. The values of trabecular bone score along the direction perpendicular to the forearm, TBSx, and along the direction parallel to the forearm, TBSy, were calculated respectively. The statistics of TBSx and TBSy were calculated, and the anisotropy of the trabecular bone, which was defined as the ratio of TBSy to TBSx and changed with subjects’ BMD and age, was reported and analyzed. Results: The results show that the correlation coefficient between TBSx and TBSy was 0.72 (p BMD and age was reported. The results showed that decreased trabecular bone anisotropy was associated with deceased BMD and increased age in the subject group. Conclusions: This study shows that decreased trabecular bone anisotropy was associated with decreased BMD and increased age.

Comparison of a Spectral Bin and Two Multi-Moment Bulk Microphysics Schemes for Supercell Simulation:Investigation into Key Processes Responsible for Hydrometeor Distributions and Precipitation

There are more uncertainties with ice hydrometeor representations and related processes than liquid hydrometeors within microphysics parameterization(MP)schemes because of their complicated geometries and physical properties.Idealized supercell simulations are produced using the WRF model coupled with“full”Hebrew University spectral bin MP(HU-SBM),and NSSL and Thompson bulk MP(BMP)schemes.HU-SBM downdrafts are typically weaker than those of the NSSL and Thompson simulations,accompanied by less rain evaporation.HU-SBM produces more cloud ice(plates),graupel,and hail than the BMPs,yet precipitates less at the surface.The limiting mass bins(and subsequently,particle size)of rimed ice in HU-SBM and slower rimed ice fall speeds lead to smaller melting-level net rimed ice fluxes than those of the BMPs.Aggregation from plates in HU-SBM,together with snow–graupel collisions,leads to a greater snow contribution to rain than those of the BMPs.Replacing HU-SBM’s fall speeds using the formulations of the BMPs after aggregating the discrete bin values to mass mixing ratios and total number concentrations increases net rain and rimed ice fluxes.Still,they are smaller in magnitude than bulk rain,NSSL hail,and Thompson graupel net fluxes near the surface.Conversely,the melting-layer net rimed ice fluxes are reduced when the fall speeds for the NSSL and Thompson simulations are calculated using HU-SBM fall speed formulations after discretizing the bulk particle size distributions(PSDs)into spectral bins.The results highlight precipitation sensitivity to storm dynamics,fall speed,hydrometeor evolution governed by process rates,and MP PSD design.

A research on the effect of plasma spectrum collection device on LIBS spectral intensity

Only a small amount of spectral information is collected because the collection solid angle of the optical fiber probe and lens is very limited when collecting spectral information.To overcome this limitation,this study presents a novel method for acquiring plasma spectral information from various spatial directions.A parabolic-shaped plasma spectral collection device(PSCD)is employed to effectively collect more spectral information into the spectrometer,thereby enhancing the overall spectral intensity.The research objects in this study were soil samples containing different concentrations of heavy metals Pb,Cr,and Cd.The results indicate that the PSCD significantly enhances the spectral signal,with an enhancement rate of up to 45%.Moreover,the signal-to-noise ratio also increases by as much as 36%.Simultaneously,when compared to the absence of a device,it is found that there is no significant variation in plasma temperature when the PSCD is utilized.This observation eliminates the impact of the spatial effect caused by the PSCD on the spectral intensity.Consequently,a concentrationspectral intensity relationship curve is established under the PSCD.The results revealed that the linear fitting R^(2)for Pb,Cr,and Cd increased by 0.011,0.001,and 0.054,respectively.Additionally,the limit of detection(LOD)decreased by 0.361 ppm,0.901 ppm,and 0.602 ppm,respectively.These findings indicate that the spectral enhancement rate elevates with the increase in heavy metal concentration.Hence,the PSCD can effectively enhance the spectral intensity and reduce the detection limit of heavy metals in soil.

Review on article of preoperative prediction in chronic hepatitis B virus patients using spectral computed tomography and machine learning

This letter comments on the article that developed and tested a machine learning model that predicts lymphovascular invasion/perineural invasion status by combining clinical indications and spectral computed tomography characteristics accurately.We review the research content,methodology,conclusions,strengths and weaknesses of the study,and introduce follow-up research to this work.

Hemiarthroplasty for irreparable distal radius fractures in the elderly:A comprehensive review

BACKGROUND Elderly patients maintaining functional independence can now be candidates for primary wrist hemiarthroplasty to manage acute irreparable distal radius fractures(DRFs).However,further investigation with long-term follow-up is required to validate these initial findings.AIM To review the literature on the outcomes of distal radius hemiarthroplasty with available implants to assess its viability as a treatment option.METHODS A comprehensive review of the literature was conducted using electronic databases,including PubMed,Medline,and Scopus.The search terms employed were"distal radius fracture","hemiarthroplasty","wrist arthroplasty",and related terminology.The search was restricted to articles published in English from 1980 until June 2023.Inclusion criteria encompassed cases or case series of DRF treated with hemiarthroplasty,providing clinical or radiographic outcomes,and published in peer-reviewed journals.RESULTS A total of 2508 articles from PubMed and 883 from Scopus were identified initially.Following screening and removal of duplicates,13 articles met the inclusion criteria.These articles,predominantly clinical retrospective studies,provided insights into hemiarthroplasty outcomes,including functional improvements and complications.Hemiarthroplasty was a treatment option for complex DRF,particularly those cases with severe comminution,intraarticular involvement,or severe osteoporosis.Functional outcomes demonstrated improvements in pain relief,wrist mobility,and grip strength,with variability across studies.Complications included implant loosening,infection,nerve injury,and stiffness,with varying incidence rates influenced by surgical techniques and implant choice.Long-term outcomes were inadequately documented,warranting further research.CONCLUSION Hemiarthroplasty is a promising treatment for irreparable DRF in the elderly.Long-term outcomes and complications require further study.

System error iterative identification for underwater positioning based on spectral clustering

The observation error model of the underwater acous-tic positioning system is an important factor to influence the positioning accuracy of the underwater target.For the position inconsistency error caused by considering the underwater tar-get as a mass point,as well as the observation system error,the traditional error model best estimation trajectory(EMBET)with little observed data and too many parameters can lead to the ill-condition of the parameter model.In this paper,a multi-station fusion system error model based on the optimal polynomial con-straint is constructed,and the corresponding observation sys-tem error identification based on improved spectral clustering is designed.Firstly,the reduced parameter unified modeling for the underwater target position parameters and the system error is achieved through the polynomial optimization.Then a multi-sta-tion non-oriented graph network is established,which can address the problem of the inaccurate identification for the sys-tem errors.Moreover,the similarity matrix of the spectral cluster-ing is improved,and the iterative identification for the system errors based on the improved spectral clustering is proposed.Finally,the comprehensive measured data of long baseline lake test and sea test show that the proposed method can accu-rately identify the system errors,and moreover can improve the positioning accuracy for the underwater target positioning.

DiriNet:An Estimation Network for Spectral Response Function and Point Spread Function

Hyper-and multi-spectral image fusion is an important technology to produce hyper-spectral and hyper-resolution images,which always depends on the spectral response function andthe point spread function.However,few works have been payed on the estimation of the two degra-dation functions.To learn the two functions from image pairs to be fused,we propose a Dirichletnetwork,where both functions are properly constrained.Specifically,the spatial response function isconstrained with positivity,while the Dirichlet distribution along with a total variation is imposedon the point spread function.To the best of our knowledge,the neural network and the Dirichlet regularization are exclusively investigated,for the first time,to estimate the degradation functions.Both image degradation and fusion experiments demonstrate the effectiveness and superiority of theproposed Dirichlet network.

Efficient simulation of spatially correlated non-stationary ground motions by wavelet-packet algorithm and spectral representation method

Although the classical spectral representation method(SRM)has been widely used in the generation of spatially varying ground motions,there are still challenges in efficient simulation of the non-stationary stochastic vector process in practice.The first problem is the inherent limitation and inflexibility of the deterministic time/frequency modulation function.Another difficulty is the estimation of evolutionary power spectral density(EPSD)with quite a few samples.To tackle these problems,the wavelet packet transform(WPT)algorithm is utilized to build a time-varying spectrum of seed recording which describes the energy distribution in the time-frequency domain.The time-varying spectrum is proven to preserve the time and frequency marginal property as theoretical EPSD will do for the stationary process.For the simulation of spatially varying ground motions,the auto-EPSD for all locations is directly estimated using the time-varying spectrum of seed recording rather than matching predefined EPSD models.Then the constructed spectral matrix is incorporated in SRM to simulate spatially varying non-stationary ground motions using efficient Cholesky decomposition techniques.In addition to a good match with the target coherency model,two numerical examples indicate that the generated time histories retain the physical properties of the prescribed seed recording,including waveform,temporal/spectral non-stationarity,normalized energy buildup,and significant duration.