Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1494 1501 (N18W44 L= 53)  1151    258        1 (1395/1475-1613/1522) -28.1 W37-W52 N13-N21
      2 1796 1397 (N02W65 L= 74)   894    239        2 (1756/1361-1849/1432)   0.0 W59-W75 N01-N03
      3 1604 1122 (S05W37 L= 46)   991    277        0 (1543/1051-1713/1225) -67.4 W33-W43 S11-N01
      4 1866 1022 (S20W56 L= 66)   452     82        2 (1844/ 987-1879/1058)  45.0 W53-W58 S21-S17
      5 1730  885 (S23W40 L= 50)  1466    322        4 (1696/ 794-1760/1000)  63.4 W36-W44 S28-S16
      6  151  976 (N23E64 L=306)   619    136        1 ( 117/ 956- 181/ 995)  -8.1 E67-E60 N21-N24
      7  208  869 (N16E60 L=309)   412     66        0 ( 186/ 844- 230/ 901) -76.0 E62-E57 N14-N18
      8 1509  718 (S25W21 L= 30)   257     67        0 (1497/ 689-1526/ 755)  78.7 W19-W21 S27-S22
      9  197  711 (N07E69 L=301)  1427    236        3 ( 150/ 677- 256/ 749) -22.6 E74-E62 N04-N09
     10  485  537 (S08E48 L=322)   292     46        0 ( 469/ 514- 503/ 563) -90.0 E49-E46 S10-S06
     11  373  483 (S10E61 L=308)   609    116        2 ( 343/ 451- 395/ 530) -90.0 E64-E58 S11-S06
  

    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