Errors with Urban Wind Field: URock v2023a

[Lalalashp]

Hello everybody!
I've been having issues with the Urban Wind Field: URock v2023a tool.
I tried to reduce the size of the area but it still gave me the same error, which is:

Versione di QGIS: 3.34.5-Prizren
Revisione codice QGIS: 4b308492
Versione di Qt: 5.15.3
Versione Python: 3.9.18
Versione di GDAL: 3.8.4
Versione di GEOS: 3.12.1-CAPI-1.18.1
Versione di PROJ: Rel. 9.3.1, December 1st, 2023
versione di PDAL: 2.6.0 (git-version: 3fced5)
Algoritmo iniziato alle: 2025-01-12T01:28:40
In avvio l'algoritmo 'Urban Wind Field: URock v2023a' …
Parametri in ingresso:
{ 'ATTENUATION_FIELD' : '', 'BUILDINGS' : 'C:/Users/Laura/Dropbox (Politecnico Di Torino Studenti)/Tesi_magistrale/UMEP - Torino/AnotherOne/19/1_UrockPrepare/19_UrP_build.shp', 'HEIGHT_FIELD_BUILD' : 'Bui_HEIGHT', 'HORIZONTAL_RESOLUTION' : 5, 'INPUT_PROFILE_FILE' : '', 'INPUT_PROFILE_TYPE' : 0, 'INPUT_WIND_DIRECTION' : 45, 'INPUT_WIND_HEIGHT' : 10, 'INPUT_WIND_SPEED' : 2, 'LOAD_OUTPUT' : True, 'OUTPUT_FILENAME' : '19', 'RASTER_OUTPUT' : None, 'SAVE_NETCDF' : False, 'SAVE_RASTER' : False, 'SAVE_VECTOR' : False, 'UROCK_OUTPUT' : 'C:\Users\Laura\Desktop\19_UrbanWindField', 'VEGETATION' : 'C:/Users/Laura/Dropbox (Politecnico Di Torino Studenti)/Tesi_magistrale/UMEP - Torino/AnotherOne/Urock_veg_clip/_19estensione_32632_UrP_veget.shp', 'VEGETATION_CROWN_BASE_HEIGHT' : '', 'VEGETATION_CROWN_TOP_HEIGHT' : 'Veg_HEIGHT', 'VERTICAL_RESOLUTION' : 5, 'WIND_HEIGHT' : '1.5' }

Writing settings for this model run to specified output folder (Filename: RunInfoURock_YYYY_DOY_HHMM.txt)
Initiating algorithm
Creates an H2GIS Instance and load data
Creates the stacked blocks used as obstacles
Rotates obstacles to the right direction and calculates geometry properties
Roughness zone properties are:

• z0: 2.3880850405465655
• d: 11.374281137024472
• Hr: 15.920566936977105
• H_ob_max: 39.0
• lambda_f: 0.19125558056004152
  Creates the 2D Röckle zones
  Creates the 2D grid
  Initializes the 3D grid within Röckle zones
  Traceback (most recent call last):
  File "JdbcPreparedStatement.java", line 262, in org.h2.jdbc.JdbcPreparedStatement.execute
  Exception: Java Exception

The above exception was the direct cause of the following exception:

Traceback (most recent call last):
File "C:\Users\Laura\AppData\Roaming\Python\Python39\site-packages\jaydebeapi_init_.py", line 534, in execute
is_rs = self._prep.execute()
org.h2.jdbc.org.h2.jdbc.JdbcSQLDataException: org.h2.jdbc.JdbcSQLDataException: Invalid value "-0.06281183352066441" for parameter "LN() argument"; SQL statement:
CREATE TABLE VEGETATION_OPEN_TEMPO_3DPOINTS
AS SELECT b.ID_POINT, a.ID_Z, CASE WHEN a.Z>b.MAX_HEIGHT
THEN LOG((a.Z - 3 * 0.05 * MAX_CANOPY_HEIGHT) / 2.3880850405465655) / LOG(a.Z / 2.3880850405465655)
ELSE CASE WHEN a.Z > b.MAX_HEIGHT OR a.Z < b.MIN_HEIGHT
THEN LOG((b.MAX_CANOPY_HEIGHT - 3 * 0.05 * MAX_CANOPY_HEIGHT) / 2.3880850405465655) / LOG(a.Z / 2.3880850405465655)
ELSE LOG((b.MAX_HEIGHT - 3 * 0.05 * MAX_CANOPY_HEIGHT) / 2.3880850405465655) / LOG(a.Z / 2.3880850405465655) * EXP(b.ATTENUATIO * (a.Z / b.MAX_CANOPY_HEIGHT - 1))
END
END
AS VEGETATION_FACTOR
FROM Z_VALUES_20250112013951 AS a, VEGETATION_OPEN_POINTS AS b
WHERE a.Z > 0 [90008-232]

During handling of the above exception, another exception occurred:

Traceback (most recent call last):
File "C:\Users/Laura/AppData/Roaming/QGIS/QGIS3\profiles\default/python/plugins\processing_umep\processor\urock_processing_algorithm.py", line 447, in processAlgorithm
MainCalculation.main(javaEnvironmentPath = javaEnvVar,
File "C:\Users/Laura/AppData/Roaming/QGIS/QGIS3\profiles\default/python/plugins\processing_umep\functions\URock\MainCalculation.py", line 540, in main
InitWindField.calculates3dVegWindFactor(cursor = cursor,
File "C:\Users/Laura/AppData/Roaming/QGIS/QGIS3\profiles\default/python/plugins\processing_umep\functions\URock\InitWindField.py", line 1892, in calculates3dVegWindFactor
cursor.execute(";".join(["""
File "C:\Users\Laura\AppData\Roaming\Python\Python39\site-packages\jaydebeapi_init_.py", line 536, in execute
handle_sql_exception()
File "C:\Users\Laura\AppData\Roaming\Python\Python39\site-packages\jaydebeapi_init.py", line 165, in handle_sql_exception_jpype
reraise(exc_type, exc_info[1], exc_info[2])
File "C:\Users\Laura\AppData\Roaming\Python\Python39\site-packages\jaydebeapi_init.py", line 57, in reraise
raise value.with_traceback(tb)
File "C:\Users\Laura\AppData\Roaming\Python\Python39\site-packages\jaydebeapi_init_.py", line 534, in execute
is_rs = self._prep.execute()
jaydebeapi.DatabaseError: org.h2.jdbc.JdbcSQLDataException: Invalid value "-0.06281183352066441" for parameter "LN() argument"; SQL statement:
CREATE TABLE VEGETATION_OPEN_TEMPO_3DPOINTS
AS SELECT b.ID_POINT, a.ID_Z, CASE WHEN a.Z>b.MAX_HEIGHT
THEN LOG((a.Z - 3 * 0.05 * MAX_CANOPY_HEIGHT) / 2.3880850405465655) / LOG(a.Z / 2.3880850405465655)
ELSE CASE WHEN a.Z > b.MAX_HEIGHT OR a.Z < b.MIN_HEIGHT
THEN LOG((b.MAX_CANOPY_HEIGHT - 3 * 0.05 * MAX_CANOPY_HEIGHT) / 2.3880850405465655) / LOG(a.Z / 2.3880850405465655)
ELSE LOG((b.MAX_HEIGHT - 3 * 0.05 * MAX_CANOPY_HEIGHT) / 2.3880850405465655) / LOG(a.Z / 2.3880850405465655) * EXP(b.ATTENUATIO * (a.Z / b.MAX_CANOPY_HEIGHT - 1))
END
END
AS VEGETATION_FACTOR
FROM Z_VALUES_20250112013951 AS a, VEGETATION_OPEN_POINTS AS b
WHERE a.Z > 0 [90008-232]

Execution failed after 671.54 secondi (11 minuti 12 secondi)

Do you know what i can do to make it work?

[Lalalashp]

Hello everybody!
I tryed to overcome the errore analysing the formula that it gave me in the error message.
Basically, the LOG needed to be >0 so i changed the parameters in a way that this would occour.
specifically a.Z>0,15*MaxHeight which means if i had veg_height of 37m the minimum height could only have been 5.5 m. In this way i made it work for my area, even tho i still have to try this on another one, due to the long computational time of the tool.
Anyway, this was my way around the error, if there's a better path that doesn't include manipulating the dataset let me know !
(btw, i have the latest version of UMEP downloaded directly from here and 26gb RAM)

Laura

[j3r3m1]

Hi Laura and thank you for reporting. Can you share the original AND the reduced size datasets that you have used and that lead to the error ?

[Lalalashp]

Hello and thanks for the support.
As I explained in the last comment I made, I manipulated the dataset to try make it work, and it didn't in fact give me the error message anymore. But when I tried to re-run the processor, it started to slow down my pc a lot to the point I couldn't work on it, so I needed to shut it down. It occurred again on another smaller area, which is not the same as this one tho. To solve this problem, I opened another issue (here: #703) where I have explained better what happen. I'm reporting it because maybe they are related, and i shouldn't have manipulated the vegetation dataset, i dunno.

Anyway, here's the data for this particular error message.
These are the two datasets used for the URock pre-processor and processor. The original output of the vegetation of the pre-processor, obtained with the input data contained in the folder '19_Input.prepare', is contained in the folder '19_UrP_veget_x', which is separated due to the size of the file. This one has all the not-modified height.
In the folder '19_Input.Wind.Field' are contained the original building output of the pre-processor and the cliped vegetation one, to which I already modify the height field, as i explained in the comment above. The data in this folder are the ones that are causing the error message. Hope i was clear enought, it is a bit tricky with all these files.
19_Input_Urock prepare.zip
19_Input_Urban Wind Field.zip

19_UrP_veget_1.zip
19_UrP_veget_2.zip