Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  882 1691 (N48E15 L=122)  2860    425        7 ( 776/1564-1023/1753)  10.3 E23-E01 N37-N54
      2 1448 1662 (N47W37 L=174)   715    116        1 (1404/1637-1487/1683)  16.7 W32-W42 N44-N48
      3 1305 1607 (N42W22 L=159)   760    140        2 (1254/1553-1343/1659)  53.1 W17-W27 N38-N46
      4 1583 1605 (N42W48 L=185)   473     81        1 (1561/1577-1611/1626) -33.7 W45-W50 N39-N43
      5 1664 1549 (N37W53 L=190)   466     69        1 (1636/1535-1693/1568) -11.3 W49-W55 N36-N38
      6  459 1425 (N27E43 L= 94)   678    207        3 ( 441/1341- 482/1512) -65.6 E45-E39 N21-N32
      7 1319 1433 (N29W20 L=157)   292     64        1 (1291/1401-1342/1466) -76.0 W18-W21 N27-N31
      8 1002 1307 (N21E02 L=135)   423     83        0 ( 967/1281-1061/1336) -31.0 E04-W01 N19-N22
      9  926 1225 (N15E07 L=130)   605    265        2 ( 797/1201-1063/1255)   3.4 E15-W02 N13-N17
     10  880 1141 (N10E09 L=128)   382    169        0 ( 832/1088- 980/1195) -31.0 E12-E02 N06-N13
     11 1597 1072 (N06W38 L=175)   379     86        1 (1566/1036-1635/1110) -45.0 W35-W40 N04-N08
     12 1294  896 (S04W17 L=154)   300     95        0 (1282/ 843-1312/ 945)  90.0 W16-W18 S07-S01
     13  317  803 (S13E50 L= 87)  1184    316        6 ( 237/ 726- 378/ 903) -42.3 E57-E44 S17-S07
     14  274  782 (S15E54 L= 83)   371     62        0 ( 254/ 760- 298/ 812) -33.7 E56-E52 S15-S12
     15 1010  685 (S18E00 L=137)   545    149        1 ( 975/ 621-1040/ 764)  84.3 E02-W01 S22-S12
     16  493  395 (S42E45 L= 93)  1514    246        2 ( 429/ 355- 541/ 460)  33.7 E52-E37 S45-S35
     17 1663  445 (S35W55 L=192)   542    145        1 (1630/ 436-1705/ 460)  -6.3 W50-W59 S36-S33
     18 1323  374 (S41W26 L=163)  1633    350        3 (1182/ 328-1498/ 429)  13.5 W15-W40 S44-S36
  

    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