Search NASA⌕ Search

SEARCH · Search NASA

Results for “Spatial correction”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 181 records · Page 10

Terra and Aqua MODIS Intercomparison Using LEO-GEO Double Difference Method

The Moderate Resolution Imaging Spectroradiometer (MODIS) onboard the Terra and Aqua satellites havesuccessfully operated since their launch in 1999 and in 2002, providing more than 18 and 16 years ofcontinuous global observations, respectively. The inter-comparison between the two MODIS instruments can bevery supportive for the instrument calibration and uncertainty assessment. Aqua and Terra MODIS have almostidentical relative spectral response, spatial resolution, and dynamic range for each band. Therefore, a sitedependent correction for a sensor spectral band pair is not necessary for their comparison. However, Terra is inthe morning orbit with an equator crossing time of 10:30 am, and Aqua is in the afternoon orbit with equatorcrossing time of 1:30 pm. Consequently, there is a dearth of simultaneous nadir overpasses (SNOs)between the two satellites. Major challenges in cross-sensor comparison of instruments on different satellitesinclude differences in observation time, solar angle, and view angle over selected pseudo-invariant sites.In this work, the inter-comparisons of thermal emissive bands are performed over a pseudo-invariant target,using the observations from a sensor onboard a geostationary satellite as a bridge. Himawari8 was launched onOctober 7, 2014. The Advanced Himawari Imager (AHI) onboard Himawari8 can be used as a reference tobridge the comparison between Terra and Aqua MODIS. AHI has 16 channels; with spatial resolutions from 0.5km to 2 km at nadir and produces a full disk observations every 10 minutes. The band spectral coveragematchup, comparable spatial resolution and near-simultaneous observation between MODIS and AHI providefeasibility to implement a double difference method. This comparison method minimizes the impact of thedifference in observation time and solar angle. The comparison results will be used as an assessment for MODISinstrument calibration and will be helpful for future enhancement of the L1B product.

Himawari8↗

Wavefront Control and Image Restoration with Less Computing

PseudoDiversity is a method of recovering the wavefront in a sparse- or segmented- aperture optical system typified by an interferometer or a telescope equipped with an adaptive primary mirror consisting of controllably slightly moveable segments. (PseudoDiversity should not be confused with a radio-antenna-arraying method called pseudodiversity.) As in the cases of other wavefront- recovery methods, the streams of wavefront data generated by means of PseudoDiversity are used as feedback signals for controlling electromechanical actuators of the various segments so as to correct wavefront errors and thereby, for example, obtain a clearer, steadier image of a distant object in the presence of atmospheric turbulence. There are numerous potential applications in astronomy, remote sensing from aircraft and spacecraft, targeting missiles, sighting military targets, and medical imaging (including microscopy) through such intervening media as cells or water. In comparison with prior wavefront-recovery methods used in adaptive optics, PseudoDiversity involves considerably simpler equipment and procedures and less computation. For PseudoDiversity, there is no need to install separate metrological equipment or to use any optomechanical components beyond those that are already parts of the optical system to which the method is applied. In Pseudo- Diversity, the actuators of a subset of the segments or subapertures are driven to make the segments dither in the piston, tilt, and tip degrees of freedom. Each aperture is dithered at a unique frequency at an amplitude of a half wavelength of light. During the dithering, images on the focal plane are detected and digitized at a rate of at least four samples per dither period. In the processing of the image samples, the use of different dither frequencies makes it possible to determine the separate effects of the various dithered segments or apertures. The digitized image-detector outputs are processed in the spatial-frequency (Fourier-transform) domain to obtain measures of the piston, tip, and tilt errors over each segment or subaperture. Once these measures are known, they are fed back to the actuators to correct the errors. In addition, measures of errors that remain after correction by use of the actuators are further utilized in an algorithm in which the image is phase-corrected in the spatial-frequency domain and then transformed back to the spatial domain at each time step and summed with the images from all previous time steps to obtain a final image having a greater signal-to-noise ratio (and, hence, a visual quality) higher than would otherwise be attainable.

Lyon, Richard G.↗

Statistical Modeling of Landscape Pattern and the Effects of Coarse Spatial Resolution

