Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  886 1749 (N57W17 L= 39)   418     65        0 ( 858/1728- 916/1769) -78.7 W13-W20 N54-N59
      2 1535 1607 (N25W55 L= 77)   965    211        1 (1457/1550-1622/1645) -59.0 W49-W59 N19-N29
      3 1367 1553 (N28W39 L= 60)   305    103        0 (1315/1519-1398/1578) -68.2 W36-W40 N25-N30
      4 1344 1455 (N23W32 L= 54)   352    127        1 (1289/1430-1404/1471) -23.2 W28-W36 N21-N24
      5 1364 1381 (N18W30 L= 52)   232     80        1 (1341/1350-1387/1402)  14.0 W28-W32 N16-N19
      6 1179 1258 (N17W15 L= 37)   602    133        2 (1122/1213-1229/1296) -63.4 W12-W17 N12-N20
      7  430 1158 (N28E37 L=345)   520    129        1 ( 377/1124- 485/1189) -16.7 E41-E31 N25-N30
      8  888 1114 (N15E06 L= 16)   280     69        1 ( 866/1079- 908/1146)  45.0 E07-E03 N13-N17
      9 1204  964 (S01W09 L= 30)   946    249        2 (1151/ 865-1269/1078)  42.0 W04-W15 S05-N04
     10 1885  972 (S21W62 L= 83)  1287    211        4 (1847/ 868-1926/1059)  45.0 W55-W70 S26-S16
     11  469  999 (N18E36 L=346)   306     67        1 ( 457/ 965- 485/1035)  63.4 E36-E34 N15-N19
     12 1885  833 (S30W62 L= 84)   580     57        1 (1856/ 811-1912/ 853) -29.1 W56-W67 S31-S26
     13 1370  710 (S20W13 L= 35)   441    187        0 (1297/ 646-1412/ 737) -54.5 W08-W15 S24-S16
     14 1265  532 (S28W02 L= 24)   590    188        0 (1187/ 497-1350/ 580)  -5.2 E03-W08 S29-S26
     15 1555  541 (S37W24 L= 46)   696     76        2 (1519/ 506-1589/ 578)  31.0 W20-W27 S40-S34
     16 1408  415 (S42W10 L= 31)   380     50        1 (1385/ 396-1438/ 445)   0.0 W07-W12 S43-S39
     17 1314  360 (S43W01 L= 22)   395     55        0 (1291/ 336-1340/ 390) -90.0 W00-W02 S45-S40
     18  692  252 (S34E51 L=331)   782     77        0 ( 667/ 222- 728/ 291)  66.8 E54-E47 S37-S30
  

    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