Development of digital interactive processing system for NOAA satellites AVHRR data

The paper discusses the digital image processing system for NOAA/AVHRR data including Land applications — configured around VAX 11/750 host computer supported with FPS 100 Array Processor, Comtal graphic display and HP Plotting devices; wherein the system software for relational Data Base together with query and editing facilities, Man-Machine Interface using form, menu and prompt inputs including validation of user entries for data type and range; preprocessing software for data calibration, Sun-angle correction, Geometric Corrections for Earth curvature effect and Earth rotation offsets and Earth location of AVHRR image have been accomplished. The implemented image enhancement techniques such as grey level stretching, histogram equalization and convolution are discussed. The software implementation details for the computation of vegetative index and normalized vegetative index using NOAA/AVHRR channels 1 and 2 data together with output are presented; scientific background for such computations and obtainability of similar indices from Landsat/MSS data are also included. The paper concludes by specifying the further software developments planned and the progress envisaged in the field of vegetation index studies.     

HS—400半自动周节检查仪的数据处理系统

ＨＳ—４００半自动周节检查仪数据处理系统是在其主机基础上研制成功的。该系统采用相对测量法，对齿轮周节参数进行测量，可直接计算出ΔＦ＿ｐ和Δｆ＿ｐ，并可打印输出测量结果和曲线。该系统的自动控制系统和安全保护系统采用了ＧＡＬ技术，使自动控制和安全保护更可靠。该系统研制成功后提高了原仪器的准确度和自动化程度。与德国的ＵＰ—６３０进行比对，其测量值之差小于２μｍ。     

Status of the T2L2/Jason2 Experiment

The T2L2 (Time Transfer by Laser Link) project, developed by CNES and OCA will permit the synchronization of remote ultra stable clocks and the determination of their performances over intercontinental distances. The principle of the experiment derives from Satellite Laser Ranging (SLR) technology with dedicated space equipment. T2L2 was accepted in 2005 to be on board the Jason2 altimetry satellite. The payload consists of both event timer and photo detection modules. The system uses the ultra-stable quartz oscillator of DORIS as on-board reference clock on one hand, and the Laser Reflector Array, making T2L2 a real two-way time transfer system on the other hand. The expected time stability of the T2L2 instrument (detection and timing), referenced by the DORIS oscillator and including all internal error sources should be at the level of 10–12 ps at 1 s and <1 ps at 1000 s. The metrological specifications of T2L2 should permit to maintain a precision of 1 to a few ps when measuring the phase of a clock during around 1000 seconds.     

Correlation of IRTAM and FPMU data confirming the application of IRTAM to support ISS Program safety

A “Real-Time” plasma hazard assessment process was developed to support International Space Station (ISS) Program real-time decision-making providing solar array constraint relief information for Extravehicular Activities (EVAs) planning and operations. This process incorporates real-time ionospheric conditions, ISS solar arrays’ orientation, ISS flight attitude, and where the EVA will be performed on the ISS. This assessment requires real-time data that is presently provided by the Floating Potential Measurement Unit (FPMU) which measures the ISS floating potential (FP), along with ionospheric electron number density (Ne) and electron temperature (Te), in order to determine the present ISS environment. Once the present environment conditions are correlated with International Reference Ionosphere (IRI) values, IRI is used to forecast what the environment could become in the event of a severe geomagnetic storm. If the FPMU should fail, the Space Environments team needs another source of data which is utilized to support a short-term forecast for EVAs. The IRI Real-Time Assimilative Mapping (IRTAM) model is an ionospheric model that uses real-time measurements from a large network of digisondes to produce foF2 and hmF2 global maps in 15?min cadence. The Boeing Space Environments team has used the IRI coefficients produced in IRTAM to calculate the Ne along the ISS orbital track. The results of the IRTAM model have been compared to FPMU measurements and show excellent agreement. IRTAM has been identified as the FPMU back-up system that will be used to support the ISS Program if the FPMU should fail.     

Fractal signatures for multiscale processing of hyperspectral image data

A multiscale approach to hyperspectral image data analysis using fractal signatures was proposed and implemented in the Interactive Data Language (IDL). For 2-D hyperspectral curves, fractal signature measures the changes in curve length with changing scale. Using NASA’s Earth Observing-1 (EO-1) Hyperion image from a study area near Denton, Texas, USA, the capabilities of fractal signatures in discriminating different land cover types were presented in three different ways: (1) fractal signature curves, (2) distances between fractal signatures, and (3) fractal signature images. The asymmetry in length measurement was found to be effective in handling hyperspectral curves obtained from Hyperion radiance data. The contribution of fractal signature images was shown through comparison of image classification results. The results from the Hyperion radiance data suggest that fractal signatures at certain scales can reveal important differences in land cover types.     

Estimation of the Love and Shida numbers: h2, l2 using SLR data for the low satellites

