Classification of multitemporal Sentinel-2 data for field-level monitoring of rice cropping practices in Taiwan

Cadastral information of rice fields is important for monitoring cropping practices in Taiwan due to official initiatives. Remote sensing based rice monitoring has been a challenge for years because the size of rice fields is small, and crop mapping requires information of crop phenology, relating to spatiotemporal resolution of satellite data. This study aims to develop an approach for mapping rice-growing areas at field level using multi-temporal Sentinel-2 data in Taiwan. The data were processed for 2018, following four main steps: (1) construct time-series Enhanced Vegetation Index (EVI) and Land Surface Water Index (LSWI), (2) noise filtering of the time-series data using wavelet transform, (3) rice crop classification using information of crop phenology, and (4) parcel-based accuracy assessment of the mapping results. The parcel-to-parcel comparisons between mapping results and ground reference data indicated satisfactory results. These findings were confirmed by close agreement between satellite-derived rice area and government’s statistics. Although some factors, including mixed-pixel issues and cloud-cover effects, lowered the mapping accuracies of townships along the coastline, this study has demonstrated the efficacy of using multitemporal Sentinel-2 data to create a reliable database of rice-growing areas over a large and heterogeneous region. Such a quantitative information was important for updating rice crop maps and monitoring cropping practices.     

Study on extraction of crop information using time-series MODIS data in the Chao Phraya Basin of Thailand

In order to acquire the crop-related information in Chao Phraya Basin, time-series MODIS data were used in this paper. Although the spatial resolution of MODIS data is not very high, it is still useful for detecting very large-scale phenomenon, such as changes in seasonal vegetation patterns. After the data processing a general crop-related LULC (land use and land cover) map, cropping intensity map and cropping patterns map were produced. Analysis of these maps showed that the main land use type in the study area was farmland, most of which was dominated by rice. Rice fields mostly concentrated in the flood plains and double or triple rice-cropping system was commonly employed in this area. Maize, cassava, sugarcane and other upland crops were mainly distributed in the high alluvial terraces. Because these area often have water shortage problem particularly in the dry season which can support only one crop in a year, the cropping intensity was very low. However, some upland areas can be cultivated twice a year with crops which have short growing seasons. The crop information extracted from MODIS data sets were assessed by CBERS data, statistic data and so on. It was shown that MODIS derived crop information coincided well with the statistic data at the provincial level. At the same time, crop information extracted by MODIS data sets and CBERS were compared with each other which also showed similar spatial patterns.     

A simplified data assimilation method for reconstructing time-series MODIS NDVI data

The Normalized Difference Vegetation Index (NDVI) is an important vegetation index, widely applied in research on global environmental and climatic change. However, noise induced by cloud contamination and atmospheric variability impedes the analysis and application of NDVI data. In this work, a simplified data assimilation method is proposed to reconstruct high-quality time-series MODIS NDVI data. We extracted 16-Day L3 Global 1 km SIN Grid NDVI data sets for western China from MODIS vegetation index (VI) products (MOD13A2) for the period 2003–2006. NDVI data in the first three years (2003–2005) were used to generate the background field of NDVI based on a simple three-point smoothing technique, which captures annual features of vegetation change. NDVI data for 2006 were used to test our method. For every time step, the quality assurance (QA) flags of the MODIS VI products were adopted to empirically determine the weight between the background field and NDVI observations. Ultimately, more reliable NDVI data can be produced. The results indicate that the newly developed method is robust and effective in reconstructing high-quality MODIS NDVI time-series.     

Object-based rice mapping using time-series and phenological data

