Преобразование вершин LWPOLYLINE, записанных в формате DXF

[Андрей Бушман]

Доброго времени суток.

Как пересчитать в глобальную систему координат координаты вершин полилинии, которая сохранена в DXF файле как LWPOLYLINE? Из DXF файла можно прочитать координаты вершин и некоторые три числа, по которым и вычисляются глобальные координаты. Как из этих трёх чисел получить матрицу направляющих косинусов, дабы можно было выполнить обозначенное преобразование?

P.S. Преобразование будет выполняться во внешнем приложении, а не в плагине AutoCAD.

Спасибо.

[Александр Ривилис]

Группа 210 у примитива LWPOLYLINE содержит вектор нормали примитива. По этому вектору используя Arbitrary Axis Algorithm можно получить матрицу преобразования: http://www.autodesk.com/techpubs/autocad/acad2000/dxf/arbitrary_axis_algorithm_dxf_ab.htm
Следует еще учесть возможность того, что у этой полилинии есть уровень (т.е.  группа 38 отлична от 0) - в этом случае потребуется дополнительный сдвиг.

[Андрей Бушман]

Для нас это чревато, т.к. мы мосты проектируем. Хотелось бы ссылку и для второго варианта.

[Александр Ривилис]

Объясняю на пальцах:
1) Получаем матрицу преобразования mat ОСК->МСК (алгоритм по ссылке, которую я дал)
2) Создаем вектор v = [0 0 Z], где Z - уровень полилинии (группа 38)
3) Преобразуем вектор v по матрице преобразования mat и получаем вектор v_wcs
4) Преобразуем точку p_ocs (группа 10) полилинии по матрице mat и получаем точку p_wcs
5) Складываем точку p_wcs и вектор v_wcs - это и будет искомая точка.
Как-то так.

[bavilla]

Уважаемые коллеги! Можете на моем примере показать, как перевести координаты точки из OCS в WCS? Без использования функций Autocad, простой арифметикой. Мне надо считать узлы полилиний из DXF в стороннем приложении, а в справке по "Arbitrary Axis Algorithm" я мало что понял... Заранее спасибо
Например, для узла LWPOLYLINE:
...
38
500.0
10
-100.0
20
430.0
...
210
-0.216
220
0.976
230
0.0

[Александр Ривилис]

А так: https://forums.autodesk.com/t5/autocad-r12-r13-r14-archive/dxf-arbitrary-axis-algorithm-in-vb6/td-p/35267 ?

[bavilla]

А так: https://forums.autodesk.com/t5/autocad-r12-r13-r14-archive/dxf-arbitrary-axis-algorithm-in-vb6/td-p/35267 ?

Спасибо, именно так я и пробовал посчитать, но где-то ошибся. Потому и прошу, если есть время и возможность, показать на конкретной точке из моего примера.

[Александр Ривилис]

Спасибо, именно так я и пробовал посчитать, но где-то ошибся. Потому и прошу, если есть время и возможность, показать на конкретной точке из моего примера.

Времени как обычно нет. Так что показывай свой код. Посмотрим в чем в нём ошибка.

[bavilla]

вот пример на С, попробовал сделать как описано в http://www.autodesk.com/techpubs/autocad/acad2000/dxf/arbitrary_axis_algorithm_dxf_ab.htm#XREF_24259_DXF_AB

Код - C++ [Выбрать]

1. struct vector
2. {
3.   double x;
4.   double y;
5.   double z;
6. };
7.  
8. void cross_product(vector v1,vector v2,vector &result)
9. {
10.   result.x = v1.y*v2.z - v1.z*v2.y;
11.   result.y = v1.z*v2.x - v1.x*v2.z;
12.   result.z = v1.x*v2.y - v1.y*v2.x;
13. }
14.  
15. void calcOCS2WCS()
16. {
17.   vector N;  //Let the given normal vector be called N
18.   vector Wy;  //Let the world Y axis be called Wy, which is always (0,1,0)
19.   vector Wz;  //Let the world Z axis be called Wz, which is always (0,0,1)
20.  
21.   vector Ax,Ay,Az;
22.  
23.   N.x = -0.1525427977445352; //код 210
24.   N.y = 0.6904568740017512;  //код 220
25.   N.z = 0.7071067811865475;  //код 230
26.  
27.   Wy.x = 0.;
28.   Wy.y = 1.;
29.   Wy.z = 0.;
30.  
31.   Wz.x = 0.;
32.   Wz.y = 0.;
33.   Wz.z = 1.;
34.  
35.   if((fabs(N.x)<1/64.)&&(fabs(N.y)<1/64.))
36.   {
37.     cross_product(Wy,N,Ax);
38.   }
39.   else
40.   {
41.     cross_product(Wz,N,Ax);
42.   }
43.  
44.   cross_product(N,Ax,Ay);
45.  
46.   Az = N; //N could also be called Az
47.  
48.   //Исходная точка в OCS
49.   double OCS_X = 100.0;
50.   double OCS_Y = 120.0;
51.   double OCS_Z = 50.0;
52.  
53.   //То что нужно вычислить
54.   double WCS_X;
55.   double WCS_Y;
56.   double WCS_Z;
57.  
58.   //Предполагаю, что Ax,Ay,Az представляют матрицу коэффициентов афинных преобразований
59.   WCS_X = OCS_X*Ax.x + OCS_Y*Ax.y + OCS_Z*Ax.z;
60.   WCS_Y = OCS_X*Ay.x + OCS_Y*Ay.y + OCS_Z*Ay.z;
61.   WCS_Z = OCS_X*Az.x + OCS_Y*Az.y + OCS_Z*Az.z;
62.  
63. }

[Александр Ривилис]

//Исходная точка в OCS
  double OCS_X = 100.0;
  double OCS_Y = 120.0;
  double OCS_Z = 50.0;

Она у тебя вроде бы была совсем другая:

10
-100.0
20
430.0

Не учтена группа 38. Я выше писал, что с ней нужно делать.

[bavilla]

//Исходная точка в OCS
  double OCS_X = 100.0;
  double OCS_Y = 120.0;
  double OCS_Z = 50.0;

Она у тебя вроде бы была совсем другая:

10
-100.0
20
430.0

Не учтена группа 38. Я выше писал, что с ней нужно делать.

Я взял произвольную точку для примера, пусть группа 38 будет 0.0, тогда ее можно и не учитывать

[bavilla]

Правильно ли я считаю Ax, Ay, Az? И как правильно их использовать для последующего пересчета координат? Я ради интереса нашел для своего примера коэффициенты афинного преобразования с помощью метода наименьших квадратов, они совсем другие чем Ax,Ay,Az

[Александр Ривилис]

Я не увидел, чтобы ты у себя в коде нормировал Ax (и Ay и Az). Ты должен разделить их Ax.x, Ax.y, Ax.z на sqrt(Ax.x*Ax.x+Ax.y*Ax.y+Ax.z*Ax.z) (аналогично для Ay и Az). В результате получаешь матрицу преобразования:

Код: [Выделить]

| Ax.x  Ay.x  Az.x |
| Ax.y  Ay.y  Az.y |
| Ax.z  Ay.z  Az.z |Произведение вектора координаты в OCS на эту матрицу даст вектор координаты в WCS

[bavilla]

Спасибо! Правда, уже сам успел разобраться.

[Александр Ривилис]

Спасибо! Правда, уже сам успел разобраться.

Получилось? Значения совпали с реальными в WCS?