Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  677 1519 (N22E28 L= 88)   293     85        1 ( 659/1483- 695/1559) -51.3 E30-E25 N19-N24
      2  600 1434 (N16E31 L= 84)   377     58        1 ( 575/1410- 625/1460) -36.9 E33-E29 N14-N17
      3 1615 1177 (N10W36 L=152)   817    250        0 (1570/1072-1680/1252)  76.0 W33-W41 N03-N15
      4 1456 1135 (N06W25 L=141)   439    104        0 (1431/1089-1490/1194) -56.3 W23-W27 N03-N09
      5  162  986 (S12E63 L= 52)  2554    317       10 ( 118/ 863- 200/1104)  70.3 E69-E58 S19-S05
      6  379  908 (S17E42 L= 74)  1071    362        1 ( 334/ 752- 441/1056)  72.5 E46-E37 S26-S07
      7  822 1009 (S08E12 L=104)   789    238        2 ( 722/ 984- 936/1055)   8.7 E18-E05 S09-S04
      8 1737  775 (S12W52 L=167)   620    248        2 (1656/ 742-1811/ 803)   3.8 W44-W59 S12-S10
      9  748  718 (S26E15 L=101)   198     75        1 ( 732/ 684- 760/ 761) -78.7 E16-E14 S28-S23
     10  912  693 (S26E04 L=112)   373    111        0 ( 864/ 658- 973/ 721)  39.8 E06-W00 S28-S23
     11 1622  565 (S26W48 L=164)   777    262        2 (1489/ 493-1715/ 626)  35.5 W40-W54 S32-S20
  

    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