Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1635 1549 (N27W59 L= 46)   428     66        0 (1598/1527-1665/1567) -31.0 W55-W61 N25-N28
      2 1779 1496 (N20W74 L= 61)   721     67        1 (1763/1466-1798/1526)  23.2 W71-W80 N18-N21
      3  726 1396 (N33E15 L=332)   911    218        3 ( 680/1293- 776/1484) -65.2 E19-E12 N26-N39
      4 1102 1383 (N27W11 L=358)   339    138        0 (1063/1330-1150/1433) -69.4 W08-W13 N22-N30
      5 1704 1155 (N03W48 L= 35)   931    110        3 (1663/1113-1746/1205) -63.4 W43-W51 N00-N06
      6  128 1092 (N18E76 L=271)   832    117        0 ( 111/1044- 154/1145) -19.7 E83-E69 N15-N20
      7 1830 1083 (S04W59 L= 47)   471     76        0 (1814/1040-1842/1117)  53.1 W57-W61 S07-S02
      8  389  695 (S07E48 L=299)   902    184        2 ( 353/ 606- 422/ 762)  39.3 E54-E43 S12-S02
      9 1136  661 (S18W02 L=349)   298    154        0 (1069/ 636-1196/ 702)  14.0 E02-W06 S18-S15
     10  918  433 (S30E17 L=330)   289     95        1 ( 906/ 386- 930/ 482)  76.0 E18-E15 S34-S26
  

    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