Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  803 1492 (N39E09 L=127)   625    119        2 ( 749/1464- 848/1523)   7.1 E13-E05 N36-N40
      2  950 1203 (N18E02 L=134)   527    228        1 ( 877/1118- 996/1270) -68.2 E06-W00 N12-N22
      3  859 1178 (N18E08 L=128)   311    135        0 ( 842/1118- 886/1237) -90.0 E08-E07 N13-N21
      4  236 1192 (N24E60 L= 76)   455     70        1 ( 211/1151- 259/1213)  63.4 E63-E56 N21-N25
      5  706  948 (N06E21 L=115)   303    143        0 ( 671/ 886- 743/1008)  45.0 E23-E17 N02-N08
      6 1884  954 (S13W65 L=201)   867    112        1 (1857/ 906-1908/1007)  35.5 W61-W69 S16-S10
      7 1710  649 (S29W46 L=182)   403     42        3 (1688/ 626-1734/ 672)  71.6 W43-W48 S30-S27
      8  566  453 (S25E42 L= 93)  1502    361        2 ( 453/ 331- 676/ 583)  23.4 E66-E29 S33-S16
      9  395  477 (S22E58 L= 78)   440     44        1 ( 375/ 456- 424/ 501) -26.6 E60-E54 S23-S20
     10  927  342 (S38E19 L=117)   487    130        2 ( 850/ 323- 981/ 369)   5.2 E25-E14 S38-S36
  

    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