[GRASS-user] r.terraflow: what's wrong?

2009-02-27 Thread FAROUX STEPHANIE

Hello,
I tried this command:

r.terraflow.short -s elev=world_dem4 filled=world_dem4_filled 
dir=world_dem4_mfdir swatershed=world_dem4_watershed 
accumulation=world_dem4_accu tci=world_dem4_tci memory=1000 
STREAM_DIR=/home/faroux/ --overwrite


(world_dem4 is the gtopo30 dem)

the first filled, dir and swatershed seems to well end; but when 
arriving at accumulation ant tci, there are a lot of messages:

warning tci has negative values
And when i draw it, there are many, many negative values everywhere.

Anybody can help me?

Thank you

Stéphanie
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] categories in raster map

2009-02-26 Thread FAROUX STEPHANIE

Hello,
I see the function r.cats which allows to read categories of a raster 
map, but how can i define these attributes? Should i pass by a vector 
map? In the database menu, it doesn't seem possible to add a table for a 
raster map.

Thank you
Stéphanie
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] categories in raster map

2009-02-26 Thread FAROUX STEPHANIE

Moritz Lennert wrote:

On 26/02/09 11:47, FAROUX STEPHANIE wrote:

Hello,
I see the function r.cats which allows to read categories of a raster 
map, but how can i define these attributes? Should i pass by a vector 
map? In the database menu, it doesn't seem possible to add a table 
for a raster map.


Raster maps have categories and labels. Categories have to be 
numerical, labels can be anything. Each pixel needs to have a category 
value (or NULL) and optionally you can have one (and only one) label 
per category value.


BTW, r.cats is now replaced by r.category.

If you want to link a cell with many attributes, you will have to use 
vector grids (v.mkgrid or r.to.vect as you suggestà).


Moritz


Thank you. And how do you put labels on categories please?
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] categories in raster map

2009-02-26 Thread FAROUX STEPHANIE

FAROUX STEPHANIE wrote:

Moritz Lennert wrote:

On 26/02/09 11:47, FAROUX STEPHANIE wrote:

Hello,
I see the function r.cats which allows to read categories of a 
raster map, but how can i define these attributes? Should i pass by 
a vector map? In the database menu, it doesn't seem possible to add 
a table for a raster map.


Raster maps have categories and labels. Categories have to be 
numerical, labels can be anything. Each pixel needs to have a 
category value (or NULL) and optionally you can have one (and only 
one) label per category value.


BTW, r.cats is now replaced by r.category.

If you want to link a cell with many attributes, you will have to use 
vector grids (v.mkgrid or r.to.vect as you suggestà).


Moritz


Thank you. And how do you put labels on categories please?
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


I use r.reclass that's it?
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] categories in raster map

2009-02-26 Thread FAROUX STEPHANIE

Moritz Lennert wrote:

On 26/02/09 13:34, FAROUX STEPHANIE wrote:

FAROUX STEPHANIE wrote:

Moritz Lennert wrote:

On 26/02/09 11:47, FAROUX STEPHANIE wrote:

Hello,
I see the function r.cats which allows to read categories of a 
raster map, but how can i define these attributes? Should i pass 
by a vector map? In the database menu, it doesn't seem possible to 
add a table for a raster map.


Raster maps have categories and labels. Categories have to be 
numerical, labels can be anything. Each pixel needs to have a 
category value (or NULL) and optionally you can have one (and only 
one) label per category value.


BTW, r.cats is now replaced by r.category.

If you want to link a cell with many attributes, you will have to 
use vector grids (v.mkgrid or r.to.vect as you suggestà).


Moritz


Thank you. And how do you put labels on categories please?
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


I use r.reclass that's it?


No, r.reclass changes category values, not labels.

I think it would be better if you explained a bit more clearly what 
you want to do, i.e. what do you have and what would you like to have.


Normally, r.category seems to be what you are looking for, using the 
rules= option.


I do see that the man page might not be particularly clear on its 
usage (the Description part only mentions printing, not managing 
labels and the part on dynamic labels is a bit cryptic). But the 
section Input from a file should give you what you need.


Moritz


