Orbit determination performances using single- and double-differenced methods: SAC-C and KOMPSAT-2

In this paper, Global Positioning System-based (GPS) Orbit Determination (OD) for the KOrea-Multi-Purpose-SATellite (KOMPSAT)-2 using single- and double-differenced methods is studied. The requirement of KOMPSAT-2 orbit accuracy is to allow 1 m positioning error to generate 1-m panchromatic images. KOMPSAT-2 OD is computed using real on-board GPS data. However, the local time of the KOMPSAT-2 GPS receiver is not synchronized with the zero fractional seconds of the GPS time internally, and it continuously drifts according to the pseudorange epochs. In order to resolve this problem, an OD based on single-differenced GPS data from the KOMPSAT-2 uses the tagged time of the GPS receiver, and the accuracy of the OD result is assessed using the overlapping orbit solution between two adjacent days. The clock error of the GPS satellites in the KOMPSAT-2 single-differenced method is corrected using International GNSS Service (IGS) clock information at 5-min intervals. KOMPSAT-2 OD using both double- and single-differenced methods satisfies the requirement of 1-m accuracy in overlapping three dimensional orbit solutions. The results of the SAC-C OD compared with JPL’s POE (Precise Orbit Ephemeris) are also illustrated to demonstrate the implementation of the single- and double-differenced methods using a satellite that has independent orbit information available for validation.     

Image restoration of the asymmetric point spread function of a high-resolution remote sensing satellite with time-delayed integration

The calibration and validation team of the Korea Aerospace Research Institute has calibrated and validated the image data of the KOMPSAT-2 since the launch of KOMPSAT-2 on July 28, 2006. The asymmetric phenomenon of the point spread function of the KOMPSAT-2 image data is evident in both the along and the across direction, most likely because KOMPSAT-2 has a 1 m ground sample distance high-resolution camera with time-delayed integration. Furthermore, because KOMPSAT-2 is in space, the KOMPSAT-2 image data has been corrected with good results by means of modulation transfer function compensation.     

卫星遥感图像信息作为姿态敏感器的应用研究

陆标敏感器可以从卫星遥感图像信息中提取卫星的姿态信息,它根据卫星实时图像与基准图像之间的偏移量,计算出卫星的姿态。对陆标敏感器的总体方案进行了设计,在图像匹配技术中采用基于区域特征的先粗匹配、后精匹配的匹配算法,在姿态确定算法中对可能达到的姿态测量精度进行了理论分析,研究表明陆标敏感器可以获得较高精度的卫星姿态信息,作为新型光学姿态敏感器具有重要的研究意义。     

卫星遥感图像数据压缩质量评价研究

针对高保真高速数据压缩算法的卫星应用,采用客观参数评价和主观评价相结合的方法,对卫星遥感图像数据压缩质量评价进行了较系统的研究。通过四次试验实践表明,主观评价具有相当的客观性,是目前遥感图像数据压缩质量评价中最基本、最有效的方法,科学合理地设计主观评价至关重要;客观参数评价是有益的参考方法,在改善某一压缩方法时,峰值信噪比(PSNR)是个较好的衡量标准;主观评价和客观评价有较好的一致性,主要体现在PSNR、相似度(XSD)和全图局域最大误差(LME)等客观评价参数上。     

Hardware-in-the-loop simulations of GPS-based navigation and control for satellite formation flying

A relative navigation and formation control algorithm for satellite formation flying was developed, and a hardware-in-the-loop (HIL) simulation testbed was established and configured to evaluate this algorithm. The algorithm presented is a relative navigation estimation algorithm using double-difference carrier-phase and single-difference code measurements based on the extended Kalman filter (EKF). In addition, a state-dependent Riccati equation (SDRE) technique is utilized as a nonlinear controller for the formation control problem. The state-dependent coefficient (SDC) form is formulated to include nonlinearities in the relative dynamics. To evaluate the relative navigation and control algorithms developed, a closed-loop HIL testbed is configured. To demonstrate the performance of the testbed, a test formation flying scenario comprising formation acquisition and keeping in a low earth orbit (LEO) has been established. The relative navigation results from the closed-loop simulations show that a 3D RMS of 0.07 m can be achieved for position accuracy. The targeted leader–follower formation flying in the along-track separation of 100 m was maintained with a mean position error of approximately 0.2 m and a standard deviation of 0.9 m. The simulation results show that the HIL testbed is capable of successful demonstration of the GPS-based satellite autonomous formation flying mission.     

