Geometric-constrained multi-view image matching method based on semi-global optimization

Targeting at a reliable image matching of multiple remote sensing images for the generation of digital surface models, this paper presents a geometric-constrained multi-view image matching method, based on an energy minimization framework. By employing a geometrical constraint, the cost value of the energy function was calculated from multiple images, and the cost value was aggregated in an image space using a semi-global optimization approach. A homography transform parameter calculation method is proposed for fast calculation of projection pixel on each image when calculating cost values. It is based on the known interior orientation parameters, exterior orientation parameters, and a given elevation value. For an efficient and reliable processing of multiple remote sensing images, the proposed matching method was performed via a coarse-to-fine strategy through image pyramid. Three sets of airborne remote sensing images were used to evaluate the performance of the proposed method. Results reveal that the multi-view image matching can improve matching reliability. Moreover, the experimental results show that the proposed method performs better than traditional methods.     

Novel fusion approach on automatic object extraction from spatial data: case study Worldview-2 and TOPO5000

The automatic extraction of information content from remotely sensed data is always challenging. We suggest a novel fusion approach to improve the extraction of this information from mono-satellite images. A Worldview-2 (WV-2) pan-sharpened image and a 1/5000-scaled topographic vector map (TOPO5000) were used as the sample data. Firstly, the buildings and roads were manually extracted from WV-2 to point out the maximum extractable information content. Subsequently, object-based automatic extractions were performed. After achieving two-dimensional results, a normalized digital surface model (nDSM) was generated from the underlying digital aerial photos of TOPO5000, and the automatic extraction was repeated by fusion with the nDSM to include individual object heights as an additional band for classification. The contribution was tested by precision, completeness and overall quality. Novel fusion technique increased the success of automatic extraction by 7% for the number of buildings and by 23% for the length of roads.     

Dynamic land use change simulation using cellular automata with spatially nonstationary transition rules

The dynamic relationships between land use change and its driving forces vary spatially and can be identified by geographically weighted regression (GWR). We present a novel cellular automata (GWR-CA) model that incorporates GWR-derived spatially varying relationships to simulate land use change. Our GWR-CA model is characterized by spatially nonstationary transition rules that fully address local interactions in land use change. More importantly, each driving factor in our GWR model contains effects that both promote and resist land use change. We applied GWR-CA to simulate rapid land use change in Suzhou City on the Yangtze River Delta from 2000 to 2015. The GWR coefficients were visualized to highlight their spatial patterns and local variation, which are closely associated with their effects on land use change. The transition rules indicate low land conversion potential in the city’s center and outer suburbs, but higher land conversion potential in the inner near suburbs along the belt expressway. Residual statistics show that GWR fits the input data better than logistic regression (LR). Compared with an LR-based CA model, GWR-CA improves overall accuracy by 4.1% and captures 5.5% more urban growth, suggesting that GWR-CA may be superior in modeling land use change. Our results demonstrate that the GWR-CA model is effective in capturing spatially varying land transition rules to produce more realistic results, and is suitable for simulating land use change and urban expansion in rapidly urbanizing regions.     

孔隙性砂岩中单条变形带与簇状变形带的特征与区别——以远安地堑上白垩统红花套组中变形带为例

孔隙性砂岩中的变形带研究近年来取得长足进展,但有关单条变形带与簇状变形带如何区分和界定仍是一个较模糊的问题。本文以远安地堑上白垩统红花套组高孔砂岩中发育的剪切增强压实带和压实剪切带为例,对单条及簇状变形带的特征和区别做出系统分析。研究表明,虽然单条变形带的长度、宽度、错断位移、凸起高度等均明显小于簇状变形带,但是不同运动学成因的变形带宏观特征、微观结构不同,因此不能单纯以规模大小区分二者,需要结合带密度将簇状变形带与变形带列、网格构造等相区分,还需要结合其微观结构特征进行准确鉴定。单条及簇状剪切增强压实带特征类似,均不具核心碎裂带,颗粒主要为中等破碎,多表现为多个不连续分布的碎裂核。簇状剪切增强压实带内部的非均质性更强。单条压实剪切带具有核心碎裂带,依剪切力大小不同,其核心碎裂带内颗粒粉碎程度及边缘过渡带宽度不同,常伴随微小剪切缝痕迹。簇状压实剪切带在镜下可清晰识别出不同的核心碎裂带、带间区域及微小剪切缝,可见S—C组构和断续分布的滑动面。单条及簇状变形带的发育会不同程度地改变储层的物性并可能会对油气运移和开采产生影响。     

1948年川西理塘 M7. 3地震地表破裂特征及Riedel剪切构造分析

