Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  798 1647 (N30E32 L= 98)   398     70        2 ( 772/1617- 823/1681) -38.7 E35-E29 N28-N32
      2  541 1516 (N16E44 L= 86)   710    211        1 ( 460/1462- 599/1555)  10.3 E50-E38 N13-N18
      3  785 1447 (N17E25 L=106)   478    165        1 ( 715/1407- 863/1465)  24.0 E29-E20 N14-N20
      4 1388 1382 (N27W15 L=145)   477     87        1 (1361/1335-1406/1418) -68.2 W12-W16 N24-N29
      5  899 1389 (N16E16 L=114)   212     61        1 ( 871/1370- 923/1413)  14.0 E17-E13 N14-N17
      6  560 1156 (S05E31 L= 99)   294     81        1 ( 531/1123- 594/1189)  53.1 E32-E29 S06-S02
      7  246 1070 (S16E56 L= 75)   336     95        1 ( 216/1038- 271/1097) -26.6 E59-E52 S17-S14
      8  292 1016 (S18E50 L= 80)   392     83        1 ( 270/ 985- 314/1067)  78.7 E52-E48 S19-S14
      9  618  972 (S14E23 L=107)   211     54        0 ( 602/ 944- 634/1003) -56.3 E24-E22 S15-S12
     10  511  829 (S25E29 L=101)   783    213        3 ( 432/ 771- 577/ 909)  53.1 E35-E24 S30-S18
     11 1369  527 (S23W35 L=166)   616     97        2 (1335/ 494-1402/ 583) -39.8 W32-W39 S25-S19
     12  717  446 (S46E07 L=123)  1209    235        4 ( 612/ 425- 833/ 473)  20.2 E17-W02 S49-S41
     13 1089  180 (S57W45 L=175)  1060    109        1 (1055/ 137-1118/ 235) -34.5 W37-W53 S62-S51
  

    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