Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1128 1878 (N59W02 L=302)   478     73        0 (1092/1863-1165/1892)   8.1 E02-W05 N57-N61
      2  758 1853 (N53E37 L=263)   673     53        1 ( 734/1823- 783/1882) -80.5 E40-E33 N49-N56
      3  972 1865 (N55E15 L=286)   439     65        1 ( 948/1840-1007/1880)  18.4 E17-E11 N52-N57
      4 1276 1596 (N34W14 L=314)   321     86        0 (1258/1553-1297/1634)  90.0 W12-W15 N30-N36
      5  739 1578 (N29E24 L=276)   337     60        1 ( 698/1561- 776/1590)   0.0 E27-E21 N27-N29
      6  652 1555 (N27E30 L=271)   316     70        0 ( 635/1518- 677/1585)  76.0 E31-E28 N24-N28
      7  880  858 (S15E08 L=293)   319    103        0 ( 848/ 818- 922/ 896)  56.3 E09-E05 S18-S12
      8  596  745 (S24E27 L=274)   752    159        1 ( 537/ 700- 686/ 779)  20.0 E31-E20 S26-S21
      9 1074  731 (S22W05 L=306)   193     85        2 (1020/ 708-1101/ 753)  31.0 W01-W07 S23-S20
     10  965  487 (S39E00 L=300)   692    105        3 ( 921/ 446-1000/ 534) -71.6 E03-W02 S41-S35
     11 1503  378 (S41W49 L=350)   455     75        0 (1478/ 359-1547/ 403)   0.0 W45-W53 S42-S39
  

    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