We introduce radial basis function (RBF) kernels for Gaussian processes that use structured hyperparameter sharing for image-valued inputs. For high-dimensional inputs such as images, a fully anisotropic RBF kernel assigns a separate lengthscale to each input dimension, which can introduce many hyperparameters and may lead to overfitting and poorly calibrated predictive uncertainty when data are limited. We propose structured relevance determination by tying lengthscales across groups of adjacent pixels and within each colour channel. This yields patch-wise and channel-wise constrained anisotropic RBF parametrisations that reduce the number of independently learned lengthscales. We further introduce channel-additive Gaussian processes for multi-channel images and combine them with structured relevance determination, resulting in a parameter-efficient decomposition of the latent function into components defined on channel subsets, enabling some non-local generalisation. The proposed kernels are evaluated on MNIST and a compact binary aircraft-detection dataset of low-resolution RGB satellite image chips. We also show how channel-wise structured relevance determination can be incorporated into convolutional Gaussian processes. We observe improved predictive performance compared to isotropic and fully anisotropic RBF baselines.
Figures & tables
Figure 1: Results for MNIST with 6,000 training points (top row) and 60,000 training points (bottom row) with isotropic RBF , RBF - ARD and RBF - APRD each with patches of shape 2×2 , 4×4 , and 7×7 . Left: evidence lower bound, middle: classification error, right: negative log predictive probability.
MNIST 6,000
MNIST 60,000
model
# lengthscales
error (%)
NLPP ( ×102 )
error (%)
NLPP ( ×102 )
isotropic
1
5.35
32.58
2.13
1.31
ARD
784
7.08
45.76
2.24
1.26
APRD (2 × 2)
196
5.88
31.49
1.94
1.15
APRD (4 × 4)
49
4.65
27.46
1.84
1.12
APRD (7 × 7)
16
4.61
27.73
1.95
1.16
Table 1: Test results for MNIST. The number of lengthscales is reported per latent GP .
Figure 2: Sample from the Planes in Satellite Imagery dataset. Top row: class 1, bottom row: class 0.
Figure 3: Results for experiment 1 on Planes in Satellite Imagery dataset. Left: evidence lower bound, middle: classification error, right: negative log predictive probability.
model
# kernel hyperparameters
error (%)
NLPP ( ×102 )
isotropic
2
3.31
9.07
ARD
1201
2.8
8.08
ACRD
4
3.16
8.67
APRD (4 × 4) + ACRD
29
2.86
7.94
channel-additive
19
2.89
8.19
channel-additive APRD (4 × 4)
44
2.55
7.32
Table 2: Test results for experiment 1 on Planes in Satellite Imagery dataset.
Figure 4: Results for experiment 2 on Planes in Satellite Imagery dataset. Left: evidence lower bound, middle: classification error, right: negative log predictive probability.
model
error (%)
NLPP ( ×102 )
Convolutional
2.85
7.82
Convolutional- ACRD
2.61
7.19
Convolutional + RBF - ARD
2.55
6.69
Convolutional- ACRD + channel-additive APRD (4x4)
2.3
5.69
Table 3: Test results for experiment 2 on Planes in Satellite Imagery dataset.
We present a theoretically grounded Gaussian process framework that leverages neural feature maps to construct expressive kernels. We show that the learned feature map can be interpreted as an optimal low-rank approximation to a Gram matrix derived from an implied RKHS, from which we establish consistency of the GP posterior. We further analyse the spectral properties of the induced kernels and introduce product feature-map kernels to address oversmoothing. This simple yet powerful approach enables fast, scalable, and accurate exact GP inference with minimal upfront work. The flexibility of kernel design supports seamless application to both regression and classification tasks across diverse data modalities, including tabular inputs and structured domains such as images. On benchmark datasets, this approach surpasses pre-existing methods in terms of accuracy and training and prediction efficiency.
Object detection in remote sensing images (RSIs) is challenged by the coexistence of geometric and spatial complexity: targets may appear with diverse aspect ratios, while spanning a wide range of object sizes under varied contexts. Existing RSI backbones address the two challenges separately, either by adopting anisotropic strip kernels to model slender targets or by using isotropic large kernels to capture broader context. However, such isolated treatments lead to complementary drawbacks: the strip-only design can disrupt spatial coherence for regular-shaped objects and weaken tiny details, whereas isotropic large kernels often introduce severe background noise and geometric mismatch for slender structures. In this paper, we extend PKINet, and present a powerful and efficient backbone that jointly handles both challenges within a unified paradigm named Poly Kernel Inception Network v2 (PKINet-v2). PKINet-v2 synergizes anisotropic axial-strip convolutions with isotropic square kernels and builds a multi-scope receptive field, preserving fine-grained local textures while progressively aggregating long-range context across scales. To enable efficient deployment, we further introduce a Heterogeneous Kernel Re-parameterization (HKR) Strategy that fuses all heterogeneous branches into a single depth-wise convolution for inference, eliminating fragmented kernel launches without accuracy loss. Extensive experiments on four widely-used benchmarks, including DOTA-v1.0, DOTA-v1.5, HRSC2016, and DIOR-R, demonstrate that PKINet-v2 achieves state-of-the-art accuracy while delivering a 3.9× FPS acceleration compared to PKINet-v1, surpassing previous remote sensing backbones in both effectiveness and efficiency.
Xinhao Cai, Liulei Li, Gensheng Pei +3
Nanjing University of Science and Technology · State Key Lab. of Intelligent Manufacturing of Advanced Construction Machinery · Zhejiang University +1
Kernels measure similarity or correlation in tasks such as regression and classification. The Gaussian kernel, other names of which include squared exponential and radial basis function kernel, is one of the most popular in Gaussian process regression. We argue that the Gaussian kernel is best avoided and should never be used as a default. The argument rests on two results demonstrating that the Gaussian kernel is extremely brittle. First, the Gaussian kernel gives rise to a conditional variance that is unrealistically small. If the variance is used to quantify predictive uncertainty, catastrophic overconfidence is almost inevitable. Second, a small variance goes hand in hand with numerical ill-conditioning, so that to use the Gaussian kernel in practice requires tricks such as nugget terms that effectively modify the underlying regression or classification model. These problems are caused by the unnatural smoothness of the Gaussian kernel, a fact we are far from the first to take notice of. The problem is not the Gaussian form itself but the analyticity of the kernel: Our argument is more broadly that analytic kernels are best avoided. For stationary kernels analyticity is essentially equivalent to an exponential decay of the spectral density.
Toni Karvonen, Chris J. Oates
School of Engineering Sciences Lappeenranta–Lahti University of Technology LUT Yliopistonkatu 34, 53850 Lappeenranta, Finland · School of Mathematics, Statistics & Physics Newcastle University Newcastle upon Tyne, NE1 7RU, United Kingdom