[R-sig-Geo] problem with writeGDAL, its reading as black box in ARcGIS

[Edzer J. Pebesma]

Debs,

what does gdalinfo sayl about the .tif map? I wouldn't trust Arcmap as 
the ultimate referee to judge whether a geotiff is sound or not. It may 
dislike your floating point values, for example.
--
Edzer

Debarchana Ghosh wrote:
> Hello,
>
> I'm using Windows with 2 GB ram.
>
> I have created a SpatialGridDataFrame object by using the functions
> GridTopology and spkernel2d. I used writeGDAL to create a grid object
> readable in ArcGIS 9.2. But the output of the writeGDAL is a black box. I
> did try to switch from "stretched" to "categories" and I'm still getting a
> black box.
>
> Here is the code:
> bird06.xy<-coordinates(splanc_bird06[1:2])
> grd1<-GridTopology(c(419983.3, 4924198), c(30,30), c(3375,3499))
> k5000<-spkernel2d(bird06.xy, splanc_bird06$poly, h0=5000, grd1)
> df<-data.frame(k5000=k5000)
> kernel5000<-SpatialGridDataFrame(grd1, data=df)
> writeGDAL(kernel5000, fname="k2006.tif", drivername="GTiff")
>
> Also,
>
> The value of k2006.tif in Arcmap shows from High: 1.#QNAN to Low: 1.#QNAN
>
> But I know k5000 is not empty
>
> summary(k5000)
>
>         Min.      1st Qu.       Median         Mean      3rd Qu.
> Max.         NA's
> 0.000000e+00 0.000000e+00 1.240000e-08 6.191000e-08 5.724000e-08
> 6.325000e-07 3.248585e+06
>
> Thanks,
> Debs.
>
>