Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  895 1896 (N52E39 L=292)   984     93        2 ( 855/1875- 944/1919)  -8.1 E46-E31 N48-N54
      2  653 1798 (N37E49 L=282)   764     97        2 ( 630/1755- 679/1848) -42.0 E56-E44 N33-N42
      3 1222 1451 (N22W04 L=335)   296     99        1 (1183/1419-1258/1477) -18.4 W01-W07 N20-N23
      4 1858 1086 (N15W58 L= 28)   558     69        2 (1828/1058-1882/1116)  38.7 W54-W60 N12-N16
      5  430  952 (S20E37 L=293)   509    127        1 ( 400/ 893- 474/ 999)  60.3 E39-E34 S23-S16
      6 1703  651 (S14W54 L= 25)   421     77        1 (1675/ 631-1735/ 680)  53.1 W52-W57 S15-S11
      7  334  581 (S45E52 L=279)   702     60        1 ( 310/ 551- 356/ 616) -90.0 E56-E47 S46-S42
      8 1699  558 (S19W59 L= 30)  1104    166        4 (1626/ 507-1765/ 604)  35.0 W52-W64 S22-S14
      9 1121  539 (S34W17 L=347)   449    152        0 (1039/ 522-1195/ 561)  21.8 W11-W21 S35-S30
     10 1635  477 (S25W59 L= 29)   729     96        1 (1604/ 444-1684/ 499)  51.3 W55-W63 S27-S22
     11 1000  254 (S57W22 L=353)   750    179        0 ( 906/ 228-1072/ 293)   0.0 W10-W30 S58-S55
  

    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