Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1280 1410 (N20W26 L= 98)   550    159        1 (1248/1353-1321/1484)  32.0 W23-W31 N16-N23
      2 1834 1276 (S05W62 L=134)  2348    220        2 (1788/1182-1873/1379)  79.7 W58-W66 S10-N01
      3 1110 1207 (N14W10 L= 81)   304    120        0 (1045/1189-1168/1239) -15.9 W06-W13 N12-N15
      4  661 1079 (N19E20 L= 52)   870    262        1 ( 551/1045- 797/1142) -38.2 E25-E11 N13-N24
      5 1625 1038 (S10W36 L=107)   400    113        0 (1609/ 970-1642/1090)  56.3 W33-W38 S13-S07
      6  190  885 (N17E62 L= 10)   492     88        0 ( 179/ 841- 201/ 931)  71.6 E63-E59 N14-N20
      7 1607  782 (S24W29 L=100)   989    188        0 (1562/ 691-1634/ 875)  48.4 W23-W31 S28-S18
      8 1357  578 (S29W07 L= 79)  1225    275        1 (1247/ 513-1429/ 693)  28.3 W00-W14 S32-S23
      9  824  626 (S11E22 L= 49)   340    146        0 ( 755/ 596- 895/ 664)   0.0 E26-E16 S11-S09
     10  903  517 (S19E22 L= 50)   564    191        0 ( 881/ 441- 931/ 611)  45.0 E25-E16 S23-S14
  

    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