Ok. I see it in the manual page of the 6.5 version.
Actually i work with the 6.2 version (i'm at office and i haven't 
control on the installed version) and i think r.cats doesn't have this 
option. But with r.reclass i achieved.
Here i have a map with integer categories values; i want to recode these 
categories but keep old values as labels.

Stéphanie
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] problem with export of raster data

2009-02-24 Thread FAROUX STEPHANIE

Markus Neteler wrote:

On Mon, Feb 23, 2009 at 12:15 PM, FAROUX STEPHANIE
stephanie.far...@meteo.fr wrote:
  

Hello
I try to export a raster integer map whose values range from - to
5880712. How can i specify that i want integer 4 bytes and not integer 2
bytes (it makes integer 2 bytes)??



If you export with r.out.gdal you can choose among

  type   File type
  options: Byte,Int16,UInt16,Int32,UInt32,Float32,Float64,CInt16,
   CInt32,CFloat32,CFloat64

This is the range which GRASS support for export to GDAL
supported formats (of course it depends if the format is compatible
with the selected type, see GDAL manual pages for details).

Markus

  

Yes i saw it.
However i wonder what is the difference between Uint32, Int32, CInt32? I 
suppose: unsigned integer 32 bits for the first, but the others?

Then, i want binary files, what should i choose among the proposed formats?
I tried with no format and with GTiff format. For example when i do: 
r.out.gdal input=world_fa2 format=GTiff2 type=Int32 output=world_fa_gdal
or r.out.gdal input=world_fa2 format=GTiff2 type=UInt32 
output=world_fa_gdal or r.out.gdal input=world_fa2 format=GTiff2 
type=CInt32 output=world_fa_gdal

it doesn't work and i get:

*** buffer overflow detected ***: g.parser terminated
=== Backtrace: =
/lib64/libc.so.6(__fortify_fail+0x37)[0x2b42e8de0157]
/lib64/libc.so.6[0x2b42e8dde7e0]
/lib64/libc.so.6[0x2b42e8b9]
/lib64/libc.so.6(_IO_default_xsputn+0x8e)[0x2b42e8d65b5e]
/lib64/libc.so.6(_IO_vfprintf+0x1633)[0x2b42e8d3b9b3]
/lib64/libc.so.6(__vsprintf_chk+0x9d)[0x2b42e8ddde5d]
/lib64/libc.so.6(__sprintf_chk+0x80)[0x2b42e8a0]
g.parser(main+0x2b9)[0x4012b9]
/lib64/libc.so.6(__libc_start_main+0xf4)[0x2b42e8d16074]
g.parser[0x400ea9]
=== Memory map: 
0040-00402000 r-xp  08:01 199895 
/usr/lib64/grass62/bin/g.parser
00602000-00603000 rw-p 2000 08:01 199895 
/usr/lib64/grass62/bin/g.parser
00603000-00624000 rw-p 00603000 00:00 0  
[heap]
2b42e81bc000-2b42e81d6000 r-xp  08:01 65538  
/lib64/ld-2.7.so

2b42e81d6000-2b42e81d8000 rw-p 2b42e81d6000 00:00 0
2b42e83d5000-2b42e83d6000 r--p 00019000 08:01 65538  
/lib64/ld-2.7.so
2b42e83d6000-2b42e83d7000 rw-p 0001a000 08:01 65538  
/lib64/ld-2.7.so
2b42e83d7000-2b42e8427000 r-xp  08:01 222338 
/usr/lib64/grass62/lib/libgrass_gis.6.2.3.so
2b42e8427000-2b42e8626000 ---p 0005 08:01 222338 
/usr/lib64/grass62/lib/libgrass_gis.6.2.3.so
2b42e8626000-2b42e8629000 rw-p 0004f000 08:01 222338 
/usr/lib64/grass62/lib/libgrass_gis.6.2.3.so

