Sintesi Tecnica: Design Poliedrico con Interpolanti n-laterali Miscelati

Definizione del Problema

Il documento affronta la sfida di generare superfici parametriche lisce che interpolino i vertici di una mesh chiusa con topologia arbitraria. Mentre i metodi esistenti basati su poliedri di controllo (come la suddivisione ricorsiva, gli spline generalizzati e gli approcci manifold) producono spesso superfici approssimative, questo lavoro si concentra sull'interpolazione diretta dei punti di input. La difficoltà principale risiede nella gestione dei vertici irregolari (quelli con valenza diversa da 4), che creano reti di controllo non quadrilaterali. I normali patch prodotto tensoriale sono insufficienti per queste regioni multi-laterali, rendendo necessaria una parametrizzazione e una tecnica di miscelazione (blending) specializzate per garantire la continuità $C^k$ su tutta la superficie.

Metodologia

1. Framework Generale

Il metodo proposto costruisce una superficie miscelando interpolanti quadratici locali multi-laterali. Si assume che l'input sia una mesh chiusa; se la mesh contiene facce non quadrilaterali, viene eseguita una divisione centrale (analogamente a un passaggio di suddivisione di Catmull–Clark) per mantenere i vertici originali convertendo tutte le facce in quadrilateri. Per ogni quad nella mesh, il vicinato 1-ring dei suoi angoli definisce una rete di controllo quadratica multi-laterale. Vengono generati quattro patch locali per interpolare questi punti di controllo, che vengono successivamente miscelati all'interno della regione parametrica associata al quad.

2. Mesh Regolari (Valenza 4)

Per le mesh regolari dove tutti i vertici hanno valenza 4, le reti di controllo sono quadrilaterali. Il metodo impiega patch di Bézier prodotto tensoriale quadratici.

• Definizione dell'Interpolante: Il patch $I(u, v)$ è definito utilizzando i polinomi di Bernstein $B^d_k(t)$. I punti di controllo $P_{ij}$ sono derivati dai punti di angolo ($C_i$), dai punti di spigolo ($E_i$) e dal punto centrale ($M$).
• Gestione degli Spigoli: I punti di controllo degli spigoli sono regolati a $\hat{E}_i = 2E_i - \frac{1}{2}(C_{i-1} + C_i)$ per garantire che la curva di Bézier quadratica interpoli il punto di spigolo originale $E_i$.
• Miscelazione (Blending): La superficie finale $S(u, v)$ per un quad è la somma pesata di quattro interpolanti locali, ciascuno valutato in uno spazio parametrico locale ruotato. La funzione di blending $\Phi(u, v)$ è il prodotto di due blend di Hermite a 1 argomento ($k=2$), garantendo la continuità $C^2$ ai bordi.

3. Vertici Irregolari (Valenza $n \neq 4$)

Per i vertici irregolari, la rete di controllo è un poligono $n$-laterale. Il documento introduce i patch di Bézier Generalizzati Quadratici (QGB).

