Generation of Landsat surface temperature product for China, 2000–2010

Land surface temperature (LST) is a key parameter in the physics of land surface processes on regional and global scales. Although there are MODIS and Landsat land surface reflectance products, there is no LST product for Landsat data due in part to many challenges in the development of an operational Landsat LST product generating system because Landsat possesses only one thermal infrared channel. The aim of this article is to describe the Landsat LST product generation project launched by the Centre for Earth Observation and Digital Earth (CEODE), Chinese Academy of Sciences. The generalized single-channel (SC) algorithm proposed by Jiménez-Muñoz et al. is used for LST retrieval. It is fully operational, requires minimal input data requirements, and has acceptable precision. Total atmospheric water vapour content is the key input parameter required by the SC algorithm. In this project, the MODIS water vapour product is employed to derive total atmospheric water vapour content. In this way, an operational Landsat LST product generation program was constructed by integration of MODIS and Landsat satellite imagery.     

Landsat8卫星影像大气校正及其植被指数与SEVI性能比较

遥感影像受大气的吸收散射以及地形起伏变化的影响,使得传感器接收到的辐射信号既包含了地物的信息,同时也包含了大气以及地形的信息。为了提高地表反射率的反演精度,需要去除遥感影像中大气和地形的影响。提出了一种基于查找表的Landsat8-OLI遥感影像的大气校正方法,该方法由6S辐射传输模型生成查找表,其中输入的参数包括大气水蒸汽含量、臭氧浓度和气溶胶光学厚度等MODIS大气参数产品。利用传统方法建立的大气参数查找表通常只考虑一部分因素,这对于以MODIS产品为输入参数的大气校正是不适用的。本文建立了一个包括大部分输入参数的高维大气校正查找表,对于Landsat-8 OLI传感器具有很高的通用性,通过进行光谱分析、与USGS地表反射率产品交叉验证等方式来验证模型的精度。验证结果表明该方法能有效地反演精确可靠的地表反射率。最后,采用目视解译、统计分析将校正结果与SEVI做对比分析,比较地形影响消减的效果。结果表明该模型与SEVI在地形消减的效果上作用相当。     

Landsat8卫星影像大气校正及其植被指数与SEVI性能比较

遥感影像受大气的吸收散射以及地形起伏变化的影响，使得传感器接收到的辐射信号既包含了地物的信息，同时也包含了大气以及地形的信息。为了提高地表反射率的反演精度，需要去除遥感影像中大气和地形的影响。提出了一种基于查找表的Landsat8-OLI遥感影像的大气校正方法，该方法由6S辐射传输模型生成查找表，其中输入的参数包括大气水蒸汽含量、臭氧浓度和气溶胶光学厚度等MODIS大气参数产品。利用传统方法建立的大气参数查找表通常只考虑一部分因素，这对于以MODIS产品为输入参数的大气校正是不适用的。本文建立了一个包括大部分输入参数的高维大气校正查找表，对于Landsat-8 OLI传感器具有很高的通用性，通过进行光谱分析、与USGS地表反射率产品交叉验证等方式来验证模型的精度。验证结果表明该方法能有效地反演精确可靠的地表反射率。最后，采用目视解译、统计分析将校正结果与SEVI做对比分析，比较地形影响消减的效果。结果表明该模型与SEVI在地形消减的效果上作用相当。     

Multi-temporal MODIS-Landsat data fusion for relative radiometric normalization, gap filling, and prediction of Landsat data

A semi-physical fusion approach that uses the MODIS BRDF/Albedo land surface characterization product and Landsat ETM+ data to predict ETM+ reflectance on the same, an antecedent, or subsequent date is presented. The method may be used for ETM+ cloud/cloud shadow and SLC-off gap filling and for relative radiometric normalization. It is demonstrated over three study sites, one in Africa and two in the U.S. (Oregon and Idaho) that were selected to encompass a range of land cover land use types and temporal variations in solar illumination, land cover, land use, and phenology. Specifically, the 30 m ETM+ spectral reflectance is predicted for a desired date as the product of observed ETM+ reflectance and the ratio of the 500 m surface reflectance modeled using the MODIS BRDF spectral model parameters and the sun-sensor geometry on the predicted and observed Landsat dates. The difference between the predicted and observed ETM+ reflectance (prediction residual) is compared with the difference between the ETM+ reflectance observed on the two dates (temporal residual) and with respect to the MODIS BRDF model parameter quality. For all three scenes, and all but the shortest wavelength band, the mean prediction residual is smaller than the mean temporal residual, by up to a factor of three. The accuracy is typically higher at ETM+ pixel locations where the MODIS BRDF model parameters are derived using the best quality inversions. The method is most accurate for the ETM+ near-infrared (NIR) band; mean NIR prediction residuals are 9%, 12% and 14% of the mean NIR scene reflectance of the African, Oregon and Idaho sites respectively. The developed fusion approach may be applied to any high spatial resolution satellite data, does not require any tuning parameters and so may be automated, is applied on a per-pixel basis and is unaffected by the presence of missing or contaminated neighboring Landsat pixels, accommodates for temporal variations due to surface changes (e.g., phenological, land cover/land use variations) observable at the 500 m MODIS BRDF/Albedo product resolution, and allows for future improvements through BRDF model refinement and error assessment.     

