Accurate lunar rock detection is vital for lunar exploration. However, the existing methods are sensitive to factors, such as the uneven lighting and terrain relief. To address these issues, a novel method combining multiscale phase feature type maps (PFTMs) and phase congruency moment maps (PCMMs) is proposed. First, rock seeds are detected through phase congruency and gradient analysis. Second, a strategy called the local scale saliency score (LSSS) is proposed to adaptively estimate the optimal scale layer for candidate rock detection. Within this layer, the specifically designed local-contextual and global-contextual (LGC) features are employed to identify the regions of interests (ROIs) for rocks. Subsequently, the process of filtering false positives (FPs) involves the utilization of geometric metrics and scale feature analysis. Finally, a specially designed edge detector named the bilateral local maximum PCMM-PFTM path is proposed to describe the edges of the rocks. Tests on Chang’E-3 and Chang’E-5 Landing Camera (LCAM) images show the proposed method’s robustness in detecting lunar rocks of varying sizes and reflectance, achieving F1-scores ranging from 0.919 to 0.947.

The mineral composition of lunar regolith influences brightness temperature (TB) as observed by the microwave radiometer (MRM); however, the large-scale spatial relationship between TB and mineral abundance has yet to be sufficiently revealed. This study aims to quantify the impact of specific mineral abundances (plagioclase and ilmenite) on TB distribution, using MRM 37-GHz data from Chang’E-2 and mineral abundance products from Kaguya. We applied hour angle correction and latitude normalization to produce high-accuracy TB maps and developed a self-adaptive moving-window method to remove strip noise to produce higher precision mineral abundance maps. Using these two types of maps, we performed comprehensive large-scale spatial analysis of TB and mineral abundance using a bivariate spatial autocorrelation model and a geographically weighted regression (GWR) approach considering spatial similarity and heterogeneity, respectively. The bivariate analysis indicates a negative spatial correlation between TB and plagioclase abundance, while a positive spatial correlation between TB and ilmenite abundance. In addition, bivariate anomalies, including both hot and cold spots of diurnal TB amplitude (i.e., noon minus nighttime), were identified through the simultaneous consideration of TB and mineral abundance. The GWR analysis reveals regional variations in the impact of mineral abundances on diurnal TB amplitudes. These correlations can be attributed to the spatial distributions and variations in TB, which arise from the unique dielectric and thermal properties of minerals across distinct regions on the Moon. These findings contribute to a better understanding of subsurface thermal behavior and regimes, enhancing our comprehension of lunar evolution.

The positioning accuracy of synthetic aperture radar (SAR) images is affected by factors, such as satellite platform instability, aging of on-board instruments, and environmental changes. Geometric calibration is a commonly employed and cost-effective method to enhance the positioning accuracy of SAR images. The classical point-based geometric calibration (PB-GC) model, however, only utilizes the location of ground control points (GCPs) and does not fully exploit the spatial relationships among the GCPs. This study introduces a high-precision geometric calibration method that builds upon the classical model for calibrating SAR imaging systems. This method incorporates the Co-Line-GC and Co-Circle-GC models, where the former uses GCPs distributed on a line while the latter uses GCPs distributed on a circle. The results reveal that, compared to the classical model, our approach enhances the positioning accuracy of Gaofen-3 and Sentinel-1A SAR images by approximately 2 m in eastern China, achieving a mean positioning accuracy of 3.02 m. In terms of calibration performance, a comparison between postcalibrated and precalibrated images indicates that the images are shifted, not distorted, and a better match of the same features between different scenes in the image mosaic is observed after calibration. The improved positioning accuracy of SAR images significantly contributes to global remote sensing mapping, land use change monitoring, and ground target detection applications.

The investigation of small lunar craters holds scientific and engineering significance. This paper presents a novel method for detecting small lunar craters. It consists of three stages: seed detection, candidate crater detection, and crater evaluation. Firstly, crater seeds are identified through morphological operations as pixels with the highest local gradient and specific gradient direction. Secondly, the optimal scale for each seed is estimated based on the maximum response of the established phase congruency maximum moment (PCMM) scale space. For detecting very small craters, a crater detector called statistical morphological constraint path-sets (SMPS), which leverages image spatial domain features, is proposed. It configures the image as a weighted directed graph, using path-sets centered on seeds to flexibly detect highlights and shadow regions of craters. For detecting craters with larger optimal scale, another crater detector named structural consistency constrained multi-paths (SCMP) is proposed, utilizing the frequency phase features. The core idea of SCMP is to configure the phase feature type (PFT) map with the optimal scale as a directed graph. Centered on the seed, the multi-path operator is designed to detect craters. Unsupervised discriminative correlation filters (UDCFs) are trained with HOG features from images or PFT maps to validate candidate craters. The results indicate that for images with a resolution of 0.5-2m/pixel, the proposed method demonstrates good detection performance for small craters with diameters of less than 5 m, 5-10 m, and greater than 10 m, with an average detection rate of 0.89, 0.91, and 0.92, respectively.

