idea Arcade projectAs Geometry Function in ArcGIS Pro Ideas

Arcade projectAs Geometry Function

JacobFougstrupNicolajsen — Fri, 06 May 2022 20:18:37 GMT

Working primarily with data in UTM projection, I often need ways to convert UTM coordinates to WGS to create deep links to various online services.

Having to create duplicate datasets, or running python scripts as regular intervals is tedious and time consuming.

And arcade equivalent of the python projectAs method would solve almost all these issues.

Re: Arcade projectAs Geometry Function - Status changed to: Under Consideration

CraigWilliams — Wed, 11 May 2022 23:45:49 GMT

We've been discussing this on the development team. If others have different use cases, please feel free to add them to the comments.

Re: Arcade projectAs Geometry Function

rsun_TQB — Thu, 12 May 2022 01:42:32 GMT

I can't agree more with this. As the largest GIS software company on earth, I was surprised that Arcade doesn't have function like proj4.js built-in.

The solution here works, but still hope to have something more built-in 🙂

Solved: Lat/Long unit conversion with Arcade - Esri Community

Reno

Re: Arcade projectAs Geometry Function

KleinschmidtBrad — Fri, 24 Jun 2022 04:13:15 GMT

I'm looking for similar functionality, my use case is to auto calculate areas in hectares or length in metres for features in a geographic coordinate system.

Re: Arcade projectAs Geometry Function

RichardHowe — Wed, 14 Sep 2022 12:51:56 GMT

Being able to have a form in a webmap automatically calculate x and y coordinates in local projections systems upon data capture/data edit via arcade would be hugely advantageous. Having to manually calculate these fields leads to a danger of them not being updated when data is edited.

Add PROJ4 support to ArcGIS Arcade

TomNeer — Mon, 30 Jan 2023 18:51:33 GMT

It would be extremely helpful to have access to a projection library (e.g., PROJ4JS | Proj4js) in ArcGIS Arcade to support the calculation of coordinates in differing systems (e.g., USNG and Latitude\Longitude). Arcade will only return the XY coordinates in the projection of the dataset.

The only solution is to translate the spherical trig into Arcade

// Calculates the latitude of a State Plane coordinate // Get Colorado State Plane North (2231) coordinates from geometry var x = Geometry($feature).x var y = Geometry($feature).y // Projection Information var linear_unit = 0.3048006096012192 var false_easting = 3000000.000316083 var false_northing = 999999.999996 var sp1 = 39.71666666666667 var sp2 = 40.78333333333333 var latitude_of_origin = 39.33333333333334 var semimajor_axis = 6378137.0 var semiminor_axis = 6356752.314140356 function degToRad(deg) { return (deg * PI / 180) } function radToDeg(rad) { return (rad * 180) / PI } function tVal(t) { return (Tan((PI / 4) - (t / 2)) / Pow((1 - (ecc * Sin(t))) / (1 + (ecc * Sin(t))), (ecc / 2))) } // Convert projection information into radians var lorR = degToRad(latitude_of_origin) var sp1R = degToRad(sp1) var sp2R = degToRad(sp2) // Spherical Trig voodoo var flat = 1 / (semimajor_axis / (semimajor_axis - semiminor_axis)) var ecc = Sqrt(2 * flat - Pow(flat, 2)) var n = (Log(Cos(sp1R) / Pow((1 - Pow(ecc, 2) * Pow(SIN(sp1R), 2)), 0.5)) - Log(Cos(sp2R) / Pow((1 - Pow(ecc, 2) * Pow(SIN(sp2R), 2)), 0.5))) / (Log(tVal(sp1R)) - Log(tVal(sp2R))) var f = (Cos(sp1R) / Pow((1 - Pow(ecc, 2) * Pow(SIN(sp1R), 2)), 0.5)) / (n * Pow(tVal(sp1R), n)) var rf = (semimajor_axis * (1 / linear_unit)) * f * Pow(tVal(lorR), n) var rz = Pow(Pow((x - false_easting), 2) + Pow((rf - (y - false_northing)), 2), 0.5) var tz = Pow((rz / ((semimajor_axis * (1 / linear_unit)) * f)), (1 / n)) // Latitude iterations in radians var latTR = (PI / 2) - (2 * Atan(tz)) var lat1R = (PI / 2) - (2 * Atan((tz * Pow((1 - (ecc * Sin(latTR))) / (1 + (ecc * Sin(latTR))), (ecc / 2))))) var lat2R = (PI / 2) - (2 * Atan((tz * Pow((1 - (ecc * Sin(lat1R))) / (1 + (ecc * Sin(lat1R))), (ecc / 2))))) var lat3R = (PI / 2) - (2 * Atan((tz * Pow((1 - (ecc * Sin(lat2R))) / (1 + (ecc * Sin(lat2R))), (ecc / 2))))) var lat4R = (PI / 2) - (2 * Atan((tz * Pow((1 - (ecc * Sin(lat3R))) / (1 + (ecc * Sin(lat3R))), (ecc / 2))))) var latFR = (PI / 2) - (2 * Atan((tz * Pow((1 - (ecc * Sin(lat4R))) / (1 + (ecc * Sin(lat4R))), (ecc / 2))))) // Latitude values in degrees var latitude = radToDeg(latFR) return latitude;

