Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1448 1746 (N55W05 L= 38)  1412    243        3 (1333/1686-1547/1788)   7.4 E06-W17 N51-N58
      2 1745  995 (N14W46 L= 78)   680    208        2 (1722/ 900-1776/1077) -68.2 W41-W48 N08-N18
      3  284  961 (S29E50 L=343)   954    204        1 ( 195/ 912- 365/1023)  26.6 E60-E42 S33-S23
      4  565  970 (S22E26 L=  7)   707    141        6 ( 526/ 928- 619/1015)  33.7 E29-E22 S24-S19
      5 1889  980 (N19W62 L= 95)   633     70        1 (1854/ 948-1909/1007) -36.9 W57-W66 N17-N20
      6  735  945 (S19E14 L= 18)   201     97        0 ( 714/ 904- 760/ 988) -45.0 E16-E11 S20-S16
      7  647  933 (S22E20 L= 13)   233     70        1 ( 628/ 903- 671/ 966) -26.6 E21-E17 S23-S19
      8  853  791 (S24E03 L= 29)   929    340        1 ( 689/ 758- 993/ 816)  34.5 E12-W04 S30-S19
      9 1591  691 (S07W44 L= 77)   358     95        1 (1563/ 656-1622/ 735) -26.6 W41-W47 S08-S05
  

    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