Latest Contributions by ChrisHills

‎08-05-2016

I would like to copy features from a file geodatabase stored on a local volume to a file geodatabase on a networked share. If I copy the file geodatabase it transfers at around 12MB/s (on a 100baseT ethernet network). If I try to copy the feature layer from within ArcMap (10.2) it transfers at around 1KB/s. This is 12,000 times slower! Why is this and is there any way to speed it up, or should I resort to copying the entire geodatabase which is much faster?

‎07-09-2015

That's what I thought Jeff. It would be nice if ArcMap (and Server) could do this on the fly!

‎07-09-2015

Thanks for the suggestion Wes. Yes this is available, but there are tens of thousands of points so I don't really want to have to create the associated label classes. Ideally I would like the label placement routine to take into account the population field when deciding which labels to remove.

‎07-09-2015

I am using the maplex labelling system, and I would like to be able to specify a population field, so that when labels need to be removed, it starts by considering the features that have the lowest population first. Is this possible in ArcMap (>= 10.2) or would this be a feature request?

‎05-07-2015

Jai, thank you very much this is exactly what I wanted! Edit: Is it possible to do this without saving to a file? I will be doing this on thousands of coordinates so this is undesirable as it will place increased load on the disk, and slow down the process.

‎05-07-2015

As part of a script I need to convert coordinates from one SRS to another. This is the code I have so far, but I am missing the bit in the middle. I have done some extensive googling but I was not able to find a solution so I would appreciate some assistance. import arcpy inputSRS = 'Projected Coordinate Systems/National Grids/Europe/British National Grid' # Projected outputSRS = 'Geographic Coordinate Systems/World/WGS 1984' # Geographic trsOut = 'OSGB_1936_To_WGS_1984_Petroleum' srIn = arcpy.SpatialReference(inputSRS) srOut = arcpy.SpatialReference(outputSRS) inX = 210000 inY = 310000 # Code goes here print outX, outY # "-4.81042685256 52.6543521582"

‎03-27-2015

I too would like to be able to maintain the same extent for DDP. All that is needed is to add an additional option on the "Extent" tab which says "Keep current extent".

‎11-13-2014

This is a great idea. I would really like to be able to create dynamic fields for when the source dataset is read-only and I do not want to create a local copy which would be a point-in-time snapshot and would quickly become out of data as the source data is changed.

‎10-09-2014

I updated my script to include a spatial reference and it now gives the expected result:- import arcpy startx = 285000 starty = 138000 boxsize = 1536 boxcount = 2 features = [] # British National Grid srid = arcpy.SpatialReference(27700) for colx in range(1,boxcount+1): for coly in range(1,boxcount+1): top = (starty+(boxsize*coly)) bottom = (starty+(boxsize*(coly-1))) left = (startx+(boxsize*(colx-1))) right = (startx+(boxsize*colx)) feature_info = [[[left, bottom], [right, bottom], [right, top], [left, top]]] for feature in feature_info: features.append(arcpy.Polygon(arcpy.Array([arcpy.Point(*coords) for coords in feature]),srid)) print left, bottom, right, top for feature in features: print feature.extent 285000 138000 286536 139536 285000 139536 286536 141072 286536 138000 288072 139536 286536 139536 288072 141072 285000 138000 286536 139536 NaN NaN NaN NaN 285000 139536 286536 141072 NaN NaN NaN NaN 286536 138000 288072 139536 NaN NaN NaN NaN 286536 139536 288072 141072 NaN NaN NaN NaN

‎10-08-2014

Thank you Dan for your helpful answer.

‎10-08-2014

When creating polygons with Polygon() the results are not accurate as expected. Take the following script:- import arcpy startx = 285000 starty = 138000 boxsize = 1536 boxcount = 2 features = [] for colx in range(1,boxcount+1): for coly in range(1,boxcount+1): top = (starty+(boxsize*coly)) bottom = (starty+(boxsize*(coly-1))) left = (startx+(boxsize*(colx-1))) right = (startx+(boxsize*colx)) feature_info = [[[left, bottom], [right, bottom], [right, top], [left, top]]] for feature in feature_info: features.append(arcpy.Polygon(arcpy.Array([arcpy.Point(*coords) for coords in feature]))) print left, bottom, right, top for feature in features: print feature.extent This produces the following output:- 285000 138000 286536 139536 285000 139536 286536 141072 286536 138000 288072 139536 286536 139536 288072 141072 285000.00012207 138000.00012207 286536.00012207 139536.00012207 NaN NaN NaN NaN 285000.00012207 139536.00012207 286536.00012207 141072.00012207 NaN NaN NaN NaN 286536.00012207 138000.00012207 288072.00012207 139536.00012207 NaN NaN NaN NaN 286536.00012207 139536.00012207 288072.00012207 141072.00012207 NaN NaN NaN NaN As you can see, instead of the resulting boundary being equal to the input points, it has changed them by 0.00012207. This is unexpected and unwanted. How do I create my polygons without this error?

