Hi,
I have a problem with i.atcorr with data aster. I tried it with grass 64
and 7, on ubuntu 15.04 64-bit.
Maybe I'm wrong to set the parameters of the txt? (Is there anyone can tell
me if they are correct? I have attached the file .met)...Otherwise I do not
know.
/10
Hi everyone.
Using the software zygrib (Also seen here http://grass.osgeo.org/wiki/GRIB)
I download the file .grb of the area of interest.
When importing the file in GRASS I get many raster ... ok but each
raster is stretched along the y
What is the problem?
peifer wrote
On 03/07/2012 14:22, Gabriele N. wrote:
Hi everyone.
Using the software zygrib (Also seen here
http://grass.osgeo.org/wiki/GRIB)
I download the file .grb of the area of interest.
When importing the file in GRASS I get many raster ... ok but each
raster is stretched
Ciao.
I.pr modules seem to work but I found a problem: i.pr.features
I saw this post because I have the same problem
i.pr.features -s training=TL_ipr_training_TL_ASTER_2006_02_01_RGB.52
features=aasasssas
ERROR: read_training- Can't open file
TL_ipr_training_TL_ASTER_2006_02_01_RGB.52 for
hamish-2 wrote
Gabriele wrote:
I installed but I have the same error.
With g.extension I can install other add-ons (eg
i.landsat.acca) but not i.pr
but are they actually on the disk somewhere,
slightly in the wrong place?
i.e. if you look in $GRASS_ADDON_PATH/ is there
a bin/ dir in
margherita wrote
I can confirm that installation of i.pr by g.extension seems not to work
properly. So, I suggest you to install it by hand following these steps:
cd your path/grass6_devel/imagery
svn co https://svn.osgeo.org/grass/grass-addons/grass6/imagery/i.pr/
cd i.pr/
make
cd ..
Thanks hamish.
I installed tcl8.5-dev and I followed this guide
http://grass.osgeo.org/wiki/Compile_and_Install_Ubuntu
and, although I got error ((Errors in:
/usr/local/src/grass6_devel/lib/ogsf
/usr/local/src/grass6_devel/lib/nviz
/usr/local/src/grass6_devel/misc/m.nviz.image
Ciao.
I have installed GRASS (6.4.1) from the repository and I can not install
i.pr (add-on).
With g.extension (with the script downloaded from
http://trac.osgeo.org/grass/browser/grass/branches/releasebranch_6_4/scripts/g.extension/g.extension)
and also g.extension present in my installation of
Markus Neteler wrote
On Mon, Jan 30, 2012 at 11:16 PM, Gabriele N. lt;gis.gn@gt; wrote:
Ciao.
I have installed GRASS (6.4.1) from the repository and I can not install
i.pr (add-on).
With g.extension
... unfortunately 6.4.1 is too old in a sense that on Ubuntu problems are
known
Ciao Chetan.
I do not know if I understand your problem.
I advise you to create a mask (which defines the area in respect of which
cut the raster) and call the raster with similar names (eg landsat_1,
landsat_2,).
At this point you should do something like this:
for i in ` g.mlist type=rast
Ciao grass users.
I have two raster.
1) raster 1:
n: 2728570
minimum: -0.791045
maximum: 0.774648
range: 1.56569
mean: 0.138333
mean of absolute values: 0.185515
standard deviation: 0.158265
variance: 0.025048
variation coefficient: 114.409 %
sum: 377452.4525757029
2) raster 2:
n: 2645080
Ciao Marcello.
I thought there was another solution but I adopted the solution proposed by
you.
Thank you very much
Gabriele
--
View this message in context:
http://osgeo-org.1803224.n2.nabble.com/third-raster-as-a-function-of-two-input-raster-r-series-tp6067501p6068331.html
Sent from the
Ciao.
Thanks.
Now, although crude, it works. However, I put exit (or exit 0 exit 1)
at the end but the script does not finish properly.
What do you say?
Thanks
http://osgeo-org.1803224.n2.nabble.com/file/n6057767/v.import_kml.sh
v.import_kml.sh
--
View this message in context:
Ciao Hamish and thanks.
It works ok :)
But I have also tried with my small script (attached) that does not work and
I did not understand why. Maybe, if you want, you can watch it.
Thank you very much.
http://osgeo-org.1803224.n2.nabble.com/file/n6052924/v.import_kml.sh
v.import_kml.sh
Ciao Grass users.
Probably my question is trivial, but I can not fix.
1) I want to import many files. kml (all in the same folder) and then I
would like a merge(with v.patch).
I tried to import v.in.ogr all together maybe I'm wrong.
2) Maybe I can do a merge of all files with ogr2ogr?
Hy Glynn, thanks for the explanation.
Gabriele
--
View this message in context:
http://osgeo-org.1803224.n2.nabble.com/r-mapcalc-3-maps-tp5929440p5952651.html
Sent from the Grass - Users mailing list archive at Nabble.com.
___
grass-user mailing list
For information:
I tried r.reclass.area under windows (before I did it from ubuntu) always
with grass 6.4. Under widows will not give me any warning but seems that
the result is the same.
Gabriele
--
View this message in context:
Ciao grass users.
I have a warning when I launch the command r.reclass.areaaborted.
The map, however, is generated but I'm not sure that everything is correct.
What do you say?
Thanks
r.reclass.area
input=fc_map_1001_int_fl1_NO_6_bin_fl2fl1ec_recode_meno_10_11___15
Excellent Maciek and Glynn :)
Although the logic is different (in the solution of Maciek I do a check
before the NULL value, while in the solution of Glynn I do a check before
non-NULL values), I believe that the result is the same.
One note: with the solution of Glynn (at least using the
Ciao grass users.
I have three maps A, B, C.
I would use r.mapcalc to combine them according to this logic:
Put the values of A and if A is null you put the values of B and if B is
null you put the values of C.
I have tried in various ways without success. .. now I'm trying with a
double IF...
Hi Carlos and thanks.
A and B do not overlap, and I should take the values of A and B also
(always). Then.in areas that are empty, I would then put the values of
C.
I tried withr.mapcalc output = if((isnull(a)),b,if(isnull(b),a),c)
but I lose the values of C
--
View
Hi Sylvain and thanks.
I had not thought of r.patch and I think that works (but I have to
carefully examine the result)
I did so
r.patch input = a, b, c = output output_patch
Now I want to do it with r.mapcalc:)
--
View this message in context:
Hi Daniel, Hi Frank.
I tried to export in different ways (as GeoTIFF, as float etc. ..) and
import it into arc in different ways. At this point, I think it is a
limitation of arc someone knows something?
Gabriele
--
View this message in context:
Hy Chetan.
Look here,
http://grass.osgeo.org/wiki/LANDSAT
maybe it can be useful
Ciao
Gabriele
--
View this message in context:
http://osgeo-org.1803224.n2.nabble.com/Importing-of-raster-data-doesn-t-seem-to-work-through-GUI-tp5805160p5805396.html
Sent from the Grass - Users mailing list
Hi.
My advice is to use linux.
In particular, I use Ubuntu (Linux human:)). However GRASS runs well on wind
but I think that using Linux has many more advantages. For example, through
the repository every time there is an updated version of GRASS you can get
it automatically. Morover (if I'm not
Hi.
I am working with colleagues that working on ArcGIS. It 's okay but now I
have a problem with esri grid.
I am using to export from GRASS r.out.arc, but the values of the cells in
arcgis are different from what I read in GRASS.
The raster has values with 13 digits.
Here's how I export:
Hy everyone.
I have a problem with a reprojection.
I have a global raster (a grid esri) that I imported (by choosing the file
w001001.adf) in one location. I created the location with a reference system
that has these parameters (in a file prj.adf):
Projection Mollweide
Units METERS
Zunits NO
I found some other posts I'm getting close to the solution (I get good
results):
gdalwarp -s_srs '+proj=goode +ellps=sphere +lon_0=30E +towgs84=0,0,0,0,0,0,0
+x_0=3335846.22854 +y_0=2' -t_srs '+proj=latlong +ellps=wgs84
towgs84=0,0,0,0,0,0,0' -te 7 33 34 59 map.tif map_wgs84.tif
I do
Hy everyone.
Two problems:
1) I can not import a file erdas why is not supported.
gdalinfo image.img
ERROR 4: `image.img' not recognised as a supported file format.
gdalinfo failed - unable to open 'image.img'.
2) Projection: Interrupted Goode Homolosine
I could not find the EPSG
Hi Martin.
Yes, strange ... I tried also with QGIS, and obviously does not work ... : /
This file is strange:) .. also because I can not understand how to create
the location with the projection Interrupted Goode Homolosine...boooh
Gabriele
--
View this message in context:
Glynn Clements wrote:
The same result can be achieved with one command:
r.mapcalc output = A * 1.0 + B
The multiply will cause A to be converted to double to match the other
argument (floating-point constants are double unless an explicit f
suffix is given), the multiply
Hy Daniel.
I do not think that r.cross can solve the problem.
With r.cross I could get a concatenation of label:
For example, with values A = 12 B = 3323 C = 1 the result with the label
category 12, category 3323, category 1
However, if you know a chance I can prove otherwise.
As a result, the
yes, maybe with grep and r.cross ... but then you think you can convert the
result of r.cross as cell value?
Gabriele
--
View this message in context:
http://osgeo-org.1803224.n2.nabble.com/Values-from-multiple-raster-to-link-together-tp5741032p5766761.html
Sent from the Grass - Users mailing
I have a new problem, this time I would like to concatenate two maps A and B.
A is with 8 digits, B is with 4 digits. So I do as usual:
r.mapcalc output = (( A * 1) ( B ))
But I get strange values with negative values. I did several tests but
without resolve. Perhaps r.mapcalc may have
Hy Glynn.
Thanks for the interesting explanation.
I did then so:
r.mapcalc A_new = (double(A)) and r.mapcalc B_new = (double(B))
and then
r.mapcalc output = ((A_new * 1)+(B_new ))
I need to check the result, but it seems that the result is right.
Thank you very much
Gabriele
--
Daniel Victoria wrote:
It's not r.mapcalc fault. I imagine it's probably a limitation due to
the number of bits in your image. For instance, 8 bit images can have
up to 256 values. 16 bit images can have 65,535. So, if you have too
many values, you will have problems. Not sure how to
Sorry, I have not explained well.
The problem was related to the value of C map. That is, I take the example
above. If the value of C map is equal to 1, I get a result - A =12 B=3323
C=1 =12332301 (8 digit number) and instead could be useful 1233231 (number
7-digit).
However I checked and I need
Hello everyone.
I have to solve a problem already discussed here:
http://osgeo-org.1803224.n2.nabble.com/r-mapcalc-concatenates-values-td1884970.html#a1884973
In summary, I have for example 3 maps A, B, C.
I need to combine the values of the overlapping cells. For example if I have
a cell in the
Excellent Alex:). The maps are as follows:
A - value 0, 1, 2, ... 19
B - value 0, 1001, 1002, 1003 3351
C - value from 0, 1, 2, 3, 4 .. 15
So, I did as you suggested:
r.mapcalc output = (( A *100)+( B *100)+ C )
right?
And if I had decimal numbers or strings instead of numbers?
Hi Markus.
I have installed grass-dev but now there is another problem (there is a
problem of permissions for the folder creation and does not solve if I
create the folder before launching g.extension):
g.extension i.landsat.toar
Fetching i.landsat.toar from GRASS-Addons SVN (be patient)...
A
Markus, still does not work.
I saw some differences compared to the old script. To understand ... which
is the part that has been changed to fix this?
GRASS 6.4.0 (utm_wgs84):~ g.extension i.landsat.toar
Fetching i.landsat.toar from GRASS-Addons SVN (be patient)...
A
I have the same problem (ubuntu 10.10, with repo ubuntugis)
GRASS 6.4.0 (utm_wgs84):~ g.extension i.landsat.toar
Fetching i.landsat.toar from GRASS-Addons SVN (be patient)...
Ai.landsat.toar/landsat_set.c
Ai.landsat.toar/local_proto.h
Ai.landsat.toar/main.c
A
Thanks Robbie.
I try to take a look.
Gabriele
--
View this message in context:
http://osgeo-org.1803224.n2.nabble.com/R-Moran-index-getis-e-ord-tp5588509p5599570.html
Sent from the Grass - Users mailing list archive at Nabble.com.
___
grass-user
Hello everyone.
I have a raster which is the result of calculating an index of remote
sensing (with decimals).
I would like to calculate the Moran's index and the index of Getis and Ord
To do this I installed R (following
http://grass.fbk.eu/statsgrass/grass6_r_install.html).
Then I installed
Hi to everyone.
I have a vector of points and I want to connect each point to a photo and a
link (to open a photo or a web page with the browser firefox).
I have the attributes in SQLite. With sqlite manager I added a column 'image'
of type blob. I uploaded an image but does not work. When I
Hi.
I'm working with grass64 on ubuntu 8.10.
I have created a sqlite db.
1)
I have a shape of municipalities that have problems of topology. In addition,
some municipalities have more non-adjacent areas. To not have the
municipalities 'split' I tried to import without topological corrections
: GRASS user list grass-user@lists.osgeo.org, Gabriele N.
gis...@libero.it
Sent: Sunday, March 1, 2009 9:10:14 PM GMT +00:00 GMT Britain, Ireland,
Portugal
Subject: Re: [GRASS-user] GRASS 7 and SQLITE o POSTGRES
On Sat, Feb 21, 2009 at 1:42 PM, Moritz Lennert
mlenn...@club.worldonline.be wrote
am
forced to import data from Postgres. right?
Thanks for any advice.
Moritz Lennert wrote:
On 21/02/09 14:25, Gabriele N. wrote:
Hi Moritz.
I have read here http://grass.osgeo.org/wiki/GRASS_7_ideas_collection
Database
establish SQLite as default DBMI driver (DBF is too limited). done
: Gabriele N. gis...@libero.it
To: grass-user@lists.osgeo.org
Sent: Friday, February 20, 2009 5:38:48 PM GMT +00:00 GMT Britain,
Ireland, Portugal
Subject: [GRASS-user] GRASS 7 and SQLITE o POSTGRES
Hello.
I work with colleagues sharing the files (shapes .. etc) that are on a NAS
server.
I
Hello.
I work with colleagues sharing the files (shapes .. etc) that are on a NAS
server.
I use ubuntu 8.10 (with GRASS locally) and my colleagues use windows.
I would do this:
- A GIS with all these data
- Building a database (sqlite or postgres/postgis/ on NAS)
- An excellent organization of
Hi.
All my scripts in bash work with tcltk.
Now, I take as example my 2 scripts: v.to.dem.sh and v.impact.sh.
With WxPython all the scripts were not working with grass 63. With Grass64,
instead, some scripts works (es v.impact.sh) and some other scripts does not
works (es.v.to.dem.sh). But I do
Screenshot-Error+in+command+execution.png
Glynn Clements wrote:
Gabriele N. wrote:
All my scripts in bash work with tcltk.
Now, I take as example my 2 scripts: v.to.dem.sh and v.impact.sh.
With WxPython all the scripts were not working with grass 63. With
Grass64,
instead, some scripts
Martin Landa-2 wrote:
Hi,
2008/7/5 Gabriele N. [EMAIL PROTECTED]:
I use grass63 with wxpython and I do not know how to use some scripts
that I
have in bash. The graphical interface opens but when I try to launch the
script I get the following error:
GRASS 6.3.0 (gb2):~ v.topo_5
Hy Glynn.
Ah ... I work on ubuntu 8.04.
But what could be the commands in the script that GRASS not 'recognize'?
Compared to version 63 with the 64 now run more scripts (always bash).
All the scripts still work well in tcltk.
?
Thanks
Gabriele
Glynn Clements wrote:
Gabriele N. wrote
Hi everyone.
I have two questions
1)
I have some raster (with the same format and are converted from vector). The
cells have the value of CAT. So the value = CAT
With r.mapcalc I would define a map of output that I give, as the value for
each cell, indications of the different CAT. For
Hello.
Following this recent post:
http://www.nabble.com/Re%3A-grass6.4-develbrunch-revision-32726%3A-undefined-reference-to-%60G_no_gisinit%27-td18947949.html,
I installed so grass64 (ubuntu 8.04):
I've just built gdal 1.5.2 (without grass)
./configure --with-ecw=yes --without-grass
and
Moritz Lennert-2 wrote:
On 09/07/08 10:42, Gabriele N. wrote:
Hi Moritz.
Moritz Lennert-2 wrote:
On 08/07/08 22:33, Gabriele N. wrote:
Hello.
I want to create a script bash (GUI) to export many vector in shapefile
having the opportunity to decide which export layer (0, 1, 2
Hello.
I want to create a script bash (GUI) to export many vector in shapefile
having the opportunity to decide which export layer (0, 1, 2) and
type=line,area ecc..
I would also like to know if you can use the characters as' * 'to be able to
export all instance vector starting with' nome_vector
Hello
1)
I use grass63 with wxpython and I do not know how to use some scripts that I
have in bash. The graphical interface opens but when I try to launch the
script I get the following error:
GRASS 6.3.0 (gb2):~ v.topo_5
$ v.topo_5 percorso=/home/gab/prova.shp nome_vett=foglio num=3
Hi Martin.
I like that.
Thank you very much
Gabriele
Martin Landa-2 wrote:
Hi,
Upgrading to grass64 (devbr6) should fix it.
generally speaking, if you want to test wxGUI and report its bugs,
upgrade to grass64.svn is necessary (wxpython code in grass63 is quite
out-dated, just
and describing small towns for you case.
Regards,
Christian.
Am Sonntag, den 15.06.2008, 08:43 -0700 schrieb Gabriele N.:
Hi all
I have 2 vector polygonal; one representing the buildings and another
that
represents the infrastructure.
I think to put them together with v.patch.
I have to create
Hi all
I have 2 vector polygonal; one representing the buildings and another that
represents the infrastructure.
I think to put them together with v.patch.
I have to create another vector polygonal, which is to enclose all the
polygons that are close to each other within 50 metres. I need to
Hello all.
I have two problems.
1) v.dissolve?
it is possible to dissolve the polygons (nonadjacent) who have the same
attribute obtaining a multipolygon?
2) v.type
I noticed that to get points (representing the centroid) from polygons must
change the order by hand ... type = point, centroid
hamish_b wrote:
Gabriele wrote:
I have a theme of a polygonal rivers. I have to locate the center line
as a line that identifies the river. Obviously the edges of the river
are not parallel and can not use v.parallel.
..
I would locate the midline of the river.
So I would like to
Hi list
--
View this message in context:
http://www.nabble.com/center-line-with-v.generalize--.-tp15062970p15062970.html
Sent from the Grass - Users mailing list archive at Nabble.com.
___
grass-user mailing list
grass-user@lists.osgeo.org
Hi list.
I have a theme of a polygonal rivers. I have to locate the center line as a
line that identifies the river. Obviously the edges of the river are not
parallel and can not use v.parallel.
I am making attempts with v.generalize but unable to solve.
Suggestions?
Thanks
Gabriele
--
I would locate the midline of the river.
So I would like to transform the polygon with the middle line.
Thanks
--
View this message in context:
http://www.nabble.com/center-line-with-v.generalize---tp15063073p15063622.html
Sent from the Grass - Users mailing list archive at Nabble.com.
Moritz Lennert-2 wrote:
On 24/01/08 14:12, Gabriele N. wrote:
Hi list.
I made a script (attached) for the definition of interactive text file
(containing the coordinates input and output) to be loaded later in
v.transform.
Could you explain how this is different from the existing
68 matches
Mail list logo