Validating MODIS land surface reflectance products using ground-measured reflectance spectra – a case study in semi-arid grassland in Inner Mongolia,China

Validation of Moderate-Resolution Imaging Spectroradiometer (MODIS) land surface reflectance products is important to effective utilization of such products for earth systems science. Ground-based measurements are normally utilized for such validation. However, the major scale mismatch between the ground ‘point’ measurement and MODIS resolution (500 m and 1 km) makes direct comparison infeasible over many land surface types. In this paper, an indirect comparison between ground ‘point’ measurements and MODIS land surface products via high-resolution remotely sensed imagery (Landsat Thematic Mapper/TM) was utilized in semi-arid grassland of Inner Mongolia in summer 2005, where ground measurements are relatively sparse in comparison with other locations around the world. Within the validation, the TM reflectance imagery was first calibrated by the ground ‘point’ measurements, and then aggregated to MODIS data resolution for determination of their accuracy. Besides common direct spectral band comparison of reflectance between TM and MODIS, empirical/indirect comparison between TM and MODIS was also implemented. Both types of validation showed that the absolute error of bidirectional reflectance from atmospheric correction (MOD09) is less than 9.4%, and for nadir bidirectional reflectance distribution function (BRDF)-adjusted reflectance (MOD43B4) it is less than 3.1%, in which the error of visible bands of two data sets is less than 1.35% and 0.95%, respectively. This validation will help improve the accuracy of MODIS products used in this area.     

Comparison between in situ and MODIS-derived spectral reflectances of snow and sea ice in the Amundsen Sea,Antarctica

The spectral albedo and directional reflectance of snow and sea ice were measured on sea ice of various types, including nilas, grey ice, pancake ice, multi-year pack ice, and land-fast ice in the Ross, Amundsen and Bellingshausen seas during a summer cruise in February through March 2000. Measurements were made using a spectroradiometer that has 512 channels in the visible and near-infrared (VNIR) region in which 16 of the 36 bands of the Moderate Resolution Imaging Spectroradiometer (MODIS) are covered. Directional reflectance is also retrieved from the MODIS radiometrically calibrated data (Level 1B) concurrently acquired from the first National Aeronautics and Space Administration (NASA) Earth Observing System (EOS) satellite, Terra. The locations of the ground ice stations are identified accurately on the MODIS images, and the spectral albedo and directional reflectance values at the 16 VNIR MODIS bands are extracted for those pixel locations. MODIS-derived reflectance is then corrected for the intervening atmosphere whose parameters are retrieved from the MODIS atmospheric profiles product (MOD07_L2) for the same granule. The corresponding spectral albedo and directional reflectance with the same viewing geometry as MODIS are derived from our ground-based spectroradiometer measurements. Because the footprint of the ground spectroradiometer is much smaller than the pixel sizes of MODIS images, the averaged spectral reflectance and albedo in the vicinity of each ice station are simulated for the corresponding MODIS pixel from the ground spectral measurements by weighting over different surface types (various ice types and open water). An accurate determination of ice concentration is important in deriving ground reflectance of a simulated pixel from in situ measurements. The best agreement between the in situ and MODIS measurements was found when the ground had 10/10 ice concentration (discrepancy range 0.2–11.69%, average 4.8%) or was oneice-type dominant (discrepancy range 0.8–16.9%, average 6.2%). The more homogeneous the ground surface and the less variable the ground topography, the more comparable between the in situ and satellite-derived reflectance is expected.     

Validation of the MODIS reflectance product under UK conditions

