Re: [Qgis-user] Projection of Shapefile export

2011-03-23 Thread Brett Adams 

Ramon,
Loaded your shape files into ARCGIS. It recognised none of the 
projections defaulting to WGS84 every time.


Attached are a couple of prj files generated by ARCGIS. These look much 
nicer.


With any luck, all these issues will be sorted in the major release of 
1.7 due at the end of this month.:-) I'm looking forward to this one.


BA



On 24/03/2011 6:48 AM, Ramon Andinach wrote:

I had a look at them when I was trying to work out why some imported text 
wasn't correctly reporting its CRS. They look like pretty strange to me, but 
since I could open the project on a linux and win7 in different computers with 
identical results I was happy at that.

It's unlikely I'll get near MI or Arc for quite some time, but I'll be curious 
what you see.

-ramon.


On 23/03/2011, at 23:38 , Brett Adams wrote:


Ramon,
I loaded the gsb files as you suggested. The grids now load into the correct 
position (correct enough for me anyway). Thats a great achievement right there.

The exported projection files (prj) look a bit odd. The segment 
DATUM["D_unknown" makes me a little nervous, however as I don't really know 
what I'm looking at, its hard to judge.  I'll need to see if they are recognised by 
Mapinfo or Arc.


Brett.



On 23/03/2011 12:58 PM, custard wrote:

Brett,

If you grab this zip file:
http://members.westnet.com.au/ramon_et_al/BrettTest.zip

In this zip file is your csv, a bundle of derived shape files, a 1.7 QGIS 
project (which will probably make 1.6 crash - but the shape files are good) and 
a folder with a few images.

The images are from the project, with projection set to MGA94 zone 51, but I 
have switched between GDA94, ADG66/84, WGS84 and the equivilent zone 51s, with 
OTF on.

The images show the points with OTF on and off, and zoomed in.

Are these close enough for your purposes?

The folder I gave for standalone should be right, but this method (for reasons 
I don't understand) doesn't seem to work on 1.6 standalone under WinXP. (I'm 
hoping it's a 1.6 thing, and not an osgeos/standalone thing)

-ramon.

- Original Message -
From: "Brett Adams"
To: qgis-user@lists.osgeo.org, "Ramon Andinach"
Sent: Wednesday, 23 March, 2011 8:20:43 AM GMT +08:00 Beijing / Chongqing / 
Hong Kong / Urumqi
Subject: Re: [Qgis-user]  Projection of Shapefile export

Ramon,
Attached CSV contains grid coordinates in 4 different projections; AGD66, 
AGD84, GDA94 and WGS84. All zone 51 and 51S

The conversions were calculated using Mapinfo. Conversions can vary up to 1-2m 
depending on the parameters used so don't be too surprised if you can't 
replicate these exactly.

These are the same points attached to Bug Ticket 3099.

Brett



On 22/03/2011 6:37 PM, Ramon Andinach wrote:
On 21/03/2011, at 19:51 , Ramon Andinach wrote:


On 21/03/2011, at 10:21 , Brett Adams wrote:


I thought this problem merely related to the failure of Australian projections. 
Projection data doesn't exist, therefore can't be exported with the shape file.

If the points are imported as WGS84 / UTMz51S, then export the shape file, the 
projection file contains the relevant information. Its all good.

So if the projection issue is fixed, shape files will then export with the 
required projection information.


As far as I can see, there is no bug here. If anyone can show otherwise, we 
need to record the bug and create a ticket (if it doesn't exist already)

The little bit of information I had before, I just tacked on to the end of your 
existing ticket.

I was going to put a bit of my spare time during this week into trying to 
replicate Roland's solution firstly on my home linux machine and then see if I 
could make that approach work on windows at work.

> From my recollection last time I was trying Ken's solution and a few others 
like just defining a CRS and couldn't quite get there.

Let you know what eventuates.

Short version:
I'm fairly sure that QGIS 1.7 is approximately 2 (fairly large) files away from 
being able to handle AGD points.

The files are here:

http://www.icsm.gov.au/icsm/gda/gdatm/national66.zip

and here:

http://www.icsm.gov.au/icsm/gda/gdatm/national84.zip


and these equate to the missing datum shift files.

In 1.7 as at this afternoon, you need to download them, unzip them rename them 
to agd66.gsb or agd84.gsb and copy them to:
/usr/share/proj  - for linux users
C:/OSGeo4W/share/proj - for OSGeo windows installs

and that means it should be something like
C:/Program Files/Quantum version/share/proj
for standalone versions (sorry forgot to check before going home - will check 
in the morning).

Paolo, I know string freeze is tomorrow, but is there any chance that these two 
files could be included in 1.7?
They should be under a CC by Atrrib (Australian variant).
(I can ring up GA tomorrow and confirm if needed. Do any of the Australians on 
the list know who migh

Re: [Qgis-user] Projection of Shapefile export

2011-03-23 Thread Ramon Andinach 
I had a look at them when I was trying to work out why some imported text 
wasn't correctly reporting its CRS. They look like pretty strange to me, but 
since I could open the project on a linux and win7 in different computers with 
identical results I was happy at that.

It's unlikely I'll get near MI or Arc for quite some time, but I'll be curious 
what you see.

-ramon.


On 23/03/2011, at 23:38 , Brett Adams wrote:

> Ramon, 
> I loaded the gsb files as you suggested. The grids now load into the correct 
> position (correct enough for me anyway). Thats a great achievement right 
> there.
> 
> The exported projection files (prj) look a bit odd. The segment 
> DATUM["D_unknown" makes me a little nervous, however as I don't really know 
> what I'm looking at, its hard to judge.  I'll need to see if they are 
> recognised by Mapinfo or Arc. 
> 
> 
> Brett.
> 
> 
> 
> On 23/03/2011 12:58 PM, custard wrote:
>> Brett, 
>> 
>> If you grab this zip file:
>> http://members.westnet.com.au/ramon_et_al/BrettTest.zip
>> 
>> In this zip file is your csv, a bundle of derived shape files, a 1.7 QGIS 
>> project (which will probably make 1.6 crash - but the shape files are good) 
>> and a folder with a few images.
>> 
>> The images are from the project, with projection set to MGA94 zone 51, but I 
>> have switched between GDA94, ADG66/84, WGS84 and the equivilent zone 51s, 
>> with OTF on.
>> 
>> The images show the points with OTF on and off, and zoomed in.
>> 
>> Are these close enough for your purposes?
>> 
>> The folder I gave for standalone should be right, but this method (for 
>> reasons I don't understand) doesn't seem to work on 1.6 standalone under 
>> WinXP. (I'm hoping it's a 1.6 thing, and not an osgeos/standalone thing)
>> 
>> -ramon.
>> 
>> ----- Original Message -
>> From: "Brett Adams" 
>> To: qgis-user@lists.osgeo.org, "Ramon Andinach" 
>> Sent: Wednesday, 23 March, 2011 8:20:43 AM GMT +08:00 Beijing / Chongqing / 
>> Hong Kong / Urumqi
>> Subject: Re: [Qgis-user]  Projection of Shapefile export
>> 
>> Ramon,
>> Attached CSV contains grid coordinates in 4 different projections; AGD66, 
>> AGD84, GDA94 and WGS84. All zone 51 and 51S
>> 
>> The conversions were calculated using Mapinfo. Conversions can vary up to 
>> 1-2m depending on the parameters used so don't be too surprised if you can't 
>> replicate these exactly. 
>> 
>> These are the same points attached to Bug Ticket 3099.
>> 
>> Brett
>> 
>> 
>> 
>> On 22/03/2011 6:37 PM, Ramon Andinach wrote:
>> On 21/03/2011, at 19:51 , Ramon Andinach wrote:
>> 
>> 
>> On 21/03/2011, at 10:21 , Brett Adams wrote:
>> 
>> 
>> I thought this problem merely related to the failure of Australian 
>> projections. Projection data doesn't exist, therefore can't be exported with 
>> the shape file.
>> 
>> If the points are imported as WGS84 / UTMz51S, then export the shape file, 
>> the projection file contains the relevant information. Its all good.
>> 
>> So if the projection issue is fixed, shape files will then export with the 
>> required projection information.
>> 
>> 
>> As far as I can see, there is no bug here. If anyone can show otherwise, we 
>> need to record the bug and create a ticket (if it doesn't exist already)
>> 
>> The little bit of information I had before, I just tacked on to the end of 
>> your existing ticket.
>> 
>> I was going to put a bit of my spare time during this week into trying to 
>> replicate Roland's solution firstly on my home linux machine and then see if 
>> I could make that approach work on windows at work.
>> 
>> >From my recollection last time I was trying Ken's solution and a few others 
>> >like just defining a CRS and couldn't quite get there.
>> 
>> Let you know what eventuates.
>> 
>> Short version:
>> I'm fairly sure that QGIS 1.7 is approximately 2 (fairly large) files away 
>> from being able to handle AGD points.
>> 
>> The files are here: 
>> 
>> http://www.icsm.gov.au/icsm/gda/gdatm/national66.zip
>> 
>> and here:
>> 
>> http://www.icsm.gov.au/icsm/gda/gdatm/national84.zip
>> 
>> 
>> and these equate to the missing datum shift files.
>> 
>> In 1.7 as at this afternoon, you need to download them, unzip them rename 
>> them to agd66.gsb or agd84.gsb and copy them to:
>>

Re: [Qgis-user] Projection of Shapefile export

2011-03-23 Thread Brett Adams 

Ramon,
I loaded the gsb files as you suggested. The grids now load into the 
correct position (correct enough for me anyway). Thats a great 
achievement right there.


The exported projection files (prj) look a bit odd. The segment 
DATUM["D_unknown" makes me a little nervous, however as I don't really 
know what I'm looking at, its hard to judge.  I'll need to see if they 
are recognised by Mapinfo or Arc.



Brett.



On 23/03/2011 12:58 PM, custard wrote:

Brett,

If you grab this zip file:
http://members.westnet.com.au/ramon_et_al/BrettTest.zip

In this zip file is your csv, a bundle of derived shape files, a 1.7 
QGIS project (which will probably make 1.6 crash - but the shape files 
are good) and a folder with a few images.


The images are from the project, with projection set to MGA94 zone 51, 
but I have switched between GDA94, ADG66/84, WGS84 and the equivilent 
zone 51s, with OTF on.


The images show the points with OTF on and off, and zoomed in.

Are these close enough for your purposes?

The folder I gave for standalone should be right, but this method (for 
reasons I don't understand) doesn't seem to work on 1.6 standalone 
under WinXP. (I'm hoping it's a 1.6 thing, and not an 
osgeos/standalone thing)


-ramon.

- Original Message -
From: "Brett Adams" 
To: qgis-user@lists.osgeo.org, "Ramon Andinach" 
Sent: Wednesday, 23 March, 2011 8:20:43 AM GMT +08:00 Beijing / 
Chongqing / Hong Kong / Urumqi

Subject: Re: [Qgis-user]  Projection of Shapefile export

Ramon,
Attached CSV contains grid coordinates in 4 different projections; 
AGD66, AGD84, GDA94 and WGS84. All zone 51 and 51S


The conversions were calculated using Mapinfo. Conversions can vary up 
to 1-2m depending on the parameters used so don't be too surprised if 
you can't replicate these exactly.


These are the same points attached to Bug Ticket 3099.

Brett



On 22/03/2011 6:37 PM, Ramon Andinach wrote:

On 21/03/2011, at 19:51 , Ramon Andinach wrote:


On 21/03/2011, at 10:21 , Brett Adams wrote:


I thought this problem merely related to the failure of Australian 
projections. Projection data doesn't exist, therefore can't be exported with 
the shape file.

If the points are imported as WGS84 / UTMz51S, then export the 
shape file, the projection file contains the relevant information. Its all good.

So if the projection issue is fixed, shape files will then export 
with the required projection information.


As far as I can see, there is no bug here. If anyone can show 
otherwise, we need to record the bug and create a ticket (if it doesn't exist 
already)

The little bit of information I had before, I just tacked on to the end 
of your existing ticket.

I was going to put a bit of my spare time during this week into trying 
to replicate Roland's solution firstly on my home linux machine and then see if 
I could make that approach work on windows at work.

> From my recollection last time I was trying Ken's solution and a few 
others like just defining a CRS and couldn't quite get there.

Let you know what eventuates.

Short version:
I'm fairly sure that QGIS 1.7 is approximately 2 (fairly large) files away 
from being able to handle AGD points.

The files are here:
http://www.icsm.gov.au/icsm/gda/gdatm/national66.zip
and here:
http://www.icsm.gov.au/icsm/gda/gdatm/national84.zip

and these equate to the missing datum shift files.

In 1.7 as at this afternoon, you need to download them, unzip them rename 
them to agd66.gsb or agd84.gsb and copy them to:
/usr/share/proj  - for linux users
C:/OSGeo4W/share/proj - for OSGeo windows installs

and that means it should be something like
C:/Program Files/Quantum version/share/proj
for standalone versions (sorry forgot to check before going home - will 
check in the morning).

Paolo, I know string freeze is tomorrow, but is there any chance that these 
two files could be included in 1.7?
They should be under a CC by Atrrib (Australian variant).
(I can ring up GA tomorrow and confirm if needed. Do any of the Australians 
on the list know who might be the right person to ask? Just to save me 
stumbling around in the dark a bit.)

Brett, I also had a problem with QGIS remembering CRS but I'll put that in 
a different thread. If you had any test files that you wanted me to check, I'd 
be happy to do so.

-ramon.



--
Brett Adams
Spinifex Geophysics <http://www.spinifexgeophysics.com.au/>
0438 861 974
SKYPE:brettadams_spinifex


___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user



--
Brett Adams
Spinifex Geophysics <http://www.spinifexgeophysics.com.au/>
0438 861 974
SKYPE:brettadams_s

Re: [Qgis-user] Projection of Shapefile export

2011-03-22 Thread custard 
Brett, 

If you grab this zip file: 
http://members.westnet.com.au/ramon_et_al/BrettTest.zip 

In this zip file is your csv, a bundle of derived shape files, a 1.7 QGIS 
project (which will probably make 1.6 crash - but the shape files are good) and 
a folder with a few images. 

The images are from the project, with projection set to MGA94 zone 51, but I 
have switched between GDA94, ADG66/84, WGS84 and the equivilent zone 51s, with 
OTF on. 

The images show the points with OTF on and off, and zoomed in. 

Are these close enough for your purposes? 

The folder I gave for standalone should be right, but this method (for reasons 
I don't understand) doesn't seem to work on 1.6 standalone under WinXP. (I'm 
hoping it's a 1.6 thing, and not an osgeos/standalone thing) 

-ramon. 

- Original Message - 
From: "Brett Adams"  
To: qgis-user@lists.osgeo.org, "Ramon Andinach"  
Sent: Wednesday, 23 March, 2011 8:20:43 AM GMT +08:00 Beijing / Chongqing / 
Hong Kong / Urumqi 
Subject: Re: [Qgis-user] Projection of Shapefile export 

Ramon, 
Attached CSV contains grid coordinates in 4 different projections; AGD66, 
AGD84, GDA94 and WGS84. All zone 51 and 51S 

The conversions were calculated using Mapinfo. Conversions can vary up to 1-2m 
depending on the parameters used so don't be too surprised if you can't 
replicate these exactly. 

These are the same points attached to Bug Ticket 3099. 

Brett 



On 22/03/2011 6:37 PM, Ramon Andinach wrote: 

On 21/03/2011, at 19:51 , Ramon Andinach wrote: 

On 21/03/2011, at 10:21 , Brett Adams wrote: 

I thought this problem merely related to the failure of Australian projections. 
Projection data doesn't exist, therefore can't be exported with the shape file.

If the points are imported as WGS84 / UTMz51S, then export the shape file, the 
projection file contains the relevant information. Its all good.

So if the projection issue is fixed, shape files will then export with the 
required projection information.


As far as I can see, there is no bug here. If anyone can show otherwise, we 
need to record the bug and create a ticket (if it doesn't exist already) The 
little bit of information I had before, I just tacked on to the end of your 
existing ticket.

I was going to put a bit of my spare time during this week into trying to 
replicate Roland's solution firstly on my home linux machine and then see if I 
could make that approach work on windows at work.

>From my recollection last time I was trying Ken's solution and a few others 
>like just defining a CRS and couldn't quite get there.

Let you know what eventuates. Short version:
I'm fairly sure that QGIS 1.7 is approximately 2 (fairly large) files away from 
being able to handle AGD points.

The files are here: http://www.icsm.gov.au/icsm/gda/gdatm/national66.zip and 
here: http://www.icsm.gov.au/icsm/gda/gdatm/national84.zip and these equate to 
the missing datum shift files.

In 1.7 as at this afternoon, you need to download them, unzip them rename them 
to agd66.gsb or agd84.gsb and copy them to:
/usr/share/proj  - for linux users
C:/OSGeo4W/share/proj - for OSGeo windows installs

and that means it should be something like
C:/Program Files/Quantum version/share/proj
for standalone versions (sorry forgot to check before going home - will check 
in the morning).

Paolo, I know string freeze is tomorrow, but is there any chance that these two 
files could be included in 1.7? 
They should be under a CC by Atrrib (Australian variant). 
(I can ring up GA tomorrow and confirm if needed. Do any of the Australians on 
the list know who might be the right person to ask? Just to save me stumbling 
around in the dark a bit.)

Brett, I also had a problem with QGIS remembering CRS but I'll put that in a 
different thread. If you had any test files that you wanted me to check, I'd be 
happy to do so.

-ramon. 


-- 
Brett Adams 
Spinifex Geophysics 
0438 861 974 
SKYPE:brettadams_spinifex ___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] Projection of Shapefile export

2011-03-22 Thread Brett Adams 

Ramon,
Attached CSV contains grid coordinates in 4 different projections; 
AGD66, AGD84, GDA94 and WGS84. All zone 51 and 51S


The conversions were calculated using Mapinfo. Conversions can vary up 
to 1-2m depending on the parameters used so don't be too surprised if 
you can't replicate these exactly.


These are the same points attached to Bug Ticket 3099.

Brett



On 22/03/2011 6:37 PM, Ramon Andinach wrote:

On 21/03/2011, at 19:51 , Ramon Andinach wrote:


On 21/03/2011, at 10:21 , Brett Adams wrote:


I thought this problem merely related to the failure of Australian projections. 
Projection data doesn't exist, therefore can't be exported with the shape file.

If the points are imported as WGS84 / UTMz51S, then export the shape file, the 
projection file contains the relevant information. Its all good.

So if the projection issue is fixed, shape files will then export with the 
required projection information.


As far as I can see, there is no bug here. If anyone can show otherwise, we 
need to record the bug and create a ticket (if it doesn't exist already)

The little bit of information I had before, I just tacked on to the end of your 
existing ticket.

I was going to put a bit of my spare time during this week into trying to 
replicate Roland's solution firstly on my home linux machine and then see if I 
could make that approach work on windows at work.

 From my recollection last time I was trying Ken's solution and a few others 
like just defining a CRS and couldn't quite get there.

Let you know what eventuates.

Short version:
I'm fairly sure that QGIS 1.7 is approximately 2 (fairly large) files away from 
being able to handle AGD points.

The files are here:
http://www.icsm.gov.au/icsm/gda/gdatm/national66.zip
and here:
http://www.icsm.gov.au/icsm/gda/gdatm/national84.zip

and these equate to the missing datum shift files.

In 1.7 as at this afternoon, you need to download them, unzip them rename them 
to agd66.gsb or agd84.gsb and copy them to:
/usr/share/proj  - for linux users
C:/OSGeo4W/share/proj - for OSGeo windows installs

and that means it should be something like
C:/Program Files/Quantum version/share/proj
for standalone versions (sorry forgot to check before going home - will check 
in the morning).

Paolo, I know string freeze is tomorrow, but is there any chance that these two 
files could be included in 1.7?
They should be under a CC by Atrrib (Australian variant).
(I can ring up GA tomorrow and confirm if needed. Do any of the Australians on 
the list know who might be the right person to ask? Just to save me stumbling 
around in the dark a bit.)

Brett, I also had a problem with QGIS remembering CRS but I'll put that in a 
different thread. If you had any test files that you wanted me to check, I'd be 
happy to do so.

-ramon.


The long version is below.
None of this is new - mostly it's that I'm slow to understand.

In this message, Ken Norris points to the relevent distortion files:
http://lists.osgeo.org/pipermail/qgis-user/2010-September/009731.html

Roland Hill picked up on this and modified the srs.db file and produced these 
messages:
http://www.mail-archive.com/qgis-user@lists.osgeo.org/msg07365.html and
http://www.mail-archive.com/qgis-user@lists.osgeo.org/msg07379.html

At least in the first message he used the following parameters to update the 
srs.db.
For AGD84:
'+ellps=aust_SA', '+a=6378160 +rf=298.2555 
+towgs84=-117.763,-51.51,139.061,-0.292,-0.443,-0.277,-0.191
and for AGD66:
'+ellps=aust_SA', '+a=6378160 +rf=298.2555 
+towgs84=-129.193,-41.212,130.73,-0.246,-0.374,-0.329,-2.955

Then placed the agd84 and agd66.gsb files into /usr/share/proj.

In attempting to replicate Roland I noticed that in 1.7 AGD84 is defined as 
ESPG:4203 and in Proj.4 format _projlonglat +elips=aust-SA 
+towards84=-117.763,-51.51,139.061,-0.443,-0.277,-0.191 +nadgrids=agd.gsb 
+wktext=no_defs which to me looked a lot like Roland's definition (please don't 
ask me to translate - I can't).

So I added in the two .gsb files and tested using this test set: 
http://members.westnet.com.au/ramon_et_al/SinglePoint2.zip

This contains a single point in my main area of operations, recorded originally 
in MGA94z51 using a Garmin GPS76. I've then converted the point into GDA94 
using GA's online converter, and into AGD66, AGD84 (and relevant UTM 
equivilants) using GA's agd-gda94.xls block calculations (see the included 
ASinglePoint.ods for calculations).

 From the spreadsheet I've loaded the point in each different projection 
(ASinglePoint.qgs).

OTF works, and the points in that project cluster within a 1.1mx0.6m area. I 
would suspect that most of that error is probably more a reflection of the 
rough methods I used to convert them.


___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user




--
Brett Adams
Spinifex Geophysics 
0438 861 974

Re: [Qgis-user] Projection of Shapefile export

2011-03-22 Thread Ramon Andinach 

On 21/03/2011, at 19:51 , Ramon Andinach wrote:

> 
> On 21/03/2011, at 10:21 , Brett Adams wrote:
> 
>> I thought this problem merely related to the failure of Australian 
>> projections. Projection data doesn't exist, therefore can't be exported with 
>> the shape file.
>> 
>> If the points are imported as WGS84 / UTMz51S, then export the shape file, 
>> the projection file contains the relevant information. Its all good.
>> 
>> So if the projection issue is fixed, shape files will then export with the 
>> required projection information.
>> 
>> 
>> As far as I can see, there is no bug here. If anyone can show otherwise, we 
>> need to record the bug and create a ticket (if it doesn't exist already)
> 
> The little bit of information I had before, I just tacked on to the end of 
> your existing ticket.
> 
> I was going to put a bit of my spare time during this week into trying to 
> replicate Roland's solution firstly on my home linux machine and then see if 
> I could make that approach work on windows at work.
> 
> From my recollection last time I was trying Ken's solution and a few others 
> like just defining a CRS and couldn't quite get there.
> 
> Let you know what eventuates.

Short version:
I'm fairly sure that QGIS 1.7 is approximately 2 (fairly large) files away from 
being able to handle AGD points.

The files are here: 
http://www.icsm.gov.au/icsm/gda/gdatm/national66.zip
and here:
http://www.icsm.gov.au/icsm/gda/gdatm/national84.zip

and these equate to the missing datum shift files.

In 1.7 as at this afternoon, you need to download them, unzip them rename them 
to agd66.gsb or agd84.gsb and copy them to:
/usr/share/proj  - for linux users
C:/OSGeo4W/share/proj - for OSGeo windows installs

and that means it should be something like
C:/Program Files/Quantum version/share/proj
for standalone versions (sorry forgot to check before going home - will check 
in the morning).

Paolo, I know string freeze is tomorrow, but is there any chance that these two 
files could be included in 1.7? 
They should be under a CC by Atrrib (Australian variant). 
(I can ring up GA tomorrow and confirm if needed. Do any of the Australians on 
the list know who might be the right person to ask? Just to save me stumbling 
around in the dark a bit.)

Brett, I also had a problem with QGIS remembering CRS but I'll put that in a 
different thread. If you had any test files that you wanted me to check, I'd be 
happy to do so.

-ramon.


The long version is below.
None of this is new - mostly it's that I'm slow to understand.

In this message, Ken Norris points to the relevent distortion files:
http://lists.osgeo.org/pipermail/qgis-user/2010-September/009731.html

Roland Hill picked up on this and modified the srs.db file and produced these 
messages:
http://www.mail-archive.com/qgis-user@lists.osgeo.org/msg07365.html and
http://www.mail-archive.com/qgis-user@lists.osgeo.org/msg07379.html

At least in the first message he used the following parameters to update the 
srs.db.
For AGD84:
'+ellps=aust_SA', '+a=6378160 +rf=298.2555 
+towgs84=-117.763,-51.51,139.061,-0.292,-0.443,-0.277,-0.191
and for AGD66:
'+ellps=aust_SA', '+a=6378160 +rf=298.2555 
+towgs84=-129.193,-41.212,130.73,-0.246,-0.374,-0.329,-2.955

Then placed the agd84 and agd66.gsb files into /usr/share/proj.

In attempting to replicate Roland I noticed that in 1.7 AGD84 is defined as 
ESPG:4203 and in Proj.4 format _projlonglat +elips=aust-SA 
+towards84=-117.763,-51.51,139.061,-0.443,-0.277,-0.191 +nadgrids=agd.gsb 
+wktext=no_defs which to me looked a lot like Roland's definition (please don't 
ask me to translate - I can't).

So I added in the two .gsb files and tested using this test set: 
http://members.westnet.com.au/ramon_et_al/SinglePoint2.zip

This contains a single point in my main area of operations, recorded originally 
in MGA94z51 using a Garmin GPS76. I've then converted the point into GDA94 
using GA's online converter, and into AGD66, AGD84 (and relevant UTM 
equivilants) using GA's agd-gda94.xls block calculations (see the included 
ASinglePoint.ods for calculations).

From the spreadsheet I've loaded the point in each different projection 
(ASinglePoint.qgs).

OTF works, and the points in that project cluster within a 1.1mx0.6m area. I 
would suspect that most of that error is probably more a reflection of the 
rough methods I used to convert them.


___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] Projection of Shapefile export

2011-03-21 Thread Ramon Andinach 

On 21/03/2011, at 10:21 , Brett Adams wrote:

> I thought this problem merely related to the failure of Australian 
> projections. Projection data doesn't exist, therefore can't be exported with 
> the shape file.
> 
> If the points are imported as WGS84 / UTMz51S, then export the shape file, 
> the projection file contains the relevant information. Its all good.
> 
> So if the projection issue is fixed, shape files will then export with the 
> required projection information.
> 
> 
> As far as I can see, there is no bug here. If anyone can show otherwise, we 
> need to record the bug and create a ticket (if it doesn't exist already)

The little bit of information I had before, I just tacked on to the end of your 
existing ticket.

I was going to put a bit of my spare time during this week into trying to 
replicate Roland's solution firstly on my home linux machine and then see if I 
could make that approach work on windows at work.

From my recollection last time I was trying Ken's solution and a few others 
like just defining a CRS and couldn't quite get there.

Let you know what eventuates.

-ramon.___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] Projection of Shapefile export

2010-09-22 Thread Roland Hill 




And if you change 20351 (AGD84
Z51) to 20251 (AGD66 Z51) in Micha's example conversion you get exactly
the same result, which shouldn't happen. Again, I believe this is
because none of the towgs84 parameters for AGD66 & AGD84 are
included in Proj4 or QGIS. The GDA94 definition is probably OK.

 
-- 
ROLAND HILL
Director
Four Winds
Technology Pty Ltd
Ph/Fax :  +61
(0)2 6366 9425
Mobile : 
+61 (0)41 880 7472
roland.h...@fourwindstechnology.com.au



Privacy and
Confidentiality Notice
The information contained
herein and any attachments are intended solely for the named
recipients. It may contain privileged confidential information. If you
are not an intended recipient, please delete the message and any
attachments then notify the sender of miss-delivery. Any use or
disclosure of the contents of either is unauthorised and may be
unlawful. All liability for viruses are excluded to the fullest extent
permitted by law.



On 22/09/10 22:50, Micha Silver wrote:
On
09/22/2010 11:59 AM, Brett Adams wrote:
  
   Gents,

Guess I'm wondering what the next step on this topic is. I don't know
what to do next.


Currently QGIS cannot separate ADG66/ AGD84 /GDA94 in either export or
import. If it can, then I haven't seen it and if its that fickle, its
way too dangerous for me.


Not keen on going to another package as they take time to learn. Is it
even an option to sponsor someone to sort this out for me?


Technically I'm way over my depth so can't offer many solutions myself.


I need some honest answers.


Thanks


  
Hi Brett:
  
  
Ken Norris indicated [1] that there might be a general problem
transforming from AGD66 or AGD84 (which use the aust_SA ellipsoid) to
GDA94 which is based on the GRS80 ellipsoid. He mentioned Oziexplorer
which does the transform correctly. To the best of my understanding,
proj4 (and thus QGIS) doesn't do any transform between the two.
  
Here's a trivial example. I took a coord pair from the file you posted
and ran it thru the proj4 cs2cs utility with epsg:20351 as the source
CRS (that's AGD84) and epsg:28351 as the target CRS (GDA94).
  
  
echo "275635 6411747" | cs2cs +init=epsg:20351 +to +init=epsg:28351
  
275635.81    6411759.44 0.00
  
  
Note that the resulting X-Y is offset by only about 12 meters. But
you're expecting a shift of nearly 200 m. So it seems proj4 is missing
here.
  
  
Normally datum transforms fit into either a 3 parameter or 7 parameter
model. It might be that the transform between those datums requires
something fancy, or localized that regular datum transforms cannot do.
So proj4 just doesn't handle it.
  
The best source for an authoritative  answer would, of course, be the
Australian govt' mapping agency.  If you learn anything more, please
let us know...
  
  
Regards,
  
Micha
  
  
  
[1]
http://osgeo-org.1803224.n2.nabble.com/Projection-of-Shapefile-export-tp5536736p5552105.html
  
  



___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] Projection of Shapefile export

2010-09-22 Thread Roland Hill 




 Brett,

Have a look at
http://blog.gmane.org/gmane.comp.gis.gdal.devel/month=20060301 for
another discussion on using AGD66 / AGD84. It basically describes how
the Proj4 definitions (which is where the Qgis CRS database comes from)
are not correct for AGD66 and AGD84. It also gives alternative
definitions a bit further down the discussion. You will notice how all
the AGD66 definitions have the same towgs84 parameters, as do the AGD84
definitions, but they are different for each datum. That's were the
shift from AGD66 to AGD84 should come from.

Note that my understanding is that the towgs84 parameters supplied are
not definitive, but are only one possible set. That is why you will
gradually build up a set of alternative, say AGD84 Z57, definitions as
you open files from various sources - MapInfo may define it with one
set of towgs84 parameters while ESRI may use a different set. You may
even be asked to supply your preferred parameters if you get data
delivered from an airborne contractor in anything other than WGS84 (I
am a geophysicist too).

I would test things by creating a custom CRS (Settings->Custom CRS)
using whichever of the definitions from the above link is appropriate
for your data. When you import/open your data (I can't quite recall
what exactly it is you where using) you should assign it to the Custom
CRS you created.

If your Custom CRS database does not appear to be working properly (ie
you can see your custom CRS, but it isn't sticking to your data
properly), then try to find your custom CRS database (qgis.db in your
home directory on Linux, I guess it is in Application Settings or
something similar in Windows) and delete it in case you have an old one
with an incorrect format.

The longer term fix is to install the correct, or at least 'a' correct,
version of the AGD66 and AGD84 definitions in the main crs database
(srs.db). While you could argue over which towgs84 parameters are
better, any would be preferable to the way it is now. The main and
custom CRS databases are SQLite3 databases. I might give this a go over
the next few days if I get a chance.

Also, this page may help to explain some of the Proj4 parameters
http://trac.osgeo.org/proj/wiki/GenParms

Hope this all makes sense - its getting late.

Regards,

Roland

 
-- 
ROLAND HILL
Director
Four Winds
Technology Pty Ltd
Ph/Fax :  +61
(0)2 6366 9425
Mobile : 
+61 (0)41 880 7472
roland.h...@fourwindstechnology.com.au



Privacy and
Confidentiality Notice
The information contained
herein and any attachments are intended solely for the named
recipients. It may contain privileged confidential information. If you
are not an intended recipient, please delete the message and any
attachments then notify the sender of miss-delivery. Any use or
disclosure of the contents of either is unauthorised and may be
unlawful. All liability for viruses are excluded to the fullest extent
permitted by law.



On 22/09/10 19:59, Brett Adams wrote:
 Gents,
  
Guess I'm wondering what the next step on this topic is. I don't know
what to do next.
  
  
Currently QGIS cannot separate ADG66/ AGD84 /GDA94 in either export or
import. If it can, then I haven't seen it and if its that fickle, its
way too dangerous for me.
  
  
Not keen on going to another package as they take time to learn. Is it
even an option to sponsor someone to sort this out for me?
  
  
Technically I'm way over my depth so can't offer many solutions myself.
  
  
I need some honest answers.
  
  
Thanks
  
  



___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] Projection of Shapefile export

2010-09-22 Thread Micha Silver 

On 09/22/2010 11:59 AM, Brett Adams wrote:

 Gents,
Guess I'm wondering what the next step on this topic is. I don't know 
what to do next.


Currently QGIS cannot separate ADG66/ AGD84 /GDA94 in either export or 
import. If it can, then I haven't seen it and if its that fickle, its 
way too dangerous for me.


Not keen on going to another package as they take time to learn. Is it 
even an option to sponsor someone to sort this out for me?


Technically I'm way over my depth so can't offer many solutions myself.

I need some honest answers.

Thanks


Hi Brett:

Ken Norris indicated [1] that there might be a general problem 
transforming from AGD66 or AGD84 (which use the aust_SA ellipsoid) to 
GDA94 which is based on the GRS80 ellipsoid. He mentioned Oziexplorer 
which does the transform correctly. To the best of my understanding, 
proj4 (and thus QGIS) doesn't do any transform between the two.
Here's a trivial example. I took a coord pair from the file you posted 
and ran it thru the proj4 cs2cs utility with epsg:20351 as the source 
CRS (that's AGD84) and epsg:28351 as the target CRS (GDA94).


echo "275635 6411747" | cs2cs +init=epsg:20351 +to +init=epsg:28351
275635.816411759.44 0.00

Note that the resulting X-Y is offset by only about 12 meters. But 
you're expecting a shift of nearly 200 m. So it seems proj4 is missing here.


Normally datum transforms fit into either a 3 parameter or 7 parameter 
model. It might be that the transform between those datums requires 
something fancy, or localized that regular datum transforms cannot do. 
So proj4 just doesn't handle it.
The best source for an authoritative  answer would, of course, be the 
Australian govt' mapping agency.  If you learn anything more, please let 
us know...


Regards,
Micha


[1] 
http://osgeo-org.1803224.n2.nabble.com/Projection-of-Shapefile-export-tp5536736p5552105.html


--
Micha Silver
Arava Development Co. +972-52-3665918
http://surfaces.co.il

___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] Projection of Shapefile export

2010-09-22 Thread Giovanni Manghi 
Hi,

>  Is it 
> even an option to sponsor someone

yes, if it is a bug or a missing feature (sorry, I didn't followed the
thread) you can fund QGIS and ask for a specific bug to be squashed or a
specific feature to be added

http://qgis.org/wiki/Bugs

cheers

-- Giovanni --


___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] Projection of Shapefile export

2010-09-22 Thread Brett Adams 

 Gents,
Guess I'm wondering what the next step on this topic is. I don't know 
what to do next.


Currently QGIS cannot separate ADG66/ AGD84 /GDA94 in either export or 
import. If it can, then I haven't seen it and if its that fickle, its 
way too dangerous for me.


Not keen on going to another package as they take time to learn. Is it 
even an option to sponsor someone to sort this out for me?


Technically I'm way over my depth so can't offer many solutions myself.

I need some honest answers.

Thanks

--
Brett Adams
Spinifex Geophysics
0438 861 974
SKYPE:brettadams_spinifex

___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] Projection of Shapefile export

2010-09-21 Thread Pierre Chevalier Géologue 
This rings a bell. I used to have similar hassles with Adindan 
projections vs. WGS84.
It works fine now, but I can't remember how and when I got this 
behaviour fixed.



Brett Adams claviota:
 Enabling "On-the-fly projection" does make a small difference but the 
offset should be approx 195-200m. I get 12m.



On 21/09/2010 4:42 AM, Jürgen E. Fischer wrote:

Hi Brett,

On Mon, 20. Sep 2010 at 22:11:17 +0800, Brett Adams wrote:
Actually, this is worse than I thought. QGIS is not recognising 
AGD84 or
GDA94. The data comes in but isn't attributed a CRS so defaults to 
WGS84.
If you import the text file as AGD84, then again as GDA94, there 
should be

approx 195m separation. QGIS plots them in exactly the same place.
But have have set a specific projection in your project properties 
and enabled

on-the-fly projection?  Otherwise QGIS won't reproject anything.

Jürgen







--

Pierre Chevalier Géologue EI
   Mesté Duran
   32100 Condom
 Tél+fax  :09 75 27 45 62
   05 62 28 06 83
06 37 80 33 64
 Émail  :   pierrechevaliergeolCHEZfree.fr
 icq#   :   10432285
 skype  :   pierre.chevalier1967
 http://pierremariechevalier.free.fr/pierre_chevalier_geologue




___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] Projection of Shapefile export

2010-09-20 Thread Brett Adams 

 Hi Roland,

Sorry but were straying well out of my depth here.

I vaguely understand the concept of what Proj4 and towgs84 do (but not 
really). I'll have to play with this idea.


Didn't know QGIS could do custom CRS definitions, although I'm not 
surprised.




On 21/09/2010 9:34 AM, Roland Hill wrote:
I had problems a while back using AGD66, AGD84 and GDA94. Like you, 
Brett, I wasn't getting the expected offset between data imported as 
AGD66/84 and GDA94. I am going from a distant memory so just ignore 
anything I say if it doesn't make sense, but there could be two issues 
here.


1) I think Proj4 only defines a CRS if it doesn't have conflicting 
definitions. There are multiple ways to define the towgs84 parameters 
for AGD66 and AGD84, so Proj4 doesn't include any of them. If you look 
at the Proj4 strings (in the QGIS CRS chooser) for say AGD66/Z55 and 
AGD84/Z55 they are the same ( +proj=utm +zone=55 +south +ellps=aust_SA 
+units=m +no_defs ). From memory, I think I got around this by 
choosing a towgs84 definition that I liked and then just used that as 
a custom CRS. You could use the custom CRS generated when importing a 
MapInfo file for example.


2) About 6 months ago a problem with the custom database structure was 
fixed. This was preventing custom CRS definitions from being used. If 
you have been using QGIS for a while it might be worth deleting the 
custom qgis.db (/.qgis/qgis.db in your home directory on Linux) so 
that Qgis recreates it. This could be the reason layers are defaulting 
to WGS84.


Roland


--
*ROLAND HILL*
*Director*
*Four Winds Technology Pty Ltd*
*Ph/Fax :  +61 (0)2 6366 9425*
*Mobile :  +61 (0)41 880 7472*
*roland.h...@fourwindstechnology.com.au 
*


*Privacy and Confidentiality Notice*
The information contained herein and any attachments are intended 
solely for the named recipients. It may contain privileged 
confidential information. If you are not an intended recipient, please 
delete the message and any attachments then notify the sender of 
miss-delivery. Any use or disclosure of the contents of either is 
unauthorised and may be unlawful. All liability for viruses are 
excluded to the fullest extent permitted by law.



On 21/09/10 10:21, Brett Adams wrote:
 Enabling "On-the-fly projection" does make a small difference but 
the offset should be approx 195-200m. I get 12m.



On 21/09/2010 4:42 AM, Jürgen E. Fischer wrote:

Hi Brett,

On Mon, 20. Sep 2010 at 22:11:17 +0800, Brett Adams wrote:
Actually, this is worse than I thought. QGIS is not recognising 
AGD84 or
GDA94. The data comes in but isn't attributed a CRS so defaults to 
WGS84.
If you import the text file as AGD84, then again as GDA94, there 
should be

approx 195m separation. QGIS plots them in exactly the same place.
But have have set a specific projection in your project properties 
and enabled

on-the-fly projection?  Otherwise QGIS won't reproject anything.

Jürgen






___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user



--
Brett Adams
Spinifex Geophysics
0438 861 974
SKYPE:brettadams_spinifex

___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] Projection of Shapefile export

2010-09-20 Thread Roland Hill 




I had problems a while back using AGD66,
AGD84 and GDA94. Like you, Brett, I wasn't
getting the expected offset between data imported as AGD66/84 and
GDA94. I am going from a distant memory so just ignore anything I say
if it doesn't make sense, but there could be two issues here.

1) I think Proj4 only defines a CRS if it doesn't have conflicting
definitions. There are multiple ways to define the towgs84 parameters
for AGD66 and AGD84, so Proj4 doesn't include any of them. If you look
at the Proj4 strings (in the QGIS CRS chooser) for say AGD66/Z55 and
AGD84/Z55 they are the same (
+proj=utm
+zone=55 +south +ellps=aust_SA +units=m +no_defs ). From memory, I
think I got around this by choosing a towgs84 definition that I liked
and then just used that as a custom CRS. You could use the custom CRS
generated when importing a MapInfo file for example.

2) About 6 months ago a problem with the custom database structure was
fixed. This was preventing custom CRS definitions from being used. If
you have been using QGIS for a while it might be worth deleting the
custom qgis.db (/.qgis/qgis.db in your home directory on Linux) so that
Qgis recreates it. This could be the reason layers are defaulting to
WGS84.

