[OSRS-PROJ] Re: [GRASS5] adding a new datum/projection?

Cliff Mugnier cjmce at LSU.edu
Wed Jul 23 10:19:58 EDT 2003


In regard to New Zealand, I have not researched their implementation in
detail, but from the little that I know; the ITRF is probably NOT the
way to go for "datum shifts" in New Zealand.  The original datum was not
recognized or utilized throughout both North Island and South Island,
and local property was referenced to local "meridional circuits" that
were perimeter surveys forming local coordinate systems and cast on
individual aphylactic Cassini-Soldner Grids.  Each county was an
"island" unto itself.  Topographic mapping at inch to the mile was on
the NZ Transverse Mercator before WWII, and the system remained that way
for decades.

They had a datum readjustment (I think) in the 70s and Lee came up with
the "NZ Map Grid" that covered both islands.

I do know that they have supposedly dropped Lee's projection for both
islands and may have gone back to Transverse Mercator.  However, NZ2000
Datum is actually changing monthly if not daily because of the enormous
speed the islands are moving towards Hawaii.  (I think maybe a
centimeter a week, not sure!)

Anyway, GIS coordinates need to be "islands" unto themselves (like the
old "meridional circuits"), and I don't think that geodetic "solutions"
are appropriate in this case.  I believe that the subject is so complex
for New Zealand, that the transformation from the ITRF_XX to a local
reference frame should be done by geodesists and maybe Land Surveyors,
but certainly not implemented into GIS software.

I may be incorrect in my assessment, but at present I believe that
"PROJ" should certainly NOT implement the current geodetic system (in an
international sense of the ITRF_XX) for New Zealand.  Leave it to the
geodesists for now.

Clifford J. Mugnier (cjmce at LSU.edu)
Chief of Geodesy and
Associate Director
Center for GeoInformatics
Department of Civil Engineering
LOUISIANA STATE UNIVERSITY
Baton Rouge, LA  70803
Voice:       (225) 578-8536
Facsimile:  (225) 578-8652
Pager:    1-(888) 365-5180
=======================================
http://www.ASPRS.org/resources.html
http://www.cee.lsu.edu/facultyStaff/mugnier/index.html
=======================================


-----Original Message-----
From: owner-osrs-proj at remotesensing.org
[mailto:owner-osrs-proj at remotesensing.org] On Behalf Of Paul Kelly
Sent: Wednesday, July 23, 2003 6:09 AM
To: H Bowman
Cc: grass5 at grass.itc.it; osrs-proj at remotesensing.org
Subject: [OSRS-PROJ] Re: [GRASS5] adding a new datum/projection?

 

 

 

 From:Clifford J Mugnier/cjmce/LSU at stjohnspoint.co.uk on 07/23/2003
