In this paper, we propose SHOPCA (Shape Operator-based Principal Component Analysis), a novel method for unsupervised metric learning and dimensionality reduction that incorporates differential geometric information into the covariance structure of classical PCA. SHOPCA regularizes the global covariance matrix using the mean shape operator, defined as the average of the absolute local shape operators estimated from the data manifold, steering principal components toward directions of both maximum variance and informative curvature. A single trace-normalized mixing coefficient α controls the regularization, recovering standard PCA at α=0 and a curvature-driven embedding as α→∞. We further introduce a fully unsupervised criterion for selecting α based on the spectral eigengap of the regularized covariance matrix, maximizing the relative separation between the top-d and remaining eigenvalues without using class labels. We evaluate SHOPCA on more than 50 real-world benchmark datasets, comparing it with PCA, ISOMAP, and UMAP using Adjusted Rand Index (ARI), Normalized Mutual Information (NMI), Fowlkes-Mallows index (FM), and V-measure. Results show that SHOPCA consistently improves clustering quality over PCA across a broad range of datasets and surpasses UMAP on small-sample settings, where iterative neighborhood-based manifold estimation can degrade. SHOPCA is computationally tractable, parameter-efficient, and applicable to domains requiring fully unsupervised, geometry-aware dimensionality reduction.
Figures & tables
Figure 1: The shape operator measures how a normal vector changes from the tail to the tip of a tangent vector v . In other words, it tells us how fast the normal changes, or roughly speaking, how fast a surface turns if you roll it along the floor in some direction.
#
Dataset
n
d
c
#
Dataset
n
d
c
1
iris
150
4
3
30
har
10299
561
6
2
wine
178
13
3
31
Olivetti_Faces
400
4096
40
3
digits
1797
64
10
32
balance-scale
625
4
3
4
diggle_table_a2
310
8
9
33
mammography
11183
6
2
5
satimage
6430
36
6
34
abalone
4177
8
3
6
mfeat-morphological
2000
6
10
35
Engine1
383
5
3
Table 1: Datasets used in the computational experiments: each row shows the dataset name, number of samples ( n ), number of features ( d ), and number of classes ( c ).
Regular PCA
Shape Operator PCA
Datasets
ARI
FM
VM
ARI
FM
VM
λ
iris
0.6254
0.7542
0.6722
0.7312
0.8222
0.7701
55.19
wine
0.6583
0.7734
0.6986
0.6764
0.7926
0.7252
12.46
digits
0.0358
0.2148
0.1322
0.3372
0.4228
0.5222
17.51
diggle_table_a2
0.2896
0.3805
0.5704
0.4708
0.5638
0.7078
127.60
satimage
0.0057
0.2840
0.0092
0.1012
0.4603
0.2627
9,95
Table 2: Comparison between standard PCA and the proposed SHOPCA on 30 real-world benchmark datasets. Both methods project the data onto a two-dimensional latent space, after which agglomerative clustering with Ward’s linkage is applied to recover the ground-truth number of classes. Clustering quality is assessed by the Adjusted Rand Index (ARI), Fowlkes–Mallows index (FM), and V-measure (VM). Bold entries indicate the best result per metric. SHOPCA achieves superior scores, with average relative improvements of 212% (ARI), 73% (FM), and 140% (VM) over regular PCA for these datasets.
Figure 2: Qualitative comparison of the two-dimensional embeddings produced by standard PCA and the proposed SHOPCA on representative datasets: cardiotocography (first row), mfeat-factors (second row) and texture (third row).
ISOMAP
Shape Operator PCA
Datasets
ARI
FM
VM
ARI
FM
VM
lambda
iris
0.7074
0.8084
0.7772
0.7312
0.8222
0.7701
55.19
digits
0.3288
0.4536
0.5855
0.3372
0.4228
0.5222
17.51
diggle_table_a2
0.4609
0.5450
0.7168
0.4708
0.5638
0.7078
127.60
mfeat-morphological
0.2873
0.3979
0.5048
0.4381
0.5242
0.6280
256.75
mfeat-fourier
0.2771
0.3610
0.4364
0.2720
0.3643
0.4760
18.34
Table 3: Comparison between ISOMAP and the proposed SHOPCA on 25 real-world benchmark datasets. Both methods project the data onto a two-dimensional latent space, after which agglomerative clustering with Ward’s linkage is applied to recover the ground-truth number of classes. Clustering quality is assessed by the Adjusted Rand Index (ARI), Fowlkes–Mallows index (FM), and V-measure (VM). Bold entries indicate the best result per metric.
Figure 3: Qualitative comparison of the two-dimensional embeddings produced by standard PCA and the proposed SHOPCA on representative datasets: cardiotocography (first row), UMIST_Faces_Cropped (second row) and texture (third row).
UMAP
Shape Operator PCA
Datasets
ARI
FM
VM
ARI
FM
VM
lambda
iris
0.6312
0.7606
0.7039
0.7312
0.8222
0.7701
55.19
diggle_table_a2
0.4073
0.4802
0.6637
0.4708
0.5638
0.7078
127.60
balance_scale
0.0942
0.4419
0.0811
0.2415
0.5376
0.2522
9.50
zoo
0.5203
0.6225
0.7311
0.8516
0.8895
0.8224
697.51
glass
0.1496
0.3383
0.2652
0.1589
0.4203
0.3024
21.38
Table 4: Comparison between UMAP and the proposed SHOPCA on 28 real-world benchmark datasets. Both methods project the data onto a two-dimensional latent space, after which agglomerative clustering with Ward’s linkage is applied to recover the ground-truth number of classes. Clustering quality is assessed by the Adjusted Rand Index (ARI), Fowlkes–Mallows index (FM), and V-measure (VM). Bold entries indicate the best result per metric.
Figure 4: Qualitative comparison of the two-dimensional embeddings produced by standard PCA and the proposed SHOPCA on representative datasets: UMIST_Faces_Cropped (first row) and mammography (second row).