So, i have a python function defined that projects coordinates using the 'projectAs()' point geometry function. When i call this function in a script from a toolbox it gives different results vs calling it in a python command prompt window. See below python function below: def ProjectCoordinates ( x , y , inputcrs , outputcrs ) : import arcpy point = arcpy . PointGeometry ( arcpy . Point ( x , y ) , inputcrs ) projectedpoint = point . projectAs ( outputcrs ) projectedcoordinates = [ projectedpoint . firstPoint . X , projectedpoint . firstPoint . Y ] return projectedcoordinates Full Toolbox Script: def ProjectCoordinates ( x , y , inputcrs , outputcrs ) : import arcpy point = arcpy . PointGeometry ( arcpy . Point ( x , y ) , inputcrs ) projectedpoint = point . projectAs ( outputcrs ) projectedcoordinates = [ projectedpoint . firstPoint . X , projectedpoint . firstPoint . Y ] return projectedcoordinates import arcpy xin = arcpy . GetParameterAsText ( 0 ) yin = arcpy . GetParameterAsText ( 1 ) srin = arcpy . GetParameter ( 2 ) srout = arcpy . GetParameter ( 3 ) xout , yout = ProjectCoordinates ( xin , yin , srin , srout ) arcpy . AddMessage ( srin . name ) arcpy . AddMessage ( srout . name ) arcpy . AddMessage ( str ( xout ) + " ," + str ( yout ) ) Results of Each: Can anyone replicate this? Does anyone have any ideas why this is happening? Also, if it helps, I know from surveyed data that the toolbox script output is correct and the command prompt output is incorrect.

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

Using 'projectAs()' Point Geometry Function Via Python Function in Command Prompt

BrandonZaitz

So, i have a python function defined that projects coordinates using the 'projectAs()' point geometry function. When i call this function in a script from a toolbox it gives different results vs calling it in a python command prompt window. See below python function below:

