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; } } }

Spatial correlation of two data sets

FrankHanssen

I have one modelled power-line path that I want to compare statistically with an existing power-line path. What I would like to measure is how the model output correlate spatially with the reality. So, no other variables than the XY-coordinates.

Best,

Frank

Finde mehr Beiträge mit dem Tag

Accepted answers

XanderBakker

I just wrote some code to see what can be done.

The code used two featureclasses with the two lines (it extracts the first feature as used the polyline) and it created two output featureclasses. The step_size (line 10) defines at what distance the points are extracted from the line:

<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> arcpy
&nbsp;&nbsp;&nbsp; arcpy<SPAN class="punctuation token">.</SPAN>env<SPAN class="punctuation token">.</SPAN>overwriteOutput <SPAN class="operator token">=</SPAN> <SPAN class="token boolean">True</SPAN>

&nbsp;&nbsp;&nbsp; fc_gis <SPAN class="operator token">=</SPAN> r<SPAN class="string token">'C:\GeoNet\compareLines\data.gdb\GIS_route'</SPAN>
&nbsp;&nbsp;&nbsp; fc_real <SPAN class="operator token">=</SPAN> r<SPAN class="string token">'C:\GeoNet\compareLines\data.gdb\Real_route'</SPAN>
&nbsp;&nbsp;&nbsp; fc_out1 <SPAN class="operator token">=</SPAN> r<SPAN class="string token">'C:\GeoNet\compareLines\data.gdb\dist_gis2real'</SPAN>
&nbsp;&nbsp;&nbsp; fc_out2 <SPAN class="operator token">=</SPAN> r<SPAN class="string token">'C:\GeoNet\compareLines\data.gdb\dist_real2gis'</SPAN>

&nbsp;&nbsp;&nbsp; step_size <SPAN class="operator token">=</SPAN> <SPAN class="number token">100</SPAN>