2b42e8629000-2b42e8641000 rw-p 2b42e8629000 00:00 0
2b42e8641000-2b42e8648000 r-xp  08:01 222312 
/usr/lib64/grass62/lib/libgrass_datetime.6.2.3.so
2b42e8648000-2b42e8848000 ---p 7000 08:01 222312 
/usr/lib64/grass62/lib/libgrass_datetime.6.2.3.so
2b42e8848000-2b42e8849000 rw-p 7000 08:01 222312 
/usr/lib64/grass62/lib/libgrass_datetime.6.2.3.so
2b42e8862000-2b42e8876000 r-xp  08:01 65654  
/lib64/libz.so.1.2.3
2b42e8876000-2b42e8a75000 ---p 00014000 08:01 65654  
/lib64/libz.so.1.2.3
2b42e8a75000-2b42e8a76000 rw-p 00013000 08:01 65654  
/lib64/libz.so.1.2.3

2b42e8a76000-2b42e8a77000 rw-p 2b42e8a76000 00:00 0
2b42e8a77000-2b42e8af6000 r-xp  08:01 65553  
/lib64/libm-2.7.so
2b42e8af6000-2b42e8cf6000 ---p 0007f000 08:01 65553  
/lib64/libm-2.7.so
2b42e8cf6000-2b42e8cf7000 r--p 0007f000 08:01 65553  
/lib64/libm-2.7.so
2b42e8cf7000-2b42e8cf8000 rw-p 0008 08:01 65553  
/lib64/libm-2.7.so
2b42e8cf8000-2b42e8e43000 r-xp  08:01 65545  
/lib64/libc-2.7.so
2b42e8e43000-2b42e9042000 ---p 0014b000 08:01 65545  
/lib64/libc-2.7.so
2b42e9042000-2b42e9046000 r--p 0014a000 08:01 65545  
/lib64/libc-2.7.so
2b42e9046000-2b42e9047000 rw-p 0014e000 08:01 65545  
/lib64/libc-2.7.so

2b42e9047000-2b42e904e000 rw-p 2b42e9047000 00:00 0
2b42e904e000-2b42e908d000 r--p  08:01 565776 
/usr/share/locale/UTF-8/LC_CTYPE
2b42e908d000-2b42e9094000 r--s  08:01 500251 
/usr/lib64/gconv/gconv-modules.cache
2b42e9094000-2b42e9095000 r--p  08:01 565931 
/usr/share/locale/fr_FR.UTF-8/LC_MESSAGES/SYS_LC_MESSAGES
2b42e9095000-2b42e90f5000 r--p  08:01 222408 
/usr/lib64/grass62/locale/fr/LC_MESSAGES/grassmods.mo
2b42e90f5000-2b42e90f6000 r-xp  08:01 500166 
/usr/lib64/gconv/ISO8859-1.so
2b42e90f6000-2b42e92f6000 ---p 1000 08:01 500166

[GRASS-user] problem with export of raster data

2009-02-23 Thread FAROUX STEPHANIE

Hello
I try to export a raster integer map whose values range from - to 
5880712. How can i specify that i want integer 4 bytes and not integer 2 
bytes (it makes integer 2 bytes)??

Thank you.
Stéphanie
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] v.digit

2009-02-20 Thread FAROUX STEPHANIE

Hamish wrote:

stephanie.faroux wrote:
  

i do submit and i have the error:

Cannot update table:
DBMI-DBF driver error:
SQL parser error in statement:
update australia set  where cat = 4
Error in db_execute_immediate()




update australia set  where cat = 4 is missing something between set
and where. set what?

Hamish

  
Yes i know it isn't right. I will see this problem next. For the moment 
grass fails to run... :-(


  



  


___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] v.digit

2009-02-20 Thread FAROUX STEPHANIE

Markus Neteler wrote:

On Thu, Feb 19, 2009 at 6:14 PM, FAROUX STEPHANIE
stephanie.far...@meteo.fr wrote:
  

Nikos Alexandris wrote:


And... what exactly does db.connect -p say?
  

driver:dbf
database:$GISDBASE/$LOCATION_NAME/$MAPSET/dbf/australia2.dbf
schema:(null)
group:toto



Please don't define a file name as database. It must be the *directory*:

db.connect -p
driver:dbf
database:$GISDBASE/$LOCATION_NAME/$MAPSET/dbf/
schema:
group:

Then it will work. So:
db.connect driver=dbf database='$GISDBASE/$LOCATION_NAME/$MAPSET/dbf/'
db.tables -p

Markus

  
Ok it works; now i have a little problem: i try to make a outline and i 
have difficulty to close it. My map is global at 1/120° resolution; 
maybe it's linked to the big distances?

I changed the threshold in the settings but it doesn't seem to be enough
Thank you


Cannot update table:

 DBMI-DBF driver error:
 SQL parser error in statement:
 update australia set  where cat = 4
 Error in db_execute_immediate()


when i want to save a new vector. But the vector is created. 



___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] v.digit

2009-02-20 Thread FAROUX STEPHANIE

Nikos Alexandris wrote:

Stephanie:
  
now i have a little problem: i try to make a outline and i 
have difficulty to close it. My map is global at 1/120° resolution; 
maybe it's linked to the big distances?

I changed the threshold in the settings but it doesn't seem to be enough



What kind of settings? g.region -p?

If you need a _frame_ (=vector) you can use v.in.region and for a grid
v.mkgrid.

Nikos


  
Nono i try to make a outline vector and i talk about the settings of the 
v.digit function


___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] to define a region which is not rectangular

2009-02-19 Thread FAROUX STEPHANIE

Hi,
In my raster map sea pixels and missing data have the same null value. 
I'd like to separate the two but i don't have a sea/land mask. Missing 
data are concentrated in some areas.

How can i do it?
I thought to define regions in which missing data would be gathered but 
if i only use the zoom to do it, it will take a long time.

Maybe using vectors?
Can anybody help me?
Thank you.
Stéphanie
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] to link two raster maps

2009-02-19 Thread FAROUX STEPHANIE

Hi,
Is it possible to link two raster maps in order to be able to get both 
pixels values by clicking on points with the mouse?

Thank you
Stéphanie
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] v.digit

2009-02-19 Thread FAROUX STEPHANIE

Hello,
I try v.digit; i want to define a new outline. When i want to close the 
outline, it opens a window form: attributes; new record was created. 
Assume data encoding as: [choice between utf-8,ascii,iso88859-1,koi8-t]; 
i do submit and i have the error:


Cannot update table:
DBMI-DBF driver error:
SQL parser error in statement:
update australia set  where cat = 4
Error in db_execute_immediate()


Could you help me?
Thank you
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] v.digit

2009-02-19 Thread FAROUX STEPHANIE

Nikos Alexandris wrote:

On Thu, 2009-02-19 at 17:28 +0100, FAROUX STEPHANIE wrote:
  

Hello,
I try v.digit; i want to define a new outline. When i want to close the 
outline, it opens a window form: attributes; new record was created. 
Assume data encoding as: [choice between utf-8,ascii,iso88859-1,koi8-t]; 
i do submit and i have the error:


Cannot update table:
DBMI-DBF driver error:
SQL parser error in statement:
update australia set  where cat = 4
Error in db_execute_immediate()


Could you help me?
Thank you



First things to try:

db.connect -p
v.db.connect -g YourVectorMap

Nikos


  

Yes i forgot to connect but it isn't enough to work...

I also try db.tables but i have the error:  


DBMI-DBF driver error:
Cannot open dbf database: $GISDBASE/$LOCATION_NAME/$MAPSET/dbf/australia.dbf

___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] v.digit

2009-02-19 Thread FAROUX STEPHANIE

Nikos Alexandris wrote:

[...]
  

Yes i forgot to connect but it isn't enough to work...



And... what exactly does db.connect -p say?


  


driver:dbf
database:$GISDBASE/$LOCATION_NAME/$MAPSET/dbf/australia2.dbf
schema:(null)
group:toto

I destroyed my vector layer and tried to make a new one...
Is there a place in the documentation where this action is described? 
For the moment i have no database in the dbf directory. I don't remember 
how i get one before (i had the database before i destroy my vector)

___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] v.digit

2009-02-19 Thread FAROUX STEPHANIE

Martin Landa wrote:

Hi,

2009/2/19 FAROUX STEPHANIE stephanie.far...@meteo.fr:

[...]

  

group:toto



What is group 'toto'?

