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
|
443
|
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
|
443
|
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
|
1270
|
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
|
516
|
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
|
384
|
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
|
384
|
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
|
384
|
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
|
813
|
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
|
177
|
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
|
692
|
POST
|
Thanks for that. I'll give it a go and get back to you. PS I'm running VBA from Access
... View more
03-07-2013
03:00 PM
|
0
|
0
|
259
|
POST
|
Hi, The following code add two arguments to the Python call. The first argument works fine but the second gives an invalid SQL error when using SELECT. The script I have attached shows the parcel string hard coded to a Parcels variable, this works OK. Can someone provide the correct string for strParcels? Note the \\ is meant to relate to a \ in the final SQL strOwner = "Lett" strParcels = "PARCEL_SPI = '1\\TP8994'" retVal = Shell("cmd.exe /K S:\MID_Owners_Database_Workspace\OwnerProperty.py " & strOwner & " " & strParcels, 1)
... View more
03-07-2013
01:01 PM
|
0
|
3
|
557
|
POST
|
That worked fine. I thank you sincerely, but what does the deep imbedding in double quotes do? Can I also pass the string in as a Python argument i.e. how would I deal with all the quotes in a variable. Could I do something like ,""""&var1&''''''"). Once again many thanks
... View more
02-26-2013
02:54 PM
|
0
|
0
|
267
|
POST
|
Hi The following script is giving a syntax errror on the arcpy.Select line. I've tried various combinations of Quotes but I still get the invalid sql statement. Any ideas? # Import arcpy module import arcpy # Local variables: PropertyNew = "M:\\GIS\\KMZ Irrigation East\\PropertyNew.shp" OwnerPropertyShp = "M:\\GIS\\KMZ Irrigation East\\OwnerProperty.shp" # Process: Select arcpy.Select_analysis(PropertyNew, OwnerPropertyShp, "[PARCEL_SPI] = '2003\\PP2679'") Traceback (most recent call last): File "M:\GIS\KMZ Irrigation East\owner.py", line 19, in <module> arcpy.Select_analysis(PropertyNew, OwnerPropertyShp, "[PARCEL_SPI] = '2003\\PP2679'") File "C:\Program Files\ArcGIS\Desktop10.0\arcpy\arcpy\analysis.py", line 88, in Select raise e ExecuteError: ERROR 999999: Error executing function. An invalid SQL statement was used. An invalid SQL statement was used. Failed to execute (Select).
... View more
02-26-2013
11:00 AM
|
0
|
3
|
2196
|
Online Status |
Offline
|
Date Last Visited |
11-11-2020
02:23 AM
|