A dump of some helpful materials
The Shapely buffer operation wasn’t working until I installed it using pip but omitting the binaries. We were also getting an assertion failure during intersect opterations. Again, omitting binaries solved this problem.
See this post
pip uninstall shapely pip install --no-binary :all: shapely
Convert ESRI GRID to GeoTIIF
Use this command:
gdalwarp -t_srs 'EPSG:26911' ./US_140EVT/us_140evt/ ./landfire.tif
Clip National Landfire to a HUC8
gdalwarp us_140evt ./test.tif -cutline ../../17040209/Shape/WBDHU8.shp -crop_to_cutline -t_srs "EPSG:26911"
Installed using terminal using the command:
pip install pygeoprocessing --no-deps --no-binary :all:
Also needed to install cython using PyCharm package manager and spatial indexing package using as per this post:
brew install spatialindex
Convert national landfire raster to compressed GeoTiff
gdal_translate -ot Int16 -of GTiff -co COMPRESS=DEFLATE ./US_140EVT_20180618/Grid/us_140evt ./test.tif
Download Jordan’s old VBET projects
aws s3 sync s3://riverscapesdata . –exclude=”” –include=”/VBET/*” –dryrun
[BLM Surface Management Agency] geodatabase was downloaded and then the “surface management layer exported to shapefile (and reprpojected) using the following command:
ogr2ogr -f "ESRI Shapefile" ~/GISData/NationalProject/ownership/surface_management_agency.shp ~/GISData/BLM_Owneership/BLM_National_Surface_Management_Agency/sma_wm.gdb SurfaceManagementAgency -t_srs 'EPSG:4269' -select 'ADMIN_AGENCY_CODE'
Running nested scripts
Use the following syntax to run scripts that are in subfolders, but run them from the parent folder in the context of that parent folder. This makes it possible for nested python scritps to refer to other nested libraries and still find these libraries when running the script directly.
lib below is the folder name but it is separated by a period and not a directory separator. And there is no
.py file suffix needed.
python -m lib.build_vrt arg1 arg2 arg3
- Downloaded the Landfire 2 REMAP dataset as multiple, overlapping tiles. Make sure to change the raster format to GTiff before commencing any of the downloads.
- Used the pyBRAT4 tool to build a single VRT file that references each GTiff.
gdal_translateto convert the VRT into a single compressed GTiff.
gdal_warpto reproject the raster.
gdal_translateto compress the reprojcted raster (because
gdal_warpproduces bloated rasters).
gdal_translate -ot Int16 -of GTiff -co COMPRESS=DEFLATE ./landfire.vrt ./landfire_vrt.tif gdalwarp -t_srs 'EPSG:4326' ./landfire_vrt_comp.tif ./landfire_vrt_warp.tif gdal_translate -ot Int16 -of GTiff -co COMPRESS=DEFLATE ./landfire_vrt_warp.if ./landfire_vrt_warp_trans.tif
Warehouse RS Context download
rscli download ~/GISData/brat/rs_context/17060204 --type RSContext --meta huc8=17060204 --verbose rscli download ~/GISData/BRAT/brat/17040201 --type BRAT --meta huc8=17040201 --verbose --file-filter "Outputs" --no-delete
Update the CLI
npm install -g @riverscapes/cli
NO_UI is an environment variable that, if set, will disable all colours and progress bars. This is mostly to create cleaner log files for headless runs.