Martin

  

nothing
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] error at launching GRASS

2009-02-19 Thread FAROUX STEPHANIE

aoutch

When i try to open grass, it doesn't work and i have the message:
D2/3: G__read_Cell_head
D2/3: G__read_Cell_head_array
D3/3: region item: proj:   3
D3/3: region item: zone:   0
D3/3: region item: north:  90N
D3/3: region item: south:  90S
D3/3: region item: east:   180E
D3/3: region item: west:   180W
D3/3: region item: cols:   43200
D3/3: region item: rows:   21600
D3/3: region item: e-w resol:  0:00:30
D3/3: region item: n-s resol:  0:00:30
D3/3: region item: top:1
D3/3: region item: bottom: 0
D3/3: region item: cols3:  1
D3/3: region item: rows3:  1
D3/3: region item: depths: 1
D3/3: region item: e-w resol3: 1
D3/3: region item: n-s resol3: 1
D3/3: region item: t-b resol:  1
D3/3: G_adjust_Cell_head: epsilon_ns: 4.62963e-08, epsilon_ew: 1e-06
D2/3: G__read_Cell_head
D2/3: G__read_Cell_head_array
D3/3: region item: proj:   3
D3/3: region item: zone:   0
D3/3: region item: north:  65:21:13.730342N
D3/3: region item: south:  56:13:13.743494N
D3/3: region item: east:   67:14:13.915718E
D3/3: region item: west:   54:08:13.934582E
D3/3: region item: cols:   1572
D3/3: region item: rows:   1096
D3/3: region item: e-w resol:  0:00:29.88
D3/3: region item: n-s resol:  0:00:29.88
D3/3: region item: top:1
D3/3: region item: bottom: 0
D3/3: region item: cols3:  71
D3/3: region item: rows3:  47
D3/3: region item: depths: 1
D3/3: region item: e-w resol3: 0:11:04.225086
D3/3: region item: n-s resol3: 0:11:39.574188
D3/3: region item: t-b resol:  1
D3/3: G_adjust_Cell_head: epsilon_ns: 9.12409e-07, epsilon_ew: 1e-06
D3/3: G_adjust_Cell_head: epsilon_ns: 9.12409e-07, epsilon_ew: 1e-06
Error in startup script: can't read monitor_zooms(1,1,n): no such variable
   while executing
lappend region $monitor_zooms($mon,1,$attr)
   (procedure MapCanvas::currentzoom line 13)
   invoked from within
MapCanvas::currentzoom $mon
   (procedure MapCanvas::coordconv line 23)
   invoked from within
MapCanvas::coordconv $mon
   (procedure MapCanvas::create line 68)
   invoked from within
MapCanvas::create
   (procedure Gm::startmon line 11)
   invoked from within
Gm::startmon
   (procedure Gm::create line 69)
   invoked from within
Gm::create
   (procedure main line 30)
   invoked from within
main $argc $argv
   (file /usr/lib64/grass62/etc/gm/gm.tcl line 521)

I was looking at the environment variables
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] v.generalize for area boundaries?

2009-02-17 Thread FAROUX STEPHANIE

Hi,
Is it possible with grass to replace missing values in a map by 
surrounding values by a sort of interpolation, knowing that sometimes 
many missing values are numerous one near another? If yes, how?

Thank you very much
Stéphanie
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] replacing missing data in a map

2009-02-17 Thread FAROUX STEPHANIE

Sorry i forgot to change the subject.

FAROUX STEPHANIE wrote:

Hi,
Is it possible with grass to replace missing values in a map by 
surrounding values by a sort of interpolation, knowing that sometimes 
many missing values are numerous one near another? If yes, how?

Thank you very much
Stéphanie
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user



___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] problem with reprojection of Hydro1k data

2009-02-16 Thread FAROUX STEPHANIE

Thank you for your answer Markus.
I didn't read the note in the r.proj page. I will try again and then we 
will see.

Stéphanie

Markus Neteler wrote:

Stéphanie,

On Fri, Feb 13, 2009 at 2:18 PM, FAROUX STEPHANIE
stephanie.far...@meteo.fr wrote:
  

