Resumo Técnico: Design Poliédrico com Interpolantes de n Lados Mistos

Declaração do Problema

O artigo aborda o desafio de gerar superfícies paramétricas suaves que interpolam os vértices de uma malha fechada com topologia arbitrária. Métodos existentes baseados em poliedros de controle (como subdivisão recursiva, splines generalizadas e abordagens de variedades) frequentemente produzem superfícies aproximativas; este trabalho foca em interpolar diretamente os pontos de entrada. A dificuldade principal reside no tratamento de vértices irregulares (aqueles com valência diferente de 4), que criam redes de controle não quadriláteras. Patches de produto tensorial padrão são insuficientes para estas regiões multilaterais, exigindo uma técnica especializada de parametrização e mistura para garantir continuidade suave ($C^k$) por toda a superfície.

Metodologia

1. Estrutura Geral

O método proposto constrói uma superfície através da mistura de interpolantes quadráticos locais de múltiplos lados. Assume-se que a entrada é uma malha fechada; se a malha contiver faces não quadriláteras, realiza-se uma divisão central (análoga a um passo de subdivisão de Catmull–Clark) para manter os vértices originais enquanto converte todas as faces em quadriláteros. Para cada quadrilátero na malha, a vizinhança de 1 anel de seus cantos define uma rede de controle quadrática de múltiplos lados. Quatro patches locais são gerados para interpolar estes pontos de controle e são subsequentemente misturados dentro da região paramétrica associada ao quadrilátero.

2. Malhas Regulares (Valência 4)

Para malhas regulares onde todos os vértices possuem valência 4, as redes de controle são quadriláteras. O método emprega patches de Bézier de produto tensorial quadrático.

• Definição do Interpolante: O patch $I(u, v)$ é definido usando polinômios de Bernstein $B^d_k(t)$. Os pontos de controle $P_{ij}$ são derivados dos pontos de canto ($C_i$), pontos de aresta ($E_i$) e do ponto central ($M$).
• Tratamento de Arestas: Os pontos de controle das arestas são ajustados para $\hat{E}_i = 2E_i - \frac{1}{2}(C_{i-1} + C_i)$ para garantir que a curva de Bézier quadrática interpole o ponto de aresta original $E_i$.
• Mistura (Blending): O surface $S(u, v)$ final para um quadrilátero é a soma ponderada de quatro interpolantes locais, cada um avaliado em um espaço paramétrico rotacionado. A função de mistura $\Phi(u, v)$ é o produto de dois blends de Hermite de 1 argumento ($k=2$), garantindo continuidade $C^2$ nas fronteiras.

3. Vértices Irregulares (Valência $n \neq 4$)

Para vérticos irregulares, a rede de controle é um polígono de $n$ lados. O artigo introduz patches Quadratic Generalized Bézier (QGB).

• Construção QGB: A superfície é definida sobre um polígono regular de $n$ lados usando coordenadas de Wachspress ($\lambda_i$). O interpolante combina contribuições das arestas de contorno (usando curvas de Bézier quadráticas) e um ponto de controle central $P_0$.
• Cálculo do Ponto Central: $P_0$ é calculado para garantir que o patch interpole o ponto central $M$. Para $n=4$, esta formulação reduz-se ao caso de produto tensorial; para $n=3$, corresponde a um patch de Bézier triangular quadrático.
• Desafio de Parametrização: Um mapeamento bilinear ingênuo do domínio do quadrilátero $[0, 1]^2$ para o domínio do polígono resulta em descontinuidades em partes adjacentes do mesmo interpolante.
• Mapeamento de Curvas Racionais: Para alcançar continuidade $C^\infty$ (exceto na origem singular), o artigo propõe um mapeamento baseado na interseção de dois feixes de curvas de Bézier racionais quadráticas.
  • As curvas $u$ e $v$ são definidas com pontos de controle específicos e pesos racionais.
  • Um termo exponencial não linear no peso do ponto de controle central é introduzido para lidar com domínios côncavos (especificamente para casos triangulares), efetivamente "puxando de volta" a curva próximo à origem.
  • A interseção destas curvas é computada via buscas de seção áurea aninhadas, um processo sugerido para ser pré-computado para eficiência.

4. Subpatches Triangulares

Para o caso específico de $n=3$, o patch QGB padrão carece de um ponto de controle central. Para interpolar o ponto central $M$, o método realiza a elevação de grau das fronteiras para cúbico e utiliza um triângulo de Bézier cúbico. Isso introduz um grau de liberdade extra ($P_{111}$) que é resolvido para satisfazer a restrição de interpolação em $M$.

5. Representações Alternativas

O artigo observa que outras representações multilaterais, como patches de Coons de Midpoint e Midpoint, são viáveis. No entanto, resultados experimentais indicaram pouca diferença visual nos modelos finais, enquanto os patches QGB propostos provaram ser computacionalmente mais simples e eficientes.

Principais Contribuições

• Superfícies Paramétricas Interpolantes: Uma mudança de esquemas de subdivisão aproximativos para uma construção que interpola estritamente os vértices da malha de entrada.
• Patches Quadratic Generalized Bézier (QGB): Uma extensão dos patches de Bézier generalizados para o caso quadrático, permitindo a interpolação de redes de controle de múltiplos lados.
• Parametrização de Curvas Racionais: Uma técnica de mapeamento de continuidade $C^\infty$ utilizando curvas de Bézier racionais para lidar com a parametrização de domínios irregulares (não quadriláteros), resolvendo os problemas de descontinuidade de mapeamentos mais simples.
• Tratamento Unificado de Topologia: Um framework que mistura perfeitamente patches quadráticos regulares com patches QGB de múltiplos lados para lidar com topologias de malha arbitrárias.

Resultados

O método foi testado em diversos modelos:

1. Toro (Malha Quadrada Regular): Demonstrou fluxo suave de linhas de isofote, embora a forma global tenha desviado de um toro perfeito devido às restrições de interpolação.
2. Modelo Trebol: Um objeto complexo com vértices de valência 3, 4, 5 e 6. As linhas de isofote confirmaram a conexão contínua e boa qualidade dos patches multilaterais. Mapas de curvatura média mostraram leve achatamento ao longo das arestas.
3. Icosaedro: Convertido para uma malha quadrática via subdivisão Catmull–Clark. A análise de curvatura mostrou valores maiores próximos a vértices de valência 5, mas flutuações pequenas fora disso, correspondendo às arestas da malha.
4. Modelo Bob: Um modelo complexo demonstrando a capacidade do método com mapeamento de ambiente, contorno e visualização de curvatura média.

Significância e Trabalho Futuro

O artigo afirma sua significância ao fornecer uma nova representação de superfície interpolante que mantém a suavidade através de topologias arbitrárias sem depender de limites de subdivisão. Os autores são modestos em suas afirmações, observando que, embora o método funcione bem para malhas fechadas, vários caminhos permanecem para pesquisas futuras:

• Generalização para malhas com bordas.
• Incorporação de restrições de contorno contínuas (interpolação de posição e derivada cruzada).
• Extensão com parâmetros de forma.
• Melhoria da qualidade da superfície através da interpolação de vetores normais.

O trabalho é apoiado pelo Fundo de Pesquisa Científica Húngaro (OTKA, No. 145970).