三维地震数据离散光滑插值的共轭梯度法

针对三维地震数据插值,提出采用Laplacian算子进行光滑约束的插值方法,并借鉴Mallet研究的离散光滑插值思路,采用预条件共轭梯度法,直接生成网格节点上的值,从而回避寻求满足插值方程的函数. 为了实现其中Laplacian算子的快速求逆,文中引入Claerbout螺旋坐标系谱因式分解理论. 在螺旋坐标系下,Laplacian算子的表示矩阵具有Toeplitz结构,其快速求逆可由谱法LU分解实现. 基于二维离散光滑插值,文中还给出共轭梯度法与NMO相结合的沿时间切片逐层处理的离散光滑插值流程. 最后,应用该方法对模型数据和实际三维地震数据进行了处理.     

Spectral identification of materials by reflectance spectral library search

Spectral library search is emerging as a viable approach for material identification and mapping by reusing spectral knowledge gained from hyperspectral remote sensing across space and time. The potential of retrieving meaningful spectral material identifications in the presence of reflectance of spectra of various material types and with various similarity metrics has been assessed in this study. Test reflectance spectra of various vegetation, minerals, soils and urban material types are identified by searching through the composite reflectance spectral library obtained by combining various institutional reflectance spectral libraries. The accuracy of material identifications under various conditions: (i) in the presence of identical, similar and dissimilar spectra; (ii) in the presence of only identical and dissimilar spectra; and (iii) in the presence of only dissimilar spectra has been assessed with several similarity metrics. Results indicate the possibility of obtaining 100% accurate material identifications by library search if the spectral library contains identical spectra. However, the presence of a large number of similar spectra, despite the presence of identical spectra, is found to increase false positives, thereby reducing the accuracy of retrievals to 82% at best. Further, the accuracy of material identifications in the presence of similar spectra is similarity metric-dependent and varied from about 52% (obtained from Binary Encoding) to 82% (obtained from Normalized Spectral Similarity Score). Overall, results support the possibility of using independent reflectance spectral libraries for material identification while calling for robust spectral similarity metrics.     

Non-destructive estimation of potato leaf chlorophyll from canopy hyperspectral reflectance using the inverted PROSAIL model

Optimizing nitrogen (N) fertilization in crop production by in-season measurements of crop N status may improve fertilizer N use efficiency. Hyperspectral measurements may be used to assess crop N status by estimating leaf chlorophyll content. This study evaluated the ability of the PROSAIL canopy-level reflectance model to predict leaf chlorophyll content. Trials were conducted with two potato cultivars under different N fertility rates (0–300 kg N ha⁻¹). Canopy reflectance, leaf area index (LAI) and leaf chlorophyll and N contents were measured. The PROSAIL model was able to predict leaf chlorophyll content with reasonable accuracy later in the growing season. The low estimation accuracy earlier in the growing season could be due to model sensitivity to non-homogenous canopy architecture and soil background interference before full canopy closure. Canopy chlorophyll content (leaf chlorophyll content × LAI) was predicted less accurately than leaf chlrophyll content due to the low estimation accuracy of LAI for values higher than 4.5.     

Optical detection of Prorocentrum donghaiense blooms based on multispectral reflectance

