Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1017 1909 (N52E45 L= 90)  1022    123        0 ( 982/1860-1043/1945) -25.5 E58-E37 N46-N56
      2 1076 1837 (N47E30 L=105)   490     62        0 (1065/1801-1091/1872) -33.7 E35-E26 N44-N50
      3 1425  918 (S00W25 L=161)  1024    346        3 (1341/ 796-1478/1062) -33.7 W16-W31 S05-N05
      4  525  858 (S29E27 L=109)   249     92        0 ( 478/ 838- 560/ 878)  38.7 E29-E24 S31-S26
      5 1712  649 (S05W56 L=192)   324     28        0 (1692/ 632-1733/ 669) -14.0 W54-W59 S05-S03
      6  382  434 (S65E49 L= 86)   486     48        0 ( 363/ 412- 398/ 464)  22.4 E58-E41 S67-S60
      7 1496  367 (S27W58 L=194)   538    137        0 (1434/ 335-1558/ 392)  49.4 W55-W62 S30-S23
      8 1255  340 (S37W43 L=178)   842    155        2 (1171/ 310-1313/ 373)  49.4 W38-W46 S41-S33
      9  760  279 (S62W14 L=149)   421     80        1 ( 737/ 245- 796/ 314) -18.4 W08-W20 S64-S59
     10  881  214 (S62W33 L=168)   579     80        2 ( 849/ 186- 922/ 244)  90.0 W28-W38 S65-S58
  

    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