ethoms-usgs / cross-section Goto Github PK

Hello!

I was really happy to find that somebody is working on making fence diagrams for ArcGIS, but I can't make it work, this message appears, and I don't undersatand what file is it about

p.s. I'm using ArcGIS 10.2.2 and I don't have other options

`Executing: fenceDiagram sectl_155 D:\three\demArcGis.tif northwest D:\three\p560506\main.gdb\baseb_pgn 1 D:\three\p560506\main.gdb
Start Time: Wed Nov 13 16:17:44 2019
Running script fenceDiagram...
Getting elevation values for cross-section line in sectl_155
sectl_155_z written to C:\Users\Georgiy_Savelyev\Documents\ArcGIS\Default.gdb
Measuring the length of the line in sectl_155_z
sectl_155_zm written to C:\Users\Georgiy_Savelyev\Documents\ArcGIS\Default.gdb
Converting D:\three\p560506\main.gdb\baseb_pgn to 3D features
Writing XYZ coordinates to D:\Current\Publications_Reviews\H\Haeussler_Kodiak\Kodiak_SIM\geologic_data\baseb_pgn_xyz.txt
File "D:\three\GeologicCrossSectionTbx\scripts\threeDsection.py", line 376, in
outF = open(xyzFile, 'w')

<type 'exceptions.IOError'>: [Errno 2] No such file or directory: u'D:\Current\Publications_Reviews\H\Haeussler_Kodiak\Kodiak_SIM\geologic_data\baseb_pgn_xyz.txt'
Failed script fenceDiagram...

Traceback (most recent call last):
File "D:\three\GeologicCrossSectionTbx\scripts\threeDsection.py", line 458, in
raise SystemError
SystemError

Failed to execute (fenceDiagram).
Failed at Wed Nov 13 16:17:52 2019 (Elapsed Time: 8,04 seconds)
`

Hi Evan-

Wondering if it's possible to tweak the "Map points to xs points" tool to allow for plotting of earthquake hypocenter depths on the cross section? Currently attempted to do this by inputting depth into the "Elevation" field with no luck.

Thanks

Ben

This is an amazing tool, can we convert it so it is compatible with ArcGIS Pro? I understand it is a side project but seems like it would be well worth the effort.

Ben

Hi Evan,

I am trying to make a cross-section from polyline SectionA in the demo cross-section layer you provided. I can't seem to get the parameters right and was wondering if you could give me some guidance as to where I am going wrong? I have attached the pictures of all 3 steps that I'm doing.

Thanks!

Brody
.

Andrea Croskrey helped me modify your "lines to surface profiles" to process multiple raster surfaces for a specific application of your code. Results as shown below:

code used:
LipanXsec.py.txt

Thanks for making this available
M.C. Robinson

I've been trying to create a surface profile for this path.

But after creating the surface profile, it results in a 2 meter in length section. It seem it is redistributing all the 25726 length to a 600x2m line. You can see in the red polygon how it is plotted. It seems to match the DEM, but not the lenght.

I'd appreciate if you could give me some light on this one. Thanks in advance.

Hi All,

First off, thank you for all of your great work on this.

More often than not, and apparently when I run any of the tools in the Toolbox, I get an ERROR 000210:

File "\psf\Home\Documents\OMSG\Active_Projects\Cross-Section-master\GeologicCrossSectionTbx\scripts\points2XsecView.py", line 251, in
arcpy.InterpolateShape_3d(dem, ptLayer, zPts)

<class 'arcgisscripting.ExecuteError'>: ERROR 000210: Cannot create output c:\ArcGIS10_DefaultGDB\Default.gdb\XS_str300-2_z
The table name is invalid.
The table was not found. [XS_str300-2_z]
Failed to execute (InterpolateShape).

Traceback (most recent call last):
File "\psf\Home\Documents\OMSG\Active_Projects\Cross-Section-master\GeologicCrossSectionTbx\scripts\points2XsecView.py", line 413, in
raise SystemError
SystemError

