[GRASS-dev] Re: [GRASS GIS] #1262: r.proj broken in all dev versions
GRASS GIS
trac at osgeo.org
Sat Apr 9 16:24:46 EDT 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:29 hamish]:
> Replying to [comment:28 cmbarton]:
> > On the Mac, the raster is not reprojected.
>
> Is that only in 7 now, or still in 6.x too?
>
> Is it now only broken in 7 on WinGrass?
> (aka is it now ok for 6.x on all platforms?)
tested in windows with
source srs
{{{
+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
+to_meter=1
}}}
destination srs
{{{
+proj=longlat +no_defs +a=6378137 +rf=298.257223563
+towgs84=0.000,0.000,0.000
}}}
'''WinGRASS-6.4.SVN-r45864-1-Setup.exe => r.proj finishes without any
problem'''
{{{
r.proj --verbose input=elev_state_500m location=nc_spm_08 mapset=PERMANENT
output=elev_state_500m_wingrass641 method=bilinear
Eingabeparameter der Projektion: +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
Einheitenfaktor der Eingabe: 1
Ausgabeparameter der Projektion: +proj=longlat +no_defs
+a=6378137 +rf=298.257223563 +towgs84=0.000,0.000,0.000
Einheitenfaktor der Ausgabe: 1
Eingabe:
Spalten: 1580 (1678)
Zeilen: 564 (669)
Norden: 266000.000000 (318500.000000)
Süden: -16000.000000 (-16000.000000)
Westen: 173000.000000 (124000.000000)
Osten: 963000.000000 (963000.000000)
EW-res: 500.000000
NS-res: 500.000000
Ausgabe:
Spalten: 9 (9)
Zeilen: 3 (3)
Norden: 36.557555 (36.557555)
Süden: 33.490551 (33.490551)
Westen: -84.225943 (-84.225943)
Osten: -75.051985 (-75.051985)
EW-res: 1.019329
NS-res: 1.022335
Allokiere Speicher und lese die Eingabekarte...
Projeziere...
r.proj komplett.
(Sat Apr 09 21:56:33 2011) Command finished (1 sec)
}}}
'''WinGRASS-6.5.SVN-r45872-1-Setup.exe => r.proj finishes without any
problem'''
{{{
r.proj --verbose input=elev_state_500m location=nc_spm_08 mapset=PERMANENT
output=elev_state_500m_wingrass65 method=bilinear_f
Eingabeparameter der Projektion: +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
Einheitenfaktor der Eingabe: 1
Ausgabeparameter der Projektion: +proj=longlat +no_defs
+a=6378137 +rf=298.257223563 +towgs84=0.000,0.000,0.000
Einheitenfaktor der Ausgabe: 1
Eingabe:
Spalten: 1580 (1678)
Zeilen: 564 (669)
Norden: 266000.000000 (318500.000000)
Süden: -16000.000000 (-16000.000000)
Westen: 173000.000000 (124000.000000)
Osten: 963000.000000 (963000.000000)
EW-res: 500.000000
NS-res: 500.000000
Ausgabe:
Spalten: 9 (9)
Zeilen: 3 (3)
Norden: 36.557555 (36.557555)
Süden: 33.490551 (33.490551)
Westen: -84.225943 (-84.225943)
Osten: -75.051985 (-75.051985)
EW-res: 1.019329
NS-res: 1.022335
Allokiere Speicher und lese die Eingabekarte...
Projeziere...
r.proj komplett.
(Sat Apr 09 22:02:14 2011) Command finished (1 sec)
}}}
'''WinGRASS-7.0.SVN-r45874-1-Setup.exe => there seems still to be
some crash nearby G_recreate_command(), but the raster is reprojected
(see https://trac.osgeo.org/grass/ticket/1262?replyto=29#comment:27)'''
{{{
GRASS 7.0.svn> g.gisenv set=DEBUG=5
GRASS 7.0.svn> r.proj --verbose input=elev_state_500m location=nc_spm_08
mapset=PERMANENT output=elevstate500mwingrass70debug method=bilinear_f
D1/5: G_find_raster(): name=elevstate500mwingrass70debug
mapset=testreproject
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: 36:33:27.198541N
D3/5: region item: south: 33:29:25.982362N
D3/5: region item: east: 75:03:07.144363W
D3/5: region item: west: 84:13:33.393728W
D3/5: region item: cols: 9
D3/5: region item: rows: 3
D3/5: region item: e-w resol: 1:01:09.583263
D3/5: region item: n-s resol: 1:01:20.405393
D3/5: region item: top: 1
D3/5: region item: bottom: 0
D3/5: region item: cols3: 9
D3/5: region item: rows3: 3
D3/5: region item: depths: 1
D3/5: region item: e-w resol3: 1:01:09.583263
D3/5: region item: n-s resol3: 1:01:20.405393
D3/5: region item: t-b resol: 1
D3/5: G_adjust_Cell_head: epsilon_ns: 0.000333333, epsilon_ew: 1e-006
D1/5: G_find_raster(): name=elev_state_500m mapset=PERMANENT
D1/5: G_find_raster2(): name=elev_state_500m mapset=PERMANENT
D1/5: G_find_raster(): name=elev_state_500m mapset=PERMANENT
D3/5: adding rule -10=-10.00 0 0 255 0=2000.00 0 191 191
D3/5: adding rule 0=-10.00 0 191 191 10=2000.00 0 255 255
D3/5: adding rule 10=-10.00 0 255 255 30=2000.00 0 255 0
D3/5: adding rule 30=-10.00 0 255 0 120=2000.00 255 255 0
D3/5: adding rule 120=-10.00 255 255 0 350=2000.00 255 127 0
D3/5: adding rule 350=-10.00 255 127 0 700=2000.00 191 127 63
D3/5: adding rule 700=-10.00 191 127 63 2000=2000.00 255 255 255
Eingabeparameter der Projektion: +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
Einheitenfaktor der Eingabe: 1
Ausgabeparameter der Projektion: +proj=longlat +no_defs +a=6378137
+rf=298.257223563 +towgs84=0.000,0.000,0.000
Einheitenfaktor der Ausgabe: 1
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: format: -1
D3/5: region item: compressed: 1
D1/5: G_find_raster(): name=MASK mapset=testreproject
D3/5: Top: xmin: 184786.522595; xmax: 919362.320937; ymin: 254813.342689;
ymax: 264840.952469
D3/5: Right: xmin: 184786.522595; xmax: 927217.047945; ymin: 33426.108872;
ymax: 264840.952469
D3/5: Bottom: xmin: 174014.338912; xmax: 927217.047945; ymin:
27992.367084; ymax: 264840.952469
D3/5: Left: xmin: 174014.338912; xmax: 927217.047945; ymin: 27992.367084;
ymax: 264840.952469
D3/5: Extra check: xmin: 174014.338912; xmax: 962750.000000; ymin:
-15750.000000; ymax: 264840.952469
D3/5: Final check: xmin: 174000.000000; xmax: 963000.000000; ymin:
-16000.000000; ymax: 265000.000000
D1/5: G_find_raster(): name=MASK mapset=testreproject
D3/5: Top: xmin: -83.844648; xmax: -75.078044; ymin: 36.049566; ymax:
36.145545
D3/5: Right: xmin: -83.844648; xmax: -75.078044; ymin: 33.547070; ymax:
36.145545
D3/5: Bottom: xmin: -83.844648; xmax: -75.078044; ymin: 33.514977; ymax:
36.145545
D3/5: Left: xmin: -83.844648; xmax: -75.078044; ymin: 33.514977; ymax:
36.145545
D3/5: Extra check: xmin: -83.716278; xmax: -75.561649; ymin: 34.001718;
ymax: 36.046388
D3/5: Final check: xmin: -84.225943; xmax: -75.051985; ymin: 33.490551;
ymax: 36.557555
D3/5: G_adjust_Cell_head: epsilon_ns: 0.000333333, epsilon_ew: 1e-006
D3/5: G_adjust_Cell_head: epsilon_ns: 0.000333333, epsilon_ew: 1e-006
D3/5: G_adjust_Cell_head: epsilon_ns: 0.000333333, epsilon_ew: 1e-006
Eingabe:
Spalten: 1580 (1678)
Zeilen: 564 (669)
Norden: 266000.000000 (318500.000000)
S³den: -16000.000000 (-16000.000000)
Westen: 173000.000000 (124000.000000)
Osten: 963000.000000 (963000.000000)
EW-res: 500.000000
NS-res: 500.000000
Ausgabe:
Spalten: 9 (9)
Zeilen: 3 (3)
Norden: 36.557555 (36.557555)
S³den: 33.490551 (33.490551)
Westen: -84.225943 (-84.225943)
Osten: -75.051985 (-75.051985)
EW-res: 1.019329
NS-res: 1.022335
D1/5: G_find_raster(): name=MASK mapset=testreproject
D1/5: G_find_raster2(): name=elev_state_500m mapset=PERMANENT
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: format: -1
D3/5: region item: compressed: 1
D1/5: G_find_raster2(): name=elev_state_500m mapset=PERMANENT
D1/5: G_find_raster2(): name=elev_state_500m mapset=PERMANENT
D1/5: G_find_raster2(): name=elev_state_500m mapset=PERMANENT
D1/5: G_find_raster2(): name=elev_state_500m mapset=PERMANENT
D1/5: G_find_raster2(): name=elev_state_500m mapset=PERMANENT
D3/5: create window mapping (1580 columns)
D1/5: G_find_raster2(): name=elev_state_500m mapset=PERMANENT
D1/5: G_find_raster2(): name=elev_state_500m mapset=PERMANENT
D1/5: G_find_raster(): name=MASK mapset=testreproject
Allokiere Speicher und lese die Eingabekarte...
100%
D3/5: G_adjust_Cell_head: epsilon_ns: 0.000333333, epsilon_ew: 1e-006
D3/5: G_adjust_Cell_head: epsilon_ns: 0.000333333, epsilon_ew: 1e-006
Projeziere...
100%
D1/5: close elevstate500mwingrass70debug compressed
D1/5: G_find_raster2(): name=elevstate500mwingrass70debug
mapset=testreproject
D1/5: G_find_raster2(): name=elevstate500mwingrass70debug
mapset=testreproject
D1/5: G_find_raster2(): name=elevstate500mwingrass70debug
mapset=testreproject
D3/5: G_recreate_command()
GRASS 7.0.svn>
}}}
Helmut
--
Ticket URL: <https://trac.osgeo.org/grass/ticket/1262#comment:31>
GRASS GIS <http://grass.osgeo.org>
More information about the grass-dev
mailing list