<div class="gmail_quote">On 11 March 2013 00:54, Alex Mandel <span dir="ltr"><<a href="mailto:tech_dev@wildintellect.com" target="_blank">tech_dev@wildintellect.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div class="im">On 03/10/2013 03:49 PM, Lester Anderson wrote:<br>
> On 10 March 2013 22:22, Alex Mandel <<a href="mailto:tech_dev@wildintellect.com">tech_dev@wildintellect.com</a>> wrote:<br>
><br>
>> On 03/10/2013 03:13 PM, Lester Anderson wrote:<br>
>><br>
>>> Hello,<br>
>>><br>
>>> The one major element missing from the otherwise great Quantum GIS, is<br>
>>> that<br>
>>> of easy setup map graticules. There is a basic way of doing simple ones in<br>
>>> the print composer which is fine for geographic (WGS84) or UTM etc type<br>
>>> projections, but will not work for conic or stereographic etc.<br>
>>><br>
>>> Is there a plugin for doing this work (that is built-in to layout mode in<br>
>>> ArcGIS) or is this an upgrade/update that will be available in a new<br>
>>> release?<br>
>>><br>
>>> Generating a vector grid does not work properly in 1.8.0, and certainly<br>
>>> only partially completes a Polar Stereographic layout. There are clearly<br>
>>> flaws that need to be addressed. If there is a reliable workaround for<br>
>>> this, it would be good to know.<br>
>>><br>
>>> Lester<br>
>>><br>
>>><br>
>> A workaround is to make the graticule in WGS84 and reproject it to the<br>
>> desired end projection. I found that it doesn't curve well so I wrote this<br>
>> script to make graticules with a higher points density that projects well.<br>
>><br>
</div>>> <a href="https://github.com/**wildintellect/pyGraticule" target="_blank">https://github.com/**wildintellect/pyGraticule</a><<a href="https://github.com/wildintellect/pyGraticule" target="_blank">https://github.com/wildintellect/pyGraticule</a>><br>
<div class="im">>><br>
>> I hope to work it into QGIS at a later point, but you can use it<br>
>> standalone with python to make what you need, import and reproject.<br>
>><br>
>> Thanks,<br>
>> Alex<br>
>><br>
>><br>
> Hi Alex,<br>
><br>
> I have tried the route you suggested. I am working on Antarctic data, so<br>
> generated a vector grid for 0-360 in X (at 10 degrees) and latitude (-90 to<br>
> -60) in 10 degrees. That all worked fine in WGS84. However, reprojecting to<br>
> Antarctic Polarstreographic did not project correctly with less than half<br>
> the meridians and no latitudes.<br>
><br>
> Cheers<br>
> Lester<br>
><br>
<br>
</div>This will probably get stripped from the list but I've attached a shp I<br>
made that should be what you need. Let me know if you need me to upload<br>
it to server for download.<br>
<br>
Note projection on the fly has some odd bug that makes 180/-180<br>
disappear, I've not spent enough time to figure out what that is.<br>
<br>
What I did:<br>
Used my pyGraticule script modified to go from -90,1<br>
Reprojected to Polar EPSG:3031 as a shapefile<br>
<br>
I'm trying to remember how I did<br>
<a href="http://geography.ucdavis.edu/image/499" target="_blank">http://geography.ucdavis.edu/image/499</a><br>
If I remember I'll post some more details.<br>
<br>
Thanks,<br>
Alex<br>
<br>
<br>
</blockquote></div>Hi Alex,<br><br>That worked very well. I had a look at the Python script and I note that it generates a GeoJSON format file. How do you convert to a shapefile?<br><br>Hopefully something along the lines of the script will be built in to QGis, the OTF reprojection of somethings does fall over.<br>
<br>Cheers<br>Lester<br>