Machine Learning · 12 of 20

DBSCAN tutorial

Groups points that are closely packed together into dense regions. Learn the concept, implement it from scratch, and apply it with scikit-learn.

Category
Unsupervised learning
Core idea
Nε(p) = {q : d(p,q) ≤ ε}, |Nε(p)| ≥ minPts
Typical use
Anomaly detection, spatial clustering, and arbitrary-shaped clusters.
Practical tool
scikit-learn

What is DBSCAN?

DBSCAN is a unsupervised learning technique. Groups points that are closely packed together into dense regions. Its central idea is summarized by Nε(p) = {q : d(p,q) ≤ ε}, |Nε(p)| ≥ minPts.

This guide connects theory to practice. You will trace the input, intermediate process, and output; decode the notation; run a dependency-light implementation; then repeat the workflow with scikit-learn and evaluate the result.

By the end of this tutorial, you will be able to

  • Explain when DBSCAN is appropriate and what assumptions it makes.
  • Read its mathematical notation or complexity statement without guessing what the symbols mean.
  • Follow and modify a from-scratch Python implementation.
  • Build a practical workflow with scikit-learn and choose useful evaluation checks.

Prerequisites and tools

You do not need an advanced software stack. Start with a recent Python environment and the fundamentals below, then install only the packages used by the practical example.

  • Comfort with Python functions, NumPy arrays, and basic descriptive statistics.
  • A clear distinction between training data, validation data, and untouched test data.
  • Familiarity with features, targets, preprocessing, and task-appropriate evaluation metrics.
LanguagePython 3
Primary libraryscikit-learn
Learning modeFrom scratch + library

How DBSCAN works

  1. 1

    Input

    Unlabelled numeric observations.

  2. 2

    Prepare and configure

    Check shapes, value ranges, ordering assumptions, missing values, and the parameters that control the algorithm’s behavior.

  3. 3

    Algorithm process

    Expand clusters from core points based on density.

  4. 4

    Output

    Cluster labels including noise points.

  5. 5

    Validate

    Report both cluster count and noise fraction while varying ε and min_samples.

Formula and intuition

Nε(p) = {q : d(p,q) ≤ ε}, |Nε(p)| ≥ minPts

Nε(p) is the ε-neighborhood of p; p is a core point when that neighborhood contains at least minPts observations.

The notation captures the main operation or complexity statement behind DBSCAN. Read it together with the symbol key above and the step-by-step process in this guide.

When to use DBSCAN

Anomaly detection, spatial clustering, and arbitrary-shaped clusters.

Learning tip

Change one parameter at a time in the interactive lesson, replay the animation, and connect the visible change to the input, process, and output described above.

DBSCAN from scratch in Python

This dependency-light example emphasizes the algorithm’s mechanics so each important step remains visible.

PythonEducational implementation
import numpy as np

rng = np.random.default_rng(0)
X = rng.uniform(0, 10, (100, 2))
eps, min_pts = 0.8, 5

# a point is a "core point" if at least min_pts other points lie within eps of it
n = len(X)
labels = np.full(n, -1)  # -1 means "noise" until a core point claims it
neighbors = [np.where(np.linalg.norm(X - X[i], axis=1) < eps)[0] for i in range(n)]
core = np.array([len(neighbors[i]) >= min_pts for i in range(n)])

cluster_id = 0
for i in range(n):
    if labels[i] != -1 or not core[i]:
        continue  # only core points can start a new cluster
    labels[i] = cluster_id
    queue = list(neighbors[i])
    while queue:  # flood-fill outward through density-connected neighbors
        j = queue.pop()
        if labels[j] == -1:
            labels[j] = cluster_id
            if core[j]:
                queue.extend(neighbors[j])  # keep expanding only from other core points
    cluster_id += 1

noise_points = np.sum(labels == -1)
How to use this example

Run it once unchanged, inspect the output, and then alter one input or parameter. The compact implementation is designed for learning; use the tested library workflow below for real projects.

Build DBSCAN with scikit-learn

DBSCAN labels density-connected clusters and noise directly, without requiring a cluster count.

Recommended environment

JupyterLab, Google Colab, or VS Code

Create an isolated virtual environment for a local project, or paste the cells into a hosted notebook. Pin package versions before deploying a reproducible application.

Install the required package python -m pip install scikit-learn
  1. Prepare the data.Unlabelled numeric observations. Validate its shape, type, range, and ordering before training or execution.
  2. Configure the algorithm.Begin with explicit, conservative parameters and a fixed random seed whenever the library supports one.
  3. Fit or execute.Expand clusters from core points based on density.
  4. Inspect the result.Cluster labels including noise points. Then apply the evaluation checks in the next section.
Pythonscikit-learn workflow
import numpy as np
from sklearn.cluster import DBSCAN
from sklearn.datasets import make_moons
from sklearn.metrics import silhouette_score
from sklearn.preprocessing import StandardScaler

X, _ = make_moons(n_samples=700, noise=0.08, random_state=7)
X = StandardScaler().fit_transform(X)
labels = DBSCAN(eps=0.22, min_samples=8).fit_predict(X)
cluster_labels = set(labels) - {-1}
mask = labels != -1
print("clusters:", len(cluster_labels))
print("noise points:", np.sum(~mask))
if len(cluster_labels) > 1:
    print("silhouette:", silhouette_score(X[mask], labels[mask]))

API details and version-specific options: official scikit-learn reference →

How to evaluate the result

A successful run is not enough. Evaluate the output against the intended use, compare it with a simple baseline, and preserve a genuinely unseen test case whenever the task involves learned parameters.

  • Report both cluster count and noise fraction while varying ε and min_samples.
  • Use a k-distance plot as a diagnostic starting point for ε.
Reproducibility check

Record the data version, package versions, parameters, random seeds, and evaluation procedure. Re-run the same workflow before publishing a benchmark or deploying a model.

Pitfalls and how to avoid them

These failure modes are common in tutorials and production systems. Treat them as review questions, not just after-the-fact debugging advice.

  • A single ε struggles when clusters have very different densities.
  • Feature scaling and the distance metric directly define neighborhoods.
  • Silhouette score can be misleading when many observations are marked as noise.

Before using DBSCAN in a project

  • Define the prediction or discovery objective before selecting the algorithm.
  • Split data before fitting preprocessing and tune only inside cross-validation.
  • Record data versions, random seeds, features, hyperparameters, and evaluation metrics.
  • Compare against a simple baseline and inspect errors by meaningful subgroups.
  • Monitor input drift and real-world performance after deployment.

Official documentation and next steps

Primary software reference scikit-learn

Use the official documentation to confirm supported parameters, current defaults, input requirements, and version changes.

Read official documentation →

This guide is an educational introduction, not a substitute for domain validation. For consequential applications, review the source documentation, test against representative data, and involve a subject-matter expert.

See DBSCAN in motion

Open the interactive lesson to adjust parameters, scrub through the process, replay the animation, and compare the explanation with the Python code.

Launch visualization →