CH-Earth / gistool

gistool: Static Data Processing tool for HPCs
GNU General Public License v3.0
2 stars 2 forks source link

soil class extraction does not generate output #44

Closed ShervanGharari closed 7 months ago

ShervanGharari commented 7 months ago

I am trying this:

rm -rf gistool
git clone https://github.com/kasra-keshavarz/gistool.git
cd gistool
./extract-gis.sh --dataset="soil_class" \
   --dataset-dir="/project/rpp-kshook/Model_Output/soil_classes/" \
   --variable="usda_mode_soilclass_250m_ll" \
   --start-date="" \
   --end-date="" \
   --shape-file="/project/6006250/CompHydCore/domain/domain_cat.shp" \
   --print-geotiff=true \
   --output-dir="$HOME/scratch/test/gistool_outputs/soil_classes" \
   --prefix="test_" \
   --stat="majority,minority,frac"

and will get this:

Attaching package: ‘dplyr’

The following objects are masked from ‘package:stats’:

    filter, lag

The following objects are masked from ‘package:base’:

    intersect, setdiff, setequal, union

Error in .local(.Object, ...) :

Error in .rasterObjectFromFile(x, band = band, objecttype = "RasterLayer",  :
  Cannot create a RasterLayer object from this file. (file does not exist)
Calls: <Anonymous> -> <Anonymous> -> .local -> .rasterObjectFromFile
Execution halted

shapefile is attached. domain_cat.zip

edit: updated the script so it can be read directly on Graham the input shapefile

kasra-keshavarz commented 7 months ago

Sorry for the confusion, the variable names are now all consistent. Previously, it was only soil_classes, but now your prompt should work fine too. I close this issue as I believe it is resolved.