1948年川西理塘M7. 3地震是川滇菱形块体内部近一个世纪以来发生的震级最大的走滑型地震。在对此次同震地表破裂进行详细调查基础上,利用差分GPS对同震地表破裂带进行了精确测量与统计分析。结果揭示该地表破裂的现存长度为36 km,北端始于无量河以北,往东南沿藏坝盆地北东缘、德巫盆地东南缘,延伸至德巫乡北,分为南、北两段,而在交德附近存在约3 km长的地表破裂空区。对同震地表破裂的线密度和同震水平位错量进行分段统计后揭示,此次地震的宏观震中应位于德巫盆地中部交德东南约4~5 km处。对理塘同震地表破裂的Riedel剪切分析结果指示,该破裂带主要由R剪切组成,以发育雁列状排列的挤压鼓包(Push- up)为主要特征,伴有少量R′剪切与T裂缝,缺少P型与X型剪切。其中R剪切占95%以上,其在藏坝段(北段)的优势方向为318°,德巫段(南段)为315°,整条地表破裂带的R剪切平均方向为316°。同时发现,此次地震形成的雁列状挤压鼓包主要以平缓的“弧形”为主,这与1981年道孚MS 6. 9地震和2010年玉树MS 7. 1地震地表破裂带中出现大量反“S”形挤压鼓包有所不同,推断与走滑断裂滑动速率密切相关。沿强滑动速率走滑断层地震破裂带的Riedel剪切发育会更为完善,挤压鼓包也更发育,易形成反“S”形,反之则以平缓的“弧形”为特征。综合分析认为,1948年理塘同震地表破裂带的展布主要受藏坝盆地与德巫盆地控制,而且藏坝段(北段)与德巫段(南段)的R剪切方向存在偏差,这可能与两个拉分盆地的发育程度有关。     

复垦土地样品中石油类物质加速溶剂萃取-荧光分光光度法分析方法研究

作为土地环境质量的一项重要指标,石油类物质含量决定了土地的用途。对于复垦土地样品,在一定程度上反映了污染土地的复垦、治理效果。土壤样中的石油类物质测定方法主要有气相色谱法、红外分光光度法、紫外光谱法等。气相色谱法主要应用于测定样品中的饱和烷烃;红外分光光度法通过测定不同波数下的特征吸收值来表征样品中石油类物质,但采用四氯乙烯为溶剂,试剂环境危害大,前处理效率低,检出限高。本文针对复垦土地样品中石油类物质含量不高,组分多为难挥发且结构复杂、基质干扰大的特点,采用正己烷为溶剂、加速溶剂萃取(ASE)前处理,荧光分光光度计测定样品中的石油类物质。以海洋环境监测石油标准物质为标准样品配制工作曲线,线性相关系数r~2=0.9997;方法检出限为0.40mg/kg,检测下限为1.60mg/kg;实际样品方法精密度为1.10%～8.76%;基质加标回收率为89.0%～95.7%;检测实际样品的结果与现行有效的HJ 1051—2019红外分光光度法基本一致,高含量样品的测定结果高于红外法。本方法所用溶剂正己烷的毒性小于四氯乙烯,前处理方法ASE的自动化程度高,萃取效率高,精密度提高11.5%～67.3%。对于石油类组分结构相对复杂、难提取的样品,本方法检出限低于红外法(4mg/kg)。     

Deformation patterns of upper Quaternary strata and their relation to active tectonics,Po Basin,Italy

Despite increased application of subsurface datasets below the limits of seismic resolution, reconstructing near‐surface deformation of shallow key stratigraphic markers beneath modern alluvial and coastal plains through sediment core analysis has received little attention. Highly resolved stratigraphy of Upper Pleistocene to Holocene (Marine Isotope Stage 5e to Marine Isotope Stage 1) alluvial, deltaic and coastal depositional systems across the southern Po Plain, down to 150 m depth, provides an unambiguous documentation on the deformation of previously flat‐lying strata that goes back in time beyond the limits of morphological, historical and palaeoseismic records. Five prominent key horizons, accurately selected on the basis of their sedimentological characteristics and typified for their fossil content, were used as highly effective stratigraphic markers (M1 to M5) that can be tracked for tens of kilometres across the basin. A facies‐controlled approach tied to a robust chronology (102 radiocarbon dates) reveals considerable deformation of laterally extensive nearshore (M1), continental (M2 and M3) and lagoon (M4 and M5) marker beds originally deposited in a horizontal position (M1, M4 and M5). The areas where antiformal geometries are best observed are remarkably coincident with the axes of buried ramp anticlines, across which new seismic images reveal substantially warped stratal geometries of Lower Pleistocene strata. The striking spatial coincidence of fold crests with the epicentres of historic and instrumental seismicity suggests that deformation of marker beds M1 to M5 might reflect, in part at least, syntectonically generated relief and, thus, active tectonism. Precise identification and lateral tracing of chronologically constrained stratigraphic markers in the ¹⁴C time window through combined sedimentological and palaeoecological data may delineate late Quaternary subsurface stratigraphic architecture at an unprecedented level of detail, outlining cryptic stratal geometries at the sub‐seismic scale. This approach is highly reproducible in tectonically active Quaternary depositional systems and can help to assess patterns of active deformation in the subsurface of modern alluvial and coastal plains worldwide.     