‎10-08-2014

Thanks for the helpful answer!

‎10-03-2014

When creating python scripts, I find that there items missing from the intellisense popups in the Python window. I also see that the same items are missing when examining an object with dir(). For example:- cursor = arcpy.SearchCursor("My Feature Class") row = cursor.next() dir(row) ['__class__', '__cmp__', '__delattr__', '__dict__', '__doc__', '__format__', '__getattr__', '__getattribute__', '__hash__', '__init__', '__module__', '__new__', '__passthrough_to_ao__', '__reduce__', '__reduce_ex__', '__repr__', '__setattr__', '__sizeof__', '__str__', '__subclasshook__', '__weakref__', '_arc_object', '_go', 'getValue', 'isNull', 'setNull', 'setValue'] dir(row.SHAPE) ['JSON', 'WKB', 'WKT', '__add__', '__class__', '__cmp__', '__delattr__', '__dict__', '__doc__', '__format__', '__from_scripting_arc_object__', '__geo_interface__', '__getattribute__', '__hash__', '__init__', '__iter__', '__module__', '__new__', '__or__', '__reduce__', '__reduce_ex__', '__repr__', '__setattr__', '__sizeof__', '__str__', '__sub__', '__subclasshook__', '__type_mapping__', '__type_string__', '__weakref__', '__xor__', '_arc_object', '_fromGeoJson', '_go', '_passthrough', 'area', 'boundary', 'buffer', 'centroid', 'clip', 'contains', 'convexHull', 'crosses', 'cut', 'difference', 'disjoint', 'distanceTo', 'equals', 'extent', 'firstPoint', 'getArea', 'getLength', 'getPart', 'hullRectangle', 'intersect', 'isMultipart', 'labelPoint', 'lastPoint', 'length', 'length3D', 'measureOnLine', 'overlaps', 'partCount', 'pointCount', 'positionAlongLine', 'projectAs', 'queryPointAndDistance', 'snapToLine', 'spatialReference', 'symmetricDifference', 'touches', 'trueCentroid', 'type', 'union', 'within'] Note that SHAPE is missing from the output of dir(row). In the intellisense popup it shows getValue, isNull, setNull, setValue, but not SHAPE. This makes it a bit more difficult when writing scripts as I have to remember that the object is missing.

‎08-13-2014

Thanks for the pointer Dan, the answer was to use measured projection instead of graticule.

‎08-13-2014

The coordinate system of my data frame is set to British National Grid, however when I try to add a grid (graticule) it is always projected using WGS84. How can I add a grid using the correct projection? ArcGIS version is 10.2.1.3497.

Online Status	Offline
Date Last Visited	‎02-12-2021 12:23 AM

My Ideas

Latest Contributions by ChrisHills

Why is copying features so slow compared to simply copying the file geodatabase?

Re: Removing labels based on population size

Re: Removing labels based on population size

Removing labels based on population size

Re: Coordinate transformation in arcpy

Coordinate transformation in arcpy

Re: Constant extent with Data Driven Pages

Re: Dynamic Fields for Shape Features and Statistics

Re: Accuracy loss when creating polygons with Polygon()

Re: Accuracy loss when creating polygons with Polygon()

Accuracy loss when creating polygons with Polygon()

Re: Why are some items missing from intellisense?

Why are some items missing from intellisense?

Re: Why does my grid always use wgs84 instead of the data frame projection?

Why does my grid always use wgs84 instead of the data frame projection?

Re: What version of WGS84 is in ArcGIS

Re: Unable to use ArcGIS Pro with a personal use l...

Cannot create new python environment in ArcGIS Pro...

Why does my grid always use wgs84 instead of the d...