[GRASS-dev] Re: [GRASS GIS] #1262: r.proj broken in all dev versions

[GRASS GIS]

#1262: r.proj broken in all dev versions
----------------------------------+-----------------------------------------
 Reporter:  cmbarton              |       Owner:  grass-dev@…              
     Type:  defect                |      Status:  new                      
 Priority:  blocker               |   Milestone:  6.4.1                    
Component:  Raster                |     Version:  unspecified              
 Keywords:  r.proj, reprojection  |    Platform:  MacOSX                   
      Cpu:  OSX/Intel             |  
----------------------------------+-----------------------------------------

Comment(by cmbarton):

 I just did an updated from the svn trunk (after a make distclean). This
 has changed the results and error, but has not corrected the problem. Just
 to be sure that there is nothing weird with the map, here is what I did. I
 opened the NC_spm_07 dataset and set the region to match the
 elev_state_500m dem. Then I tried to reproject into the state plane region
 the 1 arc second dem from the nc_ll location. I set debug=5 and here is
 the result.


 {{{
 GRASS 7.0.svn (nc_spm_07):~ > g.gisenv set="DEBUG=5"
 GRASS 7.0.svn (nc_spm_07):~ > r.proj input=elev_ned_1arcsec location=nc_ll
 mapset=PERMANENT output=elev_ned_1arcsec
 D1/5: G_find_raster(): name=elev_ned_1arcsec mapset=user1
 D2/5: G__read_Cell_head
 D2/5: G__read_Cell_head_array
 D3/5: region item: proj:       99
 D3/5: region item: zone:       0
 D3/5: region item: north:      318500
 D3/5: region item: south:      -16000
 D3/5: region item: east:       963000
 D3/5: region item: west:       124000
 D3/5: region item: cols:       1678
 D3/5: region item: rows:       669
 D3/5: region item: e-w resol:  500
 D3/5: region item: n-s resol:  500
 D3/5: region item: top:        1
 D3/5: region item: bottom:     0
 D3/5: region item: cols3:      1678
 D3/5: region item: rows3:      669
 D3/5: region item: depths:     1
 D3/5: region item: e-w resol3: 500
 D3/5: region item: n-s resol3: 500
 D3/5: region item: t-b resol:  1
 D1/5: G_find_raster(): name=elev_ned_1arcsec mapset=PERMANENT
 D1/5: G_find_raster2(): name=elev_ned_1arcsec mapset=PERMANENT
 D1/5: G_find_raster(): name=elev_ned_1arcsec mapset=PERMANENT
 D1/5: G_find_raster2(): name=elev_ned_1arcsec mapset=PERMANENT
 D1/5: G_find_raster2(): name=elev_ned_1arcsec mapset=PERMANENT
 D5/5: color buf = [0% yellow]
 D5/5: color buf = [20% green]
 D5/5: color buf = [40% cyan]
 D5/5: color buf = [60% blue]
 D5/5: color buf = [80% magenta]
 D5/5: color buf = [100% red]
 Input Projection Parameters:  +proj=longlat +a=6378137 +rf=298.257223563
 +no_defs
 Input Unit Factor: 1
 Output Projection Parameters:  +proj=lcc +lat_1=36.16666666666666
 +lat_2=34.33333333333334 +lat_0=33.75 +lon_0=-79 +x_0=609601.22 +y_0=0
 +no_defs +a=6378137 +rf=298.257222101 +towgs84=0.000,0.000,0.000
 Output Unit Factor: 1
 D2/5: G__read_Cell_head
 D2/5: G__read_Cell_head_array
 D3/5: region item: proj:       3
 D3/5: region item: zone:       0
 D3/5: region item: north:      35:54:41.000547N
 D3/5: region item: south:      35:35:17.001098N
 D3/5: region item: east:       78:27:08.000616W
 D3/5: region item: west:       78:49:17.999987W
 D3/5: region item: cols:       1330
 D3/5: region item: rows:       1164
 D3/5: region item: e-w resol:  0:00:01
 D3/5: region item: n-s resol:  0:00:01
 D3/5: region item: format:     -1
 D3/5: region item: compressed: 1
 D3/5: G_adjust_Cell_head: epsilon_ns: 8.59107e-07, epsilon_ew: 1e-06
 D3/5: G_adjust_Cell_head: epsilon_ns: 8.59107e-07, epsilon_ew: 1e-06
 D1/5: G_find_raster(): name=MASK mapset=user1
 D3/5: Top: xmin: -78.452221; xmax: -78.821668; ymin: 35.911390; ymax:
 35.588055
 D3/5: Right: xmin: -78.452221; xmax: -78.821668; ymin: 35.911390; ymax:
 35.588055
 D3/5: Bottom: xmin: -78.452221; xmax: -78.821668; ymin: 35.911390; ymax:
 35.588055
 D3/5: Left: xmin: -78.452221; xmax: -78.821668; ymin: 35.911390; ymax:
 35.588055
 D3/5: Extra check: xmin: -78.821528; xmax: -78.452361; ymin: 35.588195;
 ymax: 35.911250
 D3/5: Final check: xmin: -78.821667; xmax: -78.452222; ymin: 35.588056;
 ymax: 35.911389
 D3/5: G_adjust_Cell_head: epsilon_ns: 8.59107e-07, epsilon_ew: 1e-06
 D1/5: G_find_raster(): name=MASK mapset=user1
 D3/5: Top: xmin: 625709.948811; xmax: 659030.331786; ymin: 239768.184363;
 ymax: 239890.045778
 D3/5: Right: xmin: 625709.948811; xmax: 659228.054367; ymin:
 204049.351740; ymax: 239890.045778
 D3/5: Bottom: xmin: 625709.948811; xmax: 659228.054367; ymin:
 203927.002865; ymax: 239890.045778
 D3/5: Left: xmin: 625709.948811; xmax: 659228.054367; ymin: 203927.002865;
 ymax: 239890.045778
 D3/5: Extra check: xmin: 625709.948811; xmax: 659228.054367; ymin:
 203927.002865; ymax: 239890.045778
 D3/5: Final check: xmin: 625500.000000; xmax: 659500.000000; ymin:
 203500.000000; ymax: 240000.000000
 Segmentation fault
 }}}

 I can run a shell debugger on this again if needed.

 Michael

-- 
Ticket URL: <http://trac.osgeo.org/grass/ticket/1262#comment:14>
GRASS GIS <http://grass.osgeo.org>