Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  582 1705 (N58E37 L= 63)  1918    336        0 ( 416/1680- 719/1726)   0.0 E67-E18 N55-N60
      2  496 1361 (N33E35 L= 65)  1282    341        1 ( 402/1224- 564/1462) -46.7 E45-E28 N23-N40
      3 1785 1363 (N14W63 L=162)   464     55        0 (1760/1341-1814/1388) -36.9 W60-W66 N12-N15
      4  375 1272 (N28E45 L= 55)   366    109        0 ( 352/1231- 413/1323) -40.6 E47-E40 N24-N31
      5 1736 1217 (N07W52 L=151)   296     42        0 (1717/1188-1752/1239) -76.0 W50-W53 N04-N08
      6 1892  792 (S23W68 L=168)  1396    168        1 (1847/ 744-1919/ 883)  45.0 W60-W74 S26-S17
      7 1544  816 (S14W31 L=130)   473     94        1 (1514/ 765-1573/ 856)  59.0 W28-W32 S16-S11
      8 1445  553 (S30W23 L=122)  2420    442        2 (1251/ 464-1625/ 644)  11.3 W07-W37 S33-S25
      9  543  603 (S15E38 L= 62)   709    160        0 ( 470/ 570- 623/ 639)   8.1 E44-E30 S16-S12
     10  689  206 (S48E51 L= 49)   894     83        0 ( 655/ 185- 730/ 223)  10.6 E60-E44 S50-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