&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># get lines</SPAN>
&nbsp;&nbsp;&nbsp; line_base <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>da<SPAN class="punctuation token">.</SPAN>SearchCursor<SPAN class="punctuation token">(</SPAN>fc_real<SPAN class="punctuation token">,</SPAN> <SPAN class="punctuation token">(</SPAN><SPAN class="string token">"SHAPE@"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">.</SPAN>next<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">[</SPAN><SPAN class="number token">0</SPAN><SPAN class="punctuation token">]</SPAN>
&nbsp;&nbsp;&nbsp; line_compare <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>da<SPAN class="punctuation token">.</SPAN>SearchCursor<SPAN class="punctuation token">(</SPAN>fc_gis<SPAN class="punctuation token">,</SPAN> <SPAN class="punctuation token">(</SPAN><SPAN class="string token">"SHAPE@"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">.</SPAN>next<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">[</SPAN><SPAN class="number token">0</SPAN><SPAN class="punctuation token">]</SPAN>
&nbsp;&nbsp;&nbsp; steps <SPAN class="operator token">=</SPAN> int<SPAN class="punctuation token">(</SPAN>divmod<SPAN class="punctuation token">(</SPAN>line_base<SPAN class="punctuation token">.</SPAN>length<SPAN class="punctuation token">,</SPAN> step_size<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">[</SPAN><SPAN class="number token">0</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; sr <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>Describe<SPAN class="punctuation token">(</SPAN>fc_gis<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">.</SPAN>spatialReference

&nbsp;&nbsp;&nbsp; lst_feats <SPAN class="operator token">=</SPAN> <SPAN class="punctuation token">[</SPAN><SPAN class="punctuation token">]</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">for</SPAN> i <SPAN class="keyword token">in</SPAN> range<SPAN class="punctuation token">(</SPAN>steps<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; d <SPAN class="operator token">=</SPAN> i <SPAN class="operator token">*</SPAN> step_size
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; pnt <SPAN class="operator token">=</SPAN> line_base<SPAN class="punctuation token">.</SPAN>positionAlongLine<SPAN class="punctuation token">(</SPAN>d<SPAN class="punctuation token">,</SPAN> <SPAN class="token boolean">False</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; pnt_near <SPAN class="operator token">=</SPAN> line_compare<SPAN class="punctuation token">.</SPAN>queryPointAndDistance<SPAN class="punctuation token">(</SPAN>pnt<SPAN class="punctuation token">,</SPAN> <SPAN class="token boolean">False</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">[</SPAN><SPAN class="number token">0</SPAN><SPAN class="punctuation token">]</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; polyline <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>Polyline<SPAN class="punctuation token">(</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>Array<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">[</SPAN>pnt<SPAN class="punctuation token">.</SPAN>firstPoint<SPAN class="punctuation token">,</SPAN> pnt_near<SPAN class="punctuation token">.</SPAN>firstPoint<SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">,</SPAN> sr<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; lst_feats<SPAN class="punctuation token">.</SPAN>append<SPAN class="punctuation token">(</SPAN>polyline<SPAN class="punctuation token">)</SPAN>

&nbsp;&nbsp;&nbsp; arcpy<SPAN class="punctuation token">.</SPAN>CopyFeatures_management<SPAN class="punctuation token">(</SPAN>lst_feats<SPAN class="punctuation token">,</SPAN> fc_out1<SPAN class="punctuation token">)</SPAN>


&nbsp;&nbsp;&nbsp; steps <SPAN class="operator token">=</SPAN> int<SPAN class="punctuation token">(</SPAN>divmod<SPAN class="punctuation token">(</SPAN>line_compare<SPAN class="punctuation token">.</SPAN>length<SPAN class="punctuation token">,</SPAN> step_size<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">[</SPAN><SPAN class="number token">0</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">)</SPAN>

&nbsp;&nbsp;&nbsp; lst_feats <SPAN class="operator token">=</SPAN> <SPAN class="punctuation token">[</SPAN><SPAN class="punctuation token">]</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">for</SPAN> i <SPAN class="keyword token">in</SPAN> range<SPAN class="punctuation token">(</SPAN>steps<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; d <SPAN class="operator token">=</SPAN> i <SPAN class="operator token">*</SPAN> step_size
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; pnt <SPAN class="operator token">=</SPAN> line_compare<SPAN class="punctuation token">.</SPAN>positionAlongLine<SPAN class="punctuation token">(</SPAN>d<SPAN class="punctuation token">,</SPAN> <SPAN class="token boolean">False</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; pnt_near <SPAN class="operator token">=</SPAN> line_base<SPAN class="punctuation token">.</SPAN>queryPointAndDistance<SPAN class="punctuation token">(</SPAN>pnt<SPAN class="punctuation token">,</SPAN> <SPAN class="token boolean">False</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">[</SPAN><SPAN class="number token">0</SPAN><SPAN class="punctuation token">]</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; polyline <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>Polyline<SPAN class="punctuation token">(</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>Array<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">[</SPAN>pnt<SPAN class="punctuation token">.</SPAN>firstPoint<SPAN class="punctuation token">,</SPAN> pnt_near<SPAN class="punctuation token">.</SPAN>firstPoint<SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">,</SPAN> sr<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; lst_feats<SPAN class="punctuation token">.</SPAN>append<SPAN class="punctuation token">(</SPAN>polyline<SPAN class="punctuation token">)</SPAN>

&nbsp;&nbsp;&nbsp; arcpy<SPAN class="punctuation token">.</SPAN>CopyFeatures_management<SPAN class="punctuation token">(</SPAN>lst_feats<SPAN class="punctuation token">,</SPAN> fc_out2<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>

It creates two new featureclasses by projecting lines between points created on one line (based on the step_size) to the nearest point on the other line. This is what is created (using green for short distances and red for high distances):

If you look at the attribute table of one of the results you can run some statistics on the length of the segments being created:

Apart from analyzing the variation between them, it might be more interesting to understand why the real line does not follow the one calculated in GIS. Are there some variables that were not included? Are there features on the ground that were not considered by the GIS calculation since they may not have been available in the data used for the analysis?

All comments

FrankHanssen

Thanks Dan and Xander,

for all your valuable contributions! The modelled line is calculated based on a MCA-generated cost-surface and LCP based on a subset of criterion used in the planning of the original powerline. I have control about all criterion, my main issue was just to get comparison about how well they fitted. I did a similar calculation as you Xander in modelbuilder.

DanPatterson_Retired

Xander that is good.... yes, it is what's on the ground that should be of interest/concern for the predictive model,

XanderBakker

I just wrote some code to see what can be done.

<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> arcpy
&nbsp;&nbsp;&nbsp; arcpy<SPAN class="punctuation token">.</SPAN>env<SPAN class="punctuation token">.</SPAN>overwriteOutput <SPAN class="operator token">=</SPAN> <SPAN class="token boolean">True</SPAN>

&nbsp;&nbsp;&nbsp; fc_gis <SPAN class="operator token">=</SPAN> r<SPAN class="string token">'C:\GeoNet\compareLines\data.gdb\GIS_route'</SPAN>
&nbsp;&nbsp;&nbsp; fc_real <SPAN class="operator token">=</SPAN> r<SPAN class="string token">'C:\GeoNet\compareLines\data.gdb\Real_route'</SPAN>
&nbsp;&nbsp;&nbsp; fc_out1 <SPAN class="operator token">=</SPAN> r<SPAN class="string token">'C:\GeoNet\compareLines\data.gdb\dist_gis2real'</SPAN>
&nbsp;&nbsp;&nbsp; fc_out2 <SPAN class="operator token">=</SPAN> r<SPAN class="string token">'C:\GeoNet\compareLines\data.gdb\dist_real2gis'</SPAN>

&nbsp;&nbsp;&nbsp; step_size <SPAN class="operator token">=</SPAN> <SPAN class="number token">100</SPAN>

&nbsp;&nbsp;&nbsp; <SPAN class="comment token"># get lines</SPAN>
&nbsp;&nbsp;&nbsp; line_base <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>da<SPAN class="punctuation token">.</SPAN>SearchCursor<SPAN class="punctuation token">(</SPAN>fc_real<SPAN class="punctuation token">,</SPAN> <SPAN class="punctuation token">(</SPAN><SPAN class="string token">"SHAPE@"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">.</SPAN>next<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">[</SPAN><SPAN class="number token">0</SPAN><SPAN class="punctuation token">]</SPAN>
&nbsp;&nbsp;&nbsp; line_compare <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>da<SPAN class="punctuation token">.</SPAN>SearchCursor<SPAN class="punctuation token">(</SPAN>fc_gis<SPAN class="punctuation token">,</SPAN> <SPAN class="punctuation token">(</SPAN><SPAN class="string token">"SHAPE@"</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">.</SPAN>next<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">[</SPAN><SPAN class="number token">0</SPAN><SPAN class="punctuation token">]</SPAN>
&nbsp;&nbsp;&nbsp; steps <SPAN class="operator token">=</SPAN> int<SPAN class="punctuation token">(</SPAN>divmod<SPAN class="punctuation token">(</SPAN>line_base<SPAN class="punctuation token">.</SPAN>length<SPAN class="punctuation token">,</SPAN> step_size<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">[</SPAN><SPAN class="number token">0</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; sr <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>Describe<SPAN class="punctuation token">(</SPAN>fc_gis<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">.</SPAN>spatialReference

&nbsp;&nbsp;&nbsp; lst_feats <SPAN class="operator token">=</SPAN> <SPAN class="punctuation token">[</SPAN><SPAN class="punctuation token">]</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">for</SPAN> i <SPAN class="keyword token">in</SPAN> range<SPAN class="punctuation token">(</SPAN>steps<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; d <SPAN class="operator token">=</SPAN> i <SPAN class="operator token">*</SPAN> step_size
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; pnt <SPAN class="operator token">=</SPAN> line_base<SPAN class="punctuation token">.</SPAN>positionAlongLine<SPAN class="punctuation token">(</SPAN>d<SPAN class="punctuation token">,</SPAN> <SPAN class="token boolean">False</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; pnt_near <SPAN class="operator token">=</SPAN> line_compare<SPAN class="punctuation token">.</SPAN>queryPointAndDistance<SPAN class="punctuation token">(</SPAN>pnt<SPAN class="punctuation token">,</SPAN> <SPAN class="token boolean">False</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">[</SPAN><SPAN class="number token">0</SPAN><SPAN class="punctuation token">]</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; polyline <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>Polyline<SPAN class="punctuation token">(</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>Array<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">[</SPAN>pnt<SPAN class="punctuation token">.</SPAN>firstPoint<SPAN class="punctuation token">,</SPAN> pnt_near<SPAN class="punctuation token">.</SPAN>firstPoint<SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">,</SPAN> sr<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; lst_feats<SPAN class="punctuation token">.</SPAN>append<SPAN class="punctuation token">(</SPAN>polyline<SPAN class="punctuation token">)</SPAN>

&nbsp;&nbsp;&nbsp; arcpy<SPAN class="punctuation token">.</SPAN>CopyFeatures_management<SPAN class="punctuation token">(</SPAN>lst_feats<SPAN class="punctuation token">,</SPAN> fc_out1<SPAN class="punctuation token">)</SPAN>


&nbsp;&nbsp;&nbsp; steps <SPAN class="operator token">=</SPAN> int<SPAN class="punctuation token">(</SPAN>divmod<SPAN class="punctuation token">(</SPAN>line_compare<SPAN class="punctuation token">.</SPAN>length<SPAN class="punctuation token">,</SPAN> step_size<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">[</SPAN><SPAN class="number token">0</SPAN><SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">)</SPAN>

&nbsp;&nbsp;&nbsp; lst_feats <SPAN class="operator token">=</SPAN> <SPAN class="punctuation token">[</SPAN><SPAN class="punctuation token">]</SPAN>
&nbsp;&nbsp;&nbsp; <SPAN class="keyword token">for</SPAN> i <SPAN class="keyword token">in</SPAN> range<SPAN class="punctuation token">(</SPAN>steps<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">:</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; d <SPAN class="operator token">=</SPAN> i <SPAN class="operator token">*</SPAN> step_size
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; pnt <SPAN class="operator token">=</SPAN> line_compare<SPAN class="punctuation token">.</SPAN>positionAlongLine<SPAN class="punctuation token">(</SPAN>d<SPAN class="punctuation token">,</SPAN> <SPAN class="token boolean">False</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; pnt_near <SPAN class="operator token">=</SPAN> line_base<SPAN class="punctuation token">.</SPAN>queryPointAndDistance<SPAN class="punctuation token">(</SPAN>pnt<SPAN class="punctuation token">,</SPAN> <SPAN class="token boolean">False</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">[</SPAN><SPAN class="number token">0</SPAN><SPAN class="punctuation token">]</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; polyline <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>Polyline<SPAN class="punctuation token">(</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>Array<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">[</SPAN>pnt<SPAN class="punctuation token">.</SPAN>firstPoint<SPAN class="punctuation token">,</SPAN> pnt_near<SPAN class="punctuation token">.</SPAN>firstPoint<SPAN class="punctuation token">]</SPAN><SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">,</SPAN> sr<SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; lst_feats<SPAN class="punctuation token">.</SPAN>append<SPAN class="punctuation token">(</SPAN>polyline<SPAN class="punctuation token">)</SPAN>

&nbsp;&nbsp;&nbsp; arcpy<SPAN class="punctuation token">.</SPAN>CopyFeatures_management<SPAN class="punctuation token">(</SPAN>lst_feats<SPAN class="punctuation token">,</SPAN> fc_out2<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>

If you look at the attribute table of one of the results you can run some statistics on the length of the segments being created:

DanPatterson_Retired

Association... not correlation.

what parameters/variables are used to determine your predicted path? That is the question.

There is a strong 'association' between the actual and predicted path in several areas and not so much in others.

The question becomes, what is in the vicinity of the locations where the paths are similar? What is in the vicinity of where they are not?

What is driving your 'predicted path'. Nothing moves by itself, so how do you determine the locations that form the line.

A lot is missing from your description on how the path is derived and why it moves in the direction that it does forming a similar path and length.

XanderBakker

Perhaps creating a polygon from both lines and use the area divided by the length of the path to say something about the deviation between the two lines.

FrankHanssen

Thanks Dan! I am interested to get a measure on how well the model path fit spatially with the actual path. I have derived distance differences and have max-mean-min values. Is it possible to correlate the coordinates?

DanPatterson_Retired

Show the paths you have. Are you looking for deviations of modelled path vs actual path? I suspect this isn't a correlation in the true sense of the word since causation isn't a play here.

There are a number of factors that you could derive including distance differences, directional changes etc, but there is no reason to suspect that the actual path could be used to 'predict' the location of your modelled path since you don't have an independent and a dependent variable as is required by correlation

Nimm an der Unterhaltung teil

Angemeldete Mitglieder können Beiträge verfassen, Updates folgen und mehr. Neu hier? Registriere ein kostenloses Konto.

Schnellzugriffe

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