diff --git a/src/Functions/GeoUtils.h b/src/Functions/GeoUtils.h
index 3a5e020246916b1901a7cfd2a6737b43662cb423..d084347431404ad12093f71bbd965a2d0a04621c 100644
--- a/src/Functions/GeoUtils.h
+++ b/src/Functions/GeoUtils.h
@@ -301,6 +301,9 @@ bool PointInPolygonWithGrid<CoordinateType>::contains(CoordinateType x, Coordina
     if (has_empty_bound)
         return false;
 
+    if (std::isnan(x) || std::isnan(y))
+        return false;
+
     CoordinateType float_row = (y + y_shift) * y_scale;
     CoordinateType float_col = (x + x_shift) * x_scale;
 
diff --git a/tests/queries/0_stateless/00500_point_in_polygon_nan.reference b/tests/queries/0_stateless/00500_point_in_polygon_nan.reference
new file mode 100644
index 0000000000000000000000000000000000000000..573541ac9702dd3969c9bc859d2b91ec1f7e6e56
--- /dev/null
+++ b/tests/queries/0_stateless/00500_point_in_polygon_nan.reference
@@ -0,0 +1 @@
+0
diff --git a/tests/queries/0_stateless/00500_point_in_polygon_nan.sql b/tests/queries/0_stateless/00500_point_in_polygon_nan.sql
new file mode 100644
index 0000000000000000000000000000000000000000..37ed8dbeded8374869ea0b221aff902139d08a0f
--- /dev/null
+++ b/tests/queries/0_stateless/00500_point_in_polygon_nan.sql
@@ -0,0 +1 @@
+SELECT pointInPolygon((nan, 10.000100135803223), [(39.83154, 21.41527), (2., 1000.0001220703125), (39.90033, 21.37195), (1.000100016593933, 10.000100135803223), (39.83051, 21.42553), (39.82898, 21.41382), (39.83043, 21.41432), (39.83154, 21.41527)]);