<SPAN class="keyword token">def</SPAN> <SPAN class="token function">ProjectCoordinates</SPAN><SPAN class="punctuation token">(</SPAN>x<SPAN class="punctuation token">,</SPAN>y<SPAN class="punctuation token">,</SPAN>inputcrs<SPAN class="punctuation token">,</SPAN>outputcrs<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">:</SPAN>
	<SPAN class="keyword token">import</SPAN> arcpy
	point <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>PointGeometry<SPAN class="punctuation token">(</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>Point<SPAN class="punctuation token">(</SPAN>x<SPAN class="punctuation token">,</SPAN>y<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">,</SPAN>inputcrs<SPAN class="punctuation token">)</SPAN>
	projectedpoint <SPAN class="operator token">=</SPAN> point<SPAN class="punctuation token">.</SPAN>projectAs<SPAN class="punctuation token">(</SPAN>outputcrs<SPAN class="punctuation token">)</SPAN>
	projectedcoordinates <SPAN class="operator token">=</SPAN> <SPAN class="punctuation token">[</SPAN>projectedpoint<SPAN class="punctuation token">.</SPAN>firstPoint<SPAN class="punctuation token">.</SPAN>X<SPAN class="punctuation token">,</SPAN>projectedpoint<SPAN class="punctuation token">.</SPAN>firstPoint<SPAN class="punctuation token">.</SPAN>Y<SPAN class="punctuation token">]</SPAN>
	
	<SPAN class="keyword token">return</SPAN> projectedcoordinates<SPAN class="line-numbers-rows"><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN></SPAN>

Full Toolbox Script:

<SPAN class="keyword token">def</SPAN> <SPAN class="token function">ProjectCoordinates</SPAN><SPAN class="punctuation token">(</SPAN>x<SPAN class="punctuation token">,</SPAN>y<SPAN class="punctuation token">,</SPAN>inputcrs<SPAN class="punctuation token">,</SPAN>outputcrs<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">:</SPAN>
	<SPAN class="keyword token">import</SPAN> arcpy
	point <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>PointGeometry<SPAN class="punctuation token">(</SPAN>arcpy<SPAN class="punctuation token">.</SPAN>Point<SPAN class="punctuation token">(</SPAN>x<SPAN class="punctuation token">,</SPAN>y<SPAN class="punctuation token">)</SPAN><SPAN class="punctuation token">,</SPAN>inputcrs<SPAN class="punctuation token">)</SPAN>
	projectedpoint <SPAN class="operator token">=</SPAN> point<SPAN class="punctuation token">.</SPAN>projectAs<SPAN class="punctuation token">(</SPAN>outputcrs<SPAN class="punctuation token">)</SPAN>
	projectedcoordinates <SPAN class="operator token">=</SPAN> <SPAN class="punctuation token">[</SPAN>projectedpoint<SPAN class="punctuation token">.</SPAN>firstPoint<SPAN class="punctuation token">.</SPAN>X<SPAN class="punctuation token">,</SPAN>projectedpoint<SPAN class="punctuation token">.</SPAN>firstPoint<SPAN class="punctuation token">.</SPAN>Y<SPAN class="punctuation token">]</SPAN>
	
	<SPAN class="keyword token">return</SPAN> projectedcoordinates

<SPAN class="keyword token">import</SPAN> arcpy

xin <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>GetParameterAsText<SPAN class="punctuation token">(</SPAN><SPAN class="number token">0</SPAN><SPAN class="punctuation token">)</SPAN>
yin <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>GetParameterAsText<SPAN class="punctuation token">(</SPAN><SPAN class="number token">1</SPAN><SPAN class="punctuation token">)</SPAN>

srin <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>GetParameter<SPAN class="punctuation token">(</SPAN><SPAN class="number token">2</SPAN><SPAN class="punctuation token">)</SPAN>
srout <SPAN class="operator token">=</SPAN> arcpy<SPAN class="punctuation token">.</SPAN>GetParameter<SPAN class="punctuation token">(</SPAN><SPAN class="number token">3</SPAN><SPAN class="punctuation token">)</SPAN>

xout<SPAN class="punctuation token">,</SPAN>yout <SPAN class="operator token">=</SPAN> ProjectCoordinates<SPAN class="punctuation token">(</SPAN>xin<SPAN class="punctuation token">,</SPAN>yin<SPAN class="punctuation token">,</SPAN>srin<SPAN class="punctuation token">,</SPAN>srout<SPAN class="punctuation token">)</SPAN>

arcpy<SPAN class="punctuation token">.</SPAN>AddMessage<SPAN class="punctuation token">(</SPAN>srin<SPAN class="punctuation token">.</SPAN>name<SPAN class="punctuation token">)</SPAN>
arcpy<SPAN class="punctuation token">.</SPAN>AddMessage<SPAN class="punctuation token">(</SPAN>srout<SPAN class="punctuation token">.</SPAN>name<SPAN class="punctuation token">)</SPAN>
arcpy<SPAN class="punctuation token">.</SPAN>AddMessage<SPAN class="punctuation token">(</SPAN>str<SPAN class="punctuation token">(</SPAN>xout<SPAN class="punctuation token">)</SPAN> <SPAN class="operator token">+</SPAN> <SPAN class="string token">" ,"</SPAN> <SPAN class="operator token">+</SPAN> str<SPAN class="punctuation token">(</SPAN>yout<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>

Results of Each:

Can anyone replicate this? Does anyone have any ideas why this is happening?

Also, if it helps, I know from surveyed data that the toolbox script output is correct and the command prompt output is incorrect.

Find more posts tagged with

Accepted answers

qnadeem

Well, I have struggled a bit without any luck, I was Passing the name of Spatial Reference System it did not work. Problem is displayed name in the Spatial Reference System is Different than the name used in this function projectAs().

I solved it by passing the Spatial Reference Object for example I created outPut coordinate System using

target_Ref = arcpy.SpatialReference(4326)

point.projectAs(target_Ref)

it worked fine

All comments

qnadeem

I solved it by passing the Spatial Reference Object for example I created outPut coordinate System using

target_Ref = arcpy.SpatialReference(4326)

point.projectAs(target_Ref)

it worked fine

RandyBurton

I got this to work in 10.6:

import arcpy

x = -101.415206
y = 32.537097

ptGeometry = arcpy.PointGeometry(arcpy.Point(x,y),arcpy.SpatialReference('NAD 1927')).projectAs("NAD 1983 StatePlane Texas N Central FIPS 4202 (US Feet)","NAD_1927_To_NAD_1983_NADCON")
  
print ptGeometry.firstPoint.X, ptGeometry.firstPoint.Y  
# result: 1070116.54144 6890849.30626
<SPAN class="line-numbers-rows"><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN><SPAN>‍</SPAN></SPAN>

It seems that when you use a "transformation_name", the associated "spatial_reference" must be a name (not a factory code). This document ArcGIS 10.7.1 and ArcGIS Pro 2.4 Geographic and Vertical Transformation Tables was helpful.

DanPatterson_Retired

Of course! The help file topic has the required info. This is for the generic Geometry Class... it trickles down to the other classes

Geometry—ArcPy classes | ArcGIS Desktop

projectAs (spatial_reference, {transformation_name})

Projects a geometry and optionally applies a geotransformation.

To project, the geometry needs to have a spatial reference, and not have an UnknownCoordinateSystem. The new spatial reference system passed to the method defines the output coordinate system. If either spatial reference is unknown the coordinates will not be changed. The Z- and measure values are not changed by the ProjectAs method.

RandyBurton

I did some testing and can confirm that the tool script and code run in the Python window produce different results. I believe this is due to an automatic transformation that is occurring when the tool is being run. The GCS NA 1927 coordinates appear to be transformed into GCS NA 1983 using NAD 1927 To NAD 1983 NADCON.

As an experiment I created two layers: the Texas State Plane that you are working with, and a GCS NA 1927. The state plane was populated with two points located at the positions you highlighted in your results. The GCS 1927 was populated with one point using the lon/lat in your script. The state plane was placed in a new map. When adding the GCS 1927 layer a transformation dialog was presented. If no transformation was selected, the GCS 1927 point was placed over the state plane point located at 1070240.5212705112, 6890805.12915604. If the NAD 1297 To NAD 1983 NADCON transformation was used, the GCS 1927 point was located over the state plane point located at 1070116.54144, 6890849.30626.

While experimenting, I noticed that the Add Geometry Attributes tool also automatically included the transformation.

I suppose the next question is: Can the transformation be added to the projectAs() function?

DanPatterson_Retired

datum shift... I am thinking that the plain script doesn't deal with it, but the toolbox script is 'smart' enough to detect it.

Worth a check anyway

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.