Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1295 1643 (N44W26 L=345)   889    203        2 (1207/1603-1383/1687) -15.9 W17-W32 N41-N48
      2  467 1453 (N31E43 L=277)   332     77        1 ( 449/1410- 485/1492) -56.3 E45-E40 N28-N34
      3 1762 1207 (N12W55 L= 14)  1410    213        3 (1682/1138-1827/1292) -42.5 W48-W61 N07-N18
      4  891 1137 (N11E08 L=312)   447    150        1 ( 850/1092- 943/1197) -45.0 E10-E04 N08-N14
      5 1600 1154 (N10W39 L=359)   382     96        0 (1580/1106-1615/1197) -78.7 W37-W40 N07-N12
      6 1652 1001 (N00W42 L=  2)   338     76        0 (1611/ 986-1691/1012)   0.0 W38-W45 N00-N01
      7 1576  754 (S15W37 L=357)   369    136        0 (1554/ 690-1590/ 804)  81.9 W35-W37 S18-S11
      8 1466  735 (S15W29 L=348)   259     87        1 (1441/ 706-1499/ 773) -45.0 W26-W31 S17-S12
      9  455  582 (S24E43 L=276)   370     74        1 ( 432/ 549- 486/ 615)  45.0 E46-E40 S26-S21
  

    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