Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  432 1542 (N30E53 L=144)  1626    178        4 ( 365/1482- 510/1589) -15.3 E64-E42 N26-N32
      2  314 1021 (S04E49 L=148)   322     60        1 ( 294/ 998- 350/1042)  21.8 E51-E45 S05-S02
      3  287  879 (S13E52 L=145)   549     97        0 ( 253/ 837- 315/ 918)  39.8 E56-E48 S15-S10
      4  222  833 (S16E60 L=137)   434     62        0 ( 205/ 802- 237/ 867) -63.4 E61-E58 S18-S14
      5  482  745 (S20E35 L=162)   430    158        0 ( 420/ 717- 518/ 794)  20.6 E40-E32 S21-S16
      6 1229  739 (S14W15 L=212)   634    218        0 (1129/ 708-1323/ 769)  -4.4 W08-W21 S15-S12
      7 1711  649 (S17W55 L=252)  1022    184        2 (1640/ 604-1761/ 678)  29.1 W48-W60 S20-S15
      8 1639  568 (S23W50 L=248)   616     89        4 (1609/ 520-1674/ 597) -35.5 W47-W55 S26-S21
      9  784  360 (S46E15 L=183)   496     73        1 ( 747/ 334- 816/ 387)  35.5 E18-E11 S47-S42
     10  550  342 (S50E41 L=156)   411     59        1 ( 515/ 323- 577/ 371) -18.4 E44-E37 S52-S47
  

    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