topic Re: Modify map coordinate system with arcpy.mp? in ArcGIS Pro Questions

Modify map coordinate system with arcpy.mp?

BryanChastain — Mon, 31 Oct 2016 15:31:11 GMT

Is there any way to modify a map's coordinate system in ArcGIS Pro using Python?

Re: Modify map coordinate system with arcpy.mp?

DanPatterson_Retired — Mon, 31 Oct 2016 17:41:59 GMT

Define Projection or the Project tool are your options ... what are you looking to modify?

Re: Modify map coordinate system with arcpy.mp?

JoshuaBixby — Mon, 31 Oct 2016 17:50:14 GMT

If by "map" you mean MapFrame, at ArcGIS Pro 1.3 there is no spatialReference method like there is for a DataFrame in ArcMap. If you mean Map by "map," same issue, no method to access spatial reference. Time to submit an enhancement request and ArcGIS Idea....

Re: Modify map coordinate system with arcpy.mp?

RebeccaStrauch__GISP — Mon, 31 Oct 2016 18:19:58 GMT

My guess is you know you can already manually make the change

Specify a coordinate system—Properties of maps | ArcGIS for Desktop

(I think if you have a default basemap with a different spatial reference, that will also change to the basemaps spatial reference as default....one way to change it)

as Dan asked, what are you looking to modify? And to add, are you trying to do this once, or on-the-fly. Remember that all your layers may need to be reprojected on the fly as well.

some links to look at

Output Coordinate System (Environment setting)—Geoprocessing | ArcGIS for Desktop

Using the spatial reference class—ArcPy Get Started | ArcGIS for Desktop

Specify a coordinate system—Properties of maps | ArcGIS for Desktop

SpatialReference—ArcPy Classes | ArcGIS for Desktop

Re: Modify map coordinate system with arcpy.mp?

BryanChastain — Mon, 31 Oct 2016 18:28:46 GMT

Rebecca Strauch, GISP & Dan_Patterson, sorry I meant change the spatial reference for either an arcpy.mp.Map or arcpy.mp.MapFrame.

According to bixb0012‌, looks like it's not possible. If that's true, that's a shame. It will be hard to come up with a Python-based workaround for that limitation.

Re: Modify map coordinate system with arcpy.mp?

DanPatterson_Retired — Mon, 31 Oct 2016 18:34:57 GMT

add an empty data frame, add a layer in the coordinate system you want... the dataframe coordinate system is now set... add/copy your other layers over and remove the 'seed' layer... done... depends how hard you need to have it done

Re: Modify map coordinate system with arcpy.mp?

JoshuaBixby — Mon, 31 Oct 2016 18:42:04 GMT

Unlike most ArcPy modules, ArcPy Mapping (arcpy.mp) had to be rewritten from the ground up since the UI and structure of ArcGIS Pro is completely different than ArcMap. They will be playing catch up for quite some time, which is why it is important for people to open enhancement requests and create ArcGIS Ideas. There is a lot to port over from ArcMap to ArcGIS Pro, and without user feedback, Esri is left to either do what is easiest first or guess at what users want.

Re: Modify map coordinate system with arcpy.mp?

BryanChastain — Mon, 31 Oct 2016 18:44:07 GMT

Pro doesn't have data frames. Do you mean a Map? The only way I see to programmatically create a new map is to call the arcpy.mp.ArcGISProject.importDocument() function. However, that would require you to have an mxd pre-authored for any coordinate system the user might want to change to. Then you'd still have to copy over all the existing layers into the new map. Also, if the map had already been used in a MapFrame in a Layout, then you'd really run into some issues...

Re: Modify map coordinate system with arcpy.mp?

BryanChastain — Tue, 22 Jun 2021 08:44:04 GMT

OK, vote me up!
https://community.esri.com/ideas/12767

Re: Modify map coordinate system with arcpy.mp?

JoshuaBixby — Mon, 31 Oct 2016 18:51:55 GMT

Done. It looked like I was the first vote. Can you vote for your own idea?

Re: Modify map coordinate system with arcpy.mp?

KeithMiller4 — Tue, 26 Jun 2018 22:12:22 GMT

I'm trying to do the same thing - set the coordinate system of an ArcGIS Pro Map object using arcpy. Just downloaded ArcGIS Pro 2.2 and there still doesn't seem to be a way. Anybody come across any solutions / workarounds in the last couple of years?

Re: Modify map coordinate system with arcpy.mp?

DavidClark1 — Fri, 29 Jun 2018 04:17:57 GMT

It now appears that there is a "Camera" object associated to Map and MapFrame objects. Within the Camera class there is a getExtent() and setExtent() method that returns or sets and Extent object, which has a spatialReference property.