Roland

 
-- 
ROLAND HILL
Director
Four Winds
Technology Pty Ltd
Ph/Fax :  +61
(0)2 6366 9425
Mobile : 
+61 (0)41 880 7472
roland.h...@fourwindstechnology.com.au



Privacy and
Confidentiality Notice
The information contained
herein and any attachments are intended solely for the named
recipients. It may contain privileged confidential information. If you
are not an intended recipient, please delete the message and any
attachments then notify the sender of miss-delivery. Any use or
disclosure of the contents of either is unauthorised and may be
unlawful. All liability for viruses are excluded to the fullest extent
permitted by law.



On 21/09/10 10:21, Brett Adams wrote:
 Enabling
"On-the-fly projection" does make a small difference but the offset
should be approx 195-200m. I get 12m.
  
  
  
On 21/09/2010 4:42 AM, Jürgen E. Fischer wrote:
  
  Hi Brett,


On Mon, 20. Sep 2010 at 22:11:17 +0800, Brett Adams wrote:

Actually, this is worse than I thought.
QGIS is not recognising AGD84 or
  
GDA94. The data comes in but isn't attributed a CRS so defaults to
WGS84.
  
If you import the text file as AGD84, then again as GDA94, there should
be
  
approx 195m separation. QGIS plots them in exactly the same place.
  

