mpas_tools.mesh.creation.signed_distance.distance_from_geojson

mpas_tools.mesh.creation.signed_distance.distance_from_geojson(fc, lon_grd, lat_grd, earth_radius, nn_search='kdtree', 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

  • nn_search ({'kdtree'}, optional) – The method used to find the nearest point on the shape boundary

  • 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:

distance (numpy.ndarray) – A 2D field of distances to the shape boundary