In this paper we present results for the global elastic parameters: Love number h₂ and Shida number l₂ derived from the analysis of Satellite Laser Ranging (SLR) data. SLR data for the two low satellites STELLA (H = 800 km) and STARLETTE (H = 810 km) observed during 2.5 years from January 3, 2005 until July 1, 2007 with 18 globally distributed ground stations were analyzed. The analysis was done separately for the two satellites. We do a sequential analysis and study the stability and convergence of the estimates as a function of length of the data set used.     

The X-TRACK/ALES multi-mission processing system: New advances in altimetry towards the coast

In the context of the ESA Climate Change Initiative project, a new coastal sea level altimetry product has been developed in order to support advances in coastal sea level variability studies. Measurements from Jason-1,2&3 missions have been retracked with the Adaptive Leading Edge Subwaveform (ALES) Retracker and then ingested in the X-TRACK software with the best possible set of altimetry corrections. These two coastal altimetry processing approaches, previously successfully validated and applied to coastal sea level research, are combined here for the first time in order to derive a 16-year-long (June 2002 to May 2018), high-resolution (20-Hz), along-track sea level dataset in six regions: Northeast Atlantic, Mediterranean Sea, West Africa, North Indian Ocean, Southeast Asia and Australia. The study demonstrates that this new coastal sea level product called X-TRACK/ALES is able to extend the spatial coverage of sea level altimetry data ~3.5 km in the land direction, when compared to the X-TRACK 1-Hz dataset. We also observe a large improvement in coastal sea level data availability from Jason-1 to Jason-3, with data at 3.6 km, 1.9 km and 0.9 km to the coast on average, for Jason-1, Jason-2 and Jason-3, respectively. When combining measurements from Jason-1 to Jason-3, we reach a distance of 1.2–4 km to the coast. When compared to tide gauge data, the accuracy of the new altimetry near-shore sea level estimations also improves. In terms of correlations with a large set of independent tide gauge observations selected in the six regions, we obtain an average value of 0.77. We also show that it is now possible to derive from the X-TRACK/ALES product an estimation of the ocean current variability up to 5 km to the coast. This new altimetry dataset, freely available, will provide a valuable contribution of altimetry in coastal marine research community.     

Upgrading CCIR's foF2 maps using available ionosondes and genetic algorithms

We have developed a new approach towards a new database of the ionospheric parameter $f_{o}F2$. This parameter, being the frequency of the maximum of the ionospheric electronic density profile and its main modeller, is of great interest not only in atmospheric studies but also in the realm of radio propagation. The current databases, generated by CCIR (Committee Consultative for Ionospheric Radiowave propagation) and URSI (International Union of Radio Science), and used by the IRI (International Reference Ionosphere) model, are based on Fourier expansions and have been built in the 60s from the available ionosondes at that time. The main goal of this work is to upgrade the databases by using new available ionosonde data. To this end we used the IRI diurnal/spherical expansions to represent the $f_{o}F2$ variability, and computed its coefficients by means of a genetic algorithm (GA). In order to test the performance of the proposed methodology, we applied it to the South American region with data obtained by RAPEAS (Red Argentina para el Estudio de la Atmósfera Superior, i.e. Argentine Network for the Study of the Upper Atmosphere) during the years 1958–2009. The new GA coefficients provide a global better fit of the IRI model to the observed $f_{o}F2$ than the CCIR coefficients. Since the same formulae and the same number of coefficients were used, the overall integrity of IRI’s typical ionospheric feature representation was preserved. The best improvements with respect to CCIR are obtained at low solar activities, at large (in absolute value) modip latitudes, and at night-time. The new method is flexible in the sense that can be applied either globally or regionally. It is also very easy to recompute the coefficients when new data is available. The computation of a third set of coefficients corresponding to days of medium solar activity in order to avoid the interpolation between low and high activities is suggested. The same procedure as for $f_{o}F2$ can be perfomed to obtain the ionospheric parameter M(3000)F2.     

The accuracy of data from ionosondes for the estimation of hmF2 and the identification of global change in the ionosphere

A long temporal series of simulated ionograms was generated with a superimposed secular variation of −14 km/century on the h_mF2 parameter. These ionograms were interpreted by the automatic scaling program Autoscala. By applying four different empirical formulas, four artificial series of h_mF2 were generated and then processed with the same methods used by other authors for real data sets. Data analysis of the simulated ionograms revealed the artificially imposed long-term trend. These results lead to the conclusion, that regardless of the empirical formula used, the accuracy of h_mF2 from ionosonde measurements would be adequate to observe a long-term trend of −14 km/century.     

Open framework for objective evaluation of crater detection algorithms with first test-field subsystem based on MOLA data

