Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  593 1732 (N41E42 L=296)   386     54        0 ( 565/1709- 623/1751) -18.4 E46-E37 N39-N42
      2  491 1707 (N39E52 L=287)   585     66        0 ( 471/1667- 519/1737) -80.5 E54-E49 N35-N41
      3 1562 1564 (N34W38 L= 16)   312     51        0 (1546/1532-1576/1592) -90.0 W36-W39 N31-N35
      4  337 1547 (N27E58 L=281)   597     49        1 ( 316/1520- 359/1572)  90.0 E61-E54 N24-N28
      5  327  846 (S17E48 L=291)   486    130        0 ( 312/ 777- 353/ 903)  58.0 E49-E44 S21-S13
      6  215  761 (S22E62 L=277)   946     87        1 ( 164/ 738- 264/ 787)   3.8 E70-E55 S23-S20
      7 1653  519 (S31W54 L= 32)  1238    152        0 (1585/ 479-1703/ 555)  35.5 W47-W61 S33-S28
      8  524  509 (S39E38 L=301)   925    133        2 ( 472/ 479- 576/ 553)  20.0 E43-E32 S40-S35
      9 1200  447 (S40W17 L=356)   356    146        0 (1138/ 418-1262/ 496)  41.2 W13-W21 S42-S35
     10 1502  455 (S37W43 L= 21)   297     55        0 (1492/ 419-1516/ 489) -59.0 W40-W45 S39-S34
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y