Prorocentrum donghaiense is one of the most common red tide causative dinoflagellates in the Changjiang （Yangtze） River Estuary and the adjacent area of the East China Sea. It causes large-scale blooms in late spring and early summer that lead to widespread ecologic and economic damage. A means for distinguish- ing dinoflagellate blooms from diatom （Skeletonema costatum） blooms is desired. On the basis of measure- ments of remote sensing refectance [Rrs（λ）] and inherent optical parameters, the potential of using a mul- tispectral approach is assessed for discriminating the algal blooms due to P. donghaiense from those due to S. costatum. The behavior of two reflectance ratios [R1 = Rrs（560）/Rrs（532） and Re = Rrs（708）/Rrs（665）], suggests that differentiation of P. donghaiense blooms from diatom bloom types is possible from the current band setup of ocean color sensors. It is found that there are two reflectance ratio regimes that indicate a bloom is dominated by P. donghaiense; （1） R1 〉 1.55 and R2 〈 1.0 or （2） R1 〉 1.75 and R2 ≥ 1.0. Various sensitivity analyses are conducted to investigate the effects of the variation in varying levels of chlorophyll concentration and colored dissolved organic matter （CDOM） as well as changes in the backscattering ratio （bbp/bp） on the efficacy of this muitispectral approach. Results indicate that the intensity and inherent op- tical properties of the algal species explain much of the behavior of the two ratios. Although backscattering influences the amplitude of Rrs（λ）, especially in the 530 and 560 nm bands, the discrimination between P. donghaiense and diatoms is not significantly affected by the variation of bbp/bp. Since aCDOM （440） in coastal areas of the ECS is typically lower than 1.0 m-1 in most situations, the presence of CDOM does not interfere with this discrimination, even as SCDOM varies from 0.01 to 0.026 nm-1. Despite all of these effects, the dis- crimination of P. donghaiense blooms from diatom blooms based on multispectral measurements of Rrs（λ） is feasible.     

黄浦江水质指标与反射光谱特征的关系分析

利用地物光谱测量技术及同步配套的常规水质采样分析实验，对上海市黄浦江全河段水体进行调查研究。共选取31个典型站位进行了光谱测量和同步水质取样，每个测点分析了9个水质指标，包括TP，TN，CODMn，CODCr，NH3-N，BOD5，DO，悬浮物浓度和浊度。各水质参数之间存在较大的相关性，以总氮（TN）和悬浮物浓度两个水质指标为例，分析了这两个水质指标与单波段归一化反射率、反射率的一阶微分、不同波段之间反射率的比值以及反射率取对数等之间的相关关系，给出了以上两个水质指标的单波段归一化反射率、一阶微分反射率识别的特征波段，以及两个波段比值的最佳波段组合。建立了常规水质参数与水体光谱反射率之间的关系模型，为利用遥感技术监测水环境提供了基础。     

极端干旱区荒漠稀疏河岸林遥感分类研究

研究以位于极端干旱区的塔里木河干流中下游地区为例,基于Landsat TM影像,结合决策树分类、几何光学模型与光谱角匹配,解决混合像元信息分解,实现干旱区稀疏荒漠河岸林类别识别。首先从遥感视角的角度,将地物分解为目标和背景,提出塔里木河干流荒漠河岸林植被分类系统;其次以多变量决策树法将非荒漠植被信息剔除,采用几何光学模型模拟各类荒漠植被的像元光谱,最后以光谱角匹配的方法将荒漠植被进一步进行分解,得到塔里木河干流中下游地区典型研究区的植被分类专题图,分类精度结果表明:基于混合像元分解与几何光学模型的分类方法总精度达到了79.43%,Kappa系数为0.718,表明分类质量良好。     

基于叶片光谱的森林叶绿素浓度反演研究

准确估测森林叶片叶绿素浓度有助于进一步理解和模拟森林生态系统。基于江西省千烟州试验区主要树种光谱数据和相应叶绿素浓度等化学参量数据,使用高光谱指数和偏最小二乘(PLS)回归方法进行森林叶绿素浓度的反演,对具代表性的几种叶绿素相关的光谱指数进行反演能力评价,构造出表征叶绿素吸收特征的叶绿素吸收面积指数(CAAI),发现TCAR I(改进型的叶绿素吸收比率指数)和CAAI能较好地估测试验区森林主要树种的叶绿素浓度。针对叶片生化参量之间强相关现象,首先使用了叶片生化参量吸收特征分析方法,选定特定波段。然后对叶绿素浓度进行PLS回归估测,并从PLS获得估测模型的结果来评价和解释几种高光谱指数的叶绿素反演能力。     

渭河盆地中土层场地对地震动的放大作用

