
<html><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><br><div><div>On 3. juli. 2008, at 23.02, Dane Springmeyer wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite"><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><div>...</div><div><blockquote type="cite"><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><div><div><span class="Apple-style-span" style="-webkit-text-stroke-width: -1; ">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:</span></div><div><br></div><div><a href="http://www.ia-stud.hiof.no/~kjellare/misc/gdal_fail1.png">http://www.ia-stud.hiof.no/~kjellare/misc/gdal_fail1.png</a></div><div><a href="http://www.ia-stud.hiof.no/~kjellare/misc/gdal_fail2.png">http://www.ia-stud.hiof.no/~kjellare/misc/gdal_fail2.png</a></div><div><br></div></div></div></blockquote><div><br></div><div>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). </div><div><br></div><div>I would try: </div><div><br></div><div><a href="http://www.google.com/search?hl=en&q=site%3Ahttp%3A%2F%2Fspatialreference.org+utm+middle+east&btnG=Search">http://www.google.com/search?hl=en&q=site%3Ahttp%3A%2F%2Fspatialreference.org+utm+middle+east&btnG=Search</a></div></div></div></blockquote><div><br></div><div>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:</div><div><span class="Apple-style-span" style="font-family: Times; font-size: 16px; "><pre style="word-wrap: break-word; white-space: pre-wrap; "><font class="Apple-style-span" face="Monaco" size="3"><span class="Apple-style-span" style="font-size: 12px;">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",500000],PARAMETER["false_northing",0],UNIT["Meter",1]]</span></font></pre></span></div><div>Please let me know if I got this completely wrong and please excuse my lack of knowledge here.</div><div><br></div><blockquote type="cite"><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><div><div><span class="Apple-style-span" style="-webkit-text-stroke-width: -1; ">UTM is a likely candidate if the data came from a GPS unit, but its just a guess.</span></div></div></div></blockquote><div><br></div><div>It did. Garmin GPSMap 60 CS.</div><div><br></div><div>But, still, even with the new projection file, this thing still looks like this:</div><div><br></div><div><font class="Apple-style-span" color="#144FAE"><span class="Apple-style-span" style="text-decoration: underline;"><a href="http://www.ia-stud.hiof.no/~kjellare/misc/gdal_fail3_lebanon.png">http://www.ia-stud.hiof.no/~kjellare/misc/gdal_fail3_lebanon.png</a></span></font></div><div><br></div><div>I even tried reversing the coordinates in the input file, but with almost the same result:</div><div><br></div><div><span class="Apple-style-span" style="color: rgb(20, 79, 174); "><span class="Apple-style-span" style="text-decoration: underline; -webkit-text-decorations-in-effect: underline; "><a href="http://www.ia-stud.hiof.no/~kjellare/misc/gdal_fail3_lebanon2.png">http://www.ia-stud.hiof.no/~kjellare/misc/gdal_fail3_lebanon2.png</a></span></span></div><div><br></div><div>Hm... Not sure how to proceed from here.</div><br><blockquote type="cite"><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><div><div>...</div>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).</div></div></blockquote><div><br></div><div>I am in need of a open source command-line style software here, in order to automate this entire endevour. </div><div>Both uDIG and QGIS appears to be GUI style applications. </div><div><br></div><div>Are there a way around this to let me export the shapefiles I have to a .png raster file to complete my mission?</div><div><br></div></div>Kjell Are<div><br></div></body></html>