Select points and cross sections

cross_section_line(data, start, end)	Obtain an interpolated cross-sectional slice through gridded data.
cross_section_linestring(data, linestring)	Obtain an interpolated cross-sectional slice through gridded data.
points_in_bounds(da, **points)	Returns whether points specified by keyword arguments fall within the bounds of da.
points_values(da[, out_of_bounds])	Get values from specified points.
points_set_values(da, values[, out_of_bounds])	Set values at specified points.
points_indices(da[, out_of_bounds])	Get the indices for points as defined by the arrays x and y.
get_lower_active_grid_cells(active)	Returns grid of booleans designating location of the lowermost active grid cell.
get_lower_active_layer_number(active)	Returns planar grid of integers with the layer number of the lower most active cell.
get_upper_active_grid_cells(active)	Returns grid of booleans designating location of the uppermost active grid cell.
get_upper_active_layer_number(active)	Returns planar grid of integers with the layer number of the lower most active cell.
grid_boundary_xy(grid)	Return grid boundary on the xy plane.
active_grid_boundary_xy(active)	Return active boundary cells on the xy plane.