Urban Wind Field: URock v2023a - Missing Data in Raster File

[AlexandraStankulova]

I am performing the wind analysis using the "Urban Wind Field: URock v2023a". The area is quite large and complex, but after 3hours the analysis finishes with success, loading as always the created raster file with the wind speed.
The problem is that the raster resullts with no data, while the vector file has the vaues in it.
What could be the problem?

Here I upload the input data and the output data created by UMEP:

https://we.tl/t-TZrhbqXpdl

[biglimp]

I have the same issue, but only on my windows machine. @AlexandraStankulova , are you on windows also?

I am simulating on 4 meter resolution, both in horizontal and vertical dimensions, then I have a 2 meter raster grid that I want to use for raster output.

I have QGIS 3.40, UMEP 2.1.3.

@j3r3m1 , how do you interpolate?

Below. you find my data and settings:

processing.run("umep:Urban Wind Field: URock", {'BUILDINGS':'D:/GRUPAS/THERMAL_COMFORT_FIVE_CITIES/MALMO/UROCK/UROCKIN/build_urock44.shp','HEIGHT_FIELD_BUILD':'ROOF_H','VEGETATION':'D:/GRUPAS/THERMAL_COMFORT_FIVE_CITIES/MALMO/UROCK/UROCKIN/veg_urock44.shp','VEGETATION_CROWN_TOP_HEIGHT':'VEG_H','VEGETATION_CROWN_BASE_HEIGHT':'trunk','ATTENUATION_FIELD':'','INPUT_PROFILE_FILE':'','INPUT_PROFILE_TYPE':0,'INPUT_WIND_HEIGHT':10,'INPUT_WIND_SPEED':2,'INPUT_WIND_DIRECTION':244,'RASTER_OUTPUT':'D:/GRUPAS/THERMAL_COMFORT_FIVE_CITIES/MALMO/UROCK/UROCKIN/temp_ClippedDSM.tif','HORIZONTAL_RESOLUTION':2,'VERTICAL_RESOLUTION':2,'WIND_HEIGHT':'1.5','UROCK_OUTPUT':'C:\\Users\\xlinfr\\Desktop\\Rockle_zones','OUTPUT_FILENAME':'urock_output','SAVE_RASTER':True,'SAVE_VECTOR':False,'SAVE_NETCDF':False,'LOAD_OUTPUT':True})

ResolutionErrorURock.zip

[AlexandraStankulova]

Yes, I also work with Windows, I have Windows 10.

[j3r3m1]

The interpolation depends on the resolution of your calculation and the resolution of the raster used as output constraint. However, most of the time (and this is the case for your data) this function is used:
https://github.com/UMEP-dev/UMEP-processing/blob/65b105bb04326870f306a693e5414c30bcd6baf2/functions/URock/saveData.py#L482

And then at the end the no data values are filled using the few lines after this one:
https://github.com/UMEP-dev/UMEP-processing/blob/65b105bb04326870f306a693e5414c30bcd6baf2/functions/URock/saveData.py#L452

I have not faced any issue with Ubuntu for your case @biglimp and I will try the @AlexandraStankulova case. The issue probably comes from one or the other of the code provided above.

I now have a Windows OS at home. It may help to fix Windows issues but I need to find time for that.

[j3r3m1]

@AlexandraStankulova can you please share a smaller file ? Your is more than 1Gb.

[AlexandraStankulova]

Yes sure. I will upload now only the input data and the output raster windspeed file that actually results with no data.

Here the files:

InputData.zip

urock_outputWS.zip

[j3r3m1]

Interesting, I could reproduce the error. I will investigate what is the issue

[j3r3m1]

Well the issue was that I did not use the right DSM for interpolating the results at the end... Thus I cannot reproduce the error on Ubuntu. I have to debug on a Windows machine, not sure when.

[lgzdbf]

I am encountering the same error on my side regarding outputs without data when running the process on Ubuntu.

@j3r3m1 Here is a sample to reproduce the error.

urock_grenoble_test.zip

Let me know if you need any additional information.

Thanks in advance for your help!

[j3r3m1]

OK there was quite a big problem. When the vegetation was not intersected at all by building wake zones (which is often the case in rural area obviously...), there vegetation was removed from the analysis... It is now solved.

[j3r3m1]

UMEP-dev/UMEP-processing@af8b178

[j3r3m1]

@lgzdbf it should solve your issue.
@AlexandraStankulova or @biglimp is this fix also solves your issue ? I guess not since I do not think it is related to Windows issue but we never know...

[AlexandraStankulova]

Hello. I guess it will solve it definetely because I was doing some different tests, and I noticed that sometimes when I increase the buffer of the data to be inserted also around the area of analysis, sometimes this problem was not appearing (but I coulnd't explain nmyself why with e certain buffer are it works and with another buffer area it doesnt).

[j3r3m1]

OK the TIN interpolation does not work in some specific cases. I have reported the bug in the QGIS repo. We will see how they can help (qgis/QGIS#60416).

Two alternatives can be used:

• using the IDW interpolation seems to work but not sure it can solve the problem in all situations (at least it solves it for the dataset provided by @lgzdbf) and I have set an arbitrary distance coefficient (I have tried several and this one results in something I feel comfortable with).
• using scipy interpolation methods and writing the raster using rasterio results in slightly faster writing but also a sligth shifting of the result and this means adding a new dependency with rasterio which seems not to be a good idea in terms of long-term maintenance.

@biglimp what do you think ?