Improved kHz-SLR tracking techniques and orbit quality analysis for LEO missions

Satellite gravity field missions such as CHAMP, GRACE and GOCE are designed as low Earth orbiting spacecraft (LEO) with orbit heights of about 250–500 km. The challenging mission objectives require a very precise knowledge of the satellite orbit position in space. For these missions precise orbit information is typically provided by GPS satellite-to-satellite tracking (SST) observations supported by satellite laser ranging (SLR).     

Comparison of fragments created by low- and hyper-velocity impacts

This paper summarizes two new satellite impact experiments. The objective of the experiments was to investigate the outcome of low- and hyper-velocity impacts on two identical target satellites. The first experiment was performed at a low-velocity of 1.5 km/s using a 40-g aluminum alloy sphere. The second experiment was performed at a hyper-velocity of 4.4 km/s using a 4-g aluminum alloy sphere. The target satellites were 15 cm × 15 cm × 15 cm in size and 800 g in mass. The ratios of impact energy to target mass for the two experiments were approximately the same. The target satellites were completely fragmented in both experiments, although there were some differences in the characteristics of the fragments. The projectile of the low-velocity impact experiment was partially fragmented while the projectile of the hyper-velocity impact experiment was completely fragmented beyond recognition. To date, approximately 1500 fragments from each impact experiment have been collected for detailed analysis. Each piece has been weighed, measured, and analyzed based on the analytic method used in the NASA Standard Breakup Model (2000 revision). These fragments account for about 95% of the target mass for both impact experiments. Preliminary analysis results will be presented in this paper.     

Altimetry,bathymetry and geoid variations at the Gavdos permanent Cal/Val facility

The aim of this work has been to examine the relationship of steep bathymetry in the coastal areas around the permanent Cal/Val facility of Gavdos, and their influence on the produced calibration values for the Jason-2 satellite altimeter. The paper describes how changes in seafloor topography (from 200 to 3500 m depth over a distance of 10 km) are reflected on the determined altimeter parameters using different reference surfaces for satellite calibration. Finally, it describes the relation between these parameter trends and the region’s local characteristics.     

Spin axis orientation of Ajisai determined from Graz 2 kHz SLR data

The Graz 2 kHz Satellite Laser Ranging (SLR) measurements allow determination of the spin axis orientation of the geodetic satellite Ajisai. The high repetition rate of the laser makes it possible to determine the epoch time when the laser is pointing directly between two corner cube reflector (CCR) rings of the satellite. Identification of many such events during a few (up to 3) consecutive passes allows to state the satellite orientation in the celestial coordinate system. Six years of 2 kHz SLR data (October 2003–October 2009) delivered 331 orientation values which clearly show precession of the axis along a cone centered at 14^h56^m2.8^s in right ascension and 88.512° in declination (J2000.0 celestial reference frame) and with an half-aperture angle θ of 1.405°. The spin axis precesses with a period of 117 days, which is equal to the period of the right ascension of the ascending node of Ajisai’s orbit. We present a model of the axis precession which allows prediction of the satellite orientation – necessary for the envisaged laser time transfer via Ajisai mirrors.     

22 Years of AJISAI spin period determination from standard SLR and kHz SLR data

Satellite Laser Ranging (SLR) is a powerful and efficient technique to measure spin parameters of satellites equipped with corner cube reflectors. We obtained spin period determination of the satellite AJISAI from SLR data only: 17246 pass-by-pass estimates from standard 1–15 Hz SLR data (14/Aug/1986–30/Dec/2008) and 1444 pass-by-pass estimates (9/Oct/2003–30/Dec/2008) from data of the first 2 kHz SLR system from Graz, Austria. A continuous history of the slowing down of AJISAI spin is derived from frequency analysis, and corrected for the apparent effects. The apparent corrections, elaborated here, allowed very accurate determination of AJISAI initial spin period: 1.4855 ± 0.0007 [s]. The paper identifies also non-gravitational effects as a source of the periodical changes in the rate of slowing down of the satellite.     

