D. V. SEGUNDO, M. A. R.

SCHMIDT - REEC – Revista Eletrônica de Engenharia Civil Vol 11 - nº 1 ( 2016)

67

AJUSTAMENTO DE POLIGONAIS DE CONTROLE PARA LOCAÇÃO DE

OBRAS DE ENGENHARIA
Adjustments of control polygonals for civil building location
Denilson Viegas Segundo 1, Marcio Augusto Reolon Schmidt 2
Recebido em 16 de fevereiro de 2015; recebido para revisão em 20 de junho de 2015; aceito em 17 de agosto de
2015; disponível on-line em 02 de fevereiro de 2016.
RESUMO: A realização de um levantamento topográfico para fins de locação de
grandes obras de engenharia consiste em determinar com precisão adequada as
feições chave de projeto para que seja possível a execução conforme requisitos de
projeto. Os pilares, eixos de rodovias, alinhamento de represas e qualquer outro
ponto de ancoragem da superestrutura precisa estar na posição correta no esforço
de manter a execução da obra conforme o projeto de implantação. Entretanto, as
condições de campo como variações de umidade, temperatura, pressão, problemas
com equipamentos topográficos e de operação dão origem a diversos erros e
contribuem para a degradação da qualidade do processo de levantamento e locação
PALAVRAS CHAVE: topográfica. Por isso, esta pesquisa avalia as metodologias topográficas propostas
NBR13.133 para o controle de qualidade das coordenadas dos pontos das poligonais
Ajustamento de usadas como controle e, a partir das quais, essas grandes obras serão locadas no
observações; terreno para sua execução. Os resultados mostram que a avaliação por mínimos
Modelagem de erros quadrados de uma poligonal por caminhamento com diagonais internas ao polígono
topográficos; convexo apresenta maior rigidez na estrutura e minimiza os erros de leituras em
Locação de obra; campo, ao passo que o uso de leituras angulares apresentou os piores resultados.
Controle de estruturas ABSTRACT: The realization of a topographic survey for location of big engineering
geodésicas. building should be run with adequate precision the key features to be able to run
according to the project requirements. The pillars, road axes, dams alignment and
any other anchor point of the superstructure must be in the correct position in an
KEYWORDS: effort to maintain the execution of the work as the deployment project. However,
Adjustment computing; the field conditions such as humidity variations, temperature, pressure, problems
with topographical equipment and operation give rise to many errors and contribute
Topographic error to the degradation of its quality of the topographic location. Therefore, this research
modelling; evaluates topographical methodologies suggest in NBR13.133 proposals for the
Civil building location; quality control of the coordinates of the points of the polygonal used as control, and
Geodetic structures from which, these building works will be located for its implementation. The
control. evaluation results show that the least-squares method of a traditional polygonal
with the internal diagonal of the convex polygon has increased rigidity in structure
and minimizes reading errors in the field, while the use of angular readings showed
the worst results.
* Contato com os autores:
1 e-mail : denilsonsegundo@hotmail.com ( D. V. Segundo )
Discente do curso de Engenharia Civil da Universidade Federal de Uberlândia (UFU).
2 e-mail : marcioschmidt@feciv.ufu.br ( M. A. R. Schmidt )

Professor Dr da Faculdade de Engenharia Civil da Universidade Federal de Uberlândia (UFU ).

ISSN: 2179-0612 D.O.I. 10.5216/reec.V11i1.33819 © 2016 REEC - Todos os direitos reservados.

D. V. SEGUNDO, M. A. R. SCHMIDT - REEC – Revista Eletrônica de Engenharia Civil Vol 11 - nº 1 ( 2016)
68

1. INTRODUÇÃO Admitindo a existência de erros

