Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1287 1801 (N58E07 L=218)   869     90        2 (1264/1752-1315/1843) -65.6 E10-E03 N52-N63
      2 1014 1641 (N34E19 L=205)   691    231        1 ( 972/1545-1074/1728) -33.7 E25-E10 N28-N39
      3 1538 1548 (N44W24 L=248)   688     97        0 (1505/1503-1570/1584) -38.7 W20-W27 N40-N46
      4 1137 1481 (N27E05 L=219)   227     44        2 (1116/1454-1158/1503) -18.4 E06-E03 N24-N28
      5 1544 1333 (N28W26 L=250)   400    121        0 (1524/1273-1573/1390) -50.2 W23-W28 N25-N31
      6 1820 1208 (N29W54 L=278)  1600    306        0 (1729/1086-1902/1304) -22.2 W43-W66 N23-N33
      7  214  770 (S37E61 L=163)   811     83        1 ( 187/ 738- 245/ 803)  26.6 E66-E55 S39-S34
      8  276  577 (S50E60 L=164)   633     69        1 ( 252/ 550- 299/ 608)  18.4 E68-E53 S52-S47
      9 1340  445 (S29W39 L=263)  2497    458        5 (1129/ 381-1534/ 512)  39.3 W26-W50 S37-S19
     10  425  425 (S59E42 L=182)   440     29        1 ( 408/ 401- 441/ 450)  32.0 E46-E38 S61-S56
     11  776  255 (S63W09 L=233)   371     83        0 ( 740/ 234- 816/ 282)  -4.8 W03-W15 S63-S60
  

    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