Crater Detection Algorithms (CDAs) applications range from estimation of lunar/planetary surface age to autonomous landing on planets and asteroids and advanced statistical analyses. A large amount of work on CDAs has already been published. However, problems arise when evaluation results of some new CDA have to be compared with already published evaluation results. The problem is that different authors use different test-fields, different Ground-Truth (GT) catalogues, and even different methodologies for evaluation of their CDAs. Re-implementation of already published CDAs or its evaluation environment is a time-consuming and unpractical solution to this problem. In addition, implementation details are often insufficiently described in publications. As a result, there is a need in research community to develop a framework for objective evaluation of CDAs. A scientific question is how CDAs should be evaluated so that the results are easily and reliably comparable. In attempt to solve this issue we first analyzed previously published work on CDAs. In this paper, we propose a framework for solution of the problem of objective CDA evaluation. The framework includes: (1) a definition of the measure for differences between craters; (2) test-field topography based on the 1/64° MOLA data; (3) the GT catalogue wherein each of 17,582 craters is aligned with MOLA data and confirmed with catalogues by N.G. Barlow et al. and J.F. Rodionova et al.; (4) selection of methodology for training and testing; and (5) a Free-response Receiver Operating Characteristics (F-ROC) curves as a way to measure CDA performance. The handling of possible improvements of the framework in the future is additionally addressed as a part of discussion of results. Possible extensions with additional test-field subsystems based on visual images, data sets for other planets, evaluation methodologies for CDAs developed for different purposes than cataloguing of craters, are proposed as well. The goal of the proposed framework is to contribute to the research community by establishing guidelines for objective evaluation of CDAs.     

径向圆跳动误差的最小二乘评定数学模型和微机数据处理

用传统的测量方法只能得到被测零件的径向圆跳动误差的近似值。为了得到该项误差的准确值，需要研究新的测量方法。建立了径向圆跳动误差最小二乘评定数学模型，并编制了高级语言数据处理程序。在万能工具显微镜上获得了采样数据，给出了计算结果。     

Evaluating the statistical performance of less applied algorithms in classification of worldview-3 imagery data in an urbanized landscape

In recent decade, analyzing the remotely sensed imagery is considered as one of the most common and widely used procedures in the environmental studies. In this case, supervised image classification techniques play a central role. Hence, taking a high resolution Worldview-3 over a mixed urbanized landscape in Iran, three less applied image classification methods including Bagged CART, Stochastic gradient boosting model and Neural network with feature extraction were tested and compared with two prevalent methods: random forest and support vector machine with linear kernel. To do so, each method was run ten time and three validation techniques was used to estimate the accuracy statistics consist of cross validation, independent validation and validation with total of train data. Moreover, using ANOVA and Tukey test, statistical difference significance between the classification methods was significantly surveyed. In general, the results showed that random forest with marginal difference compared to Bagged CART and stochastic gradient boosting model is the best performing method whilst based on independent validation there was no significant difference between the performances of classification methods. It should be finally noted that neural network with feature extraction and linear support vector machine had better processing speed than other.     

Satellite ocean colour algorithm for Prochlorococcus, Synechococcus, and picoeukaryotes concentration retrieval in the South China Sea

An algorithm for retrieval of surface waters cell concentrations (in cell/ml) for three picophytoplankton components, Prochlorococcus (Pro), Synechococcus (Syn), and picoeukaryotes (Peuk) in the South China Sea (SCS), from ocean colour satellite data was developed and tested. Level 3 merged multisensor Ocean Colour Climate Change Initiative satellite data is used. Training is performed using in situ data on abundances of the three phytoplankton components. Several predictors derived from satellite reflectance data were tested. The regression form that assures the highest accuracy of the algorithm was chosen based on cross-validation (CV). According to the CV on test data subset, the algorithm performance is characterized by the r value 0.89, 0.72, and 0.73 and MAPD 38, 71 and 51% for Peuk, Pro, and Syn respectively. This is one of the few studies aimed at the Peuk, Pro, and Syn distribution research in the northern SCS using ocean colour satellite data. This is the only research providing algorithm with accuracy estimates of the Peuk, Pro, and Syn concentrations retrieval from the ocean colour data. Analysis of the developed algorithm allows us to conclude that both mechanisms (specific spectral features caused by pigments composition and spectrum features sensitive to general primary productivity, e.g. band ratios in 443–510?nm range and spectrum absolute values) are important for getting accurate information on the picophytoplankton composition.     

MultiWaveLink: An interactive data base for the coordination of multiwavelength and multifacility observations

MultiWaveLink is an interactive, computerized data base that was developed to facilitate a multiwavelength approach to studying astrophysical sources. It can be used to access information about multiwavelength resources (observers, telescopes, data bases and analysis facilities) or to organize observing campaigns that require either many telescopes operating in different spectral regimes or a network of similar telescopes circumspanning the Earth.     

