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

2010-09-22 Thread Brett Adams

 Ken,
That's a legendary effort if your getting sensible results with the icsm 
files.  Do you think the format is suitable that it could be used in a 
sponsored bug fix?


I gave this a shot myself but couldn't get it run at all, but then I 
have no idea what I'm doing so thats no real indication.


I've had no success in getting parameters to stick (AGD66, AGD84 or 
GDA94) on shape or mapinfo files.


Following solution was proposed by Roland Hill (22/09/2010  8:59pm)

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.



Brett

On 23/09/2010 12:47 PM, Ken Norris wrote:

As Micha suggested, the National AGD66 distortion grid at
http://www.icsm.gov.au/icsm/gda/gdatm/national66.zip
and AGD84 distortion grid at
http://www.icsm.gov.au/icsm/gda/gdatm/national84.zip
seem to be the way forward.

When unzipped, put into the proj directory (for me usr/share/proj) and
coupled via the +nadgrids= parameter in a Custom CRS, it seems to yield the
correct transform when tested.

+proj=utm +zone=55 +south +ellps=aust_SA +nadgrids=AGD66.gsb +units=m
+no_defs

(I renamed the unzipped distortion file 'A66 National (13.09.01).gsb' file
to 'AGD66.gsb' for convenience)

But I can't make the parameters stick in the CRS of a vector file using the
override facility. Any ideas?





--
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


[Qgis-user] Re: Projection of Shapefile export

2010-09-22 Thread Ken Norris

As Micha suggested, the National AGD66 distortion grid at
http://www.icsm.gov.au/icsm/gda/gdatm/national66.zip
and AGD84 distortion grid at
http://www.icsm.gov.au/icsm/gda/gdatm/national84.zip
seem to be the way forward.

When unzipped, put into the proj directory (for me usr/share/proj) and
coupled via the +nadgrids= parameter in a Custom CRS, it seems to yield the
correct transform when tested.

+proj=utm +zone=55 +south +ellps=aust_SA +nadgrids=AGD66.gsb +units=m
+no_defs

(I renamed the unzipped distortion file 'A66 National (13.09.01).gsb' file
to 'AGD66.gsb' for convenience)

But I can't make the parameters stick in the CRS of a vector file using the
override facility. Any ideas?


-- 
View this message in context: 
http://osgeo-org.1803224.n2.nabble.com/Projection-of-Shapefile-export-tp5536736p5561744.html
Sent from the qgis-user mailing list archive at Nabble.com.
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user


[Qgis-user] Re: Reproject WMS

2010-09-22 Thread Tim Michelsen
>> * Why not use the vector in WGS84?
>>   my result map needs to be published projected
>> * I first load the vector in arc and then the WMS tile.
> 
> Then do all your work in WGS84 and just before publication, reproject
> your result maps with the Raster -> Projection tool of QGIS (which is a
> wrapper to the gdalwarp command).
I think not the user should adapt to the software but vice-versa!

Nearly all maps on the world prepared for interchange with authorities
have some sort of projected SRS while most WMS deliver ggr. / WGS84.

So you don't wanna always spend time for reprojecting but rather do the
actual work.

Anyway, thanks for the answer.

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


[Qgis-user] Re: Reproject WMS

2010-09-22 Thread Tim Michelsen
> if you mean that arc* does raster reprojection I believe we all know
> that is true as it is also true that at the moment QGIS doesn't
> reproject on-the fly rasters. I believe also that this feature is on the
> QGIS roadmap.
Any idea who this support will be or where I could find more information?

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


[Qgis-user] standolone shapefile editor

2010-09-22 Thread Mohammed Rashad

Hi all,How can i make a standalone shapefile editor using qgis code. which 
part/file is responsible for editing shapefiles and viewing them.>Is it 
possible to make a shapfile editor from scratch. because i need an application 
to create/edit/view shapefiles  ___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user


Re: [Qgis-user] Re: Reproject WMS

2010-09-22 Thread Jean-Claude Repetto

Le 22/09/2010 15:02, Timmie a écrit :

Hello,
thanks for all your answers.

* Why not use the vector in WGS84?
  my result map needs to be published projected
* I first load the vector in arc and then the WMS tile.


Then do all your work in WGS84 and just before publication, reproject 
your result maps with the Raster -> Projection tool of QGIS (which is a 
wrapper to the gdalwarp command).

___
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


[Qgis-user] Re: Reproject WMS

