Asymptotic analysis of higher-order perturbations of the Perona--Malik functional

1. Il Problema e il Contesto

Il lavoro si concentra sull'analisi asintotica (tramite $\Gamma$-convergenza) di perturbazioni singolari di ordine superiore del funzionale di Perona-Malik, fondamentale nell'elaborazione delle immagini.
Il funzionale di Perona-Malik originale, $\int_\Omega \log(1 + |\nabla u|^2) dx$, è non convesso e formalmente mal posto: i suoi minimizzanti tendono a sviluppare singolarità. Sebbene l'involucro semicontinuo inferiore sia banale (identicamente nullo), la struttura del funzionale suggerisce un comportamento interessante per le sequenze minimizzanti: effetto di smoothing per gradienti piccoli e formazione di interfacce per gradienti grandi.

L'obiettivo è studiare il limite asintotico di una versione regolarizzata di ordine $k$-esimo, progettata per bilanciare gli effetti di volume (bulk) e di superficie (interfacce). Il funzionale considerato è:
$$F_\varepsilon(u) = \int_\Omega \frac{1}{\varepsilon |\log \varepsilon|} \log(1 + \varepsilon |\log \varepsilon| |\nabla u|^2) dx + \varepsilon^{2k-1} \int_\Omega |\nabla^k u|^2 dx$$
definito su $H^k(\Omega)$. La sfida principale risiede nel fatto che il termine logaritmico ha una crescita sub-lineare, permettendo energie finite su funzioni con discontinuità di salto, mentre il termine di regolarizzazione di ordine $k$ previene la concentrazione su derivate di tipo Cantor.

2. Metodologia

Gli autori adottano un approccio rigoroso basato sul calcolo delle variazioni e sulla teoria della $\Gamma$-convergenza. La metodologia si articola in diversi passaggi chiave:

• Riduzione alla dimensione uno: Utilizzando un argomento di "slicing" (fette) e densità, il problema multidimensionale viene ridotto allo studio di funzioni unidimensionali su $(0,1)$. Questo permette di analizzare il comportamento delle interfacce in modo più diretto.
• Stime di Coercività e Limiti Inferiori: Viene dimostrata una disuguaglianza di limite inferiore che lega il termine logaritmico a un potenziale quadratico troncato. Nello specifico, si mostra che per $\varepsilon$ sufficientemente piccolo:
  $$\frac{1}{\varepsilon |\log \varepsilon|} \log(1 + \varepsilon |\log \varepsilon| z^2) \geq (1-\eta) \max\left\{ \varepsilon^{1-2p_\varepsilon} z^2, \frac{1}{\varepsilon} \right\}$$
  dove $p_\varepsilon \to 0$. Questa stima debole è cruciale per gestire la bassa coercività rispetto ai casi precedenti (es. potenziale quadratico troncato standard).
• Tecniche di Interpolazione Locale: Per determinare la densità di energia di superficie, gli autori estendono le tecniche di interpolazione locale sviluppate da Solci (2023) per potenziali troncati. La novità risiede nell'adattare queste tecniche a un contesto in cui il coefficiente del termine quadratico cresce solo leggermente più lentamente di $\varepsilon$.
• Analisi dei Profili Ottimali: Il cuore della dimostrazione consiste nello studiare il problema di ottimizzazione unidimensionale che determina il costo energetico di un salto. Si considera un intervallo di transizione di lunghezza $T$ e si minimizza la somma della lunghezza e dell'energia di transizione data dalla derivata $k$-esima, soggetto a condizioni al contorno omogenee sulle derivate fino all'ordine $k-1$.

3. Risultati Principali

Il Teorema di $\Gamma$-convergenza

Il risultato centrale (Teorema 2.1) stabilisce che il limite $\Gamma$ del funzionale $F_\varepsilon$ nella topologia $L^1$ (e in misura) è un funzionale di discontinuità libera definito su $SBV(\Omega)$ (Funzioni a Variazione Limitata Speciali):
$$F(u) = \int_\Omega |\nabla u|^2 dx + m_k \int_{S_u \cap \Omega} |u^+ - u^-|^{1/k} d\mathcal{H}^{N-1}$$
Dove:

