Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1105 1479 (N32W05 L=333)   519    134        0 (1054/1427-1145/1534) -53.1 W00-W07 N28-N36
      2 1728 1496 (N34W62 L= 30)   503     70        0 (1692/1474-1759/1516) -14.0 W57-W65 N31-N35
      3 1394 1448 (N31W26 L=354)   645    198        0 (1307/1392-1464/1483) -20.0 W19-W30 N26-N32
      4 1682 1421 (N29W51 L= 19)   878    138        5 (1644/1372-1725/1480)  66.8 W48-W57 N25-N32
      5 1365 1396 (N27W23 L=351)   335    103        1 (1315/1368-1412/1415) -15.9 W19-W26 N25-N28
      6  985 1280 (N19E03 L=325)   235     94        0 ( 941/1267-1022/1296)  11.3 E06-W00 N18-N19
      7 1695 1179 (N13W47 L= 15)   781    115        2 (1644/1154-1756/1209)   5.7 W42-W52 N11-N15
      8 1597  845 (S07W39 L=  7)   577    166        1 (1570/ 788-1647/ 914) -45.0 W36-W43 S10-S03
      9 1730  780 (S12W51 L= 19)   516     66        1 (1701/ 751-1756/ 809) -38.7 W48-W54 S13-S09
     10 1409  606 (S23W28 L=356)  2992    555       11 (1162/ 455-1551/ 744)  39.8 W11-W36 S34-S13
     11 1189  676 (S18W11 L=339)   200     65        0 (1160/ 655-1221/ 694)  26.6 W09-W13 S19-S17
     12 1558  648 (S20W38 L=  7)   315     44        0 (1543/ 621-1574/ 675)  90.0 W37-W39 S21-S18
     13 1079  643 (S21W05 L=333)   216     69        1 (1051/ 624-1112/ 663) -26.6 W02-W06 S22-S19
     14  321  590 (S28E56 L=272)  1206    147        5 ( 292/ 522- 364/ 650)  33.7 E64-E48 S33-S23
     15  507  530 (S31E38 L=290)  2027    275        3 ( 389/ 491- 636/ 587)  17.1 E52-E26 S34-S26
     16 1392  452 (S34W30 L=358)   635    210        0 (1301/ 434-1489/ 479)   7.6 W22-W37 S35-S31
  

    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