def METHOD_NAME(lhaaso1):
    m = lhaaso1[0].spatial_model("KM2A")
    assert isinstance(m, PointSpatialModel)
    assert m.lon_0.unit == "deg"
    assert_allclose(m.lon_0.value, 1.86, atol=1e-2)
    assert_allclose(m.lon_0.error, 0.09, atol=1e-2)
    assert m.lat_0.unit == "deg"
    assert_allclose(m.lat_0.value, 57.00, atol=1e-2)
    assert_allclose(m.lat_0.error, 0.049, atol=1e-2)
    assert m.frame == "fk5"
    m = lhaaso1[11].spatial_model("WCDA")
    assert isinstance(m, GaussianSpatialModel)
    assert m.lon_0.unit == "deg"
    assert_allclose(m.lon_0.value, 75.01, atol=1e-10)
    assert m.lat_0.unit == "deg"
    assert_allclose(m.lat_0.value, 44.92, atol=1e-10)
    assert m.frame == "fk5"
    assert m.sigma.unit == "deg"
    assert_allclose(m.sigma.value, 0.41, atol=1e-3)
    model = lhaaso1["1LHAASO J0007+5659u"].spatial_model("KM2A")
    pos_err = model.position_error
    assert_allclose(pos_err.height.value, 2 * 0.12, rtol=1e-4)
    assert_allclose(pos_err.width.value, 2 * 0.12, rtol=1e-4)
    assert_allclose(model.position.ra.value, pos_err.center.ra.value)
    assert_allclose(model.position.dec.value, pos_err.center.dec.value)