nav[aria-label="Primary Navigation"] { padding: 0; & ul { list-style: none; width: 100%; display: flex; flex-direction: row; justify-content: start; align-items: start; gap: 30px; padding: 0; & li { margin: 0; } & ul li { list-style: none; } } }

Find 3d distance between 2 sets of points

SB5

Hi,

I'd like to find the 3d distance between all points in two point feature classes. One is 12 000 building points, the other 22 000 point locations of a construction machine (along a line), constructing many km of tunnel. My aim is to be able to select the building points within x meters when the construction machine is at a specific location. How can I do this without running Near 3d several thousand times?

I am using arcgis 10.3 with basic license. I am not very familiar with scripting.

Find more posts tagged with

Accepted answers

XanderBakker

As an alternative I generated a script that will calculate the 3D distance between two sets of 3D point featureclasses.

The result is a table that contains the OID of the buildings, the OID of the tunnel points and the 3D distance:

By defining two relates in ArcMap you will be able to do what you are looking for.

Relate 1: based on the tunnel points, relating to the new table.

Relate 2: based on the table relating to the building points:

Make sure to define some attribute indexes in the table to increase the performance:

How does it work. In ArcMap you select a tunnel point and open the attribute table of the tunnel point featureclass:

Click on the relate button and select the relate:

This will open the table with the related combinations selected:

Then select the relate to propagate the selection to the building points:

... which will open the selected building points:

The map will show the selection:

The script I used for this is listed below. In my case I only had 6.5M combinations and it took 2 minutes to process the combinations on an average computer. This will probably take all night processing the 264M combinations.

Script:

<SPAN class="keyword token">import</SPAN> math
<SPAN class="keyword token">import</SPAN> arcpy
<SPAN class="keyword token">from</SPAN> time <SPAN class="keyword token">import</SPAN> localtime<SPAN class="punctuation token">,</SPAN> strftime

<SPAN class="keyword token">def</SPAN> <SPAN class="token function">main</SPAN><SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">import</SPAN> os
&nbsp;&nbsp;&nbsp; ShowProgress<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"Start process..."</SPAN><SPAN class="punctuation token">)</SPAN>

&nbsp;&nbsp;&nbsp; fc_buildings <SPAN class="operator token">=</SPAN> r<SPAN class="string token">'D:\Xander\GeoNet\Generate3DNearTable\gdb\data.gdb\building_points_3D'</SPAN>
&nbsp;&nbsp;&nbsp; fc_tunnel_points <SPAN class="operator token">=</SPAN> r<SPAN class="string token">'D:\Xander\GeoNet\Generate3DNearTable\gdb\data.gdb\tunnel_points_3D'</SPAN>

&nbsp;&nbsp;&nbsp; tbl_out <SPAN class="operator token">=</SPAN> r<SPAN class="string token">'D:\Xander\GeoNet\Generate3DNearTable\gdb\data.gdb\Near3D_table_v01'</SPAN>
&nbsp;&nbsp;&nbsp; fld_building_id <SPAN class="operator token">=</SPAN> <SPAN class="string token">'OID_Building'</SPAN>
&nbsp;&nbsp;&nbsp; fld_tunnel_point_id <SPAN class="operator token">=</SPAN> <SPAN class="string token">'OID_TunnelPoint'</SPAN>
&nbsp;&nbsp;&nbsp; fld_dist3D <SPAN class="operator token">=</SPAN> <SPAN class="string token">'Distance3D'</SPAN>

&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># max distance</SPAN>
&nbsp;&nbsp;&nbsp; max_dist <SPAN class="operator token">=</SPAN> <SPAN class="number token">250</SPAN>

&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># sr = arcpy.Describe(fc_buildings).spatialReference</SPAN>

