Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  793 1688 (N40E19 L=170)   368     76        1 ( 756/1666- 831/1707) -15.9 E22-E15 N37-N41
      2  480 1684 (N40E46 L=143)   692    127        1 ( 435/1657- 517/1705) -16.7 E51-E41 N37-N41
      3  224 1218 (N09E57 L=132)   864    196        4 ( 199/1162- 261/1292) -76.0 E60-E52 N05-N13
      4  171 1128 (N03E63 L=127)  2057    383       11 ( 125/1024- 213/1216) -73.3 E70-E57 S02-N08
      5 1867 1049 (N01W61 L=250)   881    185        1 (1827/ 972-1910/1109) -37.9 W56-W67 S03-N04
      6 1000  863 (S13E01 L=188)   430    138        0 ( 966/ 791-1045/ 930)  63.4 E03-W01 S17-S09
      7  331  736 (S21E49 L=141)   428    147        0 ( 290/ 686- 366/ 778)  32.0 E53-E44 S23-S18
      8  497  557 (S33E39 L=151)   416    137        1 ( 453/ 529- 552/ 580)  12.5 E42-E33 S34-S31
  

    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