Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  420 1541 (N30E55 L=145)  1608    156        2 ( 356/1493- 492/1588) -12.8 E66-E44 N27-N32
      2  283  945 (S09E52 L=147)  1123    266        1 ( 239/ 847- 325/1053)  66.8 E57-E48 S15-S01
      3  459  755 (S19E37 L=162)   536    151        1 ( 398/ 728- 498/ 798)  24.0 E42-E33 S21-S16
      4  763  770 (S16E15 L=185)   224     29        0 ( 741/ 755- 789/ 787)  33.7 E16-E13 S16-S14
      5 1197  742 (S14W13 L=212)   615    230        1 (1093/ 707-1289/ 774)  -8.7 W06-W19 S15-S12
      6 1691  644 (S18W53 L=252)  1316    238        4 (1617/ 597-1744/ 685)  24.4 W46-W58 S21-S15
      7 1330  576 (S24W25 L=224)   235     92        0 (1298/ 552-1357/ 592)  26.6 W22-W26 S25-S23
      8  771  373 (S44E16 L=184)   384     49        0 ( 743/ 349- 803/ 400)  33.7 E18-E12 S46-S41
      9  528  346 (S50E44 L=155)   530     94        0 ( 480/ 326- 565/ 367) -12.5 E50-E39 S51-S48
  

    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