Remote sensing techniques are often used in mapping rice, but high quality time-series remote sensing data are difficult to obtain due to the cloudy weather of rice growing areas and long satellite revisit interval. As such, rice mapping is usually based on mono-temporal Landsat TM/ETM+ data, which have large uncertainties due to the spectral similarity of different vegetation types. Moreover, conventional pixel-based classification method is unable to meet the required accuracy for rice mapping. Therefore, this study proposes a new strategy for mapping rice in cloud-prone areas using fused data of Landsat-8 OLI time-series and phenological parameters, based on the object-based method. We determine the critical growth stages of paddy rice from observed phenological data and MODIS-NDVI time-series data. The spatial and temporal adaptive reflectance fusion model (STARFM) is used to blend the MODIS and Landsat data to obtain a multi-temporal Landsat-like dataset for classification. Finally, an object-oriented algorithm is used to extract rice paddies from the Landsat-like, time-series dataset. The validation experiments show that the proposed method can provide high accuracy rice mapping, with an overall accuracy of 92.38% and a kappa coefficient of 0.85.     

Reconstructing pathfinder AVHRR land NDVI time-series data for the Northwest of China

Vegetation typically elicits dynamics at the seasonal and annual level. Time-series of normalized difference vegetation index (NDVI) datasets, such as the pathfinder AVHRR land (PAL) NDVI dataset, have proven to be appropriate for the detection of long-term vegetation cover changes. It has been applied in modelling experiments for terrestrial ecosystems at the global, continental, and regional scales. But some PAL NDVI time series remain significant residual effects and noise levels. A simple method, the mean-value iteration filter (MVI), has been developed to reduce the noise and to enable the reconstruction of high quality NDVI time-series. A comparison between the newly developed method and other existing methods (the modified BISE algorithm and a fast Fourier transform algorithm) indicates that the newly developed method is an effective tool for reconstructing high-quality time series of PAL NDVI time series.     

基于CCD的单缝衍射条纹光强测量与小波去噪技术研究

通过CCD的数据采集对单缝衍射条纹的光强进行测量试验,简要分析了试验中的噪声来源。提出使用小波技术对CCD的采集信号进行噪声滤除,详细讨论了小波的分解与重构以及在去除噪声中的应用。试验结果表明,小波技术能完全滤除噪声。     

Modified gradiometer technique applied to Double Star (TC-1)

On TC-1 (Tan Ce 1), the equatorial spacecraft of the Double Star mission, a strong spin-synchronized magnetic interference from the solar panels was observed. In-flight correction techniques for spinning spacecraft that are based on minimizing spin tones in the spin-aligned component and in the magnitude of the ambient magnetic field are therefore not possible in this case. However, due to the fortunate situation that the spacecraft carries two flux-gate magnetometers on the same boom (at 0.5 m distance from each other), the spacecraft field effects could be removed from the spin-averaged data to achieve 0.2 nT relative accuracy, by using a gradiometer technique. Methodology and results are presented. The obtained accuracy allows the use of the data in multi-spacecraft studies together with the Cluster satellites.     

基于CCD的单缝衍射条纹光强测量与小波去噪技术研究

通过CCD的数据采集对单缝衍射条纹的光强进行测量试验,简要分析了试验中的噪声来源.提出使用小波技术对CCD的采集信号进行噪声滤除,详细讨论了小波的分解与重构以及在去除噪声中的应用.试验结果表明,小波技术能完全滤除噪声.     

基于经验模分解的PPP时间比对数据消噪方法

精密单点定位(PPP)时间比对数据会受到观测噪声的影响,因此对时间比对数据进行消噪是一项重要工作。提出一种基于经验模分解(EMD)的PPP时间比对数据消噪方法,并将该方法与Vondrak滤波方法的消噪效果进行对比。结果表明,两种方法的消噪效果相当,均能有效滤除PPP时间比对数据中的随机噪声,明显改善时间比对的稳定度。     

基于经验模分解的PPP时间比对数据消噪方法

精密单点定位(PPP)时间比对数据会受到观测噪声的影响,因此对时间比对数据进行消噪是一项重要工作。提出一种基于经验模分解(EMD)的PPP时间比对数据消噪方法,并将该方法与Vondrak滤波方法的消噪效果进行对比。结果表明,两种方法的消噪效果相当,均能有效滤除PPP时间比对数据中的随机噪声,明显改善时间比对的稳定度。     

