grass_script.txt

# Currently submitted line by line
r.in.gdal input=/home/ubuntu/data/Monarques/CAN_NALCMS_2015_LC_30m_LAEA_mmu5pix_.tif output=CEC memory=150000 --overwrite
g.region raster=CEC
v.in.ogr input='/home/ubuntu/data/Monarques/canada_hwy.shp' output=Roads
v.in.ogr input='/home/ubuntu/data/Monarques/canada_rails_cec.shp' output=Rails
#v.in.ogr input='/home/ubuntu/data/Monarques/transmission_lines_cec.shp' output=Transmission
v.in.ogr input='/home/ubuntu/data/Monarques/transmission_lines_OSM_CEC.gpkg' output=Transmission_OSM
v.in.ogr input='/home/ubuntu/data/Monarques/protected_areas_canada_qc.shp' output=Protected snap=1e-09
v.to.rast input=Roads output=Roads use=val value=20 memory=150000 --overwrite
v.to.rast input=Rails output=Rails use=val value=21 memory=150000 --overwrite
v.to.rast input=Transmission_OSM output=Transmission use=val value=22 memory=150000 --overwrite
v.to.rast input=Protected output=Protected use=val value=100 memory=150000 
r.patch input=Roads,Rails,Transmission,CEC output=CEC2 --overwrite
r.null null=0 map=CEC2
r.null null=0 map=Protected
r.mapcalc expression='CEC_pa = CEC2+Protected' --overwrite
r.out.gdal input=CEC_pa output=/home/ubuntu/data/Monarques/cec_pa_lines.tif format=GTiff -createopt='COMPRESS=DEFLATE' --overwrite
r.reclass input=CEC_pa output=CEC_suitable rules='/home/ubuntu/data/Monarques/cec_reclass.txt' --overwrite

v.in.ogr input='/home/ubuntu/data/Monarques/canada_mastersample_in_range_poly_cec.shp' output=gtrs_10 --overwrite
v.to.rast input=gtrs_10 output=gtrs_10 use=attr attribute_column=GTRS_ID memory=150000 --overwrite
r.stats.zonal base=gtrs_10 cover=CEC_suitable method=sum output=gtrs_suitable_zonal --overwrite

r.out.gdal input=gtrs_suitable_zonal output=/home/ubuntu/data/Monarques/gtrs_suitable_zonal.tif format=GTiff createopt='COMPRESS=DEFLATE'

r.to.vect input=gtrs_suitable_zonal output=gtrs_suitable_zonal type=area column=suitable
v.out.ogr input=gtrs_suitable_zonal type=area output=/home/ubuntu/data/Monarques/gtrs_suitable_zonal.gpkg