geoclimate
geoclimate copied to clipboard
Problem generating LCZ classification for BDTopo version 3.3
The aim is to generate an LCZ classification for the Nancy Metropole area, with data from the BDTopo 3.3. The coordinates of the zone are as follows : 6836343.6284,924280.9312,6855225.3440,944570.1370 (EPSG : 2154) However, we obtain this error: "Caught: org.h2.jdbc.JdbcSQLSyntaxErrorException: Table "ZONE_D_ACTIVITE_OU_D_INTERET_TMP" non trouvée Table "ZONE_D_ACTIVITE_OU_D_INTERET_TMP" not found; SQL statement: CREATE TABLE zone_d_activite_ou_d_interet AS SELECT ID, st_setsrid(the_geom, 2154) as the_geom, CATEGORIE, NATURE, FICTIF FROM zone_d_activite_ou_d_interet_tmp WHERE the_geom && 'SRID=0;POLYGON ((923280.9312 6835343.6284, 923280.9312 6856225.344, 945570.137 6856225.344, 945570.137 6835343.6284, 923280.9312 6835343.6284))'::GEOMETRY AND ST_INTERSECTS(the_geom, 'SRID=0;POLYGON ((923280.9312 6835343.6284, 923280.9312 6856225.344, 945570.137 6856225.344, 945570.137 6835343.6284, 923280.9312 6835343.6284))'::GEOMETRY) AND (CATEGORIE='Administratif' OR CATEGORIE='Enseignement' OR CATEGORIE='Santé') [42102-224]"
I'm on it
To be sure about the issue. You set a bbox as input to process the Nancy area right ? So you have locations : [[6836343.6284,924280.9312,6855225.3440,944570.1370]]
Yes, exactly.
I'm on it
So I run a test case with the bbox [[6836343.6284,924280.9312,6855225.3440,944570.1370]] and the BDTopo 3.3 (version BDTOPO_3-3_TOUSTHEMES_SHP_LAMB93_D054_2023-09-15) downloaded on IGNF website and I'm not able to reproduce the issue.
See the result LCZ
The problem isn't with version 3.3 of BDTopo. The error probably stems from the number of characters in the path where it fetches the data. The test was performed with two different paths:
- First path : /home/emmanuelle/Documents/Mission_CDD/Comparaisons_LCZ/AZERtx50 (306 characters)
- Second path : /home/emmanuelle/Documents/Mission_CDD/Données/Nancy/BDTopo3/BDTOPO_3-3_TOUSTHEMES_SHP_LAMB93_D054_2023-09-15/BDTOPO/1_DONNEES_LIVRAISON_2023-09-00196/BDT_3-3_SHP_LAMB93_D054-ED2023-09-15 (187 characters).
In the first case, the error occurs, but not in the second.
This issue is related to the H2GIS FILE_TABLE function