Coverage for  / home / casatest / venv / lib / python3.12 / site-packages / casatasks / private / task_msuvbinflag.py: 32%

1 

2from casatasks.private.imagerhelpers.input_parameters import saveparams2last 

3from casatasks.private.imagerhelpers.msuvbinflag_algorithms import UVGridFlag 

4 

5from casatasks import casalog 

6import time 

7 

8@saveparams2last(multibackup=True) 

9def msuvbinflag( 

10 binnedvis, # input binned vis name (output of msuvbin) 

11 method, #= 'radial_mean_annular', 'radial_per_plane' 

12 nsigma, #=5 

13 doplot=False) -> None: 

14 

15 flagger=UVGridFlag(binnedvis, doplot) 

16 if method == 'radial_mean_annular': 

17 tic = time.time() 

18 flagger.flagViaBin_radial(sigma=nsigma) 

19 toc = time.time() 

20 casalog.post("msuvbinflag running time : %.2f " % (toc - tic), 

21 'DEBUG', 

22 'task_msuvbinflag') 

23 

24 elif method == 'radial_per_plane': 

25 tic = time.time() 

26 flagger.flag_radial_per_plane(sigma=nsigma) 

27 toc = time.time() 

28 

29 casalog.post("msuvbinflag running time : %.2f " % (toc - tic), 

30 'DEBUG', 

31 'task_msuvbinflag') 

32 else: 

33 print("At this time only 'radial_mean_annular' and 'radial_per_plane' are supported") 

34 

35 return 

36 

37 #elif method == 'regionalMean': 

38 # print("No bueno") 

39 # #tic = time.time() 

40 # #flagViaBin_regionalMean(binnedvis=binnedvis, sizeRegion=20, sigma=5,ignorPoint=True) 

41 # #toc = time.time() 

42 # #print("MSuvbinflag other method Running time:", toc - tic) 

43 #elif method == 'gradient': 

44 # print("No bueno") 

45 # #tic = time.time() 

46 # #flagViaBin_gradient(binnedvis=binnedvis, radius=1, sigma=0.3) 

47 # #toc = time.time() 

48 # #print("MSuvbinflag pead method Running time:", toc - tic) 

49 #elif method == 'median': 

50 # print("No bueno") 

51 # #tic = time.time() 

52 # #flagViaBin_median(binnedvis=binnedvis, sigma=5) 

53 # #toc = time.time() 

54 # #print("MSuvbinflag median method Running time:", toc - tic) 

55 #else: 

56 # print("No bueno")