[GRASS-dev] r.proj produces null map

Hi,

on G7, I made several tentatives trying to reproject raster maps using r.proj, it doesn’t yield error but the map is completely empty (null). I made several tentatives on several maps, then also exported and reprojected with gdalwarp, to test if the problem was related to the data. Can someone reproduce this, before i open a ticket.

Thank you

···

Best regards,

Dr. Margherita DI LEO
Scientific / technical project officer

European Commission - DG JRC
Institute for Environment and Sustainability (IES)
Via Fermi, 2749
I-21027 Ispra (VA) - Italy - TP 261

Tel. +39 0332 78 3600
margherita.di-leo@jrc.ec.europa.eu

Disclaimer: The views expressed are purely those of the writer and may not in any circumstance be regarded as stating an official position of the European Commission.

Hi,

2015-05-18 17:16 GMT+02:00 Margherita Di Leo <diregola@gmail.com>:

on G7, I made several tentatives trying to reproject raster maps using
r.proj, it doesn't yield error but the map is completely empty (null). I
made several tentatives on several maps, then also exported and reprojected
with gdalwarp, to test if the problem was related to the data. Can someone
reproduce this, before i open a ticket.

are you sure that you have set in destination location region
correctly? See r.proj's -g flag.

Martin

--
Martin Landa
http://geo.fsv.cvut.cz/gwiki/Landa
http://gismentors.cz/mentors/landa

[forgot the list ]

···

On Mon, May 18, 2015 at 8:35 AM, Margherita Di Leo <diregola@gmail.com> wrote:

On Mon, May 18, 2015 at 8:18 AM, Martin Landa <landa.martin@gmail.com> wrote:

are you sure that you have set in destination location region
correctly? See r.proj’s -g flag.

I think so, unless I can’t see the obvious… here’s an example

r.proj -g input=accum location=GMTED2010-wgs84 mapset=res075
n=3036289.28490875 s=1327842.48956189 w=3859987.52509802 e=5240836.36134497 rows=7200 cols=7200

the region in the destination mapset is
g.region -p
projection: 99 (Lambert Azimuthal Equal Area)
zone: 0
datum: etrs89
ellipsoid: grs80
north: 3036500
south: 1327500
west: 3859750
east: 5241000
nsres: 250
ewres: 250
rows: 6836
cols: 5525
cells: 37768900

r.proj input=accum location=GMTED2010-wgs84 mapset=res075 output=acc_prova method=nearest

Input:
Cols: 7202 (7200)
Rows: 7200 (7200)
North: 49.999861 (49.999861)
South: 34.999861 (34.999861)
West: 364.995694 (4.999861)
East: 379.999861 (19.999861)
EW-res: 0.002083
NS-res: 0.002083

Output:
Cols: 5524 (5525)
Rows: 6834 (6836)
North: 3036250.000000 (3036500.000000)
South: 1327750.000000 (1327500.000000)
West: 3859750.000000 (3859750.000000)
East: 5240750.000000 (5241000.000000)
EW-res: 250.000000
NS-res: 250.000000

Allocating memory and reading input raster map…
Projecting…
r.proj complete.

r.info acc_prova
±---------------------------------------------------------------------------+
| Map: acc_prova Date: Mon May 18 08:26:37 2015 |
| Mapset: wetness-idx Login of Creator: user |
| Location: GMTED2010-LAEA |
| DataBase: /home/user/GIS/grassdata |
| Title: ( acc_prova ) |

Timestamp: none
Type of Map: raster Number of Categories: 0
Data Type: DCELL
Rows: 6834
Columns: 5524
Total Cells: 37751016
Projection: Lambert Azimuthal Equal Area
N: 3036250 S: 1327750 Res: 250
E: 5240750 W: 3859750 Res: 250
Range of data: min = -nan max = -nan
Data Description:
generated by r.proj
Comments:
r.proj input=“accum” location=“GMTED2010-wgs84” mapset=“res075” outp\
ut=“acc_prova” method=“nearest” memory=300
±---------------------------------------------------------------------------+

Best regards,

Dr. Margherita DI LEO
Scientific / technical project officer

European Commission - DG JRC
Institute for Environment and Sustainability (IES)
Via Fermi, 2749
I-21027 Ispra (VA) - Italy - TP 261

Tel. +39 0332 78 3600
margherita.di-leo@jrc.ec.europa.eu

Disclaimer: The views expressed are purely those of the writer and may not in any circumstance be regarded as stating an official position of the European Commission.

Best regards,

Dr. Margherita DI LEO
Scientific / technical project officer

European Commission - DG JRC
Institute for Environment and Sustainability (IES)
Via Fermi, 2749
I-21027 Ispra (VA) - Italy - TP 261

Tel. +39 0332 78 3600
margherita.di-leo@jrc.ec.europa.eu

Disclaimer: The views expressed are purely those of the writer and may not in any circumstance be regarded as stating an official position of the European Commission.