Tengo una pregunta sobre el ajuste cuadrático de un conjunto de puntos y las normales correspondientes (o, de forma equivalente, tangentes). La exploración de las superficies cuadráticas para señalar datos está bien explorada. Algunos trabajos son los siguientes:
Adaptación directa con restricciones de tipo de superficies cuadráticas , James Andrews, Diseño y aplicaciones asistidas por computadora Carlo H. Sequin , 10 (a), 2013, bbb-ccc
Adaptación algebraica de superficies cuadráticas a datos , I. Al-Subaihi y GA Watson , Universidad de Dundee
La adaptación a contornos proyectivos también está cubierta por algunos trabajos, como este .
De todos estos trabajos, creo que el método de Taubin para la adaptación cuadrática es bastante popular:
- G. Taubin, "Estimación de curvas planas, superficies y curvas espaciales no planas definidas por ecuaciones implícitas, con aplicaciones a la segmentación de imagen de borde y rango ", IEEE Trans. PAMI, vol. 13, 1991, pp1115-1138.
Déjame resumir brevemente. Una cuadrática se puede escribir en forma algebraica:
Ajuste algebraico En principio, nos gustaría resolver los parámetros que minimizan la suma de las distancias geométricas cuadradas entre los puntos y la superficie cuadrática. Desafortunadamente, resulta que este es un problema de optimización no convexo sin soluciones analíticas conocidas. En cambio, un enfoque estándar es resolver un ajuste algebraico, es decir, resolver los parámetros que minimizan:
Observe que tal minimización directa produciría la solución trivial con en el origen. Esta pregunta ha sido estudiada ampliamente en la literatura. Una resolución que funciona bien en la práctica es el método de Taubin (citado anteriormente), que introduce la restricción:
Esto se puede resolver de la siguiente manera: Sea:
donde los subíndices denotan las derivadas. La solución viene dada por la descomposición generalizada de Eigen, . El vector de parámetro de mejor ajuste es igual al vector propio correspondiente al valor propio más pequeño.
Pregunta principal En muchas aplicaciones, las normales de la nube de puntos están disponibles (o calculadas). Las normales del quadric también se pueden calcular diferenciando y normalizando la superficie implícita:
Sin embargo, el método de Taubin utiliza solo la geometría del punto, y no el espacio tangente. Y no conozco muchos métodos, que son adecuados para ajustar cuádricos de modo que las tangentes del cuadrático también coincidan con las tangentes de la nube de puntos subyacente. Estoy buscando posibles extensiones del método anterior, o cualquier otra para cubrir estos derivados de primer orden.
Lo que me gustaría lograr quizás se aborde parcialmente en espacios de dimensiones inferiores, con tipos de superficie (curva) más primitivos. Por ejemplo, ajustando líneas a los bordes de la imagen, teniendo en cuenta la información del gradiente, aquí se trata . El ajuste de planos (un tipo simple de cuadrático) a nubes 3D es muy común ( enlace 1 ) o las esferas o cilindros de ajuste se pueden ajustar a conjuntos de puntos orientados ( enlace 2 ). Entonces, lo que me pregunto es algo similar, pero la primitiva ajustada es un cuádruple.
También agradecería el análisis del método propuesto, como por ejemplo:
- ¿Cuál es el número mínimo de puntos orientados requeridos?
- ¿Cuáles son los casos degenerados?
- ¿Se puede decir algo sobre la robustez?
Actualización : me gustaría presentar una dirección a seguir. Formalmente, lo que deseo lograr: