Re: [GRASS-user] Attribute transfer from one table to another based on position of vector items

2012-06-09 Thread Michael Perdue
I have a similar goal, but the solution in the case below will not work for me. I have three points A,B C that I have linked together through Delaunay triangulation; A /\ / \ /\ / \ /\ / \ /\ B C I want to

Re: [GRASS-user] Attribute transfer from one table to another based on position of vector items

2012-06-09 Thread Michael Perdue
Lee, Nepomuk Reinhard und Nils Räder Softwarecenter 3 35037 Marburg Festnetz: +49 6421 379 6256 Mobil: +49 176 6127 7269 E-Mail: l...@isi-solutions.org Web: http://www.isi-solutions.org 2012/6/9 Michael Perdue michael_per...@yahoo.ca I have a similar goal, but the solution

Re: [GRASS-user] Determining Datum From Projection Coordinates [SOLVED]

2011-03-02 Thread Michael Perdue
Hate to be a wet blanket, but both NAD27 and NAD83 UTM coordinates are supposed to be in meters. The False Easting for the UTM zones are defined to be 50m regardless of datum. So, the presence of meters in your DB is not a solid indication that it is in NAD83. As the lat/long in you're file

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

2011-02-28 Thread Michael Perdue
In GRASS, one thing you can try is to normalize your data so the variations from scene to scene are minimized. ie, run a linear regression between overlapping maps and then use the coefficients to normalize one of the data sets. However, it is tedious and doesn't work well if the two images

Re: [GRASS-user] Re: i.atcorr with MSS band 5,6,7

2009-12-01 Thread Michael Perdue
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

Re: [GRASS-user] i.atcorr with ASTER L1B

2009-07-29 Thread Michael Perdue
. Will it be ok to post a txt file on the list with the interpretation results? Cheers Andrew On 7/29/09, Michael Perdue michael_per...@yahoo.ca wrote: Hi Andy, Yes, it looks like the right info for the SPOT sensor platforms. Iwave.cpp expects the data to be in 2.5nm increments so you will have

Re: [GRASS-user] i.atcorr with ASTER L1B

2009-07-28 Thread Michael Perdue
, Michael Perdue michael_per...@yahoo.ca wrote: Hi Markus, Here is a text file with everything organized as I think it should be. I didn't add the this Iwave.cpp to submit a patch as I wasn't really sure how to modify GeomCond.cpp. Cheers, Mike On 26-Jul-09, at 1:34 PM, Markus Neteler wrote: Great

Re: [GRASS-user] i.atcorr with ASTER L1B

2009-07-27 Thread Michael Perdue
Hi Markus, Here is a text file with everything organized as I think it should be. I didn't add the this Iwave.cpp to submit a patch as I wasn't really sure how to modify GeomCond.cpp. Cheers, Mike On 26-Jul-09, at 1:34 PM, Markus Neteler wrote: Great - the header of Iwave.cpp contains

Re: [GRASS-user] i.atcorr with ASTER L1B

2009-07-26 Thread Michael Perdue
I believe these are what you're looking for. http://www.science.aster.ersdac.or.jp/en/about_aster/sensor/excel/SPEC_ROFILES.LZH Cheers, Mike On 26-Jul-09, at 10:56 AM, Markus Neteler wrote: Hi Stefan, On Sun, Jul 26, 2009 at 1:27 PM, Stefan Syllastefansy...@gmx.de wrote: Hi Markus,

Re: [GRASS-user] i.atcorr with ASTER L1B

