geoclimate icon indicating copy to clipboard operation
geoclimate copied to clipboard

Problem generating LCZ classification for BDTopo version 3.3

Open EmmanuelleKerjouan opened this issue 1 year ago • 7 comments

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]"

EmmanuelleKerjouan avatar Dec 08 '23 10:12 EmmanuelleKerjouan

I'm on it

ebocher avatar Dec 08 '23 13:12 ebocher

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]]

ebocher avatar Dec 08 '23 14:12 ebocher

Yes, exactly.

EmmanuelleKerjouan avatar Dec 08 '23 14:12 EmmanuelleKerjouan

I'm on it

ebocher avatar Dec 08 '23 15:12 ebocher

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 geoms

ebocher avatar Dec 08 '23 16:12 ebocher

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.

EmmanuelleKerjouan avatar Dec 21 '23 16:12 EmmanuelleKerjouan

This issue is related to the H2GIS FILE_TABLE function

ebocher avatar Jan 09 '24 14:01 ebocher