randômicos e observações superabundantes,
A realização de um levantamento utiliza-se o Método dos Mínimos Quadrados
topográfico consiste em determinar com precisão (M.M.Q) para a determinação da estimativa única
coordenadas polares, através ângulos e distâncias, dos parâmetros ajustados (coordenadas das
que naturalmente e inevitavelmente contém erros feições) e avaliação da sua variância (resíduos). O
de medida que podem ser de natureza sistemática, princípio fundamental do M.M.Q, é que a soma dos
grosseira e aleatória (Comastri, 1994). Todo o quadrados dos resíduos seja minimizado (Gemael,
levantamento envolve a obtenção de medições 1994; Dalmolin 2002). Este método, constante na
com procedimentos adequados e rigorosos em NBR 13133 (ABNT, 1994), permite estimar a
campo. As redes de triangulações e trilaterações distribuição de erros e resíduos das coordenadas
topográficas são essenciais para fornecer dos pontos da poligonal ou da rede. O MMQ, em
coordenadas de pontos na superfície terrestre com seus diferentes métodos de ajustamento, permite
alta precisão, necessárias para o projeto, definir a qualidade e o grau de confiança para as
construção e monitoramento de obras importantes coordenadas e medidas envolvidas.
de engenharia como pontes, túneis ou barragens Por isto, esta pesquisa apresenta uma
(MELO et al, 2012). avaliação do controle de propagação de variâncias
Mesmo com as tecnologias disponíveis em poligonais fechadas por caminhamento,
como o posicionamento GPS-RTK, durante a trilateração e triangulação, com variação das suas
construção de grandes obras de engenharia como observações, na estimativa dos resíduos e
pontes, rodovias, túneis ou barragens, as estações indicação da mais adequada para utilização como
totais ainda são os equipamentos com mais ampla densificação do referencial nos trabalhos de
utilização. Nestas situações, as técnicas locação de obra.
empregadas para determinação das poligonais de
referência e as locações são as mesmas de muitos 2. REVISÃO BIBLIOGRÁFICA
anos atrás. São utilizadas as poligonais por
caminhamento e redes por trilateração e por Podem ser utilizadas diferentes
triangulação. metodologias de levantamento para a
Independente de qual técnica seja determinação das poligonais de apoio que servirão
utilizada, a determinação das coordenadas dos seus como referência para a locação das obras de
vértices deve ter precisão compatível com as engenharia. As técnicas topográficas para
especificações de projeto de forma que seja estabelecimento de pontos de alta precisão com
possível o uso destes pontos como apoio à tarefa uso de estações totais usadas são as poligonais por
de locação de obras. Segundo Carvalho e Leandro caminhamento, por triangulação e por trilateração.
(2007), os erros de observação não permitem que Na poligonal por caminhamento, as
haja um fechamento geométrico adequado da distâncias e direções são determinadas entre os
poligonal e, portanto, é necessário que se faça um vértices da poligonal ocupando-se os vértices
controle da qualidade dos erros nas observações de sucessivamente com a estação total, ou teodolito,
ângulos e distâncias. Estes erros ocorrem devido às e determinando os ângulos relativos entre si. O
falhas humanas, imperfeições do equipamento e procedimento é repetido até que todos os vértices
por fatores atmosféricos. Por estas razões é tenham sido ocupados. Trata-se da forma mais
imprescindível a determinação da qualidade comum de densificação de referencial utilizada em
posicional dos vértices para se estimar a levantamentos topográficos, entretanto o número
propagação de variância para os pontos a serem de equações obtidas é igual ao de variáveis. A
locados na obra. A ferramenta indicada para esta inserção de equações advindas da repetibilidade
análise é o ajustamento de observações. das leituras, incorre no estabelecimento de
D. V. SEGUNDO, M. A. R. SCHMIDT - REEC – Revista Eletrônica de Engenharia Civil Vol 11 - nº 1 ( 2016)
69
equações linearmente dependentes e tornar a função dos modernos distanciometros a laser das
solução do sistema impossível (Gemal 1994, pag estações totais. O princípio é formar triângulos que
246). Isto pode gerar um sistema indeterminado atendam ao critério de fechamento dos lados a
devido a singularidade da matriz. Por isto a partir da determinação do comprimento dos lados.
repetibilidade das leituras é utilizada para melhoria A Figura 1 apresenta a configuração dos métodos
das observações conforme as classes IIP, IIIP e IVP de trilateração (a) e triangulação (b) para a
da NBR13.133 (ABNT, 1994). Para haver determinação de coordenadas na planimetria.
superabundância de equações de forma a atender Na Figura 1a são realizadas as medidas
o M.M.Q deve-se, portanto, fixar as coordenadas das distâncias (l1 a l9) entre os pontos de forma que
de um ou mais pontos. as equações resultantes sejam não correlacionadas.
A triangulação consiste num método de Os pontos EPS4 e EPS7 foram fixados. Na figura 1b
observação de direções a partir de um alinhamento são apresentadas as leituras angulares da
de coordenadas com a distância entre os pontos triangulação (medidas α1 a α6) e a base l1 entre os
deste alinhamento conhecida. Esta base, ou seja, a pontos fixados P1 e P3.
distância tomada entre dois pontos fixados, apesar Existem problemas específicos no uso de
de usualmente ser pequena, é usada para triangulação ou trilateração isoladamente, como a
determinar a escala da rede, pois a partir dela são influência da configuração da rede em si nos
realizadas as demais observações. O importante é resultados da distribuição de erros e também obter
que o rigor do levantamento desta base seja um número de equações de observações não
aplicado a todo levantamento para minimizar a correlacionadas que possibilitem um certo grau de
propagação de variâncias na determinação dos liberdade para que o sistema seja possível. Por isso,
ângulos. uma estratégia é procurar aumentar a rigidez da
A trilateração é o método de rede através da inserção de pontos de controle.
levantamento topográfico baseado exclusivamente Estes pontos devem ter alta precisão e estar em um
na medição de distâncias horizontais (Melo et al, referencial compatível com aquele da rede ou da
2012). Esse método apresenta a vantagem de ser poligonal. Por isso, o levantamento deve ser
mais rapidamente executável que a triangulação e pensado para combinar as técnicas de campo com
produzir os resultados aceitáveis em precisão em a rigidez da rede.

[a] Trilateração [b] Triangulação

