Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  893 1783 (N61W15 L=288)   694    140        1 ( 837/1754- 967/1803)  -7.1 W06-W22 N59-N63
      2 1371 1709 (N40W49 L=323)   440    103        0 (1335/1678-1412/1746)   8.1 W43-W57 N39-N41
      3  437 1623 (N58E49 L=225)   481     56        1 ( 409/1601- 467/1641)  15.9 E53-E44 N56-N59
      4 1053 1582 (N41W16 L=290)   633    176        1 (1018/1522-1122/1616)  11.3 W12-W22 N36-N42
      5 1221 1578 (N36W28 L=302)   258     36        0 (1191/1562-1251/1594) -11.3 W25-W30 N34-N36
      6 1724 1515 (N18W67 L=341)  1212    117        0 (1681/1478-1782/1543)   0.0 W59-W80 N16-N19
      7  313  831 (N07E52 L=222)   475     86        1 ( 270/ 802- 355/ 863)   6.3 E57-E48 N05-N08
      8 1408  673 (S22W17 L=291)   285     84        0 (1384/ 636-1433/ 714)  45.0 W15-W19 S24-S20
      9 1594  663 (S28W31 L=305)   350     65        1 (1574/ 631-1609/ 695)  45.0 W29-W32 S30-S25
  

    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