Os avanços tecnológicos têm permitido a obtenção de uma enorme quantidade de dados da superfície terrestre. Nesse contexto, a busca pela construção de um modelo digital de terreno fiel a partir desses dados tem ganhado grande importância devido a sua utilidade em diversos campos do conhecimento. Ao contrário da triangulação de Delaunay clássica, a triangulação de Delaunay de ordem superior pode trazer a flexibilidade necessária para a inclusão de critérios que possam melhorar a representação de terrenos usando redes triangulares irregulares (TINs), de acordo com sua utilização final. Este trabalho tem o objetivo de implementar um algoritmo para a criação de triangulações de Delaunay de ordem superior tomando como base os estudos feitos por Rodríguez e Silveira (2016) e Wang et. al. (2001). Foi implementado um algoritmo baseado em otimização local utilizando a linguagem C++ e a Computational Geometry Algorithms Library (CGAL). Foram utilizados um terreno real, o Kinzel Springs no Tennessee, EUA, o mesmo utilizado por Rodríguez e Silveira (2016), e uma superfície sintética definida por f(x,y) = xy, com x e y entre 0 e 1. Devido à facilidade de implementação e sua relação com a suavidade da superfície do terreno, escolheu-se como critério de otimização o ângulo entre os vetores normais às faces (ABN). Conforme esperado, as TINs geradas pela triangulação de ordem superior para o terreno sintético apresentaram um menor erro quadrático médio, isto é, uma melhor fidelidade ao terreno natural em relação à triangulação de Delaunay clássica. Contudo, houve divergência dos resultados para o Kinzel Spring, quando comparado àqueles obtidos por Rodríguez e Silveira (2016). Provavelmente, esta divergência está relacionada à diferença nos pontos amostrados a partir das imagens de satélite desta região.