FIGURA1: Exemplos de [a] Trilateração e [b] Triangulação. Fonte: Adaptado de Melo et al (2012).
D. V. SEGUNDO, M. A. R. SCHMIDT - REEC – Revista Eletrônica de Engenharia Civil Vol 11 - nº 1 ( 2016)
70
2.1 AJUSTAMENTO PELO M.M.Q. - CASO Dentre os métodos de ajustamento
PARAMÉTRICO existentes no MMQ, esta pesquisa selecionou o
método paramétrico, no qual os parâmetros a
O ajustamento de observações é um
serem ajustados e os resíduos a serem estimados
ramo da matemática aplicada que, com a
são as coordenadas planas dos vértices em si e a
estatística, objetiva estimar valores únicos e mais
determinação das suas variâncias. No método
prováveis para um conjunto de observações
paramétrico, para cada observação de campo
(GEMAEL, 1994). Em toda a coleta de observações,
(ângulo e distância) formula-se ao menos uma
os valores mensurados estão sujeitos a erros
equação baseada no modelo matemático adotado
grosseiros, sistemáticos e acidentais. Gemael
para o ajustamento paramétrico. É importante que
(1994) afirma que os erros grosseiros e sistemáticos
essas equações sejam não correlacionadas para
são eventualmente detectados com relativa
que o sistema de resolução em (Equação 1) seja
facilidade, os erros sistemáticos podem ser
possível e determinável. Para isso, o número de
modelados por diferentes técnicas, porém os erros
graus de liberdade deve ser maior que 1, ou seja, o
acidentais ou aleatórios não podem ser modelados
número de equações de observações do método
dada a sua variabilidade. Devido às incertezas
paramétrico deve ser maior que o número de
desses erros, usa-se o ajustamento de observações
incógnitas. Considerando-se que se tenha n
para encontrar um valor estatisticamente
observações e u parâmetros, então obtém-se um
representativo do valor real da observação.
conjunto de n equações com u parâmetros com a
O ajustamento por M.M.Q. permite
seguinte representação matricial (Equação 2):
verificar a precisão e variância das coordenadas e
das observações como forma de identificar e
separar erros de observação de reais La=F(Xa) Eq. [2]
deslocamentos. Segundo Gemael (1994)
Em que:
ajustamento de observações diz respeito à
La: vetor (n x 1) de observações ajustadas;
minimização dos erros acidentais das medidas, com
Xa: vetor (u x 1) de parâmetros desconhecidos
a utilização do método dos mínimos quadrados, o
ajustados;
qual tem como fundamento minimizar a soma dos
F(Xa): função representativa do modelo
quadrados dos resíduos. O objetivo é resolver os
matemático que relaciona as observações
valores de parâmetros e observações através da
ajustadas com os parâmetros ajustados.
superabundância de observações, extraindo um
resultado único que representa com maior
A resolução é dada pela equação do
confiança a grandeza medida e os parâmetros
método paramétrico na sua forma explicita
estimados. As observações são relacionadas com os
pela Equação 3. A matriz (n x u) das derivadas
parâmetros por meio de modelos matemáticos
devidamente escolhidos. O princípio do Método parciais da função F é comumente chamada de
dos Mínimos Quadrados é que a soma dos matriz Jacobiana ou matriz A. Nela os modelos
quadrados dos resíduos das observações deve ser matemáticos adotados são derivados em
“mínimo”, como na Equação 1. relação a cada um dos parâmetros que se
deseja determinar para a criação da matriz A e
∑ 𝑉 𝑇 . 𝑃. 𝑉 = 𝑚𝑖𝑛 Eq. [1] os valores iniciais do levantamento são
determinados aplicando-se os valores
Em que:
encontrados em campo, chamado de valores
V: vetor dos resíduos das observações tomadas em
iniciais.
campo;
P: matriz dos pesos a priori; X = (ATPA)-1 . (ATPLb) = N-1.U Eq. [3]
D. V. SEGUNDO, M. A. R. SCHMIDT - REEC – Revista Eletrônica de Engenharia Civil Vol 11 - nº 1 ( 2016)
71

Em que: Melo et al (2012) resgatam que para

