Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1481 1736 (N39W52 L=191)   431     63        1 (1443/1717-1509/1759)  11.3 W47-W57 N37-N40
      2  439 1663 (N50E47 L= 91)  1059    142        2 ( 400/1605- 477/1727) -24.9 E66-E38 N44-N57
      3 1165 1601 (N34W18 L=157)   469    102        0 (1126/1569-1213/1642) -51.3 W16-W20 N31-N37
      4 1033 1423 (N24W06 L=144)   237     73        0 ( 995/1395-1065/1442) -45.0 W03-W07 N21-N25
      5  592 1309 (N22E24 L=114)   580    136        0 ( 541/1266- 654/1338)  23.2 E28-E20 N19-N23
      6  406 1254 (N21E39 L=100)   648    189        2 ( 326/1226- 503/1280) -14.9 E46-E31 N17-N23
      7 1740  909 (S16W47 L=186)  1637    303        2 (1706/ 766-1763/1057)  86.8 W44-W49 S24-S06
      8  237  696 (S10E63 L= 76)  1796    199        4 ( 187/ 646- 286/ 770)  18.4 E72-E54 S13-S06
      9 1737  557 (S39W55 L=193)  1637    244        2 (1648/ 501-1825/ 631)  14.0 W44-W65 S42-S34
  

    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