Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1382 1786 (N37W53 L= 97)   497     84        1 (1357/1759-1408/1804) -26.6 W50-W55 N35-N38
      2 1747 1504 (N10W65 L=109)  1573    258       10 (1683/1417-1788/1588) -85.2 W60-W70 N04-N16
      3 1189 1521 (N28W25 L= 69)   405    151        1 (1133/1487-1264/1537) -12.5 W20-W29 N25-N29
      4 1281 1333 (N15W23 L= 67)   199     51        0 (1255/1314-1300/1358)   0.0 W21-W24 N13-N15
      5 1061 1294 (N19W09 L= 53)   279     45        2 (1033/1264-1093/1320) -45.0 W07-W11 N16-N20
      6  217 1282 (N39E53 L=351)   681    217        1 ( 165/1262- 271/1305) -18.4 E59-E47 N37-N41
      7 1804  958 (S21W49 L= 92)   540    112        0 (1782/ 908-1825/1014)  80.5 W47-W50 S23-S17
      8 1484  717 (S25W19 L= 63)   798    227        1 (1452/ 600-1537/ 816)  59.7 W16-W24 S32-S20
      9 1157  595 (S23E04 L= 40)   814    304        1 (1069/ 492-1214/ 742)  38.7 E12-W03 S27-S15
     10  595  301 (S27E52 L=351)   636    145        2 ( 571/ 263- 621/ 342)  32.0 E55-E47 S30-S24
  

    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