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

Euclidean Distance and Extent

anTonialcaraz

Hi,

I need to calculate some Euclidean Distances out of a list of FCs in a GDB.

I notice that when using an extent (in this case a global extent - see script) it does not work. It does work without specifying an extent but I do need to use one otherwise the resulting raster is a bit cut out at the edges.

When using the tool in ArcMap and specifying the extent in the Environments it does work.

Any help/suggestions will be greatly appreciated.

<SPAN class="keyword token">import</SPAN> arcpy
<SPAN class="keyword token">from</SPAN> arcpy <SPAN class="keyword token">import</SPAN> env
arcpy<SPAN class="punctuation token">.</SPAN>env<SPAN class="punctuation token">.</SPAN>overwriteOutput<SPAN class="operator token">=</SPAN><SPAN class="token boolean">True</SPAN>
<SPAN class="keyword token">from</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>sa <SPAN class="keyword token">import</SPAN><SPAN class="operator token">*</SPAN>
arcpy<SPAN class="punctuation token">.</SPAN>CheckOutExtension<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"Spatial"</SPAN><SPAN class="punctuation token">)</SPAN>

arcpy<SPAN class="punctuation token">.</SPAN>env<SPAN class="punctuation token">.</SPAN>extent <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>Extent<SPAN class="punctuation token">(</SPAN><SPAN class="operator token">-</SPAN><SPAN class="number token">180.0</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="operator token">-</SPAN><SPAN class="number token">90.0</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="number token">180.0</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="number token">90.0</SPAN><SPAN class="punctuation token">)</SPAN>

env<SPAN class="punctuation token">.</SPAN>workspace <SPAN class="operator token">=</SPAN> r<SPAN class="string token">"H:\OUTPUT.gdb"</SPAN>

fcdist <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>ListFeatureClasses<SPAN class="punctuation token">(</SPAN><SPAN class="string token">""</SPAN><SPAN class="punctuation token">)</SPAN>

<SPAN class="keyword token">for</SPAN> fc <SPAN class="keyword token">in</SPAN> fcdist<SPAN class="punctuation token">:</SPAN>

