[GRASS-user] loading kc50 palsar datasets in GRASS
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/ -- ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] New page IGNF register and shift grid NTF/RGF93
Thibault, closely concerned by this item, I'll soon try to have a glance at your page. I have a set of test-data available in both NTF and RGF93 systems, and Circé to validate back & forth shiftings. Thank you for this contribution ! Vincent. Le lundi 22 novembre 2010 à 21:10 +0100, Thibault Lemaitre a écrit : > Hello, > > I wrote a new page on the wiki (url : > http://grass.osgeo.org/wiki/IGNF_register_and_shift_grid_NTF/RGF93). > It's a tutorial to well use the French projection systems. Could you > look and make some feedbacks, comments, modifications if it's necessary? > People who want to do a French translation are welcome (I will begin > this soon). > > Grassly, > Thibault > ___ > 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] Not possible to display vector data
Greetings I have a vector map with the following characteristics (v.info) v.info test01_LU ++ | Layer: test01_LU | | Mapset: Regional | | Location:Localize | | Database:D:\data | | Title: | | Map scale: 1:1 | | Map format: native | | Name of creator: monica | | Organization: | | Source date: Tue Nov 23 11:19:55 2010 | || | Type of Map: vector (level: 2) | | | | Number of points: 0 Number of areas: 76832 | | Number of lines:0 Number of islands:7949 | | Number of boundaries: 172475 Number of faces: 0 | | Number of centroids:76816 Number of kernels:0 | | | | Map is 3D: No | | Number of dblinks: 1 | | | | Projection: Transverse Mercator | | N: 276083.7674S: -300404.804 | | E: 162129.0811W: -119191.4075 | | | | Digitization threshold: 0 | | Comments: | | | ++ (Tue Nov 23 11:41:42 2010) Command finished (0 sec) When I try to display it in a Map Window my wxGUI freezes and it gets non responsive. What might be happen? DO you suggest to "wait"? ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: Values from multiple raster -
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 value of cell in the raster output should be 12332301 Ciao Gabriele -- View this message in context: http://osgeo-org.1803224.n2.nabble.com/Values-from-multiple-raster-to-link-together-tp5741032p5766350.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
Re: [GRASS-user] Re: Values from multiple raster -
Ok, I misunderstood your problem. You want the concatenated values, not just the map of the different crossings. Maybe use r.cross and then some grep magic in the raster category file Daniel On Tue, Nov 23, 2010 at 10:02 AM, Gabriele N. wrote: > > 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 value of cell in the raster output should be 12332301 > > Ciao > > Gabriele > -- > View this message in context: > http://osgeo-org.1803224.n2.nabble.com/Values-from-multiple-raster-to-link-together-tp5741032p5766350.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 mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Not possible to display vector data
Hi, 2010/11/23 Monica Buescu : > When I try to display it in a Map Window my wxGUI freezes and it gets non > responsive. What might be happen? DO you suggest to "wait"? version? OS? Martin -- Martin Landa * http://geo.fsv.cvut.cz/~landa ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Re: Help/manual on Georectifier 6.4.0
Jenny Turner wrote: >> >> > active mapset shall be which one? the one where I have my reference >> > images >> > (already georrectified) or mapset where I have my not-georreferenced >> > images)? >> >> The active location is the target location with a mapset with your >> reference images (indicated by the very first question of the >> georectifier where you have to specify the source and not the target >> location). > > My question was not related with Lcoation (because both maps > (not-georreference and georreferenced) are in the same Location but in > different mapsets > In my example: Mapset:Local: not-georreference and to be georreferenced (How > do you call this?) I chose 'source' and > Mapset: Global: georreferenced (how do you call this?) ... 'target' because the source map can be georeferenced as is probably the case in your example and you want to improve georeferencing of the source map, i.e. you want to georectify some images. Therefore the terms 'un-georeferenced' and 'georeferenced' are IMHO a bit misleading. > In the first panel it's asked about Source Location and source Mapset. You > mean that here I shall define where is my not-georreferenced image? In this > case mapset:Local? Yes. > Second panel asks for a group of images to be goerrectified. > Third asks for an image to be displayed. in this case is a georreferenced or > not-georreferenced image? Yes, one of the images you want to georectify. In the new Georectifier, the selection is restricted to the images in the chosen group to avoid confusion, and the fields are labelled source and target where appropriate. Also note that you can change the source image to be displayed through the settings dialog in the new Georectifier. > So, in the normal Map diSplay window I will put an already georreferenced > image right? in this case from Mapset Global? Yes. >> >> Source and target location can be identical. The manual of >> the new georectifier [1] may be able to answer further questions since >> all the functionality of the 6.4.0 georectifier is included in the >> georectifier in 6.4.1 and above. >> >> Markus M >> >> [1] >> http://grass.osgeo.org/grass64/manuals/html64_user/wxGUI.GCP_Manager.html >> >> > Thanks >> > Jenny >> > On Thu, Nov 18, 2010 at 11:22 PM, Jenny Turner >> > >> > wrote: >> >> >> >> Greetings >> >> I'm writting a tutorial of Georrectify tool (from 6.4.0 not latest >> >> launch) >> >> for my students. >> >> I was only able to find >> >> this:http://gama.fsv.cvut.cz/~landa/grass/swf/georect.html. Is there >> >> any >> >> other Wiki or document? >> >> Thanks >> >> Jenny >> > >> > ___ >> > 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] Re: Values from multiple raster -
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 list archive at Nabble.com. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] how to install libLAS on MAC for use in GRASS and maybe R?
Hi William, thanks for your advice and great website! I will try to compile liblas via cmake for mac osx. I have downloaded and installed cmake-2.8.3 and now I downloaded the boost package (boostI-jam-3.1.18-1-macosxx86_64.tgz). Following the http://www.boost.org/doc/libs/1_45_0/more/getting_started/unix-variants.html I just need to unzip on the right location. Which location would you recommend? Will boost install cause potential problems? Should I do something different than described here http://liblas.org/compilation.html#optional-libraries to configure the optional libraries (gdal and geotiff installed before the grass install via your website)? All tips are welcome. Having the liblas library options for grass and R would be great as I am working with terrestrial lidar data (las 1.1 and 1.2) and these files are quite large. Best regards, Jan On Mon, Nov 22, 2010 at 3:48 PM, William Kyngesburye wrote: > > On Nov 22, 2010, at 8:12 AM, Jan Verbesselt wrote: > > > Hi all, > > > > We are working with terrestrial lidar data which can be exported to LAS > > format (1.1 and 1.2). I would like to explore the options to analyse the > > data via GRASS, R, Python. > > I have seen the Lidar GRASS wiki > > (http://grass.osgeo.org/wiki/LIDAR#Micro-tutorial) which is very nice and > > also the message of Hamish > > (http://lists.osgeo.org/pipermail/grass-user/2009-February/048716.html) > > > > Now I am not a very experienced with GRASS and still need to install the > > libLAS library on my Mac (10.6.5). I have seen the following instructions: > > http://liblas.org/compilation.html#using-xcode-on-os-x on the liblas > > website to compile the package for Mac. > > Looks to me like all that is to set it up for testing and debugging. If you > just want to use it, you could probably just configure it from cmake to build > with make instead of xcode (follow the linux instructions). > > > Before I start however is there an easier way to do this? Maybe via the > > Fink (http://www.finkproject.org/) or an mac installer? > > > > I am happy to update the GRASS wiki with info about liblas install on mac > > if I manage to make it work. > > > Not much of a fan of Fink or Macports myself. Lots of overhead, and it can > interfere with compiling other software from scratch. They don't use my > frameworks, so if you are using my GRASS binaries you will be duplicating all > those libraries (GDAL, Geotiff and their dependencies). Boost could be > trouble, it looks like liblas uses boost libraries (boost is by default a > header-only library). > > Making a package has crossed my mind, though I was waiting for it to mature > (this was a while ago) and I've been distracted by other things... > > At least, when I start looking at it I will add build instructions (hopefully > simpler) to my site. > > - > William Kyngesburye > http://www.kyngchaos.com/ > > "Mon Dieu! but they are all alike. Cheating, murdering, lying, fighting, and > all for things that the beasts of the jungle would not deign to possess - > money to purchase the effeminate pleasures of weaklings. And yet withal > bound down by silly customs that make them slaves to their unhappy lot while > firm in the belief that they be the lords of creation enjoying the only real > pleasures of existence > > - the wisdom of Tarzan > > ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: Values from multiple raster -
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 problems with numbers greater than 10 digits Thanks Gabriele -- View this message in context: http://osgeo-org.1803224.n2.nabble.com/Values-from-multiple-raster-to-link-together-tp5741032p5767794.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
Re: [GRASS-user] Re: Values from multiple raster -
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 avoid this in your case but, maybe you will need to work with categories after all. Are you planing on doing any math operations in the concatenated values or is it more for displaying pourposes? Cheers Daniel http://en.wikipedia.org/wiki/Integer_(computer_science) On Tue, Nov 23, 2010 at 4:34 PM, Gabriele N. wrote: > > 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 problems with numbers greater > than 10 digits > > > Thanks > > Gabriele > -- > View this message in context: > http://osgeo-org.1803224.n2.nabble.com/Values-from-multiple-raster-to-link-together-tp5741032p5767794.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 mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: [OSGeo-Discuss] Re: Fedora Geo/GIS Spin
FYI - Thanks to Viji, there is now GRASS 6.4.0 also in Fedora available. On Tue, Nov 23, 2010 at 8:03 PM, Viji V Nair wrote: > Update: > > Grass has been update to the latest upstream version, 6.4.0, which is > currently available for rawhide (F15) > > http://koji.fedoraproject.org/koji/buildinfo?buildID=206173 > > The same will be available for EL6/F13/F14 in a while. > > Please test and report bugs, if any > > Thanks > Viji > ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Re: Values from multiple raster -
Gabriele N. wrote: > 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 problems with numbers greater > than 10 digits CELL maps contain signed, 32-bit values using sign-bit representation, so the range is -2147483647 to 2147483647 inclusive (i.e. -(2^31-1) to 2^31-1). So you aren't going to be able to store 12-digit numbers in a CELL map. DCELL maps use the C "double" type, which on most systems is IEEE-754 double precision, which has a 53-bit mantissa, plus a separate sign bit. That can exactly represent any integer between -9007199254740992 and 9007199254740992 inclusive (i.e. almost 16 decimal digits). However, depending upon what you want to do with the result, this may not be an option, as not all modules can work with floating-point data. -- Glynn Clements ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: Values from multiple raster -
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 -- View this message in context: http://osgeo-org.1803224.n2.nabble.com/Values-from-multiple-raster-to-link-together-tp5741032p5768432.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: Values from multiple raster -
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 avoid this in > your case but, maybe you will need to work with categories after all. > Are you planing on doing any math operations in the concatenated > values or is it more for displaying pourposes? > Hi Daniel... Maybe it would be useful only concatenation but I would also have the possibility to apply tools such as r.statistics and other similar that work on the value of cells. Gabriele -- View this message in context: http://osgeo-org.1803224.n2.nabble.com/Values-from-multiple-raster-to-link-together-tp5741032p5768618.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] GRASS user map - added twice
Hi everyone! Today I tried adding myself to the Worldwide GRASS GIS usermap and succeeded in doing so. However due to internet connection issues I had to refresh the page. It resulted in my name appearing and counting twice(earlier there were 1421 users, now it is 1423). I find no ways to delete the extra point. Therefore I would like to know how to go about regarding this issue. I am even ready to contact the administrator for the map site if someone could provide me the email id. Regards, Chethan S. -- Chethan S. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Re: Values from multiple raster -
Gabriele N. wrote: > 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" ))" 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 will yield a double result, and the addition will convert B to double to match this. -- Glynn Clements ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] GRASS user map - added twice
Hi Chethan S., I can look in to this. Markus On 11/24/10, Chethan S wrote: > Hi everyone! > > Today I tried adding myself to the Worldwide GRASS GIS usermap and succeeded > in doing so. However due to internet connection issues I had to refresh the > page. It resulted in my name appearing and counting twice(earlier there were > 1421 users, now it is 1423). I find no ways to delete the extra point. > Therefore I would like to know how to go about regarding this issue. I am > even ready to contact the administrator for the map site if someone could > provide me the email id. > > Regards, > > Chethan S. > > -- > Chethan S. > ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user