Dear experts, I ma trying to reproject my raster file using pyqgis script. The below script runs but there is no output in the folder created for the run.
##Settlements generation SOP
from qgis.core import *
from osgeo import gdal
import processing
import qgis.utils
import os # This is is needed in the pyqgis console also
#Reproject raster to Lat/Lon
def reproject_raster (in_fn,out_fn):
processing.run("gdalogr:warpreproject", {'INPUT':in_fn,'SOURCE_CRS':"EPSG:4326",'TARGET_CRS':'EPSG:32736','RESAMPLING':0,'NODATA':-9999,'TARGET_RESOLUTION':None,'OPTIONS':'','DATA_TYPE':0,'TARGET_EXTENT':None,'TARGET_EXTENT_CRS':None,'MULTITHREADING':False,'EXTRA':'','OUTPUT':out_fn})
in_fn = "/AQUAYA_PROJECTS/Trial_4_rwa/datasets/rwa/pop/rwa_ppp_2020_constrained.tif"
out_fn = "/AQUAYA_PROJECTS/Trial_4_rwa/datasets/rwa/comms/rwa_ppp_2020_con_repro.tif"
Please I need your help on correcting the script to obtain an output
Thank you
- Badges
- Science topic
- Similar topics
- Filing