import sys\nfrom osgeo import ogr\n\ndef get_default_scale(geopackage_path, default_enc):\n    # Ouvrir le GeoPackage en mode lecture seule\n    geopackage = ogr.Open(geopackage_path, 0)\n    default_enc = default_enc +'.000'\n\n    if geopackage is not None:\n        # R\u00e9cup\u00e9rer la couche DSID\n        dsid_layer = geopackage.GetLayerByName('DSID')\n\n        # V\u00e9rifier si la couche DSID existe\n        if dsid_layer:\n            # Rechercher la valeur de DSID_CSCL correspondant \u00e0 DSID_DSNM\n            dsid_layer.SetAttributeFilter(f\"DSID_DSNM = '{default_enc}'\")\n            feature = dsid_layer.GetNextFeature()\n            if feature:\n                default_scale = feature.GetField('DSPM_CSCL')\n                \n\n    return default_scale\ndef get_default_purpose(geopackage_path, default_enc):\n    # Ouvrir le GeoPackage en mode lecture seule\n    geopackage = ogr.Open(geopackage_path, 0)\n    default_purpose = '0'\n    default_enc = default_enc +'.000'\n\n    if geopackage is not None:\n        # R\u00e9cup\u00e9rer la couche DSID\n        dsid_layer = geopackage.GetLayerByName('DSID')\n\n        # V\u00e9rifier si la couche DSID existe\n        if dsid_layer:\n            # Rechercher la valeur de DSID_INTU correspondant \u00e0 DSID_DSNM\n            dsid_layer.SetAttributeFilter(f\"DSID_DSNM = '{default_enc}'\")\n            feature = dsid_layer.GetNextFeature()\n            if feature:\n                default_purpose = feature.GetField('DSID_INTU')\n                \n\n    return default_purpose\n\n\ndef delete_empty_tables_in_geopackages_dsid(geopackage_path, default_enc):\n    # R\u00e9cup\u00e9rer la valeur de scale \u00e0 partir de la table DSID\n    default_scale = get_default_scale(geopackage_path, default_enc)\n    default_purpose = get_default_purpose(geopackage_path, default_enc)\n    # Continuer avec le reste du code\n    if default_scale is None:\n        print(\"La valeur par d\u00e9faut de l'\u00e9chelle n'a pas \u00e9t\u00e9 trouv\u00e9e dans la table DSID.\")\n        return\n    default_scale_str = str(default_scale)\n    default_purpose_str = str(default_purpose)\n    # Ouvrir le GeoPackage en mode \u00e9dition\n    geopackage = ogr.Open(geopackage_path, 1)\n\n    if geopackage is not None:\n        # R\u00e9cup\u00e9rer le nombre de tables dans le GeoPackage\n        num_tables = geopackage.GetLayerCount()\n\n        # Liste des tables \u00e0 supprimer\n        tables_to_delete = []\n        tables_to_update = []\n\n        # Identifier les tables \u00e0 supprimer\n        for i in range(num_tables):\n            table = geopackage.GetLayerByIndex(i)\n\n            # V\u00e9rifier si la table est vide (aucun enregistrement)\n            if table.GetFeatureCount() == 0:\n                tables_to_delete.append(table.GetName())\n            else:\n                tables_to_update.append(table.GetName())\n\n        # Supprimer les tables\n        for table_name in tables_to_delete:\n            geopackage.DeleteLayer(table_name)\n            print(f\"Table supprim\u00e9e dans {geopackage_path}: {table_name}\")\n\n        # Mettre \u00e0 jour les tables restantes\n        for table_name in tables_to_update:\n            table = geopackage.GetLayerByName(table_name)\n\n            # V\u00e9rifier si le champ 'enc_chart' existe\n            enc_chart_index = table.FindFieldIndex('enc_chart', 1)\n            if enc_chart_index < 0:\n                # Ajouter le champ 'enc_chart' s'il n'existe pas\n                champ1 = ogr.FieldDefn('enc_chart', ogr.OFTString)\n                champ1.SetWidth(50)\n                champ1.SetNullable(True)\n                champ1.SetDefault(default_enc)\n                table.CreateField(champ1)\n\n            # V\u00e9rifier si le champ 'scale' existe\n            scale_index = table.FindFieldIndex('scale', 1)\n            if scale_index < 0:\n                # Ajouter le champ 'scale' s'il n'existe pas\n                champ2 = ogr.FieldDefn('scale', ogr.OFTReal)\n                champ2.SetWidth(10)\n                champ2.SetPrecision(2)\n                champ2.SetDefault(default_scale_str)  # Convertir en nombre flottant\n                table.CreateField(champ2)\n                \n            # V\u00e9rifier si le champ 'purpose' existe\n            purpose_index = table.FindFieldIndex('purpose', 1)\n            if purpose_index < 0:\n                # Ajouter le champ 'purpose' s'il n'existe pas\n                champ3 = ogr.FieldDefn('purpose', ogr.OFTReal)\n                champ3.SetWidth(10)\n                champ3.SetPrecision(2)\n                champ3.SetDefault(default_purpose_str)  # Convertir en nombre flottant\n                table.CreateField(champ3)\n                \n            # Mettre \u00e0 jour les valeurs dans les champs seulement si les champs n'ont pas d\u00e9j\u00e0 de valeurs\n            for feature in table:\n                enc_chart_value = feature.GetField('enc_chart')\n                scale_value = feature.GetField('scale')\n                purpose_value = feature.GetField('purpose')\n\n                # V\u00e9rifier si les champs n'ont pas d\u00e9j\u00e0 de valeurs d\u00e9finies\n                if enc_chart_value is None or scale_value is None or purpose_value is None:\n                    feature.SetField('enc_chart', default_enc)\n                    feature.SetField('scale', default_scale)\n                    feature.SetField('purpose', default_purpose)\n                    table.SetFeature(feature)\n                    \n\n        # Fermer le GeoPackage\n        geopackage = None\n\n    else:\n        print(f\"Impossible d'ouvrir le GeoPackage {geopackage_path} en mode \u00e9dition.\")\n\nif __name__ == \"__main__\":\n    if len(sys.argv) != 3:\n        print(\"Usage: python script.py geopackage_path default_enc\")\n        sys.exit(1)\n    \n    geopackage_path = sys.argv[1]\n    default_enc = sys.argv[2]\n\n    delete_empty_tables_in_geopackages_dsid(geopackage_path, default_enc)<\/code><\/pre>\n\n\n\n<\/div><\/div>\n\n\n\n
<\/p>\n\n\n\nYou can download the script here<\/a>
Copy the script to the same location as the scripts in the original article.<\/p>\n\n\n\n
You can either replace the lines in the .bat files<\/p>\n\n\n\n
python c:\/testgpkg\/update_geopackage_dsid.py \u00ab\u00a0%output_geopackage%\u00a0\u00bb \u00ab\u00a0!file!\u00a0\u00bb<\/em><\/strong><\/p>\n\n\n\n
by<\/p>\n\n\n\n
python c:\/testgpkg\/update_geopackage_dsid_prp.py \u00ab\u00a0%output_geopackage%\u00a0\u00bb \u00ab\u00a0!file!\u00a0\u00bb<\/em><\/strong><\/p>\n\n\n\n
Either delete the update_geopackage_dsid.py<\/em><\/strong> script and rename update_geopackage_dsid_prp.py<\/em><\/strong> to update_geopackage_dsid.py<\/em><\/strong>, which saves you having to modify the .bat files.<\/p>\n\n\n\n
Add a Python script to filter layers by purpose<\/h2>\n\n\n\nTo filter the display of all layers loaded in the QGis project, you can use the following script:<\/p>\n\n\n\n<\/div>filter_purpose.py<\/div><\/div><\/div> <\/p>\n\n\n\n# Importer le module QGIS\nfrom qgis.core import QgsProject, QgsMapLayerType\n\n# D\u00e9finir le valeur de l'attribut \"purpose\" \u00e0 filtrer\nvaleur_purpose = 3\n\n# Obtenir le projet actif\nprojet = QgsProject.instance()\n\n# Obtenir la liste des couches charg\u00e9es dans le projet\ncouches = projet.mapLayers().values()\n\n# Parcourir toutes les couches\nfor couche in couches:\n    # V\u00e9rifier si la couche est de type vecteur\n    if couche.type() == QgsMapLayerType.VectorLayer:\n        # V\u00e9rifier si la couche a un champ nomm\u00e9 \"purpose\"\n        if couche.fields().indexFromName('purpose') != -1:\n            # V\u00e9rifier si le nom de la couche commence par l'un des pr\u00e9fixes sp\u00e9cifi\u00e9s\n            if couche.name().startswith(('pt_', 'li_', 'pl_')):\n                # D\u00e9finir le filtre sur l'attribut \"purpose\" \u00e9gal \u00e0 la valeur sp\u00e9cifi\u00e9e\n                filtre = f'\"purpose\" = {valeur_purpose}'\n                couche.setSubsetString(filtre)\n                print(f\"Filtre d\u00e9fini pour la couche {couche.name()}\")\n            else:\n                print(f\"La couche {couche.name()} ne commence pas par 'pt_', 'li_', ou 'pl_'\")\n        else:\n            print(f\"La couche {couche.name()} n'a pas d'attribut 'purpose'\")\n    else:\n        print(f\"La couche {couche.name()} n'est pas une couche vecteur\")\n<\/code><\/pre>\n\n\n\n<\/div><\/div>\n\n\n\n
<\/p>\n\n\n\nYou can download the script here.<\/a><\/p>\n\n\n\n
To eliminate Purpose filters, use the same script as for scaling (unfilter.py).<\/p>\n\n\n\n
The S57 database project with Geopackage has come to an end. We’re now working on setting up an equivalent procedure using a PostgreSQL\/Postgis database. Help us bring this project to a successful conclusion<\/h3>\n\n\n\n\n\t\t\t\t\t\t\n\t\t\t\t\t\t\tCe formulaire n\u2019est pas publi\u00e9.\n\t\t\t\t\t\t<\/p>\n\t\t\t\t\t\t\n\t\t\t\t\t<\/div>\n","protected":false},"excerpt":{"rendered":"
Following the publication of the first part of our work on Geopackages, our work with PostgreSQL\/Postgis has enabled us to make a series of further advances. In this article, you’ll find an updated procedure for importing…<\/p>\n","protected":false},"author":1,"featured_media":12654,"comment_status":"open","ping_status":"closed","sticky":false,"template":"","format":"standard","meta":{"give_campaign_id":0,"_bbp_topic_count":0,"_bbp_reply_count":0,"_bbp_total_topic_count":0,"_bbp_total_reply_count":0,"_bbp_voice_count":0,"_bbp_anonymous_reply_count":0,"_bbp_topic_count_hidden":0,"_bbp_reply_count_hidden":0,"_bbp_forum_subforum_count":0,"sfsi_plus_gutenberg_text_before_share":"","sfsi_plus_gutenberg_show_text_before_share":"","sfsi_plus_gutenberg_icon_type":"","sfsi_plus_gutenberg_icon_alignemt":"","sfsi_plus_gutenburg_max_per_row":"","_monsterinsights_skip_tracking":false,"_monsterinsights_sitenote_active":false,"_monsterinsights_sitenote_note":"","_monsterinsights_sitenote_category":0,"jetpack_post_was_ever_published":false,"_jetpack_newsletter_access":"","_jetpack_dont_email_post_to_subs":false,"_jetpack_newsletter_tier_id":0,"_jetpack_memberships_contains_paywalled_content":false,"_jetpack_memberships_contains_paid_content":false,"footnotes":""},"categories":[1377],"tags":[1995,2817,1386,2721],"class_list":["post-12648","post","type-post","status-publish","format-standard","has-post-thumbnail","hentry","category-qgis-en","tag-geopackage-en","tag-import-en","tag-qgis-en","tag-s57-en"],"aioseo_notices":[],"jetpack_featured_media_url":"https:\/\/i0.wp.com\/www.sigterritoires.fr\/wp-content\/uploads\/2024\/03\/purpose5-1.jpg?fit=1200%2C860&ssl=1","jetpack_shortlink":"https:\/\/wp.me\/p6XU0A-3i0","jetpack_sharing_enabled":true,"_links":{"self":[{"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/posts\/12648","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/comments?post=12648"}],"version-history":[{"count":0,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/posts\/12648\/revisions"}],"wp:featuredmedia":[{"embeddable":true,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/media\/12654"}],"wp:attachment":[{"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/media?parent=12648"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/categories?post=12648"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/tags?post=12648"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}