Re: Arcade projectAs Geometry Function

MattRingel — Mon, 27 Feb 2023 23:40:00 GMT

+1 from me. This would be useful for data interoperability.

Re: Arcade projectAs Geometry Function

CarolineLindhe — Tue, 28 Feb 2023 13:47:11 GMT

+1 from me! My whole organization is longing for this!

Re: Arcade projectAs Geometry Function

JasonClemis — Thu, 09 Mar 2023 16:19:05 GMT

I also have a need for this based on a recent request to convert from State Plane to Decimal Degrees in an attribute rule

Re: Arcade projectAs Geometry Function

DavidSolari — Thu, 16 Mar 2023 17:53:35 GMT

My team is constantly asked to show a point's lat/lon as a standard attribute, this idea would simplify those requests and ease maintenance.

Re: Arcade projectAs Geometry Function

DavidColey — Wed, 12 Apr 2023 20:33:02 GMT

This is great @TomNeer thanks for posting this - Forgive me for asking, but are you able to point me in the right direction to get build the same calculation for a state plane coordinate translation to longitude?

Thanks--

Re: Arcade projectAs Geometry Function

TomNeer — Thu, 13 Apr 2023 21:03:56 GMT

@DavidColey - Here is the Arcade expression that can return either latitude or longitude. So at first, I thought my code was **bleep** because it was returning slightly different values than Esri's Calculate Geometry tool (even double-checking my transformations), but I ran several test SPC through NGS Coordinate Conversion and Transformation Tool (NCAT) (noaa.gov) and this code is more accurate than Esri's. Go figure...

// Calculates the latitude or longitude of a State Plane coordinate // Source: https://community.esri.com/t5/coordinate-reference-systems-questions/formula-for-state-plane-to-lat-lon-conversion/td-p/870543/page/2 // Get Colorado State Plane North (2231) coordinates from geometry var x = Geometry($feature).x; var y = Geometry($feature).y; // Projection Information (this must be changed if you are not using EPSG 2231) var linear_unit = 0.3048006096012192; var false_easting = 3000000.000316083; var false_northing = 999999.999996; var central_meridian = -105.5; var sp1 = 39.71666666666667; var sp2 = 40.78333333333333; var latitude_of_origin = 39.33333333333334; var semimajor_axis = 6378137.0; var semiminor_axis = 6356752.314140356; function degToRad(deg) { return (deg * PI / 180) } function radToDeg(rad) { return (rad * 180) / PI } function tVal(t) { return (Tan((PI / 4) - (t / 2)) / Pow((1 - (ecc * Sin(t))) / (1 + (ecc * Sin(t))), (ecc / 2))) } // Convert projection information into radians var lorR = degToRad(latitude_of_origin); var sp1R = degToRad(sp1); var sp2R = degToRad(sp2); // Spherical Trig voodoo var flat = 1 / (semimajor_axis / (semimajor_axis - semiminor_axis)); var ecc = Sqrt(2 * flat - Pow(flat, 2)); var n = (Log(Cos(sp1R) / Pow((1 - Pow(ecc, 2) * Pow(SIN(sp1R), 2)), 0.5)) - Log(Cos(sp2R) / Pow((1 - Pow(ecc, 2) * Pow(SIN(sp2R), 2)), 0.5))) / (Log(tVal(sp1R)) - Log(tVal(sp2R))) var f = (Cos(sp1R) / Pow((1 - Pow(ecc, 2) * Pow(SIN(sp1R), 2)), 0.5)) / (n * Pow(tVal(sp1R), n)) var rf = (semimajor_axis * (1 / linear_unit)) * f * Pow(tVal(lorR), n) var rz = Pow(Pow((x - false_easting), 2) + Pow((rf - (y - false_northing)), 2), 0.5) var tz = Pow((rz / ((semimajor_axis * (1 / linear_unit)) * f)), (1 / n)) var zz = Atan((x - false_easting) / (rf - (y - false_northing))); // Latitude iterations in radians var latTR = (PI / 2) - (2 * Atan(tz)) var lat1R = (PI / 2) - (2 * Atan((tz * Pow((1 - (ecc * Sin(latTR))) / (1 + (ecc * Sin(latTR))), (ecc / 2))))) var lat2R = (PI / 2) - (2 * Atan((tz * Pow((1 - (ecc * Sin(lat1R))) / (1 + (ecc * Sin(lat1R))), (ecc / 2))))) var lat3R = (PI / 2) - (2 * Atan((tz * Pow((1 - (ecc * Sin(lat2R))) / (1 + (ecc * Sin(lat2R))), (ecc / 2))))) var lat4R = (PI / 2) - (2 * Atan((tz * Pow((1 - (ecc * Sin(lat3R))) / (1 + (ecc * Sin(lat3R))), (ecc / 2))))) var latFR = (PI / 2) - (2 * Atan((tz * Pow((1 - (ecc * Sin(lat4R))) / (1 + (ecc * Sin(lat4R))), (ecc / 2))))) // Longitude calculations in radians var lonR = ((zz / n) + degToRad(central_meridian)); // Latitude and Longitude values in degrees var latitude = radToDeg(latFR) var longitude = radToDeg(lonR) // Change me return longitude;

