Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  856 1694 (N32E27 L=249)  1053    219        0 ( 807/1603- 896/1788) -40.9 E35-E20 N25-N38
      2  548 1752 (N31E56 L=219)   919     85        0 ( 507/1723- 588/1783)  38.7 E59-E53 N28-N34
      3 1462 1729 (N51W19 L=295)   737     67        2 (1433/1700-1498/1755) -38.7 W16-W23 N47-N53
      4 1484 1551 (N36W21 L=296)   338     74        0 (1463/1509-1505/1590) -59.0 W18-W22 N33-N38
      5 1862 1175 (N21W57 L=333)  1370    142        4 (1801/1124-1909/1218)   8.1 W50-W64 N18-N25
      6 1151  697 (S23W15 L=290)   451    127        0 (1103/ 673-1198/ 737)  38.7 W11-W16 S24-S19
      7 1213  509 (S33W25 L=300)   758    282        0 (1106/ 477-1361/ 560)  35.5 W17-W32 S36-S26
      8 1698  522 (S20W62 L=338)   496     58        0 (1675/ 492-1720/ 551) -76.0 W60-W64 S21-S17
      9  822  478 (S44E02 L=273)   604     98        1 ( 790/ 438- 859/ 529) -56.3 E05-W01 S46-S40
  

    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