• Il primo termine è l'energia di Dirichlet (parte di volume).
• Il secondo termine è l'energia di superficie, proporzionale all'ampiezza del salto elevata alla potenza $1/k$.
• $S_u$ è l'insieme delle discontinuità di $u$.
• $m_k$ è una costante positiva data da un problema di ottimizzazione unidimensionale:
  $$m_k = \inf_{T>0} \min \left\{ T + \int_0^T |v^{(k)}|^2 dt : v \in H^k(0,T), v(0)=0, v(T)=1, v^{(\ell)}(0)=v^{(\ell)}(T)=0 \forall \ell \in \{1, \dots, k-1\} \right\}$$

Caratteristiche Chiave dei Risultati

1. Potenza del Salto ($1/k$): La presenza dell'esponente $1/k$ nell'energia di superficie è una conseguenza diretta della scalatura del termine di regolarizzazione di ordine $k$. Questo conferma una congettura precedente di Gobbino e Picenni.
2. Condizioni al Contorno Omogenee: Un contributo tecnico significativo è la dimostrazione che, per determinare il profilo ottimale che minimizza l'energia di salto, le derivate di ordine $\ell \geq 2$ devono annullarsi ai bordi dell'intervallo di transizione. Questo è ottenuto dimostrando che, per salti sufficientemente grandi, gli intervalli sopra la soglia possono essere estesi leggermente in modo che tutte le derivate siano piccole ai bordi.
3. Indipendenza dalla Crescita Logaritmica: La costante $m_k$ risulta essere la stessa che si otterrebbe per un potenziale quadratico troncato, indicando che la crescita logaritmica non influenza il valore della costante di superficie, ma solo la struttura del limite di volume.

Applicazione a un'altra Scalatura

Gli autori applicano i risultati ottenuti per determinare il limite di una diversa scalatura del funzionale di Perona-Malik (analizzata da Gobbino e Picenni per $k=2$):
$$\tilde{F}_\varepsilon(u) = \int_\Omega \frac{1}{2\varepsilon |\log \varepsilon|} \log(1 + |\nabla u|^2) dx + \varepsilon^{2k-1} \int_\Omega |\nabla^k u|^2 dx$$
In questo caso, il limite $\Gamma$ è puramente superficiale:
$$\tilde{F}(u) = m_k \int_{S_u \cap \Omega} |u^+ - u^-|^{1/k} d\mathcal{H}^{N-1}$$
definito su funzioni in $SBV$ con gradiente nullo quasi ovunque. Questo risultato è rilevante per lo studio del fenomeno dello "staircasing" (formazione di gradini) nell'elaborazione delle immagini.

4. Significato e Impatto

• Teoria Matematica: Il lavoro estende la teoria dei funzionali di discontinuità libera a perturbazioni di ordine superiore con potenziali non convessi a crescita sub-lineare. Fornisce una caratterizzazione rigorosa del comportamento asintotico di modelli di regolarizzazione di ordine superiore.
• Elaborazione delle Immagini: La comprensione del limite $\Gamma$ aiuta a spiegare matematicamente fenomeni osservati numericamente, come la formazione di bordi netti e il fenomeno dello "staircasing". La dipendenza dell'energia di salto dalla potenza $1/k$ suggerisce come l'ordine della regolarizzazione influenzi la morfologia delle discontinuità ottimali.
• Innovazione Tecnica: La capacità di lavorare in un quadro di "bassa coercività" (dove il termine logaritmico non fornisce un controllo forte come un potenziale quadratico) e di derivare condizioni al contorno omogenee per derivate di ordine superiore rappresenta un avanzamento metodologico significativo rispetto alla letteratura esistente.

In sintesi, il paper fornisce una mappa completa del comportamento asintotico delle perturbazioni di Perona-Malik di ordine superiore, collegando la struttura analitica del funzionale alla geometria delle soluzioni limite, con implicazioni sia teoriche che applicative.