[GRASS-user] 0 points loaded in v.surf.idw

[Markus Neteler]

On Tue, Jul 13, 2010 at 10:27 PM, Chang, Elizabeth
<Elizabeth.Chang at noblis.org> wrote:
> Hi,
>
> I’m trying to import some points from Postgres and convert it into raster
> using v.in.ogr and v.surf.idw
>
> v.in.ogr dsn="PG:host=tsf.noblis.org dbname=gis port=5433"
> layer=regcm.huss_monavg_2038 output=huss_monavg_2038_395_875
> where='timestamp=395.875'

Please post a snippet of the output or check it yourself of

 v.db.select huss_monavg_2038_395_875

just to be sure that you got the data actually in which you want to be in.

Furthermore,

.... map extent:
 |               N:   70:15:29.96521N    S:  19:27:00.387268N                 |
 |               E:  35:43:20.669862W    W: 157:42:49.097902W                 |

... computational (current) region extent;
north:      70:15:29.96521N
south:      69N
west:       120W
east:       110W
nsres:      0:15:05.993042
ewres:      0:20

The "south" value of the computational region does not quite cover
the map (intentionally?). Likeweise east and west. (also the pixels
are non-square).

Perhaps you want to first use (set resolution as appropriate):
 g.region vect=huss_monavg_2038_395_875 res=0:20 -a

to better fit the computational region to the map, then run v.surf.idw.

hope this helps,
Markus