Surface reflectance obtained from remote-sensing data is the main input to almost all remote-sensing applications. The availability and special characteristics of Moderate Resolution Imaging Spectroradiometer (MODIS) products have led to their use worldwide. Validation of the MODIS reflectance product is then crucial to provid information on uncertainty in the reflectance data, and in other MODIS products and in the applied surface–atmosphere models. Compact Airborne Spectrographic Imager (CASI) and Système Pour l'Observation de la Terre (SPOT) data, collected during the Network for Calibration and Validation in Earth Observation (NCAVEO) 2006 Field Campaign, were applied to validate daily MODIS reflectance data over a site in the southern UK. The difference in the view geometry of at-nadir CASI and SPOT data and off-nadir MODIS data was dealt with using a semi-empirical bidirectional reflectance distribution function (BRDF) model. The validation results showed that for our particular study site, the absolute errors in the MODIS reflectance product were too large to allow the albedo data to be used directly in climate models. The errors were mainly related to the uncertainties in the MODIS atmospheric variables, the BRDF model, and undetected clouds and cloud shadows. More generally, the study highlights the extreme difficulty of achieving pixel-level validation of coarse spatial resolution satellite sensor data in an environment in which the atmosphere is constantly changing, and in which the landscape is characterized by high space–time heterogeneity.     

An Albedo Downscaling Method based on Stratified Linear Regression

In recent years,remote sensing technology has been widely used in the field of surface energy,thermal environment,and climate research in small and medium-sized regions.The demand for high resolution,high precision albedo products has increased.In view of this,a physical-based downscaling method is proposed for efficiently and accurately generating high-resolution albedo results.First,under the Lambert hypothesis,the primary albedo of Landsat 8 could be obtained based on Landsat 8 reflectance at 30m resolution.On the 500 m scale,it is found that after classification,the primary albedo of Landsat 8 has a better correlation with the broad-band albedo of MODIS MCD43A3 product.Therefore,a linear regression function based on surface classification is established to integrate the high-resolution primary albedo of Landsat 8 with high-precision MCD43A3 albedo to obtain downscaled albedo.Compared to MODIS albedo,the downscaled albedo provided more rich details.Verification experience based on SURFRAD observation data shows that Bias of the albedo downscaling is 0.01,and standard deviation is 0.012,which has good adaptability for different surface categories.It shows that the algorithm has great application value for producing high-resolution albedo products.     

Radiometric correction of multi-temporal Landsat data for characterization of early successional forest patterns in western Oregon

Detecting and characterizing continuous changes in early forest succession using multi-temporal satellite imagery requires atmospheric correction procedures that are both operationally reliable, and that result in comparable units (e.g., surface reflectance). This paper presents a comparison of five atmospheric correction methods (2 relative, 3 absolute) used to correct a nearly continuous 20-year Landsat TM/ETM+ image data set (19-images) covering western Oregon (path/row 46/29). In theory, full absolute correction of individual images in a time-series should effectively minimize atmospheric effects resulting in a series of images that appears more similar in spectral response than the same set of uncorrected images. Contradicting this theory, evidence is presented that demonstrates how absolute correction methods such as Second Simulation of the Satellite Signal in the Solar Spectrum (6 s), Modified Dense Dark Vegetation (MDDV), and Dark Object Subtraction (DOS) actually make images in a time-series somewhat less spectrally similar to one another. Since the development of meaningful spectral reflectance trajectories is more dependant on consistent measurement of surface reflectance rather than on accurate estimation of true surface reflectance, correction using image pairs is also tested. The relative methods tested are variants of an approach referred to as “absolute-normalization”, which matches images in a time-series to an atmospherically corrected reference image using pseudo-invariant features and reduced major axis (RMA) regression. An advantage of “absolute-normalization” is that all images in the time-series are converted to units of surface reflectance while simultaneously being corrected for atmospheric effects. Of the two relative correction methods used for “absolute-normalization”, the first employed an automated ordination algorithm called multivariate alteration detection (MAD) to statistically locate pseudo-invariant pixels between each subject and reference image, while the second used analyst selected pseudo-invariant features (PIF) common to the entire image set. Overall, relative correction employed in the “absolute-normalization” context produced the most consistent temporal reflectance response, with the automated MAD algorithm performing equally as well as the handpicked PIFs. Although both relative methods performed nearly equally in terms of observed errors, several reasons emerged for preferring the MAD algorithm. The paper concludes by demonstrating how “absolute-normalization” improves (i.e., reduces scatter in) spectral reflectance trajectory models used for characterizing patterns of early forest succession.     

