POST
|
Along each transect, line each cell value, will increase/decrease linearly using the two end point elevations.. In between transect lines, each cell value, will increase/decrease linearly using the corresponding value in the adjacent line. The shape of the resulting raster, in 3D, would resemble a ribbon with no folds or bends across it, but wavy in appearance.
... View more
10-05-2016
08:38 PM
|
0
|
1
|
16
|
POST
|
Hi Dan, I'm new at this raster creation thingy but wouldn't creating a polygon and than a raster 50cm below would mirror the same values in the DEM, just 50 cm lower?
... View more
10-05-2016
08:02 PM
|
0
|
3
|
16
|
POST
|
Hi, Arcmap raster creation for cut/fill operation. We have a irrigation channel that we plan to pipe. We need to use the cut/fill tool to show us the volumes of soil needed to be moved. We have a DEM of the area. I need to create another raster to use with the current DEM. The new raster will follow the natural surface along the pipeline. i.e. The surface after the pipe is laid and the channel filled in. The attached image shows cross section lines every 50m down the pipeline. The elevation at the ends of each cross section will somehow be used in the new raster. The new raster elevations along each cross section will be interpolated in a straight line between the end points. How can a raster be produced that mimics a natural surface along the pipeline? Cheers Bob
... View more
10-05-2016
02:47 PM
|
0
|
5
|
843
|
POST
|
Hi, I have two adjacent DEMs with different projections (see attached) I suspect the difference will be minimal but has anyone have an idea to the scale of the difference? i.e. Should I be concerned?
... View more
06-25-2014
04:41 PM
|
0
|
1
|
313
|
POST
|
Hi, I placed the appropiate prj file in the workspace. I also placed the xyz data files directly in the workspace. It then was asked for an average_point_spacing. I used 2.0 and it all worked. The figures I had to set the point_spacing; Avg. Point Density: 0.5 pts/m Footprint Size: 0.32m I used the first one to assume the value of 2.0m, Am I correct? Cheers and many thanks '''**************************************************************************** Name: ASCII3D_to_Feature_Class Example Description: This script demonstrates how to use the ASCII3D_to_Feature_Class tool to create a point feature class from a set of text files in the specified workspace. ****************************************************************************''' # Import system modules` import arcpy, os from arcpy import env import exceptions, sys, traceback try: #prjFile = os.path.join(arcpy.GetInstallInfo()["InstallDir"], # "M:/GIS/ArcGIS 10 Coordinate System Files and Projections/ArcGIS 10 Coordinate System Files and Projections/Coordinate Systems/Projected Coordinate Systems/National Grids/Australia/GDA 1994 MGA Zone 55.prj") #CS = arcpy.SpatialReference(prjFile) # Obtain a license for the ArcGIS 3D Analyst extension arcpy.CheckOutExtension("3D") # Set environment settings env.workspace = "M:/GIS/data" # Set Local Variables inFormat = "XYZ" # Name of the output file outFC = "90MileBeachPts_from_ASCII.shp" # Geometry of the output feature class outType = "MULTIPOINT" zFactor = 1 # Coordinate system of the output feature class CS = "GDA 1994 MGA Zone 55.prj" fileSuffix = "mga55.xyz" decSep = "DECIMAL_POINT" # Specifies the decimal delimeter # Create list of ASCII files txtList = arcpy.ListFiles("*" + fileSuffix) # Verify the presence of TINs in the list if len(txtList) > 0: # Execute ASCII3D_to_Feature_Class arcpy.ASCII3DToFeatureClass_3d(txtList, inFormat, outFC, outType, zFactor,CS, 2.0, fileSuffix, decSep) except arcpy.ExecuteError: print arcpy.GetMessages() except: # Get the traceback object tb = sys.exc_info()[2] tbinfo = traceback.format_tb(tb)[0] # Concatenate error information into message string pymsg = 'PYTHON ERRORS:\nTraceback info:\n{0}\nError Info:\n{1}'\ .format(tbinfo, str(sys.exc_info()[1])) msgs = 'ArcPy ERRORS:\n {0}\n'.format(arcpy.GetMessages(2)) # Return python error messages for script tool or Python Window arcpy.AddError(pymsg) arcpy.AddError(msgs)
... View more
06-23-2014
09:46 PM
|
0
|
0
|
5
|
POST
|
Hi First of all I'm using 10.0. When I run the code at the bottom I get the following result; PythonWin 2.6.5 (r265:79096, Mar 19 2010, 21:48:26) [MSC v.1500 32 bit (Intel)] on win32. Portions Copyright 1994-2008 Mark Hammond - see 'Help/About PythonWin' for further copyright information. Traceback (most recent call last): File "C:\Python26\Lib\site-packages\Pythonwin\pywin\framework\scriptutils.py", line 326, in RunScript exec codeObject in __main__.__dict__ File "M:\GIS\LidarProd\XYZToFeature\XYZ to Feature Missing 90Ml Beach 20062014\XYZ_3DfeatureclassArcPyVers.py", line 8, in <module> import arcpy, os File "C:\Program Files\ArcGIS\Desktop10.0\arcpy\arcpy\__init__.py", line 17, in <module> from geoprocessing import gp File "C:\Program Files\ArcGIS\Desktop10.0\arcpy\arcpy\geoprocessing\__init__.py", line 14, in <module> from _base import * File "C:\Program Files\ArcGIS\Desktop10.0\arcpy\arcpy\geoprocessing\_base.py", line 568, in <module> env = GPEnvironments(gp) File "C:\Program Files\ArcGIS\Desktop10.0\arcpy\arcpy\geoprocessing\_base.py", line 565, in GPEnvironments return GPEnvironment(geoprocessor) File "C:\Program Files\ArcGIS\Desktop10.0\arcpy\arcpy\geoprocessing\_base.py", line 521, in __init__ self._refresh() File "C:\Program Files\ArcGIS\Desktop10.0\arcpy\arcpy\geoprocessing\_base.py", line 523, in _refresh envset = (set(env for env in self._gp.listEnvironments())) RuntimeError: NotInitialized >>> What isn't initialised? I also noticed that arcpy.SpatialReference() didn't have SpatialReference appear in the list that is normally there after typing the dot (Whatever that is called) Cheers '''**************************************************************************** Name: ASCII3D_to_Feature_Class Example Description: This script demonstrates how to use the ASCII3D_to_Feature_Class tool to create a point feature class from a set of text files in the specified workspace. ****************************************************************************''' # Import system modules` import arcpy, os from arcpy import env import exceptions, sys, traceback try: prjFile = os.path.join(arcpy.GetInstallInfo()["InstallDir"], "M:/GIS/ArcGIS 10 Coordinate System Files and Projections/ArcGIS 10 Coordinate System Files and Projections/Coordinate Systems/Projected Coordinate Systems/National Grids/Australia/GDA 1994 MGA Zone 55.prj") CS = arcpy.SpatialReference(prjFile) # Obtain a license for the ArcGIS 3D Analyst extension arcpy.CheckOutExtension("3D") # Set environment settings env.workspace = "M:/GIS/LidarProd/XYZToFeature/Missing Lidar from DSE 21022013/elevation/cep2007-08_srw-90mb/regular_surface/dem_ground/xyz_ascii" # Set Local Variables inFormat = "XYZ" # Name of the output file outFC = "90MileBeachPts_from_ASCII.shp" # Geometry of the output feature class outType = "MULTIPOINT" zFactor = 1 # Coordinate system of the output feature class #CS = "Geographic Coordinate Systems/GDA 1994 MGA Zone 55.prj" fileSuffix = "mga55.xyz" decSep = "DECIMAL_POINT" # Specifies the decimal delimeter # Create list of ASCII files txtList = arcpy.ListFiles("*" + fileSuffix) # Verify the presence of TINs in the list if len(txtList) > 0: # Execute ASCII3D_to_Feature_Class arcpy.ASCII3DToFeatureClass_3d(txtList, inFormat, outFC, outType, zFactor, CS, fileSuffix, decSep) except arcpy.ExecuteError: print arcpy.GetMessages() except: # Get the traceback object tb = sys.exc_info()[2] tbinfo = traceback.format_tb(tb)[0] # Concatenate error information into message string pymsg = 'PYTHON ERRORS:\nTraceback info:\n{0}\nError Info:\n{1}'\ .format(tbinfo, str(sys.exc_info()[1])) msgs = 'ArcPy ERRORS:\n {0}\n'.format(arcpy.GetMessages(2)) # Return python error messages for script tool or Python Window arcpy.AddError(pymsg) arcpy.AddError(msgs)
... View more
06-23-2014
05:10 PM
|
0
|
0
|
5
|
POST
|
I tried it with and without .prj but alas no success. The sample code states that it is for Output but I beleive the syntax is for input. I'm confused '''**************************************************************************** Name: ASCII3D_to_Feature_Class Example Description: This script demonstrates how to use the ASCII3D_to_Feature_Class tool to create a point feature class from a set of text files in the specified workspace. ****************************************************************************''' # Import system modules` import arcpy from arcpy import env import exceptions, sys, traceback try: # Obtain a license for the ArcGIS 3D Analyst extension arcpy.CheckOutExtension("3D") # Set environment settings env.workspace = "M:/GIS/LidarProd/XYZToFeature/Missing Lidar from DSE 21022013/elevation/cep2007-08_srw-90mb/regular_surface/dem_ground/xyz_ascii" # Set Local Variables inFormat = "XYZ" # Name of the output file outFC = "90MileBeachPts_from_ASCII.shp" # Geometry of the output feature class outType = "MULTIPOINT" zFactor = 1 # Coordinate system of the output feature class CS = "Geographic Coordinate Systems/GDA 1994 MGA Zone 55.prj" fileSuffix = "mga55.xyz" decSep = "DECIMAL_POINT" # Specifies the decimal delimeter # Create list of ASCII files txtList = arcpy.ListFiles("*" + fileSuffix) # Verify the presence of TINs in the list if len(txtList) > 0: # Execute ASCII3D_to_Feature_Class arcpy.ASCII3DToFeatureClass_3d(txtList, inFormat, outFC, outType, zFactor, CS, fileSuffix, decSep) except arcpy.ExecuteError: print arcpy.GetMessages() except: # Get the traceback object tb = sys.exc_info()[2] tbinfo = traceback.format_tb(tb)[0] # Concatenate error information into message string pymsg = 'PYTHON ERRORS:\nTraceback info:\n{0}\nError Info:\n{1}'\ .format(tbinfo, str(sys.exc_info()[1])) msgs = 'ArcPy ERRORS:\n {0}\n'.format(arcpy.GetMessages(2)) # Return python error messages for script tool or Python Window arcpy.AddError(pymsg) arcpy.AddError(msgs)
... View more
06-22-2014
07:38 PM
|
0
|
0
|
5
|
POST
|
Hi, I'm trying to run; ASCII3D_to_Feature_Class The geographic coordinate system I require is 'GDA 1994 MGA Zone 55' I've tried CS = "Coordinate Systems/Geographic Coordinate Systems/World/GDA 1994 MGA Zone 55.prj" But, I get the error; >>> ERROR 000622: Failed to execute (ASCII 3D to Feature Class). Parameters are not valid. ERROR 000628: Cannot set input into parameter input_coordinate_system. Some help with the correct syntax would be greatly appreciated Cheers
... View more
06-22-2014
04:03 PM
|
0
|
5
|
434
|
POST
|
I have a way around it. The reason I needed to do a spatial join was to be able to do a SUM over the point feature values in each parcel. The stat was done on a table of the joined data. I then rejoined the table to the already joined data with the multiple records for each Target_Fid. What I ended up doing is joining to the original parcel layer and hence only one record per parcel. Silly me. Thanks to all viewers!
... View more
04-28-2014
04:34 PM
|
0
|
0
|
2
|
POST
|
Hi, I have done a spacial join with a land parcel layer as the target. The join is one to many using point features. The output table consequently can have multiple records with the one Target_Fid. For example, if there are three point features within a land parcel, three output records are produced. When I put 'Labels on' for the output I can have multiple values appearing. How can I limit the label for a particular parcel? The label field will always have the same value for each output record for the one parcel.
... View more
04-28-2014
03:05 PM
|
0
|
1
|
517
|
Online Status |
Offline
|
Date Last Visited |
11-11-2020
02:23 AM
|