Dear Alexandre -
Apologies if your question has been answered already, I read r-sig-geo
as a daily digest.
I'm not sure what is "obviously' wrong with the humberside example.
Because the observation window is an irregular polygon the quadrats on
the edge have an irregular shape. Their cou
Jonathan,
I used the file you send me. This shows that the problem is already
present in rgdal (i.e. before it gets passed to raster)
> x <- readGDAL("stand_bm.gis")
stand_bm.gis has GDAL driver LAN
and has 839 rows and 1758 columns
> summary(x)
(...)
Data attributes:
Min. 1st Qu. MedianM
Dear Everton,
> Regarding the rotine that identifies the possible cluster locations, I do
> not know which value should I consider for the fractpop. The spatial
> clustering is very associated to this value and, depending of the considered
> value for the fractpop, I get a different spatial confi
On Thu, Aug 20, 2009 at 9:57 PM, Jonathan
Thompson wrote:
> gdalinfo:
> Band 1 Block=1758x1 Type=Int16, ColorInterp=Undefined
> Min=0.000 Max=48499.000
>> Arcgis:
>> Pixel Type = unsigned integer
>> Pixel Depth = 16 bit
Arcgis is saying 'unsigned integer' but gdalinfo says 'Int16'. But
then gd
Thank you Barry and Robert
This is all very interesting. When I use the FWTools Shell and the
gdalinfo command it says that the Min value is zero (see below) and max
is 48499. This is different than the summary in Rgdal, which produced
the negative values. So I guess the problem is in the
> What are these .gis files? What driver is rgdal using to read them?
Ah, I've just noticed this:
> gdalDrivers()[38,]
name long_name create copy
38 LAN Erdas .LAN/.GIS FALSE FALSE
>
The output from GDALinfo will still be interesting, as will the
output from the command line g
On Thu, Aug 20, 2009 at 2:58 PM, Jonathan
Thompson wrote:
> Dear List,
>
> I use a forest landscape model, which produces 16 bit unsigned ".gis" raster
> files. I often read and manipulate them using rgdal and the raster package
> in R.
>
> However, sometimes, when I read in a unsigned 16 bit .gis,
Jonathan,
Do you only get negative values with files that should have values >
32767 ? If so, that would suggest that the .GIS files are always
interpreted as signed. I believe that values in .GIS are by definition
unsigned (as is common with remote sensing formats), and that there is
no flag for
Dear all,
I am doing some spatial analysis using the library DCluster. In order to
identify possible "cluster locations" I apply the Kulldorff's method. I am
dealing with some technical problems, however.
Regarding the rotine that identifies the possible cluster locations, I do
not know which val
Dear all,
I am doing some spatial analysis using the library DCluster. In order to
identify possible "cluster locations" I apply the Kulldorff's method. I am
dealing with some technical problems, however.
Regarding the rotine that identifies the possible cluster locations, I
do not know which val
Dear List,
I use a forest landscape model, which produces 16 bit unsigned ".gis"
raster files. I often read and manipulate them using rgdal and the
raster package in R.
However, sometimes, when I read in a unsigned 16 bit .gis, I get the
following:
stnd <- raster("stand_bm.gis", values = T
11 matches
Mail list logo