12:09   
 PM CET

 




                                                                      
                                                                      
 Hello                                                                
                                                                      
 On Wed, 23 Jul 2003, H Bowman wrote:                                 
                                                                      
 > Hi again --                                                        
 >                                                                    
 >                                                                    
 > I've got some new data in yet another new projection, which uses   
 its own                                                              
 > new datum.                                                         
 >                                                                    
 >                                                                    
 > The datum is New Zealand Geodetic Datum 2000 (NZGD2000), which is  
 based                                                                
 > on the International Terrestrial Reference Frame 1996 (ITRF96).    
 This is                                                              
 > pretty much the same as WGS84 but uses the GRS80 ellipsoid         
 (slightly                                                            
 > different flattening parameter).                                   
 >                                                                    
 > So as far as I understand: NZGD2k (==ITRF96) is the WGS84 datum    
 combined                                                             
 > with the GRS80 ellipsoid.                                          
                                                                      
 We are pushing at the limits of the accuracy of the current          
 GRASS/PROJ                                                           
 system here. I looked at the document and to me it seems that        
 NZGD2000 is                                                          
 based on the ITRS datum. It is 'realised' (which I think means a     
 description of the current location of the earth's plates based on   
 measurements from a number of stations around the world) by several  
 ITRF                                                                 
 defintions, and each can be realised at any epoch (point in time).   
                                                                      
 For high accuracy you would define a shift that related each epoch   
 to WGS84                                                             
 (would change because of plate tectonics) and each would have a      
 separate                                                             
 datum entry in GRASS. It would get out of hand very quickly if these 
 were                                                                 
 all included but it may be possible to use GRASS in this way if you  
 could                                                                
 get figures for the shift (for New Zealand you need ITRF96 at 2000.0 
 epoch). http://www.iers.org/iers/products/itrf/ looks relevant but I 
 don't                                                                
 really know.                                                         
                                                                      
 > Except g.setproj won't let me do that, either with a new location  
 or                                                                   
 > retroactively. The old 5.0.2 version asks for ellipsoid first,     
 then                                                                 
 > datum, but complains when you try to mix wgs84 and grs80 and quits 
 > without writing anything. The latest g.setproj never gives you a   
 chance                                                               
 > to set the ellipsoid, it is chosen from your datum.                
                                                                      
 No it doesn't seem logical to do that. GRASS uses a very limited     
 definition                                                           
 of a datum which really just comes down to the ellipsoid used, so by 
 changing the ellipsoid you are really taking away everything the     
 datum                                                                
 stands for.                                                          
                                                                      
 > I'm not sure if I get the correct result if I just change the      
 ellps:                                                               
 > line to grs80 in the PROJ_INFO file?                               
                                                                      
 ellps: grs80                                                         
 towgs84: 0,0,0                                                       
                                                                      
 and make sure there is no 'datum' line.                              
                                                                      
 > Is anything else based on ITRF96? Would it be worth adding that as 
 a                                                                    
 > common denominator datum? I think the Map Grid of Australia might  
 be in                                                                
 > a similar situation (??).                                          
                                                                      
 It might be worth adding ITRS (International Terrestrial Reference   
 *System*), but to be compatible with other GIS we should probably    
 explicitly add New_Zealand_Geodetic_Datum_2000 and the others. But   
 as                                                                   
 they are the same as far as the level of accuracy GRASS/PROJ         
 currently                                                            
 provides goes, it seems like we could end up with a lot of clutter.  
 However I have already done that for several countries in Europe     
 which use                                                            
 the international ellipsoid but have different datum names, so there 
 is no                                                                
 point in discriminating....                                          
                                                                      
                                                                      
 >                                                                    
 > see:  (102k)                                                       
 >                                                                    
 http://www.linz.govt.nz/rcs/linz/17880/difference_wgs84_nzgd2000.pdf 
 >                                                                    
 >                                                                    
 >                                                                    
 > As for the new projection, it is New Zealand Transverse Mercator   
 (NZTM),                                                              
 > which is in terms of that new NZGD2k datum.                        
 > Details:                                                           
 > Datum: NZGD2k                                                      
 > Origin Latitude: 0° South                                          
 > Origin Longitude: 173° East                                        
 > False Northing: 10 000 000m N                                      
 > False Easting: 1 600 000m E                                        
 > Scale Factor: 0.9996                                               
 >                                                                    
 > Setting projection as tmerc and entering those terms goes          
 smoothly.                                                            
 > Should NZTM get its own projection entry or is entering by hand    
 with the                                                             
 > tmerc projection every time the way to go?                         
                                                                      
 Yes that is the way to do it until GRASS has co-ordinate system      
 support,                                                             
 which I don't see happening any time soon.. you have it very easy    
 with                                                                 
 New Zealand Map Grid which puts in all the parameters for you.       
                                                                      
 >                                                                    
 > for more details see:  (44k)                                       
 > http://www.linz.govt.nz/rcs/linz/5684/nztransverse_mercator.pdf    
 >                                                                    
 > That PDF also has a page of useful ellipsoid to grid formulae      
 which                                                                
 > might be useful to someone.                                        
 >                                                                    
 >                                                                    
 > It doesn't help much that half the country is moving 5cm/year in   
 the                                                                  
 > opposite direction to the other half..                             
                                                                      
 It would be useful to keep an eye out for any Free Software that     
 handles                                                              
 that level of accuracy                                               
 http://www-gpsg.mit.edu/~simon/gtgk/ has some interesting links      
                                                                      
 >                                                                    
 >                                                                    
 > thanks for any insight,                                            
                                                                      
 Well hopefully someone will correct the inevitable errors and        
 misunderstandings in what I've written above                         
                                                                      
 Paul                                                                 
 ----------------------------------------                             
 PROJ.4 Discussion List                                               
 See http://www.remotesensing.org/proj for subscription,              
 unsubscription                                                       
  and other information.                                              
                                                                      






More information about the grass-dev mailing list