Skip to Content

GeoRezo

Syndiquer le contenu GeoRezo.net
GeoRezo.net, le Portail de la Géomatique et des SIG (Systèmes d'Information Géographique)
Mis à jour : il y a 3 heures 1 min

[Géomatique] Re: variogramme croisé

ven 03-07-2015
Nombre de message(s) précédent(s) : 1

vous devriez poster votre question sur 52 nord plutôt que sur Georezo.

Catégories: Sites Francophones

[QGIS] Reconstitution du relief avec QGIS

jeu 02-07-2015

Bonsoir,

Je suis novice dans la cartographie donc je m'excuse par avance si mes questions vous paraissent naïves. Je viens de télécharger la banque de données Litto3D du shom afin de visualiser le relief de nos archipels, mon problème c'est que les données sont segmentées dans des carrés de 1km de côté et pour certains archipels, les map sont réparties dans des milliers de MNT au format ASC. Je me retrouve donc à devoir reconstituer un relief à partir d'un peu plus de 4000 fichiers asc d'environ 300 ko chacun. Est ce que QGIS est capable de faire ce travail et si non, pouvez vous m'aiguiller vers un logiciel adapté? Bien sûr les données ne sont pas complètes il y a pas mal de trous, j'imagine que QGIS est capable de moyenner l'information pour combler les trous? Dans un premier temps j'ai essayé d'afficher l'information avec Matlab mais ce n'est pas un logiciel de cartographie et je me retrouve avec des problèmes de performance pour visualiser mes reliefs. Pour QGIS je n'ai pas bien compris comment fusionner mes fichiers, si vous avez des informations je suis preneur.

Par avance, merci de l'aide que vous pourrez m'apporter.

Julien

Catégories: Sites Francophones

[Webmapping] Re: Géorigami

jeu 02-07-2015
Nombre de message(s) précédent(s) : 9

Extraordinaire. Félicitation pour l'interface et le travail.

Catégories: Sites Francophones

[Géomatique] Re: Portage d'Apic vers Elyx - documents méthodologiques et couts

jeu 02-07-2015
Nombre de message(s) précédent(s) : 1

Aloha

Les utilisateurs de chez un Spatial ont en effet tendance à la discrétion.

Cordialement.

Catégories: Sites Francophones

[Mapinfo] Re: Problème code MAPBASIC

jeu 02-07-2015
Nombre de message(s) précédent(s) : 4

Bonjour
De mon côté j'aimerai bien apprendre MapBasic mais j'ai du mal à trouver une doc "débutant".
Comment as-tu fait pour apprendre ?
Merci
Aurélie

Bonjour Aurelie,

Alors pour ma part j'ai commencé à apprendre MapBasic depuis deux semaines environ, et je n'y avait jamais touché avant (c'est pour cela que je demande parfois de l'aide =) ).

Pour débuter, j'ai utilisé la fenetre MapBasic dans Mapinfo pour afficher le code qui émanait des différentes manipulations effectuées sous Mapinfo pour comprendre un petit peu comment cela marchait. Ensuite, j'ai approfondi avec de la documentation (en anglais) pour apprendre d'autres éléments et je continue a y jeter un oeil régulièrement.

Cela viendra avec l'entrainement je suppose, il faut juste apprendre tout d'abord les bases du langage et savoir ce que l'on peut/ce que l'on ne peut pas faire avec cet outils ;) . Avec un peu de bonne volonté on devrait pouvoir s'en sortir =)

Catégories: Sites Francophones

[Esri] Re: ArcGIS - Precondition Model Builder

jeu 02-07-2015
Nombre de message(s) précédent(s) : 1

Bonjour,

Si j'ai bien compris

non seulement le nom du fichier change à chaque itération
mais le nom du dossier sous "Climate Data Folder" change également

Vous voudriez que
EC25publicaProyecto_amanezadasFloraAscii_Msnom_1 change
en
EC25publicaProyecto_amanezadasFloraAscii_Msnom_1
nom_2