A: matriz jacobiana ou matriz das derivadas das redes topográficas o ajustamento pode ser
equações paramétrica em função dos parâmetros; realizado por diferentes programas como o Starnet
P: matriz dos pesos definidos a priori; da Starplus Software, o SURPAC Surveying
Lb: vetor das observações de campo; Software, entre outros. Existem as opções ainda de
N: matriz dos coeficientes dos parâmetros programas livres como o ADJUST e o WOLFPACK de
desconhecidos da equação normal; Wolf e Ghilani (1997, 2002, 2010), disponibilizados
U: matriz das estimativas dos parâmetros gratuitamente na internet.
desconhecidos a partir das observações.
2.2 ELIPSE DOS ERROS
A avaliação dos resíduos é dada pela
Equação 4: O uso da técnica do círculo de Mohr para
V= A. X + L Eq.[4] avaliar a precisão de ajustamento, baseia-se nos
Em que: valores da matriz variância-covariância dos
A: matriz jacobiana ou matriz das derivadas das parâmetros desconhecidos ajustados (∑xy)
equações paramétrica em função dos parâmetros; (Amorim, 2005). Segundo este autor, os valores de
X: vetor dos parâmetros estimados; variância e covariância dos pontos são
L: vetor das observações ajustadas. representados em posições diametralmente
opostas no círculo de Mohr, através de dois pontos
A partir dos resíduos é possível se avaliar
X e Y. As coordenadas desses pontos correspondem
a qualidade da trasnformação através da análise do
aos valores da variância (σ) e covariância (σxy).
grau de aderência entre a variância a priori, defiida
Como as variâncias em X e Y são diferentes, o
genéricamene como 1, e a variância a posteriori
círculo de Mohr se degenera em uma elipse. Para o
calculada a partir dos resíduos encontrados e do
caso desses valores serem iguais, a forma do círculo
número de graus de lliberdade (Equação 5). Esta
é mantida.
análise é realizada através do teste Qui-quadrado.
Na hipótese de distribuição normal dos
Eq.[5] erros, a interpretação geométrica da elipse em sua
forma padrão pode ser entendida como uma
Em que:
probabilidade da posição estimada para o ponto
: Sigma a posteriori; esteja dentro da elipse, centrada na posição
V: vetor dos resíduos das observações tomadas ajustada. Em outras palavras, a posição ajustada é
em campo; o valor mais provável e a precisão desta posição é
P: matriz dos pesos a priori;
dada pelo tamanho da elipse. Quanto menor a
n: número de equações de observação;
u: número de parâmetros incógnitas. elipse de erros mais preciso é o ponto, uma vez que
as variâncias medidas nos eixos desta elipse são
A matriz Variância-Covariância é: menores quando comparados aos demais vértices
determinada pela equação 6: da poligonal. Para obtenção uma região de
probabilidade de 95%, basta multiplicar o semieixo
∑ 𝑋𝑎 = 𝜎02 . 𝑁 −1 maior e menor, obtido em um processo de
Eq.[6]
ajustamento, por um fator de 2,447 (Carvalho e
Em quee: Leandro, 2007).
∑Xa: Matriz variância- covariância dos parâmetros A orientação da elipse é dada pela
ajustados;
variação do sinal da covariância σxy e pela análise da
: Sigma a posteriori; variância em X em relação à de Y. A Figura 2 mostra
N: matriz dos coeficientes dos parâmetros os elementos do círculo de Mohr para a análise da
desconhecidos da equação normal.
elipse de erros.
D. V. SEGUNDO, M. A. R. SCHMIDT - REEC – Revista Eletrônica de Engenharia Civil Vol 11 - nº 1 ( 2016)
72

FIGURA 2: Círculo de Mohr adaptado para avaliação das precisões do ajustamento. Fonte: Amorim, (2005).

3. MATERIAIS E MÉTODOS Uberlândia de forma a serem visíveis entre si e

apresentadas na figura 3. As coordenadas
Para simular uma situação real de foram levantadas pelos três métodos
implantação de uma grande obra, foi estabelecida discutidos anteriormente, caminhamento,
uma poligonal de 07 estações, distribuídas no trilateração e triangulação. Neste último caso
entorno de um novo bloco de salas de aula no de a base foi determinada no mesmo alinhamento
campus Santa Mônica da Universidade Federal de de azimute.

FIGURA 3: Localização da poligonal utilizada nos experimentos. Fonte: Autores.

