Substitution of satellite-based land surface temperature defective data using GSP method

Land surface temperature (LST) as an important environmental variable provides valuable information for earth environmental system modelling. Currently, LST is obtained through satellite thermal sensors at various spatial and temporal resolutions. Although spatially continuous satellite-based LST measurements are intended to overcome the shortcomings of sparse ground-based LST measurements, LST images often contain anomalous values due to the existence of clouds or sensor malfunctioning. The problem becomes more serious where the users deal with high spatial resolution characterized by low temporal resolution. This study examines the capability of a newly developed graph signal processing (GSP) method using two-dimensional single-date thermal data. For this purpose, four Landsat/TIRS datasets are analyzed. The data of five elliptical regions on thermal images are eliminated and then reconstructed through the GSP method and using the LST values of the enclosing rectangles containing the ellipsoids. The results indicate that the temperature variation determined by the GSP method generally conforms to the original image LST values. According to a correlation test conducted on the original image LST and those obtained through the GSP method, the values vary from 58% to 95%, which is an above-the-average rate (RMSE from 0.69 to 2.27). The statistical analysis of the original image LST in both the elliptical regions and the enclosing rectangles containing the ellipsoids indicates that an increase in the variance of LST data causes an increased error in the calculation of temperature by the GSP method, and vice versa. The results of the analysis of variance (ANOVA) and Duncan test indicated that an increase in the number of the non-zero spectral bins would result in increased RMSE values for all the dates and the regions. Moreover, the model errors were significant at the 0.05 level across all the image date and five elliptical study regions. Based on the results, the use of this method is recommended for the reconstruction of LST missing values, where dissimilarity of atmospheric conditions limits the use of other methods that depend on the time series data of various dates and a great deal of data calculation.     

A robust object-based shadow detection method for cloud-free high resolution satellite images over urban areas and water bodies

Unwanted contrast in high resolution satellite images such as shadow areas directly affects the result of further processing in urban remote sensing images. Detecting and finding the precise position of shadows is critical in different remote sensing processing chains such as change detection, image classification and digital elevation model generation from stereo images. The spectral similarity between shadow areas, water bodies, and some dark asphalt roads makes the development of robust shadow detection algorithms challenging. In addition, most of the existing methods work on pixel-level and neglect the contextual information contained in neighboring pixels. In this paper, a new object-based shadow detection framework is introduced. In the proposed method a pixel-level shadow mask is built by extending established thresholding methods with a new C₄ index which enables to solve the ambiguity of shadow and water bodies. Then the pixel-based results are further processed in an object-based majority analysis to detect the final shadow objects. Four different high resolution satellite images are used to validate this new approach. The result shows the superiority of the proposed method over some state-of-the-art shadow detection method with an average of 96% in F-measure.     

Evaluation and improvement of coastal GNSS reflectometry sea level variations from existing GNSS stations in Taiwan

Global sea level rise due to an increasingly warmer climate has begun to induce hazards, adversely affecting the lives and properties of people residing in low-lying coastal regions and islands. Therefore, it is important to monitor and understand variations in coastal sea level covering offshore regions. Signal-to-noise ratio (SNR) data of Global Navigation Satellite System (GNSS) have been successfully used to robustly derive sea level heights (SLHs). In Taiwan, there are a number of continuously operating GNSS stations, not originally installed for sea level monitoring. They were established in harbors or near coastal regions for monitoring land motion. This study utilizes existing SNR data from three GNSS stations (Kaohsiung, Suao, and TaiCOAST) in Taiwan to compute SLHs with two methods, namely, Lomb–Scargle Periodogram (LSP)-only, and LSP aided with tidal harmonic analysis developed in this study. The results of both methods are compared with co-located or nearby tide gauge records. Due to the poor quality of SNR data, the worst accuracy of SLHs derived from traditional LSP-only method exceeds 1?m at the TaiCOAST station. With our procedure, the standard deviations (STDs) of difference between GNSS-derived SLHs and tide gauge records in Kaohsiung and Suao stations decreased to 10?cm and the results show excellent agreement with tide gauge derived relative sea level records, with STD of differences of 7?cm and correlation coefficient of 0.96. In addition, the absolute GNSS-R sea level trend in Kaohsiung during 2006–2011 agrees well with that derived from satellite altimetry. We conclude that the coastal GNSS stations in Taiwan have the potential of monitoring absolute coastal sea level change accurately when our proposed methodology is used.     

Multidisciplinary Design and Optimization of Multistage Ground-launched Boost Phase Interceptor Using Hybrid Search Algorithm

This article proposes a multidisciplinary design and optimization (MDO) strategy for the conceptual design of a multistage ground-based interceptor (GBI) using hybrid optimization algorithm, which associates genetic algorithm (GA) as a global optimizer with sequential quadratic programming (SQP) as a local optimizer. The interceptor is comprised of a three-stage solid propulsion system for an exoatmospheric boost phase intercept (BPI). The interceptor's duty is to deliver a kinetic kill vehicle (KKV) to the optimal position in space to accomplish the mission of intercept. The modules for propulsion, aerodynamics, mass properties and flight dynamics are integrated to produce a high fidelity model of the entire vehicle. The propulsion module comprises of solid rocket motor (SRM) grain design, nozzle geometry design and performance prediction analysis. Internal ballistics and performance prediction parameters are calculated by using lumped parameter method. The design objective is to minimize the gross lift off mass (GLOM) of the interceptor under the mission constraints and performance objectives. The proposed design and optimization methodology provide designers with an efficient and powerful approach in computation during designing interceptor systems.     

