topic Re: Path to .prj files in ArcGIS v10.1? in Python Questions

Path to .prj files in ArcGIS v10.1?

ChrisSnyder — Sun, 12 Dec 2021 16:33:14 GMT

In v10.0 and before there was always a folder that contained all the ESRI .prj definition files... Usually something like C:\Program Files\ArcGIS\Desktop10.0\Coordinate Systems.

Maybe I messed up my install somehow, but now in v10.1 it seems that 'Coordinate Systems' folder is no longer there, and to boot, there are no .prj files at all in the ArcGIS Program Files directory... at all!

I always used the .prj files to set the arcpy.env.outputCoordinateSystem variable (and all its various past incarnations).

What gives? How can we now access the official ESRI projection definitions in a script if there are no .prj definition file anymore?

The help says you can use .prj files, and the provided example in the arcpy.env.outputCoordinateSystem help topic yields:

arcpy.env.outputCoordinateSystem = "Coordinate Systems/Projected Coordinate Systems/UTM/WGS 1984/Northern Hemisphere/WGS 1984 UTM Zone 18N.prj"
Traceback (most recent call last):
  File "<interactive input>", line 1, in <module>
  File "C:\Program Files\ArcGIS\Desktop10.1\arcpy\arcpy\geoprocessing\_base.py", line 515, in set_
    self[env] = val
  File "C:\Program Files\ArcGIS\Desktop10.1\arcpy\arcpy\geoprocessing\_base.py", line 567, in __setitem__
    ret_ = setattr(self._gp, item, value)
RuntimeError: Object: Error in accessing environment <outputCoordinateSystem>

Am I totally missing something?

Re: Path to .prj files in ArcGIS v10.1?

ChrisSnyder — Wed, 08 Aug 2012 23:48:23 GMT

Thanks Kim...

http://forums.arcgis.com/threads/43644-UTM-Conversion-in-Python?p=222580&viewfull=1#post222580

Good ole' .prj... you will be missed!

Re: Path to .prj files in ArcGIS v10.1?

MathewCoyle — Sat, 11 Dec 2021 22:04:58 GMT

Well that throws a wrench in things, thanks for posting this so I have time to make some much needed changes before I even start the switch over to 10.1. How does this impact the project tool?

In 10.0 I would do essentially something like this.

outCS = os.path.join(installDir,r"Coordinate Systems\Projected Coordinate Systems\National Grids\Canada\NAD 1983 CSRS UTM Zone 12N.prj")
inCS = os.path.join(installDir,r"Coordinate Systems\Geographic Coordinate Systems\World\WGS 1984.prj")

arcpy.Project_management(in_fc, out_feature, outCS, transform_method, inCS)

In 10.1 would look like this?

inCS = arcpy.SpatialReference()
inCS.factoryCode = some_code1
inCS.create()

outCS = arcpy.SpatialReference()
outCS.factoryCode = some_code2
outCS.create()

arcpy.Project_management(in_fc, out_feature, outCS, transform_method, inCS)

Re: Path to .prj files in ArcGIS v10.1?

ChrisSnyder — Sun, 12 Dec 2021 16:33:17 GMT

Ya - i guess that's how it would work via scripting.

Sure is a pain to find the "codes" since the .pdf is not well structured like the Coordinate Systems directory tree.

Only way I could find the right code was by using v10.0, getting the sr object, and spitting out the "factory code" such as:

>>> arcpy.env.outputCoordinateSystem = r"C:\Program Files\ArcGIS\Desktop10.0\Coordinate Systems\Projected Coordinate Systems\State Plane\NAD 1983 HARN (US Feet)\NAD 1983 HARN StatePlane Washington South FIPS 4602 (US Feet).prj"
>>> arcpy.env.outputCoordinateSystem.factoryCode
2927

So, my new syntax in v10.1 is this:

arcpy.env.outputCoordinateSystem = 2927 #WA_SPS_NAD83HARN_USFEET

Note that (contrary to the Help topic) this syntax DOES NOT appear to work for me:

>>> arcpy.env.outputCoordinateSystem = arcpy.SpatialReference('NAD_1983_HARN_StatePlane_Washington_South_FIPS_4602_Feet')Traceback (most recent call last):
  File "<interactive input>", line 1, in <module>
  File "C:\Program Files\ArcGIS\Desktop10.0\arcpy\arcpy\geoprocessing\_base.py", line 505, in set_
    self[env] = val
  File "C:\Program Files\ArcGIS\Desktop10.0\arcpy\arcpy\geoprocessing\_base.py", line 557, in __setitem__
    ret_ = setattr(self._gp, item, value)
RuntimeError: <unprintable RuntimeError object>

Note sure why the .prj files went away... Dawning of a new age of fully integrated seamless geospatial technology I guess... Sigh...

COMPLAINT TO ESRI:The v10.1 Help topics that cover to projection-related stuff via Python scripting could use some further explanation about the new code/name based syntax. Script example syntax relating to using spatial reference names to set the SR appears to be incorrect.

