卫星热红外图像与震兆异常——澜沧地震前热红外图像的启示

本文基于1988年11月6日21时澜沧7.6级地震前的气象台站地面观测温度与相近时刻卫星热红外资料,分析了热红外图像与构造、岩性、地貌及气象等的关系,阐述了震兆增温异常与正常情况下的区别,初步总结了利用热红外图像的时、空动态变化特征捕捉地震临震前兆的方法和体会。     

Modelling and Remote Sensing of Land Surface Temperature in Turkey

This study introduces artificial neural networks (ANNs) for the estimation of land surface temperature (LST) using meteorological and geographical data in Turkey (26?C45°E and 36?C42°N). A generalized regression neural network (GRNN) was used in the network. In order to train the neural network, meteorological and geographical data for the period from January 2002 to December 2002 for 10 stations (Adana, Afyon, Ankara, Eski?ehir, ?stanbul, ?zmir, Konya, Malatya, Rize, Sivas) spread over Turkey were used as training (six stations) and testing (four stations) data. Latitude, longitude, elevation and mean air temperature are used in the input layer of the network. Land surface temperature is the output. However, land surface temperature has been estimated as monthly mean by using NOAA-AVHRR satellite data in the thermal range over 10 stations in Turkey. The RMSE between the estimated and ground values for monthly mean with ANN temperature(LST_ANN) and Becker and Li temperature(LST_B-L) method values have been found as 0.077?K and 0.091?K (training stations), 0.045?K and 0.003?K (testing stations), respectively.     

Accuracy assessment of high resolution satellite imagery orientation by leave-one-out method

Interest in high-resolution satellite imagery (HRSI) is spreading in several application fields, at both scientific and commercial levels. Fundamental and critical goals for the geometric use of this kind of imagery are their orientation and orthorectification, processes able to georeference the imagery and correct the geometric deformations they undergo during acquisition. In order to exploit the actual potentialities of orthorectified imagery in Geomatics applications, the definition of a methodology to assess the spatial accuracy achievable from oriented imagery is a crucial topic.In this paper we want to propose a new method for accuracy assessment based on the Leave-One-Out Cross-Validation (LOOCV), a model validation method already applied in different fields such as machine learning, bioinformatics and generally in any other field requiring an evaluation of the performance of a learning algorithm (e.g. in geostatistics), but never applied to HRSI orientation accuracy assessment.The proposed method exhibits interesting features which are able to overcome the most remarkable drawbacks involved by the commonly used method (Hold-Out Validation — HOV), based on the partitioning of the known ground points in two sets: the first is used in the orientation–orthorectification model (GCPs — Ground Control Points) and the second is used to validate the model itself (CPs — Check Points). In fact the HOV is generally not reliable and it is not applicable when a low number of ground points is available.To test the proposed method we implemented a new routine that performs the LOOCV in the software SISAR, developed by the Geodesy and Geomatics Team at the Sapienza University of Rome to perform the rigorous orientation of HRSI; this routine was tested on some EROS-A and QuickBird images. Moreover, these images were also oriented using the world recognized commercial software OrthoEngine v. 10 (included in the Geomatica suite by PCI), manually performing the LOOCV since only the HOV is implemented.The software comparison guaranteed about the overall correctness and good performances of the SISAR model, whereas the results showed the good features of the LOOCV method.     

Calibration facility for airborne imaging spectrometers

A new facility designed to perform calibration measurements of airborne imaging spectrometers was established at the German Aerospace Center (DLR) in Oberpfaffenhofen. This Calibration Home Base (CHB) is optimized to characterize radiometrically, spectrally, and geometrically the APEX (Airborne Prism Experiment) imaging spectrometer, which is currently being developed under the authority of the European Space Agency (ESA). It however can be used for other optical sensors as well. Computer control of major laboratory equipment allows automation of time consuming measurements. In APEX configuration (wavelength range: 380 to 2500 nm, instantaneous field of view: 0.48 mrad, field of view: ±14) spectral measurements can be performed to a wavelength uncertainty of ±0.15 nm, geometric measurements at increments of 0.0017 mrad across track and 0.0076 mrad along track, and radiometric measurements to an uncertainty of ±3% relative to national standard. The CHB can be adapted to similar sensors (including those with thermal infrared detectors) by exchanging the monochromator’s lamp, the gratings and the filters, and by adjusting the distance between the sensor and folding mirror.     

