Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  529 1731 (N71E32 L=108)   908     86        2 ( 487/1699- 574/1750) -12.7 E51-E20 N66-N74
      2 1087 1638 (N41W26 L=166)   277     68        2 (1048/1619-1116/1659)  -9.5 W22-W28 N39-N42
      3 1476 1562 (N24W47 L=188)   556    194        0 (1442/1507-1524/1616)  69.4 W45-W53 N19-N27
      4 1648 1502 (N15W57 L=198)   744     62        3 (1620/1464-1674/1535)  14.0 W53-W61 N13-N16
      5 1252 1367 (N20W23 L=163)   653    319        2 (1100/1319-1374/1403) -38.2 W14-W28 N14-N25
      6  707 1264 (N29E13 L=127)   271    118        0 ( 672/1237- 753/1309) -56.3 E15-E10 N26-N32
      7  235 1279 (N40E59 L= 81)   530     63        2 ( 207/1252- 264/1302) -16.7 E64-E54 N37-N41
      8  587 1178 (N26E23 L=117)   318     79        2 ( 554/1152- 630/1209) -45.0 E25-E20 N24-N28
      9  220  958 (N20E60 L= 81)   833     94        1 ( 176/ 923- 248/ 999)  26.6 E66-E56 N18-N22
     10  229  876 (N15E60 L= 80)   406     70        1 ( 202/ 850- 246/ 906)  18.4 E64-E57 N13-N17
     11  362  870 (N13E47 L= 94)   668    133        2 ( 304/ 850- 429/ 903) -21.8 E51-E41 N10-N15
     12  573  809 (N05E32 L=108)   325    118        1 ( 548/ 756- 605/ 858)  45.0 E34-E28 N02-N07
     13  857  805 (S01E15 L=125)   261     48        0 ( 830/ 786- 886/ 828) -45.0 E16-E13 S02-N00
     14  693  586 (S10E31 L=109)   262     79        2 ( 665/ 552- 727/ 610)  18.4 E34-E28 S11-S09
     15  540  530 (S11E44 L= 96)   398     52        1 ( 524/ 495- 559/ 561)  63.4 E46-E42 S12-S08
     16 1344  494 (S34W06 L=146)   567     75        2 (1309/ 463-1376/ 528) -36.9 W03-W07 S36-S30
     17  376  460 (S12E64 L= 77)  1936    210        3 ( 314/ 381- 427/ 520) -78.7 E69-E58 S18-S08
  

    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