Analysis of classified satellite imagery was conducted to characterize errors in estimates of area based on coarse resolution satellite imagery which are due to distortions in sizes of small fragments, and to explore the feasibility of correcting for these errors using a statistical modeling approach. Sizes of bodies of open water on European Remote-Sensing (ERS-1) synthetic aperture radar (SAR) and fire scars on Landsat Multispectral Scanner (MSS) imagery were measured using geographic information system (GIS) software. Statistical analysis of the smaller scars and ponds as observed with this imagery of relatively fine resolution demonstrated that the distribution of the sizes could be modeled by either of two types of statistical distributions-a power distribution related to fractal processes or a simple exponential distribution. A comparison of the size distribution of small burn scars as observed with Landsat with the distribution observed with the much coarser Advanced Very High Resolution Radiometer (AVHRR) imagery revealed effects of coarse resolution, each of which were associated with a particular size range, and whose net effect accounted for the bias in total burn area due to coarse spatial resolution. Size dependent area correction procedures are therefore feasible. A statistical modeling approach to area estimation, based models of the underlying and observed size distributions of fragments in the scene, may lead to a procedure which requires little or no ancillary data.

Hlavka, C. A.↗

Vorticity and Vertical Motions Diagnosed from Satellite Deep-Layer Temperatures

Spatial fields of satellite-measured deep-layer temperatures are examined in the context of quasigeostrophic theory. It is found that midtropospheric geostrophic vorticity and quasigeostrophic vertical motions can be diagnosed from microwave temperature measurements of only two deep layers. The lower- ( 1000-400 hPa) and upper- (400-50 hPa) layer temperatures are estimated from limb-corrected TIROS-N Microwave Sounding Units (MSU) channel 2 and 3 data, spatial fields of which can be used to estimate the midtropospheric thermal wind and geostrophic vorticity fields. Together with Trenberth's simplification of the quasigeostrophic omega equation, these two quantities can be then used to estimate the geostrophic vorticity advection by the thermal wind, which is related to the quasigeostrophic vertical velocity in the midtroposphere. Critical to the technique is the observation that geostrophic vorticity fields calculated from the channel 3 temperature features are very similar to those calculated from traditional, 'bottom-up' integrated height fields from radiosonde data. This suggests a lack of cyclone-scale height features near the top of the channel 3 weighting function, making the channel 3 cyclone-scale 'thickness' features approximately the same as height features near the bottom of the weighting function. Thus, the MSU data provide observational validation of the LID (level of insignificant dynamics) assumption of Hirshberg and Fritsch.

Spencer, Roy W.↗

Cloud-to-Ground Lightning Estimates Derived from SSMI Microwave Remote Sensing and NLDN

