Message-ID: <d8ad40b50907130904w51d71338me80fa97f09aa0cc1@mail.gmail.com>
Date: 2009-07-13T16:04:40Z
From: Barry Rowlingson
Subject: kml generation question
In-Reply-To: <6d4aad460907130832t218c0ee7ne6f9af0492e76ce0@mail.gmail.com>
On Mon, Jul 13, 2009 at 4:32 PM, Matt Oliver<moliver at udel.edu> wrote:
> I tried to polygon.....you are right...too? slow
>
> I'm not sure how to resample to get a regular grid...using interp or interpp
> would not be wise because of the size of the grid...
A quick bodge would be to ignore the slight change in lat long
spacing and create some new lat-long coords that are regularly spaced:
> lonr=seq(min(lon),max(lon),len=length(lon))
> latr=seq(min(lat),max(lat),len=length(lat))
can you tell the difference between:
> image(lonr,latr,mcsst)
and:
> image(lon,lat,mcsst)
?
I'm sure the difference in longitude isn't worth worrying about:
> range(diff(lon))
[1] 0.01184845 0.01185608
but the lat difference is possibly significant:
> range(diff(lat))
[1] 0.008270264 0.009750366
I'll see if I can dig up my resampling code - basically it created new
coordinates for lat and lon (as above) and then mapped those into
indexes in the array.
Barry