Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1501 1577 (N24W49 L= 92)   320     77        0 (1468/1548-1532/1604) -78.7 W48-W50 N21-N26
      2 1661 1447 (N12W54 L= 97)   482     57        1 (1626/1423-1694/1466) -21.8 W51-W57 N10-N13
      3  186 1225 (N37E64 L=339)  1091    143        0 ( 122/1194- 244/1255)   0.0 E79-E55 N35-N38
      4  335 1134 (N29E45 L=358)   350     67        0 ( 302/1113- 367/1160) -26.6 E47-E41 N27-N31
      5  721  887 (N07E21 L= 22)   762    240        2 ( 685/ 787- 802/ 948)  33.7 E24-E14 N01-N09
      6 1684  793 (S25W36 L= 79)   443    114        0 (1643/ 772-1736/ 817) -26.6 W32-W40 S27-S22
      7 1479  769 (S19W20 L= 63)   501     70        1 (1454/ 732-1499/ 815)  53.1 W18-W22 S21-S16
      8  711  563 (S12E31 L= 12)   306     71        1 ( 679/ 548- 749/ 588) -14.0 E32-E28 S13-S10
      9  436  507 (S10E53 L=350)   850    142        2 ( 389/ 457- 469/ 569)  32.5 E60-E48 S12-S05
     10  676  525 (S13E34 L=  8)   254     78        0 ( 639/ 504- 719/ 540) -21.8 E36-E31 S15-S11
     11  728  487 (S17E33 L= 10)   476    108        2 ( 685/ 463- 766/ 529)  15.9 E35-E28 S18-S14
     12  345  420 (S15E71 L=332)   817     59        1 ( 325/ 385- 361/ 448)  51.3 E75-E66 S17-S12
     13  623  395 (S20E44 L=358)   365    110        0 ( 583/ 372- 677/ 433) -51.3 E46-E40 S22-S17
  

    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