Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1619 1539 (N20W57 L=187)   557    124        0 (1566/1519-1681/1556) -29.1 W52-W62 N17-N22
      2 1240 1485 (N28W26 L=156)   326     74        2 (1209/1456-1272/1515)  18.4 W23-W29 N27-N29
      3 1619 1458 (N16W52 L=182)   312     62        0 (1594/1432-1645/1478) -45.0 W50-W53 N14-N17
      4  191 1043 (N23E63 L= 67)  1623    210        5 ( 148/ 956- 234/1116) -51.3 E70-E57 N17-N27
      5 1848  944 (S19W57 L=187)   678     89        2 (1821/ 913-1881/ 989)  26.6 W53-W61 S21-S17
      6 1634  882 (S16W35 L=165)   403     92        1 (1608/ 845-1674/ 911)  18.4 W32-W38 S17-S14
      7  822  498 (S20E25 L=105)  1550    401        4 ( 688/ 363- 937/ 649)  25.7 E40-E13 S27-S12
      8  587  360 (S25E50 L= 80)  1803    354        1 ( 520/ 295- 692/ 489)  82.4 E55-E42 S31-S15
      9 1064  433 (S30E12 L=118)   273     98        0 (1018/ 417-1103/ 457)   0.0 E15-E09 S30-S29
  

    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