D. V. SEGUNDO, M. A. R. SCHMIDT - REEC – Revista Eletrônica de Engenharia Civil Vol 11 - nº 1 ( 2016)
73
Para fins de avaliação metodológica, os Os elementos da matriz jacobiana A, quando não
pontos de azimute, CIV 24 e CIV25 da poligonal, nulos são obtidos por:
foram fixados com peso de 0,0001 m em X e em Y. a) Para ângulos horizontais (Amorim, 2005):
Portanto a matriz A foi derivada em relação aos 𝜕𝐴𝑧𝑖𝑗 (𝑦𝑗 −𝑦𝑖 ) 𝜕𝐴𝑧𝑖𝑗 (𝑦𝑗 −𝑦𝑖 )
demais pontos não fixados. Entretanto, para =− 2 ; = 2 ;
𝜕𝑥𝑖 𝑑𝑖𝑗 𝜕𝑥𝑗 𝑑𝑖𝑗
𝜕𝐴𝑧𝑖𝑗 (𝑥𝑗 −𝑥𝑖 ) 𝜕𝐴𝑧𝑖𝑗 (𝑥𝑗 −𝑥𝑖 )
Eq.[10]
avaliação da propagação sobre as distâncias foi = ; =− ;
2 2
𝜕𝑦𝑖 𝑑𝑖𝑗 𝜕𝑦𝑗 𝑑𝑖𝑗
realizado a avaliação com a fixação somente do
ponto CIV24. Com isso, a matriz jacobiana teve as b) Para distâncias (Melo et al, 2012):
linhas e colunas relativas ao ponto CIV25 𝜕𝑑𝑖𝑗 (𝑥𝑖 −𝑥𝑗 ) 𝜕𝑑𝑖𝑗 (𝑥𝑗 −𝑥𝑖 )
adicionadas. A matriz dos pesos a priori (P) foi = ; = ;
𝜕𝑥𝑖 𝑑𝑖𝑗 𝜕𝑥𝑗 𝑑𝑖𝑗
determinada como a matriz diagonal cujos valores Eq.[11]
𝜕𝑑𝑖𝑗 (𝑦𝑖 −𝑦𝑗 ) 𝜕𝑑𝑖𝑗 (𝑦𝑖 −𝑦𝑗 )
são o inverso das precisões de ângulo e distância. = ; = ;
𝜕𝑦𝑖 𝑑𝑖𝑗 𝜕𝑦𝑗 𝑑𝑖𝑗
Nesta etapa não foi realizada estudos de
Para a realização do ajustamento foi
propagação de variância a priori e, portanto, o So foi
utilizado o programa ADJUST©, programa livre e
considerado como o inverso das precisões nominais
disponibilizado na internet pela Universidade da
das observações dados pela estação total utilizada.
Pensilvânia, e que exige uma configuração de
Os valores dos pesos adotados para as observações
entrada como mostra a Quadro 1.
foram determinados em função da estação total
Para cada uma das alternativas de rede
utilizada Ruide 822A de precisão angular de 05” e
(caminhamento, trilateração e triangulação) foram
linear de 2mm+ 2ppm.
realizados diversos experimentos alterando as
O modelo matemático adotado para o
configurações básicas como número de vértices
ajustamento paramétrico foi o de ângulos
fixados, número de distâncias lidas, número de
horizontais e distâncias. Os ângulos horizontais em
ângulos lidos, e assim por diante. Foram testados as
uma poligonal podem ser entendidos como uma
poligonais por caminhamento com 2 e 3 ângulos
diferença entre o azimute de vante e o de ré,
adicionais internos e variações de 3 a 6 distâncias
considerando o sentido de desenvolvimento da
internas com fixação de 1 e 2 pontos. Na
poligonal (Equações 8 e 10) e as distâncias são
triangulação foram avaliadas variações de 2 a 4
diferenças em coordenadas euclidianas (Equações
ângulos adicionais e na trilateração foram avaliadas
9 e 11). A adoção do plano topográfico como
variações de 12 a 16 distâncias com fixação de 1 e 2
referência é uma simplificação que permite esta
pontos.
interpretação desde que as leituras sejam curtas. As
equações genéricas abaixo consideram que os
pontos i, j e k pertencem a uma mesma poligonal 4. RESULTADOS E DISCUSSÕES
sendo i o ponto ocupado pelo teodolito, k o ponto
Foi realizado o ajustamento das observações pelo
de ré e j o ponto de vante à i nesta poligonal.
método paramétrico da poligonal, onde foi adotado
𝐴𝑧𝑖𝑗 = 𝐻𝑗 − 𝐻𝑖
como critério de convergência para parada da
(𝑥𝑗 − 𝑥𝑖 ) (𝑥𝑖 − 𝑥𝑘 )
𝐴𝑧𝑖𝑗 = tan−1 − tan−1 Eq.[8] iteração, um valor de δ ≤ 0,01mm para as correções
(𝑦𝑗 − 𝑦𝑖 ) (𝑦𝑖 − 𝑦𝑘 )
aos parâmetros aproximados em relação à iteração
2 2
𝑑𝑖𝑗 = √(𝑥𝑗 − 𝑥𝑖 ) + (𝑦𝑗 − 𝑦𝑖 ) Eq.[9] anterior. Segundo Dalmolin (2002, pag. 96) este
Em que: valor pode ser arbitrário analisadas as condições do
Azij: azimute do alinhamento entre os pontos i e j, levantamento. No contexto desta pesquisa,
genéricos; variações abaixo de 1cm no vetor das coordenadas
X, Y: coordenadas planas cartesianas dos pontos i, ajustadas tornam-se não relevantes. Em geral
j, k, genéricos;
foram necessárias dez iterações em cada uma das
H: ângulos horizontais entre as direções do
alinhamento entre pontos i e j, i e k; poligonais para que se atingisse a especificação
d: distância euclidiana entre os pontos i e j. de δ.
D. V. SEGUNDO, M. A. R. SCHMIDT - REEC – Revista Eletrônica de Engenharia Civil Vol 11 - nº 1 ( 2016)
74

Quadro 1: Organização do arquivo no programa Adjust©.

Exemplo de preenchimento Legenda do dado