&nbsp;&nbsp;&nbsp; output <SPAN class="operator token">=</SPAN> EucDistance<SPAN class="punctuation token">(</SPAN>fc<SPAN class="punctuation token">,</SPAN> <SPAN class="number token">250000</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="number token">10000</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">""</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; output<SPAN class="punctuation token">.</SPAN>save<SPAN class="punctuation token">(</SPAN>env<SPAN class="punctuation token">.</SPAN>workspace <SPAN class="operator token">+</SPAN> <SPAN class="string token">"\\"</SPAN> <SPAN class="operator token">+</SPAN> fc <SPAN class="operator token">+</SPAN> <SPAN class="string token">"_250km"</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>

Find more posts tagged with

Accepted answers

anTonialcaraz

It works like this:

arcpy.env.extent = arcpy.Extent(-18040095.1961317, -9020047.34806998, 18039904.8038683, 9019952.65193002)

Mollweide projection

All comments

anTonialcaraz

It works like this:

arcpy.env.extent = arcpy.Extent(-18040095.1961317, -9020047.34806998, 18039904.8038683, 9019952.65193002)

Mollweide projection

anTonialcaraz

Just to add:

The script works with the FCs in Geographic Coordinates (degrees) using the global extent (below)

I would expect the same in Mollweide Projected (meters) especially as you are measuring distances...

arcpy<SPAN class="punctuation token">.</SPAN>env<SPAN class="punctuation token">.</SPAN>extent <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>Extent<SPAN class="punctuation token">(</SPAN><SPAN class="operator token">-</SPAN><SPAN class="number token">180.0</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="operator token">-</SPAN><SPAN class="number token">90.0</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="number token">180.0</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="number token">90.0</SPAN><SPAN class="punctuation token">)</SPAN>

env<SPAN class="punctuation token">.</SPAN>workspace <SPAN class="operator token">=</SPAN> r<SPAN class="string token">"H:\OUTPUT.gdb"</SPAN>

fcdist <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>ListFeatureClasses<SPAN class="punctuation token">(</SPAN><SPAN class="string token">""</SPAN><SPAN class="punctuation token">)</SPAN>

<SPAN class="keyword token">for</SPAN> fc <SPAN class="keyword token">in</SPAN> fcdist<SPAN class="punctuation token">:</SPAN>

&nbsp;&nbsp;&nbsp; output <SPAN class="operator token">=</SPAN> EucDistance<SPAN class="punctuation token">(</SPAN>fc<SPAN class="punctuation token">,</SPAN> <SPAN class="number token">5</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="number token">0.5</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">""</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; output<SPAN class="punctuation token">.</SPAN>save<SPAN class="punctuation token">(</SPAN>env<SPAN class="punctuation token">.</SPAN>workspace <SPAN class="operator token">+</SPAN> <SPAN class="string token">"\\"</SPAN> <SPAN class="operator token">+</SPAN> fc <SPAN class="operator token">+</SPAN> <SPAN class="string token">"_250km"</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>

anTonialcaraz

To be more precise when using the tool in ArcMap session:

It works when I use the "same as layer..." (in this case I use a Bathymetry raster layer) in the Processing Extent option within Environment Settings.

If I use " As specified Below" option and I enter the extent manually it does not work.

I don't quite get it to be honest

XanderBakker

If I look at the domain property of the Mollweide (world) spatial reference I see this:

-18040200 -9020200 900701885274,099 900710905274,099

When you work inside a session of for instance ArcMap, a number of environment setting will be set based on your current dataframe settings. When you run a standalone script these settings are not there, so behavior will be different.

anTonialcaraz

As I previously mentioned, no problem at all when using the tool in ArcMap and setting the global Extent

anTonialcaraz

Thanks Dan.

The Euclidean Distance calculation is the end part of a larger script. In the previous step I do project the FCs to Mollweide. It is just when using the Extent that does not work but I guess it may still be related with the problem?

<SPAN class="comment token"># reproject to Mollweide and create 250km buffer</SPAN>

fcprj <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>ListFeatureClasses<SPAN class="punctuation token">(</SPAN><SPAN class="punctuation token">)</SPAN>

<SPAN class="keyword token">for</SPAN> fc <SPAN class="keyword token">in</SPAN> fcprj<SPAN class="punctuation token">:</SPAN>

&nbsp;&nbsp;&nbsp; mollweide <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>SpatialReference<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"Mollweide (world)"</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; arcpy<SPAN class="punctuation token">.</SPAN>BatchProject_management<SPAN class="punctuation token">(</SPAN>fc<SPAN class="punctuation token">,</SPAN> OUTworkspace<SPAN class="punctuation token">,</SPAN> mollweide<SPAN class="punctuation token">,</SPAN> <SPAN class="string token">""</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">""</SPAN><SPAN class="punctuation token">)</SPAN>

fcdist <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>ListFeatureClasses<SPAN class="punctuation token">(</SPAN><SPAN class="string token">"*1"</SPAN><SPAN class="punctuation token">)</SPAN>

<SPAN class="keyword token">for</SPAN> fc <SPAN class="keyword token">in</SPAN> fcdist<SPAN class="punctuation token">:</SPAN>

&nbsp;&nbsp;&nbsp; output <SPAN class="operator token">=</SPAN> EucDistance<SPAN class="punctuation token">(</SPAN>fc<SPAN class="punctuation token">,</SPAN> <SPAN class="number token">250000</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="number token">10000</SPAN><SPAN class="punctuation token">,</SPAN> <SPAN class="string token">""</SPAN><SPAN class="punctuation token">)</SPAN>
&nbsp;&nbsp;&nbsp; output<SPAN class="punctuation token">.</SPAN>save<SPAN class="punctuation token">(</SPAN>OUTworkspace <SPAN class="operator token">+</SPAN> <SPAN class="string token">"\\"</SPAN> <SPAN class="operator token">+</SPAN> fc <SPAN class="operator token">+</SPAN> <SPAN class="string token">"_250km"</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>

DanPatterson_Retired

Euclidean distance requires projected coordinates, it performs no projection during the process, you have to project the input data before using the tool, if the data are in geographic. Setting your extent to geographic is going to cause issues as you have found.

Euclidean - planar - projected …. otherwise things are .... spherical or geodesic

anTonialcaraz

Projected Coordinate System: World_Mollweide
Projection: Mollweide
False_Easting: 0.00000000
False_Northing: 0.00000000
Central_Meridian: 0.00000000
Linear Unit: Meter

Geographic Coordinate System: GCS_WGS_1984
Datum: D_WGS_1984
Prime Meridian: Greenwich
Angular Unit: Degree

XanderBakker

What coordinate system are you working with?

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.