Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  591 1803 (N45E53 L=265)  1210    125        4 ( 540/1779- 643/1829)  -9.0 E64-E44 N42-N47
      2  841 1733 (N40E23 L=295)   380    108        1 ( 798/1707- 895/1752)  23.2 E25-E17 N37-N41
      3 1296 1707 (N43W14 L=332)   480     97        1 (1252/1680-1347/1733) -14.0 W10-W18 N41-N45
      4 1683 1560 (N37W49 L=  6)  1442    185        1 (1613/1496-1742/1614) -35.0 W43-W53 N33-N41
      5 1593 1582 (N37W39 L=357)   388     76        0 (1560/1560-1646/1599)  -7.1 W36-W44 N36-N38
      6  937 1346 (N14E08 L=309)   846    305        1 ( 889/1241-1018/1480) -57.3 E12-E02 N07-N21
      7  228 1427 (N15E66 L=251)   718    112        3 ( 200/1391- 261/1454) -13.0 E73-E60 N12-N17
      8 1292 1038 (S02W16 L=333)   937    232        1 (1230/ 937-1372/1119)  61.4 W13-W20 S08-N03
      9 1725 1007 (N01W46 L=  4)   303     38        1 (1698/ 991-1747/1029)   0.0 W44-W48 N00-N02
     10 1616  965 (S02W38 L=356)   344     88        1 (1589/ 925-1640/1005)  90.0 W36-W39 S05-N00
     11 1550  681 (S20W39 L=356)   280    103        0 (1492/ 668-1596/ 700)  20.6 W34-W42 S21-S18
     12  333  625 (S34E52 L=265)   567    119        1 ( 316/ 569- 352/ 682)  41.2 E56-E48 S37-S30
     13  915  552 (S35E03 L=314)  1440    387        0 ( 755/ 481-1036/ 682)  41.6 E15-W04 S41-S25
     14  513  609 (S35E35 L=283)   401     65        2 ( 477/ 589- 550/ 625)   0.0 E37-E31 S35-S33
     15  968  389 (S47W03 L=320)   253     41        0 ( 943/ 370- 994/ 407)  36.9 W01-W05 S48-S45
     16 1497  316 (S44W56 L= 13)   596     78        0 (1454/ 298-1533/ 334)  23.2 W51-W61 S46-S42
  

    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