Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1150 1645 (N37W23 L= 52)   566    178        1 (1074/1617-1233/1663) -21.8 W16-W27 N33-N38
      2  483  478 (S22E49 L=340)  1234    185        4 ( 421/ 421- 558/ 536) -60.9 E53-E45 S26-S16
      3  575  489 (S22E41 L=348)   316     75        1 ( 546/ 458- 608/ 518) -76.0 E42-E38 S24-S20
      4  773  367 (S34E33 L=356)  1185    257        4 ( 701/ 294- 864/ 474) -86.4 E35-E30 S41-S25
      5 1365  443 (S42W13 L= 42)   254     63        0 (1345/ 416-1389/ 473) -63.4 W11-W14 S44-S38
      6 1474  440 (S45W23 L= 52)   316     61        0 (1453/ 417-1497/ 471) -78.7 W20-W25 S47-S42
      7  626  209 (S43E63 L=327)   736     77        0 ( 605/ 176- 652/ 252)  32.5 E69-E56 S46-S39
  

    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