Failed to execute (points2XsecView).

Permissions are set for read/write on all files and folders, including said default gdb, and everything is in the same projection, and the source files are all in the same geodatabase.

Any tips?

Hallo
I am Huy, from Vietnam. I am doing a thesis on this topic and i wonder if you can help me to show me a bit ( with demo data) how the tool work?
Thank you!

Error running lineintersect.py on ArcMap 10.7.1

"line 304, in
arcpy.CreateRoutes_lr(xsecLayer, 'ORIG_FID', mLine, 'LENGTH', '#', '#', cp)

<class 'arcgisscripting.ExecuteError'>: ERROR 999999: Error executing function.
The table was not found. [A4_line.shp_m]
Could not create the output feature class
The table name is invalid.
The table was not found. [A4_line.shp_m]
Failed to execute (CreateRoutes)."

Plz step by step guide me to how can I install the all python scripts in arcgis 10 Toolbox
Because I add the scripts in newly created under GeologicCrossSectionTbx. After adding the scripts how can I define script parameters and run the scripts file plz guide me because I have not software background skill

The tool needs both a name for intermediate files and a pointer to the layer in the map. When a layer is chosen by the parameter input form, the pointer includes the name of any groups the layer may nested within, eg 'GEOLOGY\CrossSectionLines' which can't be used as an intermediate file name prefix. My first attempt to solve this was to get the Describe.name property of the layer object which works when the data source is a featureclass in a file gdb, but does not work when the source is a shapefile because the name comes back with a .shp extension.

Hi,

I am working on demo data provided along with the Cross Section master tool. During execution of Map points to Cross Section points tool, error is shown as in the screenshot attached. I am using ArcGIS 10.3 Advanced version. I don't know where I am making the mistake. Please help me out!!

Thanks and Regards,
Praveen

Hi!
Thanks for updating this tool and making it public!

I would like to say that I ran the "Lines to surface profiles" tool on a surface topography DEM and the "Borehole locations to borehole sticks tool" and they appear to have worked great.

The next thing I attempted was to add the surface profiles of buried formation tops and bottoms from rasters surfaces we interpolated from well log points. This may not be an intended use of the "Lines to surface profiles" tool. The lines look like they plotted along the Y axxis of the cross section correctly but not so much along the X axis. Some of the formations are not continuous and instead of spacing out the breaks in the formation contact line along the X axis the lines produced are abutted to each other and all start all the way to the left. As an example I have attached the same cross section three times. One image was created by a co-worker in ViewLog, one was created in the trial version of Strater, and then the third is only partially created using tools from this toolbox. Do you see in the first two examples how the green Cretaceous unit should be in two patches and the contacts are drawn right after the other in the third? I thought I could just edit this by making points along the topographic surface making the extend of the Cretaceous units and then sliding the lines over but I haven't been able to get the "Line intersections to cross-section features" tool to run.

Thoughts?
Andrea

Thanks for these tools ethoms-usgs! I am trying to run the "Map points to cross section points- Beta" tool and it is erroring with the message here. Wondering if there is a work around. Everything is in a GDB and in the same projection.
Thanks for any insight.

"
Traceback (most recent call last):
File "C:\Users...Tools\GIS_tools\Cross-Section-master\GeologicCrossSectionTbx\scripts\points2XsecView.py", line 212, in
arcpy.SelectLayerByLocation_management(ptLayer, 'WITHIN_A_DISTANCE', mLine, buff)
File "c:\arcgis\desktop10.7\arcpy\arcpy\management.py", line 7890, in SelectLayerByLocation
raise e
ExecuteError: Failed to execute. Parameters are not valid.
ERROR 000368: Invalid input data.
Failed to execute (SelectLayerByLocation).
"

Hello. Hope you're well
Tried using the cross section tool to created segmented surface profiles but after a success dialog nothing appears on my screen even if i zoom to layer. please how do i fix this?. Thank You