Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  902 1670 (N48E10 L=336)   287     37        0 ( 873/1655- 930/1685)   9.5 E13-E07 N46-N48
      2 1052 1611 (N43W03 L=349)   676    146        1 ( 975/1573-1113/1650) -21.8 E03-W08 N39-N45
      3 1649 1419 (N27W49 L= 35)   822     78        1 (1620/1369-1671/1464) -90.0 W46-W51 N23-N29
      4  791 1273 (N19E15 L=331)   998    279        0 ( 719/1164- 886/1368) -47.3 E20-E08 N12-N25
      5 1437 1166 (N12W27 L= 13)   997    216        0 (1361/1098-1521/1231) -38.7 W22-W32 N07-N15
      6  494 1166 (N12E35 L=311)   471     75        1 ( 470/1129- 520/1206) -76.0 E37-E33 N10-N14
      7 1215 1131 (N10W12 L=358)   217     64        0 (1195/1102-1238/1165) -53.1 W10-W13 N08-N12
      8 1186  761 (S13W10 L=356)   297    138        4 (1157/ 689-1205/ 805)  90.0 W08-W11 S17-S09
      9 1295  779 (S12W17 L=  3)   279     42        1 (1270/ 753-1316/ 804) -45.0 W15-W18 S13-S10
  

    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