Lightning observations are collected using ground-based and satellite-based sensors. The National Lightning Detection Network (NLDN) in the United States uses multiple ground sensors to triangulate the electromagnetic signals created when lightning strikes the Earth's surface. Satellite-based lightning observations have been made from 1998 to present using the Lightning Imaging Sensor (LIS) on the NASA Tropical Rainfall Measuring Mission (TRMM) satellite, and from 1995 to 2000 using the Optical Transient Detector (OTD) on the Microlab-1 satellite. Both LIS and OTD are staring imagers that detect lightning as momentary changes in an optical scene. Passive microwave remote sensing (85 and 37 GHz brightness temperatures) from the TRMM Microwave Imager (TMI) has also been used to quantify characteristics of thunderstorms related to lightning. Each lightning detection system has fundamental limitations. TRMM satellite coverage is limited to the tropics and subtropics between 38 deg N and 38 deg S, so lightning at the higher latitudes of the northern and southern hemispheres is not observed. The detection efficiency of NLDN sensors exceeds 95%, but the sensors are only located in the USA. Even if data from other ground-based lightning sensors (World Wide Lightning Location Network, the European Cooperation for Lightning Detection, and Canadian Lightning Detection Network) were combined with TRMM and NLDN, there would be enormous spatial gaps in present-day coverage of lightning. In addition, a globally-complete time history of observed lightning activity is currently not available either, with network coverage and detection efficiencies varying through the years. Previous research using the TRMM LIS and Microwave Imager (TMI) showed that there is a statistically significant correlation between lightning flash rates and passive microwave brightness temperatures. The physical basis for this correlation emerges because lightning in a thunderstorm occurs where ice is first present in the cloud and electric charge separation occurs. These ice particles efficiently scatter the microwave radiation at the 85 and 37 GHz frequencies, thus leading to large brightness temperature depressions. Lightning flash rate is related to the total amount of ice passing through the convective updraft regions of thunderstorms. Confirmation of this relationship using TRMM LIS and TMI data, however, remains constrained to TRMM observational limits of the tropics and subtropics. Satellites from the Defense Meteorology Satellite Program (DMSP) have global coverage and are equipped with passive microwave imagers that, like TMI, observe brightness temperatures at 85 and 37 GHz. Unlike the TRMM satellite, however, DMSP satellites do not have a lightning sensor, and the DMSP microwave data has never been used to derive global lightning. In this presentation, a relationship between DMSP Special Sensor Microwave Imager (SSMI) data and ground-based cloud-to-ground (CG) lightning data from NLDN is investigated to derive a spatially complete time history of CG lightning for the USA study area. This relationship is analogous to the established using TRMM LIS and TMI data. NLDN has the most spatially and temporally complete CG lightning data for the USA, and therefore provides the best opportunity to find geospatially coincident observations with SSMI sensors. The strongest thunderstorms generally have minimum 85 GHz Polarized Corrected brightness Temperatures (PCT) less than 150 K. Archived radar data was used to resolve the spatial extent of the individual storms. NLDN data for that storm spatial extent defined by radar data was used to calculate the CG flash rate for the storm. Similar to results using TRMM sensors, a linear model best explained the relationship between storm-specific CG flash rates and minimum 85 GHz PCT. However, the results in this study apply only to CG lightning. To extend the results to weaker storms, the probability of CG lightning (instead of the flash rate) was calculated for storms having 85 GHz PCT greater than 150 K. NLDN data was used to determine if a CG strike occurred for a storm. This probability of CG lightning was plotted as a function of minimum 85 GHz PCT and minimum 37 GHz PCT. These probabilities were used in conjunction with the linear model to estimate the CG flash rate for weaker storms with minimum 85 GHz PCTs greater than 150 K. Results from the investigation of CG lightning and passive microwave radiation signals agree with the previous research investigating total lightning and brightness temperature. Future work will take the established relationships and apply them to the decades of available DMSP data for the USA to derive a map of CG lightning flash rates. Validation of this method and uncertainty analysis will be done by comparing the derived maps of CG lightning flash rates against existing NLDN maps of CG lightning flash rates.

Winesett, Thomas↗

Optical tomography for flow visualization of the density field around a revolving helicopter rotor blade

In this paper, a tomographic procedure for reconstructing the density field around a helicopter rotor blade tip from remote optical line-of-sight measurements is discussed. Numerical model studies have been carried out to investigate the influence of the number of available views, limited width viewing, and ray bending on the reconstruction. Performance is measured in terms of the mean-square error. It is found that very good reconstructions can be obtained using only a small number of views even when the width of view is smaller than the spatial extent of the object. An iterative procedure is used to correct for ray bending due to refraction associated with the sharp density gradients (shocks).

Snyder, R.↗

Optical communication for space missions

Activities performed at NASA/GSFC (Goddard Space Flight Center) related to direct detection optical communications for space applications are discussed. The following subject areas are covered: (1) requirements for optical communication systems (data rates and channel quality; spatial acquisition; fine tracking and pointing; and transmit point-ahead correction); (2) component testing and development (laser diodes performance characterization and life testing; and laser diode power combining); (3) system development and simulations (The GSFC pointing, acquisition and tracking system; hardware description; preliminary performance analysis; and high data rate transmitter/receiver systems); and (4) proposed flight demonstration of optical communications.

Firtmaurice, M.↗

On sound generation by turbulent convection: A new look at old results

