Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  353 1282 (N21E48 L=356)  1000    187        1 ( 263/1252- 427/1316)   0.0 E58-E41 N18-N22
      2 1810 1268 (N15W61 L=106)   786     86        2 (1774/1227-1831/1304) -38.7 W57-W64 N12-N16
      3 1653 1025 (N01W42 L= 86)   353    127        0 (1641/ 959-1667/1089) -90.0 W41-W42 S02-N05
      4 1320 1039 (N04W18 L= 63)   191     82        0 (1280/1013-1355/1062)  31.0 W15-W20 N02-N05
      5  432  973 (N02E39 L=  5)   543    189        0 ( 362/ 924- 513/1026) -20.0 E44-E33 S01-N05
      6 1404  951 (S02W23 L= 68)   243     87        0 (1391/ 908-1421/ 997) -78.7 W22-W24 S04-N01
      7 1831  929 (S06W59 L=103)   842    107        1 (1770/ 912-1887/ 958)  -8.1 W52-W66 S07-S04
      8  972  911 (S03E04 L= 41)   307     75        0 ( 944/ 880-1005/ 949)  53.1 E05-E01 S04-N00
      9  617  846 (S06E26 L= 18)   537    177        0 ( 547/ 813- 662/ 907)  32.0 E31-E23 S07-S02
     10 1311  782 (S12W17 L= 62)  1031    421        1 (1136/ 677-1433/ 896)  32.3 W06-W25 S17-S05
     11  306  713 (S15E54 L=351)   824    205        1 ( 256/ 621- 347/ 796) -70.0 E57-E50 S20-S09
     12  392  648 (S19E47 L=358)   316     89        0 ( 372/ 609- 419/ 682)  39.8 E49-E43 S21-S16
  

    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