Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  555 1521 (N48E27 L= 97)   353     36        1 ( 534/1497- 577/1543)  33.7 E29-E24 N46-N49
      2  660 1487 (N43E17 L=107)   493     96        0 ( 615/1459- 700/1524) -40.6 E21-E13 N40-N46
      3 1711 1395 (N12W57 L=181)   702    148        2 (1670/1344-1753/1436)  26.6 W52-W62 N09-N13
      4 1769 1375 (N09W62 L=186)   644     92        1 (1746/1340-1790/1427)  26.6 W58-W66 N07-N11
      5  156 1185 (N30E72 L= 52)   609     90        2 ( 124/1157- 189/1205)   4.1 E79-E65 N27-N30
      6 1850 1144 (S06W62 L=186)   895     98        1 (1805/1111-1888/1181) -33.7 W57-W67 S08-S02
      7  218 1049 (N21E60 L= 64)   792    113        1 ( 195/1009- 260/1102) -45.0 E63-E55 N18-N24
      8  573  779 (N00E33 L= 91)   777    170        1 ( 539/ 701- 630/ 834) -71.6 E34-E30 S05-N04
      9 1562  647 (S28W28 L=152)   493     60        1 (1542/ 608-1585/ 684)  68.2 W26-W30 S30-S25
     10  347  608 (S07E57 L= 67)   388     70        0 ( 332/ 573- 362/ 639)  45.0 E58-E54 S09-S05
     11  787  545 (S18E25 L= 99)   274     61        1 ( 773/ 519- 806/ 577) -90.0 E26-E23 S19-S15
  

    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