SpatialObject box = BoxLatLon.newBox(lat.subtract(r).doubleValue(),
lat.add(r).doubleValue(),
lon.subtract(r).doubleValue(),
lon.add(r).doubleValue());
long[] zValues = new long[box.maxZ()];
space.decompose(box, zValues);
for (int i = 0; i < box.maxZ(); i++) {
long z = zValues[i];
if (z != -1L) {
ExpressionNode lo = new ConstantExpression(Space.zLo(z), InternalIndexTypes.LONG.instance(true));
ExpressionNode hi = new ConstantExpression(Space.zHi(z), InternalIndexTypes.LONG.instance(true));