Pour cela je vous propose d'essayer l'outil Calculate Value (Model only tools)
Avec
            Data type : Folder
            Expression EC25publicaProyecto_amanezadasFloraAscii_Ms%Name%

Vous utiliserez alors "output_value" pour remplacer "Climate Data Folder"

Cela devrait pouvoir résoudre votre problème

--
Damien ROY
IgE MSHE ledoux

Catégories: Sites Francophones

[Géo'BD] Trigger postgresql : valeur NEW qui ne remonte pas

jeu 02-07-2015

Bonjour.
J'ai un souci avec un trigger.
Celui-ci se déclenche sur un update de la géométrie ou d'un champs nommé 'libelle'.
Lorsque je lance un update avec la valeur de libelle = 'Conforme' , le trigger se déclenche bien, fait ce qu'il faut, mais par contre le NEW.libelle n'est pas remontée dans ma base, la valeur du champ reste en l'état (avec l'ancienne valeur).

CREATE TRIGGER t_test_trigger
  BEFORE UPDATE OF the_geom, libelle
  ON sig_conformite
  FOR EACH ROW
  EXECUTE PROCEDURE test_trigger();

....
BEGIN

nb_dossier = 0;
nb_dossier = count(nom_dossier)from sig_conformite where nom_dossier = NEW.nom_dossier;

IF NEW.libelle = 'Conforme' THEN
    NEW.alerte_date = 1;
    raise notice 'nb dossiers :%', nb_dossier;
    IF nb_dossier > 1 THEN
        update sig_conformite set alerte_date = 1 where nom_dossier = NEW.nom_dossier;
    END IF;
END IF;

return NEW;



En simplifiant le trigger je me suis rendu compte que c'est la ligne update sig_conformite set alerte_date = 1 where nom_dossier = NEW.nom_dossier  qui met le bazar.
Je ne comprend pas pourquoi. J'ai supprimé tous les autres trigger (en cas d'interaction) mais le phénomène perdure. Quand je fais un raise notice sur NEW.libelle la valeur est bien 'Conforme', mais en final dans ma table, elle reste à 'Non conforme'.
Si quelqu'un à une idée sur la source du problème je suis preneur.
Cordialement
JP

Catégories: Sites Francophones

[Esri] Re: Comment s'y prendre? Automatisation d'un processus

jeu 02-07-2015
Nombre de message(s) précédent(s) : 4

J'ai le vague souvenir qu'avec l'outil de mise à jour quelque chose comme cela est possible le long d'une ligne.

Méthodologie:

Dans un premier temps, place ton .shp de linéaire dans un dossier
Puis creer un .shp de point vierge dans ce même dossier en passant par ArcCatalog
Finalement passe en mode mise à jour
Sélectionne ton linéaire
puis vas dans Editeur / Construct Points tu devrais pouvoir rentrer un critère de distance le long du linéaire
Un point est alors généré tout les 500m (10km ou autres) le long de ton linéaire

Si tu as une ligne qui représente ton cours d'eau alors c'est gagné.
Sinon tu peux toujours en généré une à partir de l enveloppe externe de tes polygones.

Damien Roy
IgE MSHE Ledoux

Catégories: Sites Francophones

[Géomatique] Re: elongation, formes des polygones

jeu 02-07-2015
Nombre de message(s) précédent(s) : 15

