Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1013 1865 (N51E25 L=293)   751     75        2 ( 974/1841-1057/1889) -11.3 E30-E20 N48-N53
      2  745 1728 (N33E36 L=282)  1093    239        2 ( 690/1614- 803/1816) -67.6 E40-E33 N23-N40
      3  952 1686 (N33E19 L=299)   359     87        1 ( 934/1642- 968/1735) -49.4 E22-E16 N30-N37
      4 1488 1517 (N33W22 L=340)   333     69        1 (1449/1500-1527/1534)  26.6 W19-W25 N31-N34
      5  597  934 (S19E25 L=293)   762    209        1 ( 539/ 853- 663/ 991)  55.0 E28-E21 S24-S14
      6  408  584 (S44E42 L=276)  1268    174        2 ( 359/ 541- 462/ 639) -26.6 E47-E35 S46-S40
      7 1715  485 (S22W67 L= 25)  2574    317        3 (1595/ 414-1817/ 552)  41.2 W58-W78 S29-S15
      8 1286  488 (S33W30 L=348)   264     66        0 (1252/ 470-1326/ 504)  26.6 W28-W32 S34-S31
      9 1098  217 (S57W36 L=354)   887    123        2 (1034/ 183-1146/ 248)  -6.3 W26-W45 S59-S53
  

    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