Korrektur eines DEM, Bearbeitung von Z-Werten, Z-Werte aus Attributen

Question

<\/HEAD>Ich habe versucht, eine Einzugsgebietsabgrenzung auf einem kleinen Stück eines 10m DEM im zentralen Bundesstaat Washington durchzuführen. Ich erhielt einige seltsame Ergebnisse und als ich das Akkumulationsraster untersuchte, stellte ich fest, dass eine Naht im DEM dazu führte, dass der gesamte Abfluss eines angrenzenden Einzugsgebiets von meinem Einzugsgebiet erfasst wurde. <\/SPAN> Ich machte mich daran, dieses DEM zu korrigieren. Ich dachte, ich würde einfach Raster zu Multipoint machen, eine Punkt-Feature-Class mit Z-Werten erstellen, die wie kleine Dämme und Kanäle in der Oberfläche wirken, ein neues Terrain aus den Multipunkten und der neuen Punkt-Feature-Class erstellen, dann Terrain zu Raster und das Einzugsgebiet erneut ausführen. <\/SPAN> Aber die Bearbeitung der Punkt-Feature-Class ist zu meinem Problem geworden. Als ich die Punkteschicht erstellte, habe ich sie so angelegt, dass sie Z-Werte enthält, und ich habe ein Attribut namens VALUES hinzugefügt. Beim Bearbeiten habe ich die gewünschte korrigierte Höhe in das VALUES-Attribut eingetragen. Die einzige Möglichkeit, die korrigierte Höhe während der Bearbeitung in die Z-Werte zu bekommen, war die Bearbeitung der Scheitelpunkte unter Sketch-Eigenschaften, und das konnte ich nur einzeln machen. Ich dachte, ich würde einfach "Populate Zs from attributes" verwenden. Aber nachdem ich so weit gekommen war, scheint "Populate Zs from attributes" nur über Production Z Management und Production Editing verfügbar zu sein, aber ich glaube, Production Editing ist nur in den Nautical- oder Aeronautical-Erweiterungen verfügbar, die ich nicht habe. Bei der Suche im Support fand ich <\/SPAN>VBA-Lösungen<\/A>, um Z-Werte aus Attributen zu aktualisieren, aber das scheint nur für 9.x anwendbar zu sein. Ich versuche herauszufinden, wie man das in 10.0, SP3 macht.<\/SPAN><\/BODY><\/HTML>

PaulHuffman · Answer

Thanks Richard.

I didn't figure out a way to easily change a group of vertexes, but I figured out how most people update DEMs.  I've been creating a new point or line feature class with a item called "value", then editing in new elevation values into this "value" item,  then running this feature class through the "Feature to Raster" tool using the "value" item as the field input,  then slipping the new cell values into the old DEM with a Con statement in Map Algebra.  In the python window the syntax is like:

demfix = Con(IsNull(Raster("newvalues")),"olddem","newvalues")

RichardFairhurst · Answer

I tried to run a watershed delimitation on a small piece of 10m DEM in Central Washington State.  I was getting some strange results and when I examined the accumulation grid, I found that a seam in the DEM was causing all the flow from an adjacent drainage to get captured by my drainage.

I set about trying to correct this DEM.  I figured I'd just go Raster to Multipoint,  create a point FC with Z values to act like little dams and channels in the surface, create a new terrain from the multipoints and the new point feature class, then terrain to raster, and rerun the watershed.

But editing the point feature class has become my problem.  When I created the points layer, I created it to contain Z-values, and I added an attribute called VALUES, and when editing, I stuck the desired corrected elevation into the VALUES attribute.  The only way I could figure out to get the corrected elevation into the Z-values while editing was by editing the vertices under Sketch Properties, and I could only do one at a time. I thought I would just "Populate Zs from attributes".  But after I got this far, it seems that "Populate Zs from attributes" is only available from Production Z Management and Production Editing,  but I think Production Editing is only available in the Nautical or Aeronautical extensions which I don't have.  Searching support I found VBA solutions to updating Z values from attributes, but that seems to be only applicable to 9.x.  I'm trying to figure out how to do this in 10.0, SP3.

Based on the ArcGIS 10 Field Calculator help you can adjust the coordinates of a point by using the field calculator on the Shape field of your point class.  Here is the example given:

"For a point feature class, shift the x coordinate of each point by 100.

Parser:
Python

Expression:
shiftXCoordinate(!SHAPE!)

Code Block:def shiftXCoordinate(shape):
   shiftValue = 100
   point = shape.getPart(0)
   point.X += shiftValue
   return point

"

Adapting this to what you need it should be something like this (where the !VALUES! field should match the name of the attribute field containing the Z values you want to assign to your points and no values of that field should be Null):

Parser:
Python

Expression:
assignZCoordinate(!SHAPE!, !VALUES!)

Code Block:def assignZCoordinate(shape, z_Value):
   point = shape.getPart(0)
   point.Z = z_Value
   return point

I hope this helps.