Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  464 1523 (N18E45 L=324)   885    206        3 ( 396/1494- 558/1564)   9.5 E51-E37 N15-N20
      2  870 1408 (N15E15 L=355)  1865    365       11 ( 760/1304- 992/1499)  36.9 E21-E07 N07-N19
      3  647 1380 (N10E28 L=342)   254     90        2 ( 626/1336- 665/1414) -68.2 E29-E25 N07-N12
      4 1278 1070 (N00W14 L= 24)   185     50        0 (1255/1045-1299/1094) -33.7 W12-W15 S01-N01
      5  634 1028 (S11E24 L=346)   377    126        0 ( 601/ 983- 658/1080)  80.5 E25-E22 S13-S07
      6 1122 1070 (S02W05 L= 15)   178    119        1 (1056/1050-1169/1080)  26.6 W01-W07 S04-S01
      7  806 1032 (S08E13 L=357)   249     96        1 ( 782/ 995- 829/1079) -78.7 E14-E11 S10-S05
      8 1510 1002 (S00W29 L= 39)   223     71        0 (1498/ 962-1523/1043)  90.0 W28-W30 S02-N02
      9  159  900 (S22E65 L=305)  1696    218        6 (  88/ 841- 217/ 964)  16.3 E81-E56 S24-S17
     10  585  788 (S25E25 L=344)  1198    319        2 ( 436/ 743- 729/ 840)  10.8 E36-E15 S29-S21
     11 1312  685 (S22W23 L= 33)  1097    401        4 (1200/ 569-1482/ 813)  57.7 W18-W31 S30-S11
     12  998  770 (S21W02 L= 12)   679    288        3 ( 866/ 745-1138/ 796)  17.4 E06-W10 S24-S18
     13 1543  635 (S20W41 L= 50)   829    166        3 (1506/ 569-1574/ 728) -66.0 W36-W43 S24-S15
     14  980  557 (S35W04 L= 14)   896    251        1 ( 859/ 537-1103/ 578)   6.7 E04-W13 S36-S33
     15 1660  523 (S24W56 L= 65)   337     35        1 (1636/ 508-1686/ 540)  11.3 W53-W59 S25-S23
     16 1404  289 (S45W51 L= 60)  1760    212        5 (1315/ 245-1482/ 341)  -2.4 W38-W62 S46-S41
     17  949  287 (S55W10 L= 19)   333     78        0 ( 923/ 262- 986/ 315) -20.6 W06-W14 S57-S53
     18 1547  260 (S42W73 L= 82)   769     69        0 (1533/ 222-1562/ 299) -15.5 W65-W85 S44-S39
  

    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