Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  623 1758 (N28E54 L=353)   617     71        1 ( 593/1721- 647/1789) -15.3 E59-E48 N25-N30
      2  459 1720 (N22E67 L=340)  1947    155        3 ( 394/1675- 515/1758)  53.1 E74-E62 N17-N25
      3  583 1603 (N17E45 L=  1)  1357    269        3 ( 516/1515- 649/1691) -22.4 E54-E37 N11-N22
      4 1673 1409 (N37W33 L= 80)  2906    466        2 (1558/1202-1788/1575) -34.2 W21-W46 N27-N46
      5  378 1431 (N02E53 L=354)  1039    201        1 ( 302/1400- 485/1464)  23.2 E60-E45 N00-N06
      6  905 1303 (N07E14 L= 32)  1090    284        2 ( 784/1246-1017/1372)  52.4 E19-E09 N01-N14
      7 1579 1261 (N24W27 L= 74)   626    168        0 (1509/1208-1643/1321)  -9.5 W21-W33 N22-N25
      8  424 1170 (S11E40 L=  6)   335     75        2 ( 403/1128- 441/1205) -63.4 E42-E38 S13-S09
      9  270  994 (S25E51 L=356)  1064    200        2 ( 203/ 951- 345/1052) -10.6 E58-E42 S26-S22
     10 1062  995 (S06W03 L= 49)   270     79        0 (1038/ 964-1095/1031)  78.7 W01-W03 S07-S02
     11 1828  684 (N00W67 L=113)  1744    259        2 (1770/ 608-1857/ 793) -84.3 W60-W76 S04-N06
     12 1622  444 (S19W61 L=108)   562    102        1 (1587/ 405-1651/ 480) -68.2 W58-W63 S21-S16
  

    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