Re: Modify map coordinate system with arcpy.mp?

JosephMcGehan — Wed, 23 Oct 2019 06:24:58 GMT

I need to read the map objects spatialReference properties. Picture that it's 2019 and Esri have GIS Software where you can't read spatial reference of a map.

aprx_object = arcpy.mp.ArcGISProject(aprx)
m = aprx_object.listMaps()
for mp in m:
   ID = 0
   mapname = mp.name
   mp_projection = mp.spatialReference() #

Re: Modify map coordinate system with arcpy.mp?

MaxSquires — Sun, 12 Dec 2021 05:29:00 GMT

you can 'get' the spatial reference from a map in the project like this:

aprx = arcpy.mp.ArcGISProject("c:\esri\project.aprx")

maps = aprx.listMaps()

map = maps[0]

print(dir(map)) 

#...

print(dir(map.getDefinition("V2").spatialReference))

‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍

and, while you cannot 'set' the spatial reference using 'arcpy', you can with a few additional libraries:

import arcpy
import xml.etree.ElementTree as ET
import tempfile
import os
import shutil
from zipfile import ZipFile

print("setting spatial reference")
output_project = "MASTER"
out_dir = os.path.dirname(__file__)

new_aprx_path = os.path.join(out_dir, output_project + "_Project", output_project) + ".aprx"
tempDirMod = tempfile.mkdtemp()
print("tempDirMod:")
print(tempDirMod)
with ZipFile(new_aprx_path, 'r') as zipObj:
   # Extract all the contents of aprx (zip file) in temporary directory
   zipObj.extractall(tempDirMod)

aprx = arcpy.mp.ArcGISProject(new_aprx_path)
maps = aprx.listMaps()