Example 13-3 of Elementary Surveying: An Linha do título, ou seja, identificação do trabalho
Introduction in Geomatics, 10th Ed. text book
# de distâncias # de ângulos horizontais lidos # de azimutes
55115 conhecidos # de vértices fixados # total de vértices da poligonal
ou rede
A 10000.00 10000.00 Estações com coordenadas X e Y fixadas (vértices de controle)
B 10125.66 10255.96 Lista das coordenadas X e Y aproximadas dos vértices calculadas
C 10716.31 10102.44 pelo método convencional
D 10523.62 9408.37
E 10517.55 9611.34
A B 285.10 0.05 Lista das distâncias aproximadas entre os vértices lidos em campo
B C 610.45 0.05 e respectivos desvios-padrão (em m)
C D 720.48 0.05
D E 203.00 0.05
E A 647.02 0.05
B A E 100 44 30 15 Lista de ângulos horizontais lidos em campo e respectivos desvios-
C B A 101 35 00 15 padrão (em segundos)
D C B 89 05 30 15
E D C 17 12 00 15
A E D 231 24 30 15
A B 26 10 00 1 Azimute observado e desvio-padrão (em segundos)

Os valores obtidos pelo método resíduo. Isto pode ser devido a ser o ponto mais
tradicional tiveram as variâncias calculadas distante do ponto CIV 24, fixado.
baseadas nas séries obtidas em campo, sendo que Para a triangulação variou-se o número
a maior parte dos pontos teve 4 séries lidas e alguns de ângulos lidos entre os vértices de 12 a 16 leituras
pontos somente 3 séries. Mesmo que a média seja com a fixação de 2 vértices (uma base) a 3 vértices
um valor mais provável para uma observação (duas bases adjacentes). Este método
(Gemael, 2004) os seus resíduos parecem estar apresentou os piores resultados, sendo que exigiu
superestimados. Quando comparados aos um grande número de observações angulares.
resultados obtidos por outros métodos. Em outras palavras, para a solução do sistema
Para a trilateração foram variadas o foi necessário criar diversas equações de
número de distâncias internas lidas entre 11 e 16 e observações internas para ser possível a
o foram fixados somente o vértice CIV24 e o CIV24 distribuição por MMQ. O problema, desta
e CIV25. A variação da distribuição dos resíduos técnica está em incorrer em equações
teve maior impacto relacionado à fixação do ponto correlacionadas. A leitura de ângulos entre os
CIV25 do que ao número de distâncias medidas. Em vértices da poligonal pode gerar equações
outras palavras os valores das coordenadas finais correlacionadas, o que levaria o determinante do
sofreram maior variação quando comparadas à sistema paramétrico (matriz A) tender a zero. Isto
técnica combinada. Além disso, os valores do teste ocorreu com os ângulos medidos no vértice CIV33
de qui-quadrado que aceitaram o resultado do que geraram resíduos acima de 3’. Porém ao se
ajustamento foram com 11 e 12 medidas lineares, retirar uma das observações deste ponto, os
sendo 7 da poligonal de contorno e as demais resíduos diminuíram. Mesmo assim, os valores
tomadas internamente. Em ambos os casos, o obtidos por este método estão muito acima do
vértice CIV33 apresentou os mais altos valores de desejável, como mostrado na Quadro 2.
D. V. SEGUNDO, M. A. R. SCHMIDT - REEC – Revista Eletrônica de Engenharia Civil Vol 11 - nº 1 ( 2016)
75
Na análise combinada os dados da sistema, pois no teste de qui-quadrado o
poligonal sofreram variações de distâncias de 7, valor ficou em 0,07 com limites entre 0,0459 e
número de alinhamentos da poligonal, a 13, os 0,502.
alinhamentos da poligonal acrescidos das diagonais Como as variâncias e covariâncias são
internas; variação no número de pontos fixados: utilizadas nos cálculos das elipses, percebe-se que
somente CIV 24 ou fixados os pontos CIV24 e CIV25; os dados angulares estão eivados de erros, com
e variações no número de ângulos internos lidos de destaque para o ponto CIV33. Este ponto foi o que
7 a 11, sendo estes os números entre os apresentou maior dificuldade para sua
alinhamentos ou os alinhamentos e as diagonais determinação em campo, por ter sofrido com
inseridas posteriormente. vandalismo e a construção de um painel de
A qualidade do ajustamento foi analisada propagando nas suas proximidades. Percebe-se
com base no teste estatístico qui-quadrado a um além disso, que nas trilaterações os maiores
nível α= 5% de significância. Os valores encontrados resíduos ocorrem também neste ponto. Os erros
para as variâncias e covariâncias entre as neste caso não se limitam ao ponto, mas foram
coordenadas estimadas por cada um dos métodos distribuídos entre os demais vértices não fixados
foram utilizados para o desenho das elipses dos conforme o desenvolvimento dos MMQ.
erros. A figura 3 apresenta a formulação da rede e É interessante notar que na análise
das elipses de erros computadas no ajustamento. A combinada, as variações no número de ângulos não
triangulação por ter somente uma base não contribuíram com a redução dos resíduos. A fixação
permite que o programa gere a estimativa da rede do segundo ponto (CIV25) gerou um ligeiro
antes do ajustamento. aumento no resíduo em X, com o consequente
Na Figura 4a, percebe-se que a deslocamento das coordenadas dos pontos nesse
trilateração apresenta uma variância das distâncias eixo. Isso pode indicar que exista um erro
no sentido esquerda-direita de maneira sistemático que deixa de ser modelado neste
semelhante. Isto pode sugerir a presença de erros vértice. Os resultados obtidos indicam os menores
sistemáticos ao longo do eixo X. Na Figura 4b, o resíduos do que quando os métodos são aplicados
método combinado redistribui essas tendências de isolados. Além disso, dos métodos aplicados em
variâncias de modo que as elipses se aproximam ao relação à coordenada calculada pelo método de
desenho de um círculo. Bowditch NBR13.133 (ABNT, 1994) a combinação
A combinação de três diagonais resultou de ângulos e distâncias apresentou o menor
nos valores com menor resíduo (Quadro 2). Esse deslocamento de coordenadas, em especial
número de distâncias adicionais possibilitou um daqueles vértices mais distantes dos pontos
maior número de equações sem gerar tendência no fixados.

