Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1612 1719 (N24W73 L= 90)   951     88        0 (1571/1701-1663/1733) -11.3 W66-W82 N22-N26
      2  352 1219 (N34E41 L=337)   481    131        0 ( 299/1167- 393/1258) -41.6 E46-E37 N29-N37
      3  755 1087 (N16E13 L=  4)   386    129        0 ( 705/1040- 789/1137) -74.1 E15-E11 N13-N20
      4 1525 1050 (S07W29 L= 46)   323    108        2 (1491/1002-1559/1095)  23.2 W25-W32 S08-S04
      5  129 1052 (N28E70 L=308)   986    100        1 (  78/1022- 173/1074)  -5.4 E83-E62 N25-N29
      6  725 1026 (N14E17 L=  1)   194     51        1 ( 699/1003- 747/1045) -33.7 E18-E15 N11-N14
      7 1594 1000 (S12W32 L= 50)   252     79        0 (1562/ 965-1612/1034)  26.6 W29-W34 S12-S09
      8 1403  940 (S09W19 L= 36)   303    104        1 (1367/ 919-1457/ 974) -45.0 W16-W21 S11-S06
      9 1340  818 (S14W12 L= 30)   662    238        2 (1251/ 764-1390/ 906)   6.3 W06-W15 S16-S08
     10 1491  831 (S18W21 L= 39)   215     48        0 (1480/ 803-1505/ 866)  76.0 W20-W22 S19-S15
     11 1790  838 (S28W46 L= 63)   457    108        0 (1753/ 798-1837/ 862)   0.0 W41-W51 S29-S26
     12  889  721 (S07E15 L=  2)   864    279        0 ( 853/ 597- 932/ 834)  42.5 E21-E09 S13-S02
     13  642  373 (S22E45 L=333)  2357    428        6 ( 567/ 224- 786/ 531)  22.6 E64-E28 S31-S14
     14  384  481 (S10E59 L=318)   602     75        2 ( 362/ 445- 409/ 511)  23.2 E63-E55 S11-S07
     15 1426  465 (S39W09 L= 27)   444     73        0 (1406/ 431-1446/ 504) -81.9 W07-W10 S42-S35
     16  524  437 (S15E50 L=328)   495     58        2 ( 499/ 408- 546/ 472)  68.2 E52-E47 S17-S12
     17 1328  345 (S45E02 L= 15)  1181    150        3 (1265/ 301-1392/ 388)   0.0 E06-W04 S48-S41
     18 1244  274 (S48E12 L=  6)   312     31        0 (1229/ 254-1260/ 299)  78.7 E13-E09 S50-S45
  

    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