[GRASS-user] Memory problems with v.in.ogr

2010-11-25 Thread Pierre Roudier
Dear list,

I am trying to load a pretty big shapefile on GRASS (6Mb +, 211621
features)  using v.in.ogr, but each attempt fails at the 'break
boundaries'' stage (at 98%!):

GRASS 7.0.svn (NZTM2000):~ > v.in.ogr
dsn=~/Documents/DATA/LCDB2/ni_nzmg.shp out=ni_lcdb2 --o
Projection of input dataset and current location appear to match
Layer: ni_nzmg
Counting polygons for 211621 features...
Importing map 211621 features...
 100%
-
Building topology for vector map ...
Registering primitives...
334611 primitives registered
18459167 vertices registered
Number of nodes: 211681
Number of primitives: 334611
Number of points: 0
Number of lines: 0
Number of boundaries: 334611
Number of centroids: 0
Number of areas: -
Number of isles: -
-
WARNING: Cleaning polygons, result is not guaranteed!
-
Break polygons:
 100%
 100%
-
Remove duplicates:
 100%
-
Break boundaries:
ERROR: G_calloc: unable to allocate 50 * 8 bytes of memory at
   allocation.c:82

The shp loads without any problems on QGIS.

Here is my config:

GRASS 7.0.svn (NZTM2000):~ > g.region -p
projection: 99 (New Zealand Map Grid)
zone:   0
datum:  nzgd49
ellipsoid:  international
north:  6115136.31665518
south:  6114206.6264583
west:   2703766.53250003
east:   2705045.64941879
nsres:  1.99933376
ewres:  1.99862019
rows:   465
cols:   640
cells:  297600

GRASS 7.0.svn (NZTM2000):~ > g.version
GRASS 7.0.svn44411 (2010)

OS: Opensuse 11.3, with 7Gb RAM.

Any pointers appreciated,

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


Re: [GRASS-user] Grid label format in psmap

2010-11-25 Thread Hamish
John wrote:
> I have a UTM dataset that I am plotting with psmap.  I
> would like to have grid squares at 1000 m intervals.  I
> can do this easily with the following:
> 
> grid 1000
>    color black
>    numbers 1
>    font helvetica
>    fontsize 12
> end*
> 
> *However, on the resulting plot, the gridline labels are
> truncated into kilometers e.g. 498, 7114, instead of 498000,
> 7114000.  Is there a parameter that I can modify
> somewhere so that I can control output format and display
> the full number.  (Similar to D_FORMAT /
> OUTPUT_DEGREE_FORMAT in GMT, I suppose).

right, it automatically truncates off extra zeros. I usually
use the map_info instruction which shows the full bounds.


Unfortunately there's no way to switch off that feature without
modifying the source code. if you wish to do that, edit
ps/ps.map/do_grid.c and in the do_grid_numbers() function
edit this line -from, +to:

-sprintf(num_text, "%s", format_northing(g, rounded_grid));
+sprintf(num_text, "%s", format_northing(g, 1));

and then recompile. (untested)


you may also try out ps.output from the Addons collection, it
has more advanced border controls:
  http://grass.osgeo.org/screenshots/images/probabilidad_de_infeccion_crop.jpg

(which at some point will be merged into ps.map too)


Hamish




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


Fwd: [GRASS-user] i.gensig comment- Signature 1 not positive definite

2010-11-25 Thread Nick Jachowski
I had this same error message sometime in the past few weeks, and I can't
remember exactly what I did, but it was one of the following:

1) Data has to be transformed to 0-255 positive integers
2) You need to apply a mask to the null values
3) Possibly your training area is not large enough

This is skirting your question a bit, as I can't tell you what the error
message means.  But hopefully you can find your way around the problem with
one of the above steps!  And, if you figure it out please let me know so I
can write it down this time...

Best of luck,
Nick

-- Forwarded message --
From: António Rocha 
Date: 2010/11/25
Subject: [GRASS-user] i.gensig comment- Signature 1 not positive definite
To: GRASS user list 


Greetings

For a few training areas I'm obtaining the following comment (which means
that my signature files are created but empty):
Signature 1 not positive definite
What  does this means?

Thanks
Antonio


__ Information from ESET NOD32 Antivirus, version of virus signature
database 5647 (20101125) __

The message was checked by ESET NOD32 Antivirus.

http://www.eset.com


___
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] Help with viewer

2010-11-25 Thread Micha Silver

On 11/25/2010 10:03 AM, Giacomo Piva wrote:

Ciao a tutti,
Hi all,
Does someone know a good GeoTIFF (and other image format) free and 
multi-platform viewer software?




Maybe OpenEV is what you want:
http://openev.sourceforge.net/


Thank you





--
Micha Silver
Arava Development Co. +972-52-3665918
http://surfaces.co.il


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


Re: [GRASS-user] Re: r.water.outlet

2010-11-25 Thread Rich Shepard

On Thu, 25 Nov 2010, skualos wrote:


I think I had the same problem. I was running r.water.outlet and getting a
plain raster. I found out that you must choose your point exactly on a
stream generated by r.watershed. Doing so solve my problem.


  I think many of us learned this lesson the hard way. :-( I found that I
needed to greatly enlarge the raster to place the outlet exactly on the
stream.

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


Re: [GRASS-user] Grid label format in psmap

2010-11-25 Thread Rich Shepard

On Thu, 25 Nov 2010, John A Stevenson wrote:


grid 1000
  color black
  numbers 1
  font helvetica
  fontsize 12
end*


John,

  Have you tried specifying grid 1 (assuming kilometers rather than meters)?
If so, what results did you see in the pdf?

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


[GRASS-user] Re: r.water.outlet

2010-11-25 Thread skualos

Hello.
I think I had the same problem.
I was running r.water.outlet and getting a plain raster.
I found out that you must choose your point exactly  on a stream generated
by r.watershed.
Doing so solve my problem.
-- 
View this message in context: 
http://osgeo-org.1803224.n2.nabble.com/r-water-outlet-tp1865650p5774891.html
Sent from the Grass - Users mailing list archive at Nabble.com.
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] Re: r.water.outlet

2010-11-25 Thread skualos

Hello.
I think I had the same problem.
I was running r.water.outlet and getting a plain raster.
I found out that you must choose your point exactly  on a stream generated
by r.watershed.
Doing so solve my problem.
-- 
View this message in context: 
http://osgeo-org.1803224.n2.nabble.com/r-water-outlet-tp1865472p5774887.html
Sent from the Grass - Users mailing list archive at Nabble.com.
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


[GRASS-user] Grid label format in psmap

2010-11-25 Thread John A Stevenson

Hi,

I have a UTM dataset that I am plotting with psmap.  I would like to 
have grid squares at 1000 m intervals.  I can do this easily with the 
following:


grid 1000
   color black
   numbers 1
   font helvetica
   fontsize 12
end*

*However, on the resulting plot, the gridline labels are truncated into 
kilometers e.g. 498, 7114, instead of 498000, 7114000.  Is there a 
parameter that I can modify somewhere so that I can control output 
format and display the full number.  (Similar to D_FORMAT / 
OUTPUT_DEGREE_FORMAT in GMT, I suppose).


Cheers
John

--
Dr John A Stevenson
RSE/Scottish Government Personal Research Fellow
(co-funded by Marie Curie Actions)
School of GeoSciences
The University of Edinburgh
Grant Institute
West Mains Road
Edinburgh EH9 3JW
Scotland
PHONE: (+44) 131 650 7526
FAX:   (+44) 131 668 3184
email: john.steven...@ed.ac.uk
web: http://www.geos.ed.ac.uk/homes/jsteven5 



The University of Edinburgh is a charitable body, registered in
Scotland, with registration number SC005336.

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


[GRASS-user] i.gensig comment- Signature 1 not positive definite

2010-11-25 Thread António Rocha

Greetings

For a few training areas I'm obtaining the following comment (which 
means that my signature files are created but empty):

Signature 1 not positive definite
What  does this means?

Thanks
Antonio


__ Information from ESET NOD32 Antivirus, version of virus signature 
database 5647 (20101125) __

The message was checked by ESET NOD32 Antivirus.

http://www.eset.com


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


[GRASS-user] Re: loading kc50 palsar datasets in GRASS

2010-11-25 Thread maning sambale
Follow-up on the, from the website:
*2. HH and HV polarization data are kept in BSQ format. Color
composition of the image data is R=HH, G=HV, B=HH/HV and enhanced.

In gdal docs reference to bsq is in the envi format:
http://www.gdal.org/frmt_various.html#ENVI

On Wed, Nov 24, 2010 at 6:21 PM, maning sambale
 wrote:
> Anyone successfully loaded kc50 palsar raw imagery with gdal?
>
>
> -- Forwarded message --
> From: maning sambale 
> Date: Tue, Nov 23, 2010 at 4:00 PM
> Subject: loading kc50 palsar datasets in GRASS
> To: grass-us...@lists.osgeo.org
>
>
> This site provides PALSAR 50m Orthorectified Mosaic Product.
> http://www.eorc.jaxa.jp/ALOS/en/kc_mosaic/kc_mosaic.htm
>
> I tried running gdalinfo to the directory, img and header files, but
> gdal doesn't recognize the format.  There is pdf (below the webpage)
> but I can't find any reference to the file formats that is supported
> by GDAL and GRASS.  Has anyone used this dataset in GRASS?
>
> PS.  This is the only dataset without clouds I can find for my study
> area (for obvious reasons :)).
> --
> cheers,
> maning
> --
> "Freedom is still the most radical idea of all" -N.Branden
> wiki: http://esambale.wikispaces.com/
> blog: http://epsg4253.wordpress.com/
> --
>
>
>
> --
> cheers,
> maning
> --
> "Freedom is still the most radical idea of all" -N.Branden
> wiki: http://esambale.wikispaces.com/
> blog: http://epsg4253.wordpress.com/
> --
>



-- 
cheers,
maning
--
"Freedom is still the most radical idea of all" -N.Branden
wiki: http://esambale.wikispaces.com/
blog: http://epsg4253.wordpress.com/
--
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user


Re: [GRASS-user] Help with viewer

2010-11-25 Thread Achim Kisseler

Well:
Q-GIS?

Achim

Am 25.11.2010 09:03, schrieb Giacomo Piva:

Ciao a tutti,
Hi all,
Does someone know a good GeoTIFF (and other image format) free and
multi-platform viewer software?

Thank you



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


[GRASS-user] Help with viewer

2010-11-25 Thread Giacomo Piva

Ciao a tutti,
Hi all,
Does someone know a good GeoTIFF (and other image format) free and 
multi-platform viewer software?

Thank you


--
G.

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