____________
From: Kilpatrick, Katherine A
Sent: Thursday, November 29, 2018 3:40 PM
To: Michael Sumner
Cc: Elizabeth Webb; r-sig-geo@r-project.org
Subject: Re: [R-sig-Geo] Reproject MODIS data using R (results in NAs or no
spatial extent)
FYI The link that Mik
FYI The link that Mike provided is for ocean color products the GLASS land
products use a different grid…see this link.
https://modis-land.gsfc.nasa.gov/MODLAND_grid.html
K
On Nov 29, 2018, at 3:16 PM, Michael Sumner
mailto:mdsum...@gmail.com>> wrote:
To answer the "Any ideas on why
To answer the "Any ideas on why reprojecting this MODIS data is so
difficult?" - it's that the sinusoidal projection is actually pretty
exotic, it matches the daily aggregation of L2 geophysical variables from
individual swaths (multiple per day) into L3 bins (aggregated daily, then
into longer
gdalUtils has a get_subdatasets() which is really helpful here, it's
just a gdalinfo with grep.
Here's what my code looks like for reprojecting, what I found is that
even if the hdf has the proj definition, gdal tools don't seem to use it
right, so I specify the proj string as the s_srs or a_srs.
Ah, never mind - it's the subdataset discovery that's probably not easy
with rgdal.
Sorry for the noise.
Mike.
On Fri, 30 Nov 2018 at 06:38 Michael Sumner wrote:
> Fwiw there shouldn't be any need to convert from hdf to tif - could you
> please try this?
>
> x <- readGDAL(
Fwiw there shouldn't be any need to convert from hdf to tif - could you
please try this?
x <- readGDAL( ".../GLASS02B05.V04.A1990161.2018062.hdf")
If that works it at least removes some steps from your process.
Cheers, Mike.
On Fri, Nov 30, 2018, 05:03 Elizabeth Webb wrote:
> I am using