Latest Contributions by DanPatterson_Retired

http://ideas.arcgis.com/ post there as a requested feature, otherwise you will have to create your own in Python ... View more

Are both in the same coordinate system?  To check, add each one to a new dataframe and mouse around on the screen to check the coordinates in the bottom right (ie don't rely on the meta data) ... View more

the values are sorted first by the first input raster, then the second, that pair is assigned the first class...and so on ... View more

Try some of the other suggestions then, there have been several others but you appear to be fixed on one, is there a reason? ... View more

e input raster can have any cell size and must be a valid integer raster dataset. from the version 10 help files...check other versions if you aren't using ArcGIS 10 ... View more

convert the tif to an ESRI grid using the IPRIM field ... View more

use "combine" to produce the unique conditions if you don't want to do the nested "con" statements, then query for your cases.  It is all documented in the help files ... View more

Con http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#//009z00000005000000.htm ... View more

run the "repair geometry" tool on the file to check and correct any geometry errors ... View more

The following function requires a shape from your search cursor, the type of shape (ie multipoint, polygon or polyline and the arcpy object.  You can glean the necessary lines for your code. def shapeToPoints(a_shape,theType,arcpy):   '''   pnts = shapeToPoints(a_shape, shape type, geoprocessor)   Purpose:  Converts a shape to points, the shape and its type   are passed by the calling script   Requires:  def pntXY(pnt)   '''   outList=[]   part_num = 0   part_count = a_shape.partCount   if theType == "Multipoint":    #Multipoints     while part_num < part_count:       pnt = a_shape.getPart(part_num)       XY = pntXY(pnt)       if XY not in outList:         outList.append(XY)       part_num += 1   else:                          #Poly* features     while part_num < part_count: #cycle through the parts       a_part = a_shape.getPart(part_num)       pnt = a_part.next()       while pnt:                 #cycle through the points         XY = pntXY(pnt)         if XY not in outList:           outList.append(XY)         pnt = a_part.next()              if not pnt:              #null point check (rings/donuts)           pnt = a_part.next()           if pnt:             XY = pntXY(pnt)             if XY not in outList:               outList.append(XY)       part_num += 1   return outList ... View more

are you using point density or kernel density?  What output units did you specify? If none are chosen and you are using projected data (assume meters as unit) then your outputs will be relative to 1 m^2.  When you say the results are strange, what are they? ... View more

I would project all the data first prior to running all the point density analysis.  Once you are done, use the reclassify tool in the Spatial Analayst toolbox to reclass your data ... View more

Is there any way you could edit your post and enclose the code in a code (the # symbol in the html editor) to make it easier to read? ... View more

"doesn't seem to work" isn't providing enough information...are you getting errors? are the input files in the same projection?  more info needed ... View more

I am sure all of the ArcInfo license holders are ecstatic...but don't get too excited for the ArcView of ArcEditor holders...(please ESRI...qualify license levels) ... View more