VTKファイルをつくる7

 投稿者：しばっち  投稿日：2018年 3月10日(土)21時01分52秒
  座標を計算して、その点と点を線で結んでいます。
ファイルを読み込んでも、ただの線しか表示されません。
まず、FiltersメニューのAlphabeticalから「Tube」を選択すると、
線をチューブ化することができます。
Radiusでチューブの太さを調整できます。

OPTION BASE 0
OPTION ANGLE DEGREES
LET N=720
DIM X(N),Y(N),Z(N)
RANDOMIZE
FILE GETSAVENAME F$,"vtkファイル|*.vtk"
IF F$="" THEN STOP
IF POS(UCASE$(F$),".VTK")=0 THEN LET F$=F$&".vtk"
OPEN #1:NAME F$
ERASE #1
PRINT #1:"# vtk DataFile Version 4.0"
PRINT #1:"lissajous"
PRINT #1:"ASCII"
PRINT #1:"DATASET POLYDATA"
PRINT #1:"POINTS";N;"double"
LET A=INT(RND*20)+1
LET B=INT(RND*20)+1
LET C=INT(RND*20)+1
LET D=INT(RND*360)
FOR I=0 TO N-1
   LET X(I)=COS(A*I/2) !'リサジュー曲線
   LET Y(I)=SIN(B*I/2)
   LET Z(I)=COS(C*I/2+D)
   PRINT #1:X(I);Y(I);Z(I) !'座標値
NEXT I
PRINT #1:"LINES";N;3*N
FOR I=0 TO N-1   !'線でつないでいく(meshデータ)
   IF I=N-1 THEN
      PRINT #1:2;N-1;0
   ELSE
      PRINT #1:2;I;I+1
   END IF
NEXT I
PRINT #1:"POINT_DATA";N
PRINT #1:"SCALARS distance float"
PRINT #1:"LOOKUP_TABLE default"
FOR I=0 TO N-1
   PRINT #1:SQR(X(I)^2+Y(I)^2+Z(I)^2)
NEXT I
CLOSE #1
END
 

戻る