from qgis.core import QgsProject, QgsRasterLayer

import processing

import os

raster_layer = QgsProject.instance().mapLayersByName("9.28.25.DEM")[0]

output_path = "C:/Users/bayan/Documents/L+U/EVLU 3004 - Ecology and Design 2/ReadingLandscapesInTime/GIS/9.28.25_MASK.tif"

result = processing.run("gdal:rastercalculator", {

'INPUT_A': raster_layer.source(),

'BAND_A': 1,

'FORMULA': "A > 246",

'OUTPUT': output_path,

'RTYPE': 5, # Float32

'NO_DATA': None,

'OPTIONS': '',

'EXTRA': ''

})

if os.path.exists(output_path):

binary_layer = QgsRasterLayer(output_path, "9.28.25_MASK", "gdal")

if binary_layer.isValid():

QgsProject.instance().addMapLayer(binary_layer)

print(" Binary mask created and loaded successfully.")

else:

print("Raster layer loaded but is invalid.")

else:

print("Output file not found. Check permissions or formula.")

BAYAN SHAERI