Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  541 1606 (N15E48 L=271)   334     74        2 ( 519/1589- 565/1632) -21.8 E50-E45 N13-N16
      2 1366 1136 (N08W16 L=335)   486    119        1 (1322/1072-1394/1192) -39.8 W12-W19 N05-N10
      3  597 1136 (S11E28 L=291)   371     85        0 ( 553/1113- 631/1160)  21.8 E30-E25 S13-S09
      4  464 1087 (S17E35 L=283)   315     74        0 ( 428/1070- 500/1103)   0.0 E38-E32 S17-S15
      5 1709  973 (N10W43 L=  1)   425     93        1 (1692/ 919-1726/1015) -80.5 W40-W44 N06-N12
      6 1563  899 (N01W34 L=353)   829    264        0 (1481/ 779-1605/ 997) -29.7 W26-W40 S04-N04
      7  844  733 (S28E03 L=316)   873    311        3 ( 708/ 699-1007/ 767)  17.5 E11-W08 S31-S25
      8 1710  656 (S06W55 L= 14)   331     65        0 (1683/ 643-1739/ 670)   9.5 W52-W58 S07-S04
  

    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