topic Re: DTM editing in Data Management Questions

DTM editing

PeterMichalak — Fri, 11 Jul 2014 08:57:05 GMT

Dear,

I'm having problems editing my DTM (Tiff format). I need to assign a certain elevation for a river. Is it possible to draw a polygon covering the river, with a certain elevation value and somehow merge it with the DTM raster? Or is there an easier way to cope with it?

Thanks for your help

Peter

Re: DTM editing

808707 — Fri, 11 Jul 2014 11:54:15 GMT

Does the river change elevation or is it the same for the entire length?

Have a look at raster calculator.

You could buffer your river to the desired pixel width and elevation

Convert this to a raster then subtract the river raster from the DTM within Raster calculator.

Hope this points you in the right direction.

Re: DTM editing

PeterMichalak — Fri, 11 Jul 2014 12:08:01 GMT

Thanks for you answer 808 707

Yes the river has got a constant elevation. However a buffer would not be accurate enough, I need to set the river boundaries more accurately. The ideal way would be to draw a polygon representing the river set a 'z value' to the polygon and then convert it to raster and somehow superseed DTM, information for river areas with constant value.

Any help appreciated,

Cheers,

Peter

Re: DTM editing

IanMurray — Fri, 11 Jul 2014 13:43:27 GMT

Peter,

Have you looked at the ArcHydro downloadable extension? It should allow you to burn streams into your DTM, overriding the existing pixel values.

Arc Hydro Overview | ArcGIS Resource Center‌

Problems with Burn Stream slope and DEM Reconditioning‌

http://www.google.com/url?sa=t&rct=j&q=&esrc=s&frm=1&source=web&cd=1&ved=0CB0QFjAA&url=http%3A%2F%2Fdownloads.esri.com%2…

Hope this gets you in the right direction

‌

Re: DTM editing

PeterMichalak — Fri, 11 Jul 2014 15:05:56 GMT

Ian,

Thanks for your answer. I have installed ArcHydro and tried to use Level DEM tool but...it did not work

I have drawn a polygon around river and added a field representing its constant elevation. Then I used Level DEM tool but it did not changed the input DTM at all (or actually changed it in a way I cannot really understand).

Is there something I'm doing wrong here?

Cheers,

Peter

Re: DTM editing

DanPatterson_Retired — Fri, 11 Jul 2014 15:15:12 GMT

If you have converted the River to raster with an appropriate elevation simply use the Con statement

Con("River" >= 0, "River","DTM")

Assuming those are the names of the river and dtm respectively

Re: DTM editing

PeterMichalak — Mon, 14 Jul 2014 08:22:52 GMT

Dan,

Thanks for you answer, however that won't solve my problem. The Con statement you've proposed is based on cell value which is not the case here.

Let me explain:

I've drawn a polygon representing a river, converted it to a raster with constant value of 87 (approx. elevation of the river in local projection). I've also got a DTM model but due to water's reflectivity and other factors the river parts of my DTM are far from OK. All I want to do is replace the water bits of the original DTM with my river raster, so the output raster will have all the values of the original DTM, except for river which will have a value of 87 (constant elevation take from river raster)

Hope I made it clear

Cheers,

Peter

Re: DTM editing

DanPatterson_Retired — Mon, 14 Jul 2014 08:27:30 GMT

I guess I am not sure if follow. assuming that you want to keep the values of the river as 87 and classify all the nodata areas within that file to the values found in the DTM, then the con expression would be

Con("River" = 87, "River","DTM")

Otherwise, I have missed something simple in your explanation. Sorry...

Re: DTM editing

PeterMichalak — Mon, 14 Jul 2014 09:09:21 GMT

Dan,

Please see attached screenshoots representing DTM (values raging from 30 - 128), River (raster of constant value of 87) and DTM+River, hope that will help you understand the problem.

DTM

River

Output = DTM + River (values of DTM where DTM does not overlap River and River values (87) within River raster boundaries)

I have tried to use your syntax but got an error 000539: Runtime error SyntaxError: Keyword can't be an expression.

I'm not sure what have caused the problem. My guess it's the Syntax or the fact that River raster is strictly limited as shown on the above picture and doesn't have any nodata areas...

Still it's my guess, but I'm still unable to generate the output raster

Any suggestions much appreicated

Cheers,

Peter

Re: DTM editing

DanPatterson_Retired — Mon, 14 Jul 2014 09:18:17 GMT

Sorry Peter, trying to help without Arc* running.

Use the Con tool in arctoolbox. It may be a double equal sign ( == for the equality check), try not to type the expression, select your rasters, the equality check etc when constructing the Con statement. Also, in the Environments tab, set the combined extent of the rasters to be the union of the inputs (or that of the DTM), make sure that the cell size is the same. Save the resultant raster to a folder (temp or otherwise) for testing purposes.

If this doesn't work, I hope that someone with Arc* running will be able to wade in. Regards Dan

Re: DTM editing

PeterMichalak — Mon, 14 Jul 2014 13:57:48 GMT

Right,

I have managed to sort the problem out, if anyone has the same issue I'd suggest the following:

1) Draw a shp representing river > Convert shp to Raster with a predefined value of 87 (elevation)

2) Cut the DTM raster using Mask tool in Image Analysis Window

3) Clip 2 rasters using Mosaic to New Raster Tool

4) I've used SetNull tool to Set Null values for the background

Everyting worked fine except I've got small gaps where there 2 rasters meet (after using Mosaic to New Raster):

I'm not sure what have caused this, probably the fact that rasters were cut using polygon which represents irregular shp, comparing to raster's regular cell grid. Anyway I'm not sure how to get rid of these gaps, I have tried Snap to Raster setting but did not help.

Any suggestions?

Regards,

Peter

Re: DTM editing

JimHall — Tue, 15 Jul 2014 16:41:55 GMT

The problem may be in your not setting the "snap raster" for your outputs. Set environment setting to your original raster for all your steps.