mpas_tools.mesh.creation.signed_distance.signed_distance_from_geojson

mpas_tools.mesh.creation.signed_distance.signed_distance_from_geojson(fc, lon_grd, lat_grd, earth_radius, max_length=None, workers=-1)[source]

Get the distance for each point on a lon/lat grid from the closest point on the boundary of the geojson regions.

Parameters:
  • fc (geometrics_features.FeatureCollection) – The regions to be rasterized

  • lon_grd (numpy.ndarray) – A 1D array of longitude values

  • lat_grd (numpy.ndarray) – A 1D array of latitude values

  • earth_radius (float) – Earth radius in meters

  • max_length (float, optional) – The maximum distance (in degrees) between points on the boundary of the geojson region. If the boundary is too coarse, it will be subdivided.

  • workers (int, optional) – The number of threads used for finding nearest neighbors. The default is all available threads (workers=-1)

Returns:

signed_distance (numpy.ndarray) – A 2D field of distances (negative inside the region, positive outside) to the shape boundary