Re: [GRASS-user] Re: Is it possible to export a Raster to KML

2011-10-28 Thread Jamie Adams
Roger is correct. GE will just drape your raster data over it's terrain mesh regardless of what vertical parameters your data has. It is essentially 2D with a z value and is interpreted as such. That is why EPSG:4326 works, because GE doesn't care about any other information. No

[GRASS-user] wiki link to addon r.bilateral broken

2011-07-18 Thread Jamie Adams
Hello all, I'd like to install the addon r.bilateral from the add-ons page ( http://les-ejk.cz/files/programs/grass/r.bilateral.tgz) but the link is broken. Does anyone have the source tgz they could share with me? Thanks, Jamie ___ grass-user

Re: [GRASS-user] Processing Atlas Lidar Data - Problem(s)

2011-03-10 Thread Jamie Adams
Hi Tim, When dealing with lidar data, I always import using r.in.xyz first. You can easily vary the resolution to generate a quick preview and it will give you a sense of the data source and what information it contains. If you then find you need to do more advanced processing, you can move

Re: [GRASS-user] How to merge rasters using feathering/blending ?

2011-02-28 Thread Jamie Adams
I wrote a script to do something like this at one point, though I can't find it now. It was intended for terrain, but could be used per band on imagery. It used r.buffer to create a distance gradient and then r.recode to convert this gradient to a percentage. I then used r.mapcalc to merge the

Re: [GRASS-user] Opensource and pretty maps design anybody?

2010-08-17 Thread Jamie Adams
I generally open GRASS data in QGIS, create a basic map and export to svg (or pdf), and then do the finishing touches in Inkscape. This workflow has worked very well for me. - Jamie On Tue, Aug 17, 2010 at 8:56 AM, kapo coulibaly kmcou...@gmail.com wrote: I've been using opensource GIS for a

Re: [GRASS-user] dems from coordinate lists

2010-05-13 Thread Jamie Adams
You want the points to represent cell centers so you need to expand your region 1/2 pixel in all four directions. Also look at using r.in.xyz, it works directly on this type of data. Jamie On Thu, May 13, 2010 at 7:26 AM, Hanlie Pretorius hanlie.pretor...@gmail.com wrote: Hi, I've obtained

[GRASS-user] Re: dems from coordinate lists

2010-05-13 Thread Jamie Adams
Did you adjust the region also? To do this properly, you need to first scan the file using r.in.xyz to get the extents. I'd also use the shell style output. r.in.xyz -sg input=2628cc.ORT.xyz output=dem_2628cc_25m_xyz Should return something like: n=1000 s=500 w=500 e=1000 t=### b=### (you can

Re: [GRASS-user] can't ingest an E00 file using r.in.gdal, or anything else?

2009-10-17 Thread Jamie Adams
Hey Jeff, This is a valid file, though not well supported by software other than ESRI. E00 files most commonly hold vector formats, but can be used to contain raster data also. AFAIK, there is no way to import these directly into GRASS, it only supports the vector version with v.in.e00 For this

Re: [GRASS-user] nnbathy

2009-09-04 Thread Jamie Adams
Reviving this old thread - I just noticed the other day that the nn code has been checked into Google Code using an MIT license. Would it be possible to integrate directly now? I find r.surf.nnbathy to be one of the more useful interpolation routines. - Jamie On Tue, Sep 23, 2008 at 8:50 AM,

Re: [GRASS-user] mosaic raster files

2009-07-02 Thread Jamie Adams
If you have GDAL 1.6 installed, you can use gdalbuildvrt to create a VRT of the files, and then link to this using r.external. This way you have a mosaic without duplicating the data. On Thu, Jul 2, 2009 at 7:53 AM, Milton Cezar Ribeiro miltinho.astrona...@gmail.com wrote: Hi Guys, I will

Re: [GRASS-user] python commands on grass.

2009-05-05 Thread Jamie Adams
Try using help() from within the python shell. python import grass help(grass) -Jamie On Tue, May 5, 2009 at 11:52 AM, Milton Cezar Ribeiro miltinho.astrona...@gmail.com wrote: Dear all, is there a way of I know the commands available on python grass library? Checking the codes available

Re: [GRASS-user] problem in projection of shape files

2009-02-23 Thread Jamie Adams
I would use ogr2ogr http://www.gdal.org/ogr2ogr.html from GDAL to reproject using a proj4 string. For example to wgs84: ogr2ogr -s_srs +proj=tmerc +lat_0=0 +lon_0=90 +k=0.9996 +x_0=50 +y_0=-200 +a=6377276.345 +b=6356075.41314024 +units=m +no_defs -t_srs wgs84 out.shp in.shp -Jamie On

Re: [GRASS-user] saving patch areas on raster maps

2008-11-14 Thread Jamie Adams
You can use r.recode to do this, but you'll need to reformat the output from r.stats. Example: *r.stats -a -n raster* 0 1234 1 2345 2 3456 .. This needs to be: 0:0:1234:1234 1:1:2345:2345 2:2:3456:3456 .. Save this to a txt file and use it with the rules flag of r.recode. If you

Re: [GRASS-user] Re: r.resamp.stats max method not working as expected

2008-10-16 Thread Jamie Adams
I think 7 is fine. My current approach is working well enough. Thanks! On Thu, Oct 16, 2008 at 2:34 AM, Markus Neteler [EMAIL PROTECTED] wrote: On Wed, Oct 15, 2008 at 5:39 PM, Glynn Clements [EMAIL PROTECTED] wrote: Glynn Clements wrote: ... Implemented in r33888 (in 7.0). The

