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

GRASS GIS trac at osgeo.org
Fri Jan 28 20:20:29 EST 2011


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

Comment(by hellik):

 Replying to [comment:26 cmbarton]:
 >
 > It still is very important to fix for 7.0. Since it is broken on the
 Mac, I'm betting it's still broken for Windows in 7.0 too, but this needs
 testing.
 >

 tested with WinGRASS-7.0.SVN-r45224-1-Setup.exe (nightly build of today).

 again with source: +proj=tmerc target: +proj=longlat

 it's very interesting, something is crashing during reprojecting, but the
 raster itself is reprojected.


 {{{
 r.proj input=DemEuropeFPp100nn location=ak_noe3 mapset=data
 output=test3debug
 D1/3: G_find_raster(): name=test3debug mapset=user1
 D2/3: G__read_Cell_head
 D2/3: G__read_Cell_head_array
 D3/3: region item: proj:       3
 D3/3: region item: zone:       0
 D3/3: region item: north:      60N
 D3/3: region item: south:      30N
 D3/3: region item: east:       30E
 D3/3: region item: west:       0
 D3/3: region item: cols:       3600
 D3/3: region item: rows:       3600
 D3/3: region item: e-w resol:  0:00:30
 D3/3: region item: n-s resol:  0:00:30
 D3/3: region item: top:        1
 D3/3: region item: bottom:     0
 D3/3: region item: cols3:      3600
 D3/3: region item: rows3:      3600
 D3/3: region item: depths:     1
 D3/3: region item: e-w resol3: 0:00:30
 D3/3: region item: n-s resol3: 0:00:30
 D3/3: region item: t-b resol:  1
 D3/3: G_adjust_Cell_head: epsilon_ns: 2.77778e-007,
 epsilon_ew: 1e-006
 D1/3: G_find_raster(): name=DemEuropeFPp100nn mapset=data
 D1/3: G_find_raster2(): name=DemEuropeFPp100nn mapset=data
 D1/3: G_find_raster(): name=DemEuropeFPp100nn mapset=data
 D3/3: adding rule 79=79.00 0 191 191  498=498.20 0 255 0
 D3/3: adding rule 498=498.20 0 255 0  917=917.40 255 255 0
 D3/3: adding rule 917=917.40 255 255 0  1336=1336.60 255 127
 0
 D3/3: adding rule 1336=1336.60 255 127 0  1755=1755.80 191
 127 63
 D3/3: adding rule 1755=1755.80 191 127 63  2175=2175.00 20
 20 20
 D2/3: G__read_Cell_head
 D2/3: G__read_Cell_head_array
 D3/3: region item: proj:       99
 D3/3: region item: zone:       0
 D3/3: region item: north:      432254.46875
 D3/3: region item: south:      253707.0625
 D3/3: region item: east:       804844.1875
 D3/3: region item: west:       610026.4375
 D3/3: region item: cols:       1948
 D3/3: region item: rows:       1785
 D3/3: region item: e-w resol:  100.00911191
 D3/3: region item: n-s resol:  100.02655812
 D3/3: region item: format:     -1
 D3/3: region item: compressed: 1
 D1/3: G_find_raster(): name=MASK mapset=user1
 D3/3: Top: xmin: 804844.187501; xmax: 610026.437499; ymin:
 432254.468751; ymax: 253707.062499
 D3/3: Right: xmin: 804844.187501; xmax: 610026.437499; ymin:
 432254.468751; ymax: 253707.062499
 D3/3: Bottom: xmin: 804844.187501; xmax: 610026.437499;
 ymin: 432254.468751; ymax: 253707.062499
 D3/3: Left: xmin: 804844.187501; xmax: 610026.437499; ymin:
 432254.468751; ymax: 253707.062499
 D3/3: Extra check: xmin: 610076.442056; xmax: 804794.182944;
 ymin: 253757.075779; ymax: 432204.455471
 D3/3: Final check: xmin: 610026.437500; xmax: 804844.187500;
 ymin: 253707.062500; ymax: 432254.468750
 D1/3: G_find_raster(): name=MASK mapset=user1
 D3/3: Top: xmin: 14.419322; xmax: 17.081289; ymin:
 49.011643; ymax: 49.027566
 D3/3: Right: xmin: 14.419322; xmax: 17.081289; ymin:
 47.420479; ymax: 49.027566
 D3/3: Bottom: xmin: 14.419322; xmax: 17.081289; ymin:
 47.407697; ymax: 49.027566
 D3/3: Left: xmin: 14.419322; xmax: 17.081289; ymin:
 47.407697; ymax: 49.027566
 D3/3: Extra check: xmin: 14.419322; xmax: 17.081289; ymin:
 47.407697; ymax: 49.027566
 D3/3: Final check: xmin: 14.416667; xmax: 17.083333; ymin:
 47.400000; ymax: 49.033333
 D3/3: G_adjust_Cell_head: epsilon_ns: 2.77778e-007,
 epsilon_ew: 1e-006
 D3/3: G_adjust_Cell_head: epsilon_ns: 5.10204e-006,
 epsilon_ew: 1e-006
 D3/3: G_adjust_Cell_head: epsilon_ns: 5.10204e-006,
 epsilon_ew: 1e-006

 Eingabe:
 Spalten: 1948 (1948)
 Zeilen: 1785 (1785)
 Norden: 432254.468750 (432254.468750)
 Süden: 253707.062500 (253707.062500)
 Westen: 610026.437500 (610026.437500)
 Osten: 804844.187500 (804844.187500)
 EW-res: 100.009112
 NS-res: 100.026558

 Ausgabe:
 Spalten: 320 (3600)
 Zeilen: 196 (3600)
 Norden: 49.033333 (60.000000)
 Süden: 47.400000 (30.000000)
 Westen: 14.416667 (0.000000)
 Osten: 17.083333 (30.000000)
 EW-res: 0.008333
 NS-res: 0.008333

 D1/3: G_find_raster(): name=MASK mapset=user1
 D1/3: G_find_raster2(): name=DemEuropeFPp100nn mapset=data
 D2/3: G__read_Cell_head
 D2/3: G__read_Cell_head_array
 D3/3: region item: proj:       99
 D3/3: region item: zone:       0
 D3/3: region item: north:      432254.46875
 D3/3: region item: south:      253707.0625
 D3/3: region item: east:       804844.1875
 D3/3: region item: west:       610026.4375
 D3/3: region item: cols:       1948
 D3/3: region item: rows:       1785
 D3/3: region item: e-w resol:  100.00911191
 D3/3: region item: n-s resol:  100.02655812
 D3/3: region item: format:     -1
 D3/3: region item: compressed: 1
 D1/3: G_find_raster2(): name=DemEuropeFPp100nn mapset=data
 D1/3: G_find_raster2(): name=DemEuropeFPp100nn mapset=data
 D1/3: G_find_raster2(): name=DemEuropeFPp100nn mapset=data
 D1/3: G_find_raster2(): name=DemEuropeFPp100nn mapset=data
 D1/3: G_find_raster2(): name=DemEuropeFPp100nn mapset=data
 D3/3: create window mapping (1948 columns)
 D1/3: G_find_raster2(): name=DemEuropeFPp100nn mapset=data
 D1/3: G_find_raster2(): name=DemEuropeFPp100nn mapset=data
 D1/3: G_find_raster(): name=MASK mapset=user1
 Allokiere Speicher und lese die Eingabekarte...
 D3/3: G_adjust_Cell_head: epsilon_ns: 5.10204e-006,
 epsilon_ew: 1e-006
 D3/3: G_adjust_Cell_head: epsilon_ns: 5.10204e-006,
 epsilon_ew: 1e-006
 Projeziere...
 D1/3: close test3debug compressed
 D1/3: G_find_raster2(): name=test3debug mapset=user1
 D1/3: G_find_raster2(): name=test3debug mapset=user1
 D1/3: G_find_raster2(): name=test3debug mapset=user1
 D3/3: G_recreate_command()
 (Sat Jan 29 02:14:39 2011) Command finished (10 sec)
 }}}

 somewhere around here
 {{{
 D1/3: close test3debug compressed
 D1/3: G_find_raster2(): name=test3debug mapset=user1
 D1/3: G_find_raster2(): name=test3debug mapset=user1
 D1/3: G_find_raster2(): name=test3debug mapset=user1
 D3/3: G_recreate_command()
 }}}
 there seems to be a crash (a windows-warning-messages pops up)

 Helmut

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



More information about the grass-dev mailing list