黑河及汉江流域MODIS叶面积指数产品质量评价

叶面积指数（LAI）是MODIS地面队伍生产的一系列标准产品之一，对其进行独立的质量评价有助于用户了解数据的适用性。本文用近同时相的Landsat影像及野外实测LAI数据获得了黑河及汉江两个研究区高分辨率的Landsat LAl分布图。基于此，对MODIS LAI数据进行了质量评价，评价指标包括统计特征和空间特征。分析结果表明，就统计特征而言，MODIS LAI数据值一般低于Landsat的LAI值，在植被覆盖较好的汉江区低估约10％，在植被覆盖稀疏的黑河区，LAI值低估达58％；就空间特征而言，两个研究区的结果都表明MODIS LAI数据无法很好地体现植被空间分布信息，在黑河区存在大量低槽被覆盖像元被归类为非植被覆盖区的情况。     

The 250 m global land cover change product from the Moderate Resolution Imaging Spectroradiometer of NASA's Earth Observing System

One of the products to be derived from data acquired by the Moderate Resolution Imaging Spectroradiometer (MODIS) onboard the National Aeronautics and Space Administration (NASA) Earth Observing System (EOS) will identify locations where land cover changes attributable to human activities are occurring. The product aims to serve as an alarm where rapid land cover conversion can subsequently be analysed with higher resolution sensors such as Landsat 7. This paper describes the production procedure and change detection algorithms for the 250 m MODIS land cover change product. Multiple change detection algorithms, including three spectral methods and two texture methods, are utilized to generate the product from the two 250 m MODIS bands. The change detection methods are implemented with look-up tables (LUTs) which are initially generated using data from the Advanced Very High Resolution Radiometer (AVHRR) and Landsat Thematic Mapper (TM) until MODIS data become available. The results from the five methods are combined to improve confidence in the product. We test the performance of each method against several sets of simulated MODIS data derived from Landsat TM image pairs. The test data represent tropical deforestation, agricultural expansion and urbanization. The commission errors of the five methods and the combination are approximately 10%, with reasonable omission errors less than 25%.     

Use of Landsat TM and EOS MODIS imaging technologies for estimation of winter wheat yield in the North China Plain

This study focuses on the methodologies of winter wheat yield prediction based on Land Satellite Thematic Map (TM) and Earth Observation System Moderate Resolution Imaging Spectroradiometer (MODIS) imaging technologies in the North China Plain. Routine field measurements were initiated during the periods when the Landsat satellite passed over the study region. Five Landsat TM images were acquired. Wheat yields of the experimental sites were recorded after harvest. Spectral vegetation indices were calculated from TM and MODIS images. The correlation analysis among wheat yield and spectral parameters revealed that TM renormalized difference water index (RDWI) and MODIS near-infrared reflectance had the highest correlation with yield at grain-filling stages. The models from the best-fitting method were used to estimate wheat yield based on TM and MODIS data. The average relative error of the root mean square error (RMSE) of the predicted yield was smaller from TM than from MODIS.     

A consistency analysis of surface reflectance and leaf area index retrieval from overlapping clear-sky Landsat ETM+ imagery

In this study, the consistency of systematic retrievals of surface reflectance and leaf area index was assessed using overlap regions in adjacent Landsat Enhanced Thematic Mapper-Plus (ETM+) scenes. Adjacent scenes were acquired within 7-25 days apart to minimize variations in the land surface reflectance between acquisition dates. Each Landsat ETM+ scene was independently geo-referenced and atmospherically corrected using a variety of standard approaches. Leaf area index (LAI) models were then applied to the surface reflectance data and the difference in LAI between overlapping scenes was evaluated. The results from this analysis show that systematic LAI retrieval from Landsat ETM+ imagery using a baseline atmospheric correction approach that assumes a constant aerosol optical depth equal to 0.06 is consistent to within ±0.61 LAI units. The average absolute difference in LAI retrieval over all 10 image pairs was 26% for a mean LAI of 2.05 and the maximum absolute difference over any one pair was 61% for a mean LAI of 1.13. When no atmospheric correction was performed on the data, the consistency in LAI retrieval was improved by 1%. When a scene-based dense, dark vegetation atmospheric correction algorithm was used, the LAI retrieval differences increased to 28% for a mean LAI of 2.32. This implies that a scene-based atmospheric correction procedure may improve the absolute accuracy of LAI retrieval without having a major impact on retrieval consistency. Such consistency trials provide insight into the current limits concerning surface reflectance and LAI retrieval from fine spatial resolution remote sensing imagery with respect to the variability in clear-sky atmospheric conditions.     