Re: Arcade projectAs Geometry Function

DavidColey — Thu, 13 Apr 2023 21:46:41 GMT

Way to go @TomNeer - I am starting to recall from my geomatics days that the longitude conversion is a simpler calc. Thanks for posting this.... I'll parameterize for Florida State Plane West wkid 2882 and let you know how it goes....

Re: Arcade projectAs Geometry Function

DavidColey — Fri, 14 Apr 2023 21:36:51 GMT

Alas I thought this might be a more-straightforward task, but as State Plane Coordinate systems using a Transverse Mercator projection do not contain parameters for standard parallels (since they're not conformal projections and the false northing is always 0), it is not a simple adjustment to @TomNeer 's Arcade above....

Re: Arcade projectAs Geometry Function

TomNeer — Fri, 14 Apr 2023 21:58:42 GMT

That is a limitation of this formula. It is possible, but my spherical
trigonometry is pretty rusty.

Would be easier if Esri just supported projections in Arcade.

Re: Arcade projectAs Geometry Function

DavidColey — Fri, 21 Apr 2023 16:30:35 GMT

Agreed. I'd have to dig back into old geomatics class work to find my spherical trig, and it has been a while that is for sure. But, for those out there working in any east-west extent state plane coord systems, your algorithm works well.

Only need to make the adjustments as you indicate.

In the meantime @TomNeer , I'd say that you have done half of Esri's work for them (as far as state plane conversions go.

Re: Arcade projectAs Geometry Function

MiltonSolano — Sat, 20 May 2023 17:42:34 GMT

Team,

+1 here! We capture new trees on the field using Field Maps. Because our trees are 1-2 meters apart, we use our local UTM Zone 17 North projection and to better improve tree geolocation, we developed a local grid reference system.

Having the option to calculate UTM positions directly on Field Maps will be a great functionality, because right now we need to update the UTM fields using ArcGIS Pro or AGOL notebook.

Best,

Mjss25

Re: Arcade projectAs Geometry Function

SeanLyons — Tue, 27 Jun 2023 20:51:28 GMT

Another use case for a built in projection is trying to deal with the difference between the control point actual location and the control point geometry location for the Pro Parcel fabric. This was fine in ArcGIS Desktop as point and the control were separate layers now that they are the same layer we need to show the correct absolute location of the control.

The snag is the fabric is on the fly projected to a UTM Zone but is stored in Albers. This means we can only get the labelling to be correct if it is in Albers. The instant we change the label to albers we get a distance that is giant. Example below.

Albers

UTM Zone

This label was made with some simple arcade for reference.

var point = Geometry($feature) var x1 = point.x; var y1 = point.y; var x2 = $feature.X; var y2 = $feature.Y; var offset = round(sqrt(pow(x2-x1,2)+pow(y2-y1,2)),2); if (offset > 0) return "(" + offset + ")";

Re: Arcade projectAs Geometry Function

shodge_fsu — Tue, 19 Sep 2023 12:51:17 GMT

Seems like a no brainer to add this functionality.

My users want to see state plane coordinates. Florida has (at least for now) three state plane zones and we use the Web Mercator projection for all of the layers in the map. The strategy for obtaining the Lat / Lon from these coordinates mentioned elsewhere in this thread works great, but now we need to project those coordinates to the appropriate state plane zone.

Hope to see this functionality in 1.13!

Re: Arcade projectAs Geometry Function

Adam_Coburn_NCC — Thu, 21 Sep 2023 11:45:20 GMT

Another +1

This functionality would make such a difference!

Being able to use an out of the box Arcade function to go from 27700 - OSGB36 (BNG) to WGS84 (4236) would allow us to pre-populate Geopoint questions in Survey123 via custom URL parameters.