Near real-time PPP-based monitoring of the ionosphere using dual-frequency GPS/BDS/Galileo data

Ionosphere delay is very important to GNSS observations, since it is one of the main error sources which have to be mitigated even eliminated in order to determine reliable and precise positions. The ionosphere is a dispersive medium to radio signal, so the value of the group delay or phase advance of GNSS radio signal depends on the signal frequency. Ground-based GNSS stations have been used for ionosphere monitoring and modeling for a long time. In this paper we will introduce a novel approach suitable for single-receiver operation based on the precise point positioning (PPP) technique. One of the main characteristic is that only carrier-phase observations are used to avoid particular effects of pseudorange observations. The technique consists of introducing ionosphere ambiguity parameters obtained from PPP filter into the geometry-free combination of observations to estimate ionospheric delays. Observational data from stations that are capable of tracking the GPS/BDS/GALILEO from the International GNSS Service (IGS) Multi-GNSS Experiments (MGEX) network are processed. For the purpose of performance validation, ionospheric delays series derived from the novel approach are compared with the global ionospheric map (GIM) from Ionospheric Associate Analysis Centers (IAACs). The results are encouraging and offer potential solutions to the near real-time ionosphere monitoring.     

A proportional counter/image intensifier system for the detection of low energy X-rays

Proportional counters (PC's) and gas scintillation proportional counters (GSPC's) currently used for detection of low energy X-rays provide information on event position and energy. Although at 1.5 keV PC's have good position resolution (～ 200 μm FWHM) they have relatively poor energy resolution (～ 40% FWHM). Conversely GSPC's have reasonable energy resolution (～ 20% FWHM), but poor spatial resolution (～ 1mm FWHM).We describe a scheme in which a parallel plate PC with a transparent anode deposited on a fibre optic substrate has been used. This allows the light emitted by electron avalanches caused by X-ray events in the PC to be detected by an image intensifier with electronic readout. Using this scheme spatial resolution better than that of conventional PC's should be attainable. In addition avalanches induced by single electrons can be resolved through observation of the time structure of the light flash. Using the ability to count the number of primary electrons created by each X-ray event, it is shown that energy resolution can be achieved which is comparable to that of the GSPC.     

Sensitivity analysis for air temperature profile estimation methods around the tropopause using simulated Aqua/AIRS data

A sensitivity analysis is performed to investigate potential improvements to the accuracy of air temperature profile retrievals near the tropopause. A simple inversion method is employed to identify and remove redundant spectral channels from the retrievals using simulated data for the high-spectral resolution sounder AIRS (Atmospheric Infrared Sounder) on the Aqua satellite. Bayesian optimal theory and inverse technique are applied for the atmospheric temperature profile retrievals, and the 15 μm CO₂ absorption bands (620–750 cm⁻¹) are chosen for this study. Sequentially elimination of redundant channels is directly integrated into the inverse scheme for the temperature profile, in order to accurately retain the valuable channels and remove all the redundant channels, for accurate retrieval of the temperature profile. Also, the tropopause and troposphere are treated differently in the inverse scheme to improve the retrieval accuracy in the tropopause. Results of a sensitivity analysis based on this method, for the Tropical and Middle-Latitude Summer models simulated by MODTRAN4.0, show that the estimated accuracies are improved by 2 K around the tropopause, and are only changed by less than 0.2 K in the troposphere.     

Chlorella vulgaris culture as a regulator of CO2 in a bioregenerative life support system

It is the primary task for a bioregenerative life support system (BLSS) to maintain the stable concentrations of CO₂ and O₂. However, these concentrations could fluctuate based on various factors, such as the imbalance between respiration/assimilation quotients of the heterotrophic and autotrophic components. They can even be out of balance through catastrophic failure of higher plants in the emergency conditions. In this study, the feasibility of using unicellular Chlorella vulgaris of typically rapid growth as both “compensatory system” and “regulator” to control the balance of CO₂ and O₂ was analyzed in a closed ecosystem. For this purpose, a small closed ecosystem called integrative experimental system (IES) was established in our laboratory where we have been conducting multi-biological life support system experiments (MLSSE). The IES consists of a closed integrative cultivating system (CICS) and a plate photo-bioreactor. Four volunteers participated in the study for gas exchange by periodical breathing through a tube connected with the CICS. The plate photo-bioreactor was used to cultivate C. vulgaris. Results showed that the culture of C. vulgaris could be used in a situation of catastrophic failure of higher plant under the emergencies. And the productivity could recover itself to the original state in 3 to 5 days to protect the system till the higher plant was renewed. Besides, C. vulgaris could grow well and the productivity could be affected by the light intensity which could help to keep the balance of CO₂ and O₂ in the IES efficiently. Thus, C. vulgaris could be included in the design of a BLSS as a “compensatory system” in the emergency contingency and a “regulator” during the normal maintenance.