Opened 13 years ago

Last modified 6 years ago

#649 new defect

Nonsensical results from r.resamp.interp with certain region settings

Reported by: sorokine Owned by: grass-dev@…
Priority: normal Milestone: 6.4.6
Component: Raster Version: 6.3.0
Keywords: r.resamp.interp Cc:
CPU: OSX/Intel Platform: MacOSX

Description

I am resampling a raster to higher resolution using bilinear or bicubic methods but I am getting erroneous result.

Steps to reproduce and more information:

r.info rast1
...
 |   Type of Map:  raster               Number of Categories: 255             |
 |   Data Type:    FCELL                                                      |
 |   Rows:         128                                                        |
 |   Columns:      256                                                        |
 |   Total Cells:  32768                                                      |
 |        Projection: Latitude-Longitude                                      |
 |            N: 89:37:40.633467N    S:        90S   Res: 1:24:12.036199      |
 |            E: 179:17:48.75E    W: 179:17:48.75E   Res: 1:24:22.5           |
 |   Range of data:    min = 0.000000  max = 0.000000                         |
...

r.region rast=rast1
g.region nsres=0.2 ewres=0.2
r.resamp.interp --o input=rast1 output=rast1resampled

rast1resampled does not look like anything close to expected.

Workaround

Shove 3 degrees from the region before resampling:

g.region n=87 s=-87 w=-177 e=177

Resampled raster will be correct.

Change History (2)

comment:1 Changed 13 years ago by hamish

Keywords: r.resamp.interp added
Milestone: 6.4.0

comment:2 Changed 6 years ago by neteler

Milestone: 6.4.06.4.6
Note: See TracTickets for help on using tickets.