Attenuation correction artifact about pyart HOT 7 CLOSED

Can you add a plot of Phi_dp and K_dp?

from pyart.

@josephhardinee please see updated image.

from pyart.

Hmm.. might be that thin slice of 2nd trip right at the receiver.. Did you hardset the system phase? Or did you let the algorithm try to find it?

BTW this is one of the issues we hope to address with CMAC2.0

from pyart.

Does Kdp go negative in that ray?

from pyart.

@scollis, here's the phase processing code snippet:

pdp_corr, kdp_corr = pyart.correct.phase_proc_lp(
                       radar, -2.0, self_const=60000.0, low_z=0.0,
                       high_z=55.0, min_phidp=0.01, min_ncp=0.3, min_rhv=0.7, 
                       fzl=3900.0, overide_sys_phase=False, sys_phase=332.0,
                       nowrap=291, LP_solver='cylp_mp', window_len=35, proc=8,
                       debug=False, really_verbose=False)

@josephhardinee, no, please see the updated image for clarity. That being said, there certainly are negative pixels...

from pyart.

Hmm.. interesting..

One workaround to keep you going is to set the ray to _FillValue and
then it will not be used in the gridding.. Is it badvalue from first
gate? Could you add a rayplot?

On 6/30/14, 2:45 PM, Kirk North wrote:

@scollis https://github.com/scollis, here's the phase processing
code snippet:

pdp_corr, kdp_corr = pyart.correct.phase_proc_lp(
radar, -2.0, self_const=60000.0, low_z=0.0,
high_z=55.0, min_phidp=0.01, min_ncp=0.3, min_rhv=0.7,
fzl=3900.0, overide_sys_phase=False, sys_phase=332.0,
nowrap=291, LP_solver='cylp_mp', window_len=35, proc=8,
debug=False, really_verbose=False)

@josephhardinee https://github.com/josephhardinee, no, please see
the updated image for clarity.

—
Reply to this email directly or view it on GitHub
#176 (comment).

from pyart.

@scollis, yup it's _FillValue for the entire ray(s). The May 24th case suffered from second trip echoes much more than the other cases I've processed, so my first guess is these artifacts are linked to that. That being said, these artifacts decrease the more I lower ncp_min towards 0, i.e. the more I allow second trip pixels in. Furthermore, I see this artifact on several tilts at different azimuths.

from pyart.