Evaluation of approaches for AWiFS multi-date registration

AWiFS onboard IRS-P6 belongs to the category of high-repetivity sensors based on large swath, but with ground trace based on narrow-swath sensor (LISS-III). This is useful for cloud removal as well as vegetation phenology studies. Such multi-date analysis has a prerequisite of accurate multi-date registration. This study investigates the accuracy of multi-date registration over a mixed plain and hilly terrain in northern India (29–31°N latitude and 77.5–79.5°E longitude; 200–4000 m.a.s.l.). Simple polynomial rectification, multi-date registration using ortho-correction technique on standard product (level-2) and radiometric product (level-1) as a function of number of ground control points (GCPs) and external Digital Elevation Model (DEM) were investigated. The results indicated that ortho-rectification on level-1 product provided better accuracy in comparison to simple rectification and ortho-rectification on level-2 product.     

Determination of the bulk temperature from NOAA/AVHRR satellite data in a midlatitude lake

The spatial and temporal variability of the bulk temperature gives important insights into biological and hydrodynamic processes. However, standard algorithms for satellite data only provide information of the surface temperature. The comparison of current and new split-window coefficients applied on NOAA-14/AVHRR brightness temperatures of Lake Constance showed that a regional adaption was most promising. To derive the bulk temperature information, a priori progression from a weather station was included into the AVHRR analysis. Among the weather is data, the mean temperature of the three preceding days and the day of the year were the most relevant additional information. By a multiple regression approach the bulk temperature in the upper 4 m of Lake Constance could be determined with an accuracy of ±1.20 °C. The training of a neural network improved the determination of the bulk temperature to ±1.04 °C.An extended field campaign demonstrated that the algorithm is also applicable to other sensors with the same spectral band settings (in this case NOAA-16/AVHRR) with an acceptable error and that it is equally accurate over the entire lake.     

Estimation and validation of snow surface temperature using modis data for snow-avalanche studies in NW-Himalaya

Snow avalanche studies require different snow-meteorological parameters for avalanche forecasting. Snow surface temperature is one of the major parameters, which is responsible for the evolution of snow pack characteristics. In the present paper, the snow surface temperature was estimated using TERRA satellite based — Moderate resolution imaging spectroradiometer (MODIS) sensor for NW-Himalayas. Ground data observed by automatic weather stations (AWS) was used to calibrate the brightness temperature obtained by MODIS thermal bands data into the actual snow surface temperature data through regression analysis. A split window technique has been implemented for the estimation of snow surface temperature. The multi-date satellite derived snow surface temperature was validated with ground data of winter 2004–05 and 2005–06 collected at various observation stations located in different ranges of NW-Himalaya. Good correlations were observed for Upper Himalaya (0.98, 0.98), Middle Himalaya (0.92, 0.96) and Lower Himalaya (0.88, 0.82) for 2004–05 and 2005–06 winter respectively. Further, estimated snow surface temperature was also verified with snow-cover information collected by manned observatories and area delineated by thematic maps of snow surface temperature was validated with the different snow climatic zones of NW-Himalaya.     

Bidirectional reflectance spectrometry of gravel at the Sjökulla test field

The Sjökulla test site is used for testing and calibrating aerial images. The permanent test field is made of four types of gravel (dark gabbro, grey granite, red granite, white limestone) in two sizes (diameters 8–16 mm and 4–8 mm) set in various patterns. The bidirectional reflection properties of the targets together with their temporal changes must be known in order to carry out radiometric and spectral evaluation and calibration.The bidirectional reflectance distribution functions (BRF) of the gravel have been measured several times in the test fields using portable field goniospectrometers belonging Finnish Geodetic Institute (FGI), and once using the European Goniometic Facility (EGO) of the Joint Research Centre (JRC) at Ispra, Italy. Detailed BRFs have been obtained, showing features typical to particulate media, e.g. a small bowl shape, strong backscattering, and smooth wavelength dependence. Temporal range measurements over several years show that the black gabbro and red granite are fairly stable, while the grey granite has changed somewhat over the years and the white limestone has experienced dramatic darkening effects, requiring action to be taken.The measured BRF data have increased the usability of the test field considerably. The results are also useful in the development and validation of scattering models for particulate media. The site has proved to be a good test bench for goniospectrometric instruments, too.     