We have revisited the problem of acoustic wave generation by turbulent convection in stellar atmospheres. The theory of aerodynamically generated sound, originally developed by Lighthill and later modified by Stein to include the effects of stratification, has been used to estimate the acoustic wave energy flux generated in solar and stellar convection zones. We correct the earlier computations by incorporating an improved description of the spatial and temporal spectrum of the turbulent convection. We show the dependence of the resulting wave fluxes on the nature of the turbulence, and compute the wave energy spectra and wave energy fluxes generated in the Sun on the basis of a mixing-length model of the solar convection zone. In contrast to the previous results, we show that the acoustic energy generation does not depend very sensitively on the turbulent energy spectrum. However, typical total acoustic fluxes of order F(sub A) = 5 x 10(exp 7) ergs/sq cm/s with a peak of the acoustic frequency spectrum near omega = 100 mHz are found to be comparable to those previously calculated. The acoustic flux turns out to be strongly dependent on the solar model, scaling with the mixing-length parameter alpha as alpha(exp 3.8). The computed fluxes most likely constitute a lower limit on the acoustic energy produced in the solar convection zone if recent convection simulations suggesting the presence of shocks near the upper layers of the convection zone apply to the Sun.

Musielak, Z. E.↗

Applications of High Etendue Line-Profile Spectro-Polarimetry to the Study of the Atmospheric and Magnetospheric Environments of the Jovian Icy Moons

Electrodynamic effects play a significant, global role in the state and energization of the Earth's ionosphere/magnetosphere, but even more so on Jupiter, where the auroral energy input is four orders of magnitude greater than on Earth. The Jovian magnetosphere is distinguished from Earth's by its rapid rotation rate and contributions from satellite atmospheres and internal plasma sources. The electrodynamic effects of these factors have a key role in the state and energization of the ionosphere-corona- plasmasphere system of the planet and its interaction with Io and the icy satellites. Several large scale interacting processes determine conditions near the icy moons beginning with their tenuous atmospheres produced from sputtering, evaporative, and tectonic/volcanic sources, extending out to exospheres that merge with ions and neutrals in the Jovian magnetosphere. This dynamic environment is dependent on a complex network of magnetospheric currents that act on global scales. Field aligned currents connect the satellites and the middle and tail magnetospheric regions to the Jupiter's poles via flux tubes that produce as bright auroral and satellite footprint emissions in the upper atmosphere. This large scale transfer of mass, momentum, and energy (e.g. waves, currents) means that a combination of complementary diagnostics of the plasma, neutral, and and field network must be obtained near simultaneously to correctly interpret the results. This presentation discusses the applicability of UV spatial heterodyne spectroscopy (SHS) to the broad study of this system on scales from satellite surfaces to Jupiter's aurora and corona.

Harris, Walter M.↗

Pulmonary blood flow distribution in sheep: effects of anesthesia, mechanical ventilation, and change in posture

BACKGROUND: Recent studies providing high-resolution images of pulmonary perfusion have questioned the classical zone model of pulmonary perfusion. Hence the present work was undertaken to provide detailed maps of regional pulmonary perfusion to examine the influence of anesthesia, mechanical ventilation, and posture. METHODS: Pulmonary perfusion was analyzed with intravenous fluorescent microspheres (15 microm) in six sheep studied in four conditions: prone and awake, prone with pentobarbital-anesthesia and breathing spontaneously, prone with anesthesia and mechanical ventilation, and supine with anesthesia and mechanical ventilation. Lungs were air dried at total lung capacity and sectioned into approximately 1,100 pieces (about 2 cm3) per animal. The pieces were weighed and assigned spatial coordinates. Fluorescence was read on a spectrophotometer, and signals were corrected for piece weight and normalized to mean flow. Pulmonary blood flow heterogeneity was assessed using the coefficient of variation of flow data. RESULTS: Pentobarbital anesthesia and mechanical ventilation did not influence perfusion heterogeneity, but heterogeneity increased when the animals were in the supine posture (P < 0.01). Gravitational flow gradients were absent in the prone position but present in the supine (P < 0.001 compared with zero). Pulmonary perfusion was distributed with a hilar-to-peripheral gradient in animals breathing spontaneously (P < 0.05). CONCLUSIONS: The influence of pentobarbital anesthesia and mechanical ventilation on pulmonary perfusion heterogeneity is small compared with the effect of changes in posture. Analysis of flow gradients indicate that gravity plays a small role in determining pulmonary blood flow distribution.

NASA Discipline Cardiopulmonary↗

Daytime adaptive optics for deep space optical communications

