f you don't want to use NASA's website. I
> am using R 4.3.1 and RStudio 2023.06.2+561.
>
> --
> Tziokas Nikolaos
> Cartographer
>
> Tel:(+44)07561120302
> LinkedIn <http://linkedin.com/in/nikolaos-tziokas-896081130>
>
> [[alternative HTML version deleted]]
>
>
ing wrong that could
> result in inaccurate location data? Thanks!
> ___
> R-sig-Geo mailing list
> R-sig-Geo@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>
--
Michael Sumner
Software and Database Engineer
Aust
play it, it says:
>
> > rw_base_blocks
> stars object with 3 dimensions and 1 attribute
> attribute(s), summary of first 1e+05 cells:
> Min. 1st Qu. MedianMean 3rd Qu. Max.
> X 0 217224 219.915 238 255
> dimension(s):
> from to offsetdel
(171.7781 24.66314)
> >
>
> I expected the POINTS to look like:
> 1 POINT (-76.61246 39.35010)
>
> Can anyone suggest what I'm doing wrong? Thanks so much in advance. I've
> worked on it all day today, without making much progress.
>
> -Kevin
>
>
> >
>
> [[alternative HTML version deleted]]
>
> ___
> R-sig-Geo mailing list
> R-sig-Geo@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>
--
Michael Sumner
Software and Database Engineer
Australia
proving the documentation, fixing bugs - if you're interested in
this and want to discuss I'm also happy to help.
The file I used in the example above is posted in the related GDAL PR (pull
request):
https://github.com/OSGeo/gdal/issues/6195
All the best!
Mike
--
Michael Sumner
Software and Data
23 at 9:36 AM Michael Sumner wrote:
> One way is to write it to a (COG) GeoTIFF with overviews ("pyramid", or
> zoom levels) - these are pre-calculated copies of the highest resolution
> data saved as lower resolution versions.
>
> Do you want this to generate a set
Institutional website: ICOMVIS
> <http://www.icomvis.una.ac.cr/index.php/manuel>
> Blog sobre Ciencia de Datos: https://mspinola-ciencia-de-datos.netlify.app
>
> [[alternative HTML version deleted]]
>
> ___
> R-sig-Geo mail
have a look at the output of printing nc_data, you are looking for
dimensions that belong to a 2D (or higher) variable of interest .
They can be called anything, but longitude, lon, x, etc are common. you can
use ncdf4 functions to find them (without reading the output), but it's all
stored in a
; Thanks,
> Alex
>
> [[alternative HTML version deleted]]
>
> ___
> R-sig-Geo mailing list
> R-sig-Geo@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>
--
Michael Sumner
Software and Database Engineer
Australian Antarctic
9nJNwPKlJaPlkHn_NGzmjMHMaoUevH04
>
> Request if someone can guide me and thanking you in advance
> sownalc
>
> [[alternative HTML version deleted]]
>
> ___
> R-sig-Geo mailing list
> R-sig-Geo@r-project.org
> https
lol, as a matrix - does it just?
very clever
On Wed, 27 July 2022, 22:36 Ben Tupper, wrote:
> Hi,
>
> I have never observed st_extract() fail to preserve row order.
>
> Your workflow works well. The documentation for st_extract() mentions
> that if you might see better performance for large
ify.app/
>
> [[alternative HTML version deleted]]
>
> _______
> R-sig-Geo mailing list
> R-sig-Geo@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
--
Michael Sumner
Software and Database Engineer
Australian Antarctic Di
ot see how to get to it and apply it to the layer.
> >>> There is also another file with complementary data related to
> >>> geolocation
> >>> but how can I query it, and get the data for geolocation and apply it
> to
> >>> the raster layer?.
> >>&
try terra::rast() or stars::read_stars() on the file, both use gdal to
interrogate and read. What you get depends on the structure of the files,
variables interpreted as subdatasets. stars is more general but sees
variables as 2D arrays with bands, terra calls these bands layers.
Alternatively
On Tue, Sep 21, 2021 at 8:04 PM Roger Bivand wrote:
> On Sun, 19 Sep 2021, Michael Sumner wrote:
>
> > Hi Jean-Luc, if the packages mentioned don't yet cover it, you can get
> > ogrInfo and ogrFIDs with a combination of vapour's vapour_read_names()
> > (for FID), vapour_
vance,
>
> Jean-Luc Dupouey
>
> --
> INRAE
> UMR Silva
> F-54280 Champenoux
> France
>
> ___
> R-sig-Geo mailing list
> R-sig-Geo@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>
--
Michael Su
> f1 <- ppm(out.ppp~1)
> E <-envelope(f1, Kinhom, nsim = 19, global = TRUE, correction = "border")
> plot(E)
>
> Best wishes,
>
> Alexandre
>
>
>
>
>
>
>
>
> Em sexta-feira, 30 de julho de 2021 15:41:08 BRT, Michael Sumner <
mdsum..
Don't use UTM, it's a legacy of a bygone era. Use a local custom projection
such as laea:
https://twitter.com/mdsumner/status/1136794870113218561?s=19
I know it's common advice, utm for projection choice but it's trash
propagated by folks who should know better. Traversing zones for projection
> work on a SpatExtent either to convert a SpatExtent to a SpatVector.
>
> Thanks,
>
> Tim
>
> ___
> R-sig-Geo mailing list
> R-sig-Geo@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>
--
Michael Sumner
Softw
j=aeqd +R_A +lat_0=", YMiddle,
> " +lon_0=", XMiddle,
> " +x_0=0 +y_0=0 +ellps=WGS84",sep="")
>
> might not be future-proof. The second option could be a local UTM based
> CRS, getting the zone from
Roger
>
>
> Thanks again,
>
> Regards, Micha
>
>
>
>
> --
> Micha Silver
> Ben Gurion Univ.
> Sde Boker, Remote Sensing Lab
> cell: +972-523-665918
>
> ___
> R-sig-Geo mailing list
> R-sig-Geo@r-project.
ta from layers for the shape? If yes, what I need to change to reproject
> the shape for the correct CRS?
>
> Thanks.
>
> Pietro Andre Telatin Paschoalino
> Doutorando em Ciências Econômicas da Universidade Estadual de Maringá -
> PCE.
>
> [[alternative HTML
gt; or
>
> Which(r>0,cells=T)->p
>
> cbind(xyFromCell(r,p),extract(r,p))
>
>
>
> what’s the terra alternative?
>
>
>
>
>
>
>
>
>
> *From:* Michael Sumner
> *Sent:* Tuesday, 14 July 2020 2:41 PM
> *To:* Herr, Alexander (L, Black M
: object of type 'S4' is not subsettable
>
> What am I doing wrong?
>
> Thanks
> herry
>
>
> in terra
>
>
> ___________
> R-sig-Geo mailing list
> R-sig-Geo@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>
--
Michael S
Hi, your longlat coordinates are in the western hemisphere (-119, 39) but
EPSG 27700 is designed for (and appropriate for some scenarios) somewhere
in the UK. You need a local projection suitable for the area you are in.
Look for an EPSG for your region, or maybe Albers conic or Lambert
Azimuthal
You need to also setExtent() of the raster. That's the coordinates that the
crs refers to xmin,xmax,ymin,ymax values in metres. Sometimes visual cues
can be used to derive these from an image, sometimes the grid is a commonly
used convention (a certain dimension, I.e. extent +resolution in a
> the next few days...
>
> So it definitely seems to be a problem of rgdal. It would be great if it
> could still be solved.
>
> Best,
>
> Thorsten
>
>
>
> Am 27.04.2020 um 15:58 schrieb Michael Sumner:
> > Try stars it worked for me on a t
g this thread could also benefit;
> if I'm wrong, that would also be good to know.
>
> Roger
>
>
> >
> > Best,
> >
> > Thorsten
> >
> >
> >
> > Am 27.04.2020 um 15:58 schrieb Michael Sumner:
> >> Try stars it worked for me on a
Try stars it worked for me on a test
On Mon., 27 Apr. 2020, 23:54 Thorsten Behrens,
wrote:
> Roger,
>
> thanks a lot for your reply!
>
> I have 256GB RAM installed (mentioned it somewhere). And there, all is
> fine when I run:
>
> rDemTest = raster(nrow = 48000, ncol = 72000, ext = extent(c(0,
l resolution is supposed to be 750m but this shows 1°. What am I
> doing wrong? Any help will be greatly appreciated. Thank you.
>
> Sincerely,
>
> Millu
>
> [[alternative HTML version deleted]]
>
> ___
> R-sig-Geo
;
> The produces what looks like a raster with values. Is that all I need to
> do or should I be reading these data another way?
> Thanks
> Julie
>
>
> [[alternative HTML version deleted]]
>
> ___
> R-sig-Geo mailin
In raster try cellFromXY() and adjacent() which separates out the two
tasks, find a cell and find cell neighbours
On Wed., 6 Nov. 2019, 17:18 Bede-Fazekas Ákos, wrote:
> Dear Cristabel,
> function focal() of package raster is what you are searching for.
> focal(x, w = matrix(1/9, ncol = 3,
Cool, thanks for clarifying
On Fri., 18 Oct. 2019, 02:48 Marta Rufino, wrote:
> Thank you very much to all :)
>
> Barry+Roger contributions sorted the problem for the moment, although in
> my mind there are still things to be understood.
> Mike, yep, thank you for alerting me of
rnaturalearth is really dirty data, there's no way forward when things are
so bad, only backing up can help.
Do you need the whole world? Where do you need, and what scale? Your
question suggests that finding the right data is the solution rather than
battling with software.
Cheers, Mike
On
really
motivated to pursue before.
Cheers, Mike.
On Sat, Sep 14, 2019 at 7:49 PM Michael Sumner wrote:
> You should set a target raster with the extent and dimensions required.
> There are inherent limits in reprojection and heuristics won't always work.
> Generally using a target r
You should set a target raster with the extent and dimensions required.
There are inherent limits in reprojection and heuristics won't always work.
Generally using a target raster is much more efficient anyway.
But, this is an intensive remodeling of the data, delivered in a projection
for good
namespace (and not attached):
> [1] compiler_3.5.1 rgdal_1.4-3 Rcpp_1.0.1 codetools_0.2-15
> ncdf4_1.16.1
> [6] grid_3.5.1 lattice_0.20-35
>
> Best regards,
>
> Frederico
>
> [[alternative HTML version dele
3-15 DBI_1.0.0
> tools_3.5.3units_0.6-2Rcpp_1.0.1 grid_3.5.3 e1071_1.7-0.1
> [10] classInt_0.3-1
>
> [[alternative HTML version deleted]]
>
> ___
> R-sig-Geo mailing list
> R-sig-Geo@r-project.org
> htt
There's code in trip for the angles, but arranging the underlying coords
and expressing those in the right form will be the crux. st_coordinates
gives what you need, or pure data frame from spbabel::sptable is what I'd
start with.
I'll take a look when I can.
Cheers, Mike
On Sat., 1 Jun. 2019,
>
>
> I tried switching "ncfname" with "Model1", but I then receive this error:
>
> Error in (function (classes, fdef, mtable) :
> unable to find an inherited method for function ‘brick’ for signature
> ‘"ncdf4"’
>
>
> Is there a rea
obtain a raster of land use and land cover of
> Europe. Which site to choose?
>
> Thanks,
>
> Lara
>
> [[alternative HTML version deleted]]
>
> ___
> R-sig-Geo mailing list
> R-sig-Geo@r-project.org
> https://stat.e
h
> lon Size:128
> units: degree East 3 global attributes:
> description: Annual global maximum precipitation from the CanESM2
> 1pctCO2 scenario
> history: Created Mon Jun 4 11:24:02 2018
> contact: rain1...@aim.com
> So, in ge
7))
> # remove the background and default gridlines
> theme_void()
>
> ## coordinates dataframe
> dput(df)
> structure(list(lon = c(2.67569724621467, 17.5766416259819,
> 28.4126232192772,
> 23.8147674538232, 29.8917589327105), lat = c(28.1503115976162,
> -12.338878738020
add = TRUE)
## convince yourself the extent is correct with things like:
abline(h = c(-90, 90), col = "red")
That use of "datum=WGS84" may well be incorrect, but compared to trying to
project data from a completely mis applied projection it's an improvement.
HTH
Cheers,
Ah, never mind - it's the subdataset discovery that's probably not easy
with rgdal.
Sorry for the noise.
Mike.
On Fri, 30 Nov 2018 at 06:38 Michael Sumner wrote:
> Fwiw there shouldn't be any need to convert from hdf to tif - could you
> please try this?
>
> x <- readGDAL( &quo
xmax, ymin, ymax)
> #coord. ref. : +proj=longlat +ellps=clrk66 +no_defs
> #values : -32768, 32767 (min, max)
>
> Any ideas on why reprojecting this MODIS data is so difficult?
>
>
> [[alternative HTML version deleted]]
>
> ___
an you; it is explained here:
>
> https://r-spatial.github.io/stars/articles/data_model.html
>
> stars uses some functions directly from GDAL which it doesn't expose to
> the user, but there is no magic going on here.
>
> >
> >
> > Thanks in advance for any tips.
> > Tim
> >
> > ___
___
> > R-sig-Geo mailing list
> > R-sig-Geo@r-project.org
> > https://stat.ethz.ch/mailman/listinfo/r-sig-geo
> >
>
> --
> Edzer Pebesma
> Institute for Geoinformatics
> Heisenbergstrasse 2, 48151 Muenster, Germany
> Phone: +49 251 8333081 &
d]]
>
> ___
> R-sig-Geo mailing list
> R-sig-Geo@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>
--
Dr. Michael Sumner
Software and Database Engineer
Australian Antarctic Div
elegant
> > mask <- crop(polys,e)
> > polys2 <- polys[mask,]
> > plot(polys2,add=T,col="green")
> >
> > Thanks a lot. Best regards
> >
>
> ___
> R-sig-Geo mail
That doesn't help though!
Sorry, I thought that was it, will sleep on it.
Cheers, Mike
On Fri, 24 Aug 2018 at 22:00 Michael Sumner wrote:
> Sorry, I was wrong - the actual issue is that the writeRaster argument is
> "datatype", it gets conflated with the function that is d
overwrite=T)
> s<-raster("test_int.tif")
> dataType(s)
>
> On Fri, 24 Aug 2018 at 18:05, Roger Bivand wrote:
>
>> Thanks, Mike. I agree that the lower level interface in rgdal is flexible
>> enough, but as you say non-trivial. Nevil: could you please provide a small
&
uot;INT2U",overwrite=T)
> > s<-raster("test_int.tif")
> > dataType(s)
> [1] "FLT8S"
> >
>
>
> Can you suggest how I ensure the values are stored as integer?
>
> Many thanks
>
> [[alternative HTML version deleted]]
>
> ___
-32 -23 1
> >>> 3 -32 -22 1
> >>> 4 -33 -22 1
> >>> 5 -32 -23 2
> >>> 6 -31 -23 2
> >>> 7 -31 -22 2
> >>> 8 -32 -22 2
> >>> 9 -31 -23 3
> >>> 10 -30 -23 3
> >>> 11 -3
was 50 km2 smaller.
>
> What could be the cause?
>
> Thanks in advance,
> N.
>
> [[alternative HTML version deleted]]
>
> ___
> R-sig-Geo mailing list
> R-sig-Geo@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-si
ace.Temperature
> > # zvar: sst
> >
> > S = raster::raster(basename(URI), varname = "sst", crs = '+proj=eqc
> +lat_ts=0 +lat_0=0 +lon_0=0 +x_0=0 +y_0=0 +ellps=WGS84 +datum=WGS84
> +units=m +no_defs')
> > # Warning message:
> > # In .getProj(prj,
ão Paulo
> Rua do Matão, 1226
> São Paulo-SP - Brasil -
> 05508-090
> +55-11-97425-3791
> Skype: sergio_ibarra1
>
> [[alternative HTML version deleted]]
>
> ___
> R-sig-Geo mailing list
> R-sig-Geo@r-project.org
>
ent : -180, 180, -90, 90 (xmin, xmax, ymin, ymax)
> > coord. ref. : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0
> > data source : in memory
> > names : clumps
> > values : 1, 1 (min, max)
> >
> > __
e
> Old Dominion University
> Department of Biological Sciences
> Mills Godwin Building, Room 312
> Norfolk VA 23529 USA
>
> [[alternative HTML version deleted]]
>
> ___
> R-sig-Geo mailing list
> R-sig-Geo@r-project.org
>
is, I haven't found anything promising.
>
> Thanks,
> Kent
>
> [[alternative HTML version deleted]]
>
> _______
> R-sig-Geo mailing list
> R-sig-Geo@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>
___
> R-sig-Geo mailing list
> R-sig-Geo@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>
--
Dr. Michael Sumner
Software and Database Engineer
Australian Antarctic Division
203 Channel Highway
Kingston Tasmania 7050 Au
ntinuous line in this example? Many
> thanks, Andy
>
> [[alternative HTML version deleted]]
>
> ___
> R-sig-Geo mailing list
> R-sig-Geo@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-g
rsity
> Department of Biological Sciences
> Mills Godwin Building, Room 312
> Norfolk VA 23529 USA
> ___
> R-sig-Geo mailing list
> R-sig-Geo@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>
--
Dr. Michael Sumner
S
___
> > R-sig-Geo mailing list
> > R-sig-Geo@r-project.org
> > https://stat.ethz.ch/mailman/listinfo/r-sig-geo
> >
>
> --
> Roger Bivand
> Department of Economics, Norwegian School of Economics,
> Helleveien 30, N-5045 Bergen, Norway.
> voice: +47 55
sig-Geo mailing list
> R-sig-Geo@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>
--
Dr. Michael Sumner
Software and Database Engineer
Australian Antarctic Division
203 Channel Highway
Kingston Tasmania 7050 Australia
[[alternative HTML version deleted]]
_
may contain information
> > that is confidential and legally privileged. Please do not read, copy,
> > forward, or store this message unless you are an intended recipient of
> > it. If you have received this message in error, please forward it to
> > the sender and delete it completely from your computer system.
> >
> > --
> > This message was scanned by Libra ESVA and is believed to be clean
_
> r-h...@r-project.org mailing list -- To UNSUBSCRIBE and more, see
> 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 cod
= ) calls into translations for rgdal's output.dim, but it
would be powerful to be easy to do that.
Cheers, Mike.
On Fri, 9 Feb 2018 at 10:32 Michael Sumner <mdsum...@gmail.com> wrote:
> Oh, I had a typo in the xml_specification due to blithe copy/paste from
> email. It works fine
t;https://maps.google.com/?q=4210+University+DriveAnchorage,+AK+99508=gmail=g>
> Anchorage, AK 99508
> <https://maps.google.com/?q=4210+University+DriveAnchorage,+AK+99508=gmail=g>
>
> voice: (907) 786-7145
>
> http://alaska.usgs.gov/science/biology/walrus/
>
> O
t; (projWinParameters = paste('-projwin', bb[1, 1], bb[2,2], bb[1,2],
> bb[2,1], sep = " "))
> aqua <- readGDAL(fname = xml_specification) # , projWinParameters) ?? how
> to implement the projwin constraints??
> (aqua_stack <- stack(aqua))
> #
> plotRGB(aqua_st
> plot(int,add=TRUE,col="red") # identifies more than just the crossings
> >
> >
> >
> >
> > [[alternative HTML version deleted]]
> >
> > ___
> > R-sig-Geo mailing list
> > R-
ossing.psp(PSP)
> >
> > plot(int,add=TRUE,col="red") # identifies more than just the crossings
> >
> >
> >
> >
> > [[alternative HTML version deleted]]
> >
> > ___
> > R-sig-Geo mai
Ugh, and you must think I'm crazy given that it's clearly HDF5 not HDF4 as
per my opening line (the file naming is more the latter). I'll stop now.
Hope it's useful.
Cheers, Mike.
On Sat, 2 Dec 2017 at 07:44 Michael Sumner <mdsum...@gmail.com> wrote:
> Ah thanks, and sorry I didn't
quot;Q20151522015181.L3m_MO_SCIA_V5.0_SSS_1deg"
>>
>> nc = ncdf4::nc_open(filename)
>> m = ncdf4::ncvar_get(nc, 'l3m_data')
>> ncdf4::nc_close(nc)
>> r = raster::raster(t(m))
>> rasterVis::levelplot(r)
>>
>> ### END
>>
>> >
gt; and variations with band and layer.
>
> I would greatly appreciate any suggestions to solve this issue.
>
> Thanks
>
> --
> Antônio Olinto Ávila da Silva
> Fisheries Institute
> São Paulo, Brasil
>
> [[alternative HTML version deleted]]
>
> ___
t; > > > >
> > > > > >
> > > > > > Many of us rely on the raster package as a crucial basis for
> their
> > > > work.
> > > > > Big thank to Robert. Nevertheless to contribute needs to focus
> needs
> > > and
> > >
t;> Date: Mon, 23 Oct 2017 22:53:22 +
>> From: "Mortimer, Linsey Anne" <l.mortimer...@aberdeen.ac.uk>
>> To: Michael Sumner <mdsum...@gmail.com>
>> Cc: "r-sig-geo@r-project.org" <r-sig-geo@r-project.org>
>
>
>> Subject: Re: [R-si
_United Kingdom.1252 LC_NUMERIC=C
>
> [5] LC_TIME=English_United Kingdom.1252
>
>
>
> attached base packages:
>
> [1] stats graphics grDevices utils datasets methods base
>
>
>
> other attached packages:
>
> [1] maptools_0.9-2 rgdal_1.2-13 sp_1.
_
> R-sig-Geo mailing list
> R-sig-Geo@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>
--
Dr. Michael Sumner
Software and Database Engineer
Australian Antarctic Division
203 Channel Highway
Kingston Tasmania 7050 Australia
[[alternative HTML vers
ellowDepartment of Climate and Space Science and
> EngineeringUniversity of Michigan
> [[alternative HTML version deleted]]
>
> ___
> R-sig-Geo mailing list
> R-sig-Geo@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r
ter_2.5-8 sp_1.2-5
>
> loaded via a namespace (and not attached):
> [1] rgdal_1.2-11parallel_3.3.3 tools_3.3.3 Rcpp_0.12.12
> grid_3.3.3
> [6] lattice_0.20-35
>
> ___________
> R-sig-Geo mailing list
> R-sig-Geo@r-project.org
>
6707267, 9941473, 9950408 (xmin, xmax, ymin, ymax)
> coord. ref. : +proj=laea +lat_0=-100 +lon_0=6370997 +x_0=45 +y_0=0
> +datum=WGS84 +units=m +no_defs +ellps=WGS84 +towgs84=0,0,0
>
>
> c<-crop(tmp, ex)
>
> Error in .local(x, y, ...) : extents do not overlap
>
>
>
Son <kkyon...@hotmail.com> wrote:
> I did the exact same way, but the output of the reprojection of 'a' raster
> did not properly align with 'b' raster.
>
>
>
>
>
> *Kyongho Son*
> Postdoctoral Fellow
> Research Foundation of The City University of New York
>
>
>
&g
ersion deleted]]
>
> _______
> R-sig-Geo mailing list
> R-sig-Geo@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>
--
Dr. Michael Sumner
Software and Database Engineer
Australian Antarctic Division
203 Channel Highway
> I’d be happy to trial it on data here if that would help.
>
>
>
Thanks! I'll take you up on that when I look at it again.
Cheers, Mike.
> Cheers
>
> Terry
>
>
>
>
>
> *From:* Michael Sumner [mailto:mdsum...@gmail.com]
> *Sent:* Thursday, 7 September 2017 6:
his message and any copies of
> this message from your computer and/or your computer system network.
> --
>
> [[alternative HTML version deleted]]
>
> ___
> R-sig-Geo mailing list
> R-sig-Geo@r-pr
to be more expert on this cross-over topic, still working on it :)
Cheers, Mike.
On Sat, 2 Sep 2017 at 20:27 Rolf Turner <r.tur...@auckland.ac.nz> wrote:
> On 02/09/17 19:00, Michael Sumner wrote:
>
>
>
> >
> > It's not clear to me if you need polygons in spatstat or a r
ooth_2.23-15 stringi_1.1.5
> [73] osmar_1.1-7 Rcpp_0.12.10 sf_0.5-3
> [76] png_0.1-7leaflet_1.1.0coda_0.19-1
> >
>
> ---
> This email has been checked for viruses by AVG.
> http://www.avg.com
>
> ___
mm_day_year <- cbind(joined_mm_day_year, joined_mm_day)}
>
> }
>
> #Join mm_day with dams codes
> joined_mm_day_year_names=cbind(dams_codes,joined_mm_day_year)
>
> newnames <-t(as.data.frame(c(paste ("day",
> seq(1,(ncol(joined_mm_day_year_names)-3),1)
;- as.data.frame(data_dams_size)
>
> Unfortunately, I get mostly NAs for the data, it seems that there's an
> issue with the CRS projections for the city coordinates. Is there a
> specific projection for city level coordinates? Or am I doing something
> completely wrong? Thank you!
&g
100, 180, 0, 80))
>
> # rotate the extent
> eExtent <- raster::extent(eR)
> raster::extent(eR) <- c(-260, -180, eExtent[3:4])
>
> # western section
> wR <- raster::crop(R, c(-180, -90, 0, 80))
>
> # merge
> newR <- raster::merge(eR, wR)
> newR
>
>
o-cgal-or-equivalent-toolset-in-r
>
> Thanks in advance,
>
> Robin
>
> [[alternative HTML version deleted]]
>
> ___________
> R-sig-Geo mailing list
> R-sig-Geo@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-
from a general framework, which I'm working on - that's
some kind of excuse for why this and related projects are quite raw and
unfinished.)
Cheers, Mike.
Cheers,
> Ben
>
> On Jun 22, 2017, at 4:46 AM, Michael Sumner <mdsum...@gmail.com> wrote:
>
>
> It used to do
extent and resolution, and the spatial points are in the same
> coordinate system.
>
> Any help would be greatly appreciated.
> Wade
>
> ___
> R-sig-Geo mailing list
> R-sig-Geo@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-si
low Drive, P.O. Box 380
> East Boothbay, Maine 04544
> http://www.bigelow.org
>
> ___
> R-sig-Geo mailing list
> R-sig-Geo@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>
--
Ds (many to one relationship) with something like merge() or
> functions in libraries like diplyr or data.table
>
> Cheers
> Herry
>
> -Original Message-
> From: R-sig-Geo [mailto:r-sig-geo-boun...@r-project.org] On Behalf Of
> Patrick Schratz
> Sent: Wednesday, 24 M
Do also see the raster vignettes, there is one on performance and using
line- or block- based chunks for large jobs.
Cheers, Mike
On Thu, 25 May 2017, 08:32 Michael Sumner <mdsum...@gmail.com> wrote:
> Try crop() (or more abstractly extract(raster, cells) once you've
> determined w
vice is welcome.
>
> Thank you very much for your help.
> Megan
>
> [[alternative HTML version deleted]]
>
> ___
> R-sig-Geo mailing list
> R-sig-Geo@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>
--
Dr. Michael
de
>
> ___
> R-sig-Geo mailing list
> R-sig-Geo@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
--
Dr. Michael Sumner
Software and Database Engineer
Australian Antarctic Division
203 Channel Highway
Kingston Tasmania 7050 Au
1 - 100 of 394 matches
Mail list logo