- R-br - DINF UFPR

Remoção de coordinates() em SpatialGridDataFrame de diferentes dimensões
by ASANTOS 01 Out '15

01 Out '15

Caros Listeiros, Eu criei dos raster r1 e r2, respectivamente com dimensões 10x10 e 15x15, porém eu converti os rasters em SpatialGridDataFrame e gostaria de remover as coordenadas e respectivos valores do objeto 15x15 para que ele fique com o mesmo tamanho e nas mesmas coordenadas que o objeto 10x10, porém não estou conseguindo determinar quais pixels estão contidos em no objeto r1.grid que não existe em r2.grid, segue CRM abaixo, alguém poderia me ajudar? ### <code r> require(raster); require(rgdal); require(dplyr) ## Criando dois raster GTiff #10x10 r1 <- raster(nc=10, nr=10) r1 <- setValues(r1, round(runif(ncell(r))* 255)) writeRaster(r1, filename="r1.tif", format="GTiff", overwrite=TRUE) #15x15 r2 <- raster(nc=15, nr=15) r2 <- setValues(r2, round(runif(ncell(r))* 255)) writeRaster(r2, filename="r2.tif", format="GTiff", overwrite=TRUE) #Converter para SpatialGridDataFrame r1.gtiff<-readGDAL('r1.tif') r1.grid<- as(r1.gtiff, 'SpatialGridDataFrame') r2.gtiff<-readGDAL('r2.tif') r2.grid<- as(r2.gtiff, 'SpatialGridDataFrame') #Detecção de quais pixels estão contidos em r1.grid e r2.grid df.pixels<-anti_join(coordinates(r1.grid),coordinates(r2.grid), copy=TRUE, replace = TRUE) #Remoção das 125 coordenadas e do valor dos pixels de r2.grid para deixar os dois SpatialGridDataFrame com mesmo tamanho new.r2.grid<-r2.grid[which(df.pixels), ] Obrigado, -- ====================================================================== Alexandre dos Santos Proteção Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Caixa Postal 244 Avenida dos Ramires, s/n Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos(a)cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 ======================================================================

2 2

NaNs produzidos
by Samuel luna de almeida 30 Set '15

30 Set '15

Boa tarde grupo, Ao lognormalizar taxas da população em estudo, recebi de retorno do R uma mensagem que desconheço. Alguém sabe identificar? segue trecho do script: > baseG04$LOGNTX_J04 = log(tx_bruta_J04 -0.001) Warning message: In log(tx_bruta_J04 - 0.001) : NaNs produzidos Agradecido, Samuel

5 7

Validação cruzada da krigagem
by Wagner Wolff 29 Set '15

29 Set '15

Olá pessoal da lista, a dúvida é a seguinte! É possível utilizar testes como Shapiro-Wilk e Breusch Pagan para verificar, respectivamente a normalidade e homocedasticidade de resíduos oriundos de uma validação cruzada da krigagem? Att Wagner Wolff

5 7

sobre gdal_translate e gdalwarp
by Antonio Silva 29 Set '15

29 Set '15

Olá Tenho algumas dúvidas sobre o uso gdal no R. Depois de converter um arquivo HDF, com dados de temperatura de superfície do mar, para geotiff eu quero reprojetar para WGS84 ou SIRGAS2000 e contar entre as lats de 22 a 29 S e longs de 40 a 50 W. o HDF para exemplo está em https://app.box.com/s/16cf7qv6af6gsz1v66staori2mtneu0r Estou seguindo basicamente https://scottishsnow.wordpress.com/2014/08/24/many-rastered-beast/ Bem, converti o HDF com: gdal_translate("A20080012008031.L3m_MO_SST_4","georef.tif",sd_index=1,a_ullr=c(0,4320,8640,0), a_srs="+proj=eqc +lat_ts=0 +lat_0=0 +lon_0=0 +x_0=0 +y_0=0 +a=6371007 +b=6371007 +units=m +no_defs") Sem a_ullr e a_srs quando utilizava o gdalwarp eu recebia mensage: "ERROR 1: Unable to compute a transformation between pixel/line and georeferenced coordinates" Copiei os valores de a_ullr and a_srs values de GDALinfo("georef.tif") map <- raster("georef.tif") plot(map) Meu problema agora é reprojetar e cortar a imagem georef.tif: gdalwarp("georef.tif", "georef2.tif", s_srs="+proj=eqc +lat_ts=0 +lat_0=0 +lon_0=0 +x_0=0 +y_0=0 +a=6371007 +b=6371007 +units=m +no_defs", t_srs="+proj=longlat +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +no_defs") map2 <- raster("georef2.tif") plot(map2) Map2 não está em SIRGAS2000 e a opção de corte te=c(-50,-29,-40,-22) não funciona. Onde será que estou errando? Talvez em a_ullr e a_srs, não sei. Agradeço muito qualquer ajuda. Sds Antonio Olinto