But have have set a specific projection in your project properties and
enabled

on-the-fly projection?  Otherwise QGIS won't reproject anything.


Jürgen


  
  
  



___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] Projection of Shapefile export

2010-09-20 Thread Brett Adams 
 Enabling "On-the-fly projection" does make a small difference but the 
offset should be approx 195-200m. I get 12m.



On 21/09/2010 4:42 AM, Jürgen E. Fischer wrote:

Hi Brett,

On Mon, 20. Sep 2010 at 22:11:17 +0800, Brett Adams wrote:

Actually, this is worse than I thought. QGIS is not recognising AGD84 or
GDA94. The data comes in but isn't attributed a CRS so defaults to WGS84.
If you import the text file as AGD84, then again as GDA94, there should be
approx 195m separation. QGIS plots them in exactly the same place.

But have have set a specific projection in your project properties and enabled
on-the-fly projection?  Otherwise QGIS won't reproject anything.

Jürgen




--
Brett Adams
Spinifex Geophysics
0438 861 974
SKYPE:brettadams_spinifex

___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] Projection of Shapefile export

2010-09-20 Thread Jürgen E . Fischer 
Hi Brett,

On Mon, 20. Sep 2010 at 22:11:17 +0800, Brett Adams wrote:
> Actually, this is worse than I thought. QGIS is not recognising AGD84 or
> GDA94. The data comes in but isn't attributed a CRS so defaults to WGS84.
> If you import the text file as AGD84, then again as GDA94, there should be
> approx 195m separation. QGIS plots them in exactly the same place.

