# Sample points
dataset = far.sample(nsamp=10000, adapt=True, seed=1234, csize=10,
                     var_dir="data",
                     input_forest_raster="forest2019.tif",
                     output_file="output/sample.txt",
                     blk_rows=0)


Sample 2x 10000 pixels (deforested vs. forest)
Divide region in 289 blocks
Compute number of deforested and forest pixels per block
100%
Draw blocks at random
Draw pixels at random in blocks
100%
Compute center of pixel coordinates
Compute number of 10 x 10 km spatial cells
... 1 cells (1 x 1)
Identify cell number from XY coordinates
Make virtual raster with variables as raster bands
Extract raster values for selected pixels
100%
Export results to file output/sample.txt
---------------------------------------------------------------------------
ValueError                                Traceback (most recent call last)
~\AppData\Local\Temp/ipykernel_19528/4146077674.py in <module>
      4                      input_forest_raster="forest2019.tif",
      5                      output_file="output/sample.txt",
----> 6                      blk_rows=0)

~\miniconda3\envs\conda-far\lib\site-packages\forestatrisk\build_data\sample.py in sample(nsamp, adapt, seed, csize, var_dir, input_forest_raster, output_file, blk_rows)
    309     # Convert to pandas DataFrame and return the result
    310     colname.extend(["X", "Y", "cell"])
--> 311     val_DF = pd.DataFrame(val, columns=colname)
    312     return val_DF
    313 

~\miniconda3\envs\conda-far\lib\site-packages\pandas\core\frame.py in __init__(self, data, index, columns, dtype, copy)
    676                     dtype=dtype,
    677                     copy=copy,
--> 678                     typ=manager,
    679                 )
    680 

~\miniconda3\envs\conda-far\lib\site-packages\pandas\core\internals\construction.py in ndarray_to_mgr(values, index, columns, dtype, copy, typ)
    322     )
    323 
--> 324     _check_values_indices_shape_match(values, index, columns)
    325 
    326     if typ == "array":

~\miniconda3\envs\conda-far\lib\site-packages\pandas\core\internals\construction.py in _check_values_indices_shape_match(values, index, columns)
    391         passed = values.shape
    392         implied = (len(index), len(columns))
--> 393         raise ValueError(f"Shape of passed values is {passed}, indices imply {implied}")
    394 
    395 

ValueError: Shape of passed values is (20000, 5), indices imply (20000, 7)