[R] extract function extracting only NA values

Fabio Berzaghi fabe at dmu.dk
Wed Apr 24 16:30:18 CEST 2013


ok so i imported the raster as a dataframe and everything is fine
    Min.   :-69826220
    Max.   :167780500

so it's something with values lower than -9999 that are interpreted as 
NA by the raster function


On 4/24/2013 13:52, Fabio Berzaghi wrote:
> So I think I might have found what is causing this problem. The values 
> I have in this raster
>
> > summary(ps0011yme)
>                 layer
> Min.        -9458.911
> 1st Qu.   1955256.000
> Median   10618870.000
> 3rd Qu.  79577730.000
> Max.    167780500.000
> NA's            0.000
>
> From ArcMap though I get different values
> min -69826224 and max 167780496
>
> So maybe when I am importing the rest are in R something goes wrong 
> and all the values below a certain threshold are considered NA. Is 
> this a bug or do I get to use a specific parameter for the raster 
> function?
>
>
> On 4/24/2013 13:10, Fabio Berzaghi wrote:
>> Hello,
>>
>> I have five raster files in ASCII format. With four of them I have no 
>> problem extracting values based on a set of X and Y coordinates. 
>> Unfortunately with one of the files all I managed to extract is NA 
>> values. To verify the problem I have opened the raster with ArcMap 
>> and there are no NA values where I am extracting. I have also plotted 
>> the spatial point class on top of the raster in R and it does 
>> correspond to the correct locations.
>>
>> These are some of the commands I am using, and as I already pointed 
>> out that works perfectly with other raster files.
>>
>> sp<-SpatialPoints(xysp)
>> xy$rasterimg<-extract(rasterimg,sp)
>>
>> Can anyone help? At this point I am rather clueless about this.
>>
>> Thanks
>>
>> ______________________________________________
>> R-help at r-project.org mailing list
>> https://stat.ethz.ch/mailman/listinfo/r-help
>> PLEASE do read the posting guide 
>> http://www.R-project.org/posting-guide.html
>> and provide commented, minimal, self-contained, reproducible code.
>
> ______________________________________________
> R-help at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide 
> http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.



More information about the R-help mailing list