But have have set a specific projection in your project properties and enabled
on-the-fly projection?  Otherwise QGIS won't reproject anything.

Jürgen 

-- 
Jürgen E. Fischer norBIT GmbH   Tel. +49-4931-918175-20
Dipl.-Inf. (FH)   Rheinstraße 13Fax. +49-4931-918175-50
Software Engineer D-26506 Norden   http://www.norbit.de

-- 
norBIT Gesellschaft fuer Unternehmensberatung und Informationssysteme mbH
Rheinstrasse 13, 26506 Norden
GF: Jelto Buurman, HR: Amtsgericht Emden, HRB 5502

___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] Projection of Shapefile export

2010-09-20 Thread Brett Adams 


/Actually, this is worse than I thought. QGIS is not recognising AGD84 
or GDA94. The data comes in but isn't attributed a CRS so defaults to WGS84.
If you import the text file as AGD84, then again as GDA94, there should 
be approx 195m separation. QGIS plots them in exactly the same place. /



 Micha,
Sorry but I can't get your suggestions to work. I tried this a couple 
of times but it fails every time.


Settings->Options-CRS window is already set to "Prompt for CRS"
And when you next import the csv, does it indeed prompt for the CRS? 
/It does. But when I go to Settings>>Project properties the CRS is WGS84./




