New answers tagged

1

The $$RowMap, $$ColMap, $$XMap and $$YMap variables (and some others) are no longer (directly) supported as of ArcGIS 10.0. You can use it via python: arcpy.env.extent = arcpy.Extent(0, 0, 10, 10) #Change to suit # Or arcpy.env.extent = "path to raster" arcpy.env.cellSize = 1 #Change to suit arcpy.gp.SingleOutputMapAlgebra_sa("$$ROWMAP + 1"), "...


3

If you create an empty tile cache, then use the Manage Tile Cache tool, you should be able to check off and generate only the tiles in your scales of interest before executing the tool.


0

A key to one aspect of this is determining the proper name for the non-user-defined fields of object id and geometry. The type of the geometry field is Double, which is unhelpful in this case. Using the describe function, one can determine the proper name for these fields across file types (i.e. shapefile v file gdb, etc; alleviating a lot of grief as oid ...


0

To have ArcGIS for Desktop, or ArcGIS Pro, generate a GeoTIFF that meets your requirements, I suspect that you will need to submit an ArcGIS Idea.


0

We have found that when using the X,Y export that the points were always on the street center line, but when using DisplayX/DisplayY, they (as documented) are in the parcel (for PointAddress matches, anyway). Then it's simply a matter of matching your parcel geometry to the lat/long point.


0

This is a documented bug: Bug BUG-000083280 Using in_memory workspaces causes the PATHCOST field to result in zero values.


0

I think you are describing features of Web Coverage Service (WCS) so maybe see if this data is being served that way too. Using WCS with ArcMap is described at Adding a WCS service to ArcMap - it's just like adding a WMS but gives access to grid cells themselves rather than just a rendered picture of them.


0

The Technical Article entitled Error: The coordinate system is invalid says that: An error message, "The coordinate system is invalid," is returned when accessing an ArcSDE version 9.3.1 feature class from a later ArcGIS version.


0

AnnotationProperties is a get/put property meaning a reference is not handed out. You need to set the collection back into the FeatureLayer after modifying it.


1

You can do this geoprocessing in QGIS by first making the image bounday of all DTED. This will convert all DTED to polygon shapefiles. Then run the spatial query for touches/intersects between your AOI and DTED polygon. You need to download the Image boundary plugin for this.


0

The function you're using to auto increment IDs is meant for use in the field calculator; in the field calculator you have to define the global variable to be used by other rows in the table, otherwise all your values would be 1. For your purposes, you necessarily need to define rec as a global variable, and you don't need to create a function to add 1 (...


0

If I understand your question, you want the same result you get from PointToRaster, but you'd prefer to have zeroes in all the NoData cells? Simply follow PointToRaster with Spatial Analyst > Reclassify, converting the NoDatas to zero and the data values to whatever makes sense for your project. Then set the symbology to get the visual contrast you are ...


1

Because a raster is continuous and points are discrete, you will need to use an interpolation method to create the kind of raster you would like. There are several different interpolation methods that may be better or worse for a specific application - but, here is a pretty brief overview of interpolation: http://pro.arcgis.com/en/pro-app/tool-reference/3d-...


1

If I get you correctly, You need to use the below steps- Use point to raster, it will give you a raster of points but blank where no points(i.e. NoData area rendered with nocolor)- maintain same pixel size with layer 2 e.g. 5*5 To get the background, convert polygon of UK into raster after setting the value for this layer 0- maintain same pixel size with ...


0

You should also check that you are not in an editing session.


-2

Have you tried: try: from osgeo import ogr except ImportError: import ogr


0

As a comment you provided a line of ArcPy code that I just tested using: import arcpy arcpy.CreateFishnet_management("C:/Temp/test.shp","223122.902847 4244274.58422","266744.967575 4410255.40647","0","0","6","6","447614.673105 4372719.14612","LABELS","wrs_dec_va_proj15033","POLYLINE") This reported an error: Traceback (most recent call last): File "C:\...


2

To dissolve down to Single-Part features your input features need to intersect (touching edges or overlapping), even when dissolving on an attribute. If the features don't intersect at all then the output will be the same as the input. If your features do not intersect you will need to set your dissolve to output Multi-Part features. This will dissolve ...


1

I'm not sure if you're asking to just create mile markers at a set distance or if you want to create the points from a distance attribute. My answer covers both methods: To use the distance field from your shapefile, export your shapefile to csv then use 'Make Route Event Layer' and point to the csv. If you only want points, select points - you'll have to ...


0

I don't know of any other way to do this unless you could access the image's EXIF data. Test "datetime" or "datetimeoriginal" and check the month the image was taken. The Python image library (PIL) module can do this like below. import arcpy, os, glob, re from arcpy import env from arcpy.sa import * from PIL import Image arcpy.env.overwriteOutput = True ...


0

I'm not familiar with the MODIS file naming conventions, but maybe you could try the regular expression module (re) to add rasters matching a pattern in your root folder, add them to "dekad" based on the pattern, then perform CellStatistics. Below, I've assumed the month is the 3 digits before the "Fpar_1km" (??) You could adjust it to suit your needs. ...


1

It sounds like you just need to iterate over each month's folder(?), add each month's tif files to the list "dekad", perform the CellStatistics command and save result as the monthly average raster, correct? Instead of the "for i in range... dekad.append[]..." method, you could do something like below using "glob"... import arcpy, os, glob from arcpy ...


1

The alternative to doing this in the raster calculator is to use the Cell Statistics tool in ArcToolbox. You can select all of your rasters at once and then pick which ever summary statistic you are interested in. It saves a lot of typing!


0

In Raster Calculator: ( "raster_1"+"raster_2"+...+"raster_n" )/n Where n is the number of rasters.


0

I created distances rasters to all 3 types of water bodies (cell size = 1000m*1000m) and used cell statistics (min) to find nearest distance. After that I ran zonal statistics as table and joined results back to polygons. This is how max distance looks like: What you call nuts is a set overlapping polygons, make sure to remove all with status < 3.


1

In response to @crld I have done some tests, Im not sure if im doing the timings right but, the first method doesnt actually work, it breaks out of the loop after 1 iteration import arcpy import itertools fc = r'D:\randomstuff\test2.gdb\Point0' #5146 points fc1 = r'D:\randomstuff\test2.gdb\points2' #484 points def printCount(start,count): end = time....


2

After looking at the code and the error you're getting, I think this can be pinned down to exporting the DDP PDFs. Python is trying to parse the outLoc path and failing because the backslashes are escaping and the full path isn't being read. The way to fix this is to us the os.path module. Insert the following in the top of your code: import os Then, ...


2

Both the HDF-EOS To GeoTIFF Conversion Tool (HEG) and the MODIS Reprojection Tool (MRT) can do that. You can download them from the USGS website: https://lpdaac.usgs.gov/tools/data_manipulation


4

Logically, they do effectively the same thing. For every row in one data set, loop through every line in one dataset (presumably until some criteria is satisfied, then break). Without a break on an if statement, if both datasets are 10,000 rows long, you would have to iterate through 100,000,000 rows. EDIT: However, the With With For For methodology doesn'...


0

As it looks like you're using MapInfo, here is a MapBasic function I wrote a while back for an in-house tool I was working on. It takes a source node (your centroid point) and a region object (the polygon) as arguments and returns a point object at the furthest node in the polygon from the source point. Function GetFurthest(ByVal oNode1 as Object, ByVal ...


2

I run the process on a server where both the 64-bits and 32-bits versions of Python are installed, as the server holds both ArcGIS Desktop and ArcGIS Server. By default, the 64-bit version is launched, and my script fails. I forced the 32-bit version use by launching the script from the command line, and the update went through successfully. C:\>C:...


1

I don't believe it is possible to do it this way, and the reason is because service layers do not support workspace paths. You can test this on your layer by lyr.supports("WORKSPACEPATH"), which should return False. However you could use UpdateLayer to replace the old service layer with a new one you create with the correct URL. I'm pretty sure to do it ...



Top 50 recent answers are included