def FixStupidNoData():
DataDirectory = "T://analysis_for_papers//Indus//"
DEMName = "indus_utm44.bil"
newDEMName = "Indus_ND.bil"
RasterFilename = DataDirectory + DEMName
NewRasterFilename = DataDirectory + newDEMName
LSDP.SetNoDataBelowThreshold(RasterFilename,
NewRasterFilename,
threshold=0,
driver_name="ENVI",
NoDataValue=-9999)
def TestNewMappingTools2():
DataDirectory = "T://analysis_for_papers//Manny_idaho//"
Filename = "TestIdaho.bil"
NewFilename = "TestIdaho_after2.bil"
ThreshFname = "ThreshIdaho.bil"
ConstFname = "ConstEros.bil"
ThisFile = DataDirectory + Filename
NewFilename = DataDirectory + NewFilename
ThreshFname = DataDirectory + ThreshFname
ConstFname = DataDirectory + ConstFname
#FigFormat = 'svg'
#FigFileN= 'Sorbas_chi.svg'
#FigFileName= DataDirectory+FigFileN
#Plot the basin over the elevation
#tcmapcolorbarlabel = "Elevation (m)"
#tcmap = 'jet'
#tcmapcolorbarlabel='Chi'
#clim_val = (50,300)
#LSDP.BasicDensityPlotGridPlot(ThisFile,tcmap,tcmapcolorbarlabel,clim_val)
#get the nodata values
NData = LSDP.getNoDataValue(ThisFile)
print "NoData is: " + str(NData)
# get the data as an array
array = LSDP.ReadRasterArrayBlocks(ThisFile)
driver_name = "ENVI"
LSDP.array2raster(ThisFile, NewFilename, array, driver_name, -9999)
#get the nodata values
NData2 = LSDP.getNoDataValue(NewFilename)
#print "NoData 2 is: " + str(NData2)
LSDP.SetNoDataBelowThreshold(ThisFile, ThreshFname)
# now print the constant value file
constant_value = 0.001
LSDP.SetToConstantValue(ThreshFname, ConstFname, constant_value)
def FloodThenHillshade():
DataDirectory = "M:\\students\\kunkelova\\"
DEMName = "bk_10m_dem.bil"
newDEMName = "bk_10m_dem_updated.bil"
HillshadeName = "bk_10m_dem_HS.bil"
RasterFilename = DataDirectory + DEMName
NewRasterFilename = DataDirectory + newDEMName
HillshadeNameFile = DataDirectory + HillshadeName
LSDP.SetNoDataBelowThreshold(RasterFilename,
NewRasterFilename,
threshold=0,
driver_name="ENVI",
NoDataValue=-9999)
LSDP.GetHillshade(NewRasterFilename,
HillshadeNameFile,
azimuth=315,
angle_altitude=45,
driver_name="ENVI",
NoDataValue=-9999)