a) b)
FIGURA 4: Elipses de erros propagados na rede (a) com trilateração, e (b) rede combinada de ângulos e distâncias.
FONTE: Autores a partir do programa Adjust.
D. V. SEGUNDO, M. A. R. SCHMIDT - REEC – Revista Eletrônica de Engenharia Civil Vol 11 - nº 1 ( 2016)
76

Quadro 2: Resultados obtidos com a variação dos métodos combinado, trilateração e triangulação.

COORDENADAS OBTIDAS

Método tradicional Triangulação (base com 2

Combinada (1 ponto fixo) Trilateração (1 ponto fixo) Trilateração (2 pontos fixos)
NBR13.133/1994 pontos fixos)
Nome do X Y X Y X Y X Y X Y
ponto (m) (m) (m) (m) (m) (m) (m) (m) (m) (m)
CIV24 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
CIV25 903,831 835,665 903,829 835,661 903,831 835,665 903,831 835,665 904,77 837,268
CIV26 1078,13 892,924 1078,127 892,922 1078,05 893,167 1078,073 893,097 1077,406 893,921
CIV32 1185,327 808,769 1185,327 808,767 1184,372 807,909 1185,11 808,739 1185,22 809,226
CIV33 1402,327 791,186 1402,328 791,187 1401,281 789,23 1402,099 791,013 1401,137 791,676
CIV20 1435,74 948,033 1435,74 948,029 1435,523 945,897 1435,534 947,856 1435,762 947,798
CIV27 1328,778 1072,051 1328,781 1072,047 1329,131 1070,405 1328,778 1070,051 1328,778 1072,051

RESÍDUOS CALCULADOS

Nome do σ²x σ²y σ²x σ²y σ²x σ²y σ²x σ²y σ²x σ²y
ponto (m) (m) (m) (m) (m) (m) (m) (m) (m) (m)
CIV24 * --- --- --- --- --- --- --- --- --- ---
CIV25 0,0007 0,0007 0,0034 0,0057 0,000 0,000 0,0010 0,0016 28,7228 49,0822
CIV26 0,0043 0,0043 0,0064 0,0071 0,052 0,055 0,0738 0,2237 22,9910 23,5449
CIV32 0,0033 0,0033 0,0080 0,0091 0,107 0,128 0,1892 0,2198 29,4930 16,9248
CIV33 0,0024 0,0024 0,0096 0,0137 0,118 0,243 0,1933 0,1393 86,5468 39,8534
CIV20 0,0047 0,0047 0,0083 0,0141 0,036 0,266 0,1490 0,1281 21,3272 25,1320
CIV27 ** 0,0056 0,0056 0,0089 0,0104 0,047 0,221 --- --- --- ---
* Ponto Civ24 fixado com variância de 0,001m respectivamente a σ²x e σ²y
** Ponto Civ27 fixado somente nos dois últimos casos da tabela com precisão igual a CIV24.
D. V. SEGUNDO, M. A. R. SCHMIDT - REEC – Revista Eletrônica de Engenharia Civil Vol 11 - nº 1 ( 2016)
77

As Figuras 4, 5 e 6, apresentam pontos mais distantes do CIV24, início da poligonal.

as variações dos resíduos das coordenadas dos A Figura 6, apresenta as variações dos
vértices por método. Nos gráficos são resíduos quando consideradas somente as
apresentadas as variações em X e em Y de distâncias observadas. No caso da triangulação
acordo com a inclusão de novas observações de com um ponto fixado, o número mínimo de
campo na matriz A. Na Figura 5, o acréscimo de observações foi 12 (7 da poligonal e mais
observações de distâncias adicionais internas à 5 observações internas) para o sistema ser
poligonal manteve as coordenadas estáveis. possível de resolução. Percebe-se que com 13
Entretanto, o acréscimo de ângulos internos distâncias os resíduos sofrem uma redução
aumentou os deslocamentos com destaque aos importante e depois se estabilizam.

FIGURA 5: Variação de coordenadas X e Y de cada vértice pelo método combinado.

Fonte: Autores.

