Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  917 1739 (N51W19 L=264)   371     89        1 ( 896/1692- 933/1784)  49.4 W16-W23 N47-N54
      2  288 1483 (N52E56 L=190)   501     58        0 ( 255/1471- 323/1497)  -4.4 E62-E49 N50-N52
      3  720 1353 (N31E10 L=235)   206     56        0 ( 686/1342- 752/1364) -14.0 E12-E08 N29-N32
      4  437 1246 (N31E33 L=212)   549    106        0 ( 393/1200- 485/1277)  12.5 E37-E28 N29-N32
      5 1589  855 (S20W30 L=275)   352    105        0 (1543/ 822-1626/ 884)   8.1 W26-W33 S21-S18
      6 1517  319 (S58W19 L=264)   671    104        0 (1467/ 297-1569/ 341)  -8.1 W12-W26 S60-S56
      7 1293  219 (S58E09 L=236)  2136    273        2 (1163/ 186-1415/ 255) -15.1 E22-W05 S62-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