Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1301 1810 (N51W54 L=170)  1619    218        3 (1203/1760-1395/1854)  -3.8 W39-W70 N48-N52
      2 1543 1703 (N36W64 L=181)   457     49        1 (1520/1679-1568/1723)  -7.1 W60-W69 N34-N37
      3 1638 1336 (N13W47 L=164)  1008    309        2 (1528/1293-1769/1390) -24.2 W40-W60 N08-N18
      4  334 1313 (N34E49 L= 68)   426    109        1 ( 306/1270- 365/1365)  50.2 E51-E46 N31-N37
      5 1044 1227 (N18W05 L=121)   461    207        1 ( 960/1197-1134/1255) -11.3 W00-W10 N16-N19
      6  656 1128 (N18E22 L= 95)   327     69        0 ( 634/1098- 682/1166) -63.4 E22-E19 N16-N21
      7 1693 1149 (N01W46 L=163)   295     98        0 (1649/1128-1731/1165) -15.9 W42-W49 S01-N01
      8 1803  982 (S12W53 L=170)  4284    459       14 (1727/ 797-1891/1105)  43.2 W46-W65 S21-S04
      9  716 1035 (N12E19 L= 98)   194     53        0 ( 688/1015- 747/1053) -26.6 E20-E16 N10-N13
     10  890  635 (S15E15 L=102)  1784    524        6 ( 753/ 480- 973/ 818) -78.2 E19-E11 S25-S01
     11 1213  729 (S14W06 L=123)   699    226        3 (1185/ 619-1241/ 815)  67.4 W03-W08 S21-S09
     12  502  700 (S05E39 L= 77)   293     95        0 ( 473/ 658- 538/ 731)  29.7 E42-E35 S07-S03
     13  392  555 (S14E54 L= 63)  1458    362        2 ( 343/ 417- 451/ 699)  78.7 E59-E48 S23-S03
     14 1071  430 (S32E08 L=108)   768    131        1 (1010/ 402-1133/ 461)   0.0 E13-E03 S33-S30
     15  660  407 (S27E38 L= 79)   781     81        2 ( 625/ 362- 696/ 445) -90.0 E40-E34 S30-S23
     16 1320  275 (S52W08 L=125)   813    117        2 (1275/ 247-1367/ 307)   5.7 W03-W13 S55-S48
  

    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