[GRASS-user] Re: r.resamp.stats max method not working as expected

2008-10-14 Thread Jamie Adams
at 12:54 PM, Jamie Adams [EMAIL PROTECTED] wrote: Hello all, When I use r.resamp.stats to resample a grid of 0s 1s using the maximum option and the -w flag, it doesn't seem to consider the edge pixels that overlap. It returns the same result as without the -w flag. I've gotten around

Re: [GRASS-user] Re: r.resamp.stats max method not working as expected

2008-10-14 Thread Jamie Adams
Very cool. I went ahead and filed a enhancement ticket for this in trac: http://trac.osgeo.org/grass/ticket/337 On Tue, Oct 14, 2008 at 2:34 AM, Glynn Clements [EMAIL PROTECTED]wrote: Jamie Adams wrote: Okay, I did find verification that this doesn't work in r.resamp.stats: *[1

Re: [GRASS-user] r.in.gdal global raster incorrect extents

2008-10-14 Thread Jamie Adams
but why then does it work fine for me? I downloaded a 180W tile, that looks fine too: ftp://e0srp01u.ecs.nasa.gov/srtm/version2/SRTM3/Australia/S20W180.hgt.zip Origin = (-180.0004166,-18.9995834) Pixel Size = (0.0008333,-0.0008333) I grabbed the same

[GRASS-user] r.in.gdal global raster incorrect extents

2008-10-13 Thread Jamie Adams
Hello all, When I import a global elevation dataset with extents 180w 90n 180e 90s using r.in.gdal, the GRASS raster shows up with w-e extents both showing 180E. While this doesn't seem to be affecting raster calculations, afaik, it's causing issues with QGIS. It displays it with it's origin at

Re: [GRASS-user] r.in.gdal global raster incorrect extents

2008-10-13 Thread Jamie Adams
at 1:35 PM, Nikos Alexandris [EMAIL PROTECTED] wrote: On Mon, 2008-10-13 at 12:41 -0700, Jamie Adams wrote: Hello all, When I import a global elevation dataset with extents 180w 90n 180e 90s using r.in.gdal, the GRASS raster shows up with w-e extents both showing 180E. While

Re: [GRASS-user] r.in.gdal global raster incorrect extents

2008-10-13 Thread Jamie Adams
13, 2008 at 4:18 PM, Nikos Alexandris [EMAIL PROTECTED] wrote: On Mon, 2008-10-13 at 16:13 -0700, Jamie Adams wrote: Thanks for the suggestions. I opened the image up in python using gdal, looked at the geotransform and found that the image is ever-so-slightly shifted west

Re: [GRASS-user] r.in.gdal global raster incorrect extents

2008-10-13 Thread Jamie Adams
On Mon, Oct 13, 2008 at 6:08 PM, Hamish [EMAIL PROTECTED] wrote: Jamie Adams wrote: In my case, the source GeoTiff had a geotransform of: (-180.012497, 0.0083297, 0.0, 90.0, 0.0, -0.0083297) One strange thing was, the original grass import (west

[GRASS-user] get modified date for raster layer

2008-09-02 Thread Jamie Adams
, but that essentially greps it out of a full list. Is there a more efficient way of getting the modified date of a raster layer? If not, getting the full path instead? Thanks, Jamie Adams ___ grass-user mailing list grass-user@lists.osgeo.org http

Re: [GRASS-user] get modified date for raster layer

2008-09-02 Thread Jamie Adams
On Tue, Sep 2, 2008 at 6:42 PM, Hamish [EMAIL PROTECTED] wrote: Jamie Adams [EMAIL PROTECTED] I've got over 14k rasters in several mapsets that I need to store modified date in a db for tracking purposes. I've got the storage and retrieval worked out, but I don't see an clean (easy) way

Re: [GRASS-user] Re: make a raster index in GRASS

2008-07-16 Thread Jamie Adams
to the cellhd files. 2008/7/16 Ivan Shmakov [EMAIL PROTECTED]: Jamie Adams [EMAIL PROTECTED] writes: Hello all, I'd like to generate a raster index polygon file, just like gdaltindex does, but using GRASS rasters as input. Is there a command I'm overlooking? It may

[GRASS-user] make a raster index in GRASS

2008-07-15 Thread Jamie Adams
Hello all, I'd like to generate a raster index polygon file, just like gdaltindex does, but using GRASS rasters as input. Is there a command I'm overlooking? If not, any ideas on how to do this? -Jamie ___ grass-user mailing list

[GRASS-user] SRTM and r.fillnulls - using neighboring tiles

2008-07-03 Thread Jamie Adams
Hello all, I'm trying to devise an efficient way to fill voids in SRTM, using neighboring tiles when there is edge nodata. In the end, I want to have tiles matching the input, just lacking any voids. My initial thought is to write a simple gdal python script to determine which tiles have edge

Re: [GRASS-user] Orthorectification of Historic Aerial Photographs

2008-04-10 Thread Jamie Adams
Would it be possible to calculate these parameters with a panoramic creation package like panotools? http://wiki.panotools.org/Lens_correction_model It calculates a value for field of view, but I didn't notice an actual focal length parameter. Maybe this is easy to associate? One benefit is