Latest Contributions by TungNg

Hi everyone, I had these errors while running Build Walls and DEM Reconditioning tools in either HEC-GeoHMS or ArcHydro  System info: ArcGIS version 10.6.1 HEC-GeoHMS 10.6: downloaded from here  Filled DEM (Fil) in C:/GIS/Layers Inner/Outer Wall layer (wall_test_buf.shp) in C:/GIS/Layers: is a buffer (30 feets) created from a line feature Errors: Executing: BuildWalls Fil 500 0 0 C:\GIS\Layers\WalledDEM # wall_test_buf # Start Time: Tue Jan 28 12:58:16 2020 Building internal walls Executing: CreateFeatureclass C:\Users\user007\AppData\Local\Temp\arc2C59\AH20200128095816\AHTemp.gdb tmpAD66 POLYLINE # DISABLED DISABLED "PROJCS['NAD_1983_UTM_Zone_10N',GEOGCS['GCS_North_American_1983',DATUM['D_North_American_1983',SPHEROID['GRS_1980',6378137.0,298.257222101]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]],PROJECTION['Transverse_Mercator'],PARAMETER['False_Easting',1640416.666666667],PARAMETER['False_Northing',0.0],PARAMETER['Central_Meridian',-123.0],PARAMETER['Scale_Factor',0.9996],PARAMETER['Latitude_Of_Origin',0.0],UNIT['Foot_US',0.3048006096012192]];-16800800 -32802000 137296494.950628;-100000 10000;-100000 10000;3.28083333333333E-03;0.001;0.001;IsHighPrecision" # 0 0 0 Start Time: Tue Jan 28 12:58:17 2020 Succeeded at Tue Jan 28 12:58:18 2020 (Elapsed Time: 0.48 seconds) Executing: PolylineToRaster GPL0 OBJECTID C:\Users\user007\AppData\Local\Temp\arc2C59\AH20200128095816\wtemp MAXIMUM_LENGTH NONE 10 Start Time: Tue Jan 28 12:58:18 2020 Succeeded at Tue Jan 28 12:58:18 2020 (Elapsed Time: 0.50 seconds) System.Runtime.InteropServices.COMException (0x80041098): ERROR 010429: Error in GRID IO: CellLyrCreateInternal: Unable to create cell layer directory c:\windows\system32\ras810 . ERROR 010302: Unable to create the output raster: C:\WINDOWS\SYSTEM32\RAS810 ERROR 010067: Error in executing grid expression. at ESRI.ArcGIS.SpatialAnalyst.RasterMathOpsClass.IsNull(IGeoDataset geoDataset) at ESRI.APWR.ApHydro.ArcHydroOp.BuildWalls(ApLayers apLayers, Double innerWallHeight, Double wallBuffer, Double breachLineBuffer, ITrackCancel& trackcancel, IGPMessages& messages, String& exMessage) Failed to execute (BuildWalls). Failed at Tue Jan 28 12:58:19 2020 (Elapsed Time: 2.11 seconds) Executing: DEMReconditioning Fil wall_test_dis 5 -500 -500 C:\GIS\Layers\AgreeDEM NEGATIVE_NO Start Time: Tue Jan 28 13:00:33 2020 Executing: PolylineToRaster C:\GIS\Layers\wall_test_dis.shp FID C:\Users\user007\AppData\Local\Temp\arc2C59\AH20200128100033\AgreeStrGeo MAXIMUM_LENGTH NONE C:\GIS\Layers\Fil Start Time: Tue Jan 28 13:00:33 2020 Succeeded at Tue Jan 28 13:00:33 2020 (Elapsed Time: 0.38 seconds) Executing: Con C:\Users\user007\AppData\Local\Temp\arc2C59\AH20200128100033\AgreeStrGeo C:\GIS\Layers\fil C:\Users\user007\AppData\Local\Temp\arc2C59\AH20200128100033\AgreeStrDEM # # Start Time: Tue Jan 28 13:00:34 2020 Succeeded at Tue Jan 28 13:00:34 2020 (Elapsed Time: 0.44 seconds) Executing: Delete C:\Users\user007\AppData\Local\Temp\arc2C59\AH20200128100033\AgreeStrGeo RasterDataset Start Time: Tue Jan 28 13:00:34 2020 Succeeded at Tue Jan 28 13:00:34 2020 (Elapsed Time: 0.05 seconds) Computing smooth drop/raise grid .. Executing: Minus C:\Users\user007\AppData\Local\Temp\arc2C59\AH20200128100033\AgreeStrDEM -500 C:\Users\user007\AppData\Local\Temp\arc2C59\AH20200128100033\smoothGeo Start Time: Tue Jan 28 13:00:35 2020 Succeeded at Tue Jan 28 13:00:35 2020 (Elapsed Time: 0.70 seconds) Computing vector distance and allocation grids... Executing: Int C:\Users\user007\AppData\Local\Temp\arc2C59\AH20200128100033\smoothGeo C:\Users\user007\AppData\Local\Temp\arc2C59\AH20200128100033\smoothGeoInt Start Time: Tue Jan 28 13:00:35 2020 Succeeded at Tue Jan 28 13:00:36 2020 (Elapsed Time: 0.36 seconds) Executing: Delete C:\Users\user007\AppData\Local\Temp\arc2C59\AH20200128100033\AgreeStrDEM RasterDataset Start Time: Tue Jan 28 13:00:36 2020 Succeeded at Tue Jan 28 13:00:36 2020 (Elapsed Time: 0.04 seconds) Executing: EucAllocation C:\Users\user007\AppData\Local\Temp\arc2C59\AH20200128100033\smoothGeoInt C:\Users\user007\AppData\Local\Temp\arc2C59\AH20200128100033\vecEucAlloGeo 45 # 10 Value C:\Users\user007\AppData\Local\Temp\arc2C59\AH20200128100033\vecEucDistGeo # PLANAR Start Time: Tue Jan 28 13:00:36 2020 ERROR 999999: Error executing function. Failed to execute (EucAllocation). Failed at Tue Jan 28 13:00:37 2020 (Elapsed Time: 0.58 seconds) System.Runtime.InteropServices.COMException (0x80004005): Error HRESULT E_FAIL has been returned from a call to a COM component. at ESRI.ArcGIS.Geoprocessing.GeoProcessorClass.Execute(String Name, IVariantArray ipValues, ITrackCancel pTrackCancel) at ESRI.APWR.ApHydro.ArcHydroOp.DEMReconditioningGP(ApLayers apLayers, Double numCellsBuff, Double sharpDrop, Double smoothDrop, Boolean raiseNegative, Boolean bUsingStreamGridInput, ITrackCancel& trackcancel, IGPMessages& messages, String& exMessage) Failed to execute (DEMReconditioning). Failed at Tue Jan 28 13:00:37 2020 (Elapsed Time: 3.65 seconds) Any clue about what went wrong?  Thanks! ... View more

