根据经纬度左边判断网格

根据经纬度判断网格的方法，通常会将地球分成若干个网格，然后根据给定的经纬度坐标来判断该坐标位于哪个网格中。

下面是一个简单的方法，将地球分成 $n$ 行 $m$ 列的网格，每个网格的宽度为 $\Delta\lambda$，高度为 $\Delta\varphi$。则可以通过以下公式来计算一个给定的经纬度所在的网格编号：

$$
\begin{aligned}
i &= \left\lfloor \frac{\lambda - \lambda_{min}}{\Delta\lambda} \right\rfloor \\
j &= \left\lfloor \frac{\varphi - \varphi_{min}}{\Delta\varphi} \right\rfloor \\
grid\_id &= i + j \times n
\end{aligned}
$$

其中，$\lambda$ 和 $\varphi$ 分别是给定的经度和纬度，$\lambda_{min}$ 和 $\varphi_{min}$ 分别是网格的最小经度和纬度，$\Delta\lambda$ 和 $\Delta\varphi$ 分别是每个网格的宽度和高度，$i$ 和 $j$ 分别是该经纬度所在网格的横纵坐标，$n$ 和 $m$ 分别是网格的列数和行数，$grid\_id$ 是该经纬度所在网格的编号。

下面是一个示例代码，其中假设将地球分成 $10$ 行 $10$ 列的网格，每个网格的宽度和高度都是 $1$ 度：

const earthRadius = 6371; // 地球半径，单位：千米

// 将角度转换为弧度
function toRadians(angle) {
  return angle * (Math.PI / 180);
}

// 计算两个经纬度坐标之间的距离
function distance(lat1, lon1, lat2, lon2) {
  const φ1 = toRadians(lat1);
  const φ2 = toRadians(lat2);
  const Δφ = toRadians(lat2 - lat1);
  const Δλ = toRadians(lon2 - lon1);

  const a = Math.sin(Δφ / 2) * Math.sin(Δφ / 2) +
          Math.cos(φ1) * Math.cos(φ2) *
          Math.sin(Δλ / 2) * Math.sin(Δλ / 2);
  const c = 2 * Math.atan2(Math.sqrt(a), Math.sqrt(1 - a));

  return earthRadius * c;
}

// 根据给定的经纬度坐标计算所在网格编号
function getGridId(lat, lon) {
  const latMin = -90; // 最小纬度
  const latMax = 90;  // 最大纬度
  const lonMin = -180; // 最小经度
  const lonMax = 180;  // 最大经度
  const n = 10; // 列数
  const m = 10; // 行数
  const Δλ = (lonMax - lonMin) / n; // 经度宽度
  const Δφ = (latMax - latMin) / m; // 纬度高度

  const i = Math.floor((lon - lonMin) / Δλ);
  const j = Math.floor((lat - latMin) / Δφ);
  const gridId = i + j * n;

  return gridId;
}

// 测试
const lat = 31.2304; // 纬度
const lon = 121.4737; // 经度
const gridId = getGridId(lat, lon);
console.log(gridId); // 输出：34

上面的代码中，`distance` 函数用来计算两个经纬度坐标之间的距离，`getGridId` 函数用来计算给定经纬度所在的网格编号。在实际应用中，可以根据具体情况来调整网格数量和大小，以达到更好的效果。