大型壁板数控钻铆的三点快速调平算法

 针对传统机翼壁板铆接过程中托架调平所存在的问题，提出一种基于有限曲面加工区域的三点调平算法。利用机翼表面上待钻铆点附近三点的坐标来表征该铆接区域的空间姿态，通过钻铆系统的几何结构模型所确定的坐标变化矩阵，逆解出调平所需的各运动副位移增量。在钻铆系统结构模型和机翼曲面方程的基础上建立三点选取的优化模型，并利用外点罚函数法确定其优化解。西安飞机工业(集团)责任有限公司ARJ21新支线飞机的机翼壁板试验件数控钻铆结果表明，该算法能够实现较高的加工精度和调平效率。     

Investigation of the role of aziridine bonding agents on the aging of the composite solid rocket propellant(CSRP)

The role of bonding agents on the aging characteristics is one of the important research topics.Aging program of the prepared propellant samples was conducting as follows:Five samples,two free of bonding agents,and three containing an aziridine based bonding agents(MAPO,HX-752,MAT4),four samples based on different bonding and curing agents all were aged at 70 ℃.The prepared bonding agent "MAT4" gave remarkable improvements and resulted in highly stable mechanical properties comparing with HX-752 or MAPO.The selected bonding agents family inhibited the rate of decomposition of the propellants during the aging periods and supported the propellant matrix against decomposition at the elevated temperatures.      

Systematic examination of the geomagnetic storm sudden commencement using multi resolution analysis

Geomagnetic storm sudden commencements (SSC) contain a wealth of information which is useful in many applications. It is important to point out that the SSCs used in this study are sampled at the rate of one sample per second in order get use of such high resolution data. In this paper, two studies are made on the geomagnetic SSC and an SSC onset automatic detection algorithm is introduced. The first study is about finding the relationship between the SSC rise time and its amplitude. Where it is found that there is a positive correlation between the amplitude and the amplitude gradient which is the amplitude divided by rise time. The second study is checking the spectrum of the SSC, starting from its onset until the end of the SSC rise time. This check had proved that the SSC contains both low and high frequency regions. This led us to introduce a new term, namely the SSC variation rate (VR). This VR is defined as the maximum rate of change of the field in the higher-frequency region of the SSC. These two studies were the guide to build an SSC automatic detector of one sample per second data using multi resolution analysis (MRA) of the discrete wavelet transform (DWT). The data set contains 134 SSCs with different VRs that were collected from the Circum-pan Pacific Magnetometer Network (CPMN). It is found that the standard deviation of the detection error is 41 s and that the average error is 9 s. From the calculated error distribution function, it is found that the detection error is within the range of −1.5 to 3 min. The detection process, as will be shown in the article, takes 70 s for one station and 3 min if the decision is related to the detection(s) of other stations. These results demonstrate the superiority of the proposed algorithm over other algorithms in which the detection error ranges between −8 and 5 min and the detection process takes 2–10 min. In addition to being faster and more accurate than the other algorithms, the proposed algorithm is the first algorithm that automatically detects the SSC onset times from high-resolution data unlike previous studies that focused on determining the SSC times automatically using one-minute resolution data.     

Automatic detection of geomagnetic sudden commencement via time–frequency clusters

One of global processes in ionosphere–thermosphere–magnetosphere system is the geomagnetic storms. It is of great importance to develop an algorithm that auto-detects sudden commencement because it could be an indicator of onset of the geomagnetic storm. Automatic detection of geomagnetic sudden commencement is based on time–frequency clusters generated by spectrogram. Proposed algorithm is tested on data set collected from stations belong to the international real-time magnetic observatory network (INTERMAGNET). Maximum standard deviation of algorithm detection times is observed to be one minute of the corresponding arrival times published by National Geophysical Data Center (NGDC).     

Asymptotic theory for thrusting,spinning-up spacecraft maneuvers

We consider the problem of a spacecraft subjected to constant body-fixed forces and moments about all three axes during a spinning-up, thrusting maneuver. In applications, undesired forces and moments can arise due to thruster imbalances and misalignments and to center-of-mass offset. In previous works, approximate analytical solutions have been found for the attitude motion, and for the change in inertial velocity and inertial position. In this paper we find asymptotic and limiting-case expressions which we derive from the analytic solutions, in order to obtain simplified, practical formulas that lend insight into the motion. Specifically, we investigate how the motion evolves (1) as time grows without bound and (2) for geometric cases of the sphere, the thin rod, and the thin plate. Closed-forms or upper-bound limits are provided for angular velocities, Eulerian angles, angular momentum pointing error, transverse and axial velocities, and transverse and axial displacements. Summaries for the asymptotic limits (for zero initial conditions) are provided in tabular form. Results are verified by numerical simulations.