The deep space optical communications subsystem offers a higher bandwidth communications link in smaller size, lower mass, and lower power consumption subsystem than does RF. To demonstrate the benefit of this technology to deep space communications NASA plans to launch an optical telecommunications package on the 2009 Mars Telecommunications orbiter spacecraft. Current performance goals are 30-Mbps from opposition, and 1-Mbps near conjunction (-3 degrees Sun-Earth-Probe angle). Yet, near conjunction the background noise from the day sky will degrade the performance of the optical link. Spectral and spatial filtering and higher modulation formats can mitigate the effects of background sky. Narrowband spectral filters can result in loss of link margin, and higher modulation formats require higher transmitted peak powers. In contrast, spatial filtering at the receiver has the potential of being lossless while providing the required sky background rejection. Adaptive optics techniques can correct wave front aberrations caused by atmospheric turbulence and enable near-diffraction-limited performance of the receiving telescope. Such performance facilitates spatial filtering, and allows the receiver field-of-view and hence the noise from the sky background to be reduced.

adaptive optics↗

The Meteoroid Fluence at Mars Due to Comet C/2013 A1 (Siding Spring)

Long-period comet C/2013 A1 (Siding Spring) will experience a close encounter with Mars on 2014 Oct 19. A collision between the comet and the planet has been ruled out, but the comet's coma may envelop Mars and its man-made satellites. By the time of the close encounter, five operational spacecraft will be present near Mars. Characterizing the coma is crucial for assessing the risk posed to these satellites by meteoroid impacts. We present an analytic model of cometary comae that describes the spatial and size distributions of cometary dust and meteoroids. This model correctly reproduces, to within an order of magnitude, the number of impacts recorded by Giotto near 1P/Halley [1] and by Stardust near comet 81P/Wild 2 [2]. Applied to Siding Spring, our model predicts a total particle fluence near Mars of 0.02 particles per square meter. In order to determine the degree to which Siding Spring's coma deviates from a sphere, we perform numerical simulations which take into account both gravitational effects and radiative forces. We take the entire dust component of the coma and tail continuum into account by simulating the ejection and evolution of dust particles from comet Siding Spring. The total number of particles simulated is essentially a free parameter and does not provide a check on the total fluence. Instead, these simulations illustrate the degree to which the coma of Siding Spring deviates from the perfect sphere described by our analytic model (see Figure). We conclude that our analytic model sacrifices less than an order of magnitude in accuracy by neglecting particle dynamics and radiation pressure and is thus adequate for order-of-magnitude fluence estimates. Comet properties may change unpredictably and therefore an analytic coma model that enables quick recalculation of the meteoroid fluence is highly desirable. NASA's Meteoroid Environment Office is monitoring comet Siding Spring and taking measurements of cometary brightness and dust production. We will discuss our coma model and nominal fluence taking the latest observations into account.

Moorhead, A.↗

Update on EPIC MAIAC V2 Algorithm: Atmospheric Correction

Alexei Lyapustin presented a new v2 MAIAC atmospheric correction algorithm for EPIC. The v2 algorithm introduced new projection removing spatial distortions, a more accurate retrieval of spectral surface reflectance using ancillary MAIAC MODIS BRDF dataset, and a new technique for simultaneous retrieval of aerosol optical depth and spectral absorption. AERONET validation shows that v2 significantly improved AOD accuracy compared to v1, while the retrieved single scattering albedo in the Blue shows a very good accuracy comparable to that of AERONET for both biomass burning smoke and mineral dust aerosols.

Alexei Lyapustin↗

Rapid Glacier Retreat Rates Observed in West Antarctica

The Pope, Smith, and Kohler glaciers, in the Amundsen Sea Embayment of West Antarctica, have experienced enhanced ocean-induced ice-shelf melt, glacier acceleration, ice thinning, and grounding line retreat in the past thirty years. Here we present observations of the grounding line retreat of these glaciers since 2014 using a constellation of interferometric radar satellites combined with precision surface elevation data. We find that the grounding lines develop spatially-variable, kilometre-scale, tidally-induced migration zones. After correction for tidal effects, we detect a sustained pattern of retreat coincident with high melt rates of un-grounded ice, marked by episodes of more rapid retreat. In 2017, Pope Glacier retreated 3.5 km in 3.6 months, or 11.7 km/yr. In 2016-2018, Smith West retreated at 2 km/yr and Kohler at 1.3 km/yr. While the retreat slowed down in 2018-2020, these retreat rates are faster than anticipated by numerical models on yearly time scales. We hypothesize that the rapid retreat is caused by un-represented, vigorous ice-ocean interactions acting within newly-formed cavities at the ice-ocean boundary.

