Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  874 1866 (N42E47 L=149)  1162    141        3 ( 821/1821- 937/1911) -11.9 E58-E38 N39-N45
      2 1131 1636 (N32E13 L=183)  1043    208        2 (1072/1546-1177/1714) -32.7 E20-E06 N27-N36
      3  773 1566 (N17E31 L=164)   406    123        0 ( 743/1527- 820/1613) -20.6 E34-E26 N15-N20
      4  718 1497 (N12E31 L=164)   438    131        0 ( 676/1454- 770/1543) -18.4 E35-E26 N10-N13
      5  405 1327 (S05E47 L=148)   833    225        1 ( 366/1226- 433/1420) -79.7 E50-E44 S11-N00
      6 1556 1299 (N24W24 L=220)   780    212        0 (1519/1215-1578/1404) -61.4 W19-W26 N18-N30
      7 1731 1290 (N30W39 L=234)   627    117        1 (1693/1234-1755/1331) -51.3 W35-W41 N27-N32
      8  987 1278 (N06E09 L=186)   415    113        0 ( 953/1232-1035/1316)  71.6 E09-E07 N02-N09
      9  501 1238 (S08E37 L=158)   438    157        0 ( 453/1174- 526/1301)  90.0 E38-E35 S12-S04
     10  175  861 (S37E65 L=130)   805    123        0 ( 136/ 825- 217/ 908)  20.6 E75-E57 S39-S33
     11  992  675 (S27W09 L=204)   388     97        1 ( 948/ 655-1041/ 700)  18.4 W05-W11 S28-S25
     12 1468  352 (S29W58 L=253)   907    105        1 (1418/ 330-1527/ 383)  15.3 W52-W63 S30-S27
  

    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