Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1255 1800 (N51E01 L=289)   815    127        0 (1206/1762-1313/1826) -16.7 E06-W04 N48-N53
      2  676 1782 (N36E46 L=245)   415     57        0 ( 661/1743- 692/1814) -36.9 E50-E42 N32-N38
      3 1104 1638 (N33E08 L=283)   394     88        0 (1091/1591-1118/1684) -66.8 E09-E06 N29-N36
      4 1170 1380 (N17W02 L=293)   286     83        1 (1155/1332-1184/1424) -63.4 W00-W03 N13-N19
      5 1051 1248 (N07E03 L=288)   539    225        0 ( 963/1199-1149/1308)  48.4 E06-W02 N02-N11
      6 1746 1277 (N24W44 L=335)   411     73        2 (1713/1259-1782/1301)  29.7 W41-W48 N22-N26
      7  934  773 (S23E01 L=290)   322    126        1 ( 885/ 749- 986/ 811)  31.0 E03-W02 S24-S20
      8  976  558 (S35W07 L=298)   445    136        0 ( 921/ 542-1055/ 574)  32.0 W03-W11 S37-S32
      9  653  532 (S44E17 L=273)   573     96        0 ( 623/ 495- 687/ 572) -38.7 E19-E14 S46-S40
     10 1193  247 (S51W41 L=332)   354     21        0 (1175/ 225-1213/ 265) -53.1 W38-W43 S53-S49
     11 1225  155 (S57W61 L=352)  1002    135        1 (1167/ 118-1270/ 188)  -3.9 W48-W78 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