Climate change↗

Rotational wind indicator enhances control of rotated displays

Rotation by 108 deg of the spatial mapping between a visual display and a manual input device produces large spatial errors in a discrete aiming task. These errors are not easily corrected by voluntary mental effort, but the central nervous system does adapt gradually to the new mapping. Bernotat (1970) showed that adding true hand position to a 90 deg rotated display improved performance of a compensatory tracking task, but tracking error rose again upon removal of the explicit cue. This suggests that the explicit error signal did not induce changes in the neural mapping, but rather allowed the operator to reduce tracking error using a higher mental strategy. In this report, we describe an explicit visual display enhancement applied to a 108 deg rotated discrete aiming task. A 'wind indicator' corresponding to the effect of the mapping rotation is displayed on the operator-controlled cursor. The human operator is instructed to oppose the virtual force represented by the indicator, as one would do if flying an airplane in a crosswind. This enhancement reduces spatial aiming error in the first 10 minutes of practice by an average of 70 percent when compared to a no enhancement control condition. Moreover, it produces adaptation aftereffect, which is evidence of learning by neural adaptation rather than by mental strategy. Finally, aiming error does not rise upon removal of the explicit cue.

Cunningham, H. A.↗

Very High Resolution, Altitude-Corrected, TMPA-Based Monthly Satellite Precipitation Product Over the CONUS

The Tropical Rainfall Measuring Mission (TRMM) Multisatellite Precipitation Analysis (TMPA) product provided over 17 years of gridded precipitation datasets. However, the accuracy and spatial resolution of TMPA limits the applicability in hydrometeorological applications. We present a dataset that enhances the accuracy and spatial resolution of the TMPA monthly product (3B43). We resample the TMPA data to a 1 km grid and apply a correction function derived from the Parameter-elevation Regressions on Independent Slopes Model (PRISM) to reduce bias in the data. We confirm a linear relationship between bias and elevation above 1,500 meters where TMPA underestimates measured precipitation, providing a proof-of-concept of how simple linear scaling can be used to augment existing satellite datasets. The result of the correction is the High-Resolution Altitude-Corrected Precipitation product (HRAC-Precip) for the CONUS. Using 9,200 precipitation stations from the Global Historical Climatology Network (GHCN), we compare the accuracy of TMPA 3B43 versus the new HRAC-Precip product. The results show an improvement of the mean absolute error of 12.98% on average.

Hossein Hashemi↗

Temporal and spatial variability of aerosol optical depth in the Sahel region in relation to vegetation remote sensing

In order to monitor the aerosol characteristics needed for atmospheric correction of remotely sensed data, a network of sun photometers was established in the Sahel region of Senegal, Mali, and Niger. Data analysis suggests that there is a high spatial variability of the aerosol optical thickness tau(a) in the western Sahel region. At a 67 percent confidence level the instantaneous values of tau(a) can be extrapolated approximately 270-400 km with an error tolerance of 50 percent. Spatial variability in the dry season is found to be of a similar magnitude. The ranges of variations in the NDVI in the Sahel region are shown to be approximately 0.02 and 0.01, respectively, due to commonly observed fluctuations in the aerosol optical thickness and aerosol size distribution.

Holben, B. N.↗

Aperiodicity Correction for Rotor Tip Vortex Measurements

The initial roll-up of a tip vortex trailing from a model-scale, hovering rotor was measured using particle image velocimetry. The unique feature of the measurements was that a microscope was attached to the camera to allow much higher spatial resolution than hitherto possible. This also posed some unique challenges. In particular, the existing methodologies to correct for aperiodicity in the tip vortex locations could not be easily extended to the present measurements. The difficulty stemmed from the inability to accurately determine the vortex center, which is a prerequisite for the correction procedure. A new method is proposed for determining the vortex center, as well as the vortex core properties, using a least-squares fit approach. This approach has the obvious advantage that the properties are derived from not just a few points near the vortex core, but from a much larger area of flow measurements. Results clearly demonstrate the advantage in the form of reduced variation in the estimated core properties, and also the self-consistent results obtained using three different aperiodicity correction methods.

Ramasamy, Manikandan↗