Tried changing the global to AGD84 / AMG zone 51, when I reopen the 
window its set to AGD66 so this setting is not holding.
Weird. There's no permissions problems? /Not that I'm aware of. It 
gives me no complaints

/




When I reimport the the CSV file the default is still WGS84, despite 
the global been set to AGD84 (or is it AGD66, I can't be sure). Again 
I select AGD84.


Layer > Properties > Specify CRS comes up as AGD66. Changed it to AGD84

When I export the layer, I'm again asked to specify the projection, 
AGD84


Exported *.prj file still has a datum of "D_unknown". Upon import of 
the shape file, my version of QGIS will now default this to AGD84 
however as the *.prj file is wrong it no good to anyone else.

Here's the prj that I'm getting:

PROJCS["UTM Zone 51, Southern Hemisphere",GEOGCS["Australian Natl & S. 
Amer. 1969",DATUM["D_unknown",SPHEROID["aust_SA",6
378160,298.25]],PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]],PROJECTION["Transverse_Mercator"],PARAMETER["l
atitude_of_origin",0],PARAMETER["central_meridian",123],PARAMETER["scale_factor",0.9996],PARAMETER["false_easting",50
],PARAMETER["false_northing",1000],UNIT["Meter",1]]


What should it be? Even tho' it's called "D_unknown" all the 
parameters seem to be there.

