Krzysztof Dyba

In `write.jl` NoData Value should be set to 0, not 65535. In `terra` the output is saved to the temporary file.

Why are the results in `raster.jl` different? ``` terra::extract(ras, buffers, fun = mean, exact = FALSE)[, -1] #> B1 B10 B11 B2 B3 B4 B5 B6 B7 B9 #> 1...

Hey @asinghvi17 and @rafaqz! I ran the scripts in `rasters_jl` and compared the output with `terra` and noticed a lot of differences. Could you please check this and respond?

Thanks! > For crop, what are you trying to measure? If its data loading time then let's call this something else maybe? I actually wrote "load", but what I mean...

Thanks! If you have any fixes, please let me know and I will run it again. I'm unsure about the output structure for point extraction and zonal statistics tasks, because...

From what I can see, converting to data frame has no impact: ``` julia> @be zonal($(Statistics.mean), $rstack; of=$(buffer_df.geom), progress=false) seconds = 10 Benchmark: 76 samples with 1 evaluation min 122.162...

Setting `missingval=0` in `write` doesn't work for me: ``` julia> @be write($stack_file, $ras; missingval=0, force=$true, options = $(Dict("COMPRESS" => "LZW", "INTERLEAVE" => "BAND"))) seconds=60 ERROR: GDALError (CE_Failure, code 1): GetNoDataValue()...

Ok, thanks! Maybe I should add information to the README that Julia functions are compiled and take longer to run first time. Later, during the benchmark, from what I noticed,...

Have you tried using the minimum or maximum instead of the mean in the zonal statistics? I get an error: ``` zonal((Statistics.minimum), rstack; of=(buffer_df.geom), progress=false) #> ERROR: TaskFailedException #> nested...

What is surprising is that the mean works, but the minimum/maximum doesn't. This problem may be partly related to the smaller extent of thermal bands (B10 and B11) compared to...