Re: [GRASS-user] reproject from lat/long to utm
try this: r.proj input=dem dbase=/work/projects/ location=Willamette mapset=PERMANENT On Wed, Dec 2, 2009 at 2:45 PM, maning sambale wrote: > On Wed, Dec 2, 2009 at 2:28 PM, Janet Choate wrote: >> I tried using the full path (i.e. the PERMANENT mapset lives within the >> /work/projects/Willamette location), but get the 'not found' error. >> is there any other command i need to run prior to r.proj, or any additional >> options to include, or how to list the path other than the way i am? > > just to check from your lat/long mapset, run g.gisenv > It should show you the location, mapset and GRASS database. > > This is what you use in r.proj parameters > >> the original data i imported was a dem from usgs national elevation dataset >> (ned) seamless data distribution system (sdds). i had to export it from ned >> as a geotiff to bring it into grass. i tried exporting the data from ned as >> an arcgrid and importing into grass with r.in.gdal, but was unable to bring >> it into grass. >> >> so if anyone can tell me the steps to bringing ned seamless data into grass, >> then reprojecting to utm, i would be grateful. or at this point, just how >> to use r.proj to reproject to utm. >> thanx, janet >> >> On Tue, Dec 1, 2009 at 10:14 PM, maning sambale >> wrote: >>> >>> Try to include the full path of your location >>> ([/work/projects/Willamette]) >>> >>> On Wed, Dec 2, 2009 at 1:47 PM, Janet Choate wrote: >>> > HI, >>> > thanx for any help. seems this should be a pretty straight forward >>> > thing to >>> > do, but am not having any success. >>> > i am trying to reproject a raster from lat/long to utm. >>> > i created a new location with a utm projection, then tried to use r.proj >>> > to >>> > reproject a raster from another location. >>> > >>> > r.proj input=dem location=/work/projects/Willamette mapset=PERMANENT >>> > >>> > however, i get the following error: >>> > Mapset [PERMANENT] in input location [/work/projects/Willamette] - not >>> > found >>> > >>> > thanx for any assistance, >>> > Janet >>> > >>> > ___ >>> > grass-user mailing list >>> > grass-user@lists.osgeo.org >>> > http://lists.osgeo.org/mailman/listinfo/grass-user >>> > >>> > >>> >>> >>> >>> -- >>> 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] reproject from lat/long to utm
On Wed, Dec 2, 2009 at 2:28 PM, Janet Choate wrote: > I tried using the full path (i.e. the PERMANENT mapset lives within the > /work/projects/Willamette location), but get the 'not found' error. > is there any other command i need to run prior to r.proj, or any additional > options to include, or how to list the path other than the way i am? just to check from your lat/long mapset, run g.gisenv It should show you the location, mapset and GRASS database. This is what you use in r.proj parameters > the original data i imported was a dem from usgs national elevation dataset > (ned) seamless data distribution system (sdds). i had to export it from ned > as a geotiff to bring it into grass. i tried exporting the data from ned as > an arcgrid and importing into grass with r.in.gdal, but was unable to bring > it into grass. > > so if anyone can tell me the steps to bringing ned seamless data into grass, > then reprojecting to utm, i would be grateful. or at this point, just how > to use r.proj to reproject to utm. > thanx, janet > > On Tue, Dec 1, 2009 at 10:14 PM, maning sambale > wrote: >> >> Try to include the full path of your location >> ([/work/projects/Willamette]) >> >> On Wed, Dec 2, 2009 at 1:47 PM, Janet Choate wrote: >> > HI, >> > thanx for any help. seems this should be a pretty straight forward >> > thing to >> > do, but am not having any success. >> > i am trying to reproject a raster from lat/long to utm. >> > i created a new location with a utm projection, then tried to use r.proj >> > to >> > reproject a raster from another location. >> > >> > r.proj input=dem location=/work/projects/Willamette mapset=PERMANENT >> > >> > however, i get the following error: >> > Mapset [PERMANENT] in input location [/work/projects/Willamette] - not >> > found >> > >> > thanx for any assistance, >> > Janet >> > >> > ___ >> > grass-user mailing list >> > grass-user@lists.osgeo.org >> > http://lists.osgeo.org/mailman/listinfo/grass-user >> > >> > >> >> >> >> -- >> 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] reproject from lat/long to utm
I tried using the full path (i.e. the PERMANENT mapset lives within the /work/projects/Willamette location), but get the 'not found' error. is there any other command i need to run prior to r.proj, or any additional options to include, or how to list the path other than the way i am? the original data i imported was a dem from usgs national elevation dataset (ned) seamless data distribution system (sdds). i had to export it from ned as a geotiff to bring it into grass. i tried exporting the data from ned as an arcgrid and importing into grass with r.in.gdal, but was unable to bring it into grass. so if anyone can tell me the steps to bringing ned seamless data into grass, then reprojecting to utm, i would be grateful. or at this point, just how to use r.proj to reproject to utm. thanx, janet On Tue, Dec 1, 2009 at 10:14 PM, maning sambale wrote: > Try to include the full path of your location ([/work/projects/Willamette]) > > On Wed, Dec 2, 2009 at 1:47 PM, Janet Choate wrote: > > HI, > > thanx for any help. seems this should be a pretty straight forward thing > to > > do, but am not having any success. > > i am trying to reproject a raster from lat/long to utm. > > i created a new location with a utm projection, then tried to use r.proj > to > > reproject a raster from another location. > > > > r.proj input=dem location=/work/projects/Willamette mapset=PERMANENT > > > > however, i get the following error: > > Mapset [PERMANENT] in input location [/work/projects/Willamette] - not > found > > > > thanx for any assistance, > > Janet > > > > ___ > > grass-user mailing list > > grass-user@lists.osgeo.org > > http://lists.osgeo.org/mailman/listinfo/grass-user > > > > > > > > -- > 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] reproject from lat/long to utm
Try to include the full path of your location ([/work/projects/Willamette]) On Wed, Dec 2, 2009 at 1:47 PM, Janet Choate wrote: > HI, > thanx for any help. seems this should be a pretty straight forward thing to > do, but am not having any success. > i am trying to reproject a raster from lat/long to utm. > i created a new location with a utm projection, then tried to use r.proj to > reproject a raster from another location. > > r.proj input=dem location=/work/projects/Willamette mapset=PERMANENT > > however, i get the following error: > Mapset [PERMANENT] in input location [/work/projects/Willamette] - not found > > thanx for any assistance, > Janet > > ___ > grass-user mailing list > grass-user@lists.osgeo.org > http://lists.osgeo.org/mailman/listinfo/grass-user > > -- 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
[GRASS-user] reproject from lat/long to utm
HI, thanx for any help. seems this should be a pretty straight forward thing to do, but am not having any success. i am trying to reproject a raster from lat/long to utm. i created a new location with a utm projection, then tried to use r.proj to reproject a raster from another location. r.proj input=dem location=/work/projects/Willamette mapset=PERMANENT however, i get the following error: Mapset [PERMANENT] in input location [/work/projects/Willamette] - not found thanx for any assistance, Janet ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: Georectify problem
Ross Benisch hi Ross, please open a bug report at https://trac.osgeo.org/grass/wiki by specifying Windows-Installer, steps to reproduce the error, Windows-version etc. so important information can't get lost in the mailing list best regards Helmut ___ Pop-Diva Sarah Kreuz und ihr beeindruckendes Debütalbum "One Moment in Time". Jetzt reinhören! http://produkte.web.de/go/05/ -- next part -- A non-text attachment was scrubbed... Name: smime.p7s Type: application/pkcs7-signature Size: 1747 bytes Desc: S/MIME Cryptographic Signature Url : http://lists.osgeo.org/pipermail/grass-user/attachments/20091201/0ed6f741/smime-0001.bin -- Message: 4 Date: Tue, 01 Dec 2009 01:45:50 +0100 From: Helmut Kudrnovsky Subject: [GRASS-user] Trouble setting location and projection To: grass-user@lists.osgeo.org Message-ID: <9735556.17178.1259628351250.javamail.fm...@fmcert01.dlan.cinetic.de> Content-Type: text/plain; charset="iso-8859-15" I am trying to set the project location and cordinate system using the 6.4.0svn version for windows but I keep getting the projection error, 'g.proj.exe-p' when I try and start Grass. I have a GIS background and I am frustrated that I can't get it to work. I want to use the projection nad83UTM >Zone >14N. I have read through all of the tuturiols and help sheets and looked at the FAQs with no such luck. I really want to use this program to help with >my job so >if you would help me get started I would greatly appreciate it. Sincerely, Ross Benisch hi Ross, please open a bug report at https://trac.osgeo.org/grass/wiki by specifying Windows-Installer, steps to reproduce the error, Windows-version etc. so important information can't get lost in the mailing list best regards Helmut maybe related to following bug report: http://trac.osgeo.org/grass/ticket/827 please add your comments best regards Helmut ___ Pop-Diva Sarah Kreuz und ihr beeindruckendes Debütalbum "One Moment in Time". Jetzt reinhören! http://produkte.web.de/go/05/ -- next part -- A non-text attachment was scrubbed... Name: smime.p7s Type: application/pkcs7-signature Size: 1747 bytes Desc: S/MIME Cryptographic Signature Url : http://lists.osgeo.org/pipermail/grass-user/attachments/20091201/c3375ab2/smime-0001.bin -- Message: 5 Date: Tue, 01 Dec 2009 01:48:12 +0100 From: Marco Tuckner Subject: [GRASS-user] i.atcorr with MSS band 5,6,7 To: grass-user@lists.osgeo.org Message-ID: Content-Type: text/plain; charset=ISO-8859-15 Dear Grass users, how can I add spectral conditions of the bands 5,6,7 of the Landsat MSS sensor? I have seen that these are defined in: http://trac.osgeo.org/grass/browser/grass/branches/develbranch_6/imagery/i.atcorr/Iwave.cpp But is is unclear to me how they are added. As far as the manual concerned, the pre-defined spectral conditions are only available to LS MSS band 1-4: http://grass.itc.it/grass64/manuals/html64_user/i.atcorr.html 31 mss (landsat5) band 1 (0.475-0.640) 32 mss (landsat5) band 2 (0.580-0.750) 33 mss (landsat5) band 3 (0.655-0.855) 34 mss (landsat5) band 4 (0.785-1.100) Thanks in advance for any hints, Marco -- Message: 6 Date: Mon, 30 Nov 2009 21:14:48 -0800 From: mh...@berkeley.edu Subject: Re: [GRASS-user] i.atcorr with MSS band 5,6,7 To: "Marco Tuckner" Cc: grass-user@lists.osgeo.org Message-ID: <65b87b606fc185067cf55224e63558f9.squir...@calmail.berkeley.edu> Content-Type: text/plain;charset=utf-8 I thought the MSS sensor only had 4 bands, and the TMS sensor sensor 7? Best, MEH Dear Grass users, how can I add spectral conditions of the bands 5,6,7 of the Landsat MSS sensor? I have seen that these are defined in: http://trac.osgeo.org/grass/browser/grass/branches/develbranch_6/imagery/i.atcorr/Iwave.cpp But is is unclear to me how they are added. As far as the manual concerned, the pre-defined spectral conditions are only available to LS MSS band 1-4: http://grass.itc.it/grass64/manuals/html64_user/i.atcorr.html 31 mss (landsat5) band 1 (0.475-0.640) 32 mss (landsat5) band 2 (0.580-0.750) 33 mss (landsat5) band 3 (0.655-0.855) 34 mss (landsat5) band 4 (0.785-1.100) Thanks in advance for any hints, Marco ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user -- Message: 7 Date: Tue, 01 Dec 2009 01:32:53 -0500 From: Peng Du Subject: [GRASS-user] reading a column/vector of a raster image To: grass-us
[GRASS-user] memory leak in calc_covariance
Hi There In grass6_devel/imagery/i.pca/main.c, calc_covariance(...), rowbuf1 and rowbuf2 are each allocated bands times, but only freed bands-1 times. Isn't this a leak? Thanks, Peng ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Re: i.atcorr with MSS band 5,6,7
i.atcorr supports the MSS sensor aboard Landsat 5 and is not applicable to Landsat 1 data (nor landsat's 2,3 or 4). In-order to atmospherically correct landsat 1 data with i.atcorr you'll need to add support for it. It's not hard and I would be willing to spare what free time I can to help you but the first step, and possibly the most difficult, is to find the spectral response curves for the RBV and MSS sensors aboard Landsat 1. Cheers, Mike On 1-Dec-09, at 3:29 PM, Marco Tuckner wrote: I thought the MSS sensor only had 4 bands, and the TMS sensor sensor 7? Sorry, I heven't expressed myself clearly. I am working with old images from the 1970ies recorded by Landsat 1. According to http://landsat.gsfc.nasa.gov/about/mss.html http://eros.usgs.gov/#/Guides/landsat_mss the MSS used the bands 4-8. I am mislead somehow or wrong? If not I would be glad to receive some advice on how to add these additional spectral characteristsics to the module. Best regards, Marco ___ 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: i.atcorr with MSS band 5,6,7
> I thought the MSS sensor only had 4 bands, and the TMS sensor sensor 7? Sorry, I heven't expressed myself clearly. I am working with old images from the 1970ies recorded by Landsat 1. According to http://landsat.gsfc.nasa.gov/about/mss.html http://eros.usgs.gov/#/Guides/landsat_mss the MSS used the bands 4-8. I am mislead somehow or wrong? If not I would be glad to receive some advice on how to add these additional spectral characteristsics to the module. Best regards, Marco ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Finding Missing Point
I have a map that's supposed to have 3 points (sites) on it, all relatively close together. The geographic coordinates were projected to lcc using v.proj and the region's boundaries and parameters adjusted to the map using 'g.region vect=sites'. My problem is that only two sites display, and the relative x values are switched (that is, the one on the left should be on the right). How should I go about determining what went wrong? What information is needed? Rich ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] r.horizon problems with large DEM
Joshua, On Mon, Nov 30, 2009 at 5:50 PM, Joshua Campbell wrote: > Markus, > > Thanks for the input. I reviewed the r.sun page and unfortunately r.sun > doesn't provide the specific horizon angle provided by r.horizon. My goal is > not to compute solar irradiance. ok. > I was able to run r.horizon on a subset that was ~4900x9900. The DEM I'm > trying to use in 3x78000 (but is has quite a bit of null values). > > Do you know if there is pixel dimension limit to r.horizon? I suspect that you have hit the 2GB limit but... > I am using the 64-bit package available through the Ubuntu site -- would it > make a difference if I built the code from source? ... if it is *really* a 64 bit binary that should not matter. If 32bit, then there is the --enable-largefile flag to configure the GRASS source code before compilation. > Do anyone know of a script that would iteratively set a region, compute > r.horizon (with a buffer), move the region, and r.patch the resulting > horizon files? > If anyone is interested in writing this script, we could discuss > compensation (contact directly at jcampbell.geo at gmail.com) Sounds good :) Unfortunately I have no time... best Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] *Direction* to nearest vector?
Just a guess... Maybe v.distance? Create a point raster, on point for each cell in your raster. Run v.distance and look for the to_angle value calculated. Then convert the points back to raster v.distance man page... http://www.ces.iisc.ernet.in/grass/grass64/manuals/html64_user/v.distance.html Or you could try something with the r.distance module and the east1,north1 values Good luck daniel On Mon, Nov 30, 2009 at 4:44 PM, Jonathan Greenberg wrote: > GRASSers: > > Is there any way to generate a raster in which the pixel values are the > directions (in radians or degrees) to the nearest vector object, e.g. a > point, line, or polygon ("nearest" being calculated as Euclidean distance)? > --j > > -- > > Jonathan A. Greenberg, PhD > Postdoctoral Scholar > Center for Spatial Technologies and Remote Sensing (CSTARS) > University of California, Davis > One Shields Avenue > The Barn, Room 250N > Davis, CA 95616 > Phone: 415-763-5476 > AIM: jgrn307, MSN: jgrn...@hotmail.com, Gchat: jgrn307 > ___ > 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] reading a column/vector of a raster image
Peng Du wrote: > I am dealing with a satellite image which has N bands, each band of size > row x col. Is there any way to read one column of a certain band or all > the bands of a certain point? I only found G_get_raster_row which works > for row. You can use G_set_window() (Rast_set_window() in 7.0) to control which samples are returned for each row. E.g. for column "n": struct Cell_head cellhd; ... G_get_set_window(&cellhd); cellhd.west += cellhd.ew_res * n; cellhd.east = cellhd.west + cellhd.ew_res; G_set_window(&cellhd); Then, each call to G_get_raster_row() etc will store a single value in the buffer. However, reducing the number of columns won't have any noticeable effect upon speed, as the library reads and decompresses entire rows. If you want to process maps in column-major order, you should create a transposed copy in a temporary file or in memory. -- Glynn Clements ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] using GDAL to view netcdf files in GRASS
Carbonari, Katie (IS) wrote: Hello. I'm trying to use GDAL to convert my netcdf file into a GTiff so it can be read into GRASS. I tried using "gdal_translate -of GTiff -b 1 NETCDF:FILE.nc:Outfile Outfile.tiff" You probably need to specify a coordinate system, too. Do this with the -a_srs option. See http://www.gdal.org/gdal_utilities.html for details. The easiest way is using the epsg codes. On my system, I can find them by: cat /usr/share/proj/epsg | grep -i lambert -A 1 There is also a list at (http://www.epsg-registry.org/) - look under Projected CRS type. Choose the one that you need, and amend your code. e.g. gdal_translate -of GTiff -a_srs epsg:12345 -b 1 NETCDF:FILE.nc:Outfile Outfile.tiff You might also need to warp your map into the new projection gdalwarp -t_srs epsg 12345 Outfile.tiff Outfile_warped.tiff See the first parts of this tutorial for details http://blog.thematicmapping.org/2008/03/generating-map-tiles-with-gdal2tiles.html Hope that helps, John -- Dr John Stevenson Postdoctoral Research Associate School of Earth, Atmospheric and Environmental Sciences Williamson Building (Room 2.42) University of Manchester Manchester M13 9PL, UK tel. +44(0)161 306 6585; fax. +44(0)161 306 9361; john.steven...@manchester.ac.uk ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Execution failed g.proj.exe
>Dear List, > >I just installed WinGRASS-6.4.0svn-r39740-2-Setup.exe on Windows XP >professional and when starting the system using the >spearfish database I get >the message Error in command execution: Execution failed g.proj.exe -p . > >Any ideas? > >Bes regards, > >Klaus please add your comments to following bug report: http://trac.osgeo.org/grass/ticket/827 by specifying the situation when the error comes up. best regards Helmut ___ Preisknaller: WEB.DE DSL Flatrate für nur 16,99 Euro/mtl.! http://produkte.web.de/go/02/ smime.p7s Description: S/MIME Cryptographic Signature ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user