[GRASS-user] again Raster names
Hi, in a script, which I executed under GRASS 6.4 under Ubunto, I included commands like r.in.gdal input=/media/.../filename output=surfac_ereflection. 1. An error message illegal filename surfac_ereflection.1, illegal character showed up. When I typed this command into the window underneath the output window, it worked fine. I also run the script using data on a FAT32 formatted partition as well as data in the /home directory. The files to be imported are MODIS data. Best regards Niels - Dr. Niels Thevs Institute of Botany and Landscape Ecology Greifswald University Grimmer Strasse 88 17487 Greifswald Germany Tel.: +49-3834-86-4132 Fax: +49-3834-86-4114 === WEB-Mailer der Uni-Greifswald ( http://www.uni-greifswald.de/ ) === ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] v.extract - sqlite driver problem
Perfect. Setting the bad point null, and exporting without table or topology (3D) proved to be the fastest solution. Thanks for the suggestion. Much thanks, Mark On Fri, Feb 5, 2010 at 2:27 AM, Hamish hamis...@yahoo.com wrote: Mark: I have a bunch of points (4+ million) with one outlier that needs removed. perhaps use r.mapcalc or r.reclass to filter out the bad point? r.mapcalc clean = if(map 99, null(), map) or r.reclass 99 thru 999 = NULL * = * and finally r.null setnull=bad_value may be the simplest. ... or d.rast.edit ... Hamish ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Radiometric Correction @GRASS6.4
Thanks. I have a portuguese friend that can help me :) And regarding i.atcorr code, is based on Mauro's code or the original one? Nikos On Thu, Feb 4, 2010 at 5:19 PM, Daniel Victoria daniel.victo...@gmail.comwrote: On Thu, Feb 4, 2010 at 10:36 AM, Nikos Dumakis nikosdu1...@gmail.com wrote: Hello Markus. Thanks for your answers As far as I saw, there are no GRASS functions to perform Radiometric (by the way, at i.atcorr manual one of the links (INPE) is dead) The INPE link was for Mauro A. Homem Antunes home page and now he is at Federal Rural University of Rio de Janeiro. He still maintains a web page where you can download his 6s fortran code (web page in Portuguese though). The link is: http://www.ufrrj.br/institutos/it/deng/mauro/downloads.htm Cheers Daniel ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Region Definition in GRASS
Ok but for instance. If my region is defined with a Spatial Resolution of 1kilometer and, if I do some r.mapcalc with two LANDSAT images, the output image will have which one of the resolutions? the Region or the 2 landsat images? Thank you Glynn for the overview. And thank you Hamish for the weblinks On Thu, Feb 4, 2010 at 7:14 PM, Glynn Clements gl...@gclements.plus.comwrote: Pedro Roma wrote: I'm a kind of a new GRASS user and I'm a bit lost on regarding something. REGION (or spatial extent). For instance the Location North-Carolina. When the location was created a certain extent, projkection a Spatial Resolution was defined. But, inside this location there is a lot of data with different extent and regional cover. And, inside GRASS I can re-define region properties (g.region) and use a certain data to define region. And, it's possible to have a region with a small extent than a vectorial that exists in the Location. My question is, what is the region and what is its purpose? Because I don't see any impact of region over my image/data processing or image visualization. Or is it? The region is mostly used for raster processing and graphical display. Most commands which read or write raster maps will use the bounds and resolution from the current region. Any input maps will be cropped, padded and/or resampled according to the region settings. Display commands translate geographic or cartographic coordinates to display coordinates such that the current region just fits inside the display window. Vector commands are normally unaffected by the current region. -- Glynn Clements gl...@gclements.plus.com ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Region Definition in GRASS
Pedro, If you run a r.mapcalc, the extent and resolution will follow the current region. So, check g.region -p before run r.mapcalc. You can use g.region rast=MapLandsat to set the current region to use MapLandsat dimensions as reference. After you can reset your extention and resolution using g.region -d . By the way, read g.region --help. best milton 2010/2/5 Pedro Roma pedroroma1...@gmail.com Ok but for instance. If my region is defined with a Spatial Resolution of 1kilometer and, if I do some r.mapcalc with two LANDSAT images, the output image will have which one of the resolutions? the Region or the 2 landsat images? Thank you Glynn for the overview. And thank you Hamish for the weblinks On Thu, Feb 4, 2010 at 7:14 PM, Glynn Clements gl...@gclements.plus.comwrote: Pedro Roma wrote: I'm a kind of a new GRASS user and I'm a bit lost on regarding something. REGION (or spatial extent). For instance the Location North-Carolina. When the location was created a certain extent, projkection a Spatial Resolution was defined. But, inside this location there is a lot of data with different extent and regional cover. And, inside GRASS I can re-define region properties (g.region) and use a certain data to define region. And, it's possible to have a region with a small extent than a vectorial that exists in the Location. My question is, what is the region and what is its purpose? Because I don't see any impact of region over my image/data processing or image visualization. Or is it? The region is mostly used for raster processing and graphical display. Most commands which read or write raster maps will use the bounds and resolution from the current region. Any input maps will be cropped, padded and/or resampled according to the region settings. Display commands translate geographic or cartographic coordinates to display coordinates such that the current region just fits inside the display window. Vector commands are normally unaffected by the current region. -- Glynn Clements gl...@gclements.plus.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] again Raster names
Niels Thevs wrote: in a script, which I executed under GRASS 6.4 under Ubunto, I included commands like r.in.gdal input=/media/.../filename output=surfac_ereflection. 1. An error message illegal filename surfac_ereflection.1, illegal character showed up. When I typed this command into the window underneath the output window, it worked fine. I also run the script using data on a FAT32 formatted partition as well as data in the /home directory. The files to be imported are MODIS data. How are you creating the script? Are you sure that invalid characters (e.g. carriage return, non-breaking space, etc) aren't being introduced? Try to create a minimal script which demonstrates the problem, then post it as a binary attachment. -- Glynn Clements gl...@gclements.plus.com ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Region Definition in GRASS
Pedro Roma wrote: Ok but for instance. If my region is defined with a Spatial Resolution of 1kilometer and, if I do some r.mapcalc with two LANDSAT images, the output image will have which one of the resolutions? the Region or the 2 landsat images? The region. Most raster modules (including r.mapcalc) generate maps with the bounds and resolution of the current region. Input maps are cropped, padded and/or resampled to the current region, so all processing within the module takes place on the grid defined by the current region. The main exceptions are the r.in.* modules, which import files cell-for-cell, using the bounds specified in the file's metadata if present. The current region is ignored. Modules which perform their own resampling (e.g. r.resamp.*, r.proj) create output maps according to the current region but read their input maps without cropping or resampling. -- Glynn Clements gl...@gclements.plus.com ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Raster Map subset
Is there a way to get a subset of a larger raster map. I have zoomed into an area and would like to make a raster map only containing that area. Can I do this? I am going to use this as input into R for some topmodel runs. thanks, -- Stephen Sefick Let's not spend our time and resources thinking about things that are so little or so large that all they really do for us is puff us up and make us feel like gods. We are mammals, and have not exhausted the annoying little problems of being mammals. -K. Mullis ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Raster Map subset
Hi Stephen, I suggest you use spgrass6 package on R, set your g.region using system(g.region n= s= w= e=) on R and use readRAST6 to load your map as SGDF. Of course, if you are under windows, you know that depending the size of your grid you will get troubles :-) best milton 2010/2/5 stephen sefick ssef...@gmail.com Is there a way to get a subset of a larger raster map. I have zoomed into an area and would like to make a raster map only containing that area. Can I do this? I am going to use this as input into R for some topmodel runs. thanks, -- Stephen Sefick Let's not spend our time and resources thinking about things that are so little or so large that all they really do for us is puff us up and make us feel like gods. We are mammals, and have not exhausted the annoying little problems of being mammals. -K. Mullis ___ 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] Raster Map subset
Wonderful- thank you very much. Stephen On Fri, Feb 5, 2010 at 6:31 PM, Milton Cezar Ribeiro miltinho.astrona...@gmail.com wrote: Hi Stephen, I suggest you use spgrass6 package on R, set your g.region using system(g.region n= s= w= e=) on R and use readRAST6 to load your map as SGDF. Of course, if you are under windows, you know that depending the size of your grid you will get troubles :-) best milton 2010/2/5 stephen sefick ssef...@gmail.com Is there a way to get a subset of a larger raster map. I have zoomed into an area and would like to make a raster map only containing that area. Can I do this? I am going to use this as input into R for some topmodel runs. thanks, -- Stephen Sefick Let's not spend our time and resources thinking about things that are so little or so large that all they really do for us is puff us up and make us feel like gods. We are mammals, and have not exhausted the annoying little problems of being mammals. -K. Mullis ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user -- Stephen Sefick Let's not spend our time and resources thinking about things that are so little or so large that all they really do for us is puff us up and make us feel like gods. We are mammals, and have not exhausted the annoying little problems of being mammals. -K. Mullis ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Raster Map subset
On Fri, 5 Feb 2010, stephen sefick wrote: Is there a way to get a subset of a larger raster map. I have zoomed into an area and would like to make a raster map only containing that area. Can I do this? I am going to use this as input into R for some topmodel runs. Stephen, Yes, you can. This subject has come up several times in the past month or so so the archives have a lot of options for you. One to look at closely is v.in.region. After zooming to the area you want, run v.in.region and save that vector map. You can then use it with the g.region command to limit the displayed bounds of your raster map. You can use r.patch to make a new raster map that's the size of the newly created vector map. Are you aware of r.topmodel within GRASS? Rich ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Raster Map subset
I am aware of r.topmodel in GRASS, but I am much more familiar with R than GRASS, and the GRASS documentation is not as comprehensible as that for R -- This is just because of my own bias and does not reflect on anybody but me. I would love to try it in GRASS but I think I need a tutorial. Any help would be greatly appreciated with this. Stephen On Fri, Feb 5, 2010 at 6:36 PM, Rich Shepard rshep...@appl-ecosys.com wrote: On Fri, 5 Feb 2010, stephen sefick wrote: Is there a way to get a subset of a larger raster map. I have zoomed into an area and would like to make a raster map only containing that area. Can I do this? I am going to use this as input into R for some topmodel runs. Stephen, Yes, you can. This subject has come up several times in the past month or so so the archives have a lot of options for you. One to look at closely is v.in.region. After zooming to the area you want, run v.in.region and save that vector map. You can then use it with the g.region command to limit the displayed bounds of your raster map. You can use r.patch to make a new raster map that's the size of the newly created vector map. Are you aware of r.topmodel within GRASS? Rich ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user -- Stephen Sefick Let's not spend our time and resources thinking about things that are so little or so large that all they really do for us is puff us up and make us feel like gods. We are mammals, and have not exhausted the annoying little problems of being mammals. -K. Mullis ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Raster Map subset
On Fri, 5 Feb 2010, stephen sefick wrote: I am aware of r.topmodel in GRASS, but I am much more familiar with R than GRASS, and the GRASS documentation is not as comprehensible as that for R -- This is just because of my own bias and does not reflect on anybody but me. I would love to try it in GRASS but I think I need a tutorial. Any help would be greatly appreciated with this. Stephen, Well, don't look to me for that tutorial. I need to explore it myself for my current project. Rich ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user