Re: Path to .prj files in ArcGIS v10.1?

JasonScheirer — Thu, 09 Aug 2012 15:49:36 GMT

You can find the correct string for the spatial reference you want to create with arcpy.ListSpatialReferences.

Re: Path to .prj files in ArcGIS v10.1?

ChrisSnyder — Sat, 11 Dec 2021 22:05:01 GMT

Thanks Jason - That did work...

For example:

>>> srList = arcpy.ListSpatialReferences("*washington*south*", "PCS")
>>> for sr in srList:
...  print sr
Projected Coordinate Systems/State Plane/NAD 1927 (US Feet)/NAD 1927 StatePlane Washington South FIPS 4602
Projected Coordinate Systems/State Plane/NAD 1983 (CORS96) (Meters)/NAD 1983 (CORS96) StatePlane Washington South FIPS 4602 (Meters)
Projected Coordinate Systems/State Plane/NAD 1983 (CORS96) (US Feet)/NAD 1983 (CORS96) StatePlane Washington South FIPS 4602 (US Feet)
Projected Coordinate Systems/State Plane/NAD 1983 (Meters)/NAD 1983 StatePlane Washington South FIPS 4602 (Meters)
Projected Coordinate Systems/State Plane/NAD 1983 (US Feet)/NAD 1983 StatePlane Washington South FIPS 4602 (US Feet)
Projected Coordinate Systems/State Plane/NAD 1983 HARN (Meters)/NAD 1983 HARN StatePlane Washington South FIPS 4602 (Meters)
Projected Coordinate Systems/State Plane/NAD 1983 HARN (US Feet)/NAD 1983 HARN StatePlane Washington South FIPS 4602 (US Feet)
Projected Coordinate Systems/State Plane/NAD 1983 NSRS2007 (Meters)/NAD 1983 NSRS2007 StatePlane Washington South FIPS 4602 (Meters)
Projected Coordinate Systems/State Plane/NAD 1983 NSRS2007 (US Feet)/NAD 1983 NSRS2007 StatePlane Washington South FIPS 4602 (US Feet)
>>> arcpy.env.outputCoordinateSystem = arcpy.SpatialReference(srList[6])
>>>

This also works (not sure what I was doing wrong before):

arcpy.env.outputCoordinateSystem = arcpy.SpatialReference('NAD 1983 HARN StatePlane Washington South FIPS 4602 (US Feet)')

Re: Path to .prj files in ArcGIS v10.1?

ChrisSnyder — Thu, 09 Aug 2012 16:39:07 GMT

Oh see what's different...

The "sR.name" property appears to be slightly different than the "official" syntax compliant name that uis given by the arcpy.ListSpatialReferences() method.

Is there any sR property that yields the "syntax compliant" name?

sR.name property of 'NAD_1983_HARN_StatePlane_Washington_South_FIPS_4602_Feet'
vs.
compliant name of 'NAD 1983 HARN StatePlane Washington South FIPS 4602 (US Feet)'

Re: Path to .prj files in ArcGIS v10.1?

NeilAyres — Fri, 10 Aug 2012 11:59:55 GMT

Phew,
I am so glad I stumbled across this little gem. I was about to upgrade my system to 10.1 this w/end.
Now I am not so sure....
The old directory system of *.prj files worked and was quite flexible. ie It was easy to create a new prj for a particular project, say a custom crafted UTM like system (ie with a non-standard central meridian). Now what do you do?
It was also useful to browse the directories when searching for obscure projection info for certain regions.

Why make this info harder to find? And use. Oh and Yes, I can also memorise all xx thousand "factory codes"....

Miffed,
Neil

Re: Path to .prj files in ArcGIS v10.1?

PeterYurkosky1 — Fri, 10 Aug 2012 14:11:35 GMT

The old directory system of *.prj files worked and was quite flexible. ie It was easy to create a new prj for a particular project, say a custom crafted UTM like system (ie with a non-standard central meridian). Now what do you do?
It was also useful to browse the directories when searching for obscure projection info for certain regions.

Neil,

