Geometria Computacional - PDF
Download "Geometria Computacional"
Vinícius de Andrade Leal
1 Universidade Federal do Espírito Santo Centro de Ciências Agrárias CCA UFES Departamento de Computação Universidade Federal do Espírito Santo CCA UFES Geometria Computacional Tópicos Especiais em Programação Site:
2 Geometria Computacional É Importante em aplicações, como Gráficos computacionais, Robótica, Design, etc. No Mundo real: A maioria dos objetos não é feita por linhas que tendem ao infinito; Já os objetos no computador: Em sua maioria representam a geometria como arranjos de segmentos de linha. Curvas fechadas ou formas são representadas por coleções ordenadas de segmentos de linha ou polígonos (discretos). Geometria computacional: A geometria de segmentos de linhas discretas e polígonos. É um divertido e interessante assunto, mas não tipicamente ensinado em cursos universitários necessários. O aluno que aprende um pouco de geometria computacional tem vantagem e uma janela aberta para uma área fascinante de algoritmos ainda sob investigação nos dias de hoje. 2
3 Segmentos de Linha Um segmento de linha s é uma reta formada por dois pontos: A = (1 ; 1) B = (3 ; 2,5) Assim, podem ser representados por um par de pontos: typedef struct { point p1,p2; } segmento; 3
4 Segmentos de Linha e Interseção Primitiva mais importante na geometria sobre segmentos: Testar se um determinado par deles se cruzam. Porém: Dois segmentos podem ser linhas paralelas, o que significa que não se cruza. Eles podem se cruzar nas extremidades, ou no meio do segmento. Isso pode ser um verdadeiro problema. Leia qualquer especificação do problema com cuidado para ver se ele não terá linhas paralelas ou segmentos sobrepostos. No entanto, é melhor fazer o programa preparado para lidar com eles. 4
5 Segmentos de Linha e Interseção Nós utilizaremos nossas rotinas de interseção para encontrar se a interseção de um ponto existe. bool point_in_box(point p, point b1, point b2){ return } (p.x >= min(b1.x,b2.x)) && (p.x <= max(b1.x,b2.x)) && (p.y >= min(b1.y,b2.y)) && (p.y <= max(b1.y,b2.y)) ; 5
6 bool segments_intersect(segment s1, segment s2) { line l1,l2; /* linhas contendo os seguimentos de entrada */ point p; /* ponto de intersecao */ Universidade Federal do Espírito Santo CCA UFES points_to_line(s1.p1, s1.p2, l1); points_to_line(s2.p1, s2.p2, l2); if (same_lineq(l1,l2)) /* segmentos sobrepostos ou disjuntos */ return( point_in_box(s1.p1,s2.p1,s2.p2) point_in_box(s1.p2,s2.p1,s2.p2) point_in_box(s2.p1,s1.p1,s1.p2) point_in_box(s2.p1,s1.p1,s1.p2) ); if (parallelq(l1,l2)) return false; if (! intersection_point(l1,l2,p) ) return false; } return point_in_box(p,s1.p1,s1.p2) && point_in_box(p,s2.p1,s2.p2); 6
7 Computação de Polígonos e Universidade Federal do Espírito Santo CCA UFES Polígonos: Ângulos São um conjunto de segmentos de linha sem interseção que fecham-se em um vértice, formando uma figura geométrica. O primeiro vértice do conjunto de segmentos é igual ao último. Suas arestas não possuem interseção: Somente pares de segmentos terão encontros em um de seus pontos (ponto final do segmento). São a estrutura básica para descrever formas no plano. Representação: Não listar os segmentos (arestas) do polígono. Fazer a listagem dos n vértices ordenados que formam as extremidades do polígono typedef struct { int n; /* número de pontos no polígono */ point p[maxpoly]; /* array de pontos do polígono */ } polygon; 7
8 Polígono Convexo A Região Convexa se e somente se todo segmento de reta reside na mesma região. A região A é convexa: Pois qualquer segmento de reta que for escolhido só tem pontos na mesma região A. O polígono B já não é convexo: Pois existe pelo menos um segmento de reta que tem extremidades na região B, mas tem pontos fora da região. Repare que os pontos M e N estão em B, mas O é um ponto fora da região. Neste caso a região B é chamada de Região Não- Convexa. Assim, todos os ângulos devem ter no máximo 180º ou π radianos. Universidade Federal do Espírito Santo CCA UFES 8
9 Polígonos Calcular o ângulo definido entre três pontos ordenados é um problema complicado. Nós podemos evitar a necessidade de saber ângulos reais na maioria dos algoritmos geométricos utilizando o atributo sentido anti-horário ccw(a, b, c). Esta rotina testa se o ponto C está posicionado à direita da linha dirigida que vai do ponto A para o ponto B. Caso positivo, o ângulo formado pela varredura de A até C em sentido anti-horário em torno B é agudo. Caso negativo, o ponto ou encontra-se à esquerda de AB ou os três pontos são colineares. 9
10 Polígonos bool ccw(point a, point b, point c) { } return (signed_triangle_area(a,b,c) > EPSILON); bool cw(point a, point b, point c) { } return (signed_triangle_area(a,b,c) < EPSILON); bool collinear(point a, point b, point c) { } return (fabs(signed_triangle_area(a,b,c)) <= EPSILON); 10
11 Casco Convexo (Convex Hulls) Universidade Federal do Espírito Santo CCA UFES O casco convexo de um conjunto Q de pontos no plano é o menor polígono convexo que os envolve. Todos os pontos de Q devem estar dentro do polígono ou sobre sua borda. Para visualizar esse polígono: Imagine que os pontos são pregos saindo do plano; Imagine uma linha esticada de modo a conter todos os pontos. Essa linha, após ser colocada prendendo-se aos pregos mais externos, delimitará o casco convexo dos pontos. 11
12 Varredura de Graham Considerado o primeiro algoritmo de Geometria Computacional (1972), há quase tantos algoritmos diferentes para envoltória convexa quanto há para ordenação. O algoritmo de Graham: Ordena os pontos do conjunto: em ordem angular ou da esquerda para a direita E depois insere os pontos na envoltória nesta ordem. 12
13 Varredura de Graham Algoritmo incremental (2D) Pontos são pré-ordenados de forma conveniente; Cada ponto é adicionado ao fecho convexo e testado. Precisamos de um ponto inicial p0 que faz parte do casco convexo de forma garantida. Solução: Tomar o ponto com menor coordenada x ou y como p 0. Na verdade, um ponto extremo em qualquer direção serve. 13
14 Varredura de Graham Os pontos restantes são ordenados de forma cíclica com respeito aos ângulos formados pelas retas p 0 p i. Pontos colineares são removidos nesse processo. 14
15 Varredura de Graham Cada ponto considerado deve estar à esquerda da aresta anteriormente computada (teste de orientação) Ou o ponto anterior faz uma curva para esquerda. 15
16 Varredura de Graham O fecho convexo é mantido como uma pilha de pontos. Percorre-se os pontos enquanto o ponto no topo da pilha não fizer uma curva para à esquerda, Quando se considera o novo ponto, ele é desempilhado. Em seguida estende-se a cadeia empilhando-se o novo ponto. 16
17 Varredura de Graham Universidade Federal do Espírito Santo CCA UFES 17
18 Triangulação Achar o perímetro de um polígono é fácil: Basta calcular o comprimento de cada lado usando a fórmula da distância euclidiana e somar tudo ao final. Computar a área de manchas irregulares já é mais difícil. A ideia é dividir o polígono em triângulos não sobrepostos e somar suas áreas. Esta operação é chamada de triangulação. 18
19 Triangulação Triangular um polígono convexo é fácil: basta conectar um dado vértice v a todos os outros (n - 1) vértices. Desta forma, o polígono é dividido em P triângulos usando arestas que não se cruzam e que se situam completamente dentro de P. Mas e quando o polígono não é convexo? 19
20 Triangulação Representação da triangulação: Listando as arestas ou, como é feito aqui, Com uma lista explícita dos índices dos vértices em cada triângulo: typedef struct { int n; int t[maxpoly][3]; } triangulation; O algoritmo de triangulação mais eficiente roda em tempo linear do número de vértices. O algoritmo mais fácil de implementar é baseado no corte da orelha (Algoritmo de Van Gogh). 20
21 Algoritmo de Van Gogh Uma orelha de um polígono P é um triângulo definido por um vértice v e seus vizinhos da esquerda e da direita (l e r), tal que o triângulo (v, l, r) se situa completamente dentro de P. Como lv e vr são arestas (limites) de P, o lado que define a orelha é rl. l v r 21
22 Algoritmo de Van Gogh rl deve se situar completamente no interior do polígono P. Para isso, lvr deve definir um ângulo não-reflexo. Nenhum outro segmento do polígono pode ser cortado por este lado, Senão um pedaço será retirado do triângulo. O fato importante é que todo polígono sempre contém uma orelha. Na verdade, duas no mínimo para n 4. 22
23 Algoritmo de Van Gogh Algoritmo: teste cada um dos vértices até achar uma orelha. Adiciona-se a aresta associada e corta-se a orelha, decrementando o número de vértices. O corte da orelha necessita testar se um dado ponto se situa dentro de um triângulo. O polígono restante deve também ter uma orelha, tal que continua-se a cortar até que apenas três vértices sejam mantidos, ficando um triângulo. 23
24 Algoritmo de Van Gogh Deve-se testar se um vértice que define uma orelha tem duas partes. Para o teste do ângulo: Utiliza-se novamente ccw/cw. Assume-se que os vértices estão ordenados em uma ordem anti-horária em volta do centro virtual. A reversão da ordem do polígono nos obriga a trocar o sinal no teste do ângulo. 24
25 Triangulação de um polígono via algoritmo de Van Gogh (corte da orelha), com triângulos rotulados na ordem da inserção (A-I ) 25
26 Computação das áreas Pode-se então computar a área de qualquer polígono: Basta triangular e somar as áreas de todos os triângulos formados. Somando as áreas com sinal dos triângulos Estes definidos por um ponto arbitrário p com cada segmento do polígono P, Chega-se na área de P, porque os triângulos com sinal negativo cancelam a área fora do polígono: n 2 Área(P)= 1 2 (x i y i+1 x i+1 y i ) i=0 26
28 Problemas Fáceis: Chainsaw Massacre; Useless Tile Packers; Intermediários: Trees on My Island; Difíceis: Nice Milk. 28
Técnicas de Detecção de Colisão para Jogos
Técnicas de Detecção de Colisão para Jogos por Leandro Silva Técnicas de Detecção de Colisão para Jogos por: Gustavo Russo Zanardo Esse artigo visa mostrar as principais técnicas para detecção de colisão
Geometria Computacional: Fecho convexo
Geometria Computacional: INF2604 Geometria Computacional Prof. Hélio Lopes lopes@inf.puc-rio.br sala 408 RDC O que é uma região convexa? Uma região é convexa se para qualquer par de pontos dessa região,
Árvores Binárias e Busca. Jeane Melo
Árvores Binárias e Busca Jeane Melo Roteiro Parte 1 Árvores Relação hierárquica Definição Formal Terminologia Caminhamento em Árvores Binárias Exemplos Parte 2 Busca seqüencial Busca Binária Grafos Conjunto
Interseção de Semiplanos
Interseção de Semiplanos Claudio Esperança Paulo Roma 1 Interseção de Semiplanos Problema consiste em construir a região convexa dada pela interseção de um conjunto de n semiplanos H = { h 1, h 2..., h