In arid and semi-arid regions, it is essential to monitor the spatiotemporal variability and dynamics of vegetation. Among other provinces of Pakistan, Punjab has produced a significant number of crops. Recently, Punjab, Pakistan, has been described as a global hotspot for extremes of climate change. In this study, the soil adjusted vegetation index (SAVI), normalized vegetation difference index (NDVI), and enhanced vegetation index (EVI) were comprehensively evaluated to monitor vegetation change in Punjab, Pakistan. The time-series MODIS (Moderate Resolution Imaging Spectroradiometer) data of different periods were used. The mean annual variability of the above vegetation indices (VIs) from 2000 to 2019 was evaluated and analyzed. For each type of vegetation, two phenological metrics (i.e., for the start of the season and end of the season) were calculated and compared. The spatio-temporal image analysis of the mean annual vegetation indices revealed similar patterns and varying vegetation conditions. In the forests and vegetation areas with sparse vegetation, the EVI showed high uncertainty. The phenological metrics of all vegetation indices were consistent for most types of vegetation. However, the NDVI result had the greatest variance between the start and end of season. The lowest annual VI variability was mainly observed in the southern part of the study area (less than 10% of the study area) based on the statistical analysis of spatial variability. The mean annual spatial variability of NDVI was <20%, SAVI was 30%, and EVI ranged between 10–20%. More than 40% of the variability was observed in the NDVI and SAVI vegetation indices.

The Remote Sensing Ecological Index (RSEI) model is widely used for large-scale, rapid Ecological Environment Quality (EEQ) assessment. However, both the RSEI and its improved models have limitations in explaining the EEQ with only two-dimensional (2D) factors, resulting in inaccurate evaluation results. Incorporating more comprehensive, three-dimensional (3D) ecological information poses challenges for maintaining stability in large-scale monitoring, using traditional weighting methods like the Principal Component Analysis (PCA). This study introduces an Improved Remote Sensing Ecological Index (IRSEI) model that integrates 2D (normalized difference vegetation factor, normalized difference built-up and soil factor, heat factor, wetness, difference factor for air quality) and 3D (comprehensive vegetation factor) ecological factors for enhanced EEQ monitoring. The model employs a combined subjective–objective weighting approach, utilizing principal components and hierarchical analysis under minimum entropy theory. A comparative analysis of IRSEI and RSEI in Miyun, a representative study area, reveals a strong correlation and consistent monitoring trends. By incorporating air quality and 3D ecological factors, IRSEI provides a more accurate and detailed EEQ assessment, better aligning with ground truth observations from Google Earth satellite imagery.

3D building models are an essential data infrastructure for various applications in a smart city system, since they facilitate spatial queries, spatial analysis, and interactive visualization. Due to the highly complex nature of building structures, automatically reconstructing 3D buildings from point clouds remains a challenging task. In this paper, a Roof Attribute Graph (RAG) method is proposed to describe the decomposition and topological relations within a complicated roof structure. Furthermore, top-down decomposition and bottom-up refinement processes are proposed to reconstruct roof parts according to the Gestalt laws, generating a complete structural model with a hierarchical topological tree. Two LiDAR datasets from Guangdong (China) and Vaihingen (Germany) with different point densities were used in our study. Experimental results, including the assessment on Vaihingen standardized by the International Society for Photogrammetry and Remote Sensing (ISPRS), show that the proposed method can be used to model 3D building roofs with high quality results as demonstrated by the completeness and correctness metrics presented in this paper.

Three-dimensional (3D) building models are closely related to human activities in urban environments. Due to the variations in building styles and complexity in roof structures, automatically reconstructing 3D buildings with semantics and topology information still faces big challenges. In this paper, we present an automated modeling approach that can semantically decompose and reconstruct the complex building light detection and ranging (LiDAR) point clouds into simple parametric structures, and each generated structure is an unambiguous roof semantic unit without overlapping planar primitive. The proposed method starts by extracting roof planes using a multi-label energy minimization solution, followed by constructing a roof connection graph associated with proximity, similarity, and consistency attributes. Furthermore, a progressive decomposition and reconstruction algorithm is introduced to generate explicit semantic subparts and hierarchical representation of an isolated building. The proposed approach is performed on two various datasets and compared with the state-of-the-art reconstruction techniques. The experimental modeling results, including the assessment using the International Society for Photogrammetry and Remote Sensing (ISPRS) benchmark LiDAR datasets, demonstrate that the proposed modeling method can efficiently decompose complex building models into interpretable semantic structures.

Accurate crop mapping is vital for monitoring agricultural resources, food security, and ecosystem sustainability. Advances in high-resolution sensing technologies now enable precise, large-scale crop mapping, improving agricultural management and decision-making. However, in scenarios where balancing precision and computational resources is important, obtaining the optimal feature combination (especially newly proposed features) and strategies from the rich feature sets contained in multi-source remote sensing imagery remains one of the challenges. In this paper, we propose a hierarchical feature optimization method, incorporating a newly reported vegetation feature, for mapping crop types by combining the Sentinel-1 Synthetic Aperture Radar (SAR) and Sentinel-2 optical imagery within the Google Earth Engine (GEE) platform. The method first calculates spectral features, texture features, polarization features, vegetation index features, and crop phenological features, with a particular focus on infrared band features and the newly developed Kernel Normalized Difference Vegetation Index (kNDVI). These 126 features are then selected to construct 15 crop type mapping models based on different feature combinations and a random forest (RF) classifier. Feature selection was performed using the feature correlation analysis and random forest recursive feature elimination (RF-RFE) to identify the optimal subset. The experiment was conducted in the Linhe region, covering an area of 2333 km². The resulting 10 m crop map, generated by the optimal model (Model 15) with 34 key features, demonstrated that integrating multi-source features significantly enhances mapping accuracy. The model achieved an overall accuracy of 90.10% across five crop types (corn, wheat, sunflower, soybean, and beet), outperforming other representative feature optimization methods, Relief-F (87.50%) and CFS (89.60%). The study underscores the importance of feature optimization and reduction of redundant features while also showcasing the effectiveness of red edge and infrared features, as well as the kNDVI, in mapping crop type.

Building height is a critical parameter for urban analysis, yet accurately estimating it from ICESat-2 photon-counting LiDAR data remains challenging due to pervasive noise photons and uneven noise distribution. To address the limitations of fixed-threshold denoising methods and improve adaptability across varying density conditions, this study proposes a dual-stage denoising framework for ICESat-2 ATL03 photon data. In the first stage, local photon densities are estimated within a reliable radius, log-transformed, and stratified into multiple levels. Adaptive thresholds are then applied at each level to suppress low-density noise while minimizing over-filtering in sparse regions. In the second stage, residual feedback-driven adaptive fitting strategy is applied along the ground track, where polynomial fitting was performed in sliding windows, with the window size dynamically adjusted based on residuals to refine local structures and eliminate outliers. The experiment was conducted in South Holland and Friesland, across 84 ICESat-2 tracks, where quantitative evaluations under varying day/night and beam conditions confirmed the effectiveness of the proposed framework. For denoising, the proposed method achieved high denoising accuracy, with F1-scores exceeding 0.97 in most cases, outperforming previous methods. Furthermore, building heights derived from footprint buffering and elevation differencing are validated against airborne LiDAR, yielding coefficient of determination (R²) values of 0.7235 and 0.9487 for the two regions, with root mean square error (RMSE) values of 1.5045 m and 1.8849 m, respectively. This study confirms the effectiveness and robustness of the proposed dual-stage framework, demonstrating its strong capability for both noise suppression in ICESat-2 ATL03 photon data and the subsequent accurate estimation of building heights.