向量计算两条link的夹角

                            double x1 = (*geom1)[0][num1-1].x();

                            double y1 = (*geom1)[0][num1-1].y();

                            double x2 = (*geom1)[0][num1-2].x();

                            double y2 = (*geom1)[0][num1-2].y();

                            double x3 = (*geom2)[0][0].x();

                            double y3 = (*geom2)[0][0].y();

                            double x4 = (*geom2)[0][1].x();

                            double y4 = (*geom2)[0][1].y(); 

                            double dx1 = x2 -x1;

                            double dx2 = x4 -x3;

                            double dy1 = y2 -y1;

                            double dy2 = y4 -y3;

                            double angle = (dx1*dx2+dy1*dy2)/sqrt((dx1*dx1+dy1*dy1)*(dx2*dx2+dy2*dy2));