Gad, whoever thought that delivering gridded data in this way would be a
good idea?

Do you see any positives at all? It amazes me.

Cheers, Mike
On 15 Feb 2014 01:31, "Barry Rowlingson" <b.rowlingson@lancaster.ac.uk>
wrote:

> On Thu, Feb 13, 2014 at 5:55 AM, Steve Pickering <
> pickering@penguin.kobe-u.ac.jp> wrote:
>
> >
> > require(raster)
> > forestCanopyIGH = raster("test.vrt")
> > plot(forestCanopyIGH,maxpixels=100*100)
> >
> > I need help with the next step of displaying the unprojected version!
> >
>
> Just to sign off this thread, I managed to fix this for Steve. The
> complications were:
>
>  1. original data had no coordinates apart from row and column count, so I
> transformed a full earth raster in epsg:4326 to "+proj=igh" to get the
> coordinates used in an IGH projection.
>
>  2. original data wasn't the full earth - it was an IGH truncated to chop
> off antarctica (no trees there) and the western and eastern oceanic parts
> of the IGH projection (no trees there either). Fixed that by visually
> working out the extremes of the IGH raster on a full world IGH raster, and
> setting the extent of the clipped IGH raster accordingly.
>
>  3. Took a deep breath and projectRaster on the clipped IGH raster using
> epsg:4326, then plotting it overlaid on the epsg:4326 full world raster
> from step 1 (downloaded from Natural Earth) for comparison.
>
> The alignment was pretty good, Steve says good enough!
>
> Barry
>
>         [[alternative HTML version deleted]]
>
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>

	[[alternative HTML version deleted]]

