Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  773 1676 (N37E33 L=155)   327     63        3 ( 740/1658- 804/1695)  21.8 E34-E29 N34-N38
      2  399 1451 (N15E52 L=136)   405     73        2 ( 380/1420- 425/1481) -45.0 E54-E49 N13-N17
      3  313 1369 (N09E57 L=131)  1474    309        8 ( 258/1258- 354/1444) -45.0 E62-E50 N02-N13
      4  647 1159 (N02E25 L=162)   268     70        1 ( 624/1125- 673/1189) -31.0 E27-E22 N00-N03
      5  262 1039 (S12E53 L=134)  1658    256        3 ( 231/ 949- 327/1154) -53.1 E58-E46 S16-S04
      6  867 1072 (S00E10 L=178)   213     90        0 ( 816/1060- 913/1085)  18.4 E13-E07 S01-N01
      7 1690  741 (S06W50 L=238)  1053    269        3 (1656/ 656-1711/ 837) -65.6 W45-W54 S11-N00
      8 1290  432 (S33W31 L=219)   367    154        1 (1212/ 414-1351/ 460)   5.2 W24-W35 S34-S31
  

    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