Fix inverted land mask

AlexeyPechnikov · Feb 26, 2024 · 350327f · 350327f
1 parent 7281edf
commit 350327f
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/pygmtsar/pygmtsar/GMT.py b/pygmtsar/pygmtsar/GMT.py
@@ -230,7 +230,7 @@ def download_landmask(self, geometry, filename=None, product='1s', resolution='f
         def grdlandmask(bounds, product, resolution, filename):
             if os.path.exists(filename):
                 os.remove(filename)
-            argv = ['gmt', 'grdlandmask', f'-R{bounds[0]}/{bounds[2]}/{bounds[1]}/{bounds[3]}', f'-I{product}', f'-D{resolution}', '-N1/NaN', f'-G{filename}']
+            argv = ['gmt', 'grdlandmask', f'-R{bounds[0]}/{bounds[2]}/{bounds[1]}/{bounds[3]}', f'-I{product}', f'-D{resolution}', '-NNaN/1', f'-G{filename}']
             #print ('grdlandmask argv:', ' '.join(argv))
             p = subprocess.Popen(argv, stdout=subprocess.PIPE, stderr=subprocess.PIPE, encoding='utf8')
             stdout_data, stderr_data = p.communicate()