Mesospheric OH temperatures: Simultaneous ground-based and SABER OH measurements over Millstone Hill

We present rotational temperature measurements of the mesospheric OH emission layer using a meridional imaging spectrograph at Millstone Hill (42.6°N, 72.5°W). The system is equipped with a state-of-the-art bare-CCD detector and can yield simultaneous quasi-meridional images of the mesospheric OH and O₂ intensity and temperature fields at 87 and 94 km altitude during the course of each night. A cross-validation study of the rotational OH temperature measurements obtained on 61 nights during the autumnal months of 2005–2007 was undertaken with near-simultaneous kinetic temperature measurements made by the SABER instrument aboard the NASA TIMED satellite during overpasses of Millstone Hill. Excellent agreement was obtained between the two datasets with the small differences being attributable to differences in the spatial and temporal averaging inherent between the two datasets.     

Measuring atmospheric density using GPS–LEO tracking data

We present a method to estimate the total neutral atmospheric density from precise orbit determination of Low Earth Orbit (LEO) satellites. We derive the total atmospheric density by determining the drag force acting on the LEOs through centimeter-level reduced-dynamic precise orbit determination (POD) using onboard Global Positioning System (GPS) tracking data. The precision of the estimated drag accelerations is assessed using various metrics, including differences between estimated along-track accelerations from consecutive 30-h POD solutions which overlap by 6 h, comparison of the resulting accelerations with accelerometer measurements, and comparison against an existing atmospheric density model, DTM-2000. We apply the method to GPS tracking data from CHAMP, GRACE, SAC-C, Jason-2, TerraSAR-X and COSMIC satellites, spanning 12 years (2001–2012) and covering orbital heights from 400 km to 1300 km. Errors in the estimates, including those introduced by deficiencies in other modeled forces (such as solar radiation pressure and Earth radiation pressure), are evaluated and the signal and noise levels for each satellite are analyzed. The estimated density data from CHAMP, GRACE, SAC-C and TerraSAR-X are identified as having high signal and low noise levels. These data all have high correlations with anominal atmospheric density model and show common features in relative residuals with respect to the nominal model in related parameter space. On the contrary, the estimated density data from COSMIC and Jason-2 show errors larger than the actual signal at corresponding altitudes thus having little practical value for this study. The results demonstrate that this method is applicable to data from a variety of missions and can provide useful total neutral density measurements for atmospheric study up to altitude as high as 715 km, with precision and resolution between those derived from traditional special orbital perturbation analysis and those obtained from onboard accelerometers.     

Sodium lidar measurements of mesopause region temperatures at 23° S

A sodium lidar, capable of measuring temperature in the 80–100 km region, has been in operation at São José dos Campos (23° S, 46 W) since March 2007. Good quality data have been obtained for late autumn, winter and spring, but weather conditions make it extremely difficult to make measurements from mid-November to mid- February. We find the temperature structure to be strongly modulated by tides and gravity waves, but average profiles typically show a primary mesopause height close to 100 km with temperatures around 180 K, and a tendency for a secondary minimum of about 185 K to occur close to 90 km. Vertical temperature gradients greater than 50 K/km are sometimes seen even on profiles averaged over several hours. The strongest gradients are always positive and are frequently associated with strong gradients in sodium concentration. On the other hand, we frequently see rapid changes in the temperature profile, suggesting that models and non-local temperature measurements, as made by satellite radiometers, for example, are of little use in applications such as the analysis of gravity wave propagation seen in airglow images.     

Development and testing of a statistical texture model for land cover classification of the Black Sea region with MODIS imagery