/
PRJ file should be something like below. /

PROJCS["AGD84 / AMG zone 
51",GEOGCS["AGD84",DATUM["D_Australian_1984",SPHEROID["Australian_National_Spheroid",6378160,298.25]],PRIMEM["Greenwich",0],

UNIT["Degree",0.017453292519943295]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",117],
PARAMETER["scale_factor",0.9996],PARAMETER["false_easting",50],PARAMETER["false_northing",1000],UNIT["Meter",1]]






The only real work around I have so far is to delete the *.prj file 
and then specify the projection on import.





On 19/09/2010 2:35 PM, Micha Silver wrote:

Micha Silver wrote:


Brett Adams wrote:


Micha,
I've attached csv data for importing and all shape files exported. 
I've added comments into your prior email also.



Hi Brett:

I do indeed see a problem exporting to the AGD84 CRS. I get a 
"Latitude or Longitude exceeds limits" error. This happens when the 
global CRS is set to WGS84 (which is the default in QGIS). Attached 
is a screen shot of the error.


Here's how to overcome the problem: Change your Global CRS to AGD84 
zone 51 (BTW, that's EPSG code 20351). This is done in the 
Settings->Options->CRS dialog window. Now re-import the csv file 
and then save as shapefile and you should be able to select the 
AGD84 CRS and have it saved correctly.
Alternatively, if you don't want to change the global CRS setting, 
you can select (in the same Settings->Options-CRS window) to 
"Prompt for CRS".  Then when you import the csv you'll be asked to 
specify the correct CRS for the points, You again select AGD84 and 
all should be well.


Cheers,
Micha

And just to complete the picture, Jurgen reminded me that in recent 
versions there's a setting for each layer to set what CRS it is. 
This is found in the Layer->properties window, General tab. Click on 
the button "Specify CRS" and choose the correct coordinate system. 
You can do this also for a delimited text file layer.  Then export 
to shapefile works as you'd expect..








--
Brett Adams
Spinifex Geophysics
0438 861 974
SKYPE:brettadams_spinifex

___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] Projection of Shapefile export

