r.external: read nan as nan (OSGeo#338)

Fixes OSGeo#329. atof() on windows can not read "nan" and returns 0 (zero) instead. Using Rast_set_d_null_value() in this case to set nan directly.
petrasovaa · Feb 18, 2020 · c106583 · c106583
1 parent b979f6c
commit c106583
Showing 1 changed file with 4 additions and 1 deletion.
diff --git a/lib/raster/gdal.c b/lib/raster/gdal.c
@@ -286,8 +286,11 @@ struct GDAL_link *Rast_get_gdal_link(const char *name, const char *mapset)
     p = G_find_key_value("null", key_val);
     if (!p)
 	return NULL;
-    if (strcmp(p, "none") == 0)
+    /* atof on windows can not read "nan" and returns 0 instead */
+    if (strcmp(p, "none") == 0 ||
+        G_strcasecmp(p, "nan") == 0 || G_strcasecmp(p, "-nan") == 0) {
 	Rast_set_d_null_value(&null_val, 1);
+    }
     else
 	null_val = atof(p);