Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1455 1516 (N32W18 L=222)   470     92        1 (1427/1470-1479/1558) -71.6 W15-W20 N29-N35
      2  981 1345 (N11E09 L=196)   593    159        0 ( 940/1262-1011/1408) -63.4 E12-E06 N05-N14
      3 1421 1378 (N22W17 L=222)   211     49        1 (1402/1349-1433/1404) -90.0 W15-W18 N20-N23
      4  793 1324 (N06E19 L=185)   248     52        2 ( 769/1293- 819/1346)  90.0 E20-E17 N04-N07
      5 1753 1292 (N26W43 L=248)   528     77        2 (1729/1263-1779/1331) -45.0 W41-W45 N23-N28
      6  220 1009 (S20E56 L=149)   494    105        0 ( 188/ 979- 251/1033)  31.0 E58-E52 S22-S18
      7  303  618 (S44E49 L=155)   487     94        1 ( 289/ 576- 321/ 648)  78.7 E51-E47 S47-S42
      8  368  550 (S48E44 L=160)   434     55        1 ( 350/ 522- 388/ 577)  76.0 E46-E42 S50-S46
  

    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