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)