2010-09-22 Thread Timmie
> I am taking the maps from here:
http://msrmaps.com/webservices.aspx 
> and need to reproject to 
> NAD83 / Texas Centric Albers Equal Area (EPSG:3083)



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


[Qgis-user] Re: Reproject WMS

2010-09-22 Thread Timmie
Hello,
thanks for all your answers.

* Why not use the vector in WGS84?
 my result map needs to be published projected
* I first load the vector in arc and then the WMS tile.


Tp be more clear:

I am taking the maps from here:

and need to reproject to 
NAD83 / Texas Centric Albers Equal Area (EPSG:3083)

Best regards,
Timmie

___
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] Re: Reproject WMS

2010-09-22 Thread Nicklas Avén
Hallo You are sure it is not the overlay that is reprojected on the fly? I 
guess that just depends on if you load the wms before the overlay layer or vice 
verse. Otherwise I learned on FOSS4G that mapproxy can do things like that. If 
you do the request through mapproxy you can do transformations like that in 
mapproxy (mapproxy.org).
 But I suspect that image quality must be quite bad but I have never tried  
/Nicklas 
2010-09-22 Timmie wrote:

> > Is it possible for QGIS to add such layer and repoject them to state plane 
> > on
>> > the fly?
>> 
>> afaik it is not possible to reproject WMS layers, you have to stick with
>> the SRSs the server makes available.
>Sorry, but in the software of some well-known vendor this works.
>Overlay matches perfectly with the same service.
>
>
>
>___
>Qgis-user mailing list
>Qgis-user@lists.osgeo.org
>http://lists.osgeo.org/mailman/listinfo/qgis-user
>
> ___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user


Re: [Qgis-user] Re: Reproject WMS

2010-09-22 Thread Jean-Claude Repetto

Le 22/09/2010 14:10, Timmie a écrit :

Is it possible for QGIS to add such layer and repoject them to state plane on
the fly?


afaik it is not possible to reproject WMS layers, you have to stick with
the SRSs the server makes available.

Sorry, but in the software of some well-known vendor this works.
Overlay matches perfectly with the same service.



Why don't you reproject the overlay to WGS84 ? If overlay is vector 
data, it is much easier and quicker to reproject.

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


Re: [Qgis-user] Re: Reproject WMS

2010-09-22 Thread Paolo Cavallini
Il 22/09/2010 14:10, Timmie ha scritto:
>>> Is it possible for QGIS to add such layer and repoject them to state plane 
>>> on
>>> the fly?
>>
>> afaik it is not possible to reproject WMS layers, you have to stick with
>> the SRSs the server makes available.
> Sorry, but in the software of some well-known vendor this works.
> Overlay matches perfectly with the same service.

however, it makes no sense tyo reproject WMS cleint-side when it can be done
server-side (of course, not all servers are configured to provide theyr wms in 
the
projection you need).
all the best.
-- 
Paolo Cavallini: http://www.faunalia.it/pc
___
Qgis-user mailing list
Qgis-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/qgis-user


Re: [Qgis-user] Re: Reproject WMS

2010-09-22 Thread Giovanni Manghi

> Sorry, but in the software of some well-known vendor this works.
> Overlay matches perfectly with the same service.

if you mean that arc* does raster reprojection I believe we all know
that is true as it is also true that at the moment QGIS doesn't
reproject on-the fly rasters. I believe also that this feature is on the
QGIS roadmap.

cheers

-- Giovanni --

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


[Qgis-user] Re: Reproject WMS

2010-09-22 Thread Timmie
> > Is it possible for QGIS to add such layer and repoject them to state plane 
> > on
> > the fly?
> 
> afaik it is not possible to reproject WMS layers, you have to stick with
> the SRSs the server makes available.
Sorry, but in the software of some well-known vendor this works.
Overlay matches perfectly with the same service.



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


Re: [Qgis-user] Reproject WMS

2010-09-22 Thread Giovanni Manghi
Hi,

> Is it possible for QGIS to add such layer and repoject them to state plane on
> the fly?

afaik it is not possible to reproject WMS layers, you have to stick with
the SRSs the server makes available.

cheers

-- Giovanni --


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


[Qgis-user] Reproject WMS

2010-09-22 Thread Timmie
Hello,
I have a WMS from MSNR (http://msrmaps.com/webservices.aspx) which is served in
WGS84.

Is it possible for QGIS to add such layer and repoject them to state plane on
the fly?

What are you ideas to handle this?

Best regards,
Timmie

___
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