Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1172 1686 (N39W33 L=143)   333     48        2 (1142/1665-1199/1714) -26.6 W30-W34 N37-N40
      2 1278 1653 (N33W38 L=148)   553    126        2 (1235/1597-1318/1682)  81.9 W34-W39 N28-N35
      3 1551 1627 (N23W56 L=166)   416     97        1 (1516/1590-1580/1662)  80.5 W54-W57 N20-N26
      4  809 1626 (N47W04 L=114)   387    103        1 ( 764/1595- 856/1647) -14.0 W00-W08 N45-N47
      5 1339 1571 (N27W37 L=147)   298     44        0 (1314/1551-1359/1596)  14.0 W34-W38 N25-N28
      6  425 1301 (N35E34 L= 76)   435    128        0 ( 391/1247- 474/1348) -53.1 E37-E30 N30-N38
      7  786 1280 (N25E08 L=102)   337    131        1 ( 764/1231- 815/1346) -90.0 E09-E06 N21-N28
      8  425 1159 (N26E35 L= 75)   368    112        2 ( 380/1131- 475/1186)   0.0 E39-E30 N25-N28
      9  241 1067 (N25E54 L= 56)   794    123        2 ( 208/1027- 302/1107)   5.2 E58-E47 N22-N27
     10  139 1002 (N23E69 L= 41)   881     87        2 ( 114/ 961- 167/1041) -45.0 E73-E64 N20-N25
     11  193  800 (N10E65 L= 45)   420     96        1 ( 170/ 763- 218/ 849) -80.5 E67-E62 N07-N13
     12  948  771 (S08E11 L=100)   603    173        1 ( 910/ 690-1000/ 847)  41.2 E14-E06 S11-S03
     13 1742  647 (S39W43 L=153)   414     64        2 (1714/ 621-1770/ 670)   9.5 W39-W45 S40-S36
     14  987  462 (S27E18 L= 92)   427    100        0 ( 945/ 440-1046/ 488) -45.0 E20-E14 S29-S24
     15  746  400 (S25E36 L= 74)   406    112        0 ( 693/ 382- 802/ 420) -14.0 E39-E31 S25-S23
     16 1138  193 (S54E25 L= 85)   439     41        0 (1111/ 176-1163/ 213) -26.6 E27-E21 S55-S51
     17 1044  128 (S57E42 L= 68)   563     54        0 (1017/ 110-1065/ 149) -76.0 E47-E36 S59-S54
  

    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