topic Re: Project a geometry to another geometry in memory in Geoprocessing Questions

Project a geometry to another geometry in memory

DarrylKlassen — Wed, 22 Feb 2012 22:17:56 GMT

Hello,
Is it possible to project a geometry to another geometry in memory using ARCPY. I am trying to take a XY coordinate and create in memory features to project the coordinate system. I have used the following code to do so, but the project statement doesn't work. I get ERROR 000944: Output feature class cannot be in the in_memory workspace.

import string,sys,os,arcpy

point = arcpy.Point(x_coord,y_coord)
Pnt_geom = arcpy.PointGeometry(point, in_coord_system)
Pnt_geom_out = arcpy.Geometry()

mytestlayer = arcpy.Project_management(Pnt_geom, Pnt_geom_out, out_coord_system, "#", in_coord_system)

desc = arcpy.Describe(mytestlayer)
print desc.extent.XMin
print desc.extent.YMin

# EOF

This code works if i put in a shapefile name in the place to the Pnt_geom_out variable.
My 4 variables (x_coord, y_coord, in_coord_system, out_coord_system) are all valid.

Thanks
Darryl Klassen
GIS Analyst
darryl.klassen@gmail.com

Re: Project a geometry to another geometry in memory

GraemeBrowning — Fri, 28 Mar 2014 05:54:57 GMT

Unfortunately, this seems to be a documented software limitation of in_memory workspaces.

I say unfortunately because I just hit it too

Re: Project a geometry to another geometry in memory

BruceHarold — Fri, 28 Mar 2014 14:13:34 GMT

Hi

The workaround is to project the geometry in-place with projectAs:

http://resources.arcgis.com/en/help/main/10.2/index.html#//018z00000070000000

Regards

Re: Project a geometry to another geometry in memory

GraemeBrowning — Fri, 28 Mar 2014 22:44:21 GMT

Hi

The workaround is to project the geometry in-place with projectAs:

http://resources.arcgis.com/en/help/main/10.2/index.html#//018z00000070000000

Regards

This looks like an excellent workaround for the original poster but I am not sure that it will be able to help me. The reason for that is that the feature class that I want to Project has (at a guess, because I am not onsite today to check) 10,000 or so polygons, and in this instance I cannot use arcpy.da because that client is at ArcGIS 10.0.