1 1

R puro ou Editores (RStudio, Tin R...)
by Leonardo Monteiro 29 Set '15

29 Set '15

Pessoal, bom dia. Provavelmente o que ocorreu comigo já ocorreu com bastante mais gente também! Eu estou criando algumas rotinas no R (RStudio) e algumas vezes "não roda". Entretanto, quando eu copio o codigo para o R (puro) ele roda normalmente. Isso é erro de versão? Ou do próprio programa (RStudio)? Outro problema que frequentemente ocorre é no caso do uso da condicional "if" pq as vezes se eu faço: if(a>10) {x} else { y } isso roda mas "nao funciona" e se eu escrevo de outra maneira, como, por exemplo: { if (a>10) {x} else {y} } aí roda numa boa ("mesmo assim, em algumas situações ocorre o msm da situação acima). Qual seria a maneira mais adequada de escrever essa estrutura do if? PS: por favor, quem for responder, responda para o email: monteiroleonardo6(a)gmail.com pq nao tenho acesso a lista do meu trabalho. obrigado, Leonardo

2 1

Informação do conteúdo de uma variável como definição de outra variável
by SILVIA PATRICIA BARRETO SANTANA 28 Set '15

28 Set '15

Muito boa tarde a todos! Pessoal, Consegui resolver meu problema! Para que possa ajudar a outros, consegui resolver meu problema com uma sugestão anterior dada pelo prof. Éder Comunello. Para o problema atual utilizei a seguinte rotina res1m9_1$c<- sapply(seq(nrow(res1m9_1)), function(x) { sel1m9_1<- seq(res1m9_1$ini[x], res1m9_1$fim[x]) print(sum(DC_MHN1m9[sel1m9_1,17]), na.rm=TRUE) }) res1m9_1$chu1<-rep(0, times=length(res1m9_1$lengths)) res1m9_1$chu1[res1m9_1$c>0]=1 Muitíssimo obrigado a todos. Prof. David muito obrigada. Cordialmente, Silvia

1 0

Informação do conteúdo de uma variável como definição de outra variável
by SILVIA PATRICIA BARRETO SANTANA 28 Set '15

28 Set '15

Prof. David, boa tarde! Obrigada pela sugestão contudo, nem sempre o vetor tem comprimento três, este valor é variável, podendo ter 20 elementos de comprimento, por exemplo. Acho que a rotina a ser construída deve pegar a informação da variável "ini" e "fim" que será correspondente a posição inicial da variável "ind" e a posição posição final, este novo vetor, conterá a quantidade de números que irei testar se zero ou um. Atenciosamente, Silvia

1 0

Informação do conteúdo de uma variável como definição de outra variável
by SILVIA PATRICIA BARRETO SANTANA 27 Set '15

27 Set '15

Caro David, Tenho dois bancos dt1 e dt2, preciso forma o dt3. Por exemplo, comp=3 então devo pegar as 3 linhas do ind, em dt2, verificar se possui algum valor "1" e retornar em dt3 o valor ind_2=1, pois (1,0,1). Cordialmente, Silvia

2 1

RMarkdown - Imagem externa
by Andre Oliveira 27 Set '15

27 Set '15

Boa tarde,inseri uma imagem no meu arquivo RMarkdown. Porém ficou muito grande a figura. Qual é o parâmetro para coloca-la em 4 cm por 3 cm? ![ifes](/home/andre/Arquivos/logo34.png) André Oliveira Souza. Graduação em Matemática, mestrado em estatística aplicada.Instituto Federal de Educação, Ciência e Tecnologia do Espirito Santo. IFES

4 8

Informação do conteúdo de uma variável como definição de outra variável
by SILVIA PATRICIA BARRETO SANTANA 27 Set '15

27 Set '15

Prezados Gostaria de obter ajuda para juntar informação de um dataframe e colocar em outro dataframe, por exemplo, tenho: dt1 Comp ini fim 1 3 1 3 2 4 4 7 3 1 8 8 4 3 9 11 outro dataframe dt2 ind 1 1 2 0 3 1 4 0 5 1 6 0 7 0 8 0 9 1 10 0 11 0 Gostaria de obter um terceiro dataframe da seguinte forma: dt3 Comp ini fim ind_2 1 3 1 3 (1,0,1)=1 2 4 4 7 (0,1,0,0)=1 3 1 8 8 (0)=0 4 3 9 11 (1,0,0)=1 Obs. A variável ind_2 é indicadora da condição: se "0" ou "1" então "1", apenas "0" se todos forem iguais a "0". Os três valores dentro dos parênteses são referentes à var "Comp" do df1 com a sequência da variável "Ind" no dt2, . Grata. Silvia

2 1