Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1386 1622 (N32W43 L= 81)   454     65        2 (1357/1597-1411/1648)  78.7 W40-W45 N29-N34
      2 1541 1522 (N22W49 L= 86)   430     44        1 (1514/1495-1570/1543)   9.5 W45-W51 N20-N23
      3 1790 1440 (N09W68 L=105)  1665    155        3 (1742/1397-1854/1484)  -9.0 W59-W78 N06-N12
      4 1766 1267 (N01W55 L= 93)  1089    226        4 (1727/1165-1792/1364)  68.2 W52-W59 S04-N07
      5 1292 1289 (N16W22 L= 60)   331    160        0 (1239/1234-1346/1316)   6.3 W17-W26 N13-N16
      6  912 1138 (N16E04 L= 33)   566    174        0 ( 833/1099- 987/1177) -45.0 E08-W00 N12-N19
      7  762 1127 (N19E13 L= 24)   199     76        0 ( 733/1104- 803/1159)  11.3 E15-E10 N17-N19
      8  632 1129 (N21E22 L= 16)   399    142        0 ( 561/1103- 693/1156) -29.1 E26-E17 N18-N24
      9 1664  935 (S14W38 L= 75)   735    235        1 (1621/ 837-1695/1042)  48.0 W32-W41 S18-S08
     10  879 1009 (N10E09 L= 29)   435    164        2 ( 841/ 966- 942/1039) -31.0 E11-E05 N07-N11
     11  487  865 (N09E37 L=  1)   756    257        0 ( 419/ 762- 521/ 976) -77.9 E39-E34 N02-N16
     12 1524  795 (S18W25 L= 62)   704    229        1 (1443/ 761-1617/ 833) -26.6 W20-W32 S19-S13
     13  875  649 (S11E17 L= 20)   311     80        2 ( 835/ 630- 922/ 673) -38.7 E19-E14 S13-S09
     14 1669  541 (S41W36 L= 74)  3330    458        1 (1499/ 415-1852/ 671)   8.5 W19-W60 S45-S37
     15  442  488 (S13E53 L=344)  1789    217        6 ( 414/ 399- 486/ 608)  57.3 E59-E47 S19-S05
     16 1288  308 (S46E01 L= 36)   676    149        1 (1216/ 288-1347/ 343)  -5.2 E07-W04 S48-S44
  

    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