Remote sensing as applied to mineral exploration in the Chalkidiki peninsula, northern Greece

Landsat (TM) imagery studies combined with photogeology, heliborne EM-radiometric interpretations and field data provided the fracture pattern in the eastern Chalkidiki peninsula, northern Greece. Processing of the respective data revealed the following four main directions of lineaments trending 070°–080°, 120°–130°, N-S and E-W. This fracture pattern is persistent in all geological units including the Tertiary basins and may represent reactivated older or even Tertiary fractures superimposed on all units. The control and the distribution of the major mineralization such as Pb---Zn (Au---Ag) massive sulphides, Mn (Fe) oxides, and porphyry Cu (Au) are of similar orientation to the deduced fracture pattern.These findings support the remote sensing is an important mineral exploration tool which can become quite effective when combined with ore deposit type, geochemical and geophysical data in partly explored areas.     

Application of remote sensing technology for Impact Assessment of Watershed Development Programme

The current study was taken up to investigate the utility of remote sensing and GIS tools for evaluation of Integrated Wasteland Development Programme (IWDP) implemented during 1997–2001 in Katangidda Nala watershed, Chincholi taluk, Gulbarga district, Karnataka. The study was carried out using IRS 1C, LISS III data of December 11, 1997 (pre-treatment) and November 15, 2002 (post-treatment) covering the watershed to assess the changes in land use / land cover and biomass that have changed over a period of five years (1997–2002). The images were classified into different land use/land cover categories using supervised classification by maximum likelihood algorithm. They were also classified into different biomass levels using Normalized Difference Vegetation Index (NDVI) approach. The results indicated that the area under agriculture crops and forest land were increased by 671 ha (5.7%) and 1,414 ha (11.94%) respectively. This is due to the fact that parts of wastelands and fallow lands were brought into cultivation. This increase in the area may be attributed to better utilization of surface and ground waters, adoption of soil and water conservation practices and changes in cropping pattern. The area under waste lands and fallow lands decreased by 1,667 ha (14.07%) and 467 ha (3.94%), respectively. The vegetation vigour of the area was classified into three classes using NDVI. Substantial increase in the area under high and low biomass levels was observed (502 ha and 19 ha respectively). The benefit-cost analysis indicates that the use of remote sensing and GIS was 2.2 times cheaper than the conventional methods. Thus, the repetitive coverage of the satellite data provides an excellent opportunity to monitor the land resources and evaluate the land cover changes through comparison of images for the watershed at different periods.     

The ADS40 Vaihingen/Enz geometric performance test

This paper presents the main results of a comprehensive ADS40 performance analysis conducted at the Vaihingen/Enz test field. As such it represents one example of an independent in-flight performance study for one of the new and commercially operational digital airborne camera systems. Based on a large number of well coordinated and defined object points, which served as independent check points, the absolute geometric accuracy of the ADS40 from true operational data has been verified. Empirical analysis of data from flying heights ranging from 1500 m to 4000 m proved the ADS40 geometric accuracy to be in the range of 1–2 μm at image scale for horizontal coordinates and 0.03–0.05‰ of the flying height for vertical components. This is fully within specification for airborne imaging.     

Parameterization of air temperature in high temporal and spatial resolution from a combination of the SEVIRI and MODIS instruments

Some applications, e.g. from traffic or energy management, require air temperature data in high spatial and temporal resolution at two metres height above the ground (T_2m), sometimes in near-real-time. Thus, a parameterization based on boundary layer physical principles was developed that determines the air temperature from remote sensing data (SEVIRI data aboard the MSG and MODIS data aboard Terra and Aqua satellites). The method consists of two parts. First, a downscaling procedure from the SEVIRI pixel resolution of several kilometres to a one kilometre spatial resolution is performed using a regression analysis between the land surface temperature (LST) and the normalized differential vegetation index (NDVI) acquired by the MODIS instrument. Second, the lapse rate between the LST and T_2m is removed using an empirical parameterization that requires albedo, down-welling surface short-wave flux, relief characteristics and NDVI data. The method was successfully tested for Slovenia, the French region Franche-Comté and southern Germany for the period from May to December 2005, indicating that the parameterization is valid for Central Europe. This parameterization results in a root mean square deviation RMSD of 2.0 K during the daytime with a bias of −0.01 K and a correlation coefficient of 0.95. This is promising, especially considering the high temporal (30 min) and spatial resolution (1000 m) of the results.     

