Fixes

Author: ghiggi

gpm/etc/colorbars/xradar.yaml

@@ -11,7 +11,7 @@ reflectivity:
   cbar:
     extend: both
     extendfrac: 0.05
-    label: Reflectivity [$dBZ$]
+    label: Reflectivity [dBZ]
   auxiliary:
     category: [reflectivity, radar]
 differential_reflectivity:
@@ -21,12 +21,12 @@ differential_reflectivity:
     bad_color: none
   norm:
     name: Norm
-    vmin: -1
-    vmax: 8
+    vmin: -0.4
+    vmax: 4
   cbar:
     extend: both
     extendfrac: 0.05
-    label: Differential Reflectivity [$dB$]
+    label: $Z_{\mathrm{DR}}$ [dB]
   auxiliary:
     category: [reflectivity, radar]
 differential_phase:
@@ -39,7 +39,7 @@ differential_phase:
     vmax: 180 #
   cbar:
     extend: neither
-    label: Total Differential Phase ($\phi_{dp}$)  [$°$]
+    label: $\phi_{\mathrm{DP}}$  [°] # Total Differential Phase
   auxiliary:
     category: [phase, radar]
 propagation_differential_phase:
@@ -52,7 +52,7 @@ propagation_differential_phase:
     vmax: 360 # can be more ... unbounded
   cbar:
     extend: neither
-    label: Propagation Differential Phase ($\psi_{dp}$)  [$°$]
+    label: $\psi_{\mathrm{DP}}$ [°] # Propagation Differential Phase
   auxiliary:
     category: [phase, radar]
 specific_differential_phase:
@@ -66,7 +66,7 @@ specific_differential_phase:
   cbar:
     extend: max
     extendfrac: 0.05
-    label: Specific Differential Phase ($K_{dp}$) [$°/km$]
+    label: $K_{\mathrm{DP}}$ [° km$^{-1}$]
   auxiliary:
     category: [phase, radar]
 cross_correlation_ratio:
@@ -79,7 +79,7 @@ cross_correlation_ratio:
     vmax: 1
   cbar:
     extend: neither
-    label: Correlation between $Z_{H}$ and $Z_{V}$ ($\rho_{hv}$)
+    label:  $\rho_{\mathrm{HV}}$ [-]
   auxiliary:
     category: [radar]
 clutter_correction:
@@ -108,7 +108,7 @@ specific_attenuation:
   cbar:
     extend: max
     extendfrac: 0.05
-    label: Specific_attenuation [$dB/km$]
+    label: k [$dB/km$]
   auxiliary:
     category: [radar]
 specific_differential_attenuation:
@@ -121,7 +121,7 @@ specific_differential_attenuation:
   cbar:
     extend: max
     extendfrac: 0.05
-    label: Specific Differential Attenuation ($A_{dp}$) [$dB/km$]
+    label: $A_{\mathrm{DP}}$ [dB km$^{-1}$]
   auxiliary:
     category: [radar]
 path_integrated_attenuation:
@@ -135,7 +135,7 @@ path_integrated_attenuation:
   cbar:
     extend: max
     extendfrac: 0.05
-    label: Path Integrated Attenuation ($PIA$) [$dB$]
+    label: PIA [dB]
   auxiliary:
     category: [radar]
 path_integrated_differential_attenuation:
@@ -149,7 +149,7 @@ path_integrated_differential_attenuation:
   cbar:
     extend: max
     extendfrac: 0.05
-    label: Path Integrated Differential Attenuation ($PIDA$) [$dB$]
+    label: $\delta$PIA [dB]
   auxiliary:
     category: [radar]
 linear_depolarization_ratio:
@@ -163,7 +163,7 @@ linear_depolarization_ratio:
   cbar:
     extend: min
     extendfrac: 0.05
-    label: Depolarization Ratio [$dB$]
+    label: Depolarization Ratio [dB]
   auxiliary:
     category: [radar]
 normalized_coherent_power:
@@ -190,7 +190,7 @@ signal_to_noise_ratio:
   cbar:
     extend: both
     extendfrac: 0.05
-    label: Signal to Noise Ratio [$dB$]
+    label: Signal to Noise Ratio [dB]
   auxiliary:
     category: [radar]
 TH:

gpm/gv/xradar/accessors.py

@@ -98,6 +98,7 @@ def plot_map(
         add_colorbar=True,
         add_background=True,
         extent=None,
+        rasterized=True,
         fig_kwargs=None,
         subplot_kwargs=None,
         cbar_kwargs=None,
@@ -113,6 +114,7 @@ def plot_map(
             add_colorbar=add_colorbar,
             add_background=add_background,
             extent=extent,
+            rasterized=rasterized,
             fig_kwargs=fig_kwargs,
             subplot_kwargs=subplot_kwargs,
             cbar_kwargs=cbar_kwargs,

gpm/gv/xradar/methods.py

@@ -319,6 +319,7 @@ def plot_map(
     add_background=True,
     add_gridlines=True,
     add_labels=True,
+    rasterized=True,
     fig_kwargs=None,
     subplot_kwargs=None,
     cbar_kwargs=None,
@@ -360,6 +361,7 @@ def plot_map(
         x=x,
         y=y,
         add_colorbar=False,
+        rasterized=rasterized,
         # cbar_kwargs=cbar_kwargs,
         **plot_kwargs,
     )
@@ -371,5 +373,5 @@ def plot_map(
         _ = plot_colorbar(p=p, ax=ax, **cbar_kwargs)
     # Set extent
     if extent is not None:
-        ax.set_extent(extent)
+        ax.set_extent(extent, crs=ccrs.PlateCarree())
     return p

gpm/utils/manipulations.py

@@ -1548,14 +1548,27 @@ def extract_within_point_distance(
     earth_radius = 6371.0 * 1000
     radius_rad = radius / earth_radius
 
-    ind = tree.query_radius(query_point, r=radius_rad)[0]
+    # Retrieve indices and distance
+    ind, dist = tree.query_radius(query_point, r=radius_rad, return_distance=True)
+
+    ind = ind[0]
+    dist = dist[0]  # in radians
+
+    # Convert to meters
+    dist_m = dist * earth_radius
 
     # Map back to flattened full-grid indices
     flat_indices = np.flatnonzero(valid)[ind]
 
     # Stack spatial dims, then select matching points
     ds_points = ds.stack(points=spatial_dims)
     ds_near = ds_points.isel(points=flat_indices)
+
+    # Add distance as coordinate
+    ds_near = ds_near.assign_coords(
+        distance=("points", dist_m),
+    )
+
     return ds_near