Hello everyone, I am trying to update a map based on a pre-built template. However, after everything ran successfully (python code below), the shapefile that was updated did not show up on the map. Its symbology was also changed from "Quantities: Graduated colors" to "Features: Single symbol". I'm not sure what I did wrong? Any pointer is appreciated. Thanks! Steps that I did: Read the original shapefile in the template with preset symbology. Copied that shapefile to a new one. Added a csv file then joined it with the newly created shapefile. Updated existing field in the new feature with the values in the csv file. Change the data source from the original shapefile to the new one. Refreshed TOC and ActiveView then exported to a new mxd file. I opened the new mxd file. The data source is correctly updated to the new file. Values in the attribute table were also correct. import os import ntpath import arcpy from arcpy import env arcpy.env.overwriteOutput = True ### set project directory project_dir = os.path.expanduser('~/Projects/') mxd_dir = os.path.join(project_dir, 'GIS/mxd/') gis_dir = os.path.join(project_dir, 'GIS/shapefile') env.workspace = os.path.join(project_dir, gis_dir) year = '2019' dframe = 'Baseline' csv_file = 'sample.csv' ### open mxd template file mxd_file = "my_template.mxd" full_path = os.path.join(mxd_dir + mxd_file) mxd = arcpy.mapping.MapDocument(full_path) ### list data frame in the mxd file df = arcpy.mapping.ListDataFrames(mxd, dframe)[0] ### List layers layers = arcpy.mapping.ListLayers(df) ### This is the shapefile in the template with preset symbology inFeatures = os.path.join(gis_dir, "Q_" + dframe + ".shp") print(inFeatures) ### make a copy outFeatures = os.path.join(gis_dir, "output", "Q_" + dframe + "_" + year + ".shp") arcpy.CopyFeatures_management(in_features = inFeatures, out_feature_class = outFeatures) ### join then do field calculator on the new feature arcpy.MakeFeatureLayer_management(in_features = outFeatures, out_layer = 'new_lyr') arcpy.MakeTableView_management(in_table = csv_file, out_view = 'csv_ec') arcpy.AddJoin_management(in_layer_or_view = 'new_lyr', in_field = 'ID', join_table = 'csv_ec', join_field = 'ID') arcpy.CalculateField_management('new_lyr', "Q_" + dframe + "_" + year + ".Q", "!" + "Q_" + dframe + ".csv" + ".Q_new!", "PYTHON") ### update data source for lyr in layers: if dframe in lyr.name: lyr.replaceDataSource(os.path.join(gis_dir, "output", dframe), "SHAPEFILE_WORKSPACE", os.path.splitext(ntpath.basename(outFeatures))[0]) ### Update Layout View arcpy.RefreshTOC() arcpy.RefreshActiveView() ### save to new mxd file new_mxd = "new_Q_" + year.upper() + ".mxd" new_mxd_path = os.path.join(mxd_dir + new_mxd) mxd.saveACopy(new_mxd) ### Clean up the MapDocument object by deleting it del mxd‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍ ... View more

Happy Friday everyone, I have a multi-line shapefile which was clipped from NHDPlusFlowline feature class (picture below and attached) https://www.usgs.gov/core-science-systems/ngp/national-hydrography/nhdplus-high-resolution   The problem is there are so many tiny little lines together with the main canals. I want to find a good way to remove those small parts. I tried using Query Builder with “"LengthKM" >= 0.3” but it actually filtered out some small parts that belong to the main canals too   Any tip to get this going? Thanks! ... View more