Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  221 1207 (N34E55 L= 54)   419     65        1 ( 201/1184- 245/1233) -56.3 E57-E52 N32-N35
      2  244  978 (N20E55 L= 55)   351    114        1 ( 213/ 955- 282/1007) -26.6 E57-E51 N18-N21
      3 1793  743 (S32W47 L=156)  2188    389        3 (1684/ 621-1898/ 879)  20.6 W35-W59 S36-S26
      4  953  680 (S11E12 L= 97)   917    308        4 ( 898/ 594-1021/ 785)  42.0 E17-E06 S15-S05
      5 1595  640 (S31W28 L=137)   538     81        1 (1566/ 592-1619/ 684)  63.4 W25-W30 S34-S28
      6 1126  560 (S23E06 L=104)   223    108        0 (1090/ 531-1172/ 605) -66.8 E06-E03 S25-S18
      7  303  568 (S06E64 L= 46)   478     90        1 ( 283/ 539- 321/ 597)  38.7 E66-E61 S07-S03
      8  614  407 (S21E44 L= 66)   437    141        1 ( 588/ 361- 646/ 456) -90.0 E46-E41 S24-S17
      9  707  222 (S37E50 L= 60)   440    138        0 ( 674/ 205- 738/ 256)  24.0 E53-E44 S38-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