Base flow and Aggregate results are coming with NoData

https://drive.google.com/drive/folders/1Wroj4pKzeVsfamXgD37wrpi81nMLhZ_7?usp=sharing

Hi @clara , thanks again for sharing your data. It appears that while you did assign -9999 as the NoData value, you did not actually reclassify the nan pixels to -9999. So the nan values remain. Here is a count of the number of pixels with a value of nan for each raster:

climando.tif
0
dem_sink.tif
0
ET0_2000_nodata10.tif
85527
ET0_2000_nodata11.tif
85527
ET0_2000_nodata12.tif
85527
ET0_2000_nodata9.tif
85527
ET0_2000_nodata_1.tif
85527
ET0_2000_nodata_2.tif
85527
ET0_2000_nodata_3.tif
85527
ET0_2000_nodata_4.tif
85527
ET0_2000_nodata_5.tif
85527
ET0_2000_nodata_6.tif
85527
ET0_2000_nodata_7.tif
85527
ET0_2000_nodata_8.tif
85527
lulc2000.tif
0
PR_2000_nodata10.tif
85527
PR_2000_nodata11.tif
85527
PR_2000_nodata12.tif
85527
PR_2000_nodata7.tif
85527
PR_2000_nodata8.tif
85527
PR_2000_nodata9.tif
85527
PR_2000_nodata_1.tif
85527
PR_2000_nodata_2.tif
85527
PR_2000_nodata_3.tif
85527
PR_2000_nodata_4.tif
85527
PR_2000_nodata_5.tif
85527
PR_2000_nodata_6.tif
85527
soilgroupnodata.tif
0

For reference, I generated this list using this script: Count the number of `nan` pixels in rasters in a directory · GitHub

1 Like

@dave Thank you!

This time somehow worked, the baseflow and local recharge didnt show the Nodata issue, nethertheless the max values of Baseflow and Local recharge are a bit out of hand, with baseflow raster with the value max of 18758158, in contrast with max baseflow of 2023 that has the value of 51,6576.

The qb presents for this scnerio of 2000 also high values, like 84.26 etc

could you take a look, because im pretty sure that this is not correct.

As you said, i set the nan values in et0 and pr as -9999 through Arcgis pro using a code. And the other rasters i set also a Nodaata value as i said before.

I am adding here the link again to this new model i ran. I would be very glad if you could check looking for eventual mistakes.

https://drive.google.com/drive/folders/1Wroj4pKzeVsfamXgD37wrpi81nMLhZ_7?usp=drive_link

I used this phyton code in Arcgispro to change the nan values to -9999:

# Função para processar os rasters de uma pasta
def process_rasters(input_folder):
    for filename in os.listdir(input_folder):
        if filename.endswith(".tif"):  # Verifica se o arquivo é um raster .tif
            input_raster = os.path.join(input_folder, filename)
            output_raster = os.path.join(output_folder, f"reclassified_{filename}")

            # Reclassificar valores NaN para -9999
            print(f"Processando: {input_raster}")
            try:
                raster = Raster(input_raster)
                reclass_raster = Con(IsNull(raster), -9999, raster)  # Substituir NaN por -9999
                reclass_raster.save(output_raster)
                print(f"Salvo: {output_raster}")
            except Exception as e:
                print(f"Erro ao processar {input_raster}: {e}")

I’m glad you resolved the missing data issue. I don’t have the expertise to say whether or not these values are reasonable. Perhaps @swolny could say?

This topic was automatically closed 7 days after the last reply. New replies are no longer allowed.