• Costruzione QGB: La superficie è definita su un poligono regolare $n$-laterale utilizzando coordinate di Wachspress ($\lambda_i$). L'interpolante combina i contributi degli spigoli del bordo (utilizzando curve di Bézier quadratiche) e un punto di controllo centrale $P_0$.
• Calcolo del Punto Centrale: $P_0$ è calcolato per garantire che il patch interpoli il punto centrale $M$. Per $n=4$, questa formulazione si riduce al caso prodotto tensoriale; per $n=3$, corrisponde a un patch di Bézier triangolare quadratico.
• Sfida di Parametrizzazione: Una mappatura bilineare ingenua dal dominio del quad $[0, 1]^2$ al dominio del poligono risulta in discontinuità nelle parti adiacenti dello stesso interpolante.
• Mappatura di Curve Razionali: Per ottenere una continuità $C^\infty$ (eccetto l'origine singolare), il documento propone una mappatura basata sull'intersezione di due fasci di curve di Bézier razionali quadratiche.
  • Le curve $u$ e $v$ sono definite con specifici punti di controllo e pesi razionali.
  • Un termine esponenziale non lineare nel peso del punto di controllo centrale viene introdotto per gestire i domini concavi (specificamente per i casi triangolari), "tirando indietro" efficacemente la curva vicino all'origine.
  • L'intersezione di queste curve viene calcolata tramite ricerche a sezione aurea annidate, un processo suggerito di essere precomputato per l'efficienza.

4. Sub-patch Triangolari

Per il caso specifico di $n=3$, lo standard patch QGB manca di un punto di controllo centrale. Per interpolare il punto centrale $M$, il metodo eleva il grado dei bordi a cubico e utilizza un triangolo di Bézier cubico. Ciò introduce un ulteriore grado di libertà ($P_{111}$) che viene risolto per soddisfare il vincolo di interpolazione in $M$.

5. Rappresentazioni Alternative

Il documento nota che altre rappresentazioni multi-laterali, come i patch di Coons Midpoint e Midpoint, sono praticabili. Tuttavia, i risultati sperimentali hanno indicato poca differenza visiva nei modelli finali, mentre i patch QGB proposti si sono dimostrati computazionalmente più semplici ed efficienti.

Contributi Chiave

• Superfici Parametriche Interpolanti: Uno spostamento dai metodi di suddivisione approssimativi verso una costruzione che interpola strettamente i vertici della mesh di input.
• Patch di Bézier Generalizzati Quadratici (QGB): Un'estensione dei patch di Bézier generalizzati al caso quadratico, che permette l'interpolazione di reti di controllo multi-laterali.
• Parametrizzazione di Curve Razionali: Una nuova tecnica di mappatura a continuità $C^\infty$ utilizzando l'intersezione di curve di Bézier razionali per gestire la parametrizzazione di domini irregolari (non quadrilateri), risolvendo i problemi di discontinuità delle mappature più semplici.
• Gestione Unificata della Topologia: Un framework che miscela fluidamente patch quadrati regolari con patch QGB multi-laterali per gestire topologie di mesh arbitrarie.

Risultati

Il metodo è stato testato su diversi modelli:

1. Toro (Mesh quad regolare): Ha dimostrato un flusso fluido delle linee di isofota, sebbene la forma globale deviasse da un toro perfetto a causa dei vincoli di interpolazione.
2. Modello Trebol: Un oggetto complesso con vertici di valenza 3, 4, 5 e 6. Le linee di isofota hanno confermato la connessione continua e la buona qualità dei patch multi-laterali. Le mappe di curvatura media hanno mostrato una leggera piattezza lungo gli spigoli.
3. Icosaedro: Convertito in una mesh quad tramite suddivisione di Catmull–Clark. L'analisi della curvatura ha mostrato valori più elevati vicino ai vertici a valenza 5, ma fluttuazioni altrimenti piccole corrispondenti agli spigoli della mesh.
4. Modello Bob: Un modello complesso che dimostra la capacità del metodo con environment mapping, contornatura e visualizzazione della curvatura media.

Significato e Lavoro Futuro

Il documento rivendica l'importanza di fornire una nuova rappresentazione di superficie interpolante che mantenga la fluidità attraverso topologie arbitrarie senza fare affidamento sui limiti di suddivisione. Gli autori sono modesti nelle loro affermazioni, notando che, sebbene il metodo funzioni bene per mesh chiuse, rimangono diverse strade per la ricerca futura:

• Generalizzazione a mesh con bordi.
• Incorporazione di vincoli di bordo continui (interpolazione di posizione e derivata incrociata).
• Estensione con parametri di forma.
• Miglioramento della qualità della superficie attraverso l'interpolazione dei vettori normali.

Il lavoro è supportato dal Fondo per la Ricerca Scientifica Ungherese (OTKA, No. 145970).