Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  479 1602 (N39E48 L=237)   597    110        1 ( 454/1566- 503/1631) -29.7 E51-E44 N36-N41
      2  217 1315 (N18E64 L=220)   562    108        1 ( 198/1293- 237/1335) -18.4 E67-E61 N16-N18
      3  714 1193 (N13E20 L=264)   949    317        2 ( 592/1141- 818/1290)  31.6 E28-E13 N09-N19
      4  185  961 (S04E64 L=220)   773    110        1 ( 160/ 928- 215/1006) -33.7 E68-E60 S06-S01
      5  659  870 (S07E23 L=262)   483    157        0 ( 606/ 834- 710/ 896)   8.1 E26-E19 S09-S05
      6 1858  821 (S10W65 L=350)   800    131        1 (1823/ 786-1889/ 854)  24.0 W61-W70 S11-S07
      7  198  790 (S15E64 L=221)   482    113        1 ( 173/ 774- 227/ 811)   0.0 E67-E59 S16-S13
      8  980  640 (S21E02 L=282)   989    254        2 ( 919/ 559-1030/ 773)  70.3 E06-W00 S26-S12
      9  694  703 (S18E21 L=263)   229     45        0 ( 667/ 683- 716/ 726) -45.0 E22-E19 S19-S16
     10  678  570 (S27E23 L=261)   255     89        0 ( 661/ 530- 693/ 607)  63.4 E24-E21 S30-S24
     11  423  393 (S43E53 L=232)  1308    275        1 ( 375/ 370- 486/ 412)  14.9 E60-E45 S45-S40
  

    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