PureBoard

Code: Alles auswählen

;- Phasenraumplot
;- "Helle" Klaus Helbing, 07.04.2008, PB4.10 
;- Dieser Phasenraumplot stellt die Bytes z.B. einer Datei 3-dimensional als Raumpunkte dar
;- Die Raum-Koordinaten errechnen sich folgendermaßen:
;- Punkt1=(Byte1-Byte2, Byte2-Byte3, Byte3-Byte4), Punkt2=(Byte2-Byte3, Byte3-Byte4, Byte4-Byte5) usw.
;- Gut verschlüsselte Daten und gute Zufallszahlen ergeben im Phasenraumplot eine gleichmäßig
;- verteilte Punktewolke, ansonsten sieht man Strukturen (z.B. Linien) 
;- Z wurde aus Platzgründen um Faktor 2 gestaucht
;- Zur besseren Übersicht wurde der Koordinaten-Ursprung nach oben verschoben

Cos10.f = Cos(10*#PI/180)
Sin10.f = Sin(10*#PI/180)
Cos60.f = 0.5                                    ;Cos(60°)
Sin60.f = Sin(60*#PI/180)

AX.l = 50
AY.l = 548
BX.l = AX + 0.75*510*Sin60                       ;0.75=Verkürzungsfaktor 
BY.l = AY - 0.75*510*Cos60
CX.l = BX + 510*Cos10
CY.l = BY + 510*Sin10
DX.l = AX + 510*Cos10
DY.l = AY + 510*Sin10

Buffer = AllocateMemory(100000)                  ;100000 Bytes für die zu testende Datei

If OpenWindow(0, 0, 0, 900, 675, "Phasenraumplot", #PB_Window_SystemMenu | #PB_Window_ScreenCentered) And CreateGadgetList(WindowID(0))
  TextGadget(0,200,10,500,20,"Zu testende Datei auswählen :",#PB_Text_Center)  
  ExplorerTreeGadget(1, 10, 30, 880, 620,"",#PB_Explorer_NoDriveRequester)

  Repeat
    Event=WaitWindowEvent()
    If Event=#PB_Event_CloseWindow
      End 
    EndIf   
  Until EventType()=#PB_EventType_LeftDoubleClick And GetGadgetState(1)=#PB_Explorer_File
 
  File$=GetGadgetText(1)

  ReadFile(0,File$)                              ;nur zum Lesen öffnen
  FileL=Lof(0)
  
  If FileL > 100000                              
    FileL = 100000                               ;hier auf 100000 Bytes begrenzt
  EndIf  

  ReadData(0, Buffer, FileL)

  If StartDrawing(WindowOutput(0))
    Box(0, 0, 900, 675, RGB(255,255,255))
    ;Z-Bemaßung
    LineXY(AX, AY, AX, 170)
    LineXY(45, 170, 55, 170)                     ;Querstrich +255
    LineXY(45, 298, 55, 298)                     ;Querstrich 0
    LineXY(45, 425, 55, 425)                     ;Querstrich -255
    DrawText(35 ,150, "( Z )")                   ;ist nicht die Achse, sondern die Bemaßung!
    DrawText(5, 165, "+255")
    DrawText(25, 290, "0")
    DrawText(5, 420, "-255")
    DrawText(475, 15, "Z", RGB(0, 0, 255), RGB(255, 255, 255))
    ;Y-Bemaßung
    LineXY(AX, AY, DX, DY)                       ;unten
    LineXY(BX, BY, CX, CY)                       ;oben
    YMXU = AX + (DX-AX)/2
    YMYU = AY + (DY-AY)/2
    YMXO = BX + (CX-BX)/2
    YMYO = BY + (CY-BY)/2
    LineXY(YMXU, YMYU, YMXO, YMYO, RGB(0, 0, 255))    ;blau
    DrawText(DX-40, DY, "( Y )")                 ;ist nicht die Achse, sondern die Bemaßung!
    DrawText(DX-10, DY+10, "+255")
    DrawText(YMXU-8, YMYU+3, "0")
    DrawText(AX-10, AY+5, "-255")
    DrawText(670, 535, "Y", RGB(0, 0, 255), RGB(255, 255, 255))
    ;X-Bemaßung
    LineXY(AX, AY, BX, BY)                       ;links
    LineXY(DX, DY, CX, CY)                       ;rechts
    XMXL = AX + (BX-AX)/2
    XMYL = BY + (AY-BY)/2
    XMXR = DX + (CX-DX)/2 
    XMYR = CY + (DY-CY)/2
    LineXY(XMXL, XMYL, XMXR, XMYR, RGB(0, 0, 255))    ;blau 
    DrawText(DX+10, DY-5, "( X )")               ;ist nicht die Achse, sondern die Bemaßung!
    DrawText(XMXR+5, XMYR, "0")
    DrawText(CX-10, CY-20, "-255")
    DrawText(335, 575, "X", RGB(0, 0, 255), RGB(255, 255, 255))
    ;"echte" Z-Achse  
    MXZ = YMXU + (YMXO-YMXU)/2  
    MYZ = YMYO + (YMYU-YMYO)/2   
   ;Nullpunkt = (MXZ, MYZ-AY+170+128)  
    MX = MXZ
    MY = MYZ-AY+170+128

;--------------------------------------------------------------------------------------------------
    For R=0 To FileL - 4
      T = PeekL(Buffer + R)

      A.c = T & 255
      B.c = T >> 8
      T>>8
      C.c = T >> 8
      T>>8
      D.c = T >> 8

      PX = A-B
      PY = B-C
      PZ = C-D

      PXX = MX - 0.75*PX*Sin60
      PXY = MY + 0.75*PX*Cos60
      
      PYX = PXX + PY*Cos10 
      PYY = PXY + PY*Sin10

      PZY = PYY - PZ/2
     
      Plot(PYX, PZY, RGB(255, 0, 0))             ;roten Punkt zeichnen
    Next 

    ;"echte" Z-Achse mit Koordinaten-Ursprung zur Orientierung nachträglich einzeichnen
    LineXY(MXZ, MYZ, MXZ, MYZ-AY+70, RGB(0, 0, 255))  ;blau  
    LineXY(MXZ-5, MYZ-AY+170+128,MXZ+5, MYZ-AY+170+128)
    DrawText(MXZ+5, MYZ-AY+170+128-8, "(0,0,0)")  

    StopDrawing()
  EndIf
  Repeat : Event = WaitWindowEvent() : Until  Event = #PB_Event_CloseWindow
EndIf