Noise analysis of continuous GPS coordinate time series for CMONOC

The Crustal Movement Observation Network of China (CMONOC) is one of the major scientific infrastructures, mainly using Global Positioning System (GPS) measurements, to monitor crustal deformation in the Chinese mainland. In this paper, decade-long coordinate time series of 26 continuous GPS sites of CMONOC are analyzed for their noise content using maximum likelihood estimation (MLE). We study the noise properties of continuous GPS time series of CMONOC for the unfiltered, filtered solutions and also the common mode signals in terms of power law plus white noise model. In the spatial filtering, we remove for every time series a common mode error that was estimated from a modified stacking of position residuals from other sites within ∼1000 km of the selected site. We find that the common mode signal in our network has a combination of spatially correlated flicker noise and a common white noise with large spatial extent. We demonstrate that for the unfiltered solutions of CMONOC continuous GPS sites the main colored noise is a flicker process, with a mean spectral index of ∼1. For the filtered solutions, the main colored noise is a general power law process, indicating that a major number of the filtered regional solutions have a combination of noise sources or local effects. The velocity uncertainties from CMONOC continuous GPS coordinate time series may be underestimated by factors of 8–16 if a pure white noise model is assumed. In addition, using a white plus flicker noise model, the median values of velocity errors for the unfiltered solutions are 0.16 (north), 0.17 (east) and 0.58 (vertical) mm/yr, and the median values for the filtered solutions are 0.09 (north), 0.10 (east) and 0.40 (vertical) mm/yr.     

转轴扭转变形动态测量噪声消除方法研究

转轴的扭转变形信号中包含了大量与转轴运行状态相关的信息,但往往被噪声所淹没,合理的去噪处理能够得到更有效的信号成分,有助于对转轴的工作状态做出正确的判断。相比于传统滤波方法,基于EMD的去噪方法有较强的自适应性,更适合动态信号的处理。从转轴相互正交的两个方向采集的相位差信号,依据两路信号EMD分解后各阶本征模态函数的相关性去除噪声并恢复出有效信号。仿真验证了该方法的有效性,并将其运用到实际转轴扭转变形信号的处理中,得到了更为准确的扭转变形信息。     

An investigation of high frequency auroral plasma emissions with a tubular dipole antenna

A band of enhanced amplitudes which follows a local plasma frequency fn in raw high frequency (HF) noise spectra is usually related to plasma emissions in the upper hybrid band (fn, fu). The enhanced band in question occurs permanently in noise spectra recorded on the Intercosmos-19, APEX and CORONAS satellites in the altitude range of 500 km–3000 km. For moderately magnetized plasma with fn > 2fc (fc – electron gyro frequency), the band occurs below fn determined from the topside sounder and impedance data or from electron beam induced spectra. The simulations of an equivalent circuit composed of a dipole antenna in a cold plasma and its preamplifiers, determined the physical origin of the band as the passive circuit resonance, due to inductive character of the antenna in a frequency band (fc, fu). The resonance spectral content is highly structured due to an inflight variability of the circuit impedances. In this report we analyze the noise and impedance spectra which are the most typical in an auroral zone if fn > fc. We focus attention on determination of local electron plasma density, essential for provisional HF mode classification. We found that the natural plasma emission in the upper hybrid band does not manifest itself as the banded natural emission, which may be used for reliable determination of local plasma frequency in the altitude range of 500–3000 km. The fast magnetosonic mode predominates in the auroral emissions. The broadband and multi banded electromagnetic emissions extending from the fast magnetosonic band well above fn > fc are characteristic for the strong wave activity and are much less frequent.     

Impact of emissions from anthropogenic sources on satellite-derived reflectance