2009-07-26 Thread Michael Perdue
Organized like this? /* 1st spectral band of enhanced thematic mapper plus */ static const float sr1[35] = { .0160, .0710, .2870, .6665, .7920, .8570, .8390, .8060, .7790, .8460, .9010, .8995, .8900, .8510, .8750, .8930, .8840, .9300, .9580, .9540, .9800, .9750, .9650,

Re: [GRASS-user] i.atcorr with ASTER L1B

2009-07-26 Thread Michael Perdue
I'm guessing that the values in these two assignments represent the left and right edges of the Spectral response curves. Is that correct? static const float wli[7] = { 0.4350, 0.5000, 0.6150, 0.7400, 1.5050, 2.0150, 0.5025 }; static const float wls[7] = { 0.5200, 0.6225, 0.7025, 0.9125,

Re: [GRASS-user] i.atcorr - help, please

2009-07-22 Thread Michael Perdue
Some one beat me to it or I must be blind. The relevant bits are in the wiki. Cheers, Mike On 21-Jul-09, at 2:49 PM, Markus Neteler wrote: On Tue, Jul 21, 2009 at 8:06 PM, Michael Perduemichael_per...@yahoo.ca wrote: I thought I'd read in a previous e-mail that the module ran much

[GRASS-user] CLI syntax for gain in i.landsat.toar

2009-07-22 Thread Michael Perdue
Hi, I've been playing around with i.landsat.toar and am having problems correctly entering the gain values on the command line. The bands are labeled as: l7_049026_01_toaDN l7_049026_04_toaDN l7_049026_08_toaDN l7_049026_02_toaDN l7_049026_05_toaDN l7_049026_6h_toaDN l7_049026_03_toaDN

Re: [GRASS-user] i.atcorr - help, please

2009-07-21 Thread Michael Perdue
I thought I'd read in a previous e-mail that the module ran much faster if the user used an intergerized DEM rather than a floating point DEM. I can't seem to find the original post but if true, I think it would be a very useful piece of information for the wiki. On 21-Jul-09, at 3:29 AM,

Re: [GRASS-user] AVIRIS data in GRASS

2009-07-06 Thread Michael Perdue
I was able to generate a ENVI style header that will allow you to read the files with GDAL. Modify as necessary and you should be able to work with the *.img files. On 5-Jul-09, at 6:56 PM, Michael Perdue wrote: ... you might be able to manipulate GDAL to read it. If you can create

Re: [GRASS-user] AVIRIS data in GRASS

2009-07-05 Thread Michael Perdue
I don't think you can use r.in.bin because there is no option to specify the number of bands and how those bands are interleaved. In the past I've used Multispec and Matlab to read the data but you might be able to manipulate GDAL to read it. If you can create a header file for a generic

[GRASS-user] Testing two vectors for equality

2009-06-26 Thread Michael Perdue
Is there a way to test if two vectors are identical, including their attribute tables? i.e., I've changed some code and now I want to test that the results are the same. Is there an easy way to do this? Cheers, Mike ___ grass-user mailing list

Re: [GRASS-user] v.lidar.edgedetection very slow

2009-06-22 Thread Michael Perdue
I read through the one of the papers today (M. A. Brovelli , M. Cannata and U. M. Longoni, 2002. Managing and processing LIDAR data within GRASS, Proceedings of the Open source GIS - GRASS users conference 2002) and the authors recommended a step value or 3 to 4 time the average point

Re: [GRASS-user] v.lidar.edgedetection very slow

2009-06-21 Thread Michael Perdue
On 19-Jun-09, at 8:54 PM, Hamish wrote: Michael wrote: I actually get the same thing. For the first while the process runs in full use of one core, but once the process starts writing results to the db, the whole process bogs down to a grinding halt. CPU usage for v.lidar.edgedetection drops

Re: [GRASS-user] v.lidar.edgedetection very slow

2009-06-19 Thread Michael Perdue
I actually get the same thing. For the first whiletheprocess runs in full use of one core, but once the process starts writing results to the db, the wholeprocessbogs down to a grindinghalt. CPU usage for v.lidar.edgedetection drops down to ~1% (1% of one core) and sqlite usage rises to ~16%.

Re: [GRASS-user] v.delaunay Bug in Mac OS X

2009-06-17 Thread Michael Perdue
OOPS, forgot to hit the list with this one. On 17-Jun-09, at 2:30 AM, Hamish wrote: William wrote: I don't see any changes in v.delaunay itself from RC4, though maybe there are changes in the vector libraries. Did this work in RC4? Note since Nov 26 2008 (ie all 6.4.0RCs) v.delaunay has

[GRASS-user] v.delaunay Bug in Mac OS X

2009-06-16 Thread Michael Perdue
Hi, v.delaunay is generating incorrect results on Mac OS X on my machine (Intel Mac Pro, Leopard, W. Kyngesburye's GRASS 6.4-RC5 binaries). The following commands give the results below; echo 456490.2375|5905949.8625|1 456690.0917|5905895.5625|2 456570.9333|5905821.65416667|3

Re: [GRASS-user] v.delaunay Bug in Mac OS X

2009-06-16 Thread Michael Perdue
? A big difference in RC5 for OSX is that it's now 64bit. Does this work: arch -i386 v.delaunay in=DelaunayPoints out=DelaunayTriangles On Jun 16, 2009, at 5:50 PM, Michael Perdue wrote: Hi, v.delaunay is generating incorrect results on Mac OS X on my machine (Intel Mac Pro, Leopard, W

Re: [GRASS-user] patching images from different UTM zones (without shell script?)

2009-04-01 Thread Michael Perdue
You can try using the -n flag for r.proj to avoid region clipping or you can extend your region beyond the edges of the zone. On 1-Apr-09, at 4:10 PM, Renae Mackas wrote: Hello GRASS users- I am trying to create a mosaic of LANDSAT taken of Eastern Hudson Bay area using GRASS 6.3.0. I

Re: [GRASS-user] v.in.ascii

2009-02-03 Thread Michael Perdue
If you create a VRT file for your CSV file, then you should be able to import it with v.in.ogr directly. v.in.ogr dsn=/Users/Desktop/Data/ 081211_Mitchell_plane_bores_list_V002.vrt out=mitchell_bores You VRT file would look like: OGRVRTDataSource OGRVRTLayer

Re: [GRASS-user] Centerlines of Polygons (skeletons)

2008-06-18 Thread Michael Perdue
You could do; v.to.rast - r.thin - r.to.vect - v.clean (and maybe v.generalize) But you'll lose the link to the attribute data. Cheers, Mike - Original Message From: Aurora Geomatics [EMAIL PROTECTED] To: [EMAIL PROTECTED] Sent: Wednesday, June 18, 2008 5:24:46 PM Subject:

Re: [GRASS-user] Problem reprojecting to UTM 34 S

2008-05-15 Thread Michael Perdue
From what I see, gdalinfo hak_ser2.asc is listing all the latitudes as being positive which means that they are in the northern hemisphere. It looks to me like gdal is doing the right thing and there is something wrong with the ASCII grid. Mike On 15-May-08, at 8:01 AM, Rainer M Krug

Re: [GRASS-user] r.in.xyz

2008-01-18 Thread Michael Perdue
You can use sed to find and replace patterns within a file. So; ~ cat temp.txt 1619923.0297213|4836235.79480289|0.01571819 1619927.57223432|4836235.09835698|nan 1619910.45256599|4836241.22590928|6.82026114 To make a substitution ~ sed 's/nan/-/g' ./temp.txt

Re: [GRASS-user] v.out.ascii -- surpressing vector IDs?

2008-01-18 Thread Michael Perdue
Not exactly what your asking, but you could pipe the output through awk to remove any fields that you don't want. Cheers, Mike On 18-Jan-08, at 3:23 PM, Jonathan Greenberg wrote: I'm trying to export a 3d vector to an ascii file that just contains x y z using v.out.ascii, but the output

Re: [GRASS-user] OFF-TOPIC - Earth Curvature in profiles

2008-01-12 Thread Michael Perdue
If I understand you correctly, you want to draw a topographic profile superimposed on an arc representing the curvature of the earth? What you want to do then is show the orthogonal distance of each sample from the chord defined by the start and end point of profile. You can try as follows