Study of thermal IR phenomena associated with 27 February 2010 Chile Mw 8.8 earthquake using MODIS data

Stresses acting on fault systems before major earthquakes can produce thermal anomalies; these abnormalities can be observed using multi-sensor satellite data. Moderate resolution spectro-radiometer on board the terra and aqua satellites can provide thermal infrared (TIR) imaging data for land and ocean. These TIR data have recorded short-lived thermal anomaly prior to major earthquakes. It is suggested by others that these electromagnetic (EM) phenomena relate to stress build up before earthquakes. The objective of this study is to find an association between spatial extent and temporal evolution of thermal anomalies and known major earthquakes near the boundary of Nazca plate and South American plate. Our approach is to map the TIR transient fields from polar orbiting satellites and analysing those data using time series temperature plots to detect the abnormal thermal trends before the earthquake. This study concentrated on marine earthquakes to detect the changes in both land and ocean before seismic activity.     

2002年秋季山东省干旱遥感监测分析

利用极轨气象卫星遥感资料对2002年山东省秋季干旱进行监测和服务,其监测方法主要采用热惯量法。为了提高土壤水分模式的计算精度,将全省分为4个区域,对利用卫星资料反演的地表温度日较差进行植被指数等环境因素订正。在地理信息系统支持下,比较准确地计算出农田受旱面积,效果良好。     

An inter-comparison of soil moisture data products from satellite remote sensing and a land surface model

Significant advances have been achieved in generating soil moisture (SM) products from satellite remote sensing and/or land surface modeling with reasonably good accuracy in recent years. However, the discrepancies among the different SM data products can be considerably large, which hampers their usage in various applications. The bias of one SM product from another is well recognized in the literature. Bias estimation and spatial correction methods have been documented for assimilating satellite SM product into land surface and hydrologic models. Nevertheless, understanding the characteristics of each of these SM data products is required for many applications where the most accurate data products are desirable. This study inter-compares five SM data products from three different sources with each other, and evaluates them against in situ SM measurements over 14-year period from 2000 to 2013. Specifically, three microwave (MW) satellite based data sets provided by ESA's Climate Change Initiative (CCI) (CCI-merged, -active and -passive products), one thermal infrared (TIR) satellite based product (ALEXI), and the Noah land surface model (LSM) simulations. The in-situ SM measurements are collected from the North American Soil Moisture Database (NASMD), which involves more than 600 ground sites from a variety of networks. They are used to evaluate the accuracies of these five SM data products. In general, each of the five SM products is capable of capturing the dry/wet patterns over the study period. However, the absolute SM values among the five products vary significantly. SM simulations from Noah LSM are more stable relative to the satellite-based products. All TIR and MW satellite based products are relatively noisier than the Noah LSM simulations. Even though MW satellite based SM retrievals have been predominantly used in the past years, SM retrievals of the ALEXI model based on TIR satellite observations demonstrate skills equivalent to all the MW satellite retrievals and even slightly better over certain regions. Compared to the individual active and passive MW products, the merged CCI product exhibits higher anomaly correlation with both Noah LSM simulations and in-situ SM measurements.     

多时相遥感数据对城市热岛效应的研究——以西安市为例

热岛效应是影响我国城市环境的重要因素,本文基于多时相的Landsat 7 ETM+和Landsat 8 TIRS影像,对西安市2002年和2013年的土地进行分类以及地表温度反演,在此基础上将土地分类影像图和温度分割影像图进行叠加,分析了不同的土地利用类型温度差异、下垫面的空间结构与热场分布的关系以及热场的分布变化与城市变迁的关系。结果表明:不同的土地利用类型的物理属性和生物学属性差异是造成城市热岛效应的内在原因,下垫面空间结构的差异是造成城市热岛效应的外在原因,城市热场的动态变化方向与城市的变迁方向具有一致性。     

Identifying natural and anthropogenically-induced geohazards from satellite ground motion and geospatial data: Stoke-on-Trent,UK