The present study deals with the impact of extensive anthropogenic activities associated with festivities and agricultural crop residues burning in the Indo-Gangetic Plains (IGP) on satellite-derived reflectance during November 2007. Intense smoke plumes were observed in the IRS-P4 OCM satellite data over IGP associated with agricultural crop residue burning during the study period. Terra-MODIS AOD and CALIPSO LIDAR backscatter datasets were analysed over the region to understand the spatial and temporal variation of the aerosol properties. Ground-based measurements on aerosol optical properties and black carbon (BC) mass concentration were carried out during 7–14 November 2007 over urban region of Hyderabad, India. Top of the Atmosphere (TOA) reflectance estimated from IRS-P4 Ocean Color Monitor (OCM) data showed large variations due to anthropogenic activities associated with crop residue burning and fireworks. Atmospheric corrections to OCM satellite data using 6S radiative transfer code with inputs from ground and satellite measurements could account for the variations due to differential aerosol loading. Results are discussed in this paper.     

Hierarchical object oriented classification using very high resolution imagery and LIDAR data over urban areas

Urban land cover information extraction is a hot topic within urban studies. Heterogeneous spectra of high resolution imagery—caused by the inner complexity of urban areas—make it difficult. In this paper a hierarchical object oriented classification method over an urban area is presented. Combining QuickBird imagery and light detection and ranging (LIDAR) data, nine kinds of land cover objects were extracted. The Spectral Shape Index (SSI) method is used to distinguish water and shadow from black body mask, with 100% classification accuracy for water and 95.56% for shadow. Vegetation was extracted by using a Normalized Difference Vegetation Index (NDVI) image at first, and then a more accurate classification result of shrub and grassland is obtained by integrating the height information from LIDAR data. The classification accuracy of shrub was improved from 85.25% to 92.09% and from 82.86% to 97.06% for grassland. More granularity of this classification can be obtained by using this method. High buildings and low buildings can, for example, be distinguished from the original building class. Road class can also be further classified into roads and crossroads. The comparison of the classification accuracy between this method and the traditional pixel-based method indicates that the total accuracy is improved from 69.12% to 89.40%.     

基于双谱分析的高原灌木回波信号分类

目标的回波信号是无线电引信获得目标信息的最重要方式,为了太赫兹频段的引信前端未来能够投入高原战场,适应高原不同的地貌环境,利用双谱对高原在灌木地形下不同高度的太赫兹波回波特性进行了分析。为减少分类时间,对双谱数据进行积分,得到实采信号双谱切片的特征,进而利用最邻近算法对此进行分类。利用经验模态分解(EMD)提取原始数据内在模态函数的特征,再次分类并与前一组分类结果进行对比。通过一系列数据的分类,结果表明：利用一维的积分双谱信息可以有效提取出距离地面分别为2 m、3 m、4 m、5 m时的特征并进行分类,经验模态分解也可以有效提高分类的成功率,成功率最高可达90%以上。     

基于EMD的MEMS陀螺仪随机漂移分析方法

为了抑制微机械电子系统（MEMS）陀螺仪的随机漂移，基于经验模态分解（EMD）和模态集合选择标准，结合时间序列建模滤波法，提出了一种改进的MEMS陀螺仪随机漂移分析方法。首先，通过EMD将MEMS陀螺仪原始数据分解为多个本征模态函数（IMF），利用模态集合选择标准将IMF分为噪声IMF、噪声与信号混合IMF和信号IMF三类；然后，对混合IMF进行重构、时间序列建模及自适应卡尔曼滤波（AKF）；最后，将3类信号重构，实现MEMS陀螺仪信号去噪。实验表明:所提方法有更好的去噪效果和实时性，提高了MEMS陀螺仪的使用精度。      

Cross testing of ionosphere models IRI-2001 and IRI-2007, data from satellite altimeters (Topex/Poseidon and Jason-1) and global ionosphere maps

