国土地理院の標高APIの使い方：「-----」の返り方、hsrcの精度、大量の地点は標高タイルで

{"elevation":3.6,"hsrc":"1m（レーザ）"}

緯度と経度を逆に渡す（lat=139.767125&lon=35.681236） -> {"elevation":"-----","hsrc":"-----"}  HTTP 200
lon を省く（lat=35.681236 のみ）                      -> {"elevation":"-----","hsrc":"-----"}  HTTP 200

// dem-tile.mjs — 国土地理院の標高タイル（PNG）から1点の標高を読む。依存なし（Node.js 18 以降）
import { inflateSync } from 'node:zlib';

const UA = 'dem-tile-example/0.1 (+https://codeagent.jp/)';
// 精度の高い順。標高APIのhsrcの優先順に合わせる
export const SOURCES = [
  { id: 'dem1a_png', z: 17, hsrc: '1m（レーザ）' },
  { id: 'dem5a_png', z: 15, hsrc: '5m（レーザ）' },
  { id: 'dem5b_png', z: 15, hsrc: '5m（写真測量）' },
  { id: 'dem5c_png', z: 15, hsrc: '5m（写真測量5C）' },
  { id: 'dem_png', z: 14, hsrc: '10m' },
];

// 緯度経度 → タイル番号と、タイル内のピクセル位置（Webメルカトル）
export function locate(lat, lon, z) {
  const n = 2 ** z;
  const rad = (lat * Math.PI) / 180;
  const fx = ((lon + 180) / 360) * n;
  const fy = ((1 - Math.log(Math.tan(rad) + 1 / Math.cos(rad)) / Math.PI) / 2) * n;
  const x = Math.floor(fx), y = Math.floor(fy);
  return { x, y, px: Math.floor((fx - x) * 256), py: Math.floor((fy - y) * 256) };
}

// 8bit RGB / RGBA のPNGだけを扱う最小のデコーダー
export function decodePng(buf) {
  if (buf.readUInt32BE(0) !== 0x89504e47) throw new Error('PNGではない');
  let pos = 8, width, height, colorType, idat = [];
  while (pos < buf.length) {
    const len = buf.readUInt32BE(pos);
    const type = buf.toString('ascii', pos + 4, pos + 8);
    const data = buf.subarray(pos + 8, pos + 8 + len);
    if (type === 'IHDR') {
      width = data.readUInt32BE(0); height = data.readUInt32BE(4);
      const depth = data[8]; colorType = data[9];
      if (depth !== 8 || ![2, 6].includes(colorType) || data[12] !== 0) throw new Error(`未対応のPNG: depth=${depth} color=${colorType} interlace=${data[12]}`);
    } else if (type === 'IDAT') idat.push(data);
    else if (type === 'IEND') break;
    pos += 12 + len;
  }
  const bpp = colorType === 6 ? 4 : 3;
  const raw = inflateSync(Buffer.concat(idat));
  const stride = width * bpp;
  const out = Buffer.alloc(stride * height);
  for (let y = 0; y < height; y++) {
    const filter = raw[y * (stride + 1)];
    const line = raw.subarray(y * (stride + 1) + 1, (y + 1) * (stride + 1));
    for (let i = 0; i < stride; i++) {
      const a = i >= bpp ? out[y * stride + i - bpp] : 0;           // 左
      const b = y > 0 ? out[(y - 1) * stride + i] : 0;              // 上
      const c = y > 0 && i >= bpp ? out[(y - 1) * stride + i - bpp] : 0; // 左上
      let v = line[i];
      if (filter === 1) v += a;
      else if (filter === 2) v += b;
      else if (filter === 3) v += (a + b) >> 1;
      else if (filter === 4) { const p = a + b - c, pa = Math.abs(p - a), pb = Math.abs(p - b), pc = Math.abs(p - c); v += pa <= pb && pa <= pc ? a : pb <= pc ? b : c; }
      out[y * stride + i] = v & 0xff;
    }
  }
  return { width, height, bpp, pixels: out };
}

// 標高タイルの詳細仕様の式: x = 2^16R + 2^8G + B、x = 2^23 は無効値、x > 2^23 は負の値
export function rgbToElevation(r, g, b) {
  const x = r * 65536 + g * 256 + b;
  if (x === 8388608) return null;
  return Math.round((x < 8388608 ? x : x - 16777216)) / 100;
}

const cache = new Map();
async function fetchTile(id, z, x, y) {
  const url = `https://cyberjapandata.gsi.go.jp/xyz/${id}/${z}/${x}/${y}.png`;
  if (cache.has(url)) return cache.get(url);
  const res = await fetch(url, { headers: { 'User-Agent': UA } });
  const result = res.ok ? { url, status: res.status, png: decodePng(Buffer.from(await res.arrayBuffer())) } : { url, status: res.status, png: null };
  cache.set(url, result);
  return result;
}

// 精度の高いタイルから順に見て、最初に値があったものを返す
export async function elevationFromTiles(lat, lon, { onFetch } = {}) {
  const tried = [];
  for (const s of SOURCES) {
    const { x, y, px, py } = locate(lat, lon, s.z);
    const t = await fetchTile(s.id, s.z, x, y);
    onFetch?.(t);
    if (!t.png) { tried.push({ source: s.id, status: t.status }); continue; }
    const i = (py * t.png.width + px) * t.png.bpp;
    const h = rgbToElevation(t.png.pixels[i], t.png.pixels[i + 1], t.png.pixels[i + 2]);
    tried.push({ source: s.id, status: t.status, value: h });
    if (h !== null) return { elevation: h, hsrc: s.hsrc, tile: t.url, pixel: [px, py], tried };
  }
  return { elevation: null, hsrc: null, tried };
}

- 標高の回答は get_elevation の結果だけを根拠にし、hsrc と出典を併記する
- elevation が null の地点は「データなし」と書き、近くの地点の値で補わない
- 緯度と経度の取り違えが疑われる入力（緯度が90を超えるなど）は、呼び出さずに確認を求める
- 浸水や避難の可否を標高だけで判断しない
- 出典は「地理院タイル（標高タイル）を加工して作成」または国土地理院の標高APIのURLで残す