Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1608 1566 (N22W59 L= 36)  1445    246        3 (1571/1457-1644/1674)  30.1 W49-W68 N16-N28
      2 1408 1547 (N27W41 L= 18)   926    150        2 (1345/1487-1457/1604) -90.0 W38-W43 N22-N32
      3 1055 1370 (N26W11 L=348)   484    211        0 ( 997/1299-1090/1453) -90.0 W08-W12 N21-N32
      4 1395 1337 (N16W30 L=  8)   694    237        0 (1295/1265-1483/1394) -54.0 W26-W34 N09-N21
      5 1214 1365 (N22W20 L=357)   237     37        1 (1191/1343-1231/1390)  45.0 W18-W22 N20-N23
      6  656 1109 (N20E21 L=316)   519    146        0 ( 635/1032- 683/1167)  69.4 E22-E18 N15-N23
      7  132 1029 (N23E74 L=263)   985     75        1 (  98/1008- 169/1056)  -6.7 E84-E66 N21-N24
      8 1463 1012 (S04W26 L=  3)   258     90        0 (1431/ 975-1490/1053)  26.6 W22-W28 S05-S02
      9  216  812 (N09E64 L=273)   560    121        0 ( 182/ 761- 249/ 862) -54.5 E67-E60 N05-N12
     10 1311  711 (S18W10 L=347)   924    395        3 (1151/ 683-1454/ 778) -20.2 W00-W19 S21-S12
     11 1559  698 (S26W27 L=  4)   971    154        3 (1491/ 644-1636/ 742)   4.8 W21-W33 S27-S23
     12  762  272 (S35E43 L=295)  2471    353        4 ( 716/ 137- 816/ 458)  31.5 E71-E27 S49-S22
  

    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