[GRASS-dev] Change the resolution in the output map

Juan Diego Giraldo Osorio Juan.Giraldo at upct.es
Fri Nov 30 11:22:18 EST 2007


Hi coworkers

I'm working with grass62. I'm triying to create a raster map with finer
resolution from another with coarse resolution. The problem is the next:

- My coarse resolution map (input map) have 500 m of resolution, with
1000x1000 rows*cols. This information is the same as the current region.
- I want to create a finer resolution map (output map) with 100 m of
resolution, with 5000x5000 rows*cols.
- I have been trying with the next code:

	/* Adjust the region to new resolution */
	region.rows = nsubpixel*nrows;
	region.cols = nsubpixel*ncols;
	region.format = 0;
	region.compressed = 1;
	G_adjust_Cell_head( &region , 1 , 1 );
	G_set_window(&region);

  nsubpixel is 5 (500/100), and I'm writing a new header file with this
new information (region.rows=5000, region.cols=5000, region.ew_res=100,
region.ns_res=100) but when I read the new map, the next message error
is shows:

  WARNING: error reading null row 115
  ... a lot of times.

Thanks for your help.

Juan Diego Giraldo
--
********************************
Juan Diego Giraldo Osorio
Tel. 968 327 047
Departamento de Ingeniería Térmica y de Fluidos
Universidad Politécnica de Cartagena
España
********************************





More information about the grass-dev mailing list