¿Es posible generar un cuadro delimitador alrededor de un conjunto de coordenadas en JMP?Todos,
Si tengo un conjunto de coordenadas, ¿existe una función geoespacial que pueda generar las coordenadas de un cuadro delimitador dentro del radio especificado?
Por ejemplo
Latitud = 47.80
Longitud = -122.00
Me gustaría generar un cuadro delimitador alrededor de estas coordenadas dentro de un radio de 10 millas. ¿Tiene JMP alguna función incorporada para esto?
Pude lograr esto usando "R", excepto que preferiría hacerlo en JMP si es posible. Además, la documentación de la función para la función acepta una fracción y no necesariamente la distancia, no está claro en lo que simboliza la fracción. Para eliminar la ambigüedad, prefiero hacerlo en JMP si es posible.
Lat <- 47.80
Lon <- -122.00
Rad <- 0.05
BBox <- ggmap::make_bbox(lon = Lon,lat = Lat,f = Rad)
Re: ¿Es posible generar un cuadro delimitador alrededor de un conjunto de coordenadas en JMP? Esta publicación te señalará a este sitio que describe la fórmula de largo alcance para calcular distancias entre puntos en el globo y otras fórmulas útiles. Utilicé las fórmulas en la sección. Punto de destino dada la distancia y el rumbo desde el punto de inicio para crear un círculo de puntos alrededor de la ubicación central. En mi publicación, utilicé una distancia que se alejaba 90 grados en todas las direcciones. 10 millas serán una fracción de un grado. (Al menos así creo que lo hice. Con un poco de suerte hay algunos comentarios en la JSL.)
Projected circle of sunlightCírculo proyectado de luz solar.
Re: ¿Es posible generar un cuadro delimitador alrededor de un conjunto de coordenadas en JMP?No estoy seguro de para qué pretende utilizar el cuadro delimitador, pero he usado algo como el siguiente para controlar el "zoom" de los mapas de Graph Builder mediante JSL (adaptación rápida del sistema métrico).
Names default to here(1);
// Desired map central position and distance to view boundary
Lat = 40.782222; //Latitude
Lon = -73.965278; //Longitude
Rad = 7; // miles
// Calculate bounding box coordinates with geodetic approximation (WGS84)
a = 6378137; // Radius of earth at equator (m)
e2 = 0.00669437999014; // eccentricity squared
m = 1609.344; // 1 mile in meters
r = Pi() / 180; // convert to radians
//Distance of 1° latitude (miles)
d1 = r * a * (1 - e2) / (1 - e2 * Sin(Lat * r) ^ 2) ^ (3 / 2) / m;
//Distance of 1° longitude (miles)
d2 = r * a * Cos(Lat * r) / Sqrt(1 - e2 * Sin(Lat * r) ^ 2) / m;
//Bounding box coordinates
{minLat, maxLat} = {Lat - Rad / d1, Lat + Rad / d1};
{minLon, maxLon} = {Lon - Rad / d2, Lon + Rad / d2};
//––––––––––––––––––––––––––––––––––––––––––––––––––
//Example map
dt = Open("$SAMPLE_DATA/SAT.jmp");
gb = dt << Graph Builder(Variables(X(:Longitude), Y(:Latitude)), Elements(Points(X, Y)));
Report(gb)[FrameBox(1)] << Background Map(Images("Street Map Service")); Wait(1);
//Set bounding box (Manhattan, ~13 miles long)
Report(gb)[ScaleBox(2)] << {Min(minLat), Max(maxLat)};
Report(gb)[ScaleBox(1)] << {Min(minLon), Max(maxLon)};