土层场地对地震动的影响较大,建(构)筑物的选址及其抗震设防必须考虑土层场地的放大作用,以避免或减轻其震害.汶川地震中,布设在渭河盆地中的数字强震动台网共有27个台站(包括2个基岩台站和25个土层台站)获得良好的主震加速度时程.利用这些加速度时程,选择汤峪台做为参考场地,基于考虑几何衰减的传统谱比法分析研究了25个土层场地的放大作用.结果表明:(1)汶川地震中渭河盆地土层场地上的加速度峰值(PGA)有随震中距增加而减小的趋势,且与土层的厚度无关.(2)根据汶川地震中渭河盆地土层场地的加速度反应谱特征及其对地震动的放大作用特征,渭河盆地的土层场地可分为3类:深厚、中厚和薄土层场地.而且,深厚、中厚和薄土层场地分别对地震动的低频、1 Hz附近和高频分量放大作用显著;值得注意的是,深厚土层对地震动的高频分量也有一定程度的放大,只不过相对于低频分量的放大作用较小.(3)汶川地震中渭河盆地的土层场地对不同方向地震动的放大作用不同,且有EW＞NS＞UD的关系.(4)汶川地震中从宝鸡到眉县沿着盆地边缘分布的Ⅶ度烈度异常区是盆地边缘效应和土层场地对地震动放大作用共同作用的结果.     

利用古温标恢复四川盆地主要构造运动时期的剥蚀量

地层剥蚀量是沉积盆地埋藏史和热史重建中一个关键的参数.本文利用石油钻井的系统古温标(Ro)资料,并结合多种恢复方法,得出四川盆地主要构造运动时期的剥蚀量.研究表明:加里东期,川东南剥蚀量较大,达2000m.东吴运动时期,川西南、川东南等距二叠纪玄武岩喷发区较近地区的剥蚀量较大,分别在260～450 m和800～900m:印支早期盆地整体遭受了抬升剥蚀作用,剥蚀厚度为100～500 m.印支中、晚期龙门山地区褶皱剥蚀,H1、Y1等钻井该时期的地层剥蚀量超过2000 m.燕山期周缘山系的继续隆升造成山前大范围地区的剥蚀:喜山期盆地周缘钻井的剥蚀量较大,在2000 m左右,而早期古隆起上的钻井如GJ、J13、Z12等钻井的地层剥蚀量则较小,在1000 m左右.可见,四川盆地不同地质时期及不同构造区位的剥蚀厚度都不尽相同,这一时空差异反映了构造运动对该区的差异影响.这一研究也表明,以系统的古温标资料(R_o)为基础,针对不同地质情况选用适当的反演方法并结合多种反演方法.能有效地恢复钻井在不同时期的剥蚀量.上述四川盆地各时期剥蚀地层厚度的恢复,对研究该区的构造、沉积和油气演化提供了基础数据.     

土壤有机质含量估测及其影响因素的光谱分析

 土壤颗粒大小差异使土壤反射光谱产生相应变化,影响土壤有机质含量等属性的光谱预测精度。本研究准备了颗粒粒径分别为2、0.25和0.15mm的土样,测定土壤有机质(Soil Organic Matter,SOM)含量,并于室内模拟条件下测定其反射光谱。通过分析不同粒径土样的原始(Raw)、多次散射校正(Multiple scattering correction, Msc)、一阶微分(First derivative, Fd)、连续统去除(Continuum removal, Cr)光谱与SOM含量之间的关系,筛选出与SOM含量相关性最强的Fd光谱单波段(2250nm, r=0.82, P<0.01),并建立线性回归模型;利用全波段光谱反射率,以偏最小二乘回归(Partial least square regression, PLSR)方法,确立2mm土样Msc处理光谱的PLSR模型为最优模型(RPD=3.56、R²=0.90、RMSEP=1.96g/kg)。土壤颗粒粒径对土壤光谱反射率变化有明显影响,但二者之间并非简单的线性关系,可能存在一个转折点;单变量(单波段光谱反射率)线性回归模型的预测能力,明显低于全波段反射光谱(Msc处理)-PLSR模型;土样样本容量对SOM含量预测精度有显著影响。因此,根据样本容量大小,选择合适的土壤颗粒粒径与光谱预处理方法组合可以提高预测精度。