Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1563 1331 (N20W39 L=350)  1849    444        7 (1395/1249-1675/1413) -30.1 W28-W47 N13-N26
      2 1852 1293 (N14W67 L= 18)  1113    137        4 (1823/1244-1883/1360) -90.0 W63-W71 N11-N18
      3  889 1230 (N18E08 L=303)   203    119        0 ( 846/1171- 925/1282) -54.5 E10-E05 N14-N21
      4  584 1053 (N08E28 L=283)   236     36        1 ( 559/1035- 607/1072)  18.4 E29-E26 N06-N09
      5 1885  951 (S07W66 L= 17)  1305    292        2 (1835/ 910-1937/ 978) -16.4 W59-W76 S10-S04
      6 1563  758 (S16W35 L=346)   947    303        1 (1461/ 647-1660/ 896)  47.0 W27-W42 S22-S07
      7  178  831 (S06E68 L=243)   763    128        1 ( 159/ 799- 208/ 882) -68.2 E70-E63 S07-S02
      8 1647  559 (S30W45 L=356)  1421    241        2 (1549/ 517-1716/ 604)  21.0 W36-W51 S32-S27
      9  476  534 (S25E43 L=268)   766    154        1 ( 415/ 497- 544/ 585) -36.9 E47-E38 S27-S21
  

    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