Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1303 1646 (N45W27 L=345)   816    226        1 (1219/1605-1394/1688) -24.8 W19-W34 N41-N48
      2  480 1454 (N32E42 L=277)   405     72        1 ( 462/1416- 501/1490) -45.0 E44-E38 N29-N34
      3 1775 1208 (N12W56 L= 14)  1550    239        3 (1698/1138-1848/1298) -40.2 W49-W63 N07-N18
      4 1282 1215 (N15W17 L=335)   374    173        1 (1208/1168-1346/1256) -32.0 W12-W20 N11-N17
      5 1615 1152 (N10W40 L=359)   448    110        1 (1594/1094-1630/1199)  80.5 W39-W41 N06-N12
      6  918 1133 (N11E06 L=312)   473    146        0 ( 883/1082- 976/1193) -54.5 E08-E02 N07-N14
      7 1664 1003 (N00W43 L=  2)   376     90        1 (1624/ 986-1707/1017)   0.0 W39-W46 N00-N01
      8 1596  790 (S13W38 L=357)   818    266        0 (1555/ 659-1635/ 866)  77.0 W36-W41 S20-S07
      9 1663  795 (S13W44 L=  2)   269     28        0 (1643/ 781-1687/ 810)  18.4 W42-W46 S13-S11
     10  466  584 (S24E42 L=276)   352     70        1 ( 444/ 553- 488/ 614)  45.0 E44-E39 S26-S21
  

    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