We've tried to make it much easier to search for coordinate systems. If you're working in the desktop apps, you can search on a text string (if you only know a fragment of the coordinate system's name), or apply a spatial filter to reduce the list to only those coordinate systems applicable to the map's extent. Once you've found that obscure coordinate system, you can add it to your Favorites and never have to search for it again.

Re: Path to .prj files in ArcGIS v10.1?

ChrisSnyder — Fri, 10 Aug 2012 16:12:49 GMT

Many (most/all?) of the scripting users were accustomed to using the file-based .prj definitions. Some issues I see:

1. The SR properties .name, .PCSName and/or .GCSName do not necessarily match the allowable spatial reference token keywords�?� For example the factory code of 2927 (the projection for all WA State Govt. agencies) has a .name and .PCSName property of 'NAD_1983_HARN_StatePlane_Washington_South_FIPS_4602_Feet', however the spatial reference token keyword is actually 'NAD 1983 HARN StatePlane Washington South FIPS 4602 (US Feet)'. So the effect of this issue is that this works:

arcpy.env.outputCoordinateSystem = arcpy.SpatialReference('NAD 1983 HARN StatePlane Washington South FIPS 4602 (US Feet)')

and this DOES NOT work:

arcpy.env.outputCoordinateSystem = arcpy.SpatialReference('NAD_1983_HARN_StatePlane_Washington_South_FIPS_4602_Feet')

Perhaps my projection (WA SPS NAD83HARN USFEET) is the exception here, but it would be nice to have a SR property that would yield the correct sr token keyword.

2. The Help system topics (Project tool, Project Raster tool, arcpy.env.outputCoordinateSystem, etc) do not adequately describe the disappearance of the file-based .prj files and it's impact on script-based use of the geoprocessing tools, nor does it describe the proper syntax of the factory codes (or how to search for them in projected coordinate systems .pdf reference doc), fails to mention the arcpy.ListSpatialReferences() method, and fails to mention a way to retrieve the proper sr token keywords. Also, script examples for the arcpy.env.outputCoordinateSystem method have the wrong syntax for using the spatial reference token keywords (Help provided example of arcpy.env.outputCoordinateSystem = "Coordinate Systems/Projected Coordinate Systems/UTM/WGS 1984/Northern Hemisphere/WGS 1984 UTM Zone 18N.prj" results in an error.

3. There are no file-based .prj files anymore. This threw me for a big loop! I suppose you can make some with the .exportToString() method, but that's kinda a pain... Seems like some of the legacy tools such as Project Coverage required .prj files... Not sure why the help system mentions .prj files anymore if there aren't any to speak of!

I liked the .prj file folders because they were easy to search... Other than the directory tree in the the coordinate system parameter GUI in the projection-related toolbox tools, I don't really see an easy way to browse and retrieve the factory codes or SR token keywords (which you now need to reference when specifying coordinate systems in a script).

Getting rid of the file-based .prj definitions is a pretty big change... Like everything else though I'm sure we'll all get used to it.

Re: Path to .prj files in ArcGIS v10.1?

PeterYurkosky1 — Fri, 10 Aug 2012 18:04:46 GMT

By way of suggestions for a workaround:

Keep in mind that you can still export coordinate systems to .prj files in the dialog (right click the coordinate system in the treeview and select Save As...).

You can also get the factory code in the same dialog (and avoid searching through that .pdf file). It's shown as the "WKID" in the large text box at the bottom.

Re: Path to .prj files in ArcGIS v10.1?

NeilAyres — Mon, 13 Aug 2012 09:57:23 GMT

Okay Pete,
As I havn't yet upgraded my system, is there a way to select a coordinate system, modify it then "Save As"?
But this is quite a change and I agree with the comments above.

N

Re: Path to .prj files in ArcGIS v10.1?

PeterYurkosky1 — Mon, 13 Aug 2012 13:26:27 GMT

is there a way to select a coordinate system, modify it then "Save As"?

Yes, just right-click the coordinate system in the tree-view and select Copy and Modify. This will create a new coordinate system that will appear under a "Custom" category. Then either add it to your "Favorites" (which creates a .prj file in your %APPDATA% folder) or select Save As to put it wherever you want.

Basically, the software still reads and honors .prj files, but they're no longer the mechanism for deploying the 4000+ coordinate systems included with ArcGIS.

Re: Path to .prj files in ArcGIS v10.1?

T__WayneWhitley — Sat, 11 Dec 2021 22:05:03 GMT

This may be a little late to post here, but I followed up on Chris Snyder's suggestion (post #10) to write the prj file using the exportToString() method. I didn't find it posted anywhere else and it's simple, so here it is (below). The topic came up again yesterday with a post by Janet Brewster.

This script version simply writes a single file...

- gets the SR object from an input dataset and exports it to string.

- replaces single quotes with double quotes in the string (mine was rejected when I tried to use it with 10.0; was accepted as the coord sys param by tools such as Create Feature Class after replacing the quotes).

- writes the result to a text file (use a prj extension).

# define a few variables 1st:
dataset = r'< fc or shp pathname you want to create prj from >'
outFile = r'< text file pathname you want to write to (make sure to use prj ext) >'

print 'executing...please wait.\n'
import arcpy, time

# get Spatial Reference object
SR = arcpy.Describe(dataset).spatialReference

# get the string from the object
seeString = SR.exportToString()

# Testing w/ system ArcGIS tools at 10.0 rejected the string - the problem:  single quotes, simply replace with double quotes:
seeString = seeString.replace('\'', '\"')

# let's see it...
print 'This is the string after replacing single quotes, if any found:\n'
print seeString + '\n'

# write it and use it!
writeTo = open(outFile, 'wb')
writeTo.writelines(seeString)
writeTo.close()

print 'done.  ...check your new prj file in the directory you defined.'
time.sleep(5)

Enjoy,
Wayne