Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  603 1656 (N30E47 L=  8)   812    204        1 ( 588/1551- 624/1735) -45.0 E53-E39 N23-N36
      2  901 1481 (N24E17 L= 38)   219     60        0 ( 877/1455- 930/1500) -21.8 E19-E14 N23-N25
      3 1345 1428 (N30W14 L= 69)   483    234        0 (1248/1390-1437/1494) -11.3 W05-W20 N28-N32
      4  784 1436 (N19E23 L= 32)   512    138        2 ( 715/1411- 840/1464)  29.7 E27-E20 N16-N21
      5 1684 1245 (N25W42 L= 98)   394     32        1 (1661/1220-1705/1270) -11.3 W39-W44 N23-N26
      6  206 1219 (S04E65 L=351)   592     97        2 ( 179/1171- 230/1257) -32.0 E70-E60 S06-S01
      7 1722 1148 (N20W46 L=101)   624     80        1 (1674/1123-1760/1177)   0.0 W41-W49 N18-N21
      8  437  855 (S22E36 L= 19)  3155    698        8 ( 201/ 703- 683/1097)  36.5 E67-E21 S36-S02
      9 1773 1063 (N15W52 L=107)   415     45        0 (1737/1048-1805/1079)   8.1 W48-W55 N14-N16
     10 1208 1043 (N03W10 L= 66)   207     80        0 (1172/1014-1241/1068)  45.0 W08-W12 N01-N05
     11  181  974 (S19E64 L=352)  3327    274        7 ( 121/ 866- 229/1059)  72.9 E75-E57 S25-S12
     12 1722  929 (N06W49 L=104)   291     74        0 (1699/ 893-1751/ 959) -26.6 W46-W52 N04-N07
     13 1535  883 (S00W35 L= 90)   630    127        3 (1483/ 851-1593/ 925)  -6.3 W30-W39 S01-N01
     14 1265  700 (S16W21 L= 77)   344    117        2 (1218/ 676-1332/ 715)  29.7 W18-W25 S17-S13
     15 1726  582 (S14W62 L=118)  1658    274        2 (1639/ 495-1807/ 678)  63.4 W56-W67 S20-S06
     16  637  493 (S43E18 L= 37)   280     51        1 ( 616/ 471- 670/ 513)  31.0 E20-E15 S44-S41
     17 1501  387 (S31W53 L=109)  1185    133        2 (1449/ 356-1561/ 429)  16.7 W48-W61 S34-S28
  

    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