Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  438 1607 (N59E36 L=205)   389     54        0 ( 412/1585- 466/1624)   7.1 E40-E31 N58-N60
      2 1819  891 (S25W51 L=292)  1350    222        2 (1755/ 804-1896/ 957)  12.5 W42-W61 S29-S23
      3 1699  734 (S31W37 L=278)  1112    254        1 (1651/ 683-1772/ 820) -90.0 W32-W44 S35-S24
      4  911  704 (S09E15 L=227)   276    123        2 ( 891/ 654- 941/ 764)  81.9 E15-E12 S12-S05
      5 1425  649 (S27W14 L=255)   440    132        0 (1389/ 587-1456/ 702)  35.5 W10-W17 S30-S24
      6 1441  568 (S33W13 L=255)   305     81        1 (1394/ 554-1483/ 581) -26.6 W10-W16 S34-S31
      7  610  304 (S27E54 L=188)  1658    253        3 ( 530/ 228- 660/ 380)  67.4 E64-E49 S32-S20
  

    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