Comparison of atmospheric correction algorithms for TM image in inland waters

In order to extract quantitative water‐leaving information from the Thematic Mapper (TM) image accurately in inland waters, atmospheric correction is a necessary step. Based on former researchers' results, the paper presents two atmospheric correction algorithms based on meteorological data (MD) and on Moderate Resolution Imaging Spectroradiometer (MODIS) Vicarious Calibration (MVC) for TM image in inland waters according to the theory of radiative transfer. Studying Taihu lake, China, in this paper we derived water remote sensing reflectance from a TM image of 26 July 2004 by these two atmospheric correction algorithms and we compare the results with that of dark object subtraction (DOS) and 6S code. The results show that the effect of atmospheric correction based on meteorological data and MODIS Vicarious Calibration is much better than that of DOS and 6S code. Although the MD is more accurate, MVC may be an ideal choice for TM images in inland water because TERRA MODIS images can be acquired easily than collecting meteorological data at the time of satellites passing over.     

Monitoring on-orbit calibration stability of the Terra MODIS and Landsat 7 ETM+ sensors using pseudo-invariant test sites

The ability to detect and quantify changes in the Earth's environment depends on sensors that can provide calibrated, consistent measurements of the Earth's surface features through time. A critical step in this process is to put image data from different sensors onto a common radiometric scale. This work focuses on monitoring the long-term on-orbit calibration stability of the Terra Moderate Resolution Imaging Spectroradiometer (MODIS) and the Landsat 7 (L7) Enhanced Thematic Mapper Plus (ETM+) sensors using the Committee on Earth Observation Satellites (CEOS) reference standard pseudo-invariant test sites (Libya 4, Mauritania 1/2, Algeria 3, Libya 1, and Algeria 5). These sites have been frequently used as radiometric targets because of their relatively stable surface conditions temporally. This study was performed using all cloud-free calibrated images from the Terra MODIS and the L7 ETM+ sensors, acquired from launch to December 2008. Homogeneous regions of interest (ROI) were selected in the calibrated images and the mean target statistics were derived from sensor measurements in terms of top-of-atmosphere (TOA) reflectance. For each band pair, a set of fitted coefficients (slope and offset) is provided to monitor the long-term stability over very stable pseudo-invariant test sites. The average percent differences in intercept from the long-term trends obtained from the ETM + TOA reflectance estimates relative to the MODIS for all the CEOS reference standard test sites range from 2.5% to 15%. This gives an estimate of the collective differences due to the Relative Spectral Response (RSR) characteristics of each sensor, bi-directional reflectance distribution function (BRDF), spectral signature of the ground target, and atmospheric composition. The lifetime TOA reflectance trends from both sensors over 10 years are extremely stable, changing by no more than 0.4% per year in its TOA reflectance over the CEOS reference standard test sites.     

Spectral assessment of new ASTER SWIR surface reflectance data products for spectroscopic mapping of rocks and minerals

ASTER reflectance spectra from Cuprite, Nevada, and Mountain Pass, California, were compared to spectra of field samples and to ASTER-resampled AVIRIS reflectance data to determine spectral accuracy and spectroscopic mapping potential of two new ASTER SWIR reflectance datasets: RefL1b and AST_07XT. RefL1b is a new reflectance dataset produced for this study using ASTER Level 1B data, crosstalk correction, radiance correction factors, and concurrently acquired level 2 MODIS water vapor data. The AST_07XT data product, available from EDC and ERSDAC, incorporates crosstalk correction and non-concurrently acquired MODIS water vapor data for atmospheric correction. Spectral accuracy was determined using difference values which were compiled from ASTER band 5/6 and 9/8 ratios of AST_07XT or RefL1b data subtracted from similar ratios calculated for field sample and AVIRIS reflectance data. In addition, Spectral Analyst, a statistical program that utilizes a Spectral Feature Fitting algorithm, was used to quantitatively assess spectral accuracy of AST_07XT and RefL1b data.Spectral Analyst matched more minerals correctly and had higher scores for the RefL1b data than for AST_07XT data. The radiance correction factors used in the RefL1b data corrected a low band 5 reflectance anomaly observed in the AST_07XT and AST_07 data but also produced anomalously high band 5 reflectance in RefL1b spectra with strong band 5 absorption for minerals, such as alunite. Thus, the band 5 anomaly seen in the RefL1b data cannot be corrected using additional gain adjustments. In addition, the use of concurrent MODIS water vapor data in the atmospheric correction of the RefL1b data produced datasets that had lower band 9 reflectance anomalies than the AST_07XT data. Although assessment of spectral data suggests that RefL1b data are more consistent and spectrally more correct than AST_07XT data, the Spectral Analyst results indicate that spectral discrimination between some minerals, such as alunite and kaolinite, are still not possible unless additional spectral calibration using site specific spectral data are performed.     

