Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  608 1711 (N48E43 L=344)   562     85        2 ( 578/1686- 631/1742)  76.0 E45-E38 N45-N50
      2  854 1497 (N33E14 L= 13)   644    163        2 ( 788/1455- 922/1543) -24.4 E19-E08 N29-N35
      3  449 1387 (N23E43 L=344)   748    138        2 ( 426/1328- 468/1466) -83.7 E44-E40 N19-N28
      4 1393 1324 (N22W24 L= 51)  1715    177        3 (1223/1229-1559/1421) -25.6 W12-W35 N16-N28
      5  923 1313 (N21E08 L= 20)   866    215        4 ( 876/1220- 969/1390) -61.4 E11-E04 N14-N25
      6  727 1085 (N06E19 L=  8)  1407    371        4 ( 641/ 969- 879/1219) -34.5 E25-E08 S01-N13
      7  585  942 (S04E28 L=359)   253     78        0 ( 548/ 924- 620/ 955)  11.3 E30-E25 S04-S02
      8 1313  255 (S52W33 L= 60)   503    155        0 (1270/ 227-1355/ 279) -22.6 W27-W39 S54-S49
  

    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