As the accuracy of satellite laser altimetry is susceptible to real-time atmospheric conditions, with-in footprint topography fluctuation, and detector noise, etc., we proposed a method by comprehensively analyzing the laser ranging error and evaluation labels to extract high-accuracy elevation control points that is suitable for satellite imagery based topographic mapping applications. Using the ICESat laser altimetry data, a global high accuracy laser altimetry dataset including more than 60 million control points, based on the laser altimetry ranging error model and waveform quality analysis is proposed by the paper. For land areas, except for areas of water, snow/ice, and polar ice sheets, the dataset can provide the elevation control points for worldwide satellite topographic mapping using high spatial resolution imageries or other science researches that depend on accurate earth’s elevation information. We further used airborne lidar data from six study areas around the world to carefully validate the dataset’s accuracy. The results showed that, this dataset can meet the accuracy requirement of global mapping using high spatial resolution satellite imageries in terrains with a slope below 25°. Compared to the raw dataset, the proportion of footprint elevations that conform to the accuracy standard (0.5m@ slope<2°, 1.5m@ 2°≤slope<6° and 3m@ 6°≤slope<25°) is increased from 68.24%, 59.97% and 26.52% to 87.58%, 90.04% and 83.91% respectively. This method can assure that its extracted results’ accuracy is either very close to or better than that obtained by the methods proposed in relevant studies, with a much larger number of laser footprints have been reserved.

With more detailed spatial information being represented in very-high-resolution (VHR) remote sensing images, stringent requirements are imposed on accurate image classification. Due to the diverse land-objects with intraclass variation and interclass similarity, efficient and fine classification of VHR images especially in complex scenes is challenging. Even for some popular deep learning (DL) frameworks, geometric details of land-object may be lost in deep feature levels, so it is difficult to maintain the highly-detailed spatial information (e.g., edges, small objects) only relying on the last high-level layer. Moreover, many of the newly developed DL methods require massive well-labeled samples, which inevitably deteriorates the model generalization ability under the few-shot learning. Therefore, in this paper, a lightweight shallow-to-deep feature fusion network (SDF2N) is proposed for VHR image classification, where the traditional machine learning (ML) and DL schemes are integrated to learn rich and representative information to improve the classification accuracy. In particular, the shallow spectral-spatial features are first extracted, and then a novel triple-stage fusion (TSF) module is designed to learn the saliency and discriminative information at different levels for classification. The TSF module includes three feature fusion stages, i.e., low-level spectral-spatial feature fusion, middle-level multi-scale feature fusion, and high-level multi-layer feature fusion. The proposed SDF2N takes advantages of the shallow-to-deep features, which can extract representative and complementary information of crossing layers. It is important to note that even with limited training samples, the SDF2N still can achieve satisfying classification performance. Experimental results obtained on three real VHR remote sensing data sets including two multispectral and one airborne hyperspectral images covering complex urban scenarios confirm the effectiveness of the proposed approach compared with the state-of-the-art methods.

Crater detection from planetary images is a challenging issue due to the complicated variations in geometry shape, illumination, and scale. An automatic crater detection algorithm (CDA) that is robust to these factors is, therefore, necessary. In this article, a novel automatic CDA that is robust to these factors is proposed to detect the multiscale craters of the Moon. The proposed method consists of two main steps: 1) in the hypothesis generation (HG) step, a novel feature operator called the path-profile, which is constructed based on the self-defined adjacency graph and a path descriptor, is presented to derive the highlight-shadow feature of craters for detecting candidate craters. 2) In the hypothesis verification (HV) step, based on the idea of anomaly detection, the isolation forest algorithm which is an unsupervised learning anomaly detection method is applied to eliminate falsely detected craters. Lunar Reconnaissance Orbiter Camera Wide Angle Camera and Narrow Angle Camera images and Chang’E-4 landing camera images were used to test the accuracy and robustness of the proposed method. The experimental results indicate that: on average, the accuracy of the detection result of the HG step is about 90%, and the HV step can further improve this by 3%–4%. The proposed method is a reliable way to detect multiscale lunar craters for various resolutions images with diameters ranging from five pixels to hundreds of pixels, and it is robust to the different terrains and illumination conditions on the Moon.