FIGURA 6: Variação de coordenadas X e Y de cada vértice pelo método de trilateração com 1 vértice fixo (CIIV24).
Fonte: Autores.
D. V. SEGUNDO, M. A. R. SCHMIDT - REEC – Revista Eletrônica de Engenharia Civil Vol 11 - nº 1 ( 2016)
78
Na Figura 7, estão apresentadas as Verificou-se que o MMQ permite a
variações dos resíduos considerando a fixação de identificação da existência de erros nas
um ponto adicional, o CIV27 ré do ponto CIV24. observações de campo através da verificação das
O experimento utilizou os mesmos dados da variações das coordenadas e da determinação dos
trilateração com único ponto, porém percebe-se resíduos dessas coordenadas em cada método.
que a dimensão das variações é semelhante a O método que combina ângulos e
anterior com exceção do primeiro teste com distâncias atingiu os melhores resultados e permite
12 distâncias. Em ambas as trilaterações, o afirmar que as coordenadas apresentadas no
valor mínimo dos resíduos foi obtido com 13 Quadro 2 apresentam confiabilidade superior aos
distâncias. Este resultado pode indicar que demais, uma vez que os resíduos ficam abaixo de
algumas das observações adicionais pode 1,5 cm. Os outros métodos, envolvendo somente
conter algum tipo de inconsistência sistemática ou distâncias ou somente ângulos, apresentaram
acidental que voltou a elevar os valores. maiores variações de coordenadas e variâncias. Isto
A triangulação não apresentou resultados permite afirmar que o referencial para locação de
melhores que o apresentado no Quadro 2. obras na engenharia civil deve ser prioritariamente
Além disso, a combinação de ângulos internos definido a partir de uma poligonal de controle que
degradou a qualidade do levantamento permita visadas internas, além da condição de
indicando erros nas observações de campo intervisibilidade dos pontos ré e vante. Dessa
que se propagaram na determinação dos forma, garante-se um número maior de equações
vértices. do que incógnitas e o MMQ pode ser aplicado ao
cálculo de poligonais. Além disso, deve-se
5. CONCLUSÕES aumentar o número de pontos fixados ao longo das
maiores diagonais como forma de minimizar os
Esta pesquisa avaliou a propagação de valores de propagação.
erros na estimativa das coordenadas de uma Como continuidade dos trabalhos sugere-
poligonal de controle para locação de obras de se a variação da forma e dimensões da poligonal
engenharia e das precisões destas coordenadas em bem como avaliar os efeitos da propagação de
função dos diferentes métodos existentes. Foram variâncias a partir de equipamentos diversos e com
avaliadas diversas variações no número de maior tempo de uso, na qualidade dos dados
observações adicionais de distâncias e ângulos. obtidos.

FIGURA 7: Variação de coordenadas X e Y de cada vértice pelo método de trilateração com 2 vértices fixos (CIV24 e
CIV27). Fonte: Autores.
D. V. SEGUNDO, M. A. R. SCHMIDT - REEC – Revista Eletrônica de Engenharia Civil Vol 11 - nº 1 ( 2016)
79

6. AGRADECIMENTOS

Parte desta pesquisa tem apoio do CNPq

através do Edital Universal nº 04/2012.

7. REFERÊNCIAS BIBLIOGRÁFICAS

AMORIM, G. P. 2005. Curso de formação continuada em

Georreferenciamento Aplicado ao Cadastro Rural -
Ajustamento de Observações. MEC/SEMT/CEFET – ES.

ASSOCIAÇÃO BRASILEIRA DE NORMAS TÉCICAS.

NBR13.133 (1994). Execução de levantamento
topográfico. Rio de Janeiro. Brasil. 1994

________. ADJUST 6.0.2. 2013. Disponível em:

http://www.personal.psu.edu/cdg3/free.htm Acessado
em: 15/11/2014

CARVALHO, A. S. e LEANDRO. D. 2007. Verificação da

precisão das coordenadas ajustadas em poligonal
fechada e enquadrada. Anais: XXIII Congresso Brasileiro
de Cartografia, Rio de Janeiro, Brasil.

COMASTRI, J, 1994. Topografia: Planimetria, Imprensa

Universitária, Viçosa - MG, Brasil, 336 páginas.

GEMAEL, C. Introdução ao ajustamento de observações:

aplicações geodésicas. Curitiba: Editora da UFPR, 1994.
319p

GHILANI, C. D. e WOLF, P. R. 2007. Adjustment

Computations: Spatial Data Analysis, 10a Edição. ISBN:
9780471697282

GHILANI, C. D. e WOLF, P.R. 2010. Geomática, 13ª

EDIÇÃO. ISBN: 9788581434506

DALMOLIN, Quintino, 2002. Ajustamento por Mínimos

Quadrados, Imprensa Universitária, Curitiba, Brasil, 175
páginas.

MELO, W. D. A.; SOUZA, A. N. e DA SILVA, D. C. 2012.

Utilização do programa adjust em ajustamento de
triangulações e trilaterações. Anais: IV Simpósio
Brasileiro de Ciências Geodésicas e Tecnologias da
Geoinformação Recife – PE.