Quantitative relationships of near-surface spectra to Landsat radiometric data

The aim of our research has been to determine the quantitative relationship between the surface spectral character of a variety of geologic terrains and that sensed by the Landsat multispectral scanner. A spectral sampling and measurement program was conducted to accurately characterize the surface spectral reflectance of the Landsat resolution element and, for the first time, to establish statistically the degree of sampling required for a variety of natural terrains. Results from the study showed that for typical homogeneous and moderately heterogeneous terrains, the number of samples required to estimate the mean reflectance of a pixel is small. Only 9–20 samples are required to be within 2% reflectance at the 95% probability level. Coincident field measurements and satellite observations were used to test the equivalency and correlation of the reflectance data. Before the Landsat data could be compared with the surface measurements the satellite brightness values must be converted to absolute radiometric units, and corrected for atmospheric attenuation and scattering. A conversion method using a standard/target comparison, which indirectly compensated for atmospheric attenuation and scattering, produced a Landsat equivalent reflectance that exhibited a root-mean-square error of ± 4% reflectance, when compared with the surface measured value at 12 test sites. Although the equivalence of the surface and satellite data cannot be shown to be better than 4% reflectance, statistical study indicates that the surface and satellite data are highly correlated within defined contrast constraints. However, this correlation is present only after the satellite brightness values are corrected for between band gain differences and compensation is made for atmospheric attenuation and scattering.     

基于Landsat-5 TM的洪河湿地地表温度估算方法对比研究

利用洪河湿地2008年5月15日过境的Landsat/TM图像和实测地面数据以及MODIS 地表发射率数据,分别运用大气辐射传输模型、覃志豪的单窗算法和Jimenez\|Munoz & Sobrino 的单波段算法估算洪河湿地的地表温度,并且对比了大气校正前后的NDVI、LSE以及各种算法估算地表温度的差异。分析估算结果表明,覃志豪的单窗算法与实测地面数据估算结果非常一致。指出在没有实时探空数据的情况下,应用只有一个热红外通道的Landsat/TM数据源,采用覃志豪的单窗算法估算的精度是可以接受的。     

Simulation of hyperspectral and directional radiance images using coupled biophysical and atmospheric radiative transfer models

Future remote sensing satellite missions exploring the earth will feature advanced hyperspectral and directional optical imaging instruments. Given the complex nature of the data to be expected from these missions, a thorough preparation for them is essential and this can be accomplished by realistic simulation of the imagery data, years before the actual launch. Based on given spectral and directional capabilities of the instrument, and in combination with biophysical land surface properties obtained from existing imagery, the spectral and directional responses of several types of vegetation and bare soil have been simulated pixel by pixel using the radiative transfer models PROSPECT (for hyperspectral leaf reflectance and transmittance), GeoSAIL (for two-layer canopy bidirectional spectral reflectance), and MODTRAN4 (for atmospheric hyperspectral and directional effects). In this way, one obtains realistically simulated hyperspectral and directional top-of-atmosphere spectral radiance images, with all major effects included, such as heterogeneity of the landscape, non-Lambertian reflectance of the land surface, the atmospheric adjacency effect, and the limited spatial resolution of the instrument. The output of the image simulations can be used to demonstrate the capabilities of future earth observation missions. In addition, instrument specifications and image acquisition strategies might be optimized on the basis of simulated image analysis results, and new advanced data assimilation procedures could be validated with realistic inputs under controlled circumstances. This paper describes the applied methodology, the study area with the input images, the set-up of the actual image simulations, and discusses the final results obtained.