Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1577 1614 (N22W56 L=241)  1779    338        4 (1470/1586-1683/1639) -19.3 W47-W67 N17-N24
      2  685 1375 (N35E12 L=174)   469    147        0 ( 615/1334- 759/1408)   0.0 E17-E05 N33-N36
      3  214 1299 (N42E54 L=131)   690    199        1 ( 169/1277- 263/1314) -11.3 E59-E49 N41-N43
      4 1129 1287 (N18W13 L=198)   264     70        1 (1092/1259-1170/1311) -53.1 W10-W14 N15-N19
      5  517 1000 (N17E31 L=154)   301     59        0 ( 491/ 967- 542/1034) -76.0 E32-E29 N15-N19
      6 1207  593 (S23E01 L=184)   207     54        0 (1186/ 575-1235/ 613) -26.6 E02-W00 S24-S21
      7  634  459 (S16E40 L=145)   869    178        3 ( 583/ 403- 667/ 512)  76.0 E45-E37 S20-S12
  

    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