Re: [OSGeo-Discuss] Mapnik in windows with gdal with ecw support
On Nov 11, 2009, at 9:24 AM, Pedro-Juan Ferrer Matoses wrote: Hello Frank. Thank you for your quick response. I have gdal 1.6 running properly (I mean with ecw support) in my sistem. And I've tried the OSGeo4w version and the official windows version. With both, I have no problems managing ecw images, for example using gdalinfo with my ecw images. The problem seems to be with Mapnik and Mapnik's gdal plug in. It seems that the gdal plug in which comes with mapnik in doesn't support ecw... Pedro, That is correct, the Mapnik binaries provided for windows do not include the gdal plugin with support for formats such as ecw or MrSid. You would need to recompile the library that is the 'gdal plugin' (with is separate from Mapnik's core library) against a gdal version built with that support. Dane ___ Discuss mailing list Discuss@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/discuss
Re: [OSGeo-Discuss] Mapnik in windows with gdal with ecw support
Thanks Mateusz, fixed. - dane (pronounced dayne, or danee, or dude, what does your name stand for?) On Nov 11, 2009, at 10:45 AM, Mateusz Loskot wrote: Dane Springmeyer wrote: On Nov 11, 2009, at 9:24 AM, Pedro-Juan Ferrer Matoses wrote: Hello Frank. Thank you for your quick response. I have gdal 1.6 running properly (I mean with ecw support) in my sistem. And I've tried the OSGeo4w version and the official windows version. With both, I have no problems managing ecw images, for example using gdalinfo with my ecw images. The problem seems to be with Mapnik and Mapnik's gdal plug in. It seems that the gdal plug in which comes with mapnik in doesn't support ecw... Pedro, That is correct, the Mapnik binaries provided for windows do not include the gdal plugin with support for formats such as ecw or MrSid. You would need to recompile the library that is the 'gdal plugin' (with is separate from Mapnik's core library) against a gdal version built with that support. Dane, GDAL (pronounced Goodle) [1] Someone wants Google to fine GDAL for too similar pronunciation of project name ;-) Anyway, there is not only one way correct [2] [1] http://trac.mapnik.org/wiki/MapnikDependencies [2] http://trac.osgeo.org/gdal/wiki/FAQGeneral#WhatdoesGDALstandsfor Best regards, -- Mateusz Loskot, http://mateusz.loskot.net ___ Discuss mailing list Discuss@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/discuss ___ Discuss mailing list Discuss@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/discuss
Re: [OSGeo-Discuss] The Easy steps to get your project on the FOSS4G Live DVD]
On Aug 29, 2009, at 12:34 PM, Tyler Mitchell (OSGeo) wrote: http://wiki.osgeo.org/wiki/FOSS4G_2009_Press_Release_28 To Geospatial Open Source Project Leaders: Tyler, Thanks for the heads up on getting projects included into the live dvd. I've got Mapnik covered and have just submitted my install script to the live-d...@lists.osgeo.org: http://lists.osgeo.org/pipermail/live-demo/2009-August/000205.html Thats for everyone's work on this to make it so easy. Cheers, Dane ___ Discuss mailing list Discuss@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/discuss
Re: [OSGeo-Discuss] Call for help: GDAL and txt2shp.py
Hi Kjell, Not having worked much with GIS systems or projection/coordinate systems before, I am struggling to understand the meaning of all this. However, the following (esriwkt) was the result of the visit to the link above, and I am injected it into the .prj file before running the script all over again: You are not alone. Projected coordinate systems can confuse beginners and experts alike. My understanding is that you have gps data and you are not sure what coordinate system it is in. You ran Matt Perry's script to convert the text-based coordinates into a shapefile format, but you still need to assign a coordinate system. However, the following (esriwkt) was the result of the visit to the link above, and I am injected it into the .prj file before running the script all over again: PROJCS[Nahrwan 1967 / UTM zone 39N,GEOGCS[Nahrwan 1967,DATUM[D_Nahrwan_1967,SPHEROID[Clarke_1880_RGS, 6378249.145,293.465]],PRIMEM[Greenwich,0],UNIT[Degree, 0.017453292519943295 ]],PROJECTION[Transverse_Mercator],PARAMETER[latitude_of_origin, 0],PARAMETER[central_meridian,51],PARAMETER[scale_factor, 0.9996],PARAMETER[false_easting,50],PARAMETER[false_northing, 0],UNIT[Meter,1]] Please let me know if I got this completely wrong and please excuse my lack of knowledge here. First, you only need to assign the projection AFTER you run the txt2shp.py script. Second, I can see now (from your image links) that your data was collected in Lebanon, which means that the link I directed you to will not provide a reference to the correct UTM zone. There are two possible UTM zones for Lebanon: UTM 36N or UTM 37N. WGS 84 / UTM zone 36N http://spatialreference.org/ref/epsg/32636/ WGS 84 / UTM zone 37N http://spatialreference.org/ref/epsg/32637/ You can check for yourself by downloading a world borders shapefile (http://thematicmapping.org/downloads/world_borders.php ) and a world UTM zone shapefile (http://earth-info.nga.mil/GandG/coordsys/gislayers/zips/mgrs6x8.zip ), both in the WGS 84/EPSG 4326 coordinate system. If you have more GPS data from other countries I highly recommend getting familiar with UTM zones. It did. Garmin GPSMap 60 CS. But, still, even with the new projection file, this thing still looks like this: http://www.ia-stud.hiof.no/~kjellare/misc/gdal_fail3_lebanon.png I even tried reversing the coordinates in the input file, but with almost the same result: http://www.ia-stud.hiof.no/~kjellare/misc/gdal_fail3_lebanon2.png Hm... Not sure how to proceed from here. If you try to load data in different coordinate systems into the same QGIS project, QGIS will not (by default) make an effort to 'reproject- on-the-fly' when rendering, unless you explicitly set the Qgis project to do so. So even if you assign the correct UTM zone to your points, if you open them in Qgis along with a base layer in WGS 84 projection, they will not line up. You need to go to SETTINGS PROJECT PROPERTIES PROJECTION ENABLE ON THE FLY PROJECTION. Or use the command line program ogr2ogr to reproject your shapefile to WGS84/EPSG 4326 from your assumed UTM projection. That command would look like: $ ogr2ogr -s_srs EPSG:32636 -t_srs EPSG:4326 lebanon_points_wgs84.shp lebanon_points_utm36N.shp Then, if the assumed source projection was correct, your new shapefile should line up with other data in WGS 84. ... In Qgis you should be able to reproject to UTM and label your points. Qgis will also export to PNG format if that is what you want. However I don't think that Qgis will output an ESRI world file (the only way that I know to keep a PNG spatially reference... but the application uDIG will output a .wld world file to spatially reference a PNG). I am in need of a open source command-line style software here, in order to automate this entire endevour. Both uDIG and QGIS appears to be GUI style applications. Are there a way around this to let me export the shapefiles I have to a .png raster file to complete my mission? Kjell Are Yes, they are GUI applications... but don't you just need to automate the processing of your GPS data? If you want to create one PNG raster then QGIS or uDIG are the ideal tools to layout your png map. Just open all your shapefiles (you can even merge them all into one shapefile with ogr2ogr), label them, then zoom to each group and export a PNG file If you truly want to automate the creation of *many* PNG graphics then you'll need to look into scripting a mapping toolkit like Mapserver or Mapnik. Dane ___ Discuss mailing list Discuss@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/discuss ___ Discuss mailing list Discuss@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/discuss
Re: [OSGeo-Discuss] Call for help: GDAL and txt2shp.py
Refsvik, Try reversing your coordinates, as the graphic looks like your output shapefile has had the x/y reversed. Then if you want to display the data against other layers not in WGS 84/EPSG 4326 you'll need to add a .prj file, which can be grabbed from http://spatialreference.org/ref/epsg/4326/ Dane On Jul 3, 2008, at 12:51 AM, Kjell Are Refsvik wrote: Hi, I am a graduate student in Information Technology at Ostfold University College - Dept. of Computer Science. Nearing the end of my thesis on the challenges and opportunities related to digital photographs I would like to write some code that process pictures and makes a HTML web album (with an embedded map with a numbered legend on it) using open source tools. Inspired by GeoPerrys note on command-line style approach to using software (http://www.perrygeo.net/wordpress/?p=105) I have been looking for a basic unix toolset to process geodata. My geo-data come in the form of a pile of geo-tagged images and my own maps (shapefiles), and going back to GeoPerry, I found a python script: http://www.perrygeo.net/wordpress/?p=5 ...that appears to do more or less what I want. I should mention that my I have little experience in programming in general and that diving deep into the terminal on my Mac OS X installation feels a bit daunting. Even so, I was able to upgrade my Python-installation yesterday and install GDAL and all its dependencies (I think) to get the script (txt2shp.py) to work. When I run the script: http://www.ia-stud.hiof.no/~kjellare/misc/Terminal.png ...using my data as input: http://www.ia-stud.hiof.no/~kjellare/misc/Norway_map.zip http://www.ia-stud.hiof.no/~kjellare/misc/input.txt ...and view the output shapefile that comes out of it, I observe this strange offset between the map (or Norway) and my data (geodata from images taken in Norway) and know something must be wrong. http://www.ia-stud.hiof.no/~kjellare/misc/QGIS.png Not knowing a whole lot about GIS, I would really appreciate help on how to proceed from here to: 1. Work out the reason behind, and how to fix the offset problem 2. Change the map projection (preferably to Mercator) and have the legend dots numbered 3. Output a png of a mesh between the map and the legend for web use Best regards, Kjell Are Refsvik Graduate student Ostfold University College Norway ___ Discuss mailing list Discuss@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/discuss ___ Discuss mailing list Discuss@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/discuss
Re: [OSGeo-Discuss] Call for help: GDAL and txt2shp.py
On Jul 3, 2008, at 1:40 PM, Kjell Are Refsvik wrote: On 3. juli. 2008, at 21.36, Dane Springmeyer wrote: Refsvik, Try reversing your coordinates, as the graphic looks like your output shapefile has had the x/y reversed. Splendid! That did the trick. great. However - I have another dataset from the middle east, and approaching this data the same way, even trying to reverse the dataset from the images, it still fails: http://www.ia-stud.hiof.no/~kjellare/misc/gdal_fail1.png http://www.ia-stud.hiof.no/~kjellare/misc/gdal_fail2.png I can see from the Qgis mouse coordinates in the lower right that the shapefile you created likely has x,y coordinates in a projected coordinate system (at least different from WGS 84/4326). My first (wild) guess would be to add a .prj file from spatialreference.org for a UTM zone that covers your area of the middle east (you'll be need to find out what exact zone number is needed if it is in UTM). I would try: http://www.google.com/search?hl=enq=site%3Ahttp%3A%2F%2Fspatialreference.org+utm+middle+eastbtnG=Search UTM is a likely candidate if the data came from a GPS unit, but its just a guess. ... trying to have the dots numbered or named and outputting the two layers as a transverse mercator projected .png file are the two things left to do. I would really appreciate some input on how to proceed, as I cannot figure out if gdal_rasterize can be used for this purpose. In Qgis you should be able to reproject to UTM and label your points. Qgis will also export to PNG format if that is what you want. However I don't think that Qgis will output an ESRI world file (the only way that I know to keep a PNG spatially reference... but the application uDIG will output a .wld world file to spatially reference a PNG). I'd recommend #qgis or #udig with further questions. Dane Best regards, Kjell Are ___ Discuss mailing list Discuss@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/discuss ___ Discuss mailing list Discuss@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/discuss
