Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1259 1799 (N51E01 L=289)   699    119        0 (1212/1760-1322/1825)  -5.7 E05-W05 N48-N53
      2  682 1782 (N36E46 L=244)   455     66        0 ( 664/1742- 700/1815) -36.9 E50-E42 N32-N38
      3 1113 1640 (N33E07 L=283)   336     80        0 (1100/1592-1125/1683) -66.8 E09-E05 N29-N36
      4  519 1626 (N22E48 L=242)   319     67        0 ( 485/1606- 553/1638)  -8.1 E51-E44 N20-N22
      5 1106 1292 (N10E00 L=290)   854    350        1 ( 972/1198-1193/1411)  63.4 E06-W03 N02-N18
      6 1737 1271 (N23W43 L=334)   356     40        0 (1706/1255-1763/1289)   9.5 W40-W46 N22-N24
      7  938  769 (S23E00 L=290)   323    140        0 ( 891/ 748- 997/ 800)  31.0 E03-W02 S24-S20
      8  986  557 (S35W08 L=298)   455    120        0 ( 929/ 541-1060/ 574)  24.0 W03-W12 S37-S32
      9  659  527 (S44E17 L=273)   558     95        1 ( 624/ 493- 699/ 563)  -9.5 E19-E13 S46-S41
     10 1196  245 (S51W42 L=332)   371     17        0 (1178/ 224-1217/ 262) -76.0 W39-W44 S53-S49
     11 1246  145 (S57W66 L=356)   626     71        1 (1209/ 118-1272/ 176)  78.7 W60-W77 S58-S53
  

    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