2010-09-20 Thread Micha Silver 

Brett Adams wrote:


 Micha,
Sorry but I can't get your suggestions to work. I tried this a couple 
of times but it fails every time.


Settings->Options-CRS window is already set to "Prompt for CRS"

And when you next import the csv, does it indeed prompt for the CRS?


Tried changing the global to AGD84 / AMG zone 51, when I reopen the 
window its set to AGD66 so this setting is not holding.

Weird. There's no permissions problems?


When I reimport the the CSV file the default is still WGS84, despite 
the global been set to AGD84 (or is it AGD66, I can't be sure). Again 
I select AGD84.


Layer > Properties > Specify CRS comes up as AGD66. Changed it to AGD84

When I export the layer, I'm again asked to specify the projection, AGD84

Exported *.prj file still has a datum of "D_unknown". Upon import of 
the shape file, my version of QGIS will now default this to AGD84 
however as the *.prj file is wrong it no good to anyone else.

Here's the prj that I'm getting:

PROJCS["UTM Zone 51, Southern Hemisphere",GEOGCS["Australian Natl & S. 
Amer. 1969",DATUM["D_unknown",SPHEROID["aust_SA",6
378160,298.25]],PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]],PROJECTION["Transverse_Mercator"],PARAMETER["l
atitude_of_origin",0],PARAMETER["central_meridian",123],PARAMETER["scale_factor",0.9996],PARAMETER["false_easting",50
],PARAMETER["false_northing",1000],UNIT["Meter",1]]


What should it be? Even tho' it's called "D_unknown" all the parameters 
seem to be there.






The only real work around I have so far is to delete the *.prj file 
and then specify the projection on import.





On 19/09/2010 2:35 PM, Micha Silver wrote:

Micha Silver wrote:


Brett Adams wrote:


Micha,
I've attached csv data for importing and all shape files exported. 
I've added comments into your prior email also.



Hi Brett:

I do indeed see a problem exporting to the AGD84 CRS. I get a 
"Latitude or Longitude exceeds limits" error. This happens when the 
global CRS is set to WGS84 (which is the default in QGIS). Attached 
is a screen shot of the error.


Here's how to overcome the problem: Change your Global CRS to AGD84 
zone 51 (BTW, that's EPSG code 20351). This is done in the 
Settings->Options->CRS dialog window. Now re-import the csv file and 
then save as shapefile and you should be able to select the AGD84 
CRS and have it saved correctly.
Alternatively, if you don't want to change the global CRS setting, 
you can select (in the same Settings->Options-CRS window) to "Prompt 
for CRS".  Then when you import the csv you'll be asked to specify 
the correct CRS for the points, You again select AGD84 and all 
should be well.


Cheers,
Micha

And just to complete the picture, Jurgen reminded me that in recent 
versions there's a setting for each layer to set what CRS it is. This 
is found in the Layer->properties window, General tab. Click on the 
button "Specify CRS" and choose the correct coordinate system. You 
can do this also for a delimited text file layer.  Then export to 
shapefile works as you'd expect..





___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] Projection of Shapefile export

2010-09-19 Thread Brett Adams 

 Micha,
Sorry but I can't get your suggestions to work. I tried this a couple of 
times but it fails every time.


Settings->Options-CRS window is already set to "Prompt for CRS"

Tried changing the global to AGD84 / AMG zone 51, when I reopen the 
window its set to AGD66 so this setting is not holding.


