Despite low spatial resolutions, thermal infrared bands (TIRs) are generally more suitable for mineral mapping due to fundamental tones and high penetration in vegetated areas compared to shortwave infrared (SWIR) bands. However, the weak overtone combinations of SWIR bands for minerals can be compensated by fusing SWIR-bearing data (Sentinel-2 and Landsat-8) with other multispectral data containing fundamental tones from TIR bands. In this paper, marble in a granitic complex in Mardan District (Khyber Pakhtunkhwa) in Pakistan is discriminated by fusing feature-oriented principal component selection (FPCS) obtained from the ASTER, Landsat-8 Operational Land Imager (OLI), Thermal Infrared Sensor (TIRS) and Sentinel-2 MSI data. Cloud computing from Google Earth Engine (GEE) was used to apply FPCS before and after the decorrelation stretching of Landsat-8, ASTER, and Sentinel-2 MSI data containing five (5) bands in the Landsat-8 OLI and TIRS and six (6) bands each in the ASTER and Sentinel-2 MSI datasets, resulting in 34 components (i.e., 2 × 17 components). A weighted linear combination of selected three components was used to map granite and marble. The samples collected during field visits and petrographic analysis confirmed the remote sensing results by revealing the region’s precise contact and extent of marble and granite rock types. The experimental results reflected the theoretical advantages of the proposed approach compared with the conventional stacking of band data for PCA-based fusion. The proposed methodology was also applied to delineate granite deposits in Karoonjhar Mountains, Nagarparker (Sindh province) and the Kotah Dome, Malakand (Khyber Pakhtunkhwa Province) in Pakistan. The paper presents a cost-effective methodology by the fusion of FPCS components for granite/marble mapping during mineral resource estimation. The importance of SWIR-bearing components in fusion represents minor minerals present in granite that could be used to model the engineering properties of the rock mass.