Bonjour,
je suis désolé pour le déterrage de topic mais je me réfère à un lien récent (https://github.com/igeofr/qgis2).

Je cherche à calculer les indices de Morton pour une série de polygones mais les résultats obtenus ne sont pas satisfaisants.

Indice de Morton : Rapport de la superficie d'une entité, à celle d'un cercle de même périmètre = Imor = Area / Pi*(0,5*L)²
(Vaut 1 si le polygone est un cercle, 0 s'il est de surface nulle)


Le modèle proposé pour QGis (https://github.com/igeofr/qgis2/blob/ma … rton.model) propose cette formule :

Code:($area)/($pi*((0.5*bounds_height( $geometry ))*(0.5*bounds_height( $geometry ))))

J'ai l'impression que ce mode de calcul rend l'indice dépendant de l'orientation du polygone car pour des géométries similaires on a des indices significativement différents. Et accessoirement les indices vont bien au-delà de 1 alors que je pensais qu'il était borné à 1...

Et dernier point, les résultats sont encore différents si je les calcule manuellement.

D'où ma question :
Est-ce pertinent d'assimiler le bounds_height à L largeur maximale du polygone ?

Catégories: Sites Francophones

[QGIS] Re: QGIS : déplacement de point

jeu 02-07-2015
Nombre de message(s) précédent(s) : 4

Bonjour,

Je pense que c'est un prob de paramétrage. Je vous envoie une copie d'écran

Catégories: Sites Francophones

[Données] Re: Fichiers MAJIC 2015

jeu 02-07-2015
Nombre de message(s) précédent(s) : 8

Re,
De pm66 :

merci de m'indiquer si j'ai oublié quelque chose dans mon raisonnement

Oui, c'est tout bon !
Mais, je maintiens : la livraison PCI de Mai-Juin, est celle qui offre le meilleur taux d'appariement avec les données MAJIC III reconstituées au 01/01.
A yartostout :
Les données au 01/01/2015 sont prêtes à la diffusion, mais c'est du lourd, donc il se peut que les fichiers MAJIC III ne soient pas prêts sur tout le territoire: si vous êtes accrédités pour accéder au portail de la fonction publique, vous devriez les trouver.
Sinon, en faire la demande à votre correspondant local de la DGFiP.
Cordialement

Catégories: Sites Francophones

[Mapinfo] Re: MAPINFO 10 Mise à jour une colonne en fonction d'une localisation geo

jeu 02-07-2015
Nombre de message(s) précédent(s) : 3

Bonjour,

tiyan34 a écrit:

il me semble que ça n'existait pas dans les tres anciennes versions de Mapinfo

en fait si, cela existe depuis longtemps (au moins depuis la v4.5).

Nicolas.

Catégories: Sites Francophones

[Geolibre_web] Extraire les métadonnées d'OSM

jeu 02-07-2015

Bonjour à tous,

Je travaille sur des données OpenStreetMap récupérées en .shp et en .pbf
Le .shp étant bien pauvre en information, je me suis intéressé au .pbf

J'ai lu, ci-dessous, que le .pbf contenait des métadonnées
(http://wiki.openstreetmap.org/wiki/PBF_Format)

Savez vous si ces métadonnées sont bien renseignées? (au moins la date de mise à jours des polygones)
Sauriez vous extraire ces métadonnées?

J'ai essayé OsmPbfMetadata mais ce n'est pas ce que je recherche
J'ai lu que Osmosis pouvait être utilisé mais je ne sais pas comment...

Auriez-vous le liens pour des tutoriels d’utilisation d'Osmosis?

Merci d'avance!

Catégories: Sites Francophones

[QGIS] Re: QGIS 2.2 : plugin Cadastre

jeu 02-07-2015
Nombre de message(s) précédent(s) : 5

Bonjour j'utilise Qgis 2.8.2, avec le plugin Cadastre en 1.2.
Lorsque je lance l'import des données MAJIC III et EDIGEO, après quelque minutes de chargement j'ai ce message
Une erreur est apparue lors de l'exécution du code Python :

Traceback (most recent call last):
  File "C:/Documents and Settings/philippe/.qgis2/python/pluginscadastrecadastre_dialogs.py", line 1055, in processImport
    qi.importEdigeo()
  File "C:/Documents and Settings/philippe/.qgis2/python/pluginscadastrecadastre_import.py", line 631, in importEdigeo
    self.importAllEdigeoToDatabase()
  File "C:/Documents and Settings/philippe/.qgis2/python/pluginscadastrecadastre_import.py", line 1102, in importAllEdigeoToDatabase
    self.importEdigeoVecToDatabase(vec)
  File "C:/Documents and Settings/philippe/.qgis2/python/pluginscadastrecadastre_import.py", line 1238, in importEdigeoVecToDatabase
    self.qc.updateLog('Erreurs pendant la requête : %s' % sql)
UnboundLocalError: local variable 'sql' referenced before assignment

Version de Python :
2.7.4 (default, Apr  6 2013, 19:54:46) [MSC v.1500 32 bit (Intel)]

Version de QGIS :
2.8.2-Wien Wien, 1b929ef

Chemin vers Python : ['C:/Documents and Settings/philippe/.qgis2/python/plugins\processing', 'C:/PROGRA~1/QGISWI~1/apps/qgis/./python', u'C:/Documents and Settings/philippe/.qgis2/python', u'C:/Documents and Settings/philippe/.qgis2/python/plugins', 'C:/PROGRA~1/QGISWI~1/apps/qgis/./python/plugins', 'C:\PROGRA~1\QGISWI~1\bin\python27.zip', 'C:\PROGRA~1\QGISWI~1\apps\Python27\DLLs', 'C:\PROGRA~1\QGISWI~1\apps\Python27\lib', 'C:\PROGRA~1\QGISWI~1\apps\Python27\lib\plat-win', 'C:\PROGRA~1\QGISWI~1\apps\Python27\lib\lib-tk', 'C:\PROGRA~1\QGISWI~1\bin', 'C:\PROGRA~1\QGISWI~1\apps\Python27', 'C:\PROGRA~1\QGISWI~1\apps\Python27\lib\site-packages', 'C:\PROGRA~1\QGISWI~1\apps\Python27\lib\site-packages\PIL', 'C:\PROGRA~1\QGISWI~1\apps\Python27\lib\site-packages\jinja2-2.7.2-py2.7.egg', 'C:\PROGRA~1\QGISWI~1\apps\Python27\lib\site-packages\markupsafe-0.23-py2.7-win32.egg', 'C:\PROGRA~1\QGISWI~1\apps\Python27\lib\site-packages\python_dateutil-2.2-py2.7.egg', 'C:\PROGRA~1\QGISWI~1\apps\Python27\lib\site-packages\pytz-2014.2-py2.7.egg', 'C:\PROGRA~1\QGISWI~1\apps\Python27\lib\site-packages\win32', 'C:\PROGRA~1\QGISWI~1\apps\Python27\lib\site-packages\win32\lib', 'C:\PROGRA~1\QGISWI~1\apps\Python27\lib\site-packages\Pythonwin', 'C:\PROGRA~1\QGISWI~1\apps\Python27\lib\site-packages\Shapely-1.2.18-py2.7-win32.egg', 'C:\PROGRA~1\QGISWI~1\apps\Python27\lib\site-packages\six-1.6.1-py2.7.egg', 'C:\PROGRA~1\QGISWI~1\apps\Python27\lib\site-packages\wx-2.8-msw-unicode', u'C:/Documents and Settings/philippe/.qgis2//python', 'C:\Documents and Settings\philippe\.qgis2\python\plugins\cadastre/forms', 'C:\PROGRA~1\QGISWI~1\apps\qgis\python\plugins\fTools\tools', 'C:\Documents and Settings\philippe\.qgis2\python\plugins\cadastre/forms']

je me suis balader sur pas mal de site et forum mais je ne trouve pas la solution... J'ai bien toute les configuration au top...

Auriez-vous une solution?

Catégories: Sites Francophones

[Géomatique] variogramme croisé

jeu 02-07-2015

Bonjour,
Je travailactuellement sur une méthode de cokrigeage et je me suis penché sur le variogramme croisé qui selon moi doit représenter l'influence d'une donnée Z2(x) à la position x sur une donnée Z1(x+h) à la position x+h

Or la définition que l'on retrouve le plus souvent pour ce variogramme croisé est gamma(h) = E{Z1(x+h)-Z1(x)}{Z2(x+h)-Z2(x)} qui nécessite de connaitre chacune des deux données au même point ?

pourquoi n'a on pas plus simplement gamma(h) = E{Z1(x+h)-Z2(x)}^2 qui correspondrai plus a ma définition ?

sur quoi je me base : http://streaming.stat.iastate.edu/~stat … enotes.pdf

j'éspère que quelqu'un saura répondre a ma cette interrogation !
merci d'avance.

Catégories: Sites Francophones

[Géo'BD] Re: trigger qui boucle

jeu 02-07-2015
Nombre de message(s) précédent(s) : 6

Bonjour.
J'ai presque réussi à me dépatouiller avec mon problème.
J'ai testé les RULES mais je me suis rendu compte qu'elles avaient elles aussi tendance à se relancer en boucle avec le risque de ne pas en sortir (si on a oublié un test en cours de route ....)
J'ai donc joué sur la définition des colonnes qui déclenchent les triggers (UPDATE OF XXX) , ajouté des tests, distingué les triggers en insert et ceux en update, et enfin bloqué le déclenchement de certains triggers (disable /enable) lors d'un INSERT.
Merci donc pour vos suggestions qui m'ont bien inspiré.
Cordialement
JP

Catégories: Sites Francophones

[Mapinfo] Re: Etiquettes

jeu 02-07-2015
Nombre de message(s) précédent(s) : 2

Merci pour votre info. Je l'ai déjà fait mais je vais réessayé.
Merci, à bientôt.
Vous n'avez pas une autre solution ?
Cordialement.
Jérôme BAYLE.

Catégories: Sites Francophones

[Mapinfo] Re: MAPINFO 10 Mise à jour une colonne en fonction d'une localisation geo

jeu 02-07-2015
Nombre de message(s) précédent(s) : 1

Bonjour Yannick,

Pour cela ouvrez la table contenant vos points, et celle de vos communes.

Créez tout d'abord la colonne qui va contenir les codes insee.

Ensuite, allez dans Table > Mettre à jour la colonne.

Ici définissez d'une part votre table et la colonne mettre à jour, et d'autre part la table et la colonne qui contiennent vos codes insee.

Choisissez une jointure par objets géographiques en définissant que les objets de votre table de communes contiennent les objets de votre table points.

Et voila vous avez pour chaque point le code insee de la commune dans laquelle il est situé.

Pierre Vivet

Catégories: Sites Francophones

[Données] Re: CORINE LAND COVER - Prochaine mise à disposition ?

jeu 02-07-2015
Nombre de message(s) précédent(s) : 14

pour la fin du premier semestre 2015.

Si mes calculs sont bons, ça y est, c'est le début du second semestre 2015. ;-)
Mais pas encore la France sur http://land.copernicus.eu/pan-european/ … -2012/view
Quelqu'un a des news ?

Sylvain M.

Catégories: Sites Francophones

[Données] Re: Fichiers MAJIC 2015

jeu 02-07-2015
Nombre de message(s) précédent(s) : 6

Bonjour,

Quand je disais en début d'année, je pensais courant 1er trimestre.
Les données sont en effet mises à jour jusqu'en Avril Mai (date D4) car on comprend bien que l'on ne peut pas saisir toutes les informations pour le 1er Janvier au soir.
D'ailleurs la DGFIP communique une estimation des bases aux communes en début d'année (date D2).
Donc lorsque les fichiers fonciers annualisés dit Majic sont prêts en début d'été, la situation saisie est bien celle du 1er Janvier et pour les communes ayant un parcellaire dynamique avoir le plan en situation Juin ou Septembre N implique que des parcelles soient devenus obsolètes ou que de nouvelles se soient créées entre temps.
Les actes notariés doivent être enregistrés aux services de la publicité foncière avant un mois après la signature. Il faut tenir compte aussi du délai d'enregistrement de ces services. Cela peut aller de 1à 6 mois!!
Le lien plan-fichier se faisant généralement par le numéro de parcelle merci de m'indiquer si j'ai oublié quelque chose dans mon raisonnement

Cordialement

Catégories: Sites Francophones