A statistical model is proposed for analysis of the texture of land cover types for global and regional land cover classification by using texture features extracted by multiresolution image analysis techniques. It consists of four novel indices representing second-order texture, which are calculated after wavelet decomposition of an image and after texture extraction by a new approach that makes use of a four-pixel texture unit. The model was applied to four satellite images of the Black Sea region, obtained by Terra/MODIS and Aqua/MODIS at different spatial resolution. In single texture classification experiments, we used 15 subimages (50 × 50 pixels) of the selected classes of land covers that are present in the satellite images studied. These subimages were subjected to one-level and two-level decompositions by using orthonormal spline and Gabor-like spline wavelets. The texture indices were calculated and used as feature vectors in the supervised classification system with neural networks. The testing of the model was based on the use of two kinds of widely accepted statistical texture quantities: five texture features determined by the co-occurrence matrix (angular second moment, contrast, correlation, inverse difference moment, entropy), and four statistical texture features determined after the wavelet transformation (mean, standard deviation, energy, entropy). The supervised neural network classification was performed and the discrimination ability of the proposed texture indices was found comparable with that for the sets of five GLCM texture features and four wavelet-based texture features. The results obtained from the neural network classifier showed that the proposed texture model yielded an accuracy of 92.86% on average after orthonormal wavelet decomposition and 100% after Gabor-like wavelet decomposition for texture classification of the examined land cover types on satellite images.     

An approach to multi-temporal MODIS image analysis using image classification and segmentation

This paper discusses an approach for river mapping and flood evaluation based on multi-temporal time series analysis of satellite images utilizing pixel spectral information for image classification and region-based segmentation for extracting water-covered regions. Analysis of MODIS satellite images is applied in three stages: before flood, during flood and after flood. Water regions are extracted from the MODIS images using image classification (based on spectral information) and image segmentation (based on spatial information). Multi-temporal MODIS images from “normal” (non-flood) and flood time-periods are processed in two steps. In the first step, image classifiers such as Support Vector Machines (SVM) and Artificial Neural Networks (ANN) separate the image pixels into water and non-water groups based on their spectral features. The classified image is then segmented using spatial features of the water pixels to remove the misclassified water. From the results obtained, we evaluate the performance of the method and conclude that the use of image classification (SVM and ANN) and region-based image segmentation is an accurate and reliable approach for the extraction of water-covered regions.     

Determination of station positions and velocities from laser ranging observations to Ajisai,Starlette and Stella satellites

The positions and velocities of the four Satellite Laser Ranging (SLR) stations: Yarragadee (7090), Greenbelt (7105), Graz (7839) and Herstmonceux (7840) from 5-year (2001–2005) SLR data of low orbiting satellites (LEO): Ajisai, Starlette and Stella were determined. The orbits of these satellites were computed from the data provided by 20 SLR stations. All orbital computations were performed by means of NASA Goddard’s GEODYN-II program. The geocentric coordinates were transformed to the topocentric North–South, East–West and Vertical components in reference to ITRF2005. The influence of the number of normal points per orbital arc and the empirical acceleration coefficients on the quality of station coordinates was studied. To get standard deviation of the coordinates determination lower than 1 cm, the number of the normal points per site had to be greater than 50. The computed positions and velocities were compared to those derived from LAGEOS-1/LAGEOS-2 data. Three parameters were used for this comparison: station coordinates stability, differences from ITRF2005 positions and velocities. The stability of coordinates of LEO satellites is significantly worse (17.8 mm) than those of LAGEOS (7.6 mm), the better results are for Ajisai (15.4 mm) than for Starlette/Stella (20.4 mm). The difference in positions between the computed values and ITRF2005 were little bit worse for Starlette/Stella (6.6 mm) than for LAGEOS (4.6 mm), the results for Ajisai were five times worse (29.7 mm) probably due to center of mass correction of this satellite. The station velocities with some exceptions were on the same level (≈1 mm/year) for all satellites. The results presented in this work show that results from Starlette/Stella are better than those from Ajisai for station coordinates determination. We can applied the data from LEO satellites, especially Starlette and Stella for determination of the SLR station coordinates but with two times lower accuracy than when using LAGEOS data.     

Limb Viewing Hyper Spectral Imager (LiVHySI) for airglow measurements onboard YOUTHSAT-1