Hello,
First i'd like to know if anybody already used the Hydro1k data from USGS?
I try to make global maps of parameters from continental ones.
The initial projection is lambert azimutal equal area, with different
latitude and longitude of origins for different continents. I defined one
location by continent; then i reprojected all maps in a latlon global
location, with WGS84 datum.
There are gaps between continents, for example between Asia and Europe and
between Europe and Africa.
I wonder if the problem comes:
- from the data which don't perfectly cover the earth



... could be possible, I didn't try with Hydro1k data from USGS.

  

- from my reprojection in grass



How did you do that? Say, did you use the v.in.region trick
mentioned in the r.proj manual page which helps a lot?

Markus

  


___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] problem with reprojection of Hydro1k data

2009-02-16 Thread FAROUX STEPHANIE
I did it; results are not exactly the same but it's not enough because 
gaps between continents are still more important.


FAROUX STEPHANIE wrote:

Thank you for your answer Markus.
I didn't read the note in the r.proj page. I will try again and then 
we will see.

Stéphanie

Markus Neteler wrote:

Stéphanie,

On Fri, Feb 13, 2009 at 2:18 PM, FAROUX STEPHANIE
stephanie.far...@meteo.fr wrote:
 

Hello,
First i'd like to know if anybody already used the Hydro1k data from 
USGS?

I try to make global maps of parameters from continental ones.
The initial projection is lambert azimutal equal area, with different
latitude and longitude of origins for different continents. I 
defined one

location by continent; then i reprojected all maps in a latlon global
location, with WGS84 datum.
There are gaps between continents, for example between Asia and 
Europe and

between Europe and Africa.
I wonder if the problem comes:
- from the data which don't perfectly cover the earth



... could be possible, I didn't try with Hydro1k data from USGS.

 

- from my reprojection in grass



How did you do that? Say, did you use the v.in.region trick
mentioned in the r.proj manual page which helps a lot?

Markus

  


___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user



___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] problem with reprojection of Hydro1k data

2009-02-16 Thread FAROUX STEPHANIE
I try the v.in.region + v.proj + g.region +  r.proj process (explained 
in the r.proj page) with a region whose reprojection doesn't allow to 
get a closed boarder: the projection of the raster doesn't work. Is it 
normal?

Thank you

FAROUX STEPHANIE wrote:
I did it; results are not exactly the same but it's not enough because 
gaps between continents are still more important.


FAROUX STEPHANIE wrote:

Thank you for your answer Markus.
I didn't read the note in the r.proj page. I will try again and then 
we will see.

Stéphanie

Markus Neteler wrote:

Stéphanie,

On Fri, Feb 13, 2009 at 2:18 PM, FAROUX STEPHANIE
stephanie.far...@meteo.fr wrote:
 

Hello,
First i'd like to know if anybody already used the Hydro1k data 
from USGS?

I try to make global maps of parameters from continental ones.
The initial projection is lambert azimutal equal area, with different
latitude and longitude of origins for different continents. I 
defined one

location by continent; then i reprojected all maps in a latlon global
location, with WGS84 datum.
There are gaps between continents, for example between Asia and 
Europe and

between Europe and Africa.
I wonder if the problem comes:
- from the data which don't perfectly cover the earth



... could be possible, I didn't try with Hydro1k data from USGS.

 

- from my reprojection in grass



How did you do that? Say, did you use the v.in.region trick
mentioned in the r.proj manual page which helps a lot?

Markus

  


___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user



___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user



___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] problem with reprojection of Hydro1k data

2009-02-13 Thread FAROUX STEPHANIE

Hello,
First i'd like to know if anybody already used the Hydro1k data from USGS?
I try to make global maps of parameters from continental ones.
The initial projection is lambert azimutal equal area, with different 
latitude and longitude of origins for different continents. I defined 
one location by continent; then i reprojected all maps in a latlon 
global location, with WGS84 datum.
There are gaps between continents, for example between Asia and Europe 
and between Europe and Africa.

I wonder if the problem comes:
- from the data which don't perfectly cover the earth
- from my reprojection in grass
Thank you if somebody can help me
Regards
Stéphanie

___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user