Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  815 1858 (N55E48 L=164)   945     97        1 ( 783/1827- 848/1882) -81.9 E54-E41 N51-N58
      2  673 1724 (N40E45 L=167)   510     55        0 ( 649/1699- 703/1750) -12.5 E49-E40 N38-N41
      3 1529 1675 (N54W38 L=250)   451     69        0 (1496/1654-1565/1692)  17.1 W32-W45 N51-N56
      4 1395 1589 (N43W20 L=231)   358     99        0 (1355/1558-1431/1615) -26.6 W16-W22 N41-N44
      5 1309 1419 (N29W13 L=225)   430    142        2 (1242/1401-1377/1439)  21.8 W08-W18 N26-N31
      6 1563  455 (S27W52 L=264)   612    129        0 (1502/ 438-1615/ 472)  10.3 W46-W57 S29-S25
      7 1202  355 (S41W28 L=239)   321     63        1 (1174/ 340-1231/ 375)  -9.5 W24-W30 S42-S39
      8  757  211 (S66E11 L=200)   634     64        1 ( 733/ 187- 785/ 240) -56.3 E14-E06 S69-S62
  

    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