The Limb Viewing Hyper Spectral Imager (LiVHySI) is one of the Indian payloads onboard YOUTHSAT (inclination 98.73°, apogee 817 km) launched in April, 2011. The Hyper-spectral imager has been operated in Earth’s limb viewing mode to measure airglow emissions in the spectral range 550–900 nm, from terrestrial upper atmosphere (i.e. 80 km altitude and above) with a line-of-sight range of about 3200 km. The altitude coverage is about 500 km with command selectable lowest altitude. This imaging spectrometer employs a Linearly Variable Filter (LVF) to generate the spectrum and an Active Pixel Sensor (APS) area array of 256 × 512 pixels, placed in close proximity of the LVF as detector. The spectral sampling is done at 1.06 nm interval. The optics used is an eight element f/2 telecentric lens system with 80 mm effective focal length. The detector is aligned with respect to the LVF such that its 512 pixel dimension covers the spectral range. The radiometric sensitivity of the imager is about 20 Rayleigh at noise floor through the signal integration for 10 s at wavelength 630 nm. The imager is being operated during the eclipsed portion of satellite orbits. The integration in the time/spatial domain could be chosen depending upon the season, solar and geomagnetic activity and/or specific target area. This paper primarily aims at describing LiVHySI, its in-orbit operations, quality, potential of the data and its first observations. The images reveal the thermospheric airglow at 630 nm to be the most prominent. These first LiVHySI observations carried out on the night of 21st April, 2011 are presented here, while the variability exhibited by the thermospheric nightglow at O(¹D) 630 nm has been described in detail.     

Surface temperature estimation in Singhbhum Shear Zone of India using Landsat-7 ETM+ thermal infrared data

Land surface temperature (LST) is an important factor in global change studies, heat balance and as control for climate change. A comparative study of LST over parts of the Singhbhum Shear Zone in India was undertaken using various emissivity and temperature retrieval algorithms applied on visible and near infrared (VNIR), and thermal infrared (TIR) bands of high resolution Landsat-7 ETM+ imagery. LST results obtained from satellite data of October 26, 2001 and November 2, 2001 through various algorithms were validated with ground measurements collected during satellite overpass. In addition, LST products of MODIS and ASTER were compared with Landsat-7 ETM+ and ground truth data to explore the possibility of using multi-sensor approach in LST monitoring. An image-based dark object subtraction (DOS3) algorithm, which is yet to be tested for LST retrieval, was applied on VNIR bands to obtain atmospheric corrected surface reflectance images. Normalized difference vegetation index (NDVI) was estimated from VNIR reflectance image. Various surface emissivity retrieval algorithms based on NDVI and vegetation proportion were applied to ascertain emissivities of the various land cover categories in the study area in the spectral range of 10.4–12.5 μm. A minimum emissivity value of about 0.95 was observed over the reflective rock body with a maximum of about 0.99 over dense forest. A strong correlation was established between Landsat ETM+ reflectance band 3 and emissivity. Single channel based algorithms were adopted for surface radiance and brightness temperature. Finally, emissivity correction was applied on ‘brightness temperature’ to obtain LST. Estimated LST values obtained from various algorithms were compared with field ground measurements for different land cover categories. LST values obtained after using Valor’s emissivity and single channel equations were best correlated with ground truth temperature. Minimum LST is observed over dense forest as about 26 °C and maximum LST is observed over rock body of about 38 °C. The estimated LST showed that rock bodies, bare soils and built-up areas exhibit higher surface temperatures, while water bodies, agricultural croplands and dense vegetations have lower surface temperatures during the daytime. The accuracy of the estimated LST was within ±2 °C. LST comparison of ASTER and MODIS with Landsat has a maximum difference of 2 °C. Strong correlation was found between LST and spectral radiance of band 6 of Landsat-7 ETM+. Result corroborates the fact that surface temperatures over land use/land cover types are greatly influenced by the amount of vegetation present.     

星载合成孔径雷达的分辨特性

较全面和详细地论述了星载合成孔径雷达（ＳＡＲ）的分辨特性;导出了距离分辨率、方位分辨率和辐射分辨率的定量关系式;阐述了多种辐射分辨率定义的内涵,并对其错误概率进行了比较;介绍了ＳＡＲ图像的分辨体积和可懂度的概念,并且完成了ＳＡＲ图像分辨特性的方形等效。这些论述对于ＳＡＲ图像质量的评估和比较,以及系统参数的折衷设计具有重要意义。     

星载高速图像数据压缩原理样机的研制

针对星上高速遥感数据实时压缩要求,选用多模式自适应量化压缩算法,设计了同步并行阵列与流水线相结合的压缩系统体系结构,并以可编程序门阵列物理实现.采用四路压缩系统并行的阵列结构,研制成功高速数据压缩原理样机,数据处理速度大于1000Mbit/s,恢复图像平均峰值信噪比大于40dB.