Determining the location and nature of hazardous ground motion resulting from natural and anthropogenic processes such as landslides, tectonic movement and mining is essential for hazard mitigation and sustainable resource use. Ground motion estimates from satellite ERS-1/2 persistent scatterer interferometry (PSI) were combined with geospatial data to identify areas of observed geohazards in Stoke-on-Trent, UK. This investigation was performed within the framework of the EC FP7-SPACE PanGeo project which aimed to provide free and open access to geohazard information for 52 urban areas across Europe. Geohazards identified within the city of Stoke-on-Trent and neighbouring rural areas are presented here alongside an examination of the PanGeo methodology.A total of 14 areas experiencing ground instability caused by natural and anthropogenic processes have been defined, covering 122.35 km². These are attributed to a range of geohazards, including landslides, ground dissolution, made ground and mining activities. The dominant geohazard (by area) is ground movement caused by post-mining groundwater recharge and mining-related subsidence (93.19% of total geohazard area), followed by landsliding (5.81%). Observed ground motions along the satellite line-of-sight reach maxima of +35.23 mm/yr and −22.57 mm/yr. A combination of uplift, subsidence and downslope movement is displayed.‘Construction sites’ and ‘continuous urban fabric’ (European Urban Atlas land use types) form the land uses most affected (by area) by ground motion and ‘discontinuous very low density urban fabric’ the least. Areas of ‘continuous urban fabric’ also show the highest average velocity towards the satellite (5.08 mm/yr) and the highest PS densities (1262.92 points/km²) along with one of the lowest standard deviations. Rural land uses tend to result in lower PS densities and higher standard deviations, a consequence of fewer suitable reflectors in these regions. PSI is also limited in its ability to identify especially rapid ground motion. As a consequence the supporting geospatial data proved especially useful for the identification of landslides and some areas of ground dissolution. The mapped areas of instability are also compared with modelled potential geohazards (the BGS GeoSure dataset).     

Determination of land degradation causes in Tongyu County,Northeast China via land cover change detection

Tongyu County in Northeast China is highly prone to land degradation due to its fragile physical settings characterized by a flat topography, a semi-arid climate, and a shallow groundwater table. This study aims to determine the causes of land degradation through detecting the long-term trend of land cover changes. Degraded lands were mapped from satellite images recorded in 1992 and 2002. These land cover maps revealed that the area subject to land degradation in the form of soil salinization, waterlogging and desertification increased from 2400 to 4214 km², in sharp contrast to most severely degraded land that decreased by 122.5 km². Newly degraded land stems from productive farmland (263 km²), harvested farmland (551 km²), and grassland (468 km²). Therefore, the worsened degradation situation is attributed to excessive reclamation of grassland for farming, over cultivation, overgrazing, and deforestation. Mechanical, biological, ecological and engineering means should be adopted to rehabilitate the degraded land.     

Estimation of regional evapotranspiration over the North China Plain using geostationary satellite data

Data from the first operational Chinese geostationary satellite Fengyun-2C (FY-2C) satellite are applied in combination with Moderate Resolution Imaging Spectroradiometer (MODIS) satellite products for the assessment of regional evapotranspiration over the North China Plain. The approach is based on the improved triangle method, where the temperature–vegetation index space includes thermal inertia. Two thermal infrared channels from FY-2C are used to estimate surface temperature (Ts) based on a split window algorithm originally proposed for the MSG-SEVIRI sensor. Subsequently the high temporal resolution of FY-2C data is exploited to give the morning rise in Ts. Combined with the 16 days composite MODIS vegetation indices product (MOD13) at a spatial resolution of 5 km, evaporative fraction (EF) is estimated by interpolation in the ΔTs–NDVI triangular-shaped scatter space. Finally, regional actual evapotranspiration (ET) is derived from the evaporative fraction and available energy estimated from MODIS surface albedo products MCD43. Spatial variations of estimated surface variables (Ts, EF and ET) corresponded well to land cover patterns and farmland management practices. Estimated ET and EF also compared well to lysimeter data collected for the period June 2005–September 2007. The improved triangle method was also applied to MODIS products for comparison. Estimates based on FY-2C products proved to provide slightly better results than those based on MODIS products. The consistency of the estimated spatial variation with other spatial data supports the use of FY-2C data for ET estimation using the improved triangle method. Of particular value is the high temporal frequency of image acquisitions from FY-2C which improves the likelihood of obtaining cloud free image acquisitions as compared to polar orbiting sensors like MODIS.