云量对祁连山老虎沟12号冰川表面能量平衡的影响

为探讨云量对冰川表面能量平衡（SEB）的影响,利用架设在老虎沟12号冰川（简称12号冰川）消融区（4 550 m a.s.l.）的自动气象站资料,结合能量平衡模型计算各能量分量并分析其季节变化,通过云量参数化方案获取云量因子并量化其对冰川表面能量收支的影响。结果表明：净短波辐射为冰川表面主要的能量来源（92%）,净长波辐射为主要能量支出（61%）,二者均受云量影响,但云的短波辐射效应更强（-37 W?m^-2）。云量通过影响辐射收支和湍流通量进而影响冰川表面能量收支,随云量的增加,冰川表面获得的能量减少,冰川消融速率降低。与其他区域的冰川表面能量收支对比,除地理位置、反照率、气温等因素外,海拔和云量的影响也非常显著。     

内蒙古鄂尔多斯高原古冰缘遗迹科学考察研究进展

冰缘遗迹（特别是冷生楔形构造及融冻褶皱）是重建古气候及第四纪晚期多年冻土环境的重要证据。内蒙古鄂尔多斯高原是我国北方地区冰缘现象最为发育的地区之一。为准确了解鄂尔多斯高原冰缘遗迹类型及其分布特征、区域冻土演化历史等,中国科学院西北生态环境资源研究院和荷兰自由大学共同组成科研小组,于2018年5—6月组织了“鄂尔多斯高原冰缘遗迹科学考察”。考察区域涉及靖边—城川—乌审旗—鄂尔多斯东胜区一带约12 000 km²的范围。考察内容主要包括鄂尔多斯高原冰缘遗迹类型及特征、分布区域、各类型冰缘遗迹所指示的气候条件的初步推断等。结果表明：冻融褶皱和冷生楔体构造是鄂尔多斯高原主要存在的两大类冰缘遗迹。基于本次考察中关于冰缘遗迹的分布与特征等新发现,并综合前人研究成果,初步推断：在气温极低、多年冻土非常发育的时段,有利于形成各类冷生楔状构造,如冰楔假形和大型原生砂楔等;在气候转暖、多年冻土退化,但还没有全部融化完阶段,可能形成融冻褶皱;区域性大面积分布和成群出现的融冻褶皱一般反映较暖气候环境下,多年冻土层上部已退化到一定程度。基于光释光（OSL）年代测试结果,结合冰缘遗迹的特征及其所指示的古气候环境,初步重建了鄂尔多斯5万年以来的冻土环境变化序列。区内多年冻土在多年冻土最大期（LPM,25~19 ka BP）时最发育,以大面积连续多年冻土为主;之后,随气温转暖,总趋势呈退化状态,多年冻土分布逐渐变为片状→岛状→零星斑状,直至现今全部融完变为深季节冻土区。     

四川稻城海子山方枝柏树轮宽度与气候因子的响应关系研究

对四川稻城县海子山高原上的方枝柏进行树芯采集,经交叉定年后建立林线上限位置的方枝柏树轮宽度年表。标准化年表与气象资料的响应分析结果揭示了该地区树木生长主要受当年生长季前期温度条件的控制。轮宽指数与前一年9月到当年2月的平均气温呈显著正相关（R=0.56,P<0.01）。由此重建了该地区1850—2019年这170 a的前一年9月到当年2月平均气温的变化,转换函数方差解释量31.7%。重建结果与邻近地区树轮记录的温度变化具有较好的一致性。重建序列显示工业革命以来,研究区经历较冷的时期有：1870—1890年和1960—1980年;较暖的时期有：1884—1892年、1919—1925年、1940—1960年以及1980—2019年。1940—1960年是20世纪最暖的时期,20世纪60年代后温度呈现下降趋势,1980年开始气温保持持续升高,2010年之后气温有小幅下降。重建气温序列与太阳黑子数变化和北大西洋多年代际涛动AMO有较强的响应,揭示了太阳活动、海温等因子可能对该区温度变化产生重要影响。