Calculando a área em uma polilinha?

Question

Olá colegas tikzers.

Existe um pacote/hack para calcular a área delimitada por uma polilinha, ou seja, por exemplo, algo assim

\draw (p1) to (p2) to (p3) to cycle;

onde os pontos são os vértices. O polígono pode não ser convexo.

Imagino que o algoritmo de preenchimento poderia ser melhorado/modificado para retornar o valor da área (nas unidades utilizadas). Mas não tenho conhecimento técnico e talvez já exista uma resposta.

Obrigado.

PS: editado para acrescentar que a área é simplesmente a integral ao longo do limite

$$ \oint x dy $$

então é muito fácil de calcular. É claro que não estou interessado apenas em quadriláteros, mas em polilinhas com muitos vértices, de modo que o cálculo manual se torna irritantemente longo.

Answer 1

Aplicação da fórmula Shoelace usando Tikbiblioteca z math. Como exemplo, é tomada a aplicação complexa no link especificado por @Thruston.

\documentclass[margin=3mm]{standalone}

\usepackage{tikz}  
\usetikzlibrary{math}
\usepackage{tkz-euclide}
\tikzmath{
\x1 = 3; \y1 =4; 
\x2 = 5; \y2 =11; 
\x3 = 12; \y3 =8;
\x4 = 9; \y4 =5;
\x5 = 5; \y5 =6;
\Det = (\x1*\y2) + (\x2*\y3) + (\x3*\y4) + (\x4*\y5) + (\x5*\y1)-
(\x2*\y1) - (\x3*\y2) - (\x4*\y3) - (\x5*\y4) - (\x1*\y5);
\Area = abs (\Det / 2);
} 

\begin{document}

\begin{tikzpicture}
\tkzInit[xmax=13,ymax=12,xmin=0,ymin=0]
   \tkzGrid
   \tkzAxeXY

\draw[fill=gray!30](\x1,\y1)--(\x2,\y2)--(\x3,\y3)--(\x4,\y4)--(\x5,\y5)--cycle;

 \tkzText [below](6.5,-1){$A = \Area$} ;
\end{tikzpicture}

\end{document}

Answer 2

Isso não deve ser muito difícil de fazer na área do usuário - enquanto você espera pela equipe do TikZ, aqui está um simplesAlgoritmo de cadarçofeito de forma simplesMetapost. Compilar mpostou adaptar lualatexcom luamplib.

prologues := 3;
outputtemplate := "%j%c.%{outputformat}";

vardef shoelace_area(expr p) = 
    for i=1 upto length p: 
        + 1/2 xpart point i of p * ypart point i+1 of p 
        - 1/2 ypart point i of p * xpart point i+1 of p 
    endfor
enddef;

beginfig(1);

    path t; t = for i=0 upto 4: 100 dir 72i -- endfor cycle;
    draw origin -- point 1 of t; label.lft(btex $R=100$ etex, 1/2 point 1 of t);
    draw t; dotlabel.rt(decimal shoelace_area(t), origin) withcolor 2/3 red;

    draw thelabel(btex $\displaystyle A = {5\over4}R^2\sqrt{5+\sqrt5\over2}$ etex, origin) scaled 3/4 shifted 30 down;
    draw thelabel(btex $A = 2.37764129 R^2$ etex, origin) scaled 3/4 shifted 54 down;

endfig;
end.

A figura de exemplo aplica-o a um pentágono:

Como tentei mostrar, isso dá cerca de 6 sig. Figo. de precisão usando a aritmética em escala simples do MP, e se você aplicá-la a algo muito maior do que isso, obterá um estouro aritmético. Mas se você compilar, mpost -numbersystem doubleobterá uma resposta precisa e menos chance de estouro: