Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1527 1642 (N29W50 L=269)   410     69        0 (1486/1615-1557/1660) -36.9 W47-W52 N27-N31
      2  617 1358 (N27E21 L=197)  2168    484        4 ( 447/1242- 837/1450)  13.6 E34-E04 N22-N30
      3  287 1326 (N31E50 L=168)  1769    206        3 ( 187/1293- 394/1369) -14.5 E65-E38 N27-N35
      4 1797 1047 (S10W53 L=271)   386     69        0 (1778/1014-1814/1082) -90.0 W51-W54 S12-S08
      5 1731  923 (S17W45 L=264)   275     25        1 (1714/ 903-1752/ 940)  14.0 W43-W47 S17-S15
      6 1308  828 (S15W14 L=232)   232     69        0 (1272/ 806-1343/ 845)  11.3 W11-W16 S15-S14
      7 1490  751 (S23W25 L=244)   415     86        1 (1467/ 703-1511/ 798) -81.9 W23-W26 S26-S19
      8 1595  714 (S27W33 L=252)   578    103        0 (1542/ 685-1650/ 751)  14.0 W29-W38 S28-S25
      9 1276  688 (S23W10 L=228)   243     96        0 (1241/ 654-1324/ 705)   9.5 W07-W13 S24-S22
     10  411  494 (S21E55 L=164)   706    171        0 ( 335/ 474- 471/ 533)   3.6 E63-E47 S22-S19
  

    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