Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  354 1503 (N53E51 L=  0)   471     76        0 ( 317/1485- 392/1519)  -4.4 E57-E44 N51-N54
      2  181 1101 (N28E65 L=346)   559     67        1 ( 152/1080- 216/1126) -15.3 E70-E59 N26-N29
      3  785 1005 (N12E14 L= 37)  1296    316        5 ( 703/ 871- 826/1124) -86.4 E17-E11 N03-N19
      4 1783  879 (S22W48 L= 99)  3471    646        2 (1575/ 667-1919/1103)  19.7 W27-W69 S31-S12
      5  493  953 (N14E35 L= 16)   398    164        0 ( 468/ 901- 551/ 988)  26.6 E36-E30 N11-N16
      6  192  916 (N17E64 L=347)   607    151        1 ( 176/ 855- 228/ 985) -63.4 E66-E60 N12-N20
      7 1614  834 (S19W32 L= 83)  1176    255        2 (1550/ 723-1682/ 921)  26.6 W25-W39 S24-S15
      8  722  729 (S03E25 L= 26)   467    114        0 ( 700/ 676- 748/ 790)  63.4 E26-E22 S06-N00
      9 1827  729 (S34W54 L=105)   649     88        0 (1794/ 695-1855/ 766)  23.2 W49-W58 S35-S32
     10  941  535 (S20E17 L= 34)   309     88        0 ( 897/ 515- 981/ 549) -36.9 E20-E15 S21-S18
     11 1149  342 (S39E12 L= 39)  1498    269        2 (1066/ 260-1225/ 444)  29.1 E21-E02 S44-S33
     12  860  377 (S28E29 L= 22)   406     89        0 ( 817/ 360- 907/ 400) -38.7 E31-E25 S30-S25
  

    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