Message posté par : Christophe BADOL
----------------------------------------
-----------------
Hydrolithe écrit :
Bonjour,
En lançant gdalinfo avec l'option -approx_stats sur le .vrt de la BD Otho à l'échelle d'un département, cela prend un temps très long. J'ai laissé tourner la commande plus de 10 heures et le calcul de la bande 1 n'était pas terminé.
Est-ce qu'il y a une subtilité que je n'ai pas vue ou bien est-ce long dans tous les cas ?
Pierre
-----------------
Un gdalinfo avec l'option -approx_stats sur un VRT de la BD Ortho HR d'un département prends toujours moins d'une heure, même pour les plus gros départements.
Le problème doit être ailleurs !
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=353640#p353640
Pour y répondre : geolibre(a)ml.georezo.net ou reply de votre messagerie
Pour vous désabonner connectez-vous sur le forum puis Profil / Abonnement
--
Association GeoRezo - le portail géomatique
https://georezo.net
Message posté par : Thais781
----------------------------------------
Bonjour,
Cf. l'image ci-joint.
J'ai un polygon (en bleu) et un point (celui en orange)
Et je cherche a calculer les points de contact les plus proches entre ce point et le polygon
Pour le moment, je trace un cercle a partir du point orange et je fais évoluer le diamètre jusqu'a qu'il touche (intercepts) le polygone et j'ai donc le point de contact.
Cette solution est un peu ... bourin ... peut être qu'il existe une solution plus élégante pour obtenir ces point (les verts dans le dessin)
N'hésitez pas si vous avez une idée ;-)
Thais
----------------------------------------
Ce message est accompagné de fichiers, pour les télécharger, suivre le lien ci-dessous.
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=356620#p356620
Pour y répondre : geolibre(a)ml.georezo.net ou reply de votre messagerie
Pour vous désabonner connectez-vous sur le forum puis Profil / Abonnement
--
Association GeoRezo - le portail géomatique
https://georezo.net
Message posté par : PierreVit
----------------------------------------
Bonjour,
Je suis à mon compte et une entité publique me demande si j'ai une solution pour leur proposer un visuel SIG avec mise à jour simple sans passer par l'installation d'un logiciel interne et open source de cartographie type QGis.
Auriez vous une idée de ce que je pourrais leur proposer ? c'est la première fois que j'ai ce cas. Car à part leur proposer QGis pour que le fichier reste en interne, je sèche un peu.
Pour les détails techniques :
Il n'y a pas forcement besoin de géolocalisation
Ce sera un fichier de points tous alignés avec moins de 10 couches, et en "fond de plan" un rectangle qui entoure les points.
La mise à jour serait idéalement faite dans un tableur, rattaché aux points par leurs identifiants bien entendu.
En vous remerciant ! N'hésitez pas à me demander d'autres infos.
Pierre
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=355343#p355343
Pour y répondre : geolibre(a)ml.georezo.net ou reply de votre messagerie
Pour vous désabonner connectez-vous sur le forum puis Profil / Abonnement
--
Association GeoRezo - le portail géomatique
https://georezo.net
Message posté par : FredDuhart
----------------------------------------
Bonjour à tous,
Je cherhce à enregistrer le flux WMS du RGE alti de l'IGN, avec gdal_translate. Voici le code :
-----------------
Code :
gdal_translate -tr 1 1 -b 1 "WMS:https://wxs.ign.fr/altimetrie/geoportail/r/wms?SERVICE=WMS&VERSION=1.3.0&REQUEST=GetMap&LAYERS=RGEALTI-MNT_PYR-ZIP_FXX_LAMB93_WMS&CRS=IGNF:LAMB93&BBOX=476250,6381200,477100,6381600" rge_sortie.tiff
-----------------
La capture se passe bien mais toutes les valeurs du .tiff sont arrondies à l'unité (au mètre). Est-ce normal, et est-ce possible d'avoir une précision centimétrique comme sur le geoportail?
merci d'avance pour vos réponses.
Frédéric
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=354151#p354151
Pour y répondre : geolibre(a)ml.georezo.net ou reply de votre messagerie
Pour vous désabonner connectez-vous sur le forum puis Profil / Abonnement
--
Association GeoRezo - le portail géomatique
https://georezo.net
Message posté par : Pierre
----------------------------------------
Bonjour,
En lançant gdalinfo avec l'option -approx_stats sur le .vrt de la BD Otho à l'échelle d'un département, cela prend un temps très long. J'ai laissé tourner la commande plus de 10 heures et le calcul de la bande 1 n'était pas terminé.
Est-ce qu'il y a une subtilité que je n'ai pas vue ou bien est-ce long dans tous les cas ?
Pierre
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=351779#p351779
Pour y répondre : geolibre(a)ml.georezo.net ou reply de votre messagerie
Pour vous désabonner connectez-vous sur le forum puis Profil / Abonnement
--
Association GeoRezo - le portail géomatique
https://georezo.net
Message posté par : image95
----------------------------------------
Bonjour,
Je travaille sous windows. Postgis et gdal (j'utilise gdal/ogr2ogr depuis l'invit de commande. Version de gdal disponible sous C:\Program Files\QGIS 3.16\bin).
Jusqu'à maintenant, au sein de ma base postgis, j'effectuais des insertions de nouvelles entités au sein de table existante
via des requetes SQL du type :
-----------------
Code :
INSERT INTO SCHEMA_humain.TableDestination (
"fieldXX",
"geom")
SELECT
"fieldXX",
"geom"
FROM table_Source
-----------------
Ou bien des mises à jour d'entités (update)
grace à
-----------------
Code :
update
set
from
where
-----------------
Celà fonctionnait. Mais la grosse limite de cette approche est qu'elle necessite au prealable d'importer au sein de la bdd postgis les shp.
Mon but est donc de directement inserer de nouvelles entités / mettre à jour (depuis fichier source shp/table destinatation postgis en utilisant gdal.
Afin de s'affranchir de cette importation de shp ds la bdd.
Le fichier source est 1 fichier shapefile de POINTS (EPSG 2154)
La table cible/destination est donc une table postgis deja existante de POINTS (EPSG 2154)
(Je précise que le schema de la table de destination n'est PAS le schema public mais le schema 'humain'.)
Il semble y avoir 2 approches pour accomplir cette tache. Soit d'utiliser les options -append associé à -fieldmap.
Soit l'approche avec -sql en définissant la requete SQL.
Je privilegie si possible l'approche avec -sql.
Voici la structure à laquelle je pensais.
-----------------
Code :
ogr2ogr - progress - append D:\XXX.shp
-f PG:"dbname='dbname' host='ipXXXXX' port='5432' user='XXXX' password='XXXX' active_schema='humain'"
-sql
INSERT INTO schema_humain.TableName (
"champ1",
"champ2",
"champ3",
"geom")
--selection des champs du SHP
SELECT
"champ1",
"champ2",
"champ3",
"geom"
FROM
shp
-----------------
Une personne pourrait elle me dire si je me fourvoie ? Avec cette méthode, est il possible de specifier les champs du shp au sein de la clause SQL ? (champ non geom et champ geom?). Lorsque comme moi on a en entrée du shp et en sortie une table postgis, doit on plutot prioriser l'approche avec -append -addfield ou -update plutot que l'approche avec -sql ?
Merci beaucoup.
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=353329#p353329
Pour y répondre : geolibre(a)ml.georezo.net ou reply de votre messagerie
Pour vous désabonner connectez-vous sur le forum puis Profil / Abonnement
--
Association GeoRezo - le portail géomatique
https://georezo.net
Message posté par : Elise Henrot
----------------------------------------
Bonjour,
Re-quelques années après, le lien ne fonctionne plus... ou je suis une cruche ! :-(
Elise
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=346527#p346527
Pour y répondre : geolibre(a)ml.georezo.net ou reply de votre messagerie
Pour vous désabonner connectez-vous sur le forum puis Profil / Abonnement
--
Association GeoRezo - le portail géomatique
https://georezo.net
Message posté par : Thomas Corvé (thomas.corve(a)gmail.com)
----------------------------------------
Bonjour à tous,
Je suis développeur d'une appli qui entres autres reçoit des coordonnées WGS84 d'un GNSS et qui doit placé la position de l'utilisateur sur un plan.
Ce plan est exprimé (par exemple) en CC49, et utilise une grille de correction Raf18.
Après avoir télécharger les grilles avec l'outil projsync, la correction altimétrique fonctionne (je n'avais pas les grilles avant).
Voici comment j'utilise proj:
PJ* pj = proj_create_crs_to_crs(m_projCtx, "epsg:4326", "epsg:3949+5720", nullptr);
J'essaye maintenant de faire la même chose pour un projet en allemagne. Pour la planimétrie j'utilisais le code epsg:4647.
Pour ajouter la correction verticale, il semble que mes clients utilisent "DHHN2016". Après une recherche sur espg.io je trouve le code epsg correspondant epsg:7837.
Malheureusement cela ne fonctionne pas. Après une recherche avec projinfo, le code 5720 donne une string Proj avec le nom d'un fichier grille, mais pas le code 7837:
-----------------
Code :
projinfo epsg:5720
PROJ.4 string:
+geoidgrids=fr_ign_RAF18.tif +vunits=m +no_defs +type=crs
-----------------
-----------------
Code :
projinfo epsg:7837
PROJ.4 string:
+vunits=m +no_defs +type=crs
-----------------
Du coup je me suis dit qu'il fallait modifier la string Proj pour ajouter le nom du fichier Grille. Dans la doc, il est mention de fichiers ini, dans lesquel on peut modifier la string avec le mot clef +init. J'ai donc créé un fichier avec la ligne suivante:
-----------------
Code :
+init=epsg:7837 +geoidgrids=monfichier.tiff
-----------------
Je n'ai pas constaté d'effet. Et projinfos ne m'a pas ajouter cette ligne dans la string Proj.
Si une bonne âme pouvait me donner quelques pistes... Éventuellement on cherche un contractuel pour nous aider sur cette partie...
Bien cordialement, Thomas C.
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=352956#p352956
Pour y répondre : geolibre(a)ml.georezo.net ou reply de votre messagerie
Pour vous désabonner connectez-vous sur le forum puis Profil / Abonnement
--
Association GeoRezo - le portail géomatique
https://georezo.net
Message posté par : Bruno IRATCHET (bruno(a)georezo.net)
----------------------------------------
Merci de poster dans ce forum.
----------------------------------------
Le message est situé https://georezo.net/forum/viewtopic.php?pid=352907#p352907
Pour y répondre : geolibre(a)ml.georezo.net ou reply de votre messagerie
Pour vous désabonner connectez-vous sur le forum puis Profil / Abonnement
--
Association GeoRezo - le portail géomatique
https://georezo.net