visinski profil hrvatske - nece mi dobro izrezati

[Dragica Šalamon]

Bok!

Zeza me dobiti kartu visinskog profila HR (sa susjedima) koristeći ovaj kod:

# elevation profile map

library(maptools)

library(raster)

srtm <- getData("SRTM", lon = 17, lat = 44.5)

#extent

e1 <- extent(10, # xmin

             20, # xmax

             40, # ymin

             50) # ymax

srtm_ct <- crop(srtm, e1)

# plot slope and aspect with hill shades

slope <- terrain(srtm_ct, opt = "slope")

aspect <- terrain(srtm_ct, opt = "aspect")

hill <- hillShade(slope, aspect, angle = 45, direction = 45, normalize = TRUE)

plot(hill, col = grey(0:100/100), legend = FALSE)

plot(srtm_ct, col = rainbow(25, alpha = 0.35), add = TRUE)

Ne znam zašto, ali izreže mi kartu 15 20 40 45, umjesto 10 20 40 50. Odnosno, uopće neće na sjever i zapad. Probala sam mijenjati i lon lat u getData... ali nista. 

E, sad kad skidam podatke, dobijem Warning, al ne znam jel to ima neke veze s tim: 

trying URL 'ftp://xftp.jrc.it/pub/srtmV4/tiff/srtm_40_04.zip'

trying URL 'http://hypersphere.telascience.org/elevation/cgiar_srtm_v4/tiff/zip/srtm_40_04.ZIP'

trying URL 'http://srtm.csi.cgiar.org/SRT-ZIP/SRTM_V41/SRTM_Data_GeoTiff/srtm_40_04.ZIP'

Content type 'application/x-zip-compressed' length 26250379 bytes (25.0 MB)

downloaded 25.0 MB

Warning messages:

1: In utils::download.file(url = aurl, destfile = fn, method = "auto",  :

  InternetOpenUrl failed: 'The server name or address could not be resolved

2: In utils::download.file(url = aurl, destfile = fn, method = "auto",  :

  InternetOpenUrl failed: 'The server name or address could not be resolved'

Ima tko ideju?

Ili alternativu?

Tenks,

Dragica

[Andreja Radović]

Dragice pozdrav,

problem se dogodio tijekom downloada podataka.

ako pogledaš strukturu objekta strm_ct ili slope, aspect

.@ extent  :Formal class 'Extent' [package "raster"] with 4 slots
  .. .. ..@ xmin: num 15
  .. .. ..@ xmax: num 20
  .. .. ..@ ymin: num 40
  .. .. ..@ ymax: num 45

ako pokušaš nacrtati veći obuhvat od podataka: 
spplot(srtm_ct, xlim=c(10,20))

dobiješ ovo:

ne znam da li je to do internetske veze ili dostupnosti podataka na webu

Predlažem da uzmeš podatke putem oznake države, znam da nije optimalno ako radiš više država


alt_hrv <- getData("alt", country="HRV")

spplot(alt_hrv)

Andreja

--
Primili ste ovu poruku jer ste se pretplatili na grupu "rtuntun" Google grupa.
Da biste otkazali pretplatu na ovu grupu i prestali primati njezine e-poruke, pošaljite e-poruku na adresu rtuntun+unsubscribe@googlegroups.com.
Posjetite ovu grupu na web-lokaciji https://groups.google.com/group/rtuntun.
Da biste pogledali ovu raspravu na webu, posjetite https://groups.google.com/d/msgid/rtuntun/CAJQ%2ByLV5WFeY0H%2ByO0nq2Ky%3DKH9z9qkCAUYnzLP_exdBciLrLA%40mail.gmail.com.
Dodatne opcije potražite na web-lokaciji https://groups.google.com/d/optout.

[Dragica Šalamon]

Hvala!

mda, bojala sam se da je tome tako, ako se dobro sjećam, mislim da sam isti problem imala i kad sam druge karte skidala ranije (za isti obuhvat)

ok, može i tak, valjda

-  a kak onda pridodam susjede? ili bolje da skinem negdje nekako neki prostorno veci raster nize rezolucije pa onda izrezujem?

Dragica

On Fri, Apr 21, 2017 at 10:49 AM Andreja Radović <andreja.r...@gmail.com> wrote:

Dragice pozdrav,

problem se dogodio tijekom downloada podataka.

ako pogledaš strukturu objekta strm_ct ili slope, aspect

.@ extent  :Formal class 'Extent' [package "raster"] with 4 slots
  .. .. ..@ xmin: num 15
  .. .. ..@ xmax: num 20
  .. .. ..@ ymin: num 40
  .. .. ..@ ymax: num 45

ako pokušaš nacrtati veći obuhvat od podataka: 
spplot(srtm_ct, xlim=c(10,20))

dobiješ ovo:

ne znam da li je to do internetske veze ili dostupnosti podataka na webu

Predlažem da uzmeš podatke putem oznake države, znam da nije optimalno ako radiš više država


alt_hrv <- getData("alt", country="HRV")

spplot(alt_hrv)

Andreja

Da biste otkazali pretplatu na ovu grupu i prestali primati njezine e-poruke, pošaljite e-poruku na adresu rtuntun+u...@googlegroups.com.

[Andreja Radović]

spajanje više rastera:

 alt_cro <- getData("alt", country="HRV")
> alt_bih <- getData("alt", country="BA")
trying URL 'http://biogeo.ucdavis.edu/data/diva/msk_alt/BIH_msk_alt.zip'
Content type 'application/zip' length 121735 bytes (118 KB)
downloaded 118 KB

> k <- merge( alt_cro,  alt_bih)
> spplot(k)

2017-04-21 12:00 GMT+02:00 Andreja Radović <andreja.r...@gmail.com>:

To ti je jako jednostavno

funkcija merge iz istog paketa raster.

Andreja

Hvala!

Da biste otkazali pretplatu na ovu grupu i prestali primati njezine e-poruke, pošaljite e-poruku na adresu rtuntun+unsubscribe@googlegroups.com.

[Dragica Šalamon]

izi-pizi
Hvala!

On Fri, Apr 21, 2017 at 12:03 PM Andreja Radović <andreja.r...@gmail.com> wrote:

spajanje više rastera:

 alt_cro <- getData("alt", country="HRV")
> alt_bih <- getData("alt", country="BA")
trying URL 'http://biogeo.ucdavis.edu/data/diva/msk_alt/BIH_msk_alt.zip'
Content type 'application/zip' length 121735 bytes (118 KB)
downloaded 118 KB

> k <- merge( alt_cro,  alt_bih)
> 

spplot(k)

Hvala!

Da biste otkazali pretplatu na ovu grupu i prestali primati njezine e-poruke, pošaljite e-poruku na adresu rtuntun+u...@googlegroups.com.

[Hr Ma]

test

izi-pizi
Hvala!

Hvala!

Da biste otkazali pretplatu na ovu grupu i prestali primati njezine e-poruke, pošaljite e-poruku na adresu rtuntun+unsubscribe@googlegroups.com.