Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1508 1774 (N35W59 L=175)  1080    263        1 (1477/1733-1542/1813)  15.3 W53-W64 N33-N37
      2  826 1645 (N45W02 L=118)   440    148        0 ( 752/1613- 901/1673)  -4.8 E04-W08 N44-N46
      3  409 1582 (N52E39 L= 77)   440    102        1 ( 377/1566- 445/1597) -20.6 E42-E34 N49-N53
      4  424 1409 (N38E34 L= 82)   316     59        1 ( 399/1381- 452/1428)  11.3 E36-E31 N36-N39
      5 1021 1333 (N20W06 L=123)   432    161        1 ( 939/1308-1103/1354) -21.8 W01-W11 N17-N22
      6 1825 1212 (S05W59 L=175)  3317    619        8 (1726/1022-1905/1350)  23.4 W44-W76 S13-N01
      7 1566 1295 (N06W39 L=155)   528    160        0 (1516/1211-1606/1340) -83.7 W36-W40 N00-N09
      8  181 1016 (N18E59 L= 57)  2463    505        9 ( 106/ 919- 308/1145) -41.6 E66-E48 N10-N26
      9  528  615 (S11E40 L= 76)   302    110        0 ( 492/ 570- 559/ 650)  29.7 E43-E36 S13-S09
     10  963  496 (S27E16 L=100)   890    243        2 ( 922/ 408- 996/ 613)  53.1 E20-E11 S32-S20
     11  470  393 (S24E55 L= 61)  1622    385        5 ( 448/ 292- 504/ 506)  43.0 E63-E47 S30-S16
     12  860  233 (S44E36 L= 80)   396     77        0 ( 845/ 205- 875/ 266)  45.0 E39-E34 S46-S41
  

    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