Hello,
I'm running a series of simulations for different urban areas. When I do it for a small area that requires only 1 DEM file from de NASA ASTER collection, it works perfectly.
However, when my area is bigger than the ASTER DEM file, I have to merge those files, but when doing that, SUMO doesn't read the .tif file correctly. In fact I have the following error message:
Error: Unknown data type in Madrid_elevation.tif
where Madrid_elevation.tif has been created with the merge tool of the raster package in R from 9 DEM files. This is the full R code:
library(raster)
h1=raster("ASTGTMV003_N39W003_dem.tif")
h2=raster("ASTGTMV003_N39W004_dem.tif")
h3=raster("ASTGTMV003_N39W005_dem.tif")
h4=raster("ASTGTMV003_N40W003_dem.tif")
h5=raster("ASTGTMV003_N40W004_dem.tif")
h6=raster("ASTGTMV003_N40W005_dem.tif")
h7=raster("ASTGTMV003_N41W003_dem.tif")
h8=raster("ASTGTMV003_N41W004_dem.tif")
h9=raster("ASTGTMV003_N41W005_dem.tif")
mmerge = merge(h1, h2, h3, h4, h5, h6, h7, h8, h9)
tmp <- "D:/ALEJANDRO-TESIS/3_sumo/Topografia/Para_unir/Madrid/"
plot(mmerge)
rf <- writeRaster(mmerge, filename=file.path(tmp, "Madrid_elevation.tif"), format="GTiff", overwrite=TRUE)
The only difference between the original DEM files and the merged one is that originals are 16bit and the merged one are 32bit colors. And there is no problem reading the merged file with any GIS software or image visor.
Any clue of what it is wrong?
Best regards,
Alejandro
--
Alejandro Rodríguez Sánchez
PhD student
CIEMAT // Universidad Complutense de Madrid
28040 Madrid
email:
alrodr25@xxxxxx