Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  588 1706 (N58E36 L= 63)  1774    318        4 ( 421/1678- 724/1737)  -1.2 E67-E17 N55-N60
      2  503 1367 (N33E34 L= 65)  1225    347        1 ( 414/1222- 574/1471) -46.7 E44-E27 N23-N40
      3 1794 1365 (N14W64 L=163)   457     50        1 (1768/1344-1822/1389) -26.6 W61-W67 N12-N15
      4  384 1272 (N28E44 L= 55)   369    113        0 ( 361/1231- 419/1328) -54.5 E46-E40 N24-N31
      5 1055 1192 (N16W04 L=103)   184     55        1 (1036/1163-1068/1228)  71.6 W02-W05 N14-N17
      6 1555  816 (S14W31 L=130)   476     95        1 (1528/ 764-1583/ 862)  59.0 W29-W33 S17-S11
      7 1895  781 (S24W69 L=168)  1085    121        2 (1856/ 743-1920/ 842)  35.5 W62-W76 S26-S20
      8 1424  537 (S31W21 L=120)  2433    558        1 (1184/ 402-1634/ 645)  14.0 W02-W38 S37-S26
      9  556  605 (S15E37 L= 62)   668    159        0 ( 482/ 572- 631/ 640)   8.7 E43-E30 S16-S12
     10  834  274 (S42E29 L= 70)   315     50        2 ( 804/ 255- 864/ 301) -11.3 E31-E26 S44-S40
     11  693  206 (S48E51 L= 48)   901     83        1 ( 656/ 186- 737/ 224)   7.6 E59-E43 S49-S45
  

    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