Le blog de Jean David TECHER, un Réunionnais à Saint-Priest/Lyon

Aller au contenu | Aller au menu | Aller à la recherche


< 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 217 218 219 220 221 222 223 224 225 226 227 228 229 230 231 232 233 234 235 236 237 238 239 240 241 242 243 244 245 246 247 248 249 250 251 252 253 254 255 256 257 258 259 260 261 262 263 264 265 266 267 268 269 270 271 272 273 274 275 276 277 278 279 280 281 282 283 284 285 286 287 288 289 290 291 292 293 294 >

mercredi 16 mai 2007

Faire un peu joujou avec GoogleEarth

Ayé j'ai fait mon premier fichier kml personnel - hors contexte professionnel bien sûr -. Il est disponible à http://www.davidgis.fr/download/jeandavidtecher.kml

Le bâtiment en question représente le bâtiment où je vis à savoir la Résidence Prado Lido. Le polygone qui est dessiné a pour représentation dans le système de projection WGS 84 pour une hauteur de 12 mètres -un peu ficitf je l'avoue , je ne connais pas vraiment la hauteur de mon bâtiment -:

SELECT GeomFromEWKT('SRID=4326;POLYGON((3.896026111319043 43.63077293064538 12,3.896125412083016 43.63129020374376 12,
3.896017722403089 43.63130336362247 12,3.896023316938252 43.6313529652266 12,
3.895547802871119 43.63140560406845 12,3.895538013117296 43.63135701466003 12,
3.895430322945158 43.63136612505183 12,3.895325433631279 43.63084986305638 12,
3.895428927275387 43.6308407529833 12,3.895543607566332 43.63132462175059 12,
3.896003736551513 43.63126793394599 12,3.895911429755532 43.63078305367532 12,
3.896026111319043 43.63077293064538 12))');

lundi 14 mai 2007

Samedi soir - barbecue et petit soirée devant Oggy et les cafards

Samedi soir avec Madame ca a été soirée barbecue. On a profité qu'il fasse bon en fin d'après-midi pour sortir le barbeuk' et le mettre sur le balcon. On a grillé quelques belle côtes de viandes. Une petite salade a accompagné notre belle viande. Madame avait aussi conconcté de belles saucisses blanches de Lorraine. Hum un régal !!!

Comme tous les samedis - comme on ne change pas une équipe qui gagne - vers 20h30, j'ai squatté le clic-clac de ma salle informatique, ma télécommande ATI Remote Wonder II en main. J'allume Tvtime et zou direction gulli pour regarder «Oggy et les cafards»


Oggy et les cafards: dans cet épisode, un taureau qui va lui faire voir de toutes les couleurs.

Comme toujours, tordu de rire j'ai été! Vers 22h30, j'ai ensuite zappé pour regarder «Ken le survivant» sur Mangas. Finalement je me suis endormi devant la télé après en regardant les infos sur LCI.

samedi 12 mai 2007

QGIS 0.8.0 : test de l'exportateur de projet vers mapfile pour MapServer

Récemment, je me suis aperçu que par rapport à mon premier billet sur l'installation de QGIS 0.8.0, l'option d'exportation de QGIS ne marchait pas! Par «exportation de QGIS », j'entends

fichier de projet de QGIS en .qgs --[Exportation]--> fichier .map pour MapServer

En effet dans le sous-répertoire tools/mapserver_export/, la compilation ne passait pas.

Or celui-ci a besoin des headers de Python 2.4 (chez moi) pour être compilé. La solution a donc été

apt-get install python2.4-dev
puis comme cette option de Qgis nécessite python, c'est surtout l'option --with-python qu'il faut passer à configure pour l'activer
 ./configure --with-qtdir=/usr --with-wfs --with-python --with-geos=/usr/local/bin/geos-config --with-projdir=/usr/local \
--with-gdal=/opt/gdal-1.4.0/bin/gdal-config --with-postgresql=/usr/local/pgsql/bin
Et ensuite comme d'habitude
make
make install

Pour tester ensuite l'option, je me suis ensuite connecté au server WMS de GEOSIGNAL que j'ai chargé dans QGIS. Conformément à la directive de licence en cours sur le site de cette société, je rappelle ici l'URL de la licence du cadre d'utilisation. Je me suis amusé à zoomersur ma propre ville à savoir Castelnau-Le-Lez. En utilisant l'exportateur, j'ai eu

vendredi 11 mai 2007

Errors with PostGIS and POLYGON

Here I will try to show two examples of error with PostGIS and POLYGON.

1. Self-Intersection
The follwing POLYGON is not valid.
testgis=# SELECT IsValid('POLYGON((402 253,401 434,675 435,675 314,487 314,559 314,559 201,404 253,402 253))'::geometry);
INFO:  Self-intersection
INFO:  Self-intersection
 isvalid
---------
 f
He intersects hilmself. Have a look on the following picture!
2. Ring Self-Intersection
Suppose that you have a the following POLYGON
POLYGON ((380 340, 40 340, 40 20, 380 20, 380 340),(120 300, 300 280, 320 200, 
160 140, 200 80, 320 120, 320 200, 360 60, 120 40, 120 300))
Let's try Summary() function.
testgis=# SELECT summary('POLYGON ((380 340, 40 340, 40 20, 380 20, 380 340),
(120 300, 300 280, 320 200, 160 140, 200 80, 320 120, 320 200, 360 60, 120 40, 120 300))'::geometry);
         summary
-------------------------

 Polygon[B] with 2 rings
    ring 0 has 5 points
    ring 1 has 10 points

(1 ligne)
With
  • ring 0 = [380 340, 40 340, 40 20, 380 20, 380 340] --> My POLYGON
  • ring 1 = [120 300, 300 280, 320 200, 160 140, 200 80, 320 120, 320 200, 360 60, 120 40, 120 300] --> the hole of my polygon
If I try IsValid() function with PostGIS, I get
testgis=# SELECT IsValid('POLYGON ((380 340, 40 340, 40 20, 380 20, 380 340),
(120 300, 300 280, 320 200, 160 140, 200 80, 320 120, 320 200, 360 60, 120 40, 120 300))'::geometry);

INFO:  Ring Self-intersection
INFO:  Ring Self-intersection
 isvalid
---------
 f
(1 ligne)
That's not correct according to OGC Specification. A Polygon with a hole (ring 1) that self-intersects (i.e. POINT (320 200) )! Remember that it is just a case! This example has been taken from Geos SVN repository at ./tests/xmltester/TestValid.xml