for m in maps:
    extracted_map_CIM_path = m.getDefinition("V2").uRI #i.e. 'CIMPATH=map11/map11.xml' --> "C:\admin\create_query_layers\MASTER_Project\MASTER\map11\map11.xml"
    xml_path = extracted_map_CIM_path.split("=")[1]
    this_map_xml = os.path.join(tempDirMod, xml_path)
    
    tree = ET.parse(this_map_xml)
    root = tree.getroot()
       
    DefaultExtent = root.find("DefaultExtent")
    Xmin = ET.SubElement(DefaultExtent, "XMin")
    Xmin.text = "-180"
    Ymin = ET.SubElement(DefaultExtent, "YMin")
    Ymin.text = "-90"
    Xmax = ET.SubElement(DefaultExtent, "XMax")
    Xmax.text = "180"
    Ymax = ET.SubElement(DefaultExtent, "YMax")
    Ymax.text = "90"
    SpatialReference = ET.SubElement(DefaultExtent, "SpatialReference")
    SpatialReference.set("xsi:type", "typens:GeographicCoordinateSystem")
    WKT = ET.SubElement(SpatialReference, "WKT")
    WKT.text = r"GEOGCS[&quot;GCS_WGS_1984&quot;,DATUM[&quot;D_WGS_1984&quot;,SPHEROID[&quot;WGS_1984&quot;,6378137.0,298.257223563]],PRIMEM[&quot;Greenwich&quot;,0.0],UNIT[&quot;Degree&quot;,0.0174532925199433],AUTHORITY[&quot;EPSG&quot;,4326]]" 
    XOrigin = ET.SubElement(SpatialReference, "XOrigin")
    XOrigin.text = "-400"
    YOrigin = ET.SubElement(SpatialReference, "YOrigin")
    YOrigin.text = "-400"
    XYScale = ET.SubElement(SpatialReference, "XYScale")
    XYScale.text = "999999999.99999988"
    ZOrigin = ET.SubElement(SpatialReference, "ZOrigin")
    ZOrigin.text = "-100000"
    ZScale = ET.SubElement(SpatialReference, "ZScale")
    ZScale.text = "10000"
    MOrigin = ET.SubElement(SpatialReference, "MOrigin")
    MOrigin.text = "-100000"
    MScale = ET.SubElement(SpatialReference, "MScale")
    MScale.text = "10000"
    XYTolerance = ET.SubElement(SpatialReference, "XYTolerance")
    XYTolerance.text = "8.983152841195215e-09"
    ZTolerance = ET.SubElement(SpatialReference, "ZTolerance")
    ZTolerance.text = "0.001"
    MTolerance = ET.SubElement(SpatialReference, "MTolerance")
    MTolerance.text = "0.001"
    HighPrecision = ET.SubElement(SpatialReference, "HighPrecision")
    HighPrecision.text = "true"
    LeftLongitude = ET.SubElement(SpatialReference, "LeftLongitude")
    LeftLongitude.text = "-180"
    WKID = ET.SubElement(SpatialReference, "WKID")
    WKID.text = "4326"
    LatestWKID = ET.SubElement(SpatialReference, "LatestWKID")
    LatestWKID.text = "4326"
    
    #IMPORTANT NOTE: if you have no spatial reference in your map yet you will want to append a new Element to the xml (as below).  Otherwise do as above with DefaultExtent where you would find the SpatialRefrence Node and overwrite existing data.
    SpatialReference2 =ET.Element("SpatialReference")
    SpatialReference2.set("xsi:type", "typens:GeographicCoordinateSystem")
    WKT = ET.SubElement(SpatialReference2, "WKT")
    WKT.text = r"GEOGCS[&quot;GCS_WGS_1984&quot;,DATUM[&quot;D_WGS_1984&quot;,SPHEROID[&quot;WGS_1984&quot;,6378137.0,298.257223563]],PRIMEM[&quot;Greenwich&quot;,0.0],UNIT[&quot;Degree&quot;,0.0174532925199433],AUTHORITY[&quot;EPSG&quot;,4326]]" 
    XOrigin = ET.SubElement(SpatialReference2, "XOrigin")
    XOrigin.text = "-400"
    YOrigin = ET.SubElement(SpatialReference2, "YOrigin")
    YOrigin.text = "-400"
    XYScale = ET.SubElement(SpatialReference2, "XYScale")
    XYScale.text = "999999999.99999988"
    ZOrigin = ET.SubElement(SpatialReference2, "ZOrigin")
    ZOrigin.text = "-100000"
    ZScale = ET.SubElement(SpatialReference2, "ZScale")
    ZScale.text = "10000"
    MOrigin = ET.SubElement(SpatialReference2, "MOrigin")
    MOrigin.text = "-100000"
    MScale = ET.SubElement(SpatialReference2, "MScale")
    MScale.text = "10000"
    XYTolerance = ET.SubElement(SpatialReference2, "XYTolerance")
    XYTolerance.text = "8.983152841195215e-09"
    ZTolerance = ET.SubElement(SpatialReference2, "ZTolerance")
    ZTolerance.text = "0.001"
    MTolerance = ET.SubElement(SpatialReference2, "MTolerance")
    MTolerance.text = "0.001"
    HighPrecision = ET.SubElement(SpatialReference2, "HighPrecision")
    HighPrecision.text = "true"
    LeftLongitude = ET.SubElement(SpatialReference2, "LeftLongitude")
    LeftLongitude.text = "-180"
    WKID = ET.SubElement(SpatialReference2, "WKID")
    WKID.text = "4326"
    LatestWKID = ET.SubElement(SpatialReference2, "LatestWKID")
    LatestWKID.text = "4326"
    root.append(SpatialReference2)
    #convert xml byte string to string 
    new_xml = str(ET.tostring(root), 'utf-8')
    #write xml over top old file
    with open(this_map_xml, "w") as xmlW:
        #xml 'wrangler' replaces stubborn escaping and quotation types
        xmlW.write(new_xml.replace("&amp;", "&").replace('"',"'"))


##Delete reference to aprx
del aprx

#rezip everything.   
shutil.make_archive(os.path.join(out_dir, output_project + "_Project", output_project), "zip", tempDirMod)


new_aprx_path_with_sr = os.path.join(out_dir, output_project + "_Project", output_project + "_withSR") + ".aprx"
# rename the .zip file to .aprx
os.rename(os.path.join(out_dir, output_project + "_Project", output_project + ".zip"), new_aprx_path_with_sr)

#end set spatial ref

‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍

* i pulled this from a larger project i'm working on and it hasn't been tested in the above state and may need tweeks.

Re: Modify map coordinate system with arcpy.mp?

RyanDavis1 — Wed, 19 May 2021 19:47:51 GMT

This is an old thread, but it led me to believe this wasn't possible for a long time. After looking over documentation, I was able to use arcpy to set the spatial reference on the map object.

Something like this works for me:

mapx.spatialReference = arcpy.SpatialReference(prjFile)

where mapx is a map object. Hope that helps someone.

Re: Modify map coordinate system with arcpy.mp?

DaveTaylor011 — Sun, 12 Dec 2021 03:04:15 GMT

Thanks Ryan, your reply to an old thread saved my bacon!

Re: Modify map coordinate system with arcpy.mp?

ArtK — Tue, 21 Dec 2021 20:26:45 GMT

Thanks for the answer Ryan. This is what my code looks like. Where DEM is the filepath for a layer you want your map to match CRS of.

aprx = arcpy.mp.ArcGISProject("CURRENT")
mapx = aprx.activeMap
mapx.spatialReference = arcpy.Describe(DEM).spatialReference