When I reimport the the CSV file the default is still WGS84, despite the 
global been set to AGD84 (or is it AGD66, I can't be sure). Again I 
select AGD84.


Layer > Properties > Specify CRS comes up as AGD66. Changed it to AGD84

When I export the layer, I'm again asked to specify the projection, AGD84

Exported *.prj file still has a datum of "D_unknown". Upon import of the 
shape file, my version of QGIS will now default this to AGD84 however as 
the *.prj file is wrong it no good to anyone else.




The only real work around I have so far is to delete the *.prj file and 
then specify the projection on import.





On 19/09/2010 2:35 PM, Micha Silver wrote:

Micha Silver wrote:


Brett Adams wrote:


Micha,
I've attached csv data for importing and all shape files exported. 
I've added comments into your prior email also.



Hi Brett:

I do indeed see a problem exporting to the AGD84 CRS. I get a 
"Latitude or Longitude exceeds limits" error. This happens when the 
global CRS is set to WGS84 (which is the default in QGIS). Attached 
is a screen shot of the error.


Here's how to overcome the problem: Change your Global CRS to AGD84 
zone 51 (BTW, that's EPSG code 20351). This is done in the 
Settings->Options->CRS dialog window. Now re-import the csv file and 
then save as shapefile and you should be able to select the AGD84 CRS 
and have it saved correctly.
Alternatively, if you don't want to change the global CRS setting, 
you can select (in the same Settings->Options-CRS window) to "Prompt 
for CRS".  Then when you import the csv you'll be asked to specify 
the correct CRS for the points, You again select AGD84 and all should 
be well.


Cheers,
Micha

And just to complete the picture, Jurgen reminded me that in recent 
versions there's a setting for each layer to set what CRS it is. This 
is found in the Layer->properties window, General tab. Click on the 
button "Specify CRS" and choose the correct coordinate system. You can 
do this also for a delimited text file layer.  Then export to 
shapefile works as you'd expect..



--
Brett Adams
Spinifex Geophysics
0438 861 974
SKYPE:brettadams_spinifex

___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] Projection of Shapefile export

2010-09-18 Thread Micha Silver 

Micha Silver wrote:


Brett Adams wrote:


Micha,
I've attached csv data for importing and all shape files exported. 
I've added comments into your prior email also.



Hi Brett:

I do indeed see a problem exporting to the AGD84 CRS. I get a 
"Latitude or Longitude exceeds limits" error. This happens when the 
global CRS is set to WGS84 (which is the default in QGIS). Attached is 
a screen shot of the error.


Here's how to overcome the problem: Change your Global CRS to AGD84 
zone 51 (BTW, that's EPSG code 20351). This is done in the 
Settings->Options->CRS dialog window. Now re-import the csv file and 
then save as shapefile and you should be able to select the AGD84 CRS 
and have it saved correctly.
Alternatively, if you don't want to change the global CRS setting, you 
can select (in the same Settings->Options-CRS window) to "Prompt for 
CRS".  Then when you import the csv you'll be asked to specify the 
correct CRS for the points, You again select AGD84 and all should be 
well.


Cheers,
Micha

And just to complete the picture, Jurgen reminded me that in recent 
versions there's a setting for each layer to set what CRS it is. This is 
found in the Layer->properties window, General tab. Click on the button 
"Specify CRS" and choose the correct coordinate system. You can do this 
also for a delimited text file layer.  Then export to shapefile works as 
you'd expect..

--
Micha
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] Projection of Shapefile export

2010-09-17 Thread Micha Silver 

Brett Adams wrote:


Micha,
I've attached csv data for importing and all shape files exported. 
I've added comments into your prior email also.



Hi Brett:

I do indeed see a problem exporting to the AGD84 CRS. I get a "Latitude 
or Longitude exceeds limits" error. This happens when the global CRS is 
set to WGS84 (which is the default in QGIS). Attached is a screen shot 
of the error.


Here's how to overcome the problem: Change your Global CRS to AGD84 zone 
51 (BTW, that's EPSG code 20351). This is done in the 
Settings->Options->CRS dialog window. Now re-import the csv file and 
then save as shapefile and you should be able to select the AGD84 CRS 
and have it saved correctly.
Alternatively, if you don't want to change the global CRS setting, you 
can select (in the same Settings->Options-CRS window) to "Prompt for 
CRS".  Then when you import the csv you'll be asked to specify the 
correct CRS for the points, You again select AGD84 and all should be well.


Cheers,
Micha


Brett

On 16/09/2010 9:29 PM, Micha Silver wrote:

Hi Brett:
(cc-ing the list also)

So if I understand correctly:
1- You have coordinate locations in an ascii file which are already 
in AGD84 zone 51 CRS : /Yes/
2- You have imported the locations into QGIS using the Delimited Text 
plugin :/Yes/
3- You then export to a shapefile. In the export dialog you select 
EPSG: 20351 for the CRS, but the resulting shapefile has NO 
projection information.(i.e. no *.prj is created?) /No./
/In the export dialogue I choose UTM - AGD84 / AMG zone 51(Don't know  
EPSG:20351 )/



4- You tried to set projection info with the "Vector->Data Management 
Tools->Define Current Projection" /Yes, this will run or say 
"Identical outout spatial reference system chosen"
/ 


5- Now the new shapefile still does not have CRS ? (i.e. QGIS asks 
what CRS when you load the shapefile) : /Not quite. Both prj and qpj 
files are created but the DATUM is always "D_unknown" which defaults 
to WGS84 when imported into QGIS / Mapinfo. /


If I've got the picture right, can you check that NO *.prj or *.qpj 
is created on the disk with the other elements of the shapefile?  
Would you mind sending me a bit of the ascii file with coord 
locations so I can retry?  If this persists, I suggest to report to 
Carson Farmer, creator of the fTools plugin. There was a problem in 
the past with creating a *.prj file, and I believe it has already 
been fixed. If not, he needs to know./ I think this problem is different.

/
Regards,
Micha

On 16/09/2010 14:47, Brett Adams wrote:

 Micha,
Updated to 1.5 and have tried defining the current projection as 
suggested but my exports never carry the correct projection. 
Everything goes in as AGD84/zone51 and comes out projection free 
which is then defaults to WGS84. Even tried loading files with the 
AGD84 projection from MAPINFO then exporting, and again they export 
without projection and default to WGS84.


I'm obviously missing a fairly major step, but is it really that hard?





On 16/09/2010 3:09 PM, Micha Silver wrote:

On 16/09/2010 09:01, Brett Adams wrote:

The original file is one created by me via import of delimted text so
I'm certain of the projection.
Question is does QGIS know the projection. If there is no *.prj 
attached to the shapefile, try to add projection info with the 
Vector->Data Management->Define Currnet Projection. Then do your 
export.



I'll try upgrading to 1.5. See if it makes any difference.


2010/9/16 Micha Silver:

On 16/09/2010 03:19, Brett Adams wrote:

Hi folks,
When I export a layer as a shapfile the projection is always WGS84,
even though I specify AGD84 on export.

Am I doing something wrong or is this a bug?


Are you sure the original shapefile has correct CRS data? That 
is, does it

have a *.prj ?
Also, I'm not 100% sure that exporting to a new CRS was available 
in 1.4.

You might go ahead and upgrade to 1.5 where the feature does work.

Running Version 1.4.0

Thanks




--
Micha Silver
http://www.surfaces.co.il/
Arava Development Co.  +972-52-3665918

















--
Brett Adams
Spinifex Geophysics
0438 861 974
SKYPE:brettadams_spinifex


This mail was received via Mail-SeCure System.


<>___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] Projection of Shapefile export

2010-09-16 Thread Micha Silver 

Hi Brett:
(cc-ing the list also)

So if I understand correctly:
1- You have coordinate locations in an ascii file which are already in 
AGD84 zone 51 CRS

2- You have imported the locations into QGIS using the Delimited Text plugin
3- You then export to a shapefile. In the export dialog you select EPSG: 
20351 for the CRS, but the resulting shapefile has NO projection 
information.(i.e. no *.prj is created?)
4- You tried to set projection info with the "Vector->Data Management 
Tools->Define Current Projection"
5- Now the new shapefile still does not have CRS ? (i.e. QGIS asks what 
CRS when you load the shapefile)


If I've got the picture right, can you check that NO *.prj or *.qpj is 
created on the disk with the other elements of the shapefile?  Would you 
mind sending me a bit of the ascii file with coord locations so I can 
retry?  If this persists, I suggest to report to Carson Farmer, creator 
of the fTools plugin. There was a problem in the past with creating a 
*.prj file, and I believe it has already been fixed. If not, he needs to 
know.


Regards,
Micha

On 16/09/2010 14:47, Brett Adams wrote:

 Micha,
Updated to 1.5 and have tried defining the current projection as 
suggested but my exports never carry the correct projection. 
Everything goes in as AGD84/zone51 and comes out projection free which 
is then defaults to WGS84. Even tried loading files with the AGD84 
projection from MAPINFO then exporting, and again they export without 
projection and default to WGS84.


I'm obviously missing a fairly major step, but is it really that hard?





On 16/09/2010 3:09 PM, Micha Silver wrote:

On 16/09/2010 09:01, Brett Adams wrote:

The original file is one created by me via import of delimted text so
I'm certain of the projection.
Question is does QGIS know the projection. If there is no *.prj 
attached to the shapefile, try to add projection info with the 
Vector->Data Management->Define Currnet Projection. Then do your export.



I'll try upgrading to 1.5. See if it makes any difference.


2010/9/16 Micha Silver:

On 16/09/2010 03:19, Brett Adams wrote:

Hi folks,
When I export a layer as a shapfile the projection is always WGS84,
even though I specify AGD84 on export.

Am I doing something wrong or is this a bug?


Are you sure the original shapefile has correct CRS data? That is, 
does it

have a *.prj ?
Also, I'm not 100% sure that exporting to a new CRS was available 
in 1.4.

You might go ahead and upgrade to 1.5 where the feature does work.

Running Version 1.4.0

Thanks




--
Micha Silver
http://www.surfaces.co.il/
Arava Development Co.  +972-52-3665918














--
Micha Silver
http://www.surfaces.co.il/
Arava Development Co.  +972-52-3665918

___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user

Re: [Qgis-user] Projection of Shapefile export

2010-09-15 Thread Micha Silver 

On 16/09/2010 03:19, Brett Adams wrote:

Hi folks,
When I export a layer as a shapfile the projection is always WGS84,
even though I specify AGD84 on export.

Am I doing something wrong or is this a bug?

   
Are you sure the original shapefile has correct CRS data? That is, does 
it have a *.prj ?
Also, I'm not 100% sure that exporting to a new CRS was available in 
1.4. You might go ahead and upgrade to 1.5 where the feature does work.

Running Version 1.4.0

Thanks

   



--
Micha Silver
http://www.surfaces.co.il/
Arava Development Co.  +972-52-3665918

___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user