&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># create lists of buildings and tunnel points</SPAN>
&nbsp;&nbsp;&nbsp; ShowProgress<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"Create lists of buildings and tunnel points..."</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; flds <SPAN class="operator token">=</SPAN> <SPAN class="punctuation token">(</SPAN><SPAN class="string token">'OID@'</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">'SHAPE@'</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; lst_building <SPAN class="operator token">=</SPAN> <SPAN class="punctuation token">[</SPAN><SPAN class="punctuation token">[</SPAN>r<SPAN class="punctuation token">[</SPAN><SPAN class="number token">0</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">,</SPAN> r<SPAN class="punctuation token">[</SPAN><SPAN class="number token">1</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">]</SPAN> <SPAN class="keyword token">for</SPAN> r <SPAN class="keyword token">in</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>da<SPAN class="punctuation token">.</SPAN>SearchCursor<SPAN class="punctuation token">(</SPAN>fc_buildings<SPAN class="punctuation token">,</SPAN> flds<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">]</SPAN>
&nbsp;&nbsp;&nbsp; lst_tunnelpoints <SPAN class="operator token">=</SPAN> <SPAN class="punctuation token">[</SPAN><SPAN class="punctuation token">[</SPAN>r<SPAN class="punctuation token">[</SPAN><SPAN class="number token">0</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">,</SPAN> r<SPAN class="punctuation token">[</SPAN><SPAN class="number token">1</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">]</SPAN> <SPAN class="keyword token">for</SPAN> r <SPAN class="keyword token">in</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>da<SPAN class="punctuation token">.</SPAN>SearchCursor<SPAN class="punctuation token">(</SPAN>fc_tunnel_points<SPAN class="punctuation token">,</SPAN> flds<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">]</SPAN>

&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># create output table for relate</SPAN>
&nbsp;&nbsp;&nbsp; ShowProgress<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"Create output table..."</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; tbl_ws<SPAN class="punctuation token">,</SPAN> tbl_name <SPAN class="operator token">=</SPAN> os<SPAN class="punctuation token">.</SPAN>path<SPAN class="punctuation token">.</SPAN>split<SPAN class="punctuation token">(</SPAN>tbl_out<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; arcpy<SPAN class="punctuation token">.</SPAN>CreateTable_management<SPAN class="punctuation token">(</SPAN>tbl_ws<SPAN class="punctuation token">,</SPAN> tbl_name<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; ShowProgress<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"Add fields to output table..."</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; AddField<SPAN class="punctuation token">(</SPAN>tbl_out<SPAN class="punctuation token">,</SPAN> fld_building_id<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"LONG"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="number token">8</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; AddField<SPAN class="punctuation token">(</SPAN>tbl_out<SPAN class="punctuation token">,</SPAN> fld_tunnel_point_id<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"LONG"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="number token">8</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; AddField<SPAN class="punctuation token">(</SPAN>tbl_out<SPAN class="punctuation token">,</SPAN> fld_dist3D<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">"DOUBLE"</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="number token">8</SPAN><SPAN class="punctuation token">)</SPAN>


&nbsp;&nbsp;&nbsp; cnt_total <SPAN class="operator token">=</SPAN> len<SPAN class="punctuation token">(</SPAN>lst_building<SPAN class="punctuation token">)</SPAN> <SPAN class="operator token">*</SPAN> len<SPAN class="punctuation token">(</SPAN>lst_tunnelpoints<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; flds <SPAN class="operator token">=</SPAN> <SPAN class="punctuation token">(</SPAN>fld_building_id<SPAN class="punctuation token">,</SPAN> fld_tunnel_point_id<SPAN class="punctuation token">,</SPAN> fld_dist3D<SPAN class="punctuation token">)</SPAN>

&nbsp;&nbsp;&nbsp; ShowProgress<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"Start creating combinations..."</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">with</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>da<SPAN class="punctuation token">.</SPAN>InsertCursor<SPAN class="punctuation token">(</SPAN>tbl_out<SPAN class="punctuation token">,</SPAN> flds<SPAN class="punctuation token">)</SPAN> <SPAN class="keyword token">as</SPAN> curs<SPAN class="punctuation token">:</SPAN>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; cnt <SPAN class="operator token">=</SPAN> <SPAN class="number token">0</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; cnt_found <SPAN class="operator token">=</SPAN> <SPAN class="number token">0</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">for</SPAN> building_data <SPAN class="keyword token">in</SPAN> lst_building<SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; oid_building <SPAN class="operator token">=</SPAN> building_data<SPAN class="punctuation token">[</SPAN><SPAN class="number token">0</SPAN><SPAN class="punctuation token">]</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; pntg_building <SPAN class="operator token">=</SPAN> building_data<SPAN class="punctuation token">[</SPAN><SPAN class="number token">1</SPAN><SPAN class="punctuation token">]</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">for</SPAN> tunnel_data <SPAN class="keyword token">in</SPAN> lst_tunnelpoints<SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; oid_tunnel_point <SPAN class="operator token">=</SPAN> tunnel_data<SPAN class="punctuation token">[</SPAN><SPAN class="number token">0</SPAN><SPAN class="punctuation token">]</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; pntg_tunnel_point <SPAN class="operator token">=</SPAN> tunnel_data<SPAN class="punctuation token">[</SPAN><SPAN class="number token">1</SPAN><SPAN class="punctuation token">]</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; cnt <SPAN class="operator token">+=</SPAN> <SPAN class="number token">1</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">if</SPAN> cnt <SPAN class="operator token">%</SPAN> <SPAN class="number token">10000</SPAN> <SPAN class="operator token">==</SPAN> <SPAN class="number token">0</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ShowProgressPercentage<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"Processing combination: {0},&nbsp; found: {1}"</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN>cnt<SPAN class="punctuation token">,</SPAN> cnt_found<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">,</SPAN> cnt<SPAN class="punctuation token">,</SPAN> cnt_total<SPAN class="punctuation token">)</SPAN>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; dist3D <SPAN class="operator token">=</SPAN> get3Ddistance<SPAN class="punctuation token">(</SPAN>pntg_building<SPAN class="punctuation token">,</SPAN> pntg_tunnel_point<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">if</SPAN> dist3D <SPAN class="operator token"><=</SPAN> max_dist<SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; cnt_found <SPAN class="operator token">+=</SPAN> <SPAN class="number token">1</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; row <SPAN class="operator token">=</SPAN> <SPAN class="punctuation token">(</SPAN>oid_building<SPAN class="punctuation token">,</SPAN> oid_tunnel_point<SPAN class="punctuation token">,</SPAN> dist3D<SPAN class="punctuation token">,</SPAN> <SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; curs<SPAN class="punctuation token">.</SPAN>insertRow<SPAN class="punctuation token">(</SPAN>row<SPAN class="punctuation token">)</SPAN>

&nbsp;&nbsp;&nbsp; ShowProgressPercentage<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"Ready processing {0} combinations,&nbsp; found: {1}"</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN>cnt<SPAN class="punctuation token">,</SPAN> cnt_found<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">,</SPAN> cnt<SPAN class="punctuation token">,</SPAN> cnt_total<SPAN class="punctuation token">)</SPAN>


<SPAN class="keyword token">def</SPAN> <SPAN class="token function">AddField</SPAN><SPAN class="punctuation token">(</SPAN>fc<SPAN class="punctuation token">,</SPAN> fld_name<SPAN class="punctuation token">,</SPAN> fld_type<SPAN class="punctuation token">,</SPAN> fld_length<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">if</SPAN> len<SPAN class="punctuation token">(</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>ListFields<SPAN class="punctuation token">(</SPAN>fc<SPAN class="punctuation token">,</SPAN> fld_name<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN> <SPAN class="operator token">==</SPAN> <SPAN class="number token">0</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; arcpy<SPAN class="punctuation token">.</SPAN>AddField_management<SPAN class="punctuation token">(</SPAN>fc<SPAN class="punctuation token">,</SPAN> fld_name<SPAN class="punctuation token">,</SPAN> fld_type<SPAN class="punctuation token">,</SPAN> None<SPAN class="punctuation token">,</SPAN> None<SPAN class="punctuation token">,</SPAN> fld_length<SPAN class="punctuation token">)</SPAN>


<SPAN class="keyword token">def</SPAN> <SPAN class="token function">get3Ddistance</SPAN><SPAN class="punctuation token">(</SPAN>pntg1<SPAN class="punctuation token">,</SPAN> pntg2<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp; pnt1 <SPAN class="operator token">=</SPAN> pntg1<SPAN class="punctuation token">.</SPAN>firstPoint
&nbsp;&nbsp;&nbsp; pnt2 <SPAN class="operator token">=</SPAN> pntg2<SPAN class="punctuation token">.</SPAN>firstPoint
&nbsp;&nbsp;&nbsp; distance <SPAN class="operator token">=</SPAN> math<SPAN class="punctuation token">.</SPAN>sqrt<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">(</SPAN>pnt1<SPAN class="punctuation token">.</SPAN>X<SPAN class="operator token">-</SPAN>pnt2<SPAN class="punctuation token">.</SPAN>X<SPAN class="punctuation token">)</SPAN><SPAN class="operator token">**</SPAN><SPAN class="number token">2</SPAN> <SPAN class="operator token">+</SPAN> <SPAN class="punctuation token">(</SPAN>pnt1<SPAN class="punctuation token">.</SPAN>Y<SPAN class="operator token">-</SPAN>pnt2<SPAN class="punctuation token">.</SPAN>Y<SPAN class="punctuation token">)</SPAN><SPAN class="operator token">**</SPAN><SPAN class="number token">2</SPAN> <SPAN class="operator token">+</SPAN> <SPAN class="punctuation token">(</SPAN>pnt1<SPAN class="punctuation token">.</SPAN>Z<SPAN class="operator token">-</SPAN>pnt2<SPAN class="punctuation token">.</SPAN>Z<SPAN class="punctuation token">)</SPAN><SPAN class="operator token">**</SPAN><SPAN class="number token">2</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">return</SPAN> distance


<SPAN class="keyword token">def</SPAN> <SPAN class="token function">ShowProgressPercentage</SPAN><SPAN class="punctuation token">(</SPAN>txt<SPAN class="punctuation token">,</SPAN> cnt<SPAN class="punctuation token">,</SPAN> cnt_tot<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp; current_time <SPAN class="operator token">=</SPAN> strftime<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"%H:%M:%S"</SPAN><SPAN class="punctuation token">,</SPAN> time<SPAN class="punctuation token">.</SPAN>localtime<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">if</SPAN> cnt_tot <SPAN class="keyword token">is</SPAN> None<SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">print</SPAN> <SPAN class="string token">"{0}&nbsp;&nbsp; - {1}"</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN>current_time<SPAN class="punctuation token">,</SPAN> txt<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">else</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; percentage <SPAN class="operator token">=</SPAN> float<SPAN class="punctuation token">(</SPAN>cnt<SPAN class="punctuation token">)</SPAN> <SPAN class="operator token">*</SPAN> <SPAN class="number token">100.0</SPAN> <SPAN class="operator token">/</SPAN> float<SPAN class="punctuation token">(</SPAN>cnt_tot<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">print</SPAN> <SPAN class="string token">"{0}&nbsp;&nbsp; - {1} ({2}%)"</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN>current_time<SPAN class="punctuation token">,</SPAN> txt<SPAN class="punctuation token">,</SPAN> round<SPAN class="punctuation token">(</SPAN>percentage<SPAN class="punctuation token">,</SPAN> <SPAN class="number token">2</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>

<SPAN class="keyword token">def</SPAN> <SPAN class="token function">ShowProgress</SPAN><SPAN class="punctuation token">(</SPAN>txt<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># current_time = strftime("%Y-%m-%d %H:%M:%S", gmtime())</SPAN>
&nbsp;&nbsp;&nbsp; current_time <SPAN class="operator token">=</SPAN> strftime<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"%H:%M:%S"</SPAN><SPAN class="punctuation token">,</SPAN> time<SPAN class="punctuation token">.</SPAN>localtime<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">print</SPAN> <SPAN class="string token">"{0} - {1}"</SPAN><SPAN class="punctuation token">.</SPAN>format<SPAN class="punctuation token">(</SPAN>current_time<SPAN class="punctuation token">,</SPAN> txt<SPAN class="punctuation token">)</SPAN>

<SPAN class="keyword token">if</SPAN> __name__ <SPAN class="operator token">==</SPAN> <SPAN class="string token">'__main__'</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp; main<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">)</SPAN>
<SPAN class="line-numbers-rows"><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN></SPAN>

All comments

XanderBakker

Add a new field (double) and calculate the field using this formula:

Sqr (([POINT_X_building]-[POINT_X_constructionmachineno2850]^2 + &nbsp;[POINT_Y_building]-[POINT_Y_constructionmachineno2850]^2 + [POINT_Z_building]-[POINT_Z_constructionmachineno2850]^2 ))<SPAN class="line-numbers-rows"><SPAN>‍</SPAN></SPAN>

SB5

Hi again,

As mentioned in a comment to Xander above, I used 3D Buffer and Inside 3D plus some table joining to end up with a table with xyz coordinates for both building points and construction machine points. Can you tell me how to calculate 3D distance in a new field, either using excel or in the attribute table?

This is an example (with very many zeros for some reason)

POINT_X_building	POINT_Y_building	POINT_Z_building	POINT_X_constructionmachineno2850	POINT_Y_constructionmachineno2850	POINT_Z_constructionmachineno2850
599017.980859000000000	6640584.964650000400000	123.000000000000000	599075.213999999690000	6640639.921700000800000	-6.040000000000000
598973.976436999980000	6640633.025369999900000	121.950400000000000	599075.213999999690000	6640639.921700000800000	-6.040000000000000
598973.976436999980000	6640633.025369999900000	121.950400000000000	599075.213999999690000	6640639.921700000800000	-6.040000000000000
599111.963903000000000	6640589.985279999700000	126.000000000000000	599075.213999999690000	6640639.921700000800000	-6.040000000000000
599090.217955000000000	6640602.671959999900000	126.000000000000000	599075.213999999690000	6640639.921700000800000	-6.040000000000000

Thanks so much in advance.

SB5

Thanks, this seems interesting. I unfortunately only have the basic license and need advanced to use Generate Near Table.

SB5

The Z difference is between 5 and 100 meters.

Limiting the search radius would at the most cut the number of building points in two, leaving me with 6000 points.The building points are all within 400 meters of the line that the construction machine points run along, and I would need to select building points within 150-250 meters from each construction machine point.

I unfortunately cannot share the data.

Thanks again Xander.

DarrenWiens2

You could use Generate Near Table to calculate the 2d distances ('c' below), then Spatial Join to transfer start/end height attribute (difference is 'd' below), then calculate the 3d distance to a new field ('e' below).

XanderBakker

Just wondering... is there a lot of difference in Z between the building points and the construction machine points?

The fact that you will limit the search distance can limit a lot the number of results. Also if you apply a "Relate" and define your attribute indexes correctly, this might yield be a workable solution.

Is it possible to share your data, to see what's can be achieved?

Kind regards, Xander

SB5

Thanks Xander. I guess the easiest solution is to run Near 3D after specifying one location of the construction machine in the first feature class. My aim was to have this available for my colleagues to query from a table, but there would be too much data.

XanderBakker

The Near 3D—Help | ArcGIS for Desktop tool will get you the distance to the nearest object in the other layer, not all objects. The Generate Near Table—Help | ArcGIS for Desktop will calculate the distance between all objects in two featureclasses, but not in 3D.

I suppose you would have to script this, but this means evaluating 12000 x 22000 = 264M combinations. That will be a pretty long process. Also to work with the results will be slow, since your table will have 264M records.

To calculate the 3D distance between is not very complex. Just take the square root from (dx² + dy² + dz²).

Join the Conversation

Signed in members can post, follow updates, and more. New here? Register a free account.

Quick Links

Community Help
Find useful guides, FAQs, and documents to help you navigate and make the most of Esri Community.