In this study, we use a great body of statistical data covering the entire 23rd solar cycle to cross test data of satellite altimeters, Global Ionosphere Maps and the International Reference Ionosphere models, IRI-2001 and IRI-2007. It is revealed that experimental TEC values of the satellite altimeters regularly exceed the model ones by ∼3 TECU (1 TECU = 10¹⁶ m⁻²). The best possible value of difference between TECs obtained from altimeter and GIM-map data significantly differs for different laboratories: the maximum for CODG data falls on 2.5 TECU, ESAG – 3 TECU, JPLG – 0 TECU, UPCG – 2 TECU. The dependence of experimental and model data root-mean-square deviation on the F10.7 index is shown to be nearly linear. IRI-2001 and IRI-2007 relative errors are characterized by considerable 11-year and annual variations. Given the geomagnetic planetary index K_p under 7, IRI-2001 and IRI-2007 reproduce TEC in the ionosphere with an accuracy of ∼30% relative to measurement data from satellite altimeters. The amplitude of absolute error variations resulting from the difference in ionization enhancement between the model and the real ionosphere during the morning solar terminator transit is ∼5 TECU.     

A comparative analysis of SLR,MLR, ANN,XGBoost and CNN for crop height estimation of sunflower using Sentinel-1 and Sentinel-2

Sustainable monitoring and determining the biophysical characteristics of crops is of global importance due to the increase in demand for food. In this context, remote sensing data provide valuable information on crops. This study investigates the relationship between the variables determined from both Synthetic Aperture Radar (SAR) and optical images and crop height. For this purpose, backscatter (σVH, σVV, σVH / σVV) and coherence (?VH, ?VV) of multi-temporal dual-polarized Sentinel-1 and vegetation indices of multi-temporal Sentinel-2 data are analyzed. Two indices, namely, Normalized Difference Vegetation Index (NDVI) and NDVI with the red-edge band (NDVIred), are interpreted to identify the contribution of the red-edge band over the near-infrared band. The Zile District of Tokat province in Turkey where dominantly sunflower cultivation is carried out, was selected as the study area. In the analysis of the data, Simple Linear Regression (SLR), Multiple Linear Regression (MLR), Artificial Neural Network (ANN), EXtreme Gradient Boosting (XGBoost), and Convolutional Neural Network (CNN) were used. In the results of the study, ANN showed the lowest RMSE = 3.083 cm (RMSE%= 11.284) in the stem elongation period. The CNN followed the lowest RMSE for the Inflorescence development and flowering stages 19.223 cm (RMSE%=15.458) and 8.731 cm (RMSE%=5.821), respectively. In the ripening period, XGBoost achieved the lowest RMSE = 8.731 cm (RMSE%=6.091). All the best models in four methods were created using common variables of σVH, σVV, ?VH, ?VV and NDVIred, except ANN which exclude coherence variables. The results concluded that NDVIred contributed more than NDVI which is widely interpreted in previous studies.     

Validation of the STORM model in IRI-2001 at a high latitude station

The International Reference Ionosphere IRI-2001 model contains geomagnetic activity dependence based on an empirical storm time ionospheric correction (STORM model). An extensive validation of the STORM model for the middle latitude region has been performed. In this paper the ability of the STORM model to predict foF2 values at high latitudes is analyzed. For this, ionosonde data obtained at Base Gral. San Martin (68.1°S, 293°E) are compared with those obtained by the IRI-2001 model with or without storm correction during four geomagnetic storms that occurred in 2000 (Rz₁₂ = 117) and 2001 (Rz₁₂ = 111). The results show that predicted values with the STORM model follow the behaviour of foF2 experimental data better than without the STORM model. The relative deviation between measured and predicted foF2 reaches values of up to 24% and 43% with and without the STORM model in IRI-2001, during the main phase of the storms. In order to explain increases of electron density that occurred prior to the storm onset and also decreases of electron density observed during the first part of the recovery of the storm, possible physical mechanisms are discussed.